$G_0W_0$@HF and BSE methods in periodic systems from Hartree-Fock theory: gaussian orbital and density fitting approach

本論文は、半導体および酸化物におけるハートリー・フォック法によるバンドギャップおよび価電子帯幅の過大評価を、プラズモン極近似なしの厳密なRPAスクリーニングと仮想状態のハイブリッド収束戦略を採用することで補正するガウス型軌道および密度近似を用いた周期系向けの$G_0W_0$@HFおよびベテ・サルペーター方程式の枠組みを提示する。

要約

ダイヤモンドやシリコンの断片のような固体材料が、光を浴びたときや電気が流れたときにどのように振る舞うかを理解しようとしていると想像してください。これを行うために、科学者たちは材料内部の電子の正確なエネルギー準位を計算する必要があります。これらのエネルギー準位を、電子が住む超高層ビルの「階」と考えてみてください。階がどこにあるかを正確に知れば、その建物がどのように機能するかを知ることができます。

長年にわたり、これらの階をマッピングする標準的な方法は、密度汎関数理論（DFT）と呼ばれる手法を用いるものでした。しかし、DFT は少しぼやけた地図を使うようなもので、建物の全体的な形状は正しく捉えますが、階の高さを正確に捉えることはよくありません。より鮮明な画像を得るために、科学者たちはGW（方程式中の記号 G と W にちなんで命名）と呼ばれるより高度な技術を使用します。この手法は、ぼやけたスケッチから高解像度の 3D モデルに切り替えるようなものですが、計算コストが極めて高く、通常、特定の種類の材料では扱いにくい数学的な「格子」（平面波と呼ばれる）を必要とします。

新しいアプローチ：異なるレンズ
チャールズ・H・パターソンによって書かれたこの論文は、その高解像度の 3D モデルを構築する新しい方法を紹介します。標準的なぼやけた地図（DFT）を出発点とするのではなく、著者は**ハートリー・フォック（HF）**と呼ばれる、非常に鋭いが過度に硬直した別の地図を出発点とします。

• 出発点の問題点： ハートリー・フォック法は、厳しすぎる定規で描かれた地図のようなものです。これは階が離れすぎている（「バンドギャップ」が大きすぎる）と予測し、部屋が広すぎる（「バンド幅」が大きすぎる）と予測します。この地図だけをそのまま使えば、予測は誤ったものになります。
• 解決策： 著者は巧妙な戦略を用います。彼らはこの厳格なハートリー・フォック地図を出発点とし、誤りを修正する「補正レンズ」（GW 法）を適用します。この論文は、この補正レンズが実際には広すぎる部屋を実際のサイズに縮小させるのに非常に優れており、最終的な地図が実験的な現実と非常に良く一致することを示しています。

ツール：ガウス型軌道と密度フィット
ほとんどの GW 計算は、電子を記述するために「平面波」（無限の平らなシートで構成された格子のようなもの）を使用します。この論文では、代わりにガウス型軌道を使用します。

• アナロジー： 複雑な彫刻を記述することを想像してください。平面波アプローチは、何百万もの平らな正方形のタイルを積み重ねてそれを記述しようとするようなものです。一方、ガウス型アプローチは、彫刻の曲線に完璧にフィットするように成形できる、柔らかい丸い粘土の塊を使用するようなものです。これは複雑な分子や結晶に対して、しばしばより効率的です。
• 密度フィット： コンピュータがクラッシュすることなく、これらの粘土の塊を使って数学を機能させるために、著者は密度フィットと呼ばれる手法を使用します。これは「圧縮アルゴリズム」と考えてください。すべての粘土の塊のペア間の相互作用を計算する（これには永遠にかかるでしょう）代わりに、この方法はそれらをクラスターにグループ化し、グループ全体の相互作用を計算します。これは、一人ひとりを個別に体重を量るのではなく、代表的な数人を体重を量ってそれを乗算することで群衆の体重を推定するようなものです。

「近似なし」のトリック
これらの計算における一般的なショートカットに「プラズモンポール近似」というものがあります。

• アナロジー： ドラムの音を予測しようとしていると想像してください。ショートカット法は、「残りを無視して、ドラムが特定の 1 つの音だけを出すものと仮定しましょう」と言います。これは速いですが、ニュアンスを見逃します。
• 論文の主張： この論文はそのショートカットを回避します。それが単一の音に過ぎないと仮定することなく、電子相互作用の完全な複雑な音（完全な周波数依存性）を計算します。これはより正確ですが、材料構造のすべての点に対して巨大で複雑なパズル（ベテ・サルペッター方程式）を解くことを必要とします。

