Markov Inequality as a Tool for Linear-Scaling Estimation of Local Observables

技術要約：局所観測量の線形スケーリング推定のためのツールとしてのマルコフ不等式

問題提起
大規模な無秩序格子系における、局所状態密度（LDoS）、局所電流、トポロジカルマーカーといった空間分解された量子観測量を明らかにすることは、依然として重大な計算上の課題である。チェビシェフ展開に基づく線形スケーリング・アルゴリズムは、体積平均されたスペクトル関数（演算子のトレース）を効率的に推定できるが、フルサンプル（全領域）の空間マップへの適用には苦慮する。これらのマップの生成は、局所的な量には自己平均性が欠如しているため、通常、$O(N^2)$（$N$はシステムサイズ）のスケーリングによる計算コストを伴う。アンダーソン局在や、金属ー絶縁体転移付近でのマルチフラクタル・スケーリングを示すような、強い空間的ゆらぎを持つ系では、標準的な確率論的サンプリング手法は一様に収束しない。具体的には、素朴な（naive）推定量の相対的なサンプリング誤差は、局所的な波動関数の振幅（$|\psi_\alpha(r)|^{-1}$）に反比例してスケールするため、非一様な収束と非物理的な結果（例：負のLDoS値）をもたらす。

手法
著者らは、タイトバインディングモデルにおける任意の正定値局所スペクトル演算子の実空間マップを、制御された精度で計算するために設計された、線形スケーリングの確率論的手法を導入する。核心となる革新は、正定値確率論的推定器とマルコフ不等式の組み合わせにある。

1. 正定値推定器： 符号の変化を伴う可能性のある従来の非バイアス推定器とは異なり、提案手法は厳密に非負となる推定器を構築する。正定値演算子 $X$ に対して、ランダム場 $\phi_r = \langle \xi | \circ \sqrt{X} | \xi \rangle$ が導入される。ターゲットとなる対角要素 $X_r$ は、$|\phi_r|^2$（複素ベクトルの場合）、または特定のモーメント条件を満たす実ベクトルに対する $\text{Var}(\phi_r)$ として推定される。
2. マルコフ不等式による境界： 推定器が非負であることを保証することで、マルコフ不等式（$P(\eta_r \ge C) \le \eta_r C^{-1}$）を用いて、各サイトにおける局所平均に対するサンプリング誤差を厳密に束縛する。これにより、局所的な観測量の大きさに関わらず、相対的なゆらぎが全サンプルにわたって一定に保たれ、均質な収束が保証される。
3. 非対角観測量への拡張： この手法は、局所ユニタリ変換を通じて非対角要素（例：局所電流）へと拡張される。サイトのペア $\{r, r'\}$ に対して局所的な回転を適用することで、非対角行列要素を回転した基底の対角要素へとエンコードする。これにより、回転が互いに素なサイトペアに対して作用するため、追加のサンプリング相関を誘発することなく、正定値推定技術を非対角量に適用することが可能となる。

主な結果
本手法の有効性は、正方格子上の無秩序2D $\pi$-フラックスモデルを用いたベンチマーク計算を通じて実証されている。

• 単一空孔ケース： $4096 \times 4096$ のスーパーセル内に中央に単一の空孔が存在する場合、本手法はゼロエネルギーLDoSの計算に成功した。それは、数桁にわたる期待される代数的減衰（$\rho_r(0) \propto r^{-2}$）を忠実に捉えた。対照的に、従来の素朴な推定法は低振幅領域で失敗し、非物理的な負の値や大きな相対誤差を生じた。
• 複数空孔： $2048 \times 2048$ のスーパーセルにおいて、0.5%の空孔濃度を持つ系に対し、本手法は数桁に及ぶ収束したLDoSマップを生成した。局所相対ゆらぎの95パーセンタイルは3.5%に抑えられ、臨界ゼロエネルギー状態の存在下でも制御された精度が示された。
• 誤差解析： $512 \times 512$ のサンプルを用いた厳密なLDoS計算との比較により、正定値推定器は、素朴な推定器と比較して、局所相対誤差のモードを約2桁減少させることが示された。誤差分布は、ランダムベクトル $R$ の数に対して $R^{-1/2}$ としてスケールした。計算量については、確率論的推定に約4 CPU時間要したのに対し、厳密な計算には約171 CPU日を要した。
• 局所電流： 本手法は、長距離ポテンシャル散乱体の存在下における定常状態の局所電流場の計算に適用された。得られたマップは、複雑な渦–反渦ペアと複雑な流れのパターンを明らかにし、これまで高い空間分解能での捕捉が困難であったグラフェンの局在転移（コステリッツ＝ザウレス普遍クラス）に関連する構造を解明した。

意義および主張
本論文は、このアプローチが、計算複雑性を二次スケーリングから線形スケーリングへと低減させることで、現実的に大規模な格子における無秩序駆動型のメゾスコピック現象のシミュレーションを可能にすると主張している。マルコフ不等式によるサイトごとの確率論的誤差の厳密な束縛は、局所量における自己平均性の欠如という根本的な限界を克服するものである。

著者らは、本手法が以下のことを実現すると断言している：

• 強力な空間的ゆらぎの下でも、局所的な観測量の正確な推定を可能にする。
• 大規模系におけるマルチフラクタリティーや無秩序駆動型転移の研究を促進する。
• 効率的かつ制御された精度の局所入力を提供することで、実空間の自己整合的平均場計算を加速させる。
• 任意の正の局所スペクトル演算体に適用可能であり、ユニタリ変換を通じて非対角観測量にも拡張可能な一般的なフレームワークを提供する。

本研究は、実験的な空間分解能（STM、SQUIDなど）と、凝縮系物理学における理論的シミュレーション能力との間の溝を埋めるためのツールとして提示されている。