来源论文: https://arxiv.org/abs/2605.11470v1 生成时间: May 13, 2026 00:23

一步法相对论驱动相似重整化群多参考微扰理论：X2C-DSRG-MRPT2 深度解析

0. 执行摘要

在现代量子化学中，同时处理相对论效应（Relativistic Effects）和强电子关联（Strong Electron Correlation）是计算重元素体系（如锕系、镧系及过渡金属配合物）的核心挑战。传统的二步法（Two-step approaches）往往在处理自旋-轨道耦合（SOC）与电子关联的交织作用时力有不逮。

由 Zijun Zhao 和 Francesco A. Evangelista 开发的 X2C-DSRG-MRPT2 方法，通过将精确二分量（X2C）汉密尔顿量直接集成到驱动相似重整化群多参考微扰理论（MR-DSRG）框架中，实现了SOC效应的变分处理。该方法具有以下核心优势：

  1. 一步法处理：在 MCSCF 阶段变分地优化分子旋量，并在微扰阶段保持一致的相对论描述。
  2. 消除入侵态：继承了 DSRG 的特性，通过流参数（Flow Parameter）有效解决了传统 CASPT2 中的入侵态问题。
  3. 计算高效：采用密度泛函（Density Fitting）技术，计算标度维持在 $O(N^5)$，适用于中大型分子体系。
  4. 高精度：对第 2 至第 6 周期 p 区元素的自旋-轨道分裂（ZFS）预测精度极高，平均绝对百分比误差（MAPE）低于 7%。

1. 核心科学问题，理论基础与方法细节

1.1 核心科学问题：相对论与关联的“纠缠”

对于包含重元素的分子，其价层电子的运动速度接近光速，导致显著的标量相对论效应（收缩与膨胀）和矢量相对论效应（SOC）。同时，这些体系通常具有近简并的轨道（如 d 轨道或 f 轨道），呈现出强烈的静态关联特征。传统的处理方式通常分为：

  • 四分量方法 (4C)：基于 Dirac 方程，最精确但计算开销巨大（$O(N^{8-10})$）。
  • 二步法 (Two-step)：先计算不含自旋的态，再通过态相互作用（SISO）引入 SOC。这种方法在态密集或 SOC 极强时（如第 6 周期元素）会产生严重的系统偏差。

X2C-DSRG-MRPT2 的目标是通过“一步法”二分量框架，在降低计算复杂度的同时，提供媲美四分量方法的精度。

1.2 理论基础：X2C 汉密尔顿量

方法的起点是矩阵修改的 Dirac 方程。X2C 框架的核心在于寻找一个变换矩阵 $\mathbf{U}$，将 Dirac 汉密尔顿量解耦为正能量（电子）和负能量（正电子）分量：

$$\mathbf{U}^{\dagger}\mathbf{h}^{D}\mathbf{U} = \begin{bmatrix} \mathbf{h}^{+} & 0 \\ 0 & \mathbf{h}^{-} \end{bmatrix}$$

其中 $\mathbf{h}^{+}$ 即为包含一电子相对论修正的二分量汉密尔顿量。为了处理二电子图像变化误差（2ePC），作者采用了屏蔽核自旋-轨道 (SNSO) 近似，通过经验参数 $Q$ 缩放一电子 SOC 项：

$$h_{\mu\nu}^{+,\text{SD}} \leftarrow \left( 1 - \sqrt{\frac{Q(n_\mu, l_\mu)Q(n_\nu, l_\nu)}{Z_\mu Z_\nu}} \right) h_{\mu\nu}^{+,\text{SD}}$$

这一修正能够以极低的成本模拟二电子自旋-自旋和自旋-轨道相互作用。

1.3 关键技术：驱动相似重整化群 (DSRG)

MR-DSRG 不同于传统的微扰理论。它通过连续幺正变换 $s$（流参数）来控制 Hamilton 矩阵的对角化过程：

$$\bar{H}(s) = e^{-\hat{A}(s)} \hat{H} e^{\hat{A}(s)}$$

在 DSRG-MRPT2 中，能量修正项通过二阶微扰展开得到。其核心优势在于分母正规化：当轨道能量差趋于零时，微扰项不会发散，从而彻底解决了入侵态问题。这对于处理具有极窄能隙的重元素体系至关重要。

1.4 一步法变分实现细节

