Random phase approximation-based local natural orbital coupled cluster theory

本文在基于局域自然轨道的耦合簇（LNO-CC）框架内，将随机相位近似（RPA）作为二阶 Møller-Plesset 微扰理论（MP2）的一种稳健替代方案加以介绍，表明基于 RPA 的 LNO-CC 在保持对具有较大能隙体系精度的同时，为金属体系提供了显著更快的收敛速度。

要点

想象一下，你正在试图计算一台巨大而复杂的机器的总能量，比如一台城市规模的引擎。为了得到完全准确的答案，你需要追踪每一个运动部件，以及每个部件如何与其他每一个部件相互作用。在化学世界中，这台“机器”是一个分子或晶体，而“部件”则是电子。

对大型系统进行如此完美的计算，就像在涨潮时试图数清沙滩上的每一粒沙子——它需要如此庞大的计算能力，以至于实际上是不可能的。

问题：“足够好”的捷径
为了解决这个问题，科学家们使用了一种称为片段嵌入的技巧。他们将这台大机器分解为更小、更易处理的块（片段）。

1. 高精度区域：他们在片段的中心计算最重要的相互作用，使用极端且昂贵的高精度方法。
2. “低层级”区域：对于片段中远离中心的部分，他们使用一种“低层级”理论——一种更快、更便宜但精度较低的方法——来估算那些遥远部分的行为。

几十年来，标准的“低层级”方法被称为MP2。这就像用粗略的草图来估算背景景色。它对大多数情况都有效，但存在两个主要缺陷：

• 粘合问题：它经常高估非粘性物质（如两个分离的分子）之间的结合强度。
• 金属问题：当应用于金属（其中电子像河流一样自由流动）时，MP2 完全失效，并给出荒谬的、无限大的答案。

新解决方案：RPA 和 SOSEX
本文介绍了两种新的“低层级”方法来取代 MP2：RPA（随机相位近似）和SOSEX（二阶屏蔽交换）。

将 MP2 想象成用钝铅笔绘制的草图。它很快，但线条粗糙且有时错误。

• RPA 就像是用更细的笔绘制的草图，它理解了“电风”（屏蔽）如何平滑相互作用。它更好地处理了“粘合问题”，并且关键的是，在观察金属时不会失效。
• SOSEX 是 RPA 的更精细版本，它修正了 RPA 有时会犯的一种特定错误（自相互作用）。

作者做了什么
研究人员构建了他们计算引擎（称为LNO-CC）的新版本，该版本可以将旧的 MP2“草图”替换为这些新的 RPA 和 SOSEX 草图。他们在三类挑战上测试了这个新引擎：

1. 非粘性分子：由弱力结合的分子系统。
2. 化学反应：计算反应发生所需跨越的能量“山丘”。
3. 金属：锂和铜的块状样本。

结果

• 对于非粘性分子：新的 RPA/SOSEX 方法的表现与旧的 MP2 方法一样好。它们没有使情况变糟；它们的精度同样高。
• 对于金属：这正是新方法大放异彩的地方。虽然 MP2 难以给出金属的良好答案，但 RPA 尤其是 SOSEX 提供了更快、更准确的结果。它们以远少的计算努力就达到了“完美”的答案。
• “速度”因素：作者发现，使用 RPA 和 SOSEX 作为背景“草图”，使得计算的高精度部分能够更快地收敛（确定最终答案）。这就像拥有一张更好的背景景色地图，让你能够专注于前景的细节而不会迷失方向。

结论
本文证明，RPA 和 SOSEX 是这些复杂计算中替代旧 MP2 方法的优秀现代方案。它们对于标准分子同样有效，但对于金属以及加速整个计算过程则显著更优。它们提供了一种更可靠的方式来模拟量子世界，而无需使用城市规模的超级计算机。

技术摘要

技术摘要：基于随机相位近似（RPA）的局域自然轨道耦合簇理论

问题陈述
片段嵌入和局域关联方法（如局域自然轨道耦合簇，LNO-CC）依赖于低阶理论来定义局域嵌入子空间，并处理活性区域之外的长程静电和关联效应。二阶 Møller-Plesset 微扰理论（MP2）是这些框架中的标准低阶理论。然而，MP2 在特定体系中表现出显著局限性：它在非共价结合分子复合物和分子固体中定量地高估相互作用能，并且在体相金属的热力学极限（TDL）下，由于能隙消失而在形式上发散。这些缺陷损害了局域关联方法在这些关键体系中的准确性和收敛性。

方法论
作者提出在 LNO-CC 框架内用随机相位近似（RPA）及其密切相关的二阶屏蔽交换（SOSEX）近似替代 MP2。该实现利用了 RPA 的直接环耦合簇双激发（drCCD）公式，使得 RPA 能够无缝集成到现有的基于 MP2 的嵌入结构中。

