Full Uncertainty Quantification of Sign-Problem-Free Quantum Monte Carlo Methods and Nuclear Lattice Effective Field Theory Benchmarks

技术摘要：无符号问题量子蒙特卡洛方法的全不确定度量化与核格点有效场论基准

问题陈述
现代从头算（ab initio）核理论通过结合手征有效场论（EFT）与多体求解器，在重现轻核及无限核物质的可观测物理量方面取得了显著成功。然而，将这些方法扩展到中重质量核（最高至 $^{100}$Sn）面临着一个关键瓶颈：缺乏对系统算法误差的严格不确定度量化（Uncertainty Quantification, UQ）。若缺乏此类控制，理论与实验之间的差异往往会被归因于底层的核力，从而导致低能常数（LECs）的“重新调整”以补偿求解器缺陷。这掩盖了理论误差的真实来源，并损害了从头算方法的预测能力。

具体而言，最近的批评（文献 [1]）对核格点有效场论（NLEFT）基准的有效性提出了质疑，指出了其在变分界限、相关能以及核饱和机制方面存在的问题。这些批评通常源于将不同的正则化格点理论（例如，混合转移矩阵与哈密顿量形式）进行比较，或是由于有限体积处理不一致，而非源于多体方法本身的失效。

方法论
作者采用了无符号问题的量子蒙特卡洛（QMC）方法，具体为辅助场量子蒙特卡洛（AFQMC），并将其应用于 NLEFT 框架内。本研究重点关注了近期开发的 LAT-OPT1 作用量，这是一种极简的无符号问题格点哈密顿量，它在保持偶偶核所需的反幺正对称性的同时，纳入了非微扰自旋-轨道相互作用。

该方法涉及全面的误差预算分析和基准比较的诊断性再检查：

1. 不确定度量化： 作者系统地量化了从 $^4$He 到 $^{100}$Sn 核的统计、虚时、有限体积、时间离散、辅助场诱导以及残余符号不确定度。
2. 非微扰自旋-轨道处理： 不同于以往将自旋-轨道力作为微扰插入的处理方式，LAT-OPT1 将自旋-轨道密度直接嵌入到欧氏时间投影哈密顿量中。
3. 正则化一致性检查： 本文严格比较了不同的计算形式：
   • 转移矩阵 vs. 哈密顿量： 分析有限时间格点间距（$a_t$）下的转移矩阵计算与 $a_t \to 0$ 哈密顿量极限之间的差异。
   • 边界条件： 研究周期性边界条件（PBC）与平均扭曲边界条件（ATBC）对核物质中有限体积壳层效应的影响。
   • 重整化： 考察连续极限以及短程相互作用的重整化。

核心贡献与结果

• 系统误差的定量控制：
  研究表明，对于高达 $^{100}$Sn 的双幻核，其多体计算的累积不确定度远低于百分之一个数量级。

  • 虚时与有限体积： 向无限投影时间（$\tau \to \infty$）和无限体积（$L \to \infty$）的外推是稳定的。研究证明，对于所研究的核，盒尺寸 $L=11$（配合 $a=1.32$ fm）已足够收敛。
  • 时间离散： 使用具有有限 $a_t$ 的转移矩阵形式会诱导出与目标哈密顿量 $H$ 不同的有效哈密顿量 $H_{eff}$。作者展示了虽然 $H_{eff}$ 依赖于 $a_t$，但能量误差随 $O(a_t^2)$ 缩放，从而实现了精确控制。
  • 辅助场诱导项： 为了维持无符号问题采样所需的平方完成过程，生成了诱导的二体（$\rho_s^2$）和三体项。这些项经过量化后发现仅为微小的修正（数 MeV），证实了该构造的有效性。
  • 残余符号问题： 对于奇 $A$ 和奇奇核，存在轻微的残余符号问题。然而，平均相位仍接近于 1（0.94–0.98），且统计不确定度与偶偶核相当。

• 自旋-轨道耦合的非微扰作用：
  一个关键发现是，自旋-轨道相互作用必须在欧氏时间投影中以非微扰方式包含。仅将其作为 SU(4) 对称波函数的微扰插入，会遗漏在中质量核中达数至十 MeV 的壳层结构贡献。通过将其嵌入传播过程，可以动态生成正确的壳层结构。

• 解决基准差异：
  本文通过识别不同格点理论之间在概念和技术上的错误，解决了近期的批评（文献 [1]）：

  • 不等价比较： 将来自 $a_t \to 0$ 哈密顿量的 Hartree-Fock (HF) 界限与来自有限 $a_t$ 转移矩阵的 AFQMC 结果进行比较是不合理的。有限 $a_t$ 转移矩阵定义了一个带有诱导相互作用的截然不同的正则化理论。当使用一致的正则化理论时，变分原理成立（HF 能量界定投影能量）。
  • 相关能： “基本要素”（EE）相互作用在中子物质中较小的相关能并不意味着波函数是平凡的。对不可约密度的详细分析显示，尽管能量差异很小，但仍存在显著的多模超流性（s 波、p 波及四重态相关性）。
  • 核饱和： 高密度下相互作用能的平坦化并非格点填充的人为产物。在连续空间中使用伽利略不变正则化进行的解析推导表明，相互作用能的平坦化是动量空间正则化的自然结果。NLEFT 相互作用的饱和特性追踪的是弥散诱导的正则化标度，而非格点位点占据标度。
  • 连续极限： 关于格点相互作用导致系统在连续极限下坍缩的说法被驳回。用于批评的缩放论证依赖于错误的 Born 能级重整化；正确的向物理散射可观测量进行的重整化会产生不同的缩放行为（$C_c \propto a$），从而防止坍缩。

意义与主张
本文声称为具有全手征力的 NLEFT 模拟建立了坚实的完整误差分析基础。通过严格量化不确定度并证明无符号问题方法可以在重核中实现亚百分比精度，这项工作验证了 NLEFT 是可靠的多体方法基准。

作者断言，当在具有一致参数、体积和边界条件的等效正则化理论之间进行比较时，近期文献中的表象差异均得到了解决。LAT-OPT1 作用量不仅被呈现为一个唯象工具，更是一个能够捕捉本质壳层结构动力学的精确可解的非微扰基准模型。该框架使得区分真实的核力缺陷与多体求解器误差成为可能，为中重质量核的壳层演化、光谱学及致密物质的研究提供了路径。