One-Step Relativistic Driven Similarity Renormalization Group Multireference Perturbation Theory

本文提出了 X2C-DSRG-MRPT2，这是一种基于精确二分量哈密顿量的高效单步相对论多参考微扰理论，能够以五次方计算标度和高精度准确描述强关联体系中的自旋 - 轨道耦合效应。

要点

想象一下，你正在试图预测一个复杂舞蹈团体的行为。在化学世界中，“舞者”是电子，而它们栖身的“舞池”则是原子或分子。

长期以来，科学家在模拟含有重元素（如金、铅或铊）的分子时，主要面临两个问题：

1. “重”问题：重原子中的电子运动速度极快，其行为遵循爱因斯坦的相对论理论。这产生了一种棘手的“自旋”效应（称为自旋 - 轨道耦合），使得电子的舞步变得复杂得多。
2. “拥挤”问题：在这些重原子中，电子并非独自起舞；它们之间会产生强烈的相互影响。这被称为“强关联”。如果你试图通过一次观察一个电子来预测舞蹈，就会出错。你必须同时观察整个群体。

新解决方案：“一步到位”的舞蹈教练

本文介绍了一种名为 X2C-DSRG-MRPT2 的新计算方法。你可以将其想象为一位高效、全能且能同时解决上述两个问题的舞蹈教练。

以下是作者如何利用简单的类比来分解这一方法：

1. “精确二分量”（X2C）地图
想象一下试图在城市中导航。最准确的地图是一个 4D 全息图（代表相对论的全部复杂性），但它体积庞大、加载缓慢，且需要超级计算机。
作者使用了一张"2D 地图”（精确二分量哈密顿量）。这是一个巧妙的捷径，它捕捉了 4D 全息图的所有关键细节，但体积更小、处理速度更快。这就像使用一款高清 GPS，无需依赖一栋大楼大小的卫星，就能精准定位你的位置。

2. “驱动相似性重整化群”（DSRG）
这是处理“拥挤”电子问题的引擎。想象一个杂乱无章的房间，人们互相碰撞。

• 旧方法可能会试图通过先清理一个角落，再清理另一个角落来整理房间，但这往往会导致陷入困境或忽略整体大局。
• DSRG 方法则像一台智能清洁机器人，能够系统地抚平混乱。它不会被“闯入者”问题（即数学失效的情况）所困扰，并且具有良好的可扩展性，意味着随着房间变大，其运行速度不会呈指数级下降。

3. “一步到位”方法
这是本文最大的创新。

• “两步走”方法（旧方式）：首先，你在不考虑重相对论自旋效应的情况下计算舞步。然后，在第二步中，你将自旋效应作为修正项加入。这就像在没有音乐的情况下排练舞蹈，然后在最后试图加入节奏。这往往会导致不匹配。
• “一步到位”方法（新方式）：X2C-DSRG-MRPT2 方法在音乐（相对论）播放的同时计算舞步。它一次性优化了整个表演。论文表明，这种“一步到位”的方法要准确得多，特别是对于“音乐”最响亮的重元素而言。

他们证明了什么？

作者在各种各样的“舞者”身上测试了这一新方法：

• 单原子：从轻元素（如硼）到非常重的元素（如铊和铅）。
• 分子：如铊化氢（TlH）等原子对。

结果：

• 准确性：该方法预测的“自旋 - 轨道分裂”（不同舞步之间的能隙）与真实世界实验相比，平均误差小于 7%。对于许多体系，其准确性甚至更高。
• 效率：尽管精度极高，但其计算成本很低。其运行时间随系统规模呈合理的五次方缩放，这意味着它可以在标准计算机上运行，而无需依赖庞大的超级计算机。
• “秘密武器”：论文发现，如果你试图在主计算之后（即“两步走”或近似方法中）添加相对论效应，对于重元素而言，准确性会显著下降。你必须从一开始就将相对论和电子拥挤作为一个整体来处理。

总结

作者开发了一种新工具，使科学家能够在无需超级计算机的情况下，准确模拟重而复杂的分子。通过将“相对论自旋”和“电子拥挤”视为一个统一的单一问题，他们达到了与最昂贵方法相媲美的精度，但成本却只是其零头。

他们还指出，该方法已实现在一个名为 Forte2 的开源软件包中，这意味着其他科学家现在就可以利用它来研究重元素化学。

技术摘要

技术摘要：一步相对论驱动相似性重整化群多参考微扰理论

问题陈述
准确模拟含重元素的分子需要同时处理相对论效应（包括标量和矢量效应）与强电子关联。虽然基于狄拉克方程的四分量（4C）方法能提供高精度，但其计算成本昂贵。两分量（2C）方法，如精确两分量（X2C）哈密顿量，提供了一种更经济的替代方案，但在平衡相对论与电子关联的处理方面往往存在困难，特别是在存在强自旋–轨道耦合（SOC）的情况下。现有的单步 2C 多参考方法寥寥无几，且许多当前方法依赖于两步方案，即在轨道优化后微扰性地添加 SOC，这可能导致在具有密集低能态流形的体系中产生不准确的结果。此外，在 2C 方法中纳入双电子相对论效应（双电子图像变化误差）仍然是一个挑战。

