三核子力的旋钮账本
两个核子之间的力,被几万条散射数据钉住了。三个核子之间的力,一条"三体散射数据"都没有。
这就是核力问题上最别扭的地方。两体力你可以拿相移分析往死里拟合,一套相互作用能带上两千多个数据点、误差压到百分之一量级。等到三个核子凑一块儿,你既没法造一个"三核子散射实验",也说不清它到底该有几个自由参数、每个参数该由谁来定。
三核子力不是可有可无的补丁。没有它,核物质不会饱和。只有两体力的从头算,要么束缚不住原子核,要么把核物质压到塌缩——两头都错。所以这几十年,做核力的人一直在干两件事:把三核子力的算符清单列干净,和找新的观测量来给它定系数。
这一篇只盯住这两个问题,别的不碰。
知识库里的分工说清楚。《核力与核结构》讲力长什么样、怎么从数据里反推出来。《核有效场论手征微扰论》讲这套有效理论怎么搭架子、幂次计数怎么数。《核物质EOS与中子星结构》和《星震学与引力波约束EOS》讲状态方程本身和星体怎么算。两篇《PREX-CREX悖论与对称能》讲对称能这个量从哪来、为什么两批实验打架。《少体物理与Faddeev方程》讲三体、四体怎么解。这些都不重复。
这一篇看的是账本:三核子力在理论里到底有几个旋钮,每个旋钮站在哪一级,谁在账上动了手脚,以及——中子星这种几十公里外的东西,怎么反过来给地上的系数标价。
"三"为什么这么难
先看量级对比,才知道难在哪。
两核子这一侧:有几万条散射数据、上千条束缚态性质,理论拟合的自由参数五六十个,摊到每个参数头上几十个数据点。所以 NN 相互作用被压到百分之一左右的精度,这套活儿早就成熟了。
三个核子这一侧:没有三体散射实验这种装置。你能做的只有两种——要么用三核子系统(氚、³He、中子-氘散射)的束缚能和散射观测量去反推,要么干脆拿原子核和中子星当"巨型探测器",用宏观量去约束微观系数。前者观测量少得可怜,后者一路要跨过十几个数量级。
而且"三体力"这个概念本身就比"两体力"绕。两个核子之间交换一个π介子,这是两体力。三个核子在场的时候,π介子可以被第三个核子"截胡"——两个核子本来要交换的东西,被中间那个吃掉了。1957 年藤田和宫泽就把这种机制写出来了,中间那个可以是个 Δ 共振,所以这条也叫藤田-宫泽力 [9]。
1994 年,van Kolck 从手征拉氏量出发,把核力系统化地推了一遍,三核子力第一次有了"同一套理论里自然长出来"的地位——不用外面塞,它就是理论的产物 [7]。2002 年,Epelbaum 和合作者把 N²LO 阶的三核子力完整算了出来,这才有了可以拿去做计算的显式表达式 [8]。
还有一层麻烦是重整化。三个粒子凑在一块儿,两体那套重整化程序不够用,必须专门加一个三体接触项才能把发散收干净——这件事 1999 年就被 Bedaque、Hammer 和 van Kolck 论证清楚了 [17]。所以三核子力不是"顺手多出来的修正",它是理论自洽性的硬要求。
这里有个容易被忽略的点:三核子力出现得比两核子力晚。手征展开一级一级往上走,两核子力在领头阶就有了,三核子力要到第三级(N²LO)才露面。这不是巧合——从算符构造上,三体接触项需要更多的核子场组合,天然是更高阶的。
算符账本:每一级添几个旋钮
这是全篇最需要钉死的一张表。手征展开每往上走一级,三核子力这边到底添什么东西、添几个自由参数。
| 手征阶 | 三核子力添什么 | 新添旋钮数 |
|---|---|---|
| NLO(第二级) | 什么都不添 | 0 |
| N²LO(第三级) | 两π交换 + 一π交换 + 一个纯接触项 | 2(cD、cE) |
| N³LO(第四级) | 标准说法:不出现新的三体接触项 | 0(有争议,见下) |
| N⁴LO(第五级) | 两导数接触算符 | 13 个独立算符(同位旋投影后 10 个 LEC 组合) |
那两个旋钮是这么分工的。cD 管的是"一π交换加接触项"这一块——它同时挂在一个π介子和一个短程接触上,所以是连接长程与短程的桥。cE 管的是纯粹的核子-核子-核子接触项,完全短程。这两个系数在两核子散射里压根不出现,所以 NN 数据再准也拿不到它们。这就是为什么它们必须"另找来源"。
N³LO 那一行值得多说一句。按标准的幂次计数结论(Bernard、Epelbaum、Krebs、Meißner 在 2008 年和 2011 年的两篇文章里给出的 [10][11]),三核子力到 N³LO 为止,不出现新的三体接触项。这句话的分量在于:它意味着三核子力往上走一级,一个旋钮都不添——贡献完全由已知的东西决定。所以 N³LO 的三核子力是个"纯预言",拿它去算核物质,等于顺便检验幂次计数本身对不对。当然,这个结论并不是所有人都点头,Girlanda 一组就一直在从另一个角度提问。
再往上就是 N⁴LO 的两导数接触项,这里进入了"算符账本"真正复杂的地方。
2011 年,Girlanda、Kievsky、Viviani 三个人做了件苦活 [4]:把两导数的三核子接触拉格朗日量做了最小化——要求同时满足分立对称性、Fierz 关系、庞加莱协变性。结果是 13 个独立的低能常数,还给了它们在坐标空间里的定域形式。
这篇有个后续必须记住:2020 年他们发了一篇勘误。原因是原来的 Fierz 关系算错了几条,改正之后发现还有三个算符 [5]。所以现在说"13 个",用的是勘误后的版本。引用时如果只引 2011 年那篇、不给勘误,等于引了一个错数字。
2023 年,同一组人加上 Marcucci 把这套东西带进了 p-d 散射计算,检验它对散射观测量的影响 [6]。这篇文章之所以重要,是因为核物理里挂了几十年的"分析本领之谜"——p-d 和 n-d 散射在低能下,理论算出的矢量分析本领系统性低于实验——一直被怀疑跟更高阶的三核子力有关。
到了 2026 年,Filandri、Girlanda、Capitani 把这块推进了一步 [3]。他们做的主要是数学整理工作,但整理得很关键:把这 13 个算符按总同位旋投影,发现其中只有 11 个在 T=1/2 通道里有贡献;接着再投到 N-d 散射的分波上,最终归结为 10 个 LEC 的组合。
这件事听着像纯记账,其实它是后面所有事的前置条件。你手上有一堆散射数据想去拟合 N⁴LO 的三体力,得先知道哪些组合对应哪些分波、哪些通道根本用不上——不然就是盲拟合。所以这篇文章的定位很清楚:不是给出系数值,而是把"该拟合什么"这张表格先画出来。
这一整套账本,如果想要一份完整的地图,首选是 Hebeler 2021 年那篇综述 [16]——三核子力怎么实现、怎么用在原子核和稠密物质上,讲得最系统。再往前,手征核力的两篇经典综述是 Epelbaum、Hammer、Meißner 2009 [12] 和 Machleidt、Entem 2011 [13];三体力的普适图像(从冷原子一路到原子核)见 Hammer、Nogga、Schwenk 2013 [14];有效场论整体的现状与展望见 Hammer、König、van Kolck 2020 [15]。
2025 年冒出来的一类,和半年后的反驳
上面那套账本用了几十年,大家默认它差不多齐了。2025 年 7 月,有人往桌上拍了一个新条目。
华盛顿大学核理论研究所(Institute for Nuclear Theory)的四个人——Cirigliano、Dawid、Dekens、Reddy——在《物理评论快报》上发表了一篇文章,标题直白:《三核子力的新一类及其含义》[1]。
他们的主张是这样的。在含π介子的低能有效理论里,还存在一类三核子力,形式是"短程算符把两个π介子和两个核子拴在一起"。关键在于,按最朴素的量纲分析(NDA)去估,这类算符的系数应该很小、可以忽略。但他们认为实际上它被增强了,贡献可能和文献里一直用的那个领头三核子力一样大。
更要命的是后半句:这个力的强度由一个"π介子连到两个核子"的耦合决定,而这个耦合目前实验上没有被约束住。也就是说,如果这个力真的大到那个地步,那过去所有关于中子物质状态方程的结论、所有关于核物质饱和的推断,底下都垫了一个没人量过的自由参数。
这篇文章很短,是快报,对核物质能量和状态方程的影响只做了简短讨论。但它把一个可能被漏掉的结构指了出来。
半年后,有人不服。
2026 年 4 月,波鸿鲁尔大学理论物理二所的五个人——Epelbaum、Gasparyan、Gegelia、Hog、Krebs(Gegelia 还兼第比利斯国立大学)——在《物理评论C》上发了一篇《对三核子力新一类的推敲》[2]。
他们的反驳走的是重整化方案这条路。核心论点是:对应的那些低能常数,本身就是重整化方案依赖的——你换个方案,数值就变。所以单看某个方案里的数值大小、然后说"它比 NDA 预期大",这个推理链是断的。他们的做法是从π介子圈图里把方案依赖的短程成分剔掉,再看这些三核子力还剩多少贡献。
剔完之后的结果:贡献变得相当小,跟 Weinberg 幂次计数的预期一致。他们还把这批力和低阶π交换图里"同类"的贡献放在一起比了比,作为交叉检查。最后的结论是,这些三核子力对中子物质和对称核物质状态方程的改动都不大。
这场分歧现在怎么判?
两边都同意算符存在——这不是算符造得对不对的问题。分歧全在"到底多大",而分歧的根在重整化:一方说增强是真的,另一方说那是方案的人为效应。
这种时候,决定性的证据不会来自争论本身,只会来自第三方的、不依赖同一个方案的计算,或者来自观测——如果这个力真的大到跟领头三体力平级,它总该在某个核的性质或某条天体观测量上留痕。目前两边都还在等这个。
我的读法是:算符是真的,增强与否未决。在有人用不同方案独立算一遍之前,任何一边的结论都只能算一半。
谁给 cD 和 cE 标价
回到那两个"必须另找来源"的旋钮。它们是怎么定下来的?
最早的路子不是用能量,而是用"流"。2009 年,Gazit、Quaglioni、Navrátil 提出一个很漂亮的想法:手征有效理论里,相互作用和流必须来自同一套拉氏量、互相自洽。氚的 β 衰变率——也就是那个 Gamow-Teller 矩阵元——对这种自洽性极度敏感。拿它当约束,cD 就定住了 [18]。
但单个观测量有个老问题:cD 和 cE 的组合是简并的,你沿一条线滑过去,氚结合能都能对上。所以必须有第二个独立约束把它劈开。
2020 年,Hüther、Vobig、Hebeler、Machleidt、Roth 给出了一套后来被广泛采用的做法:拿 ³H 基态能加 ¹⁶O 基态能两个条件一起锁。他们从 Entem–Machleidt–Nosyk 的两核子力出发,配上一套手征阶一致、用非局域正规化子的三核子力,截断取 450、500、550 三个值,把 cD 和 cE 用两个基态能钉住。得到的相互作用能稳定地重现从 p 壳一直到镍同位素链的结合能和电荷半径,把之前那批相互作用的很多毛病都修掉了 [19]。
这套"³H + ¹⁶O"后来变成了行业惯例。
但惯例不等于没争议。2021 年,Wesolowski 一组(Salisbury 大学、查尔姆斯理工大学、俄亥俄州立大学、俄亥俄大学与达姆施塔特工业大学)用贝叶斯方法认认真真查了一遍这件事。他们取四个观测量:³H 结合能、⁴He 结合能、⁴He 电荷半径、³H 的 β 衰变率(Gamow-Teller 矩阵元),用本征矢量延拓做快速模拟,把实验误差、计算方法误差和手征展开被截断的误差全都放进去 [20]。
结论有三条,每一条都值得单独记下来。
第一,把截断误差排除在外,根本找不出一组能同时满足 ³H、⁴He 和 ³H β 衰变率的解。换句话说,过去很多看起来"拟合成功"的做法,其实是在噪声里挑了一个好看的点。
第二,除了 β 衰变率以外,其他观测量给出的约束本质上都是简并的。只有 Gamow-Teller 矩阵元能把简并打开。所以那个不起眼的衰变率,其实是这四个观测量里最不可替代的一个。
第三,他们顺手拟合出这套观测量对应的手征展开参数 Q = 0.33 ± 0.06。这个数告诉你有效展开的收敛速度到底有多快。
"手征展开被截断"这个误差怎么算、怎么写进似然函数,本身是两篇方法学文章立的规矩:Furnstahl 等 2015 [22] 与 Melendez 等 2019 [21]。
这一节的意思是:cD 和 cE 虽然有了公认用法,但它们真正的"值"取决于你用哪几个观测量、以及你用什么口径算截断误差。这两个前提一换,数值就换。所以谈这两个系数时,不交代前提是没法对话的。
一条新探针:让中子从 α 粒子旁边蹭过去
2025 年 10 月,出现了一个思路很干净的新办法。
北京大学核物理与核技术国家重点实验室的杨一龙、孟杰、赵鹏巍,和波鸿鲁尔大学理论物理二所的 Epelbaum、孟璐(孟璐同时在东南大学物理学院),在《物理评论快报》上提出:低能下的中子-α 擦边散射,可以当长程三核子力的探针 [46]。
道理很直接。擦边散射的时候,中子跟 α 粒子之间主要发生作用的距离比较远,短程的东西——接触项、高动量成分——插不进来。于是最长程的那部分三核子力,也就是两π交换那一块,就被单独露出来了。这比在氚里逆推要干净得多,因为氚内部所有成分都缠在一起。
他们用从头算量子蒙特卡洛(QMC),配手征有效理论的两体加三体力、展开到第三级(N²LO),算了中子-α 散射。
结果出来了:中子-α 的 D 波相移,非要有那个最长的两π交换三核子力才描述得了。
这件事的意义在于通道。原来约束长程三核子力,主要靠少数几个观测量,而且彼此之间共享同一套理论假设。现在多了一条物理上完全不同的通道——轻核擦边散射。如果两π交换三体力真有问题,这两条通道应该会给出不一致的答案。这是一种交叉检查。
(这里有个小陷阱顺便记一下:这篇文章里的两位作者都姓孟,但不是一个人。孟杰在北京大学,孟璐在波鸿鲁尔大学与东南大学。二手材料里把两人混成一个的不少见。)
矛盾:定核的旋钮,和定核物质的旋钮,对不上
前面几节讲的是"怎么把旋钮算准"。这一节说的是一个更难受的问题:这些旋钮给不出自洽的答案。
2020 年,爱达荷大学的 Sammarruca 和 Millerson 把这个矛盾写得很直白 [23]。
一边是这样:那批手征两体加三体相互作用,在 N²LO 和 N³LO 都有,它们能给出很好的有限核描述——从 p 壳一直到镍同位素,结合能和半径都对得上。
另一边是这样:同一批相互作用,压不住无限核物质。也就是说,它们不会在正确的密度和能量上出现饱和。核物质不饱和,就意味着能量曲线一直往下走,物理上等于说核物质要塌缩。
反过来呢?如果拟合 cD 和 cE 的时候,除了照例用氚结合能,再加上"核物质必须饱和"这个条件,那么结果掉到另一头:中重核束缚不足,半径系统性地偏大。
两头都错,而且错的量级不是小数点后的事,是几个 MeV 起步。
Sammarruca 和 Millerson 把这个叫"明显的不一致",然后做了一堆测试计算去找原因。他们试的方向包括:是不是局域密度近似这一步引入了偏差、是不是三体力的正规化方式不同造成的、是不是核物质和有限核对外推的敏感度本来就不同。
这篇文章没有给出终局答案。但它把问题钉在墙上:在同一个理论框架里,"能算好核"和"能算好核物质"这两件事,用同一组旋钮办不到。 这是三核子力问题真正的痛处,比"某个系数该取多少"要深一层。
反过来:让中子星给 c₁ 和 c₃ 标价
前面讲的是地面上的定值方式。2025 年 11 月,出来了目前最直接的一次"用天体给三核子力系数标价"。
Somasundaram、Svensson、De、Deneris、Dietz、Landry、Schwenk、Tews 八个人(雪城大学与洛斯阿拉莫斯国家实验室、达姆施塔特工业大学与 EMMI 与马普核物理研究所、洛斯阿拉莫斯理论部与计算统计科学部、多伦多大学理论天体物理研究所与圆周理论物理研究所)在《自然·通讯》上发了一篇,标题就是《从多信使中子星观测推断三核子耦合》[24]。
先纠正一个常见的误记:这篇文章定的不是 cD 和 cE,而是 c₁ 和 c₃。这两个是 π 介子-核子相互作用里的低能常数,它们通过两π交换进入三核子力。c₃ 是其中分量最重的那个,直接决定三核子力有多排斥。所以这是另一条通路,跟上一节讲的 cD、cE 分开。
他们建了一个贝叶斯框架,把天体观测直接接到场论的系数上。用的数据是两块:一是 GW170817 这个双中子星并合的引力波事件(取 LIGO/Virgo 给出的后验),二是 NICER 的三颗脉冲星——J0030+0451、J0740+6620、J0437−4715。
结果摆出来就很清楚了。
实验室那边,πN 散射的 Roy-Steiner 分析给出 c₁ = −0.74 ± 0.02 每 GeV、c₃ = −3.61 ± 0.05 每 GeV(单位按 GeV⁻¹)。(πN 耦合常数本身的精密确定是另一条技术线,Reinert、Krebs、Epelbaum 2021 用它把精度又推高了一档 [50]。)
中子星那边:对 c₁ 几乎没有约束——因为 c₁ 对中子物质状态方程的影响本来就不如 c₃ 大,它的后验基本还是先验的样子。对 c₃,后验中位数落在 −2.52 GeV⁻¹,90% 可信区间是 [−6.15, −0.31] GeV⁻¹。
怎么理解这两组数。它们在 90% 水平上是一致的,这是好消息——天体观测没有跟实验室打架。但天体给出的区间宽出一百倍不止。所以现在的状况是:实验室定得紧、天体定得松,但两边方向不冲突,天体这一路因此拿到了"独立复核"的资格。
还有一个有意思的倾向:观测偏爱不那么排斥的三核子力。原因是 GW170817 和 NICER 的观测都偏向更致密、也就是内部压力更低的星,而更排斥的三核子力会把压力顶上去。这是天体数据在给三核子力"往软的方向拉"。
再往前看一步。他们做了模拟:如果换成下一代引力波天文台(爱因斯坦望远镜、宇宙探险者),用一年内最响的 20 个事件(信噪比大于 225)做层级贝叶斯分析,对 c₃ 的约束精度可以逼近实验室水平。也就是说,三核子力系数这件事,可能会在下一个十年从"只有实验室说了算"变成"两边对等对话"。
这篇文章的分量在于它接的那根线:它把量子场论里的微观耦合常数和几十公里外一个天体的宏观观测量,直接连成了一条可推断的链条。
状态方程这台机器是怎么搭起来的
要说清"EOS 约束 3NF",得先知道这台机器长什么样。它有几层。
最底下是多体方法。2016 年,Drischler、Carbone、Hebeler、Schwenk 做了件标志性的事:第一次把 N³LO 的三核子力放进多体微扰论(算到三阶)和自洽格林函数(SCGF)两套框架,算纯中子物质,还顺带把正规排序框架推广到了有限温度。
为什么要提这一篇:前面说过,N²LO 的三核子力带两个新系数,N³LO 的三核子力原则上不带新系数——它是纯粹预言出来的。所以把 N³LO 加进来是件有风险的事 [29]:万一它对状态方程的贡献大得离谱,那就说明幂次计数出问题了。结果它贡献确实不小,尤其在对称核物质里。这本身就是一条关于幂次计数的信息。
再往上是饱和问题。2019 年,Drischler、Hebeler、Schwenk 把手征相互作用推到 N³LO,专门看核物质饱和能不能对上 [30]。
然后是误差怎么算。2020 年,Drischler、Furnstahl、Melendez、Phillips 用贝叶斯方法处理核物质状态方程 [31],重点是把相关误差——不同密度点之间的误差不是独立的——正确地放进去。同一年,Leonhardt 一组走了另一条路 [32]:低密度用手征 EFT(到 2n₀),高密度用泛函重整化群直接从 QCD 推,两条线在密度上接得上,而且自然冒出一个声速极大值。
2023 年,Keller、Hebeler、Schwenk 造了个高斯过程模拟器,把状态方程扩展到任意质子分数和温度 [33]。
同一条线上的系统比较也是必需的:Marino、Jiang、Novario 2024 把几种图解式从头算方法放在现代手征相互作用下横向比了一遍 [49];Tews 等 2025 把局部手征相互作用推到更大的截断值 [28]。
最新的完整版本是 2026 年 Göttling 一组那篇 [26](达姆施塔特工业大学与 EMMI,Hebeler 和 Schwenk 兼马普核物理研究所)。
他们的做法值得细看,因为它体现了目前这套机器最讲究的地方。他们造了一个二维高斯过程,两个自变量是密度 n 和质子分数 x,训练数据是手征两体加三体力(从领头一直到 N³LO)算出来的多体微扰论结果。用它来评估状态方程和热力学导数,同时把手征 EFT 截断误差带进来。
几个具体数字:
- 均匀物质一直算到两倍饱和密度,n₀ = 0.16 fm⁻³,也就是 0.32 fm⁻³。崩坏尺度取 Λ_b = 600 MeV。
- β 平衡下的质子分数,在 2n₀ 处单调升到约 7.5%。
- 中子星壳层分两段:外壳层接已知的 BPS 状态方程;内壳层用可压缩液滴模型加 Wigner-Seitz 近似(含表面与库仑修正),算中子滴落、质子滴落和两相共存,最后给出与 N³LO 均匀物质自洽的内壳层状态方程。
- 壳-核过渡密度落在 0.062 到 0.088 fm⁻³ 之间。具体位置看状态方程软硬:偏软的落到 0.077–0.088,中等 0.068–0.082,偏硬的 0.058–0.062。
有一点必须说清楚:这篇文章只做到中子星的壳层和外核,没有解 TOV 方程,所以它不报最大质量、不报半径。二手材料里把它当成"给出了中子星最大质量"来引用的,是错读。
同一条线上还有一篇量"费米液体"性质的:Alp、Dietz、Hebeler、Schwenk 在 2025 年算出朗道参数、有效质量、纯中子物质的声速 [27]。费米液体语言的好处是把准粒子之间的相互作用参数化,这些参数是超新星和并合模拟的直接输入。
从外面夹:高密度和实验室各给一道墙
手征 EFT 自己只能管到两倍饱和密度左右。再往上,得靠另外两堵墙从两头夹。
高密度那头是 pQCD。 微扰 QCD 只在 40 倍饱和密度附近才靠得住,离中子星内部差着十几倍。但 2022 年 Komoltsev 和 Kurkela 找到了把它往下压的办法 [36]:不用通常的"压强对能量密度" p(ε),改用压强对化学势 p(μ),再加上热力学稳定性、因果性、以及一个必须成立的整体约束——从低密度到高密度积分出来的压强差,要正好等于两端给定的值。这三条一摆上,40n_s 的信息就能被约束到 2.2n_s 往上,在 5n_s 处至少排除掉原本允许面积的一大半。
2025 年,Finch 一组把从壳层一直到 pQCD 密度的推断做成了统一的非参数框架 [37],也就是不再假设一个函数形式,而是让数据自己长。
实验室那头是重离子碰撞。 2025 年 Cozma 用 dcQMD 输运模型,把理论预测跟 FOPI 的核子流数据、SπRIT 的带电π多重数比数据对比,抽出来的是两个数字 [38]:
- 对称核物质的不可压缩模量 K₀ = 230 加 9 减 11 MeV
- 对称能斜率 L = 63 加 10 减 13 MeV(核子流通道);π 介子通道给 L = 79.9 ± 37.6 MeV
两个物理上独立的通道(核子流 vs 介子产生)给出一致的 L,这是这条路近年最硬的成果。同一年,Cozma 和 Trautmann 把核理论、重离子碰撞、天体观测三头合起来看,把 1.4 倍太阳质量中子星的半径收在 12 到 13 公里之间 [39],并且明确指出剩下的不确定性主要来自饱和密度到两倍饱和密度之间那一段——恰好就是手征 EFT 能管的区间。
SπRIT 的原始数据来自 Estee 等 2021 年那篇(64 位作者)[40]。这条路的总论见 Sorensen 等 2024 年的综述(136 位作者)[41]。
天空那头的三把尺子
天体观测给了三条独立的尺子。
NICER 测的是脉冲星的质量和半径。Miller 一组 2021 年给出 J0740+6620 的半径(28 位作者)[44];2026 年又给出 J0437−4715 的半径(15 位作者)[45],后者是目前半径测量里质量精度最高的一颗,质量 1.418 ± 0.044 倍太阳质量。这颗星的价值在于它质量小、精度高,能卡住状态方程在低质量端的形状。
GW170817 测的是双中子星并合时的潮汐形变。2018 年 LIGO/Virgo 那篇有 1151 位作者 [43]。潮汐形变对状态方程的软硬特别敏感,是中子星半径之外的第二把独立尺子。
对称能这条线可以看 2024 年 Lim 和 Schwenk 的工作。他们用能量密度泛函,配上手征 EFT 的约束、GW170817 和 NICER 的数据,做贝叶斯分析。95% 可信区间给出:对称能 S_v = 30.6 到 33.9 MeV,斜率 L = 43.7 到 70.0 MeV,1.4 倍太阳质量中子星半径 R₁.₄ = 11.6 到 13.2 公里。他们还有一个副产品结论:如果中子星最大质量超过 2.1 倍太阳质量,那么星内部的声速必须超过共形极限——这是对状态方程"有多硬"的一个硬性要求 [25]。
低密度那头的老锚点是 Hebeler 一组 2013 年的工作 [34],以及 Tews 一组 2018 年用状态方程约束声速的工作 [35]。
至于 PREX-II(Yue 一组 2022 年 [42])给出的那条中子皮约束,知识库里两篇专文已经讲透了,这里只当外部输入引用,不展开。
坏消息:尺子自己可能不准
上面所有推断都建立在一个前提上——我们用来算核物质的多体方法是可靠的。2026 年 8 月,这个前提被敲了一下。
北京大学核物理与核技术国家重点实验室的胡荣哲、金绍良、甄鑫、商浩宇、裴俊琛、许甫荣(后三人兼惠州中科院近代物理研究所南方核科学理论中心),加上美因茨约翰内斯古腾堡大学核物理研究所与 PRISMA+ 卓越集群的 Marino,在《物理评论快报》上发表了一篇从头算的严格计算 [48]。
他们用的是全组态相互作用量子蒙特卡洛(FCIQMC),直接算无限核物质。方法本身先做了标定:在可精确求解的 Richardson 配对模型和小模型空间的严格对角化里,用 10⁵ 个 walkers 的 FCIQMC 能重现精确能量,误差比其他方法小几个数量级。
然后是结果。
对称核物质是强关联的,纯中子物质则相对弱关联。 这句话听着平淡,含义不轻。
先看差异有多大。纯中子物质里,各方法每核子能量的差别小于 0.5 MeV——大家基本一致。对称核物质里,FCIQMC 和其他方法能差到 2 MeV。
更要紧的是被丢掉的东西有多少。用较硬的那个相互作用(不含 Δ 的 N²LO 版本),各种多体展开截断丢掉的高阶关联贡献,在饱和密度处约 10 MeV,到两倍饱和密度处约 40 MeV。
40 MeV 是个什么概念。核物质每核子能量的饱和值本身才 −16 MeV 左右。也就是说,在被丢掉的项里,藏着一个比答案本身还大的数字。所有截断方法都系统性地高估饱和密度和饱和能量。
这件事跟本文主题的接口在哪?就在前面那个矛盾上。
cD 和 cE 这两个系数,很多做法是拿"核物质必须饱和"来定的;核物质饱和性质本身,又是用截断的多体方法算出来的。如果这些方法在饱和密度附近就带着 MeV 级的系统偏差,那么"用饱和性质定出来的系数"自然也带着偏差。回过头看,Sammarruca 和 Millerson 指出的"能算好核 vs 能算好核物质"这个两头堵,很可能有一部分根本不来自相互作用——它来自我们算核物质用的方法。
这篇文章把 FCIQMC 定位成无限核物质的严格基准工具。它的用处是把两类误差分开:哈密顿量本身的误差,和多体近似的误差。以前这两团是缠在一起的。分开之后,才有可能判断"到底是力不对,还是方法不对"。
中国这条线上的两件事
第一件是给局部核力配三体伙伴。2026 年 7 月,胡荣哲、李建国、范思钦、许甫荣(北京大学物理学院与核物理与核技术国家重点实验室;李建国与许甫荣兼中科院近代物理研究所与惠州南方核科学理论中心;通讯作者许甫荣)在《物理评论C》上发了一篇 [47]。
他们解决的问题很具体:新出的 Idaho 局部位置空间两核子力已经有了,但它缺一套配套的三核子力——绝大多数现成的三核子力都是动量空间非局域形式,配不上去。他们专门造了一套,正规化子采用"局部加非局域"的混合方案,然后一路算到了 ¹³²Sn 的结合能和电荷半径。
有一个细节值得注意:他们的两个低能常数,用的是 ³H 和 ¹⁶O 的基态能来定的,并且明说这是跟着近期一项工作(就是前面 Hüther 一组 2020 年那套办法)的建议走的。
这个细节的分量在于:它说明"³H + ¹⁶O 锁 cD、cE"已经从一篇论文的做法变成了领域内的标准动作。前面讲的那个矛盾——用这套定法得到的系数配不上核物质饱和——也就跟着被继承下来了。
结论方面,他们认为混合正规化子比纯局部正规化子在同时对上能量和半径上更有优势。
第二件就是前面单独讲过的那条中子-α 探针,北大与波鸿、东南大学的合作。
现在这本账怎么记
把前面所有内容收进一张表。
| 旋钮 / 条目 | 谁来定 | 定住了没有 |
|---|---|---|
| cD、cE(N²LO 短程两个) | ³H + ¹⁶O 基态能;³H 的 β 衰变率(流守恒路线);若干贝叶斯分析 | 用法已标准化,但数值依赖观测量选择与截断误差口径。与核物质饱和之间存在已知不一致 |
| c₁、c₃(πN 两个,经两π交换进三体力) | 实验室:πN 散射的 Roy-Steiner 分析;天体:多信使中子星观测 | 实验室定得紧(c₃ = −3.61 ± 0.05 GeV⁻¹);天体给得松(90% 可信区间 [−6.15, −0.31] GeV⁻¹),但方向独立、两者一致 |
| N³LO 三体力 | 不带新旋钮,纯预言 | 拿来当幂次计数本身的检验 |
| N⁴LO 两导数接触项(13 个算符) | 尚无定论 | 2026 年才完成同位旋投影与分波分解,是拟合的前置工作 |
| "新 3N 类"(两π加两核子的短程) | 理论上存在 | 未决。一派说与领头三体力同级,另一派说剔除方案依赖后贡献很小 |
最后一段判断
三核子力的问题,从来不是"有没有"。是"有多大、谁说了算"。
看这张表,能看出两个方向都在收紧。一头是把算符清单列干净——从 13 个算符,到 11 个在 T=1/2 有贡献,到 10 个 LEC 组合。另一头是找新的、更干净的观测量:中子-α 擦边散射、氚的 β 衰变率、中子星半径,各自从不同角度切进来。
但最要紧的遗留问题,不是表里任何一格的数值。是两套标价体系之间的那道裂缝:让原子核好看的那组 cD、cE,和让核物质饱和的那组,对不上,而且差得不小。
2026 年那个 FCIQMC 结果提供了目前最有希望的一个解释方向——裂缝的一部分可能根本不在相互作用里,而在我们用来算核物质的方法里。如果对称核物质真的像那个结果显示的那么强关联,那么过去所有基于截断多体方法得到的核物质结论,都需要重新过一遍。包括用它们定的系数。
至于"新 3N 类"这场争论,我的判断是:算符是真的,增强与否还没定。在这种争论里,靠同一个重整化方案内部反复计算是分不出胜负的。得等一个不同方案的独立计算,或者等某条观测量上的痕迹。在那之前,两边的话都只能算一半。
---
参考文献
- V. Cirigliano, M. Dawid, W. Dekens, S. Reddy,"New Class of Three-Nucleon Forces and Their Implications",*Physical Review Letters* 135(2), 022501 (2025). DOI 10.1103/kn79-f5m9;arXiv:2411.00097。印刷日期 2025 年 7 月 8 日。4 位作者,全部署名于美国华盛顿大学核理论研究所(Institute for Nuclear Theory, University of Washington, Seattle 98195)。注意 arXiv 初版标题作 "A New Class of Three Nucleon Forces and their Implications",与正式刊名有三处差异(A、连字符、大小写)。
- E. Epelbaum, A. M. Gasparyan, J. Gegelia, D. Hog, H. Krebs,"Scrutiny of the new class of three-nucleon forces",*Physical Review C* 113(4), 044005 (2026). DOI 10.1103/gkjm-f8j9;arXiv:2512.14117。印刷日期 2026 年 4 月 20 日。5 位作者。单位:德国波鸿鲁尔大学物理与天文学院理论物理二所(Institut für Theoretische Physik II, Ruhr-Universität Bochum);J. Gegelia 另兼格鲁吉亚第比利斯国立大学。arXiv 页面未挂 journal-ref,仅凭预印本会误判为未发表。
- E. Filandri, L. Girlanda, Y. Capitani,"Partial-wave decomposition of isospin-projected subleading three-nucleon contact interactions",arXiv:2601.19033 (2026)。3 位作者。单位:E. Filandri——意大利特伦托的欧洲核物理及相关领域理论研究中心(European Centre for Theoretical Studies in Nuclear Physics and Related Areas)与 Bruno Kessler 基金会;L. Girlanda——意大利萨伦托大学"E. De Giorgi"数学与物理系;Y. Capitani——萨伦托大学同系与意大利国家核物理研究院(INFN)莱切分部。截至 2026 年 9 月未见正式期刊发表。
- L. Girlanda, A. Kievsky, M. Viviani,"Subleading contributions to the three-nucleon contact interaction",*Physical Review C* 84(1), 014001 (2011). DOI 10.1103/PhysRevC.84.014001;arXiv:1102.4799。3 位作者。单位:萨伦托大学与 INFN 莱切(Girlanda)、INFN 比萨(Kievsky、Viviani)。
- L. Girlanda, A. Kievsky, M. Viviani,"Erratum: Subleading contributions to the three-nucleon contact interaction [Phys. Rev. C 84, 014001 (2011)]",*Physical Review C* 102(1), 019903 (2020). DOI 10.1103/PhysRevC.102.019903。这篇勘误把原论文的 Fierz 关系修正后多出 3 个算符——引用"13 个算符"必须同时引这一篇,否则数字对应的是修订前的版本。
- L. Girlanda, E. Filandri, A. Kievsky, L. E. Marcucci, M. Viviani,"Effect of the N³LO three-nucleon contact interaction on p-d scattering observables",*Physical Review C* 107(6), L061001 (2023). DOI 10.1103/PhysRevC.107.L061001;arXiv:2302.03468。5 位作者。单位:萨伦托大学与 INFN 莱切、INFN 比萨、比萨大学物理系。
- U. van Kolck,"Few-nucleon forces from chiral Lagrangians",*Physical Review C* 49(6), 2932–2941 (1994). DOI 10.1103/PhysRevC.49.2932。单作者(署名单位:加州理工凯洛格实验室)。
- E. Epelbaum, A. Nogga, W. Glöckle, H. Kamada, Ulf-G. Meißner, H. Witała,"Three-nucleon forces from chiral effective field theory",*Physical Review C* 66(6), 064001 (2002). DOI 10.1103/PhysRevC.66.064001。6 位作者。
- J. Fujita, H. Miyazawa,"Pion Theory of Three-Body Forces",*Progress of Theoretical Physics* 17(3), 360–365 (1957). DOI 10.1143/PTP.17.360。印刷 1957 年 3 月。第一作者全名 Jun-ichi Fujita,第二作者 Hironari Miyazawa。藤田-宫泽力(即三体力中的两π交换经 Δ 共振介导)即出自此文。
- V. Bernard, E. Epelbaum, H. Krebs, Ulf-G. Meißner,"Subleading contributions to the chiral three-nucleon force: Long-range terms",*Physical Review C* 77(6), 064004 (2008). DOI 10.1103/PhysRevC.77.064004。4 位作者。
- V. Bernard, E. Epelbaum, H. Krebs, Ulf-G. Meißner,"Subleading contributions to the chiral three-nucleon force. II. Short-range terms and relativistic corrections",*Physical Review C* 84(5), 054001 (2011). DOI 10.1103/PhysRevC.84.054001。4 位作者。第 10、11 条合起来构成"N³LO 不出现新的三体接触项"这一标准结论的出处。
- E. Epelbaum, H.-W. Hammer, Ulf-G. Meißner,"Modern theory of nuclear forces",*Reviews of Modern Physics* 81(4), 1773–1825 (2009). DOI 10.1103/RevModPhys.81.1773。3 位作者。
- R. Machleidt, D. R. Entem,"Chiral effective field theory and nuclear forces",*Physics Reports* 503(1), 1–75 (2011). DOI 10.1016/j.physrep.2011.02.001。印刷 2011 年 6 月。2 位作者。单位:美国爱达荷大学(Machleidt)、西班牙萨拉曼卡大学 IUFFyM(Entem)。
- H.-W. Hammer, A. Nogga, A. Schwenk,"Colloquium: Three-body forces: From cold atoms to nuclei",*Reviews of Modern Physics* 85(1), 197–217 (2013). DOI 10.1103/RevModPhys.85.197。3 位作者。
- H.-W. Hammer, S. König, U. van Kolck,"Nuclear effective field theory: status and perspectives",*Reviews of Modern Physics* 92(2), 025004 (2020). DOI 10.1103/RevModPhys.92.025004;arXiv:1906.12122。3 位作者。
- K. Hebeler,"Three-nucleon forces: Implementation and applications to atomic nuclei and dense matter",*Physics Reports* 890, 1–116 (2021). DOI 10.1016/j.physrep.2020.08.009;arXiv:2002.09548。印刷 2021 年 1 月。单作者。本主题最完整的综述,是"三核子力怎么实现、怎么用"的首选参考。
- P. F. Bedaque, H.-W. Hammer, U. van Kolck,"Renormalization of the Three-Body System with Short-Range Interactions",*Physical Review Letters* 82(3), 463–467 (1999). DOI 10.1103/PhysRevLett.82.463。3 位作者。三体系统需要三体接触项才能重整化的经典论证。
- D. Gazit, S. Quaglioni, P. Navrátil,"Three-Nucleon Low-Energy Constants from the Consistency of Interactions and Currents in Chiral Effective Field Theory",*Physical Review Letters* 103(10), 102502 (2009). DOI 10.1103/PhysRevLett.103.102502。印刷 2009 年 9 月 1 日。3 位作者。用相互作用的流守恒(氚 β 衰变)定 cD 的开山之作。
- T. Hüther, K. Vobig, K. Hebeler, R. Machleidt, R. Roth,"Family of chiral two- plus three-nucleon interactions for accurate nuclear structure studies",*Physics Letters B* 808, 135651 (2020). DOI 10.1016/j.physletb.2020.135651;arXiv:1911.04955。文章号 135651(该刊用文章号而非页码)。印刷 2020 年 9 月。5 位作者。单位:德国达姆施塔特工业大学核物理研究所(Hüther、Vobig、Hebeler、Roth)、EMMI(Hebeler)、美国爱达荷大学(Machleidt)。"用 ³H 加 ¹⁶O 基态能同时锁定 cD、cE"的出处。
- S. Wesolowski, I. Svensson, A. Ekström, C. Forssén, R. J. Furnstahl, J. A. Melendez, D. R. Phillips,"Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables",*Physical Review C* 104(6), 064001 (2021). DOI 10.1103/PhysRevC.104.064001;arXiv:2104.04441。7 位作者。单位:美国索尔兹伯里大学数学与计算机科学系(Wesolowski)、瑞典查尔姆斯理工大学物理系(Svensson、Ekström、Forssén)、美国俄亥俄州立大学物理系(Furnstahl、Melendez)、美国俄亥俄大学物理与天文系与核与粒子物理研究所,兼德国达姆施塔特工业大学与 EMMI(Phillips)。含截断误差的贝叶斯定值;结论是 ³H 的 β 衰变率(Gamow-Teller 矩阵元)是唯一能打开简并的约束,展开参数 Q = 0.33 ± 0.06。
- J. A. Melendez, R. J. Furnstahl, D. R. Phillips, M. T. Pratola, S. Wesolowski,"Quantifying correlated truncation errors in effective field theory",*Physical Review C* 100(4), 044001 (2019). DOI 10.1103/PhysRevC.100.044001。5 位作者。
- R. J. Furnstahl, N. Klco, D. R. Phillips, S. Wesolowski,"Quantifying truncation errors in effective field theory",*Physical Review C* 92(2), 024005 (2015). DOI 10.1103/PhysRevC.92.024005。4 位作者。
- F. Sammarruca, R. Millerson,"Exploring the relationship between nuclear matter and finite nuclei with chiral two- and three-nucleon forces",*Physical Review C* 102(3), 034313 (2020). DOI 10.1103/PhysRevC.102.034313;arXiv:2005.01958。2 位作者,署名单位:美国爱达荷大学物理系。"能算好有限核的相互作用压不住核物质、加了饱和约束核又欠束缚"这一矛盾的出处。
- R. Somasundaram, I. Svensson, S. De, A. E. Deneris, Y. Dietz, P. Landry, A. Schwenk, I. Tews,"Inferring three-nucleon couplings from multi-messenger neutron-star observations",*Nature Communications* 16, 9819 (2025). DOI 10.1038/s41467-025-64756-6;arXiv:2410.00247。印刷 2025 年 11 月 6 日。8 位作者。单位:美国雪城大学物理系与洛斯阿拉莫斯国家实验室理论部(Somasundaram,并列第一作者)、德国达姆施塔特工业大学与 EMMI 与马普核物理研究所(Svensson,并列第一作者)、洛斯阿拉莫斯理论部(De、Tews)、洛斯阿拉莫斯计算计算与统计科学部(Deneris)、达姆施塔特工业大学与 EMMI(Dietz)、加拿大多伦多大学理论天体物理研究所与圆周理论物理研究所(Landry)、达姆施塔特工业大学与 EMMI 与马普核物理研究所(Schwenk)。约束对象是 c₁ 与 c₃(πN 低能常数),不是 cD、cE。
- Y. Lim, A. Schwenk,"Symmetry energy and neutron star properties constrained by chiral effective field theory calculations",*Physical Review C* 109(3), 035801 (2024). DOI 10.1103/PhysRevC.109.035801;arXiv:2307.04063。印刷 2024 年 3 月 5 日。2 位作者。95% 可信区间:S_v = 30.6–33.9 MeV,L = 43.7–70.0 MeV,R₁.₄ = 11.6–13.2 km。
- H. Göttling, L. Hoff, K. Hebeler, A. Schwenk,"Neutron star crust and outer core equation of state from chiral effective field theory with quantified uncertainties",*Physical Review C* 113(5), 055804 (2026). DOI 10.1103/t9y8-nbf4;arXiv:2512.19593。印刷 2026 年 5 月 11 日。4 位作者。单位:德国达姆施塔特工业大学物理系与 EMMI(全部四人);Hebeler 与 Schwenk 另兼马普核物理研究所。只做壳层与外核(到 2n₀),未解 TOV 方程,故不报最大质量或半径。
- F. Alp, Y. Dietz, K. Hebeler, A. Schwenk,"Equation of state and Fermi liquid properties of dense matter based on chiral effective field theory interactions",*Physical Review C* 112(5), 055802 (2025). DOI 10.1103/ls3l-dn1y;arXiv:2504.18259。4 位作者,单位同第 26 条。注意:该文 DOI 为 2025 年后 APS 启用的随机式格式,按卷期页码拼出的 10.1103/PhysRevC.112.055802 会返回 404。
- I. Tews, R. Somasundaram, D. Lonardoni, H. Göttling, R. Seutin, J. Carlson, S. Gandolfi, K. Hebeler, A. Schwenk,"Neutron matter from local chiral effective field theory interactions at large cutoffs",*Physical Review Research* 7(3), 033024 (2025). DOI 10.1103/r314-6r62。9 位作者。同为随机式 DOI。
- C. Drischler, A. Carbone, K. Hebeler, A. Schwenk,"Neutron matter from chiral two- and three-nucleon calculations up to N³LO",*Physical Review C* 94(5), 054307 (2016). DOI 10.1103/PhysRevC.94.054307。4 位作者。单位:达姆施塔特工业大学核物理研究所与 EMMI(四人)、马普核物理研究所(Schwenk)。首次把 N³LO 三核子力放进多体微扰论与自洽格林函数。
- C. Drischler, K. Hebeler, A. Schwenk,"Chiral Interactions up to Next-to-Next-to-Next-to-Leading Order and Nuclear Saturation",*Physical Review Letters* 122(4), 042501 (2019). DOI 10.1103/PhysRevLett.122.042501;arXiv:1710.08220。3 位作者。
- C. Drischler, R. J. Furnstahl, J. A. Melendez, D. R. Phillips,"How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties",*Physical Review Letters* 125(20), 202702 (2020). DOI 10.1103/PhysRevLett.125.202702;arXiv:2004.07232。4 位作者。
- M. Leonhardt, M. Pospiech, B. Schallmo, J. Braun, C. Drischler, K. Hebeler, A. Schwenk,"Symmetric Nuclear Matter from the Strong Interaction",*Physical Review Letters* 125(14), 142502 (2020). DOI 10.1103/PhysRevLett.125.142502。7 位作者。手征 EFT 低密度与泛函重整化群高密度两条线在密度上接续。
- J. Keller, K. Hebeler, A. Schwenk,"Nuclear Equation of State for Arbitrary Proton Fraction and Temperature Based on Chiral Effective Field Theory and a Gaussian Process Emulator",*Physical Review Letters* 130(7), 072701 (2023). DOI 10.1103/PhysRevLett.130.072701。3 位作者。
- K. Hebeler, J. M. Lattimer, C. J. Pethick, A. Schwenk,"Equation of State and Neutron Star Properties Constrained by Nuclear Physics and Observation",*The Astrophysical Journal* 773(1), 11 (2013). DOI 10.1088/0004-637X/773/1/11。4 位作者。
- I. Tews, J. Carlson, S. Gandolfi, S. Reddy,"Constraining the Speed of Sound inside Neutron Stars with Chiral Effective Field Theory Interactions and Observations",*The Astrophysical Journal* 860(2), 149 (2018). DOI 10.3847/1538-4357/aac267。4 位作者。
- O. Komoltsev, A. Kurkela,"How Perturbative QCD Constrains the Equation of State at Neutron-Star Densities",*Physical Review Letters* 128(20), 202701 (2022). DOI 10.1103/PhysRevLett.128.202701;arXiv:2111.05350。2 位作者。用热力学稳定性、因果性与热力学势的整体匹配,把 pQCD 的 40n_s 结果下压到 2.2n_s 以上;在 5n_s 处至少排除 65% 原本允许的面积。
- E. Finch, I. Legred, K. Chatziioannou, R. Essick, S. Han, P. Landry,"Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities",*Physical Review D* 112(10), 103023 (2025). DOI 10.1103/krc7-kz2l;arXiv:2505.13691。6 位作者。随机式 DOI。
- M. D. Cozma,"The equation of state of nuclear matter above saturation density from intermediate energy heavy-ion collisions",*Il Nuovo Cimento C* 48(2), 66 (2025). DOI 10.1393/NCC/I2025-25066-5。单作者,署名单位:罗马尼亚霍里亚胡卢贝伊国家物理与核工程研究所(IFIN-HH)。该 DOI 未收录于 CrossRef,本条以意大利物理学会期刊目次与 CNR 官方仓储记录核实。给出 K₀ = 230⁺⁹₋₁₁ MeV、L = 63⁺¹⁰₋₁₃ MeV(核子流,FOPI 数据)与 L = 79.9 ± 37.6 MeV(带电π谱比,SπRIT 数据)。
- M. D. Cozma, W. Trautmann,"Neutron star radii from laboratory experiments",*International Journal of Modern Physics E* 34(08), 2530001 (2025). DOI 10.1142/S0218301325300012;arXiv:2505.00390。印刷 2025 年 8 月。2 位作者。把核理论、重离子碰撞与天体观测合并后,1.4 倍太阳质量中子星半径落在 12–13 km。
- J. Estee, W. G. Lynch, C. Y. Tsang, J. Barney, G. Jhang, M. B. Tsang 等 64 人(SπRIT 合作组),"Probing the Symmetry Energy with the Spectral Pion Ratio",*Physical Review Letters* 126(16), 162701 (2021). DOI 10.1103/PhysRevLett.126.162701。64 位作者。装置:理研 RIBF 的 SπRIT 时间投影室。
- A. Sorensen, K. Agarwal, K. W. Brown, Z. Chajęcki, P. Danielewicz, C. Drischler 等 136 人,"Dense nuclear matter equation of state from heavy-ion collisions",*Progress in Particle and Nuclear Physics* 134, 104080 (2024). DOI 10.1016/j.ppnp.2023.104080。印刷 2024 年 1 月。136 位作者。
- T.-G. Yue, L.-W. Chen, Z. Zhang, Y. Zhou,"Constraints on the symmetry energy from PREX-II in the multimessenger era",*Physical Review Research* 4(2), L022054 (2022). DOI 10.1103/PhysRevResearch.4.L022054;arXiv:2102.05267。4 位作者。这是 PRR 不是 PRL 也不是 PRC。
- B. P. Abbott 等 1151 人(LIGO Scientific 与 Virgo 合作组),"GW170817: Measurements of Neutron Star Radii and Equation of State",*Physical Review Letters* 121(16), 161101 (2018). DOI 10.1103/PhysRevLett.121.161101。1151 位作者。
- M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Arzoumanian, K. C. Gendreau 等 28 人,"The Radius of PSR J0740+6620 from NICER and X-Newton Data",*The Astrophysical Journal Letters* 918(2), L28 (2021). DOI 10.3847/2041-8213/ac089b。印刷 2021 年 9 月 1 日。28 位作者。题名中的 XMM-Newton 是欧洲 X 射线空间望远镜,不是这里的计算或模型名。
- M. C. Miller, A. J. Dittmann, I. M. Holt, F. K. Lamb, C. Chirenti, Z. Arzoumanian, J. Berteaud, S. Bogdanov, K. C. Gendreau, W. C. G. Ho, S. M. Morsink, P. S. Ray, R. A. Remillard, Z. Wadiasingh, M. T. Wolff,"The Radius of PSR J0437-4715 from NICER Data",*The Astrophysical Journal Letters* 1000(2), L48 (2026). DOI 10.3847/2041-8213/ae5057;arXiv:2512.08790。在线 2026 年 3 月 27 日,印刷 2026 年 4 月 1 日。15 位作者。该脉冲星质量 1.418 ± 0.044 倍太阳质量。
- Y. Yang (杨一龙), E. Epelbaum, J. Meng (孟杰), L. Meng (孟璐), P. Zhao (赵鹏巍),"Chiral Symmetry and Peripheral Neutron-α Scattering",*Physical Review Letters* 135(17), 172502 (2025). DOI 10.1103/45g7-bmp6。印刷 2025 年 10 月 21 日。5 位作者。单位:北京大学物理学院与核物理与核技术国家重点实验室(杨一龙、孟杰、赵鹏巍)、德国波鸿鲁尔大学物理与天文学院理论物理二所(Epelbaum、孟璐)、东南大学物理学院(孟璐)。两位作者同姓孟但并非同一人。做法:从头算量子蒙特卡洛,手征 EFT 两体加三体力到 N²LO;结论:中子-α 的 D 波相移必须有最长的两π交换三核子力才能描述。
- R. Z. Hu, J. G. Li, S. Q. Fan, F. R. Xu,"Chiral three-nucleon forces for the new local position-space two-nucleon potential in ab initio many-body calculations",*Physical Review C* 114(1), 014313 (2026). DOI 10.1103/gjmd-fyjy;arXiv:2601.08199。印刷 2026 年 7 月 13 日。4 位作者。单位:北京大学物理学院与核物理与核技术国家重点实验室(胡荣哲、范思钦、许甫荣)、中国科学院近代物理研究所(李建国、许甫荣)、南方核科学理论中心(SCNT,惠州;李建国、许甫荣)。通讯作者许甫荣。两个三核子力低能常数用 ³H 与 ¹⁶O 基态能约束,明示沿用第 19 条的建议。计算范围从 ⁴He 到 ¹³²Sn。
- R. Z. Hu, S. L. Jin, X. Zhen, H. Y. Shang, J. C. Pei, F. R. Xu, F. Marino,"Ab Initio Exact Calculation of Strongly Correlated Nucleonic Matter",*Physical Review Letters* 137(6), 062503 (2026). DOI 10.1103/q4p2-r1bt;arXiv:2508.09252。印刷 2026 年 8 月 7 日。7 位作者。单位:北京大学物理学院与核物理与核技术国家重点实验室(胡荣哲、金绍良、甄鑫、商浩宇、裴俊琛、许甫荣)、中科院近代物理研究所南方核科学理论中心(裴俊琛、许甫荣)、德国美因茨约翰内斯古腾堡大学核物理研究所与 PRISMA+ 卓越集群(Marino)。核心结论:对称核物质强关联、纯中子物质弱关联;用不含 Δ 的 N²LO 较硬相互作用时,被截断方法丢掉的高阶关联贡献在饱和密度处约 10 MeV、在两倍饱和密度处约 40 MeV。
- F. Marino, W. G. Jiang, S. J. Novario,"Diagrammatic ab initio methods for infinite nuclear matter with modern chiral interactions",*Physical Review C* 110(5), 054322 (2024). DOI 10.1103/PhysRevC.110.054322。3 位作者。
- P. Reinert, H. Krebs, E. Epelbaum,"Precision Determination of Pion-Nucleon Coupling Constants Using Effective Field Theory",*Physical Review Letters* 126(9), 092501 (2021). DOI 10.1103/PhysRevLett.126.092501。3 位作者。第一作者是 Reinert 不是 Epelbaum——首稿曾把它写成 E. Epelbaum 打头,CrossRef 直查后改正。
---
核实与纠错记录
下表逐条列出本轮核查中发现的错误。既包括中文与英文二手材料里反复出现的说法,也包括本轮起草时凭印象写错、随后查实改正的条目。每一条都以 CrossRef 直查、arXiv 官方 API 反查、期刊官网署名页或装置所属机构官方发布为依据。
| # | 常见说法 / 起草时的写法 | 实际情况 | 依据 |
|---|---|---|---|
| 1 | Epelbaum 等《对三核子力新一类的推敲》只有预印本 arXiv:2512.14117,尚未发表 | 已正式发表:*Phys. Rev. C* 113(4), 044005 (2026),DOI 10.1103/gkjm-f8j9,2026 年 4 月 20 日 | CrossRef DOI 直查;APS 期刊页。其 arXiv 页面未挂 journal-ref,只看预印本会误判 |
| 2 | 核物质那篇多信使推断约束的是 cD 与 cE | 约束的是 c₁ 与 c₃——π 介子-核子低能常数,经两π交换进入三核子力。c₃ 后验中位数 −2.52 GeV⁻¹,90% 可信区间 [−6.15, −0.31] GeV⁻¹;实验室值 c₃ = −3.61 ± 0.05 GeV⁻¹ | *Nat. Commun.* 16, 9819 原文与摘要 |
| 3 | 胡荣哲等 FCIQMC 那篇只是"2026 年 PRL" | 准确条目:*Phys. Rev. Lett.* 137(6), 062503 (2026),DOI 10.1103/q4p2-r1bt,印刷 2026 年 8 月 7 日 | CrossRef DOI 直查与 arXiv journal-ref 双向核对 |
| 4 | Lim 与 Schwenk 的对称能约束还没发表 | 早已发表:*Phys. Rev. C* 109(3), 035801 (2024),DOI 10.1103/PhysRevC.109.035801,印刷 2024 年 3 月 5 日。arXiv 是 2307(2023) | CrossRef DOI 直查。按 arXiv 年份记成 2023 年发表是常见错法 |
| 5 | Yue 等 PREX-II 那篇是 PRL | 是 *Physical Review Research* 4(2), L022054 (2022),DOI 10.1103/PhysRevResearch.4.L022054 | CrossRef 刊名直查。PRL、PRC、PRR 三个缩写只差一两个字母,最易串 |
| 6 | 三核子接触项"13 个算符"出自 Girlanda 2011 那一篇 | 必须同时引 2020 年的勘误(*Phys. Rev. C* 102(1), 019903)。原文的 Fierz 关系有误,修正后多出三个算符,现在通用的"13 个"对应的是勘误后版本 | 勘误标题与 arXiv:1102.4799 修订说明("Three more operators found after correcting some mistaken Fierz relations") |
| 7 | 单靠卷期页码拼 DOI 就能引用 APS 2025 年后的新文章 | 不能。APS 自 2025 年起对新文章启用随机式 DOI。第 27、28、37 条的正确 DOI 分别是 10.1103/ls3l-dn1y、10.1103/r314-6r62、10.1103/krc7-kz2l;按卷期页码拼出的 10.1103/PhysRevC.112.055802 等会返回 404 | CrossRef DOI 直查(拼出来的返回 HTTP 404,随机式可正常解析) |
| 8 | Tews 等 2018 年 ApJ 860, 149 的 DOI 是 10.3847/1538-4357/aac8cb | 实为 10.3847/1538-4357/aac267。aac8cb 不存在 | CrossRef 模糊检索后直查 |
| 9 | Göttling 等那篇给出了中子星最大质量 | 没有。该文只做到壳层与外核(均匀物质上限 2n₀ = 0.32 fm⁻³),未求解 TOV 方程,因此不报最大质量、不报半径。它报的是能量、压强、化学势、β 平衡质子分数(2n₀ 处约 7.5%)与壳层结构 | arXiv:2512.19593 正文 |
| 10 | 中子-α 探针那篇是"孟杰与 Epelbaum 的合作" | 5 位作者:杨一龙(第一作者)、Epelbaum、孟杰、孟璐、赵鹏巍。其中的孟璐不是孟杰——孟璐署名于波鸿鲁尔大学与东南大学,孟杰署名于北京大学 | APS 期刊页署名与机构脚注 |
| 11 | 胡荣哲等那篇 3NF 文章的 arXiv 页面标注期刊是 Physics Letters B | 正式发表为 *Physical Review C* 114(1), 014313 (2026),DOI 10.1103/gjmd-fyjy。CrossRef、arXiv 的 journal-ref 与 INSPIRE 三处一致 | CrossRef 直查 + arXiv API journal_ref + INSPIRE 三源交叉。arXiv 前端渲染的期刊字段不可作为依据 |
| 12 | Cirigliano 等《新一类三核子力》的标题 | 刊名为 "New Class of Three-Nucleon Forces and Their Implications";arXiv 初版作 "A New Class of Three Nucleon Forces and their Implications"(多一个冠词、少一个连字符、F 小写)。引用以所引版本为准 | APS 期刊页 vs arXiv:2411.00097 初版 |
| 13 | 四名作者分属华盛顿大学物理系等不同单位 | 四人全部署名于华盛顿大学核理论研究所(Institute for Nuclear Theory),同一条机构脚注 | APS 期刊页机构脚注 |
| 14 | 局部 3NF 那篇的通讯作者是李建国 | 通讯作者是许甫荣(frxu@pku.edu.cn),署名北京大学与中科院近物所与惠州 SCNT 三家;李建国署名近物所与 SCNT 两家 | arXiv:2601.08199 署名页 |
| 15 | Cozma 那篇的 DOI 可以直接查 CrossRef | 该 DOI(10.1393/NCC/I2025-25066-5)未收录于 CrossRef,直查返回 404。本条以意大利物理学会期刊目次与 CNR 官方仓储记录核实:*Il Nuovo Cimento C* 48(2), 66 (2025) | CNR Solar 仓储记录与 CORE 记录 |
| 16 | 中子-α 探针那篇(第 46 条)用的是"手征 EFT 到 N³LO" | 原文写的是 "up to third expansion order",按手征展开计数就是第三级、即 N²LO(领头 LO 算第一次,NLO 第二次,N²LO 第三次)。不要与 N³LO 混——两者差一级,对应的三核子力内容完全不同 | 论文摘要原文 |
| 17 | πN 耦合常数那篇(第 50 条)的第一作者是 E. Epelbaum | 实为 P. Reinert:P. Reinert, H. Krebs, E. Epelbaum,*Phys. Rev. Lett.* 126(9), 092501 (2021)。本条为首稿凭印象写错、经 CrossRef 直查改正 | CrossRef DOI 直查 |
| 18 | Bernard 一组关于三核子力次领头贡献的两篇文章可以混着引 | 不能,是典型的高危"双胞胎"。第一篇是 *Phys. Rev. C* 77(6), 064004 (2008),讲长程项;第二篇是 *Phys. Rev. C* 84(5), 054001 (2011),讲短程项与相对论修正。两者标题只差副标题,卷号差 7、期号差 1、文章号只差一位(064004 对 054001) | CrossRef 两条分别直查 |