关键方法论组件包括：

1. 基于 RPA 的 LNO 构建：作者不使用 MP2 振幅构建用于生成局域自然轨道（LNO）的单粒子密度矩阵，而是利用 drCCD 振幅（$T_{iajb}$）。这些振幅通过求解非线性 drCCD 方程获得，该方程仅保留直接粒子 - 空穴环图。
2. 高效实现：作者采用了一种因子化的 drCCD 算法（基于 Heßelmann 和 Kállay 的工作），将分子的计算成本降低至 $O(N^4)$，将周期性固体的计算成本降低至 $O(N_k^2 n^4)$，这在渐近意义上优于传统的密度拟合 MP2（周期性体系为 $O(N^5)$ 标度）。该方法利用轨道能量分母的 Cholesky 分解和电子排斥积分的密度拟合。
3. 外部能量修正：该框架保留了标准的 LNO-CC 结构，即总关联能是活性空间内的局域 CCSD 能量与外部修正之和。作者定义了基于 RPA 和基于 SOSEX 的外部修正（$\Delta E^{ext}_{c}$），其类比于 MP2 修正，计算为正则 RPA/SOSEX 能量与 LNO 截断对应项之间的差值。
4. 适用范围：该方法已在 PySCF 包的开发者版本中实现，适用于分子体系和周期性固体（使用 $k$ 点采样）。

主要贡献

• 形式化推广：该工作通过利用 drCCD 公式，确立了将 RPA 和 SOSEX 纳入 LNO-CC 的严格理论途径，在初始实现中避免了为密度矩阵构建而使用复杂的拉格朗日量或解析梯度公式的需求。
• 周期性扩展：作者将四次标度的 drCCD 算法扩展至具有显式 $k$ 点采样的周期性固体，使得基于 RPA 的局域关联计算能够应用于体相金属。
• 系统性基准测试：该研究提供了针对三类 MP2 已知难以处理的体系的全面基准测试：非共价分子复合物（如鸟嘌呤三聚体、 coronene 二聚体）、分子晶体（蒽）、反应势垒高度（BHDIV10 数据集）以及体相金属（Li 和 Cu）。

结果

• 非共价体系和晶体：对于 LNO-CCSD，基于 RPA 和 SOSEX 的外部修正相比基于 MP2 的修正，能显著更快地收敛至正则极限。虽然基于 MP2 的修正系统地高估关联能，但基于 RPA 的修正缓解了该误差，而基于 SOSEX 的修正仅需更少的活性 LNO（约 100 个，而 MP2 需 >400 个）即可实现近收敛。对于 LNO-CCSD(T)，RPA/SOSEX 的性能与 MP2 相当，所有方法均通过使用约 300 个活性 LNO 的外推法达到化学精度。
• 反应势垒高度：势垒高度观察到类似趋势。基于 SOSEX 的修正为 LNO-CCSD 提供了最快的收敛速度，而基于 MP2 的修正对 LNO-CCSD(T) 显示出轻微优势，这可能是因为 MP2 的微扰性质与微扰三激发修正更契合。所有方法在少于 100 个活性轨道的情况下均收敛至化学精度。
• 体相金属：改进在金属体系中最为显著。对于面心立方（FCC）Cu，基于 RPA 的外部修正仅用 200 个活性 LNO 就将误差降低至约 2 kcal/mol，而 MP2 需要约 500 个。基于 SOSEX 的修正甚至在最宽松的阈值下（约 50–200 个 LNO）即可达到化学精度。至关重要的是，随着 $k$ 点网格变密（趋近于 TDL），基于 MP2 的修正显著恶化，而基于 RPA 和 SOSEX 的修正保持稳定且准确。
• LNO 质量：发现由 RPA 振幅构建的 LNO 质量与由 MP2 构建的相当，且由于 RPA 的屏蔽效应，基于 RPA 的 LNO 通常导致稍小的活性空间尺寸。

意义
该论文确定了 RPA 和 SOSEX 作为定义片段嵌入和局域关联方法中低阶理论的有吸引力且实用的 MP2 替代方案。结果表明，虽然 MP2 在许多情况下仍适用于 LNO-CCSD(T)，但 RPA 和 SOSEX 在 LNO-CCSD 方面表现出更优越的性能，特别是在收敛速度和 MP2 在形式上失效的金属体系稳定性方面。通过使金属和非共价体系能够进行准确且线性标度的耦合簇计算，同时避免 MP2 的发散问题，这项工作解决了将高精度量子化学应用于复杂分子和凝聚相体系的关键瓶颈。作者指出，这些发现突显了低阶理论选择的关键作用，并为将超越 RPA 的方法进一步集成到局域关联框架中开辟了途径。