格点QCD核力与多核子系统

核力从哪来？这个问题困扰了物理学八十多年。1935年汤川秀树猜是介子交换，后来发现这个图太粗糙。1970年代量子色动力学——QCD——确立了：质子和中子是三个夸克的束缚态，它们之间的力归根结底是夸克之间交换胶子的结果。但知道方程不等于知道答案。QCD在低能区是强耦合的，微扰论失效，没人能从方程直接算出两个核子之间那个势阱长什么样。直到2007年，石井晶（Ishii）、初田哲男（Hatsuda）和青木慎也（Aoki）在PRL上发表了一篇论文，第一次从格点QCD模拟中"看见"了核力——那个中程吸引、短程排斥的形状，跟核物理学家用了几十年的唯象势对得上。从此，核力不再只是猜的，而是从第一性原理算出来的。

格点QCD怎么算核力？核心思想是把连续的时空变成离散的网格，把路径积分变成蒙特卡洛求和。夸克场定义在格点上，胶子场定义在格点之间的链上。给一个格点间距a和一个时空体积L⁴，你就可以生成一批规范场构型，然后在每个构型上算夸克传播子。两个核子的关联函数——包含四个夸克传播子的收缩——随时间衰减的方式，编码了它们之间的相互作用信息。问题是怎么把这些信息提取出来。

这里出现了两条路。第一条是Lüscher方法：在有限体积L³中，两个粒子的能量E和它们在无穷大体积中的散射相移δ之间有一个精确的解析关系。你只需要从格点关联函数中提取出能量E，就能得到相移。干净，优雅，理论上无可挑剔。第二条是HAL QCD方法：从关联函数中提取出两个核子的"波函数"ψ(r)，然后定义一个势V(r)使得它满足Schrödinger方程。有了V(r)，你就可以解散射问题、算束缚态、推广到三体甚至多体。灵活，可扩展，但定义势这一步需要理论输入——势不是可观测量，不同定义给出不同的V(r)。

两条路吵了十年。NPLQCD合作组用Lüscher方法，在mπ≈800 MeV的重夸克质量下，声称在¹S₀和³S₁两个NN道都看到了浅束缚态。HAL QCD合作组用自己的方法在同一参数下什么束缚态也没看到。谁对？2020年代的一系列工作终于搞清楚了：NPLQCD早期用的wall-source关联函数有严重的激发态污染，负权重让有效质量在中间时间出现一个假的平台——假束缚态。换用Hermitian关联矩阵后，所有组都得出了同样的结论：在重夸克质量下，NN系统没有束缚态。争议画上句号。Morningstar等人在arXiv:2504.01950中详细展示了这个假平台的产生机制：非对角关联函数中A₁=-1.05的负权重让有效质量先上升再趋于真值，中间那段看起来像基态的信号完全是假的。这个教训很深：在格点QCD中，关联函数的信噪比随时间指数恶化，而激发态污染衰减得慢，两者的竞争让plateau识别变成了一个高风险操作。

争议解决了，但两条路各有自己的领域。Lüscher方法适合计算共振参数——有限体积能级的移动直接反映散射相移，不依赖势的定义。HAL QCD方法适合多体问题——有了势函数就可以用Few-body方法算三核子力、轻核结合能、核物质状态方程，这是Lüscher方法做不到的。近年来两个组也在互相借鉴：HAL QCD用改进的关联函数减少系统误差，NPLQCD开始探索三粒子有限体积形式理论。

最大的瓶颈一直是夸克质量。真实世界的π子质量是139.6 MeV，但格点计算中减小夸克质量会让计算量暴涨——轻夸克意味着更长的关联长度、更慢的蒙特卡洛采样、更严重的临界慢化。所以直到2020年代初，大多数核力计算都在mπ=300-800 MeV下做的，离真实世界太远，需要手征外推才能跟实验比较——而外推本身引入的不确定性可能比计算还大。2024-2025年，HAL QCD合作组在Fugaku超级计算机上生成了(2+1)味Wilson-Clover规范构型，π子质量降到137.1 MeV——几乎就是物理值。格点间距a=0.084 fm，体积(8.1 fm)³，温度只有15-23 MeV。这是里程碑式的突破。在近物理质量下，NN势的中程吸引明显增强——π子越轻，单π交换力越长程，双π交换力也越强——这正是手征对称性要求的。Hatsuda在arXiv:2507.08359中展示了mπ=146 MeV下的结果：双π交换特征能量E_eff(r)在r>1 fm处精确趋于2mπ=292.8 MeV的平台，这在更重的π子质量下是看不到的。换句话说，核力中那个"中程吸引"的物理来源——两π子交换——第一次在格点QCD中直接看到了。

但物理π质量下的NN相移和散射长度至今还没发表——这是整个领域最等待的结果。一旦有了它，就可以直接跟实验比较：¹S₀道的散射长度a_t≈-23.7 fm，³S₁道的a_s≈5.4 fm，氘核结合能2.22 MeV。这些数字是检验"从QCD到核力"这条链的终极标准。

格点QCD不只有NN力。超子-核子力同样重要。直接散射实验极度稀缺——Λ束流不存在，Σ束流几乎不可能——所以格点QCD几乎是唯一的第一性原理工具。HAL QCD在mπ=146 MeV下算了ΛΛ-ΞN耦合道的势：ΛΛ吸引弱，不够支撑H-dibaryon束缚态；ΞN的¹¹S₀道最强吸引，接近unitarity。ALICE femtoscopy的Ξ⁻p关联数据不需要微调就能被格点QCD解释——这是理论与实验的直接对话。2026年3月，HAL QCD在PRL上发表了物理点的KN S波相互作用——kaon-nucleon势终于有了第一性原理数值，这对理解SIDDHARTA-2的kaonic deuterium数据至关重要。

