[論文レビュー] Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

Xinguo Ren, Patrick Rinke+2 more|arXiv (Cornell University)|Dec 15, 2012

Advanced Chemical Physics Studies参考文献 1被引用数 4

ひとこと要約

本稿では、Kohn-Sham枠組み内でRPA、SOSEX、および厳密に取り扱った再正規化単一励起状態（rSE）を組み合わせた、改良された再正規化第二準位摂動理論（rPT2）を導入する。この手法は共有結合、非共有結合、弱い結合系のすべてにおいてバランスの取れた精度を達成し、特に希ガス二量体やS22/S66ベンチマークのような分散支配の相互作用において、以前の近似rSE手法よりも優れた性能を示す。

ABSTRACT

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schrödinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

研究の動機と目的

• 弱い結合系において標準RPAが示す不足結合問題を是正する、体系的に改善された電子相関法の開発。
• 以前のrPT2バージョンで近似されていた再正規化単一励起状態（rSE）寄与の厳密な実装。
• 共有結合エネルギー、反応エネルギー障壁、および分散力を含む非共有結合相互作用を含む多様な化学的環境におけるrPT2のベンチマーク。
• rPT2がバランスの取れた性能と収束性を示すことで、一般に適用可能な電子構造手法としての可能性を確立すること。

提案手法

• 本手法は、RPA（ランダム位相近似）、SOSEX（第二準位スクリーン付き交換）、rSE（再正規化単一励起状態）という3つの無限図式和を組み合わせており、Kohn-Sham基底状態からの第二準位レイノルズシュレーディンガー摂動論に基づく。
• rSE寄与は、フォック行列の準正規化変換を経て、占有および仮想軌道ブロックを別々に対角化することで、エネルギー式における逆行列の評価を効率化する。
• 最終的なrSEエネルギーは、変換された軌道に関する和として表される：$ E_{\text{c}}^{\text{rSE}} = \sum_{ia} \frac{\tilde{f}_{ai}\tilde{f}_{ia}}{\tilde{\epsilon}_{i}-\tilde{\epsilon}_{a}} $、ここで$ \tilde{f}_{ia} $および$ \tilde{\epsilon} $は再解釈された遷移振幅および固有値を表す。
• 本手法は、交換相関関数形（例：PBE）を用いたKohn-Sham基底状態を用い、相関エネルギーを後処理としてSCF計算後に評価する。
• 基底関数の収束性は、Dunningスタイルのcc-pVXZ/aug-cc-pVXZおよびFHI-aims NAOティアーベースセットを用いて評価され、BSSE補正が適用されている。
• 実装の妥当性は、原子化エネルギー、反応エネルギー障壁、非共有結合相互作用セット（S22、S66、希ガス二量体）を用いたベンチマークによって検証されている。

実験結果

リサーチクエスチョン

• RQ1 再正規化単一励起状態（rSE）寄与を厳密に取り扱うことで、希ガス二量体のような弱い結合系の結合エネルギー曲線に見られる物理的に不自然な振るまいが解消されるか？
• RQ2 RPA、SOSEX、rSEを統合的に組み合わせたrPT2フレームワークが、共有結合、遷移状態、分散支配系を含む多様な化学的環境においてバランスの取れた精度を達成するか？
• RQ3 S22およびS66テストセットに関して、rPT2は他のRPAに基づく手法と比較してどのように性能を示すか？
• RQ4 特に拡張関数を含む基底関数の選択が、rPT2の結合エネルギー予測の収束性および精度にどの程度影響を与えるか？

主な発見

• 厳密なrSE取り扱いを施した改良型rPT2は、以前の近似rSEバージョンが示していた希ガス二量体の結合エネルギー曲線における物理的に不自然な振るまいを完全に解消した。
• rPT2は共有結合エネルギー（原子化エネルギー）、反応エネルギー障壁、非共有結合相互作用において、一貫性を保ったバランスの取れた性能を示し、標準RPAおよび他のRPAに基づく手法を上回っている。
• S22およびS66テストセットにおいて、rPT2はRPAに比べて顕著に高い精度を示し、他の先進的なpost-RPA手法と同等またはそれ以上の平均絶対偏差を達成している。
• 特にメタン二量体のような分散支配系の正確な記述には、拡張関数の導入（例：aug-cc-pV5Z、またはハイブリッドティア3/4＋拡張基底（t3/4+a5Z-d））が不可欠である。
• 基底関数収束性に関しては、FHI-aims NAOティアーベースセットの方が標準的なcc-pVXZよりも速いが、同じサイズの基底関数で比較すると、aug-cc-pVXZがより優れた精度を示している。
• rPT2は多様な系において頑健で信頼性が高く、一般に適用可能な電子構造手法としての可能性を裏付けている。