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

0. 执行摘要

传统的辅助场量子蒙特卡洛（AFQMC）方法在量子化学和凝聚态物理中有着广泛的应用。然而，一个显著的局限性在于其基于构型相互作用单双激发（CISD）试探波函数（AFQMC/CISD）的版本缺乏尺寸外延性，这极大地限制了其在大分子和材料科学领域的应用。本研究引入了一种创新的方法来克服这一挑战，即通过对耦合簇单双激发（CCSD）试探波函数进行微扰处理，构建了尺寸外延的AFQMC变体：AFQMC/ptCCSD 和 AFQMC/pt2CCSD。这两种方法在计算局部能量时均实现了O(N⁵)的计算复杂度，与AFQMC/CISD相当，但显著提高了尺寸外延性。通过对从主族分子到3d过渡金属配合物的各种基准体系的广泛数值验证，研究表明，微扰处理引入的偏差可以忽略不计。对于小体系，新方法达到了与AFQMC/CISD相当甚至更优的精度和噪声水平。特别是在处理非相互作用单体和一维原子链等大体系时，其尺寸外延性优势得到了明确体现。此外，对均匀电子气（UEG）的模拟证明，与传统CCSD(T)方法不同，新方法在热力学极限（TDL）下不会出现红外发散问题。这项工作为AFQMC在复杂大尺度量子体系中的应用开辟了新途径，并且为未来结合局部关联技术以实现线性标度算法奠定了基础。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题

辅助场量子蒙特卡洛（AFQMC）是一种强大的从头算方法，用于精确计算多电子体系的基态能量。它通过虚时间传播来投影出基态，并利用无相近似（phaseless approximation）来解决所谓的“符号问题”，从而实现多项式标度计算。然而，无相AFQMC方法引入的偏差（bias）取决于所使用的试探波函数（trial wavefunction）的质量。此前，研究组开发的AFQMC/CISD方法在计算精度上通常优于被认为是量子化学“黄金标准”的CCSD(T)方法，并且计算复杂度更低（O(N⁵) 对比 O(N⁷)）。

尽管AFQMC/CISD具有高精度和相对低的计算成本，它却面临一个核心局限性：缺乏尺寸外延性（size-extensivity）。尺寸外延性是指一个由互不作用的子系统构成的总体系的能量，应该等于所有子系统能量的总和。对于AFQMC/CISD而言，这意味着当系统尺寸增大时，其计算出的能量并不严格满足这一性质。这种缺陷导致AFQMC/CISD在处理大分子、低维材料以及具有周期性边界条件的固体等体系时变得不适用，因为在这些情况下，对关联能（以及所有广泛性质）随系统尺寸的正确缩放是至关重要的。尺寸外延性缺乏也是开发簇内分子（Cluster-in-Molecule, CIM）式线性标度方法的一个根本障碍。

因此，本研究的核心科学问题是：如何在AFQMC框架内，克服AFQMC/CISD方法的尺寸外延性不足，同时保持甚至提高其计算精度和效率？ 这要求开发一种新的试探波函数处理方案，既能保证尺寸外延性，又能避免过高的计算成本。

1.2 理论基础：AFQMC的原理与尺寸外延性分析

AFQMC基本原理

AFQMC通过以下虚时间传播公式来获取基态：

$|\Psi(\tau)\rangle = \exp(-\tau(\hat{H} - E_0))|\Phi(0))$

其中 $\hat{H}$ 是哈密顿量，$E_0$ 是基态能量的估计值。通过Trotter分解和Hubbard-Stratonovich（HS）变换，复杂的多体哈密顿量被分解为一系列单体算符，从而可以通过蒙特卡洛采样来有效传播波函数。最终的AFQMC能量由投影能量公式给出：

$E = \frac{\langle\psi_T|\hat{H}|\Psi_{AF}\rangle}{\langle\psi_T|\Psi_{AF}\rangle}$

其中 $|\psi_T\rangle$ 是试探波函数，$|\Psi_{AF}\rangle$ 是AFQMC传播得到的波函数。为了克服符号问题并实现多项式标度计算，AFQMC采用了无相近似，该近似通过试探波函数 $|\psi_T\rangle$ 来限制蒙特卡洛采样路径，从而引入了一个偏差（bias）。这个偏差可以通过提高 $|\psi_T\rangle$ 的质量来系统性地减少。

哈密顿量的Cholesky分解与平均场减法

哈密顿量中的两体项通常通过Cholesky分解为一系列单体算符的平方和：

$\hat{V} = \sum_{\gamma} \hat{v}_\gamma^2$

这允许在AFQMC传播中更有效地处理两体相互作用。为了降低随机噪声，通常还会进行平均场减法（mean-field subtraction），即通过从哈密顿量中减去一个平均场项并将其加回核心能量，来重新构造哈密顿量。

尺寸外延性分析

对于AFQMC，如果试探波函数是一个乘积态（product state），即 $|\psi_T\rangle = |\psi_{T,A}\rangle|\psi_{T,B}\rangle$ 对于两个不相互作用的子系统A和B，那么在小时间步长极限下，AFQMC能量可以保持尺寸外延性。这是因为传播子$B(x)$也会因子化为$B_A(x_A)B_B(x_B)$。然而，无相近似中的一个关键组成部分是权重更新公式中的余弦项（Eq. 16）：

$w^{(n)}_i \leftarrow w^{(n-1)}_i \left| \frac{\langle\psi_T|\hat{B}^{(n-1)}(x_i)|\phi^{(n-1)}_i\rangle}{\langle\psi_T|\phi^{(n-1)}_i\rangle} \right| \times \max(0, \cos(\theta))$

