[論文レビュー] Wannier Functions Dually Localized in Space and Energy

Aaron Mahler, J. Williams+2 more|arXiv (Cornell University)|Jan 19, 2022

Advanced Chemical Physics Studies被引用数 4

ひとこと要約

本稿では、空間的分散とエネルギー分散の重み付き和を最小化することで、実空間およびエネルギー空間の両方で局在化する二重局在Wannier関数（DLWFs）を導入する。この手法により、分数占有軌道を含む特定のエネルギー準位に関連するWannier関数が得られ、占有状態と非占有状態の両方の前線軌道を一元的かつ一貫的に記述可能となる。これは、密度汎関数理論（DFT）における局在化誤差を是正するにあたり、占有状態と非占有状態を別々に取り扱う必要がないという点で重要である。

ABSTRACT

The construction of Wannier functions from Bloch orbitals offers a unitary freedom that can be exploited to yield Wannier functions with advantageous properties. Minimizing the spatial variance is a well-known choice; another, previously proposed for Wannier functions constructed from the occupied Bloch manifold, minimizes a weighted sum of spatial and energy variance. Departing from all previous work, we extend dual localization to include both valence and conduction bands together. Near the Fermi energy, these dually localized Wannier functions yield frontier (bonding and antibonding) orbitals in bulk silicon and molecular ethylene, as well as $d$-orbital character in metallic copper. Because they are both localized and retain information about the orbital energy spectrum, dually localized Wannier functions are well suited to orbital-dependent methods that associate Wannier functions with specific energy ranges. They naturally induce fractional occupations, allowing for corrections to the DFA total energy.

研究の動機と目的

• 従来の最大局在Wannier関数（MLWFs）の制限を克服し、空間的およびエネルギー的両方で局在化するWannier関数を構築する統一的枠組みを開発すること。
• 占有状態と非占有状態の両方の前線軌道を、別々の局在化手順をとらずに同一のWannier関数枠組み内で一貫して記述する課題に取り組むこと。
• 1つのWannier関数セットを用いて、シリコンやエチレンのような系における化学的結合性および反結合性の特徴を正確に表現できるようにすること。
• 軌道依存DFT法を支援するため、特定のエネルギー準位に結びついた内因的分数占有を持つWannier関数を提供すること。
• バンド構造が複雑に交差する系における逐次的脱結合化と局在化の限界を克服し、コア最適化にエネルギー局在化を統合すること。

提案手法

• 空間分散（Ω）とエネルギー分散（E）のバランスを取るためのコスト関数 F = (1−γ)Ω + γE を最小化する。γは空間的局在化とエネルギー的局在化のトレードオフを制御する。
• 各kポイントでバンドのユニタリ混合を許容する一般化されたWannier関数構築法を用い、占有状態と非占有状態の両方のバンドを同一のWannier関数セットに含めることを可能にする。
• シリコンとエチレンで妥当性が確認された、Loschmidt（LOSC）コスト関数（γ = 0.47714）を用いて、空間的局在化を保ちつつ最適なエネルギー局在化を達成する。
• k空間のBloch軌道のフーリエ変換によりWannier関数を構築し、位相およびユニタリ自由度を最適化してコスト関数Fを最小化することで、滑らかなk依存性と局在化を確保する。
• 等値面および等値面閾値（例：160, 60, 30）を用いてDLWFの空間的分布を可視化し、結合性・反結合性の特徴を特定する。
• エネルギー局在化を活用してWannier関数に分数占有を割り当て、DFTにおける局在化誤差是正法に直接応用可能にする。

実験結果

リサーチクエスチョン

• RQ1 Wannier関数を空間的およびエネルギー的両方で同時に局在化させることで、占有状態と非占有状態の両方の電子状態を統一的に記述できるか？
• RQ2 Wannier関数におけるエネルギー局在化は、σおよびπ結合性・反結合性状態といった化学的に意味のある軌道を表現する能力にどのように影響するか？
• RQ3 分数占有を持つ二重局在Wannier関数は、シリコンやエチレンのような系における前線軌道の記述をどの程度改善できるか？
• RQ4 統一的局在化アプローチは、バンドが複雑に交差する系における逐次的脱結合化と局在化の落とし穴を回避できるか？
• RQ5 占有数が1から大きく離れたWannier関数が顕著な虚数成分を示す理由は何か？また、これは物理的解釈にどのように影響を及ぼすか？

主な発見

• γ = 0.47714で構築されたDLWFsは、エチレンにおけるσおよびπ結合性・反結合性軌道を的確に回復し、それらの化学的特徴を反映した明確な等値面パターンを示した。
• 本手法により、特定のエネルギー準位に自然に結びついた分数占有Wannier関数が得られ、占有状態と非占有状態の空間を別々に取り扱わずとも、前線軌道を直接記述可能となった。
• 占有数が1に近いDLWFsは実数であり、原子的またはタイトバインディング的性質を保ち、対称性と局在化が維持される一方、非整数占有のものは顕著な虚数成分を示した。
• エネルギー局在化により、広く離れたバンド（例：占有状態と非占有状態）を同一のWannier関数セットに含めても物理的に不自然な混合が生じず、表現の物理的整合性が向上した。
• 本アプローチは、逐次的脱結合化と局在化の代替として統一的かつ一貫した手法を提供し、バンドが複雑に交差する系において、劣悪な解を避ける可能性を示した。
• LOSCコスト関数（γ = 0.47714）は、ダイヤモンド格子をとるシリコンおよびエチレンの両方において、化学的に直感的なWannier関数（π*反結合性軌道を含む）を生成し、広範な適用可能性を示した。