经济二阶方法精确计算芯能级电离能

《Journal of Chemical Theory and Computation》：Accurate Core-Level Ionization Energies from an Affordable Second-Order Approach

　　

研究人员提出了一种利用二阶方法精确计算芯能级电离势（IPs）的方法。所评估的理论框架基于迭代二阶代数图解构造（ADC(2)）和组态相互作用单激发加微扰二阶修正（CIS(D)）方法。在此，研究人员将IP-ADC(2)和IP-CIS(D)的高效实现 [Mester, D.; Kállay, M. J. Chem. Theory Comput. 2023, 19, 3982–3995] 与芯-价分离（CVS）近似相结合，从而能够从芯区进行电离。这些方法表现出非常有利的标度行为：计算成本实际上是三次方的，仅有一个依赖于活性芯轨道数量的温和前置因子。此外，所得的基于波函数的方法与自旋标度技术相结合，并成功扩展到双杂化（DH）泛函，而无需对二阶修正的实现进行任何修改。通过基准计算对所提出方法的性能进行了全面评估。研究人员的成果表明，对双激发的迭代处理至关重要，凸显了更先进的DH ansatz的必要性。此外，SOS0-PBE0-2/CVS-IP-ADC(2)方法与更昂贵的高阶耦合簇方法具有很强的竞争力。最后，二阶修正仅在总体计算时间中引入可忽略的开销——对于使用三-zeta基组的61个原子的氮杂富勒烯分子，仅需约1分钟——从而使得能够对扩展分子系统进行精确计算。

**论文解读：经济二阶方法精确计算芯能级电离能**

**研究背景与问题**
X射线光电子能谱（XPS）是材料科学和化学中的核心实验技术，能够通过芯轨道的局域化和元素特异性提供分子与电子结构的详细信息。然而，准确且高效的理论描述芯区电离态仍是量子化学中的重大挑战。现有方法中，基于波函数的方法（如耦合簇CC）精度高但计算成本高昂，随系统大小呈陡峭标度；密度泛函理论（DFT）方法成本低但精度存疑且缺乏系统改进路径。因此，亟需发展一种兼具精度与效率的方法。本文旨在提出基于二阶方法的芯能级电离能（IPs）计算方案，通过结合芯-价分离（CVS）近似、自旋标度技术和双杂化（DH）泛函，在保证精度的同时大幅降低计算成本。论文发表在《Journal of Chemical Theory and Computation》。

**关键技术与方法**
本文采用以下主要技术：
1. **二阶电离势方法**：基于IP-CIS(D)（组态相互作用单激发加微扰二阶修正）和IP-ADC(2)（迭代二阶代数图解构造）框架，直接计算电离能。
2. **芯-价分离（CVS）近似**：忽略芯激发与价激发之间的耦合，将问题限制在芯激发子空间，实现高效求解。
3. **密度拟合（DF）近似**：将四中心双电子积分因子化为三中心积分，避免存储和变换大型阵列。
4. **自旋标度技术**：对二阶项的同自旋和反自旋分量分别进行标度（如SOS方案中反自旋因子1.3）。
5. **双杂化（DH）泛函**：将DFT参考轨道与标度的二阶关联修正结合，具体采用SOS0-PBE0-2泛函。
计算中采用cc-pCVXZ-X2C基组和自旋自由X2C单电子哈密顿量（SFX2C-1e）处理标量相对论效应。基准测试使用两组文献数据集：小型分子集（14个分子，23个电离）和CORE65集（32个分子，65个芯结合能）。所有计算在Mrcc程序包中实现。

**研究结果**
**1. Small-Molecule Compilation（小分子编译）**
在小型分子基准集上，比较CVS-EOMIP-CCSD和CCSDT等高阶方法。结果表明：
- CCSDT方法表现最佳（平均绝对误差MAE=0.14 eV，标准差SD=0.17 eV），但成本高。
- ADC(2)方法平均误差（ME）小（约0.06 eV），但SD大（1.58 eV），精度不足。
- SOS-ADC(2)降低SD但引入显著蓝移（过估电离能），导致ME和MAE增大。
- 双杂化SOS0-PBE0-2/ADC(2)方法显著改进：MAE=0.68 eV，SD=0.59 eV，精度接近CCSD（MAE=1.82 eV），且优于CCSD。
- CIS(D)类方法误差更大，不具竞争力。

