[论文解读] 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中的局限性。
• 确保激发态几何结构、锥形交叉点以及垂直/激发态能量的精确计算。
• 通过高效整合虚数位移项到现有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秒 vs. 3387秒），主要源于虚数位移项评估的 $O(N_{\text{act}}^9)$ 标度。
• 在FeP中，虚数位移项 $\mathbf{d}^{(2)}_{\text{shift}}$ 占了相关密度矩阵计算时间的52%（达383秒），而在小体系中可忽略不计（<1秒）。
• 求解 $\lambda$-方程的额外开销仅增加11秒（FeP），表明对迭代求解器影响极小。
• 该方法已成功实现并集成到公开的BAGEL量子化学程序包中，推动了多参考量子化学应用的广泛应用。