彼らは何を見つけたか
著者は、この新しい方法をダイヤモンド、シリコン、酸化マグネシウム（MgO）、二酸化チタン（TiO2）の 4 つの材料でテストしました。

1. ダイヤモンドとシリコン： 標準的なハートリー・フォック法は、「部屋」（価電子帯）が約 25% 広すぎると予測していました。新しい方法はこれを修正し、実験で測定されるものと正確に一致するようにそれらを縮小しました。
2. 酸化物（MgO と TiO2）： この方法は、エネルギーギャップ（階の間の距離）と材料が光を吸収する方法を成功裡に予測しました。予測されたギャップは実験で観測されるものよりわずかに大きかった（この分野では一般的な問題ですが）、エネルギーマップの全体的な形状は非常に正確でした。
3. 光吸収： これらの材料が光を吸収する方法（その「光学スペクトル」）をシミュレートしたとき、この方法はピークの位置（吸収される色）を非常に良く再現しました。ただし、酸化物については、この方法は光吸収がわずかに強すぎると予測しました。これは、マイクが少し大きすぎる音を拾うようなものです。

結論
この論文は、厳格な「ハートリー・フォック」モデルを出発点とし、洗練された「GW」補正を適用しながら、柔軟な「ガウス」数学的言語と賢明な「圧縮」技術（密度フィット）を使用することで、固体中の電子エネルギーの非常に正確で高解像度の地図を構築できることを実証しています。優れた結果を得るために標準的な「平面波」格子を必要としないことを証明しており、実際には、この代替アプローチは出発点の手法の特定の誤りを修正し、現実の実験と一致する結果を生み出すことができます。

技術概要

技術的概要：ハートリー・フォック理論に基づく周期系における GoWo@HF および BSE 法

問題提起
ab-initio GW 法およびベテ・サルペター方程式（BSE）法は、結晶材料の準粒子（QP）エネルギーおよび光学特性を計算するための標準的なツールである。しかし、既存の実装の大部分は、平面波（PW）基底セットおよび密度汎関数理論（DFT）ハミルトニアン（通常は LDA または PBE）を出発点として依存している。これらは効果的ではあるが、これらのアプローチはしばしば大きな基底セットと、ダイ電関数の周波数依存性を効率的に処理するためのプラズモンポール近似（PPA）などの特定の近似を必要とする。さらに、DFT 出発点は、多体補正から分離することが困難なバンド構造における特定の誤差を導入しうる。量子化学で標準であり、大きな単位格子や開放体積を持つ系に有利であるガウス型軌道（GO）基底セットへのこれらの手法の検証と拡張、ならびに DFT ではなくハートリー・フォック（HF）理論を出発点とする有効性の探求が必要とされている。

手法
著者らは、ガウス型軌道基底セットおよび密度フィット（恒等式分解）アプローチを用いた周期系に対する $G_0W_0$（$GoW_0$ と表記）および BSE 計算を行うための枠組みを提示する。この実装は Exciton コードを利用し、以下の主要な技術的特徴によって特徴づけられる：

• ハートリー・フォック出発点： 固体状態 GW 応用のほとんどとは異なり、摂動を受けていないハミルトニアンおよび単一粒子軌道は HF 理論から導出される。著者らは、HF が通常バンドギャップおよび価電子帯幅を過大評価することを認めるが、その後の $GoW_0$ 自己エネルギー補正がこれらの量を再規格化して実験データと一致させることができると主張する。
• 密度フィット： 周期系における 4 中心積分の計算コストを管理するために、この手法はクーロン計量を用いた密度フィットを採用する。波動関数の積は補助ガウス基底 onto 射影され、積分の複雑さが低減される。クーロン行列はこの補助基底で表現され、その逆行列は遮蔽相互作用の構築に用いられる。
• プラズモンポール近似なしの RPA： 遮蔽クーロン相互作用 $W$ は、$W = v + v\Pi v$ の関係を用いて計算され、ここで $\Pi$ はランダム位相近似（RPA）レベルで扱われる分極率である。決定的な点として、この手法はブリルアンゾーン内の固有の有限波数ベクトル（$Q$）において RPA 方程式（大規模な RPA ハミルトニアンの対角化）を解く。これによりプラズモンポール近似を不要とし、RPA および基底セットの制限内で $W$ の周波数依存性を厳密に扱うことが可能となる。
• ハイブリッド自己エネルギー戦略： 仮想状態の数に対する自己エネルギーの収束に対処するため、ハイブリッド戦略が採用される。自己エネルギーの大部分は、RPA 遮蔽に十分な限られた仮想状態セットを用いた $GoW_0$ によって計算される。高エネルギー仮想レベルからの残りの寄与は、2 次摂動理論（$\Sigma^{(2)}$）を用いて評価される。このアプローチは計算コストと収束のバランスを取る。
• BSE および TDA： ベテ・サルペター方程式は、線形応答の運動方程式アプローチを用いて定式化される。共鳴項と反共鳴項間の結合を無視するタマン・ダンコフ近似（TDA）が採用され、これは光学吸収特性に対して十分であるとみなされる。
• 小 Q 極限： 自己エネルギーおよび遮蔽電子 - 正孔相互作用における発散を処理するため、$Q \to 0$ 極限に対して特別な処理が施され、立方晶固体における交換エネルギーに用いられる手法と類似の手法が利用される。

