[論文レビュー] Electronic excitation energies of molecular systems from the Bethe-Salpeter equation: Example of the H2 molecule

Elisa Rebolini, Julien Toulouse+1 more|arXiv (Cornell University)|Apr 4, 2013

Advanced Chemical Physics Studies被引用数 5

ひとこと要約

この論文は、H₂分子の最小基底において、静的GW近似におけるBethe-Salpeter方程式を用いて励起エネルギーを計算し、異なる初期1粒子グリーン関数の影響を検証する。正確なグリーン関数を用いることで、一重項励起状態のエネルギー曲線は定性的に正しいが、三重項状態では初期近似が不十分であるため失敗し、静的核の制限により二重励起状態も完全には捉えられていないことが判明した。

ABSTRACT

We review the Bethe-Salpeter equation (BSE) approach to the calculation of electronic excitation energies of molecular systems. We recall the general Green's function many-theory formalism and give the working equations of the BSE approach within the static GW approximation with and without spin adaptation in an orbital basis. We apply the method to the pedagogical example of the H2 molecule in a minimal basis, testing the effects of the choice of the starting one-particle Green's function. Using the non-interacting one-particle Green's function leads to incorrect energy curves for the first singlet and triplet excited states in the dissociation limit. Starting from the exact one-particle Green's function leads to a qualitatively correct energy curve for the first singlet excited state, but still an incorrect energy curve for the triplet excited state. Using the exact one-particle Green's function in the BSE approach within the static GW approximation also leads to a number of additional excitations, all of them being spurious except for one which can be identified as a double excitation corresponding to the second singlet excited state.

研究の動機と目的

• 分子系における電子励起エネルギーを計算するためのBethe-Salpeter方程式（BSE）の静的GW近似の性能を評価すること。
• H₂における励起エネルギー曲線に与える影響を検証するため、初期1粒子グリーン関数として非相互作用HFと正確な解を用いること。
• BSE-GW手法がH₂分子における電荷移動状態および二重励起状態をどの程度正しく記述できるかを検討すること。
• Bethe-Salpeter核が励起エネルギーを補正する役割を果たすメカニズム、特に解離極限における役割を分析すること。

提案手法

• 軌道基底における静的GW近似下でのBSEの作業方程式を導出し、閉殻系のスピン適応型表現を含む。
• BSE計算の初期点として、ハートリー・フォック（非相互作用）および正確な解からの1粒子グリーン関数を用いる。
• 還元可能極効率の逆行列の固値問題を解くことで、Bethe-Salpeter方程式を用いて励起エネルギーを計算する。
• 正確な結果との比較のため、フル配置相互作用（FCI）ベンチマークと不変極効率χIP(ω)の極を用いて精度を評価する。
• 核間距離Rの関数としての励起エネルギー曲線（特に一重項および三重項状態）を分析する。
• 近似された静的核に起因する誤った励起状態を特定し、物理的意味のある状態と区別する。

実験結果

リサーチクエスチョン

• RQ1 静的核を用いたBSE-GW手法は、核間距離Rの全範囲にわたりH₂の一重項および三重項励起状態を正確に記述できるか？
• RQ2 初期1粒子グリーン関数（非相互作用 vs. 正確な解）の選択が、BSE励起エネルギー曲線の精度にどのように影響するか？
• RQ3 周波数依存核が存在しない場合、BSE-GW手法は二重励起状態の特徴を捉えることができるか？
• RQ4 正確なグリーン関数を用いても、解離極限において三重項励起状態を正しく記述できないのはなぜか？
• RQ5 BSE-GW手法に生じる誤った励起状態の起源は何か？また、それらは物理的状態と区別できるか？

主な発見

• 非相互作用1粒子グリーン関数を用いると、解離極限（R ≥ 6.5 bohr）において、最初の一重項および三重項状態の両方で誤った虚数の励起エネルギーが得られる。
• 正確な1粒子グリーン関数を初期点として用いることで、一重項励起状態（1Σg+）のエネルギー曲線は解離極限まで定性的に正しい形状を示す。これは、基本ギャップが初期近似として適切であるためである。
• 同様のアプローチは三重項励起状態（3Σu+）では失敗する。解離極限において、基本ギャップは三重項励起エネルギーの初期近似として不適切である。
• 正確なグリーン関数を用いたBSE-GW手法は複数の誤った三重項励起状態を生成するが、そのうちの1つは二重励起性質を持つ第二一重項励起状態（1Σg+）に対応している。
• 静的Bethe-Salpeter核は、χIP(ω)の初期極と比較して二重励起状態の記述に顕著な改善をもたらさないが、全体的なエネルギー曲線の形状は捉えている。
• 静的核の制限により、平衡幾何近辺でさえも、二重励起状態の正確なエネルギー曲線を記述できないことが判明した。