Dilated coordinate method for solving nuclear lattice effective field theory

技术摘要：用于求解核格点有效场论的扩张坐标法

问题陈述
核格点有效场论（NLEFT）是研究低能核结构与反应的标准从头算（ab initio）框架。然而，传统的 NLEFT 实现依赖于均匀立方格点，这在空间分辨率与模拟体积之间面临着根本性的权衡。

• 分辨率与体积： 准确描述短程核相互作用需要精细的格点间距（$a \approx \pi/\Lambda$，其中 $\Lambda \sim 500$ MeV）。相反，弱束缚系统（例如氘、氚、晕核）以及低能散射过程具有向远端延伸的远程波函数尾部。捕捉这些尾部需要巨大的模拟体积（$L$）。
• 计算瓶颈： 在均匀格点中，若要在保持精细分辨率以描述短程力的情况下增加体积以捕捉远程尾部，会导致格点数量呈指数级增长。对于三体系统，希尔伯特空间的维度按 $O(L^6)$ 比例缩放，这使得针对浅束缚态或靠近连续谱阈值的激发态的计算在计算上变得难以承受。标准均匀格点通常会受到严重的有限体积效应影响，需要进行大量的外推，且这种外推可能无法高效收敛。

方法论：扩张坐标法
作者引入了一种专门的自适应网格细化（AMR）技术，称为扩张坐标法。该方法并非直接使用非均匀网格，而是采用一种解析坐标变换，将均匀的计算空间（s-空间）映射到非均匀的物理空间（r-空间）。

1. 坐标变换：
   物理相对坐标 $\mathbf{r}$ 通过以下方式与计算坐标 $\mathbf{s}$ 相关联：
   $$\mathbf{r} = \mathbf{s} f(s)$$
   其中 $s = |\mathbf{s}|$ 且 $f(s)$ 为缩放函数：
   $$f(s) = 1 + (\lambda - 1) \left[ 1 - \exp\left(-s^6/R^6\right) \right]$$
   这里 $R$ 定义了相互作用区域的半径，而 $\lambda > 1$ 是渐近扩张比例。

• 内区（$s \ll R$）： $f(s) \approx 1$，因此 $\mathbf{r} \approx \mathbf{s}$。格点间距保持精细，以解析短程核力。
• 外区（$s \gg R$）： $f(s) \approx \lambda$，因此 $\mathbf{r} \approx \lambda \mathbf{s}$。物理格点间距会扩大 $\lambda$ 倍，从而通过更少的点高效地捕捉长程波函数尾部。

2. 实现：

• 将薛定谔方程从物理 $\mathbf{r}$-空间变换到均匀的 $\mathbf{s}$-空间。
• 利用解析导出的变换雅可比矩阵对动能算符进行重新表述。
• 在均匀 $\mathbf{s}$-格点上使用标准的有限差分或快速傅里叶变换（FFT）方法求解变换后的方程，随后对少体系统（$A \le 3$）进行直接对角化（Lanczos 算法）。
• 对称化： 对于全同粒子（例如三体系统），在 $\mathbf{s}$-空间对波函数进行对称化，以确保物理结果与粒子标记无关，并减轻由坐标扩张引入的格点伪影。

关键结果
论文通过在一维和三维空间中使用玩具模型及真实的核力验证了该方法：

• 一维和三维两体系统：

  • 使用高斯势和无派生有效场论（氘），扩张格点相比均匀格点表现出显著加速的无限体积极限收敛。
  • 对于氘核，扩张格点在较小的盒尺寸下即可收敛至实验结合能（$-2.224$ MeV）。即使在更大的尺寸下，均匀格点仍表现出残余的有限体积误差。
  • 该方法有效地在计算空间中压缩了波函数尾部，使得相同数量的格点可以覆盖更大的物理体积。

• 长程相互作用（库仑力）：

  • 该方法被用于测试具有无限程吸引库仑相互作用的系统。
  • 对于靠近连续谱阈值的高激发态，均匀格点表现出明显的偏离解析解现象，并由于有限体积效应导致错误的简并模式。
  • 扩张格点成功恢复了高达 $n=5$ 的状态的正确 $n^2$ 简并度和能量谱，证明了其处理长程尾部的优越性。

• 三体系统：

  • 在一维和三维三玻色子系统中，扩张法加速了深层束缚态和浅层束缚态的收敛。
  • 对于三维三玻色子系统，扩张格点在 $L/a \approx 16$ 时即达到收敛，而均匀格点则需要 $L/a \approx 20$。鉴于 $O(L^6)$ 的计算成本缩放，这一减少意味着约 4 倍的计算效率提升。
  • 两者的外推结合能差异在 $0.36\%$ 以内，证实了扩张法的准确性。

意义与主张
论文声称，扩张坐标法为研究尺度分离至关重要的核系统奠定了稳固的基础。

• 效率： 它解决了在不增加计算成本或降低精度的情况下，同时满足精细分辨率（短程力）和大规模体积（长程尾部）的需求。
• 适用性： 该方法对于以下情况特别有利：
  • 浅束缚态： 靠近连续谱阈值的系统（如晕核、$\Lambda$-超核），其波函数具有空间扩展性。
  • 激发态： 获取对边界条件敏感的高激发态。
  • 散射与反应： 该框架可直接扩展到散射过程和核反应（例如 1 MeV 以下的恒星反应），这可能实现簇相互作用（如 $\alpha$-$\alpha$ 散射）的单格点计算，其中需要在原点附近使用密集网格，而在远距离处使用稀疏网格。
• 未来展望： 虽然目前主要使用直接对角化对少体系统（$A \le 3$）进行基准测试，但作者指出该方法对多体蒙特卡洛模拟具有潜在应用价值。他们建议，对于足够短程的相互作用，仅变换动能项可能是可行的，从而为更大规模的系统保留效率增益。

这项工作并非声称立即解决所有 NLEFT 的挑战，而是提供了一个经过验证的、高效的数值工具，用以克服弱束缚和低能核系统中有限体积效应这一特定的瓶颈。