其中 $\theta = \arg \left[ \frac{\langle\psi_T|\hat{B}(x_i)|\phi^{(n-1)}_i\rangle}{\langle\psi_T|\phi^{(n-1)}_i\rangle} \right]$。如果 $|\psi_T\rangle$ 是多行列式波函数（如CISD），$\cos(\theta)$ 项不再具有乘积可分离性，导致总能量出现尺寸外延性误差。理论分析表明，这个误差随时间步长以O(dt²)标度，并且对于整个模拟以O(dt)标度，与无相AFQMC的时间步长误差相同。因此，核心挑战在于设计一种新的试探波函数，既能保持高精度，又能确保尺寸外延性。

1.3 技术难点与方法细节：微扰耦合簇试探波函数

为了解决尺寸外延性问题，本研究提出了一种创新的方法：将耦合簇单双激发（CCSD）试探波函数以微扰方式纳入AFQMC。CCSD波函数本身是尺寸外延的，因为它具有乘积可分离的形式：

$|\Psi_{CCSD}\rangle = e^{\hat{T}}|\Phi_0\rangle$, 其中 $\hat{T} = \hat{T}_1 + \hat{T}_2$ 是簇算符。

然而，直接在AFQMC中使用CCSD波函数进行能量、重叠积分和力偏置的评估计算成本过高。因此，本研究通过引入一个微扰参数 $\lambda$ 将簇算符以微扰形式处理。

我们将AFQMC能量的分子（Numerator, N($\lambda$)）和分母（Denominator, D($\lambda$)）定义为：

$N(\lambda) = \langle\Psi_0|e^{\lambda\hat{T}^\dagger}\hat{H}|\Psi_{AF}\rangle$ $D(\lambda) = \langle\Psi_0|e^{\lambda\hat{T}^\dagger}|\Psi_{AF}\rangle$

其中 $|\Psi_0\rangle$ 通常是 Hartree-Fock (HF) 基态。通过对 $N(\lambda)$ 和 $D^{-1}(\lambda)$ 进行泰勒展开，可以得到不同阶的微扰能量 $E^{[n]}$。文章推导了直到二阶的微扰能量表达式 $E^{[0]}$、$E^{[1]}$ 和 $E^{[2]}$。重要的是，这些微扰表达式被证明是尺寸外延的，通过对不相互作用子系统的各项进行分组，非尺寸外延的贡献会在求差时精确抵消。

三种近似方法

为了平衡计算成本和精度，本研究提出了三种不同的微扰AFQMC/CCSD方法：

  1. AFQMC/ptCCSD (PT)：

    • 核心思想：将簇算符 $\hat{T}_1$ 和 $\hat{T}_2$ 都视为微扰项进行处理。因此，试探波函数为HF波函数 $|\Psi_0\rangle$。能量表达式中的所有关于 $\hat{T}$ 的项都通过一阶微扰展开来计算。
    • 计算细节：能量表达式形如： $E^{[0]} = \frac{\langle\Psi_0|\hat{H}|\Psi_{AF}\rangle}{\langle\Psi_0|\Psi_{AF}\rangle}$ $E^{[1]} = \frac{\langle\Psi_0|\hat{T}^\dagger\hat{H}|\Psi_{AF}\rangle}{\langle\Psi_0|\Psi_{AF}\rangle} - E^{[0]}\frac{\langle\Psi_0|\hat{T}^\dagger|\Psi_{AF}\rangle}{\langle\Psi_0|\Psi_{AF}\rangle}$ 在实际计算中，AFQMC的引导波函数（guiding wavefunction，用于计算力偏置和施加无相约束）是 Hartree-Fock 波函数。局部能量的评估成本为 O(N⁵)。

  2. AFQMC/pt2CCSD (PT2)：

    • 核心思想：利用Thouless定理，$\hat{T}_1$ 算符可以精确地作用到HF波函数上，因此将其非微扰处理，吸收进参考波函数，形成 $|\Psi'_0\rangle = e^{\hat{T}_1}|\Psi_0\rangle$。而 $\hat{T}_2$ 算符则以微扰方式处理。
    • 计算细节：零阶和一阶能量表达式变为： $E^{[0]} = \frac{\langle\Psi'_0|\hat{H}|\Psi_{AF}\rangle}{\langle\Psi'_0|\Psi_{AF}\rangle}$ $E^{[1]} = \frac{\langle\Psi'_0|\hat{T}_2^\dagger\hat{H}|\Psi_{AF}\rangle}{\langle\Psi'_0|\Psi_{AF}\rangle} - E^{[0]}\frac{\langle\Psi'_0|\hat{T}_2^\dagger|\Psi_{AF}\rangle}{\langle\Psi'_0|\Psi_{AF}\rangle}$ 为了避免结果质量下降，尽管 $|\Psi'_0\rangle$ 看起来更自然作为引导波函数，但在实际计算中，引导波函数仍然使用HF波函数 $|\Psi_0\rangle$，而 $|\Psi'_0\rangle$ 则作为投影能量评估的试探波函数。这种方法也被证明具有 O(N⁵) 的局部能量评估成本。

  3. AFQMC/stoCCSD (STO)：

    • 核心思想：将 $\hat{T}_1$ 和 $\hat{T}_2$ 都处理到所有阶次（非微扰）。首先，将 $\hat{T}_1$ 吸收进HF参考态形成 $|\Psi'_0\rangle$。然后，利用HS变换将 $e^{\hat{T}_2}$ 算符分解为辅助场 $y$ 的积分，从而可以随机采样整个CCSD波函数。
    • 计算细节： $e^{\hat{T}_2}|\Psi'_0\rangle = \int dy \frac{1}{(2\pi)^{N_y/2}} e^{-\frac{1}{2}y^2} e^{y\cdot\hat{v}}|\Psi'_0\rangle$ 其中 $\hat{v}$ 是从 $\hat{T}_2$ 分解得到的单体算符。这种方法能够直接采样CCSD波函数并评估局部能量。为了降低噪声，采用了关联采样技术。局部能量的评估成本在形式上与微扰方法相同，但由于需要额外的采样 $N_y$，实际成本为 O(N⁵) + O(N_y N⁴)。对于小体系，STO方法可以作为评估两种微扰方法精度的参考。

