降成本四分量相对论双电离势运动方程耦合簇方法：包含4空穴-2粒子激发和三体簇

《Journal of Chemical Theory and Computation》：Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

　　

双电离势（DIP）运动方程（EOM）耦合簇（CC）方法，在包含单、双和三激发（CCSDT）计算的基础上加入4空穴-2粒子（4h-2p）激发（缩写为DIP-EOMCCSDT(4h-2p)），连同其微扰近似DIP-EOMCCSD(T)(a)(4h-2p)，被扩展到相对论四分量（4c）框架。此外，研究人员引入并测试了一种新的计算上实用的DIP-EOMCC方法，称为DIP-EOMCCSD(T)(?)(4h-2p)，它近似处理了DIP-EOMCCSD(T)(a)(4h-2p)方法中的4h-2p关联，并将DIP-EOMCCSDT(4h-2p)和DIP-EOMCCSD(T)(a)(4h-2p)的特征性??⁸N⁸标度（scaling）降低为系统大小??的??⁷N⁷。通过使用冻结核自然旋量（frozen natural spinor, FNS）近似，根据明确定义的基于占据数的阈值减少进入DIP-EOMCC计算相关步骤的未占据旋量数目，进一步提高了计算效率。所得的4c-FNS-DIP-EOMCC方法被用于计算从氩到氡的惰性气体原子系列的双电离势，以及Cl₂、Br₂、HBr和HI的垂直双电离势，这些系统已在实验上被研究过。研究人员证明，当使用完全基组外推和10^-4.5的FNS截断阈值时，4c-FNS-DIP-EOMCCSD(T)(?)(4h-2p)计算能够预测与实验数据一致的双电离势，优于其非相对论和无自旋标量相对论对应方法，特别是在检查具有更强自旋-轨道耦合效应的双电离势时。

本文解读基于论文《Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters》，发表于《Journal of Chemical Theory and Computation》。

**研究背景与问题**：准确处理相对论效应是现代计算化学的重要议题，双电离势（DIP）的预测对理解光电子和俄歇电子能谱至关重要。自旋-轨道耦合等相对论贡献显著影响含重原子体系的双电离谱。现有标量相对论方法（如零阶正则近似ZORA、Douglas-Kroll-Hess变换DKH、无自旋精确二分量SFX2C）虽广泛应用，但更完整的四分量（4c）方法因使用旋量基组、复值哈密顿量及缺乏自旋对称性而导致计算成本极高，难以推广。此前4c-DIP-EOMCC方法局限于DIP-EOMCCSD(3h-1p)水平，无法充分处理4h-2p激发和三体簇关联。因此，研究人员旨在开发计算上实用的4c-DIP-EOMCC方法，使其能准确预测重元素体系的DIP，并适用于大基组（如四重ζ QZ质量）。

**研究内容与结论**：研究人员将原有高水平DIP-EOMCCSDT(4h-2p)及其微扰近似DIP-EOMCCSD(T)(a)(4h-2p)扩展到四分量相对论框架，并引入一种新的低成本近似方法DIP-EOMCCSD(T)(?)(4h-2p)，将计算标度从??⁸降至??⁷。结合冻结核自然旋量（FNS）技术进一步降低虚拟旋量数目。将所发展的4c-FNS-DIP-EOMCC方法应用于惰性气体原子（Ar至Rn）及分子（Cl₂、Br₂、HBr、HI），通过与实验数据对比，验证了方法的准确性。研究证实，结合完全基组外推（CBS）和FNS阈值10^-4.5时，4c-FNS-DIP-EOMCCSD(T)(?)(4h-2p)能预测与实验高度一致的DIP，特别是对自旋-轨道耦合效应较强的体系，其表现优于非相对论（NR）和无自旋标量相对论方法。