在 X2C-DSRG-MRPT2 中，分子轨道（旋量）在 X2C-CASSCF 层面进行变分优化。这意味着：

  1. 旋量优化：使用复数算术处理含有 $i$ 项的 SOC 矩阵。
  2. 态平均 (SA)：为了描述激发态，采用态平均 CASSCF 产生参考波函数。
  3. 动力学关联：在优化的旋量基础上，应用 MR-DSRG 处理动力学关联。由于旋量本身已包含 SOC 贡献，微扰过程自动捕捉了 SOC 与动力学关联的交叉效应。

2. 关键 Benchmark 体系与性能数据

2.1 p 区元素原子的零场分裂 (ZFS)

作者对第 2 周期至第 6 周期的 p 区元素进行了系统测试。这是验证相对论方法最严苛的指标之一。

  • 数据对比：对于第 6 周期元素（如 Tl, Pb, Bi, Po），传统的二步法（sf-X2C-SO）误差急剧增加，甚至给出定性错误的能级顺序。而 X2C-DSRG-MRPT2 的 MAPE 始终保持在较低水平。
  • 与 4C 方法对比：在 cc-pVTZ 基组下，X2C-DSRG-MRPT2 的平均绝对误差 (MAE) 为 $19.3 \text{ cm}^{-1}$，与极其昂贵的 4C-DSRG-MRPT2 ($15.5 \text{ cm}^{-1}$) 非常接近，但计算速度提升了数倍。
  • 动力学关联的重要性：对比 X2C-CASSCF 与 X2C-DSRG-MRPT2，发现引入动力学关联后，MAE 从 $22.3$ 降低至 $19.3 \text{ cm}^{-1}$，且误差分布显著收敛。

2.2 过渡金属原子 (Cu, Ag, Au, Sc, Y, La)

过渡金属原子的 $d$ 壳层具有强烈的关联效应。计算结果显示：

  • 对于 Au 原子（金），X2C-DSRG-MRPT2 预测的分裂值为 $12613.82 \text{ cm}^{-1}$，实验值为 $12274.01 \text{ cm}^{-1}$，误差仅为 2.8%。
  • 方法成功捕捉了从轻元素到重元素随原子序数增加而急剧增强的相对论修正。

2.3 氢化铊 (TlH) 的势能曲线 (PES)

TlH 是研究 SOC 效应的典型体系。计算展示了最低 6 个电子态的势能曲线：

  • 避越相交 (Avoided Crossings)：方法精准描述了 $^3\Pi_{0^+}$ 与 $^3\Pi_1$ 态在 $2.0$ Å 附近的交叉行为。
  • 光谱常数：TlH 的平衡键长 $R_e$ 预测值为 $1.8706$ Å，实验值为 $1.872$ Å，精度达到了 $0.001$ Å 级别。

2.4 计算性能数据

在一台普通的 MacBook Air M2 (16GB RAM) 上：

  • TlH 分子的全电子计算（含 54 个内层电子冻结）总耗时仅为 630 秒。
  • 其中 X2C-DSRG-MRPT2 微扰步骤仅占 30 秒，展示了极其优异的计算效率。
  • 内存占用约 1.2 GB，这主要归功于三索引密度泛函（Density Fitting）技术的应用。

3. 代码实现细节与复现指南

3.1 软件包支持

该方法已集成在 FORTE2 软件包中。FORTE2 是一个开源的高级多参考量子化学计算套件，主要由 Evangelista 课题组开发。

  • 代码仓库: https://github.com/evangelistalab/forte2
  • 依赖项:
    • Psi4: 提供积分引擎和基础架构支持。
    • BasisSetExchange: 用于自动获取相对论基组。
    • Numpy/opt_einsum: 优化张量收缩运算。

3.2 关键算法实现：密度泛函 (DF)

X2C-DSRG-MRPT2 的高效性在于避免了存储四中心双电子积分。在复数旋量基组下，双电子积分的存储量是实数轨道的 16 倍。作者通过以下方式优化：

  • DF 拆分：将 $(pq|rs)$ 拆分为三中心积分 $(pq|P)$ 和 $(P|rs)$。
  • On-the-fly 收缩：在计算微扰振幅 $t$ 时，直接利用三中心积分进行张量收缩，大幅降低了内存需求。