计算复杂度和噪声标度

  • 计算复杂度：所有三种新方法（AFQMC/ptCCSD, AFQMC/pt2CCSD, AFQMC/stoCCSD）在局部能量评估阶段都保持了 O(N⁵) 的标度，这与AFQMC/CISD相同。传播步骤的成本为 O(N³)，与AFQMC/HF相同，因为它们都使用HF波函数作为引导态。
  • 噪声标度：
    • AFQMC/HF：总噪声随体系尺寸以 O(N¹/²) 标度，这是由于能量可以分解为N个独立体系的贡献，方差具有可加性（Eq. 40）。
    • AFQMC/CISD, AFQMC/ptCCSD, AFQMC/pt2CCSD：这些方法的局部能量表达式耦合了独立片段，导致其方差随体系尺寸以 O(N²) 标度（对于单点能量），总噪声随体系尺寸以 O(N) 标度（Eq. 42, 43, 44）。
    • AFQMC/stoCCSD：由于其随机评估 $e^{\hat{T}_2}$ 重叠，其噪声随体系尺寸呈指数增长（Eq. 45, 46, 47），这使得它对于大体系不实用。

总结来说，本研究通过引入尺寸外延的微扰CCSD试探波函数，在AFQMC框架内解决了尺寸外延性这一核心问题，并在计算效率和精度之间取得了良好的平衡。

2. 关键 benchmark 体系，计算所得数据，性能数据

为了全面验证新方法的性能，本研究在多种关键基准体系上进行了广泛的数值测试，包括非相互作用体系、相互作用体系、主族分子和3d过渡金属配合物。以下是主要结果的详细分析：

2.1 非相互作用体系：H₂ 和 O₂ 分子链

目的：验证新方法的尺寸外延性以及其噪声特性。

2.1.1 氢分子（H₂）链

  • 体系设置：采用STO-6G基组，键长1.05835 Å，由1到50个非相互作用的H₂分子组成。计算在自旋限制轨道下进行。
  • 能量结果（图1）：
    • AFQMC/CISD：随着单体数量的增加，每个H₂分子的能量偏离CCSD参考值逐渐增大，对于50个H₂分子体系，误差达到约1 mHa/分子。这清晰地表明了AFQMC/CISD缺乏尺寸外延性。
    • AFQMC/pt2CCSD（在此体系中与AFQMC/ptCCSD相同，因为H₂不含单激发）：每个H₂分子的能量保持恒定（在随机误差范围内），与CCSD参考值非常接近。这有力地证实了新方法的严格尺寸外延性。
  • 噪声结果（图2）：
    • AFQMC/HF：噪声随体系尺寸呈O(N¹/²)增长。
    • AFQMC/CISD 和 AFQMC/pt2CCSD：噪声随体系尺寸呈O(N)线性增长。对于较小的体系（1-16个H₂），AFQMC/pt2CCSD的随机噪声显著低于AFQMC/HF。此外，AFQMC/pt2CCSD的噪声水平与AFQMC/CISD相当，没有表现出明显的更高方差。

2.1.2 氧分子（O₂）链

  • 体系设置：采用STO-6G基组，键长1.20577 Å，由1到10个非相互作用的O₂分子组成。O₂分子具有自旋三重态开壳电子基态，$\hat{T}_1$ 振幅更为显著，因此适用于比较PT和PT2。
  • 能量结果（图3）：
    • AFQMC/CISD：对于5个O₂分子的体系，能量偏离参考值约5 mHa/分子，再次印证了其尺寸外延性不足。
    • AFQMC/ptCCSD 和 AFQMC/pt2CCSD：能量保持恒定。特别是AFQMC/pt2CCSD的结果更接近FCI和AFQMC/stoCCSD参考值，表明对 $\hat{T}_1$ 的精确处理是必要的。
  • 噪声结果（图4）：AFQMC/CISD 和两种微扰方法（PT和PT2）的随机噪声彼此相当，且显著低于AFQMC/HF。AFQMC/stoCCSD的噪声随体系尺寸急剧增加，再次确认其不适用于大体系。

2.2 相互作用体系：一维氢原子链

目的：在热力学极限下，评估新方法在相互作用体系中的精度和鲁棒性，特别是其对红外发散的抵抗能力。

  • 体系设置：采用STO-6G基组的一维氢原子链，固定原子间距2玻尔，链长N=4, 8, 16, 32, 64, 100。参考值来自DMRG。
  • 能量结果（图5）：
    • AFQMC/CISD：随着链长的增加，每个氢原子的能量偏离DMRG参考值越来越大（接近0.7 mHa/原子），再次凸显尺寸外延性问题。
    • AFQMC/ptCCSD 和 AFQMC/pt2CCSD：在所有链长下，这两种方法都表现出与CCSD(T)相当的精度，且与DMRG参考值保持在约0.1 mHa/原子以内，显示出卓越的尺寸外延性和精度。
  • 噪声结果（图6）：所有方法的噪声均随体系尺寸线性增长。对于较小的链，PT/PT2/CISD的噪声水平相当，且低于AFQMC/HF。随着体系尺寸增大，PT/PT2的噪声可能超过AFQMC/CISD，但仍显著低于AFQMC/HF。

