Dynamical Hubbard approach to correlated materials: the case of transition-metal monoxides

技術要約：相関物質に対する動力学的ハバード・アプローチ

問題提起
遷移金属モノオキサイド（MnO, FeO, CoO, NiO）は、モット絶縁体および電荷移動絶縁体に分類される複雑な電子相関を研究するための典型的な系として機能する。標準的なコーン・シャム密度汎関数理論（DFT）と半局所的汎関数は、CoOおよびFeOの絶縁状態を予測できず、MnOおよびNiOにおけるバンドギャップを過小評価するが、より高度な手法は重大なトレードオフに直面している。

• 静的な平均場手法: DFT+Uやハイブリッド汎関数などの手法は、絶縁状態とギャップの大きさを回復できるが、本質的に静的である。これらは、バンド再正規化、スペクトル重みの移動、粒子寿命といった、電子相関の重要な徴候である周波数依存の自己エネルギーを捉えることができない。
• 多体摂動論 (MBPT): $G_0W_0$ や $GW_0$ といった手法は、金属的なDFT基底状態から開始する場合に苦戦することが多い。QSG$\hat{W}$ のような拡張手法はギャップの記述を改善するが、計算コストのために完全な周波数依存性を放棄することが多く、$k$ 空間分解されたスペクトル重みの移動を記述する能力を制限している。
• DFT+DMFT: 動的平均場理論（DFT+DMFT）は、相関部分空間をアンダーソン不純物モデルにマッピングすることで、局所的な動的相関と揺らぐ磁気モーメントをモデル化することに成功している。しかし、これには高い計算コストと数値的な複雑さが伴い、厳密解はシステムサイズに対して指数関数的にスケールする。さらに、近年のDFT+DMFTの結果は、4つのモノオキサイドすべてに対して広いギャップを予測しており、NiOにおける特定の分散バンドのような、角度分解光電子分光（ARPES）で観察される特定のスペクトル特徴を再現できていない。

手法
著者らは、スピン偏極したマルチサイト磁性系を扱うために、以前に導入した定式化を一般化した「動力学的ハバード（dynH）」汎関数を提案している。この手法は、動的な相関を捉えつつ、数値的に効率的かつ物理的に透明であるように設計されている。

• 定式化: この手法は、局所的な周波数依存の遮蔽相互作用ポテンシャルを導入することにより、Dudarevの回転不変なDFT+U定式化を拡張する。自己エネルギーは、局所グリーン関数に関する動力学的ハバード汎関数の関数微分から導かれる：
  $$\Sigma^{\sigma, I}_{\text{dynH}}(\omega) = -\int d\omega' U_I(\omega') G^{\sigma}_I(\omega + \omega') + \frac{U^{\infty}_I}{2}$$
  ここで、$U_I(\omega)$ は（RPAスピンDFT分極によって計算された）局所的な遮蔽クーロン相互作用の平均であり、$U^{\infty}_I$ は遮蔽されていない裸のクーロンポテンシャルである。
• 実装:
  • 局所自己エネルギーは、ディゾン方程式（$G^{-1} = G^{-1}_0 - \Sigma$）を解くために$k$空間へとアップフォールドされる。
  • 解法には、**極の和に基づくアルゴリズム的反転法（AIM-SOP）**が用いられる。この技術は、局所的な遮蔽相互作用とグリーン関数の両方を、一次の極の和として表現する。
  • 決定的なことに、AIM-SOPは実周波数軸上で直接動作するため、標準的なDMFT不純物ソルバーで必要とされる虚数松原周波数や、その後の解析接続を必要としない。これにより、計算精度が向上する。
  • 基底状態の計算には r2SCAN 汎関数が採用されており、これは4つの化合物すべてに対して絶縁相を等しく正しく予測する。
  • 有効な遮蔽相互作用は、一次摂動論で計算され、ゼロ周波数極限は $U(0) = 5.3$ eV である。

主な結果
本フレームワークは、MnO, FeO, CoO, NiOの基底状態に適用され、実験的な光電子分光（PES）および逆光電子分光（IPES）データ、ならびに最先端のDFT+DMFTおよびQSG$\hat{W}$ 計算と比較された。

• MnO: 被占有状態密度（DOS）はDFT+DMFTの結果と一致しており、実験と比較して全体的なバンド幅がわずかに収縮している。5 eVにおけるLUMO励起は、位置、幅、強度のすべてにおいてIPESデータと一致している。高エネルギーのサテライトについては依然として未解明であり、これはDFT+DMFTやQSG$\hat{W}$ が共有する限界である。
• FeO: 本手法は約2 eVの小さなバンドギャップを予測しており、これはIPESデータと一致している。これは、4つのモノオキサイドすべてに対して広いギャップを予測するDFT+DMFTおよびQSG$\hat{W}$ と対照的である。第一励起後の強度の減少もよく記述されている。
• CoO: DOSはDFT+DMFTおよびPESの両方とよく一致している。空状態の分布に関して顕著な違いがあり、DFT+dynHは$t_{2g}$ と $e_g$ 対称状態に対して同等の重みを示すが、DFT+DMFTは重みを主に高エネルギー状態へと移動させる。
• NiO: このアプローチは、実験的なARPESデータおよびPES/IPESスペクトルと極めて良好な一致を示す。
  • 全体的なバンド幅（6 eV）はARPESと一致している。
  • 本手法は、Zhang-Rice二重項束縛状態および平坦なNi(3d)/O(2p)バンドを正常に再現している。
  • 決定的な知見: 特定のバンドに対して消失するスペクトル重みを予測するDFT+DMFTとは異なり、DFT+dynHは -8 eVにおける分散バンド および -4 eV付近の反強磁性的（AFM）平坦バンド を再現している。これらの特徴は実験的なARPESデータに存在するが、以前のDMF不純物解では欠落していたものである。著者らは、AFMバンドのわずかな過小評価について、r2SCANの開始点、自己整合性の欠如、または二重計上項の近似性に起因すると考えている。

意義および主張
本論文は、DFT+dynHアプローチが、計算コストの高いDFT+DMFTのような手法に代わる「単純で透明な」選択肢を提供すると主張している。その主な貢献は以下の通りである：

1. 一般化: マルチサイト、スピン偏極磁性系に対して動力学的ハバード汎関数を拡張することに成功した。
2. 精度とコスト: スペクトル特性（バンド再正規化やスペクトル重みの移動を含む）において、最先端の理論（DFT+DMFT, QSG$\hat{W}$）と同等またはそれを上回る結果を、大幅に低い計算コストと数値的複雑さで達成した。
3. 不一致の解消: FeOにおける小さなギャップの予測や、DMFT解では欠落していたNiOのスペクトル特徴の回収など、従来の理論の失敗を修正した。
4. 物理的洞察: 厳密な不純物ソルバーの膨大なスケーリングを伴うことなく、電子相関の根本的な徴候（粒子寿命、サテライト、スペクトル重みの移動）へのアクセスを提供した。

著者らは、DFT+dynHを、対称性の破れた軌道秩序やハイブリダイゼーションを持つ現実的で複雑な相関物質をシミュレートするための有望なツールとして位置づけている。