《The Journal of Physical Chemistry A》：A Benchmark and Basis-Set Extrapolation Study of Hyperfine Coupling Constants from the Random Phase Approximation and σ-Functionals
编辑推荐：
为解决HFCCs（hyperfine coupling constants，超精细耦合常数）计算在准确性与计算成本之间难以兼顾的挑战，德国G?rling课题组的研究人员发展了σ-泛函理论框架，并与RPA（random phase approximation，随机相位近似）、CCSD(T)（coupled-cluster singles doubles and perturbative triples）等方法对比，发现性能最佳的σ-泛函精度可比肩DLPNO–CCSD，且计算需求显著更低。该工作提供了经基组外推的高精度CCSD(T)参考值，为未来方法发展或基准测试奠定了基础。
在化学、物理和材料科学的前沿，理解分子内部电子与原子核之间微妙的磁相互作用至关重要。这种相互作用，具体体现为超精细耦合常数（HFCCs），是连接理论与实验的桥梁，尤其在电子顺磁共振（EPR）光谱学中扮演着核心角色。EPR是研究顺磁性物质（如有机自由基、过渡金属催化剂）不可或缺的工具，而HFCCs和g-张量则是其解读分子自旋密度分布与结构信息的关键参数。然而，精准计算HFCCs长期以来一直是一项艰巨的挑战：高精度的波函数方法如耦合簇（CC）或多组态自洽场（CASSCF）虽然可靠，但其惊人的计算成本将其应用范围死死限制在仅包含几个原子的小体系；而计算成本更可承受的密度泛函理论（DFT），其准确性却往往不尽如人意，且高度依赖于所研究的体系和所选的泛函。尽管近年来出现的基于局域的耦合簇方法，如DLPNO–CCSD（domain-based local pair-natural orbital coupled-cluster singles doubles），为较大体系的CC计算打开了大门，但其本身仍存在计算成本较高、可能引入与正则CCSD的显著偏差以及在几何变化时可能产生数值不稳定性等问题。因此，发展一种既能保持高精度、又具备合理计算效率的理论方法，成为计算化学与光谱学领域一个亟待突破的瓶颈。
在此背景下，随机相位近似（RPA）及其扩展——σ-泛函，作为一种有潜力的理论框架进入了研究者的视野。RPA可以通过多种方式推导，例如基于Kohn-Sham与非相互作用系统间的绝热连接以及涨落-耗散定理，或者从波函数理论的角度将其等价于环状耦合簇双激发。σ-泛函则在RPA的基础上引入了一个经验性的校正项，通过优化一系列标准基准集的反应能、过渡态能、原子化能和非共价相互作用能等参数，显著提升了计算精度，同时其计算成本与标准RPA相当，展现出巨大的吸引力。前期研究已经表明，RPA在主族体系HFCCs的计算上能够超越多种全局、范围分离和局域杂化泛函，达到接近CCSD的精度，而最新的σ-泛函在核磁共振（NMR）屏蔽常数的计算上甚至表现优于RPA。这些发现激励着研究者们探索σ-泛函在HFCCs计算中的潜力。
本项研究由德国埃尔朗根-纽伦堡大学理论化学研究所的Andreas G?rling教授课题组完成，论文发表于《The Journal of Physical Chemistry A》。研究团队旨在系统评估RPA及多种σ-泛函在计算HFCCs方面的表现，并为该领域建立新的高精度基准。为此，他们首先在理论部分完整推导了基于σ-泛函能量表达式的非相对论性HFCCs计算公式。该推导从一般的能量表达式出发，通过分析微分总能量对核磁矩的导数，并利用Handy和Schaefer的Z-矢量（Z-vector）技术绕过显式计算扰动密度矩阵，最终将HFCCs的各向同性（费米接触，FC）和各向异性（自旋偶极，SD）贡献表达为有效驰豫密度矩阵Q与微扰核心哈密顿量的迹。研究表明，将现有RPA实现升级到σ-泛函水平，仅需在计算Z-矢量方程的相关能量导数时，将标准的RPA屏蔽库仑相互作用Wc(iω)替换为修改后的σ-泛函相互作用Wcσ(iω)即可，显示了良好的兼容性与可扩展性。
在计算细节上，研究者选用了由Bartlett课题组引入的两组主族化合物自由基测试集，一组包含26个小分子（少于6个原子），共65个各向同性HFCCs值；另一组包含8个较大自由基，共73个HFCCs值。他们系统地考察了原子轨道基组的影响，重点使用了Jensen课题组开发的专门用于磁性性质计算的pcH和pcJ基组系列。这些基组系列支持对所有HFCCs计算进行系统的完全基组（CBS）外推，且不同的外推方案给出了非常相似的结果。由此获得的经CBS外推的CCSD(T) HFCCs值，可以作为未来方法发展或基准测试的新参考标准。研究中，RPA和σ-泛函计算均采用了分辨率近似（RI），并利用了线性标度的ω-CDGD-RI-RPA算法。作为RPA计算的起点，研究者测试了包括PBE、TPSS、TPSSh、TPSS0、PBE0、B3LYP和ωB97X-V在内的多种密度泛函。对于σ-泛函，则考察了基于PBE、TPSS、PBE0和B3LYP的W1参数化，以及基于PBE和PBE0的S1、S2参数化。所有计算均在其课题组开发的量子化学程序包FermiONs++中完成，其中HFCCs在RPA和σ-泛函水平的解析梯度计算为本工作新实现的功能。
研究结果部分，作者通过多个小节详细展示了他们的发现：
首先，在“小分子测试集的基准结果”中，研究确立了新的CCSD(T)/CBS参考值。与之前使用aug-cc-pVTZ-J基组的结果相比，新参考值变化很小，验证了其可靠性。随后，他们评估了基于不同起点轨道（PBE、TPSS、PBE0、B3LYP、HF）的RPA方法的精度。结果显示，基于HF轨道的RPA表现最佳，平均绝对误差（MAE）和最大绝对误差（MaxAE）最小；而基于DFT轨道的RPA精度稍逊，其中基于PBE0和B3LYP轨道的RPA略优于基于PBE和TPSS轨道的RPA。总体而言，RPA的精度显著优于所测试的常规DFA方法。
其次，在“σ-泛函的基准结果”中，研究比较了多种σ-泛函的表现。其中，基于PBE0的σ-泛函（特别是PBE0-W1和PBE0-S2参数化）达到了所有测试方法中的最高精度，其MAE甚至低于基于HF轨道的RPA，并且与DLPNO–CCSD的精度相当。而基于B3LYP的σ-泛函表现则不尽如人意，误差甚至大于其对应的RPA。研究指出，σ-泛函的表现强烈依赖于其底层参考DFA的性质，成功的参数化需要参考DFA本身已能合理描述所研究的性质。
接着，在“辅助基组对RPA和σ-泛函HFCCs的影响”中，作者考察了用于RI近似的辅助基组大小对计算结果的影响。他们发现，使用从cc-pVDZ-RI到cc-pV5Z-RI的不同辅助基组，对基于HF轨道的RPA以及性能最佳的PBE0-W1 σ-泛函的HFCCs影响非常小。这表明在当前精度要求下，使用较小的辅助基组（如cc-pVDZ-RI）已足以获得收敛的结果，这为大规模计算提供了显著的效率提升。
然后，研究将分析扩展到了“较大主族化合物的测试集”。在这组包含更大、更具挑战性自由基的测试中，PBE0-W1 σ-泛函再次展现了其优越性，取得了所有方法中最小的MAE。基于HF轨道的RPA也表现良好，精度与σ-泛函相近。而基于PBE0轨道的RPA和DLPNO–CCSD的误差则稍大。这一结果进一步证实了性能最佳的σ-泛函在处理更复杂体系时的稳健性和高精度。
最后，在“计算成本考量”中，作者简要比较了不同方法的计算缩放比例。正则CCSD(T)的缩放比例高达O(N7)，而DLPNO–CCSD将其降低到近似O(N)（利用电子结构的局域性），但其计算成本依然可观。原始的RPA缩放比例为O(M6)，但通过RI方法可降至O(M4)，若进一步利用电子局域性则可实现线性标度O(M)。σ-泛函的计算成本与标准RPA相当。因此，性能最佳的σ-泛函在保持与DLPNO–CCSD相近精度的同时，具备了显著更低的计算需求，显示出巨大的应用潜力。
在结论与讨论部分，本研究系统评估了RPA及一系列σ-泛函在计算主族化合物自由基HFCCs方面的表现。研究者成功推导并实现了σ-泛函水平的HFCCs解析梯度计算，为该方法的应用提供了工具。通过使用Jensen课题组开发的pcH和pcJ基组系列，他们首次为测试集提供了系统经CBS外推的高精度CCSD(T)参考值，这将成为未来方法发展与基准测试的宝贵资源。研究表明，基于纯HF轨道的RPA能够显著且稳健地改进DFA的预测精度。而性能最佳的σ-泛函，特别是基于PBE0的W1和S2参数化，其精度甚至超越了RPA，达到了与目前最准确的、可应用于稍大体系的DLPNO–CCSD方法相媲美的水平，同时计算成本显著更低。然而，σ-泛函的成功高度依赖于其底层参考密度泛函近似（DFA）的质量，一个糟糕的参考DFA无法通过σ-泛函的修正获得好的结果。此外，研究还发现使用较小的辅助基组（如cc-pVDZ-RI）对RPA和σ-泛函的HFCCs结果影响甚微，这为大规模计算提供了重要的效率优化策略。综上所述，本工作不仅为HFCCs的计算建立了新的基准，也凸显了σ-泛函作为一种高精度、高效率的理论方法，在计算磁性性质方面的巨大潜力与应用前景，为未来在更大体系乃至周期性体系中应用高级电子相关方法计算EPR参数铺平了道路。