2.3 均匀电子气（UEG）

目的：进一步测试新方法在低能激发体系中是否受红外发散影响。

  • 体系设置：14电子UEG体系，rs=2（弱关联）和rs=5（强关联）两种密度。通过对基组函数（平面波）数量进行逆线性拟合，外推至完全基组（CBS）极限。
  • CBS外推能量（表II，图7）：
    • rs=2：所有AFQMC方法（HF, CISD, pt2CCSD）均获得了与FCIQMC能量相当的高精度结果。
    • rs=5：AFQMC/HF 和 AFQMC/CISD的能量表现出非单调性，随着基组函数数量增加，能量可能不规则变化。这限制了它们在UEG这种具有重要低能激发体系中的应用。AFQMC/pt2CCSD则没有出现这种非单调行为，并与AFQMC/stoCCSD结果几乎一致，这强烈表明微扰近似引入的偏差可以忽略不计。
  • 红外发散（图8）：当带隙缩小（同性核-最低非占据轨道（HOMO-LUMO）带隙）时，CCSD(T)方法表现出红外发散，能量急剧偏离CCSD参考值。与此形成鲜明对比的是，所有AFQMC方法（AFQMC/CISD, AFQMC/pt2CCSD）都成功捕捉了CCSD之外的额外关联，且没有表现出这种非物理发散。这表明新方法在热力学极限下对低能激发具有鲁棒性。

2.4 主族分子体系：HEAT 和 W4-MR 数据集、N₂ 解离曲线和苯分子

目的：在实际分子体系中验证新方法的精度，并与AFQMC/CISD和CCSD(T)进行比较。

2.4.1 HEAT 数据集（表III 和 表IV）

  • 体系：26个小分子，cc-pVDZ基组。
  • 基态能量 RMSD（表III）：AFQMC/CISD RMSD最低（0.75 mHa）。AFQMC/pt2CCSD略高于AFQMC/CISD（1.1 mHa），但优于AFQMC/ptCCSD（1.4 mHa）和CCSD(T)（1.79 mHa）。AFQMC/pt2CCSD 和 AFQMC/stoCCSD 结果统计上几乎相同，再次证实微扰处理的偏差可忽略。
  • 原子化能 RMSD（表IV）：AFQMC/CISD RMSD最低（0.55 mHa）。AFQMC/pt2CCSD（0.9 mHa）优于CCSD(T)（1.42 mHa）。

2.4.2 W4-MR 数据集（表V 和 表VI）

  • 体系：17个具有强多参考特征的大分子，cc-pVDZ基组。
  • 基态能量 RMSD（表V）：AFQMC/pt2CCSD（1.7 mHa）与AFQMC/CISD（1.9 mHa）精度相当，并且两者都显著优于CCSD(T)（4.1 mHa）。在此数据集中，PT2相对于PT的改进（约2.5 mHa）更为显著，主要归因于对 $\hat{T}_1$ 振幅较大的强关联分子的精确处理。
  • 原子化能 RMSD（表VI）：结果相似，AFQMC/pt2CCSD（1.5 mHa）和AFQMC/CISD（1.6 mHa）均优于CCSD(T)（3.37 mHa）。

2.4.3 N₂ 解离曲线

  • 体系：N₂分子在不同键长下的能量（cc-pVDZ基组）。
  • 能量误差（图9）：CCSD(T)在中间键长区域（d ≈ 3 玻尔）显著低估了关联能，这是由于微扰三重激发在近简并区域不可靠。AFQMC/pt2CCSD相对于AFQMC/ptCCSD有显著改进，并且精度与AFQMC/CISD相当。
  • RMSD 和 NPE（表VIII）：CCSD(T)表现最差（RMSD 12.3 mHa, NPE 18.6 mHa）。AFQMC/pt2CCSD显著优于PT，并达到了与AFQMC/CISD相当的精度（RMSD 2.3 mHa, NPE 4.4 mHa）。PT2和STO之间的差异可忽略。

2.4.4 苯分子

  • 体系：苯分子的基态能量（cc-pVDZ基组）。
  • 能量（表VII）：所有三种QMC方法（AFQMC/CISD, AFQMC/ptCCSD, AFQMC/pt2CCSD）在1 mHa内保持一致，且显著优于CCSD(T)，与DMRG和CCSDTQ参考值高度吻合。

2.5 3d 过渡金属配合物

目的：评估新方法在具有显著多参考特征的过渡金属体系中的性能。

2.5.1 过渡金属单氧化物（ScO, TiO, VO, CrO, MnO, FeO, CuO）

  • 体系：七种3d过渡金属单氧化物的CBS外推解离能。
  • RMSD（表IX，图S6）：AFQMC/ptCCSD 和 AFQMC/pt2CCSD 的性能劣于CCSD(T) 和 AFQMC/CISD。这可能归因于在无相近似中使用的HF引导波函数在处理这些强多参考体系时的局限性。

2.5.2 Cu₂O₂⁺异构化能

  • 体系：Cu₂O₂⁺的f=0和f=1几何结构之间的异构化能。
  • 能量（表X）：AFQMC/pt2CCSD与自由投影AFQMC（fp-AFQMC）和AFQMC/CISD高度吻合。而CCSD(T)高估了异构化能约6 kcal/mol，AFQMC/HF高估了约10 kcal/mol。尽管AFQMC/pt2CCSD和AFQMC/CISD都受益于误差抵消，但AFQMC/pt2CCSD的绝对能量与fp-AFQMC参考值更为接近。

