Full Configuration Interaction Quantum Monte Carlo for Accurate $\textit{Ab Initio}$ Nuclear Structure Calculations

技术摘要：用于精确从头算核结构计算的全构型相互作用量子蒙特卡洛法

问题陈述
开发针对强相互作用原子核的从头算（ab initio）多体方法是现代核理论的核心挑战。虽然手征有效场论（EFT）为从量子色动力学导出核相互作用提供了系统性框架，但求解中等质量及重核的多体薛定谔方程在计算上仍然是难以实现的。

• 无心壳模型 (NCSM)： 虽然在轻核研究中取得了成功，但 NCSM 受限于模型空间维数的指数级增长，这种增长会迅速超过当前超级计算机的计算能力（达到 $\sim 10^{12}$ 维）。因此，NCSM 高度依赖于对较大系统的外推。
• 截断多体展开： 如多体微扰理论 (MBPT)、内中类相似重整化群 (IMSRG) 和耦合簇 (CC) 理论等方法具有多项式级缩放特性，能够进行较重核的计算。然而，它们受限于截断误差。为了获得高精度，需要将这些展开推向更高阶，这增加了数学复杂性和计算成本。此外，在这些截断方案中提供可靠的多体不确定性估计非常困难，特别是对于具有强团簇或集体相关性的系统，其收敛速度可能较慢。

目前迫切需要一种既能保留全构型相互作用 (FCI) 的精度，又无需显式存储哈密顿量的算法，以此作为验证和指导截断展开方法的基准。

方法论
作者将全构型相互作用量子蒙特卡洛 (FCIQMC) 方法引入到从头算核结构领域。FCIQMC 是一种随机投影方法，用于求解虚时薛定谔方程：
$$-\frac{d}{d\tau} |\Psi(\tau)\rangle = (\hat{H} - E_0) |\Psi(\tau)\rangle$$
其实现的关键特征包括：

• 随机采样： 波函数由分布在全构型空间中不同斯莱特行列式上的“行走者”（walkers，有符号粒子）系综表示。行列式的系数由该行列式上行走者的净符号数量决定。
• 传播步骤： 行走者系综在每个虚时步 ($\Delta\tau$) 内通过三个随机步骤进行演化：
  1. 生成 (Spawning)： 采样非对角哈密顿量耦合 ($H_{ij}$)，在相连的行列式上创建新的行走者。
  2. 死亡/克隆 (Death/Cloning)： 采样由参数 $S$ 偏移后的对角项 ($H_{ii}$)，调整当前行列式上行走者的数量。
  3. 湮灭 (Annihilation)： 抵消位于同一行列式上的相反符号的行走者，这对于克服费米子符号问题至关重要，且无需使用固定节点近似或引导波函数。
• 发起器近似 (Initiator Approximation)： 为了控制大模型空间中的符号问题，作者采用了发起器近似结合自适应偏移修正。这使得随着行走者数量 ($N_w$) 的增加，可以实现对 FCI 极限的可控且可改进的估计。
• 系统设置： 计算针对 $^4$He、$^8$Be、$^{12}$C 和 $^{16}$O，使用了优化后的手征两体相互作用 (N2LOopt)，并采用球谐振荡器基底 ($\hbar\omega = 20$ MeV)。在这些 FCIQMC 计算中未应用任何构型截断。

关键结果
研究给出了所选原子核的基态能量 ($E$) 和电荷半径 ($R_{ch}$)，其多体不确定性达到了亚百分比水平。

• 小空间验证： 在小模型空间 ($e_{max}=2$) 中，FCIQMC 结果与通过直接对角化获得的精确 FCI 结果完全吻合，验证了随机采样的准确性。
• 大模型空间表现： 在大模型空间 ($e_{max}=10$) 中，其中 $^{16}$O 的希尔伯特空间维度达到 $\sim 10^{33}$（使得直接对角化变得不可能），FCIQMC 仍提供了精确的结果。
  • 对于 $^4$He，FCIQMC 结果与近乎精确的 NCSM 外推结果完美一致。
  • 对于较重的原子核（$^8$Be、$^{12}$C、$^{16}$O），在 NCSM 失效的情况下，FCIQMC 提供了精确的结果。
• 截断方法的基准测试：
  • MBPT： 显示出与 FCIQMC 结果最大的偏差。
  • IMSRG： 表现总体合理，但在 $^8$Be 中表现出显著偏差，其中高阶修正 (IMSRG(3f2)+T) 甚至比精确值偏离得更远，表明存在缺失的集体相关性。
  • 耦合簇 (CC)： 展示了从 CCSD 到 CCSDT-3 的系统性改进。然而，对于 $^{16}$O 的电荷半径，包含领先迭代三体项后，其偏差反而比 FCIQMC 更大，凸显了在某些观测物理量上收敛的难度。
• 团簇效应： 对于具有显著分子/团簇结构的 $^8$Be，结果表明基于参考态的展开（如单参考态 IMSRG）由于忽略了高阶集体相关性，难以实现收敛。FCIQMC 通过包含基底内所有激发秩的相关性，保持了对这些团簇态的准确性。

意义与主张
本文声称 FCIQMC 建立了从头算核结构计算的新标准：

1. 随机全空间求解器： 它能够通过直接采样全构型空间中的精确波函数，处理超出传统 NCSM 处理范围的系统。
2. 量化的不确定性： 与难以进行不确定性估计的截断展开方法不同，FCIQMC 提供的结果具有量化的、亚百分比级的多体不确定性。
3. 理论基准： 其结果作为评估手征相互作用以及多项式级数多体方法（MBPT、IMSRG、CC）收敛性的精确基准。研究明确展示了在截断方案中经常缺失的高阶相关性的关键作用。
4. 未来影响： 作者认为，该方法将促进下一代高精度多体方法的发展，并有助于解决长期存在的核理论难题，例如钙链中电荷半径的演化以及超允许费米衰变中的同位素破缺修正，其中高阶多体效应起着核心作用。

这项工作代表了首次将 FCIQMC 应用于强相互作用下有限核的精确从头算计算，标志着在以空前精度理解核力与基本对称性方面迈出了重要一步。