[論文レビュー] Imaginary shift in CASPT2 nuclear gradient and derivative coupling theory

Jae Woo Park, Rachael Al-Saadon+2 more|arXiv (Cornell University)|Apr 14, 2019

Spectroscopy and Quantum Chemical Studies被引用数 4

ひとこと要約

本稿では、虚数シフト正則化を用いたCASPT2における解析的核勾配理論を提示する。従来の実数シフト実装を拡張し、励起状態の性質および円錐状交差計算における精度を向上させつつ、計算コストの増加は最小限に抑えられる。p-HBDIおよびFePモデルを用いた検証でその有効性が示された。

ABSTRACT

We report the analytical nuclear gradient theory for complete active space second-order perturbation theory (CASPT2) with imaginary shift, which is commonly used to avoid divergence of the perturbation expression. Our formulation is based on the Lagrangian approach and is an extension of the algorithm for CASPT2 nuclear gradients with real shift. The working equations are derived and implemented into an efficient parallel program. Numerical examples are presented for the ground- and excited-state geometries and conical intersections of a green fluorescent protein model chromophore, $p$-HBDI$^-$. We also report timing benchmarks with adenine, $p$-HBDI$^-$, and iron porphyrin. It is demonstrated that the energies and geometries obtained with the imaginary shift improve accuracy at a minor additional cost which is mainly associated with evaluating the effective density matrix elements for the imaginary shift term.

研究の動機と目的

• イントラーダー状態の発散を回避するため、虚数シフト正則化を用いたCASPT2における解析的核勾配理論の開発。
• 実数シフトの限界を克服し、より安定で特異点のない正則化スキームを実装すること。
• 励起状態の幾何構造、円錐状交差、垂直励起状態エネルギーの正確な計算を保証すること。
• 虚数シフト項を既存のCASPT2勾配アルゴリズムに効率的に統合することで、計算コストの増加を最小限に抑えること。
• p-HBDIおよび鉄プロフィリンを含むベンチマーク系を用いて、異なる活性空間サイズでの手法の妥当性を検証すること。

提案手法

• ラグランジュ法を応用し、虚数シフトを含むCASPT2における解析的核勾配の導出を実施。従来の実数シフト形式を拡張する。
• 虚数シフト項のための作業式を導出。$\mathbf{d}^{(2)}_{\text{shift}}$および$\tilde{y}_{I,M}^{\text{shift}}$を含み、エネルギー分母における正則化を記述する。
• 並列計算を活用し、BAGEL量子化学パッケージ内に理論を実装。
• 既存コードの変更を最小限に抑えつつ、$\lambda$-方程式および$Z$-ベクトル方程式の解法に虚数シフト項を統合。
• 特異点のない挙動を保証するため、複素数シフト$1/\Delta \to \Re[1/(\Delta + i\epsilon)]$を用いる。特に円錐状交差付近で重要である。
• 新しい項と従来のCASPT2寄与項を組み合わせることで計算コストを最適化し、大規模な活性空間ではコスト増加を約15%に抑える。

実験結果

リサーチクエスチョン

• RQ1 CASPT2核勾配における虚数シフトは、特に円錐状交差付近で実数シフトと比較して、精度をどのように向上させるか？
• RQ2 CASPT2核勾配計算に虚数シフト項を含めることによる計算コストの増加はどの程度か？
• RQ3 活性空間サイズの増大に伴い、虚数シフト手法の性能はどのようにスケーリングするか？
• RQ4 虚数シフト法は、実数シフトと比較して、シフトパラメータ$\epsilon$への感度をどの程度低減するか？
• RQ5 虚数シフト形式は、励起状態の幾何最適化および円錐状交差計算において、数値的安定性と精度を維持できるか？

主な発見

• 虚数シフト形式は、$\epsilon$の変動に対して実数シフトと比較して感受性が低く、励起状態計算におけるロバスト性が向上する。
• 小規模な活性空間（例：アデニン、p-HBDI）では、虚数シフト勾配の計算コストは実数シフトとほぼ同一であり、追加時間は2–3%にとどまる。
• 大規模な活性空間（例：FeP、CAS(10e,9o)）では、ウォールタイムが約15%増加（2947 s 対 3387 s）したが、主に虚数シフト項評価の$O(N_{\text{act}}^9)$スケーリングに起因する。
• FePでは、虚数シフト項$\mathbf{d}^{(2)}_{\text{shift}}$が相関密度行列計算時間の52%を占め（383秒）、小規模系では無視できる程度（<1秒）であった。
• $\lambda$-方程式を解くための追加コストはFePでたったの11秒増加にとどまり、反復解法への影響は最小限であった。
• 本手法は、公開されたBAGEL量子化学パッケージに正常に実装・統合され、多電子状態量子化学応用分野における広範な利用を可能にした。