**2. CORE65 Compilation（CORE65编译）**
在更大、更具应用相关性的CORE65集上，趋势类似：
- CCSD*（含三重微扰）表现最佳（MAE=0.43 eV，SD=0.22 eV）。
- CCSD过估严重（MAE>2 eV）。
- ADC(2)和SOS-ADC(2)性能与小型集一致。
- SOS0-PBE0-2/ADC(2)与CCSD*相比精度相当（MAE=1.02 eV，SD=0.50 eV），且优于CCSD。
- 峰间距分析表明，ADC(2)类方法在相对误差上优于CCSD和CCSD*，凸显实际应用优势。

**3. Performance by Element（按元素性能）**
按C、N、O、F元素分别分析：
- 对所有元素，SOS0-PBE0-2/ADC(2)的精度均高于CCSD，且随电离能增加优势更明显（如F的ME差异达2 eV）。
- 两类方法的精度（SD）几乎相同（0.3–0.5 eV），但DH方法最大误差（MAX）更小。

**4. Timings（计时）**
对61原子氮杂富勒烯分子（N K-边电离）使用cc-pCVTZ-X2C基组：
- HF/KS参考计算约500分钟，分子积分约88分钟。
- 二阶修正本身仅需0.5–1.3分钟，占总时间可忽略。
- 迭代（ADC(2)）与微扰（CIS(D)）方法计算成本几乎相同。

**总结与结论**
讨论部分强调：迭代处理双激发至关重要，CIS(D)类方法不推荐使用；双杂化泛函有效补偿了SOS方案引入的蓝移，使SOS0-PBE0-2/ADC(2)达到与CCSD相当甚至更高的精度，同时计算成本极低。结论部分如下（翻译）：
在本研究中，实现了使用二阶方法计算芯结合能。为此，研究人员提出了一个理论框架，将微扰IP-CIS(D)和迭代IP-ADC(2)方法与CVS近似相结合，从而能够从芯区进行电离。提供了高效实现所需的工作方程。算法分析表明，这些方法表现出非常有利的标度行为。修正的计算成本实际上是三次方的，仅有一个依赖于活性芯轨道数量的温和前置因子。所得的CVS-IP-CIS(D)和CVS-IP-ADC(2)方法进一步扩展到DH泛函。由于XC贡献仅出现在参考轨道集的确定中，因此无需对实现进行修改；只需将二阶修正项按泛函的关联混合因子进行标度。使用两个广泛采用的基准集（一个包含小分子，另一个为著名的CORE65编译），以实验电离能为参考，评估了所提出方法的性能。总之，结果表明，CVS-IP-ADC(2)通常提供比CCSD更低的ME和MAE值，但其精度由于较大的SD和误差范围而有所限制。SOS-ADC(2)方法显著改善了这一点；然而，电离能中观察到的显著蓝移导致ME和MAE值增大，尽管SD和范围更小。通过应用DH泛函，这种缺陷得到了有效修正。SOS0-PBE0-2/ADC(2)方法提供的精度明显超过昂贵的CCSD方法，同时保持几乎相同的精度。该研究还证实了双激发迭代处理的重要性；因此，不推荐使用基于CIS(D)的方法。这一结论进一步得到如下观察的支持：微扰和迭代二阶方法的总体计算成本几乎相同。基于对使用cc-pCVTZ-X2C基组的61原子氮杂富勒烯分子进行的计算，研究人员证明了在获得参考轨道集和所需的分子轨道积分后，二阶修正的评估仅需约一分钟。此外，内存需求极低，因为高阶方法所需的分子积分列表的虚-虚块在此情形中并不需要。此外，在此形式下完全避免了四指标数组的计算和存储，且所需的所有中间体都不包含超过一个虚指标。