**主要关键技术方法**（不超过250字）：采用四分量狄拉克-库仑（DC）哈密顿量描述相对论效应，基组使用未收缩的dyall.avnz（n=2,3,4）系列。核心方法为DIP-EOMCCSD(T)(?)(4h-2p)：该法通过微扰方式近似处理三体簇（???₃）和4h-2p激发，将CCSDT步骤简化为CCSD后非迭代修正，使总标度降至??⁷。进一步结合冻结核自然旋量（FNS）近似：基于MP2单粒子密度矩阵的虚拟-虚拟块对角化，按占据数阈值（本文取10^-4.5）截断低贡献虚拟旋量，减少后续相关步骤的计算量。所有计算采用自主软件BAGH实现。样本为惰性气体原子及实验几何下的双原子分子（Cl₂、Br₂、HBr、HI），化学芯轨道在相关步骤中冻结。

**研究结果**（保留各小标题）：

**4.1 选择FNS阈值**：通过Cl₂分子在4c-FNS-DIP-EOMCCSD(3h-1p)水平下的收敛性测试，发现FNS阈值10^-4.5时DIP截断误差仅0.01 eV，且仅需约40%的虚拟空间即可达到与规范基组一致的结果，验证了该阈值在精度与成本间的合理折中。

**4.2 基组收敛性**：对Cl₂使用dyall.avnz（n=2-4）基组，DIP值随基组增大而增加，CBS外推后进一步增大0.1-0.2 eV；加入弥散函数影响可忽略。表明必须采用大基组或CBS外推才能与实验可靠比较。

**4.3 包含高激发的效应**：比较不同方法在Cl₂上的表现，发现DIP-EOMCCSD(3h-1p)/dyall.av4z高估DIP约0.46-0.66 eV，而加入4h-2p激发（CCSD(4h-2p)）导致低估约0.23 eV，表明N电子参考态与(N-2)电子目标态关联不平衡。引入三体簇修正的DIP-EOMCCSD(T)(?)(4h-2p)恢复平衡，结果与昂贵DIP-EOMCCSDT(4h-2p)和DIP-EOMCCSD(T)(a)(4h-2p)一致（偏差<0.03 eV），证实新方法能准确再现高级别结果。

**4.4 基准测试：原子与分子**：对惰性气体原子（Ar-Rn），DIP-EOMCCSD(T)(?)(4h-2p)/CBS的统计误差（MAD=0.11 eV, RMSD=0.04 eV）远优于DIP-EOMCCSD(3h-1p)/CBS（MAD=0.47 eV）。对分子（Cl₂、Br₂、HBr、HI），新方法的MAD=0.30 eV, RMSD=0.15 eV，同样优于低级别方法。特别对Br₂，新方法能准确分辨自旋-轨道耦合导致的A 0_g与A 1_g态分裂（差值0.08 eV），与实验及多参考组态相互作用（MRCI）结果一致。

**4.5 相对论处理的影响**：比较NR、SFX2C1e和4c-DC哈密顿量下的DIP。对原子，NR的MAD=0.77 eV，SFX2C1e为0.67 eV，而4c-DC降为0.11 eV；对分子，NR的MAD=0.34 eV，SFX2C1e为0.30 eV，4c-DC为0.30 eV。表明只有完整的4c-DC哈密顿量能可靠描述重原子DIP；Gaunt和Breit修正影响可忽略。

**总结与结论**：研究成功开发了一系列包含4h-2p激发和三体簇的四分量相对论DIP-EOMCC方法，其中新引入的DIP-EOMCCSD(T)(?)(4h-2p)结合FNS技术，在保持精度的同时大幅降低计算成本。基准测试表明，经CBS外推后，该方法预测的DIP与实验吻合良好，且对重元素体系的自旋-轨道耦合效应描述至关重要。结论翻译如下：“在这项工作中，研究人员开发和测试了一套包含高达4h-2p激发和三体簇的四分量相对论DIP-EOMCC方法。将原始DIP-EOMCCSDT(4h-2p)及其微扰近似扩展到完全相对论领域，并引入了新的低成本近似DIP-EOMCCSD(T)(?)(4h-2p)，能够以低得多的计算成本准确再现DIP-EOMCCSD(T)(a)(4h-2p)和DIP-EOMCCSDT(4h-2p)的结果，使该方法成为研究原子和分子系统DIP的实用工具。结合FNS截断方案，显著降低了计算成本而不损失精度。基准测试表明，经CBS外推后，该方法与实验吻合良好。DIP值对基组大小敏感，需大基组或CBS外推。仅四分量哈密顿量能提供稳健准确的DIP。”