Efficient Interpolation of the Bethe-Salpeter Hamiltonian at Arbitrary Momenta

技術要約：任意の運動量におけるベテ・サルピーター・ハミルトニアンの効率的な補間

問題提起
ベテ・サルピーター方程式（BSE）は、光学スペクトルやエキシトン結合エネルギーを含む、固体中の中性励起に関する標準的な第一原理的記述を提供する。しかし、その予測精度は、電子・正孔遷移および遮蔽クーロン相互作用を収束させるために高密度な運動量サンプリングが必要となるため、膨大な計算コストによって制限されている。標準的な実装は、規則的な電子運動量メッシュに依存しており、任意の転送運動量（$Q$）における計算を制約的なものにしている。任意の $Q$ に対しては、基礎となる高密度な計算を新たにやり直す必要があり、これは計算コストの面から極めて困難である。この制限は、$\Gamma$ 点付近の二次元系において特に顕著である。なぜなら、長距離の電子・正孔交換相互作用は、典型的な ab initio メッシュの間隔よりもはるかに小さい運動量での非解析的なエキシトン分散を誘起するためである。さらに、エキシトンを格子自由度（例：エキシトン・フォノン相互作用や非平衡ダイナミクス）に結合させるには、エキシトニック状態と相互作用バーテックス（頂点）の両方に対して共通の滑らかな電子表現が必要であるが、異なるサンプリンググリッドを用いる場合にはこれを達成することが困難である。

手法
著者らは、粗い第一原理グリッドから任意の結晶運動量における静的なBSEハミルトニアンを再構成するための逆空間フレームワークを開発した。この手法は、電子・正孔遷移運動量（$k$）および転送運動量（$Q$）を、基礎となる ab initio サンプリングからデカップル（分離）するものである。この手法はコード EPIq に実装されており、以下の4つの主要なステップで進行する。

1. ブロッホ状態のワニエ補間: ブロック状態は、最適に滑らかなワニエ部分空間（OSS）内で補間される。ブロッホ関数の周期部分は平面波に展開され、最大局在化ワニエ関数（MLWF）を用いてOSSへと変換される。これにより、逆空間における滑らかな表現が保証され、中間的な実空間への変換を経ることなく、任意の運動量における波動関数係数の直接的な補間が可能になる。補間された状態は、遷移密度を構築するためにバンド固有状態表現へと回転される。
2. 遮蔽相互作用の再構成: 逆誘電行列（$\epsilon^{-1}$）を直接補間すると、裸のクーロン相互作用による急激な変化が残ってしまうため、本手法ではより滑らかな補助関数 $f_{GG'}(q)$ を補間する。この関数は、明示的なクーロン依存性を除去するものである。その後、必要な交換運動量 $q = k - k'$ において、クーロン因子を再導入することで、遮蔽相互作用 $W_{GG'}(q)$ を解析的に再構成する。裸のクーロン相互作用は、要求される運動量 $Q$ で直接評価され、決して補間されない。
3. ハミルトニアンの構築: 補間された遷移密度と再構成された遮蔽相互作用を用いて、直接および交換電子・正孔カーネルが計算される。BSEハミルトニアンは、任意の遷移運動量 $k, k'$ および転送運動量 $Q$ に対して組み立てられる。
4. 解法: 全対角化を避け、反復ソルバー（光学スペクトルのための Lanczos–Haydock または個別のエキシトン状態のための Davidson）を用いてハミルトニアンを解く。

主要な貢献

• サンプリングのデカップリング: 本手法は、遷移運動量グリッド（$k$）と転送運動量（$Q$）を、初期の ab initio 計算（ブロッホ状態および誘電応答）で使用された粗いソースグリッドから独立して選択することを可能にする。
• 共通の滑らかな表現: ブロッホ状態と誘電応答の両方にOSSを利用することで、本フレームワークは、エキシトン波動関数と相互作用バーテックスが共通の滑らかな電子表現内に存在することを保証する。これにより、エキシトニック効果と電子・フォノン結合のような他の相互作用バーテックスとの一貫した結合が容易になる。
• 任意の運動量へのアクセス: 本フレームワークは、全ブリルアンゾーンの均一な精細化を必要とせずに、エキシトン分散の長波長領域への直接的なアクセスを可能にする。

結果およびベンチマーク
本手法は、単層 hBN、単層 MoS$_2$、およびバルク GaAs について、Yambo を用いた明示的な ab initio 計算と比較してベンチマークを行った。

• 補間精度: 補間されたブロッホ状態の平面波成分および再構成された誘電遮蔽は、大幅に粗いソースメッシュ（例：$12\times12\times1$）から導出された場合でも、高密度なターゲットグリッド（例：$36\times36\times1$ または $48\times48\times1$）を用いた明示的な計算と優れた一致を示した。ブロッホ状態の相対誤差は1%未満であり、誘電遮蔽の誤差は $\Gamma$ 点付近で最大となったが、管理可能な範囲（粗いグリッドで約5%）に留まり、ソースメッシュの精細化に伴って減少した。
• 光学スペクトル: 再構成されたBSE光学スペクトルは、エキシトンのラインシェイプ、ピーク位置、および振動子強度を含む参照スペクトルを密接に再現した。Wasserstein距離を用いた収束解析により、特に低次元系においては、遮蔽ソースメッシュの精細化が精度を決定付ける主要な要因であることが示された。
• 計算効率: 本アプローチは大幅な高速化を実現した。単層 MoS$_2$ の場合、高密度グリッド（$96\times96\times1$）に対して最大 22倍、バルク GaAs（$28\times28\times28$）に対して 10倍 の高速化が推定された。これらの利得は、高価な第一原理計算（NSCF および誘電応答）が粗いソースメッシュ上で一度だけ行われ、補間とBSE構築がより有利にスケールすることに起因する。
• 長波長エキシトン分散: 本手法は、hBN および MoS$_2$ における、実用的な一様グリッドの間隔よりも約7倍小さい運動量におけるエキシトン分散の非解析的挙動を解明した。hBN においては、長距離交換相互作用によって誘起される横方向および縦方向の枝の線形分裂の解明を可能にした。MoS$_2$ では、輝点 A-エキシトン二重線の分裂を解明した。
• 局所的に精細化されたサンプリング: バルク GaAs について、$\Gamma$ 点付近に点を集中させた局所的に精細化された遷移運動量サンプリングを用いることで、わずか数千のBSEポイントのみを使用して、3重縮退した最低エキシトンの結合エネルギーを 4.0 meV まで収束させた。これは実験値（3.8–4.2 meV）および従来のベンチマークと一致しており、一方で一様メッシュでは $10^5$ ポイント以上を必要としても未収束であった。

意義
本論文は、このフレームワークが柔軟な運動量サンプリングをBSE計算において実用的なものにし、あらゆるターゲット運動量に対して高価な ab initio 計算を繰り返す必要性を排除すると主張している。エキシトンのエネルギーと振幅を任意の結晶運動量において共通の滑らかな電子表現内で提供することにより、本手法は、非局所的な電子・フォノンバーテックスへのエキシトニック効果の組み込みや、エキシトン・フォノン相互作用および非平衡エキシトンダイナミクスへの計算の拡張への道を開くものである。著者らは、本手法が、粗いグリッド上のブロッホ状態と微視的な逆誘電行列の利用可能性のみに依存しているため、幅広い材料や電子構造コードに適用可能であることを強調している。