Real-Space Spectral Approach to Orbital Magnetization

技術要約：軌道磁化への実空間スペクトル・アプローチ

問題提起
結晶性固体における軌道磁化の正確な計算は、物性物理学における重要な課題であり続けている。現代的な$k$空間理論では、占有されたブロッホ状態に関わるベリー曲率や量子計量を用いたブリルアンゾーン積分として軌道磁化を表現するが、第一原理計算におけるその実装は技術的に困難である。それはワニエ補間や非絡み合い（disentanglement）の手続きに依存しており、特にバンドが孤立していない場合には、複雑さと不確実性が導入される。さらに、$k$空間の定式化は並進対称性を仮定しているため、無秩序、不均一性、あるいは大きな実空間スーパーセルを持つ系には適用できない。

BiancoとRestaによる理論のような既存の実空間定式化は、局所的な磁気双極子密度のマクロな平均として磁化を表現する代替案を提示している。しかし、これらの手法は通常、基底状態の一体粒子密度行列の構築を必要とし、それは厳密な固有状態および基底状態射影体の知識を必要とする。この要件は、大規模な系へのスケーラビリティを著しく制限し、無秩序な材料に対する密度汎関数理論のワークフローへの効率的な統合を妨げている。

手法
著者らは、ハミルトニアンの対角化や基底状態射影体の構築を必要とせずに、有限サイズの系の軌道磁化を計算する、新しい実空間スペクトル法を導入する。この手法は、有限の二次元系における全軌道磁化演算子 $\hat{M}_z$ の交換子形式に基づいている：
$$M_z = -\frac{ie}{2\hbar c A} \sum_{\epsilon_n < E_F} \langle \phi_n | \hat{r} \times [\hat{H}, \hat{r}] | \phi_n \rangle$$
基底状態射影体を直接評価する代わりに、著者らはエネルギー分解されたスペクトル関数、すなわち磁化スペクトル密度 $m_z(E)$ を定式化する：
$$m_z(E) = \text{Tr}[\hat{M}_z \delta(E - \hat{H})]$$
全軌道磁化は、フェルミ準位までの単純な数値積分によって得られる：$M_z(E_F) = \int^{E_F} m_z(E) dE$。

$m_z(E)$ を効率的に評価するために、本手法はスペクトル演算子 $\delta(E - \hat{H})$ に対してチェビシェフ多項式展開を用いる。ハミルトニアン $\hat{H}$ は、区間 $[-1, 1]$ 内の $\tilde{H}$ に再スケーリングされる。展開係数（モーメント $\mu_n$）は、チェビシェフ再帰手順を用いて反復的に計算される。計算をさらに加速させるため、トレースはランダムベクトルを用いた確率的（Stochastic）な手法によって評価される。すなわち、$\mu_n \approx R^{-1} \sum_{r=1}^R \langle r | \hat{M}_z T_n(\tilde{H}) | r \rangle$ となる。

このフレームワークは、位置演算子に伴う曖昧さを回避するために開境界条件を利用している。計算コストは、システムサイズ、格子サイト数、およびチェビシェフ反復回数に対して線形にスケールするため、明示的な対角化を必要とせずに最大$10^{10}$サイトの系まで扱うことが可能である。

主な貢献と結果
著者らは、軌道磁化の定式化をテストするための典型的な系として、ハニカム格子上のハルデン（Haldane）モデルを用いて本手法を検証している。主な知見は以下の通りである：

1. $k$空間理論とのベンチマーク： 本手法は、トポロジカルおよび自明な絶縁体領域の両方において、現代的な$k$空間定式化の結果と極めて良好に一致する結果を再現した。トポロジカル相においては、バルクギャップ内での $M_z$ のフェルミエネルギー ($E_F$) に対する線形依存性を正しく捉えており、これはカイラルなエッジ状態のスペクトルフローの徴候である。自明な相においては、電子・正孔対称性に応じて、消失またはプラトー（平坦化）した応答を正しく予測している。
2. 無秩序と欠陥： 本アプローチは、無相関なアンダーソン無秩序および点欠陥（空孔）の存在に対しても堅牢であることが示されている。
   • アンダーソン無秩序： トポロジカルギャップにおける $M_ز$ の線形依存性は、移動度ギャップが存続する限り強固である。自明な絶縁体においては、ギャップを埋める無秩序誘起状態が、最終的に磁化のプラトーを破壊する。
   • 空孔欠陥： 空孔は状態密度に鋭い共鳴をもたらし、磁化曲線に局所的なキンクを引き起こすが、チャーン数に関連するグローバルな線形傾向は維持される。これは、空孔共鳴が局在しているものの、マクロな磁化を決定する長距離の軌道循環に対して弱く結合していることを示唆している。
3. チャーン数の直接抽出： 著者らは、エネルギー微分とスペクトル密度を介して、軌道磁化のエネルギー微分とバルク・チャーン数 ($C$) との間の直接的な関係を確立した：$-\frac{2\hbar c}{e} m_z(E_F) = \frac{C}{2\pi}$。これにより、射影体、曲率マーカー、またはワニエ関数を用いることなく、磁化スペクトル密度からチャーン数を直接抽出することが可能となる。本手法は、トポロジカル相転移（$C=0 \to 1$）を正常に解明し、システムサイズが増大するにつれて熱力学的極限の値へと収束する。
4. ゲージ不変性と境界条件： 本定式化は、有限の試料のゲージ不変な全磁化を対象としている。著者らは、局所的な磁化の空間マップはゲージや原点の選択に依存する一方で、グローバルな観測量 $M_z$ は明確に定義されていることを明らかにしている。本手法は、人工的なバルクのみの平均化を必要とせず、バルクと境界の寄与を等しく扱うことができる。

意義
本論文は、このチェビシェフ・スペクトル法が、大規模で無秩序なトポロジカル系を調査するための、コンパクトで計算効率が高く、かつ広く適用可能なフレームワークを提供すると主張している。固有状態や基底状態射影体を回避することで、本手法は既存の実空間技術のスケール性の限界や、$k$空間理論の並進対称性の制約を克服している。

著者らは、このアプローチが、強い無秩序や有限温度を含む、現実的な複雑さを持つ量子材料の実空間モデルにおける、軌道応答やトポロジカル転移を調査するための有望な経路を開くものであると断じている。これにより、従来のメソッドではアクセスできなかったはるかに大規模で現実的な系における軌道磁化の研究が可能となり、明示的な対角化を欠く第一原理計算のワークフローへの軌道磁化計算の統合を促進する。