3.3 复现指南：以 TlH 为例

  1. 安装 FORTE2: 通过 Conda 或源代码编译安装最新的 FORTE2 版本。
  2. 准备输入文件:
    • 设置 molecule 块，指定 Tl 和 H 的坐标。
    • 指定基组：推荐使用针对相对论优化的基组，如 ANO-RCC-VTZP 或 cc-pVTZ-DK。
    • 设置 set forte_options { ... } 块：
      • hamiltonian = x2c
      • intrinsic_so = true (开启 SNSO 缩放)
      • df_basis = autoaux (启用自动辅助基组)
  3. 执行计算: 调用 energy('dsrg-mrpt2')。

4. 关键引用文献与局限性评论

4.1 关键参考文献

  1. DSRG 理论核心: Li, C.; Evangelista, F. A. J. Chem. Phys. 2016, 144, 204111. (奠定了 MR-DSRG 的微扰框架)
  2. X2C 汉密尔顿量: Iliaš, M.; Saue, T. J. Chem. Phys. 2007, 126, 064102. (X2C 理论的基石)
  3. SNSO 经验缩放: Ehrman, J. et al. J. Chem. Theory Comput. 2023, 19, 5785–5790. (本文采用的缩放方案来源)
  4. 四分量对照组: Zhao, Z.; Evangelista, F. A. J. Phys. Chem. Lett. 2024, 15, 7103–7110. (作者之前的 4C 工作)

4.2 方法局限性评论

尽管 X2C-DSRG-MRPT2 表现出色，但仍存在以下局限：

  • SNSO 的经验性：SNSO 参数 $Q$ 依赖于特定基组和拟合数据。虽然对于 p 区和过渡金属效果很好，但对于极端的锕系元素体系，其普适性仍需进一步验证。
  • 二电子图像变化误差：尽管通过缩放补偿了 2ePC 误差，但在原子序数极大（$Z > 100$）的超重元素中，原子平均场（AMFI）或分子平均场（X2C-MMF）方法可能更为稳健。
  • 基组依赖性：相对论计算需要去收缩（decontracted）基组以保持动力学平衡，这会导致基组线性相关问题。虽然作者通过丢弃小特征值解决了这一问题，但在超大体系中可能会影响数值稳定性。
  • 三阶微扰的缺失：目前仅实现到二阶（MRPT2）。在某些极强关联体系中，三阶修正（MRPT3）对于精确描述态能量间隔是必要的。

5. 其他必要补充：深度原理解析

5.1 为什么是 $O(N^5)$？

在量子化学中，$O(N^5)$ 是分水岭。传统的 CASPT2 如果不采用内部收缩（internal contraction），标度会随活性空间大小呈指数增长。DSRG-MRPT2 通过以下张量结构实现多项式标度：

$$E^{(2)} = \sum_{ijab} t_{ab}^{ij} \tilde{v}_{ij}^{ab} + \dots$$

由于它只需要三体还原密度累积量（3-body RDM cumulants），而不需要全组态相互作用的张量，因此对于活性空间较大的体系表现极为稳健。在本文的复数实现中，虽然复数运算带来了 4 倍的常数项增长，但其渐近标度依然锁定在 $O(N_C^2 N_V^2)$（$C$: 核心轨道, $V$: 虚拟轨道）。

5.2 图像变化 (Picture Change) 的直观理解

相对论量子化学中最令人头疼的是：当你将一电子汉密尔顿量从四分量变换到二分量时，双电子算符（库仑相互作用）理论上也应该发生变换。然而，直接变换双电子项会导致极其复杂的积分计算。SNSO 的巧妙之处在于，它意识到二电子 SOC 效应在物理上主要表现为对核电荷产生的强电场的“屏蔽”，因此可以通过修改一电子核吸引项的 SOC 部分来近似模拟这种屏蔽效果，这正是其高效的原因。

5.3 复数算术的工程挑战

将 DSRG 扩展到复数域并非简单的类型转换。所有的张量收缩（einsum）必须支持 Hermitian 对称性，且内存分配翻倍。FORTE2 在这方面做了深度优化，确保在处理复数旋量时，内存分配的冗余降到最低。这对于目前在 Apple M2 系列等统一内存架构芯片上运行大规模模拟具有重要的实践意义。

总结：X2C-DSRG-MRPT2 的问世，为强关联重元素体系的常规处理开辟了道路。它在精度、速度和稳定性之间找到了一个极佳的平衡点，是未来 routine 处理相对论效应的有力工具。