2.5.3 Fe(H₂O)₂⁺自旋劈裂能

  • 体系：Fe(H₂O)₂⁺的低自旋（LS）单重态和高自旋（HS）五重态之间的能量差。
  • 能量（表XI）：在随机不确定性和缺乏精确参考值的情况下，CCSD(T), AFQMC/CISD, AFQMC/ptCCSD 和 AFQMC/pt2CCSD 之间存在广泛的一致性。

总结：AFQMC/pt2CCSD在大多数基准体系上持续表现出比AFQMC/ptCCSD更准确的结果，这表明精确纳入 $\hat{T}_1$ 是至关重要的。同时，AFQMC/pt2CCSD与AFQMC/stoCCSD的结果统计上几乎相同，这表明在CCSD试探波函数框架内，微扰处理引入的误差可以忽略不计，且无需超越PT2级别来捕捉额外关联。新方法在处理大体系时展现出卓越的尺寸外延性和对红外发散的鲁棒性，使其成为解决传统方法难题的有力工具。

本研究提出的尺寸外延AFQMC方法，其核心在于将耦合簇（Coupled Cluster, CC）试探波函数以微扰方式纳入到AFQMC框架中。这一创新性的实现得益于对现有量子化学软件包的深入理解和灵活扩展。

3.1 核心实现思想

  1. 微扰处理 $\hat{T}$ 算符：AFQMC/ptCCSD 和 AFQMC/pt2CCSD 的关键在于，并非直接使用完整的 $e^{\hat{T}}|\Phi_0\rangle$ 作为试探波函数进行能量评估，而是将其中的 $\hat{T}$ 算符（或其部分 $\hat{T}_2$）视为微扰项，并通过泰勒展开的形式，计算其对AFQMC能量的修正。这样避免了直接评估 $e^{\hat{T}}$ 算符在AFQMC中带来的高昂计算成本和复杂性。
  2. 引导波函数（Guiding Wavefunction）：为了保持计算效率和避免更高的传播成本，所有新方法在AFQMC的传播阶段都采用了Hartree-Fock（HF）波函数作为引导波函数。HF波函数具有乘积可分离性，使得传播和力偏置计算的成本能够保持在O(N³)级别。
  3. 局部能量评估：
    • AFQMC/ptCCSD (PT)：能量评估涉及对 $\langle\Psi_0|\hat{H}|\Psi_{AF}\rangle$、$\langle\Psi_0|\hat{T}^\dagger\hat{H}|\Psi_{AF}\rangle$ 和 $\langle\Psi_0|\hat{T}^\dagger|\Psi_{AF}\rangle$ 等项的计算。这些项可以通过对AFQMC波函数 $|\Psi_{AF}\rangle$（它是一个行列式集合）进行适当的投影和积分来高效完成。其中，$\hat{T}^\dagger\hat{H}$ 项的计算复杂度是决定O(N⁵)标度的主要因素。
    • AFQMC/pt2CCSD (PT2)：将 $\hat{T}_1$ 吸收进HF波函数后，形成新的参考态 $|\Psi'_0\rangle = e^{\hat{T}_1}|\Psi_0\rangle$。能量评估涉及 $\langle\Psi'_0|\hat{H}|\Psi_{AF}\rangle$、$\langle\Psi'_0|\hat{T}_2^\dagger\hat{H}|\Psi_{AF}\rangle$ 和 $\langle\Psi'_0|\hat{T}_2^\dagger|\Psi_{AF}\rangle$ 等项。与PT类似，这些项的计算也保持O(N⁵)标度。
    • AFQMC/stoCCSD (STO)：这是一种非微扰方法，用于基准测试。其通过对 $\hat{T}_2$ 算符进行Hubbard-Stratonovich分解，引入辅助场 $y$ 来随机采样 $e^{\hat{T}_2}$。虽然在形式上局部能量评估成本与微扰方法相同，但由于需要额外的采样步骤，实际计算成本会更高。为减少噪声，需要实施关联采样（correlated sampling）策略。

本研究所依赖的计算框架和工具都基于开源项目，这极大地促进了研究的可复现性和进一步发展。

  • PySCF：

    • 功能：所有的平均场（Mean-Field）和耦合簇（Coupled Cluster）计算都是使用PySCF软件包完成的。PySCF是一个功能强大且高度灵活的开源量子化学平台，支持各种波函数方法和密度泛函理论计算。
    • 开源链接：https://github.com/pyscf/pyscf.git

  • AFQMC 接口：

    • 功能：本研究开发了一个定制的AFQMC接口，使其能够与PySCF无缝连接，从而利用PySCF提供的CC振幅和积分。这是实现文中描述的微扰AFQMC方法的核心组件。
    • 开源链接：https://github.com/zh-yichi/afqmc

  • 均匀电子气（UEG）计算脚本：

    • 功能：用于复现文中所有均匀电子气（UEG）计算的脚本和数据。这些脚本对于理解如何在AFQMC框架下处理周期性边界条件体系以及验证新方法在热力学极限下的鲁棒性至关重要。
    • 开源链接：https://github.com/zh-yichi/ueg

  • 数值结果和数据分析脚本：

    • 功能：包含本研究中所有数值结果、数据分析和用于复现每个计算的脚本。这个仓库是理解和验证研究发现的关键资源。
    • 开源链接：https://github.com/zh-yichi/afqmc_ptccsd_data

  • TROT（未来稳定版AFQMC功能）：

    • 功能：未来，本研究中描述的AFQMC的稳定特性将发布在TROT项目上。TROT是PySCF-Forge的一部分，致力于提供高性能的量子蒙特卡洛工具。
    • 开源链接：https://github.com/ankit76/trot.git