方法论
作者提出了 X2C-DSRG-MRPT2，这是一种高效的一步相对论二阶多参考微扰理论实现。该方法整合了三个核心组件：

1. 哈密顿量：利用**精确两分量（X2C）哈密顿量将正能电子态与负能正电子态解耦。为了考虑双电子图像变化（2ePC）误差（即被忽略的双电子自旋–自旋和自旋–轨道相互作用），作者采用了屏蔽核自旋–轨道（SNSO）**近似。这涉及基于原子序数和基函数角动量对单电子哈密顿量的自旋依赖部分进行经验缩放。
2. 参考波函数：该方法使用X2C-CASSCF（完全活性空间自洽场）在 SOC 存在的情况下变分优化分子旋量。这通过态平均（SA）形式处理激发态，并采用复数算术 L-BFGS 算法进行轨道优化。
3. 关联处理：动态关联通过**驱动相似性重整化群多参考微扰理论（DSRG-MRPT2）进行处理。该形式不受闯入态问题的困扰，且仅需最多三体约化密度累积量。实现中使用了密度拟合（DF）**近似以避免存储四指标中间量，关联步骤的渐近标度为 $O(N^2_C N^2_V)$（其中 $N_C$ 和 $N_V$ 分别为芯层和虚轨道旋量），相对于活性旋量为 $O(N^6_A N_V)$。

作者还调查并比较了几种近似方案（标记为 A、B 和 C），其中 SNSO-X2C 哈密顿量在计算的不同阶段引入（例如，仅在参考态弛豫期间或在微扰步骤期间），以确定在轨道优化级别对 SOC 进行完全变分处理的必要性。

主要贡献

• 实现：作者在开源Forte2包中提供了首个高效、密度拟合的一步相对论 MR-DSRG-MRPT2 方法实现。
• 变分处理：该工作证明，在 SOC 存在下（通过 X2C-CASSCF）变分优化轨道对于精度至关重要，其表现优于仅在轨道优化后微扰性添加 SOC 的方案。
• 效率：通过将 X2C 哈密顿量与 DSRG-MRPT2 形式及密度拟合相结合，该方法实现了与系统尺寸大致五次方的计算标度，使得对包含直至第六周期元素的强关联体系常规处理相对论效应成为可能。
• 基准测试：该方法在广泛的体系范围内，与实验数据及其他高水平理论方法（包括 4C-DSRG-MRPT2/3、4C-iCIPT2 以及各种 NEVPT2 和 PDFT 变体）进行了严格基准测试。

结果

• 精度：对于包含直至第六周期元素的体系，X2C-DSRG-MRPT2 产生的自旋–轨道分裂的平均绝对百分比误差（MAPE）始终低于 7%。对于一组 15 个第二至第四周期的 p 区元素，该方法实现了 5.4% 的 MAPE，与昂贵得多的 4C-DSRG-MRPT2（4.9%）相当。
• 近似方案：研究表明，仅在 CASSCF 优化后引入 SNSO 哈密顿量的近似方案（方案 B 和 C）会导致精度的显著损失，特别是对于重元素。完全变分方法（X2C-DSRG-MRPT2）被证明更优越。
• 动态关联：通过 DSRG-MRPT2 纳入动态关联显著提高了相对于单独 X2C-CASSCF 的精度，将 p 区元素的 MAPE 从 7.8% 降低至 6.2%（使用 ANO-RCC 基组）。
• 体系性能：
  • p 区原子：该方法准确重现了从硼到碘以及铊的零场分裂，性能优于许多现有的 2C 方法，并与 4C 方法相媲美。
  • 过渡金属：对于 Cu、Ag、Au、Sc、Y 和 La，该方法实现了 4.9% 的 MAPE，与未收缩的 X2C-MRCI 相当，但所需的活性空间更小。
  • 分子：该方法对开壳层双原子分子（如 OH、SH、IO）表现良好，并准确重现了 TlH 分子的势能面和光谱常数，与实验值及高水平 4C-ic-MRCI+Q 结果吻合。
• 计算成本：TlH 的计时分析显示，在标准笔记本电脑上，DSRG-MRPT2 步骤相对于 X2C-CASSCF 步骤（约 533 秒）仅增加了适度的计算开销（约 30 秒），证明了其高效率。

意义
本文声称，X2C-DSRG-MRPT2 为强关联分子体系中相对论效应的常规处理提供了一条有前景的途径。通过以适中的计算标度（系统尺寸的五次方）实现高精度（MAPE < 7%）并避免闯入态问题，该方法弥合了昂贵四分量理论的精度与两分量方法的效率之间的差距。作者强调，在轨道优化级别对 SOC 进行变分处理对于获得可靠结果至关重要，特别是对于强 SOC 与电子关联相互作用显著的重元素。未来的工作 noted 将侧重于将此框架扩展至三阶微扰理论（DSRG-MRPT3）。