[論文レビュー] Scalable and predictive spectra of correlated molecules with moment truncated iterated perturbation theory

Oliver J. Backhouse, Alejandro Santana‐Bonilla+1 more|arXiv (Cornell University)|Jun 7, 2021

Spectroscopy and Quantum Chemical Studies参考文献 94被引用数 19

ひとこと要約

本論文は、相関分子の全単粒子スペクトル関数を計算するためのスケーラブルで予測可能かつ自己整合的な多体摂動法—モーメント截断反復摂動理論（AGF2）—を導入する。スペクトルモーメント保存による有効自己エネルギー動的特性の圧縮により、本手法はO[N⁵]スケーリングを達成し、大規模並列計算（最大約1000軌道）を可能にするとともに、アーテミシンリンのような複雑系の電子励起スペクトルを正確に解明し、長年のその最低励起エネルギー状態の割り当てに関する矛盾を解消する。

ABSTRACT

A reliable and efficient computation of the entire single-particle spectrum of correlated molecules is an outstanding challenge in the field of quantum chemistry, with standard density functional theory approaches often giving an inadequate description of excitation energies and gaps. In this work, we expand upon a recently-introduced approach which relies on a fully self-consistent many-body perturbation theory, coupled to a non-perturbative truncation of the effective dynamics at each step. We show that this yields a low-scaling and accurate method across a diverse benchmark test set, capable of treating moderate levels of strong correlation effects, and detail an efficient implementation for applications up to $\sim1000$ orbitals on parallel resources. We then use this method to characterise the spectral properties of the artemisinin anti-malarial drug molecule, resolving discrepancies in previous works concerning the active sites of the lowest energy fundamental excitations of the system.

研究の動機と目的

• 強相関分子の全単粒子グリーン関数スペクトルを計算するためのスケーラブルで高精度な手法の開発。
• 標準的DFTおよび摂動的アプローチが相関系における励起エネルギーとバンドギャップを記述するうえで抱える限界の解消。
• 複雑な分子系のための効率的並列実装を備えた大規模シミュレーション（最大約1000軌道）の実現。
• 抗マラリア薬アーテミシンリン分子における最低エネルギー励起状態の空間的局在化に関する文献上の矛盾の解消。

提案手法

• 本手法は完全に自己整合的な補助グリーン関数アプローチ（AGF2）を採用し、補助自由度を介して2次自己エネルギーを静的固有値問題に再定式化する。
• 第一スペクトルモーメントを完全に保存するモーメント截断近似を用いることで、計算コストを低減しつつも精度を確保する。
• 有効自己エネルギーは極の和として表現される：Σpq(ω) = Σk (vpk v†qk)/(ω − Ek)，これによりダイソン方程式が行列固有値問題として定式化される。
• アルゴリズムはオープンソースのPySCFフレームワークに実装されており、効率的な並列化と大規模応用を可能にする。
• 粒子およびホール励起を分離する必要がなく、あらゆるエネルギースケールにわたる全励起スペクトルへの同時アクセスが可能である。
• 本手法は、動的自己エネルギー成分の体系的截断を施した自己整合的拡張版ADC(2)と同等であり、動的自己エネルギー成分の系統的截断により、スペクトル精度が向上している。

実験結果

リサーチクエスチョン

• RQ1 モーメント截断と自己整合性を備えた多体摂動法は、相関分子スペクトルにおいて高い精度と低い計算スケーリングを両立できるか？
• RQ2 多様な分子系において、モーメント截断AGF2手法のスペクトル精度は、ADC(2)、GW、EOM-CCSDといった既存手法と比較してどの程度か？
• RQ3 アーテミシンリンのような複雑分子における励起状態の特性割り当てに関する曖昧さを、この手法はどの程度解消できるか？
• RQ4 並列アーキテクチャ上でのこの手法を用いて、信頼性のあるスペクトル精度を維持できる最大の系サイズ（軌道数）はどの程度か？
• RQ5 相関系において、完全な動的自己エネルギー記述と比較して、モーメント截断戦略はスペクトル分解能を向上させるか？

主な発見

• AGF2手法はO[N⁵]の計算スケーリングを達成し、効率的な並列処理を可能にし、最大約1000軌道の系のシミュレーションが可能である。
• GW100ベンチマークセットにおいて、AGF2は電荷励起に対してADC(2)やGWと同等の精度を示し、EOM-CCSDよりもスケーラビリティに優れている。
• 本手法は、アーテミシンリンにおける最低エネルギー励起状態に関する文献上の矛盾した割り当てを明確に解消し、活性部位がエンドパーオキサイドブリッジであることを特定した。これは実験的Photoelectronスペクトロスコピーと整合的である。
• モーメント截断自己エネルギー近似は、厳密な保存近似を破るものの、重要なスペクトル特徴を保持し、非截断自己整合的アプローチよりも精度が向上している。
• PySCFへの実装により、再現性があり大規模な計算が可能となり、エネルギー依存の全スペクトル関数への直接アクセスが可能である。
• 完全な自己整合性のおかげで、基底状態の選択に依存しないため、強相関領域においても高いロバスト性を示している。