3.3 复现指南

复现本研究结果的步骤如下：

  1. 安装 PySCF：首先，按照PySCF官方GitHub仓库 https://github.com/pyscf/pyscf.git 上的说明安装PySCF。
  2. 克隆 AFQMC 接口和相关数据仓库：
    • 克隆 https://github.com/zh-yichi/afqmc 仓库，其中包含了AFQMC的核心接口和实现。
    • 克隆 https://github.com/zh-yichi/ueg 仓库，用于UEG计算。
    • 克隆 https://github.com/zh-yichi/afqmc_ptccsd_data 仓库，其中包含了所有基准测试的输入文件、脚本和数据分析工具。
  3. 准备环境：确保你的Python环境配置正确，并安装了所有必要的依赖库（通常PySCF的安装会处理大部分依赖）。
  4. 运行计算：
    • 进入 afqmc_ptccsd_data 仓库，查阅其中的 README 文件或相关文档，了解如何针对特定体系运行计算。
    • 通常，你需要运行PySCF脚本来生成初始的HF轨道和CC振幅。然后，这些数据将作为AFQMC计算的输入。
    • AFQMC计算将调用 zh-yichi/afqmc 中的代码，执行虚时间传播、蒙特卡洛采样和局部能量评估。对于微扰方法，这些计算将涉及到对CC振幅的微扰展开。
    • 对于AFQMC/stoCCSD，你需要指定额外的辅助场采样次数 ($N_y$)。
  5. 数据分析：使用 afqmc_ptccsd_data 仓库中提供的数据分析脚本来处理原始AFQMC输出，计算平均能量、标准差以及生成图表，与论文中的结果进行比较。
  6. 并行计算：考虑到AFQMC计算的计算密集性，建议在高性能计算集群上利用多核或多节点进行并行计算。本研究中的计算使用了ACCESS项目CHE240162提供的计算资源。

通过遵循这些步骤，研究人员应该能够复现论文中的主要结果，并进一步探索这些新方法在不同体系中的应用。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献

本研究工作在AFQMC和耦合簇理论的交叉领域取得了显著进展，因此其引用文献涵盖了这两个方向的经典和最新研究。

  1. AFQMC基础与符号问题 (Ref. 1-7):

    • [1-4] Zhang et al. 的早期工作奠定了AFQMC在多体量子系统中的基础，特别是其在真实和模型系统中的应用。
    • [5,6] Feynman 和 Loh Jr. 等人的工作揭示了虚时间传播中的符号问题，这是AFQMC面临的核心挑战。
    • [7] Zhang 和 Krakauer 的论文引入了无相近似（phaseless approximation），为克服符号问题提供了有效的多项式标度解决方案，是现代AFQMC方法的基础。

  2. AFQMC/CISD与多行列式试探波函数 (Ref. 9-13):

    • [9-11] Wouters, Mahajan, Xiao 等人的工作探讨了使用更复杂的试探波函数（如矩阵乘积态、选择构型相互作用和变分AFQMC）来提高AFQMC精度。
    • [12] Mahajan et al. (本研究作者之一的早期工作)：这篇论文首次提出了AFQMC/CISD方法，并展示了其在精度上超越CCSD(T)的潜力，是本研究的直接前驱。本研究旨在解决AFQMC/CISD的尺寸外延性问题。
    • [13] Kjønstad et al. 的工作进一步探讨了截断构型相互作用波函数作为试探波函数。

  3. 耦合簇理论 (Ref. 14-17):

    • [14,15] Paldus & Li 和 Bartlett & Musiał 的综述文章，提供了耦合簇理论的全面背景，特别是CCSD(T)作为量子化学黄金标准的地位及其O(N⁷)的计算复杂度。
    • [16,17] Lyakh et al. 和 Nooijen et al. 的工作讨论了多参考耦合簇及其在处理强关联体系中的重要性。

  4. 尺寸外延性与AFQMC (Ref. 23):

    • [23] Lee et al.：这篇论文深入探讨了AFQMC能量在小时间步长极限下的尺寸外延性，特别指出乘积态试探波函数的重要性，并为本研究中尺寸外延性的理论分析提供了基础。

  5. UEG与红外发散 (Ref. 28, 37-40):

    • [28, 37-40] Masios, Shepherd, Neufeld 等人的工作讨论了MP2和CCSD(T)方法在带隙缩小体系中（如UEG在热力学极限下）的红外发散问题，强调了开发对这一问题鲁棒的新方法的重要性。本研究的新方法成功避免了红外发散。

  6. DMRG作为1D体系的参考 (Ref. 26, 27):

    • [26,27] Motta 等人的工作展示了DMRG在提供1D系统（如氢链）近乎精确解方面的能力，被用作本研究的基准参考。

  7. 过渡金属体系的挑战 (Ref. 75, 76):

    • [75,76] Cramer 和 Truhlar 的工作讨论了过渡金属配合物在量子化学中的计算挑战，强调了其多参考特性和d轨道电子的复杂性。

  8. PySCF软件包 (Ref. 89, 90):

    • [89,90] Sun et al. 的工作介绍了PySCF，这是本研究中所有平均场和耦合簇计算所依赖的开源量子化学平台。

4.2 对这项工作局限性的评论

