Theory of local orbital magnetization: local Berry curvature

技術要約：局所軌道磁化および局所ベリー曲率の理論

問題提起
周期的な結晶固体における軌道磁化の整合的な量子論は2005年に確立されたが（軌道モーメントとベリー曲率を通じて磁化を定義）、多様な幾何学的構造において有効な局所軌道磁化の微視的理論はこれまで存在していなかった。結晶に関する現代的な理論などの既存の枠組みは、単位胞の平均値に依存しており、サブ格子や単一サイトレベルの磁気テクスチャを不明瞭にしている。さらに、Středaの公式は周期系において軌道磁化とホール伝導度を結びつけているが、トポロジカルな情報が有限系、リボン、あるいは非結晶構造においてどのように局所的にエンコードされるのかは不明である。Bianco-Restaの局所Chernマーカーのような従来の実空間アプローチは、グローバルな全磁化から局所磁化を定義しており、表面電流やヒンジ効果に関連するゲージの曖昧さを導入している。本研究は、周期結晶、リボン、および有限の開境界系にわたって、原子スケールで軌道磁化を解決する、ゲージ不変な熱力学的理論の必要性に対処するものである。

手法
著者らは、磁場依存局所状態密度（LDOS）、$\rho(r, \varepsilon, B)$ の摂動展開に基づく熱力学理論を開発した。

• 定式化: 局所グランドポテンシャル $\Xi(r, \mu, B) = -\int d\varepsilon F(\varepsilon)\rho(r, \varepsilon, B)$ から出発し、局所軌道磁化は零磁場極限 $M(r, \mu) = -\partial_B \Xi(r)|_{B=0}$ として定義される。
• モデル: 理論は、ペイヤーズ置換を通じて磁場を取り入れた、非相互作用のスピンレス電子のタイトバインディングモデルに対して定式化されている。LDOSは、ゲージ不変な摂動論を用いて評価される。
• 導出: 中心となる結果は、遅延グリーン関数 $G(\varepsilon)$ と速度演算子 $\hat{v}$ を含む局所磁化の式である：
  $$M_{orb}(r, \mu) = \frac{e}{2\pi\hbar} \int d\varepsilon F(\varepsilon) \text{Re}\langle r | G [(\hat{v}G) \times (\hat{v}G)] | r \rangle$$
  この式は、有限系（開境界）、無限周期結晶（サブ格子分解）、およびリボン幾何学という3つの異なる幾何学的構造に対して評価される。

主要な貢献と結果

1. 局所ベリー曲率の特定:
   本理論は、これまで特定されていなかった量、すなわち局所ベリー曲率 $\Omega_n(r)$（有限系の場合）、または $\Omega_{nk}(r_\alpha)$（サブ格子の場）を明らかにした。

   • 有限系において、局所磁化は、射影された軌道モーメント $m_n(r) = |\langle r|n\rangle|^2 m_n$ と局所ベリー曲率項 $F(\varepsilon_n)\Omega_n(r)$ の2つの項に分離される。
   • 物理的には、軌道モーメントが（位置に依存しない）エネルギー準位の磁場によるシフトを記述するのに対し、局所ベリー曲率は、実空間における電子スペクトル重みの磁場誘起再分配を支配する。これは $\sum_r \Omega_n(r) = 0$ という総和則を満たしており、状態の全正規化を変更することなく、重みを再分配することを示す。

2. 幾何学を横断する統一的記述:
   本定式化は、以下のすべてに対して一貫した記述を提供する：

   • 有限系: 局所磁化は、すべてのサイトにわたる総和によって従来の全軌道磁化を回収する。このとき、局所ベリー曲率の寄与は相殺される。
   • 周期結晶: 理論は磁化をサブ格子スケールで分解する。局所ベリー曲率は、従来のブロッホ・ベリー曲率の射影と、純粋に幾何学的な寄与 $\Omega_{nk}^{geom}(r_\alpha)$ に分解される。この幾何学的項は、単位胞内のサブ格子間でスペクトル重みを再分配し、$\sum_\alpha \Omega_{nk}^{geom}(r_\alpha) = 0$ を満たす。
   • リボン: 同様の軌道モーメントと局所ベリー曲率への分離が保持され、これら3つの幾何学すべてに対する統一的な枠組みを提供する。

3. トポロジカルなエンコーディングと勾配の量子化:
   著者らは、絶縁体ギャップ内における局所軌道磁化の勾配が量子化されており、リボンや有限の幾何学においても局所的にStředaの関係 $\partial_\mu M = \frac{e}{h}C$ を満たすことを示している。

   • この勾配は、局所ベリー曲率によって支配される。
   • 2バンド系において、この勾配は、サブ格子に依存しないトポロジカルな寄与（Chern数の半分）と、$\Omega_{nk}^{geom}$ から生じるサブ格子依存の幾何学的寄与で構成される。
   • 幾何学的寄与は粒子・正孔対称な系では消失するが、自明なギャップにおいては顕著であり、サブ格子レベルでの異なる軌道磁気秩序（強磁性、反強磁性、フェリ磁性）を区別する。

4. Bianco-Restaアプローチとの比較:
   本論文は、その結果をBianco-Restaの局所Chernマーカーと比較している。

   • 計算効率: 新しい定式化は、$N^3$ の計算量を要するBianco-Restaに対し、$N^2$ でスケールする。
   • 概念的な違い: Bianco-Restaのアプローチは、グローバルな全磁化から局所磁化を構築するため、表面磁化に関連するゲージの曖昧さが生じる。現在の理論は、磁場依存LDOSから直接局所応答を導出しており、一意でゲージ不変な局所磁化の定義を提供する。
   • 物理的区別: 両アプローチは著しく異なるサブ格子テクスチャを与える。決定的なことに、Bianco-Restaのアプローチは、自明な相において顕著となる局所ベリー曲率の幾何学的部分に関連する、ギャップ内の磁化の勾配を説明できない。

5. バルク・トポロジーの解剖学:
   理論は、トポロジカルな応答をバンド状態とギャップ状態の寄与に分解する。

   • 局所的（バルク内）: トポロジカルな応答は、完全にバンド状態の局所ベリー曲率に由来する。ギャップ状態は局所的なバルクのトポロジカル応答には寄与しない。
   • マクロ的: 量子化された勾配は、幾何学に応じて異なる微視的メカニズムを通じて回収される。リボンでは、それはバンド状態のベリー曲率によって運ばれる。有限の開境界系では、ベリー曲率の寄与は総和によって消失し、量子化された勾配はギャップ状態の軌道モーメントを通じて回収される。

意義
本論文は、局所軌道磁化を、単位胞平均の量から原子スケールへと拡張することで、真の微視的観測量として確立したと主張している。局所ベリー曲率を特定することにより、本研究は有限系におけるバルクのトポロジーに関する記述を提供し、トポロジカルな応答がエッジ状態とは独立してバルク内に局所的にエンコードされ得ることを明らかにしている。この枠組みは、単一サイトおよびサブ格子スケールでの軌道磁気テクスチャを解決し、異なる微視的メカニズム（軌道モーメント vs ベリー曲率による再分配）がいかに協力してマクロな量子化応答を生み出すかについて、統一的な熱力学的視点を提供する。これは、分子、非晶質材料、準結晶、モアレ構造など、並進対称性を欠く系における軌道磁気テクスチャの研究への道を開くものである。