主要な結果
これらの手法は、元素半導体（ダイヤモンド、ケイ素）および広帯域幅酸化物（MgO、アナターゼ、およびルチル $TiO_2$）に適用された。

• バンド構造の再規格化：
  • ダイヤモンドおよびケイ素： HF 理論は価電子帯幅を著しく過大評価する（Si およびダイヤモンドで約 25%）。$GoW_0@HF$ 補正はこれらの幅を成功裡に再規格化し、実験的な光電子分光データと優れた一致をもたらす（例：Si の価電子帯幅は HF での 17.00 eV から $GoW_0@HF$ で 13.03 eV に減少し、実験値 12.5 eV と一致）。
  • バンドギャップ： HF がギャップを過大評価する一方、$GoW_0@HF$ は一般的に $GoW_0@LDA$ の結果よりも大きなギャップをもたらす。ダイヤモンドの場合、間接ギャップは 5.75 eV（実験値 5.48 eV に対して）；ケイ素の場合、1.38 eV（実験値 1.17 eV に対して）である。MgO の場合、ギャップは 9.13 eV であり、零点振動再規格化を考慮しても実験値 7.77 eV よりも高い。
• 光学特性（BSE）：
  • BSE-TDA によって計算されたダイ電関数は、ダイヤモンドおよびケイ素について実験的なピーク位置と良好な一致を示す。
  • 酸化物（MgO および $TiO_2$）については、ピーク位置は妥当な精度（実験値から 0.1–0.4 eV 以内）であるが、この手法は系統的に振動子強度（強度）を過大評価する。この過大評価は、DFT 波動関数よりも局在化している HF 波動関数の使用に起因し、それが BSE カーネル内でより強い電子 - 正孔相互作用をもたらすことに帰因される。
  • $TiO_2$ については、実験と比較して計算されたスペクトルにおいて低エネルギー側への小さなエネルギーシフト（0.2–0.4 eV）が観測され、著者らはこれを遮蔽電子 - 正孔引力の小さな $Q$ 極限における近似に起因するとしている。

意義および主張
本論文は、プラズモンポール近似に依存しない、周期系における $GoW_0$ および BSE 計算のための堅牢な密度フィット、ガウス型軌道法を確立する。著者らは以下のことを実証する：

1. HF 出発点の実用性： バンドギャップおよび幅における HF の既知の欠陥にもかかわらず、$GoW_0$ 自己エネルギー補正は Si およびダイヤモンドのバンド構造を再規格化し、実験と良好な一致する結果をもたらす能力があり、バンド幅の点では $GoW_0@LDA$ と同等かそれ以上の性能を示す。
2. 密度フィットの効率性： このアプローチは周期 GW/BSE 計算の計算要求を成功裡に処理し、必要な積分数を削減し、平面波アプローチが効率的でないか、または異なる収束戦略を必要とする可能性のある大きな単位格子を持つ系（$TiO_2$ のような）に手法を適用可能にする。
3. プラズモンポール近似の不在： 有限 $Q$ ベクトルにおいて RPA ハミルトニアンを対角化することにより、この手法は遮蔽相互作用の完全な周波数依存性を捉え、PPA ベースのアプローチに対するより厳密な代替案を提供する。

著者らは、$GoW_0@HF$ アプローチは一般的に実験値（および $GoW_0@LDA$ よりも）よりもわずかに大きなバンドギャップをもたらすものの、ギャップを持つ材料における準粒子励起および光学スペクトルの一貫性のある正確な記述を提供し、固体における多体摂動理論に対するガウス型基底セットおよび HF 出発点の使用を検証していると結論づける。