尽管本研究提出的尺寸外延AFQMC/ptCCSD和AFQMC/pt2CCSD方法在克服尺寸外延性、提高精度和应对红外发散方面取得了显著进展，但仍存在一些局限性，值得进一步探讨和改进。

  1. 引导波函数对过渡金属体系的限制：

    • 问题：在对3d过渡金属单氧化物（表IX）的基准测试中，AFQMC/ptCCSD和AFQMC/pt2CCSD的性能劣于CCSD(T)和AFQMC/CISD。论文将此归因于在无相近似中使用的Hartree-Fock（HF）引导波函数。对于具有强多参考特征或存在多个对称性破缺HF解的体系，HF引导波函数可能无法充分捕捉体系的复杂性，导致力偏置（force bias）不够准确，进而影响AFQMC的最终结果。
    • 潜在影响：这限制了新方法在处理某些特定类型的强关联体系（如许多过渡金属配合物和开壳层体系）时的普适性。
    • 未来方向：需要开发更先进的引导波函数（例如，多构型或自优化引导波函数）或尺寸外延的无相近似方案，以弥补这一缺陷。

  2. 微扰处理的固有偏差：

    • 问题：尽管数值结果表明微扰处理（特别是PT2）引入的偏差可以忽略不计，且与完全非微扰的AFQMC/stoCCSD结果高度一致，但从理论上讲，微扰展开本身就是一种近似。当试探波函数恰好是精确基态时，零方差原理（zero-variance principle）保证AFQMC能量是精确的且没有随机噪声。微扰处理可能会打破这一原理，从而引入一个非零的偏差。
    • 潜在影响：对于一些极其敏感或需要极高精度的体系，这种微扰偏差可能变得不可忽略。此外，随着微扰阶数的增加，表达式会变得越来越复杂，计算成本也会随之增加。

  3. 噪声标度的限制：

    • 问题：虽然AFQMC/ptCCSD和AFQMC/pt2CCSD的噪声随体系尺寸呈O(N)线性增长，这优于AFQMC/HF的O(N¹/²)（总噪声），但在非常大的体系中，为了达到相同的统计精度，O(N)的噪声标度仍意味着需要更多的蒙特卡洛样本，最终可能使其计算成本高于AFQMC/HF。
    • 潜在影响：这可能限制了方法在处理尺寸极其庞大的体系时的效率，尤其是在需要非常低统计噪声的场景。
    • 未来方向：结合局部相关技术（如CIM）来将总体系分解为更小的、独立的局部区域，可能有助于降低实际噪声的增长速度，甚至实现线性标度。

  4. AFQMC/stoCCSD的实际应用性：

    • 问题：AFQMC/stoCCSD作为完全非微扰的CCSD试探波函数方法，在小体系上用于验证微扰方法的精度。然而，其噪声随体系尺寸急剧增加，这使得它在当前形式下对于大体系的应用变得不切实际。
    • 潜在影响：这意味着AFQMC/stoCCSD主要是一个理论基准，而不是一个通用的实际计算工具。为了使其更具实用性，可能需要进行显著的方法改进，例如结合相似变换。

  5. 计算成本的进一步优化空间：

    • 问题：尽管O(N⁵)的计算复杂度优于CCSD(T)的O(N⁷)，但对于真正的“大规模”或“线性标度”应用，O(N⁵)仍然相对较高。本研究本身并未实现线性标度。
    • 未来方向：论文中也提到了与局部关联技术（CIM）的兼容性。结合CIM等方法，将AFQMC/pt2CCSD进一步发展为线性标度算法，是其应用于超大体系和材料科学的关键。

  6. 微扰方法的普适性：

    • 问题：虽然微扰方法在概念上是通用的，可以推广到其他相关的波函数形式，但实际实现和分析每种具体情况的复杂性可能会很高。例如，更高阶的CC方法（如CCSDT, CCSDTQ）的微扰处理将涉及更复杂的项。
    • 潜在影响：这表明在将这种方法推广到更复杂的波函数或体系时，仍需进行大量的理论和计算工作。

总而言之，这项工作为AFQMC方法带来了重要的进步，但同时也指明了未来研究的几个关键方向，特别是在引导波函数选择、噪声控制和与线性标度技术的结合方面。

5. 其他必要补充

5.1 本研究的意义与影响

本研究提出的尺寸外延AFQMC/ptCCSD和AFQMC/pt2CCSD方法，在量子化学和凝聚态物理领域具有深远的意义和影响：

  1. 拓展AFQMC的应用范围：

    • 解决尺寸外延性瓶颈：通过解决AFQMC/CISD缺乏尺寸外延性的问题，新方法使得AFQMC能够可靠地应用于以前受限的大分子体系、低维材料、固体以及具有周期性边界条件体系。这对于精确预测材料性质（如晶格能、表面吸附能等）至关重要。
    • 处理热力学极限问题：新方法在均匀电子气（UEG）体系中展现出对红外发散的鲁棒性，这意味着它可以用于研究金属系统和绝缘体-金属转变等具有小带隙或无带隙的体系，而传统微扰方法如CCSD(T)则会在此类体系中失效。

  2. 高性能与高精度计算的融合：

    • 效率与精度兼备：新方法在局部能量评估中保持了O(N⁵)的计算复杂度，这远低于CCSD(T)的O(N⁷)。同时，在大多数基准测试中，其精度与AFQMC/CISD相当或更优，甚至超越了CCSD(T)，为精确量子化学计算提供了更优的平衡点。
    • 提高多参考体系描述能力：虽然在某些过渡金属体系中仍有改进空间，但在N₂解离曲线和Cu₂O₂⁺异构化能等强关联体系中，新方法展现出优于CCSD(T)的性能，这对于描述化学键断裂、激发态和磁性体系至关重要。

  3. 促进线性标度AFQMC的发展：

    • 与局部关联技术的兼容性：尺寸外延性是实现线性标度（O(N)）计算的先决条件。新方法固有的尺寸外延性使其与局部关联技术（如簇内分子（CIM）方法）天然兼容。这意味着未来可以结合这些技术，将AFQMC/pt2CCSD推广到真正的大尺度体系，实现其计算成本随体系尺寸线性增长。
    • 为材料设计提供工具：线性标度的精确量子蒙特卡洛方法将为新型材料的理论设计和模拟提供强大工具，从而加速科学发现和技术创新。

