[论文解读] Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

Xinguo Ren, Patrick Rinke+2 more|arXiv (Cornell University)|Dec 15, 2012

Advanced Chemical Physics Studies参考文献 1被引用 4

一句话总结

该论文提出了一种改进的重整化二阶微扰理论（rPT2），用于电子相关能计算，结合了RPA、SOSEX以及在Kohn-Sham框架下严格处理的重整化单激发（rSE）。该方法在共价键、非共价键及弱键合体系中均实现了平衡的精度，相较于以往近似处理的rSE方法，性能显著提升，尤其在以色散作用主导的相互作用（如稀有气体二聚体及S22/S66基准测试）中表现更优。

ABSTRACT

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schrödinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

研究动机与目标

• 开发一种系统性改进的电子相关方法，以纠正标准RPA在弱键合体系中出现的过度解离问题。
• 严格实现重整化单激发（rSE）贡献，此前在rPT2的早期版本中曾被近似处理。
• 在多种化学环境中对rPT2进行基准测试，包括共价键能、反应能垒以及非共价相互作用（如色散力）。
• 通过展示其平衡的性能和收敛行为，确立rPT2作为一种有前景且普遍适用的电子结构方法。

提出的方法

• 该方法结合了三种无限图求和：RPA（随机相 approximation）、SOSEX（二阶屏蔽交换）和rSE（重整化单激发），均源自Kohn-Sham参考态的二阶Rayleigh-Schrödinger微扰理论。
• 通过Fock矩阵的半正则变换计算rSE贡献，分别对占据轨道和虚拟轨道块进行对角化，以实现能量表达式中逆矩阵的高效求解。
• 最终的rSE能量表示为变换后轨道的求和形式：$ E_{\text{c}}^{\text{rSE}} = \sum_{ia} \frac{\tilde{f}_{ai}\tilde{f}_{ia}}{\tilde{\epsilon}_{i}-\tilde{\epsilon}_{a}} $，其中 $ \tilde{f}_{ia} $ 和 $ \tilde{\epsilon} $ 分别为重新解释的跃迁振幅和本征值。
• 该方法采用Kohn-Sham参考态，结合交换-相关泛函（如PBE），并将相关能作为后-SCF修正进行计算。
• 通过Dunning型cc-pVXZ/aug-cc-pVXZ和FHI-aims NAO层级基组评估基组收敛性，并应用BSSE校正。
• 通过原子化能、反应能垒以及非共价相互作用集（S22、S66、稀有气体二聚体）的基准测试验证了实现的正确性。

实验结果

研究问题

• RQ1 对重整化单激发（rSE）贡献进行严格处理，是否能消除弱键合体系（如稀有气体二聚体）中结合能曲线的非物理行为？
• RQ2 在统一的rPT2框架中同时包含RPA、SOSEX和rSE，是否能在包括共价键、过渡态及色散主导体系在内的多种化学环境中实现平衡的精度？
• RQ3 与其它基于RPA的方法相比，rPT2在非共价相互作用能计算中表现如何，特别是在S22和S66测试集中？
• RQ4 基组的选择，尤其是弥散函数的引入，对rPT2结合能预测的收敛性和精度有多大影响？

主要发现

• 经过严格rSE处理的改进rPT2方法，消除了早期近似rSE版本中在稀有气体二聚体结合能曲线中出现的非物理行为。
• rPT2在共价键能（原子化能）、反应能垒和非共价相互作用方面均表现出良好的平衡性能，其一致性优于标准RPA及其他基于RPA的方法。
• 在S22和S66测试集中，rPT2相较于RPA显著提升了精度，平均绝对偏差与其它先进后-RPA方法相当或更优。
• 为准确描述甲烷二聚体等色散主导体系，必须引入弥散函数，特别是使用aug-cc-pV5Z或混合层级3/4 + 弥散基组（t3/4+a5Z-d）时效果最佳。
• 与标准cc-pVXZ基组相比，FHI-aims NAO层级基组在rPT2中收敛速度更快，但相同基组尺寸下，aug-cc-pVXZ提供的精度更优。
• rPT2在各类体系中表现出强健性和可靠性，证实其作为普遍适用电子结构方法的巨大潜力。