更重的味呢？ΩΩ双超子——六个s夸克组成的dibaryon——在格点QCD中显示为束缚态。PRL 130, 111901 (2023)报道了HAL QCD的结果：全底味的Ω_bbbΩ_bbb系统强束缚，结合能远大于纯运动学效应。2025年PRL进一步确认Ω_ccc-N在物理点有吸引势。这些重味dibaryon虽然没有实验验证，但它们测试了QCD在极端夸克质量下的行为：当s夸克或b夸克变得很重，系统趋近静态极限，势的短程行为由单胶子交换主导——这是QCD最干净的区域。

三核子力呢？这是核物理的硬骨头。核物质饱和性质的定量描述、中子星状态方程的约束、轻核ab initio计算的精度——全部依赖于三体力。但三体力太弱，实验上几乎看不到它的直接效应，只能从间接线索（A=3散射、⁴He结合能、核物质饱和密度）反推。格点QCD能直接算三体力吗？原则上可以，但三核子关联函数需要九个夸克传播子的收缩，计算量比两核子大几个量级。Doi等人在2011年首次尝试：在mπ≈710 MeV的重夸克质量下，用三核子BS波函数提取3NF，但统计误差太大，只能看到3NF可能存在的微弱信号。物理质量下的3NF至今还是空白。一个间接但可行的路是"手征格点EFT"：把χEFT放到空间格点上（NLEFT），用辅助场变换把多体问题变成路径积分，然后用蒙特卡洛计算。2024年arXiv:2411.19613在N³LO上做了三核子格点EFT计算，不需要依赖特定截断的正则化——这是纯χEFT做不到的。这类方法虽然不是从QCD直接出发的，但χEFT的LEC（低能常数）可以被格点QCD的NN/YN结果约束，所以整条链是QCD→格点NN势→χEFT参数→格点EFT多体计算，间接地实现了"从QCD到原子核"。

最令人意外的应用出现在核物质状态方程上。PRL 134, 011903 (2025)——Abbott等人的工作——利用了一个QCD不等式：同位旋化学势μ_I下的配分函数Z_I给出了重子化学势μ_B=(3/2)μ_I下配分函数Z_B的上界。格点QCD在同位旋密度下没有符号问题（μ_I是实数，不引入复的费米子行列式），所以可以精确计算。他们在四个系综上做了连续极限外推，物理π质量，温度低至14.7 MeV，计算了从π子凝聚阈值到μ_I≈3250 MeV范围内的压强、能量密度和声速。声速超过共形极限c²_s/c²>1/3，确认不是格点误差。然后通过QCD不等式，把这些结果转化为对称核物质EOS的上限——这是第一个来自QCD第一性原理的核物质状态方程约束，独立于任何唯象模型。虽然目前上限还不够紧，不能排除任何主流EOS模型，但它提供了一个完全不同来源的约束，跟核结构实验和中子星观测互补。更重要的是，他们提取了同位旋超导gap：在μ_I∈[1500,3250] MeV范围内，格点结果减去不含配对的微扰QCD，残余压强与NLO微扰gap公式一致。这个gap比重子密度下的gap大——因为同位旋gap的指数前因子中1/g比baryon case小1/√2，指数增强。

还有一个里程碑：pp聚变矩阵元。北京大学冯旭、刘川团队的王紫宇等人在arXiv:2603.09554(2026)中第一次从格点QCD计算了p+p→d+e⁺+ν_e的弱作用矩阵元。这是太阳能量产生链的第一步反应，其S因子S_pp(0)的精度直接决定标准太阳模型的预言能力。之前S_pp的理论不确定性主要来自短程算符的LEC标定——现在格点QCD可以直接算这些LEC了。虽然目前还在较重的夸克质量下，物理质量外推的不确定性还大，但方法已经通了。

整个领域的驱动力是计算能力的增长。HAL QCD的物理点构型在Fugaku上生成，这台超级计算机的峰值算力超过400 PFLOPS。2025年8月，RIKEN宣布与NVIDIA合作开发FugakuNEXT，预算1100亿日元。美国的USQCD合作组也在争取FY2025-2029的计算资源，明确将核力与多核子计算列为高优先级。中国方面，华中师范大学的CCNU LQCD组在丁亨通领导下研究QCD热力学和相结构，中科院理论物理所和科大的格点组也在推进。但这些计算实在太贵了：物理π质量下一个NN关联函数需要O(10⁵)个core-hour，三核子是O(10⁷)，四核子以上还没有人敢碰。

```从QCD到原子核——链条与裂缝

裂缝标记:[!] 物理mπ下NN相移未发表[!] 3NF在格点QCD中不可行(计算量)[!] 符号问题阻断有限重子密度EOS[!] 四核子以上无直接格点结果```

回头看，格点QCD对核力的研究已经从概念验证走到了定量计算。2007年的淬火近似、重夸克质量、粗格点，到2025年的Nf=2+1全动力夸克、物理π质量、连续极限外推——每一步都是计算能力和方法学共同推动的。但距离"从QCD算出⁴He结合能"这个终极目标还很远。三核子力要靠χEFT间接推断，四核子以上连这个间接路径都没有完全打通，有限重子密度的符号问题像一堵墙挡在那里。不过，同位旋密度EOS的约束、pp聚变矩阵元的第一性原理计算、物理质量下NN势的出现——这些都是五年前还不敢想的结果。格点QCD正在一步步把核力从"猜测"变成"计算"。