5.2 与现有方法的对比

本研究提出的方法在以下几个关键方面与现有方法形成了对比：

  1. 对比 AFQMC/CISD：

    • 优势：核心在于解决了AFQMC/CISD的尺寸外延性不足问题。尽管AFQMC/CISD在小分子上表现优异，但其在处理多体效应在大体系中累积时会产生非物理误差。新方法通过微扰CCSD试探波函数，在保持O(N⁵)计算成本的同时，严格确保了尺寸外延性。
    • 精度：在小分子体系（如HEAT、W4-MR）上，新方法的精度与AFQMC/CISD相当或略低。但在非相互作用分子链和一维氢链等体系中，新方法在尺寸外延性方面表现出压倒性优势，精度保持恒定，而AFQMC/CISD则随体系尺寸增大而偏离参考值。

  2. 对比 CCSD(T)：

    • 计算效率：新方法（O(N⁵)）在计算复杂度上显著优于CCSD(T)（O(N⁷)），使得其能够处理更大尺寸的体系。
    • 精度：在许多基准测试（HEAT、W4-MR、N₂解离曲线、苯分子）中，新方法展现出与CCSD(T)相当甚至更优的精度。特别是在强关联体系（如N₂解离的中间键长区域）和低能激发体系（如UEG的红外发散问题）中，新方法表现出卓越的鲁棒性，而CCSD(T)则会失效或表现不佳。
    • 尺寸外延性：CCSD(T)本身是尺寸外延的，但其在处理红外发散问题时的失败限制了其在某些热力学极限应用。

  3. 对比 AFQMC/HF：

    • 精度：AFQMC/HF使用简单的HF波函数作为试探波函数，其精度通常远低于包含关联效应的AFQMC/CISD或新方法。新方法能够捕捉更强的电子关联效应，提供更精确的能量。
    • 噪声标度：AFQMC/HF的噪声总方差标度为O(N¹/²)，理论上优于新方法的O(N)。但在实践中，对于中等大小的体系，新方法的实际噪声水平可能更低，因为其试探波函数质量更高。

  4. 对比其他AFQMC变体（如VAFQMC, Selected CI-AFQMC, MPS-AFQMC）：

    • 这些方法通常旨在通过更复杂的试探波函数来提高精度，但可能带来更高的计算成本或不同的复杂性。本研究的优势在于通过微扰处理，以相对较低的成本实现了尺寸外延性，并在O(N⁵)的计算框架内提供了高精度。

5.3 未来研究方向

本研究为AFQMC领域开辟了激动人心的未来研究方向：

  1. 改进引导波函数和无相近似方案：

    • 针对过渡金属配合物等强关联体系，目前的HF引导波函数存在局限性。未来的工作可以探索使用多构型参考态、CASSCF波函数或自优化引导波函数来提高力偏置的准确性。
    • 开发新的尺寸外延无相近似方案，可能进一步降低或消除无相近似引入的偏差，同时保持尺寸外延性。

  2. 实现线性标度AFQMC/pt2CCSD：

    • 本研究已证明新方法具有尺寸外延性，这为结合局部关联技术（如簇内分子（CIM）方法）实现线性标度AFQMC/pt2CCSD奠定了基础。这将使得AFQMC能够处理数千甚至上万个原子的体系，从而直接应用于复杂材料体系。
    • 正在进行的将AFQMC/pt2CCSD与局部关联变体结合的工作，是这一方向的重要一步。

  3. 探索更高阶的微扰处理：

    • 尽管目前的研究表明PT2对于CCSD试探波函数已足够精确，但在某些极端情况下，更高阶的微扰修正可能仍有益处。未来的工作可以探索更高阶的微扰项对精度和计算成本的影响。

  4. 推广微扰方法到其他关联波函数：

    • 本研究的微扰方法在概念上是通用的，未来可以将其推广到其他类型的关联波函数，例如广义耦合簇（GCC）、多参考构型相互作用（MRCI）或变分量子蒙特卡洛（VMC）中的Jastrow因子，以进一步提高AFQMC的灵活性和精度。

  5. 应用于更广泛的体系：

    • 将新方法应用于更多元的体系，包括拓扑材料、催化反应机制、生物分子系统等，以进一步验证其普适性和可靠性。
    • 特别是在研究相变、带隙工程和非常规超导等材料物理问题时，新方法在热力学极限下的鲁棒性将发挥关键作用。

5.4 总结与展望

这项工作代表了AFQMC领域的一个重要里程碑，它成功地解决了现有方法在尺寸外延性和热力学极限应用方面的核心挑战。通过巧妙地结合微扰耦合簇理论，研究人员开发出了一种高效、高精度且尺寸外延的量子蒙特卡洛方法。其卓越的性能和广阔的应用前景，预示着AFQMC将在未来的量子化学和凝聚态物理研究中扮演越来越重要的角色，特别是在推动大规模、复杂体系的精确模拟方面。随着进一步的优化和与其他先进计算技术的融合，本方法有望加速新材料的发现和基础科学问题的解决。