跳到正文
返回核结构
至真探索 · 核结构

从头算走到 A≈100：双开壳、形变，和一支终于画出来的误差棒

2026-09-19

2022 年，一篇发表在 *Nature Physics* 上的论文给铅-208 的中子皮厚度报了一个数：0.14 到 0.20 飞米[28]。同一年，杰斐逊实验室的 PREX-II 实验报的是 0.283 ± 0.071 飞米[30]。两个数不算直接冲突，但也谈不上握手。

有意思的不是差多少，而是报数的方式。PREX-II 是在电子散射装置上"量"出来的；另一边的数，是在超级计算机上"算"出来的——而且算法很抠门：只允许输入"两个核子之间怎么相互作用""三个核子之间怎么相互作用"，不允许输入"铅-208 长什么样"。同位素的名字、质量数、它是不是双幻核、中子皮有多厚，全都是算完之后才出现的输出。

这就是核物理里的"从头算"（ab initio）：从核力出发，解多体薛定谔方程，中间不许偷偷调参数去凑实验值。

这件事听起来像纯粹的方法论洁癖。但它其实是一个赌注：如果你真的能从核力一路推出任何原子核，那么你就能推到那些永远做不出实验的核上去——比如中子星并合里几毫秒内一闪而过的极端缺中子核，那些核的半衰期是毫秒量级，地球上造出来还没装进探测器就没了。

问题是，这个赌注下的筹码是算力，而算力的账本非常难看。

---

一、先说清楚"从头算"在赌什么

核力的第一性原理描述来自手征有效场论（chiral EFT，简称 χEFT）。它的思路不复杂：量子色动力学（QCD）管的是夸克和胶子，但原子核的能量尺度只有 MeV 量级，远远够不着夸克解禁闭的地方。所以在低能下，可以把重得多的核子（质子和中子）和轻得多的 π 介子当成基本的自由度，把夸克胶子的细节"积分掉"，写成一个有效理论[2]。

有效理论的好处是它自带一把尺子。你可以按"阶"（order）做展开：领头阶、次领头阶、再下一阶……每一阶都贡献一个可计算的修正项，而且每一项的系数都由实验数据（主要是核子-核子散射相移）定下来。展开到第四阶（N3LO）的两体核力，精度已经能和唯象高精度势相提并论[3]。

这套相互作用里有两个具体的"零件"值得记住名字。

第一是 EM。它由 Entem 和 Machleidt 在 2003 年给出的 N3LO 两体力[3]，加上一个在 N2LO 截断的三体力，两体力再用相似性重整化群（SRG）演化到一个低分辨率尺度 λ = 1.8 fm⁻¹，三体力截断取 Λ = 2.0 fm⁻¹。合起来写作"EM 1.8/2.0"。这个名字在后面会反复出现，因为很多重核计算都在用它。

第二是带 Δ 同质异能素的相互作用。核子激发成 Δ(1232) 共振态，按手征幂次算其实是"轻"的自由度，不进展开反而不自然。Jiang 等人在 2020 年把 Δ 自由度显式做进了相互作用，得到 ΔNNLOGO：两体散射相移、A = 2 到 4 的束缚态、核物质饱和性质一起拟合，然后一口气把结合能和半径算到 A = 132[5]。另一条类似的线是 NNLOsat，2015 年给出[4]。

这里有个容易搞混的点，值得单独说清楚。

"从头算"不等于"没有参数"。 χEFT 有低能常数，得靠实验定。从头算的意思是：参数定完之后，算任何别的核都不许再动它们。你可以理解成——先用几个轻核把秤砣调准，然后拿这把秤去称所有没见过的东西，不许在称的过程中偷偷挪秤砣。

这个约束一旦立起来，预测就有了分量。也正因为有这个约束，"算得准不准"这个问题才突然变得极其重要：如果误差说不清楚，那从头算和唯象模型就没有本质区别。

---

二、A≈100 这道坎是怎么形成的

先说现在已经走到哪了。

2025 年 9 月，在法国卡昂召开的欧洲核物理会议（EuNPC 2025）上，CEA Paris-Saclay 的 Vittorio Somà 做了一场特邀报告，题目叫《可扩展的从头算方法》[1]。摘要里有一句话，是目前对行业前沿最干净的官方表述：

"在过去十年里，由于多体理论和核力建模两方面的进展，从头算技术稳步推进，现在已能覆盖数十个同位素直到质量数 A≈100 左右，以及若干选定的重核。长期目标是把这类计算最终推广到整个核素图，也就是到质量数 A≈300 为止的数千个核。"

注意这句话的措辞。"数十个同位素"，不是"数千个"；"直到 A≈100"，不是"到 A=200"；"以及若干选定的重核"——"选定"两个字是重点，说明重核那一片不是成片覆盖，是点状突破。

为什么是 A≈100 这个位置？

要理解这个坎，得先理解从头算方法是怎么"骗"过指数爆炸的。

精确对角化（比如无核芯壳模型）的代价随核子数指数增长，A 到 16 左右就吃不消了。所以 A 更大的核必须换思路。主流做法的核心是"选一个参考态，然后做相似变换，把高能自由度和低能物理解耦"。耦合簇（CC）和介质内相似重整化群（IMSRG）都走这条路[6][7][22][23]，计算代价随质量数只是多项式增长，一下就打开了中等质量区。

但这个技巧有一个隐含前提：参考态得选对。

具体分三档：

```第一档 双关壳（doubly closed shell）--------------------------------质子数和中子数都填满壳层，基态几乎是单个斯莱特行列式。参考态天然就对。CC / IMSRG 一路推到 132Sn、208Pb 都没问题。

第二档 单关壳 / 半幻核（semi-magic）--------------------------------一头的壳层是满的，另一头是开的 → 要处理配对。允许粒子数不守恒（Bogoliubov 变换），把配对塞进参考态。代价可控，能横扫整条同位素链。

第三档 双开壳 + 形变（doubly open shell + deformed）--------------------------------两头都开着，基态是一个"球形参考态"描述不了的东西。核子在集体地、相干地往一个方向堆，堆成一个椭球——参考态自己就错了，而且错得很彻底。```

第三档就是那道墙。不是因为算力不够，是因为起点错了。 你可以把参考态想成一张草图：前两档里草图的轮廓基本对，后面只是把细节画细；第三档里草图把圆的画成了方的，后面再怎么描都描不回来。

雪上加霜的是成本。三体力在计算里必须处理，而三体矩阵元的存储量随单粒子基的截断 E3max 急剧膨胀。长期以来能用的上限是 E3max = 18，这个基太小，重核根本收敛不了。

2022 年，Miyagi、Stroberg、Navrátil、Hebeler、Holt 提出了一种新存储方案，只存真正进入 normal-ordered two-body（NO2B）近似的线性组合[9]。内存需求直接砍掉大约两个数量级，E3max 的上限从 18 推到 28。他们用这个方案把 ¹³²Sn 的基态能量算到对 E3max 截断只差几个 MeV——这是"重核的从头算终于能收敛"的标志性一步。

顺便说一句 2019 年的另一条线：Stroberg、Holt、Schwenk、Simonis 用价空间 IMSRG 把滴线预测一路做到了从氦到铁、将近 700 个同位素[8]。注意目标是"滴线"——不是挑几个双幻核报数，是整片区域的稳定边界。那次工作还第一次给滴线预言配上了由实验数据标定的理论误差棒。

从"能算几个核"到"能画一张带误差的地图"，这是两个不同阶段的能力。

---

三、拆墙（一）：让参考态自己变形

既然问题是参考态错了，最直接的办法就是：换个会变形的参考态。

这一思路下面有好几条并行的技术路线，互相之间不是竞争关系，而是各有各的适用面。理解它们的分工，要先区分两种截然不同的关联：

```动力学关联（dynamical correlations）-----------------------------------少数几个核子在高能、短程上互相推搡。典型的：多粒子-多空穴激发。IMSRG 这一族方法抓这个特别有效。对集体运动（整块核一起变形）却几乎无能为力。

静态/集体关联（static / collective correlations）-----------------------------------整个核作为一个整体在形变、配对、转动。典型的：大振幅集体运动。生成坐标方法（GCM）抓这个特别有效。但球形核的低激发态它又对付不了，因为非集体激发根本没有被显式包含进去。```

这两件事是互补的，不是一回事。所以最自然的想法就是：把两套东西拼起来。

路线一：IM-GCM（介质内生成坐标方法）。 用多参考 IMSRG（MR-IMSRG）管动力学关联，用带量子数投影的生成坐标方法（PGCM）管集体关联，拼成一个统一框架。具体做法是：先生成一组不同形变的参考态（形变参数 β₂ 从扁到长，像一串珠子），每个都做粒子数和角动量投影，然后解 Hill-Wheeler-Griffin 方程把它们混合起来，得到既有正确角动量又有形状混合的波函数。

中山大学和密歇根州立大学 FRIB 合作的 Zhou、Ding、Yao、Bally、Hergert、Jiao、Rodríguez 在 *Physics Letters B* 上把这条路用到了镁同位素链[11]。他们看的不是别的，正是"形变共存"这个最刁钻的现象：同一个核里，接近球形的状态和强变形的状态能量几乎一样低。²⁸Mg 到 ³³Mg 一路算下来，他们能把这种共存关系复现出来，而且能追踪 N=20 岛反转（island of inversion）是怎么在 IMSRG 演化过程中自然浮现的。³³Mg 的结果尤其漂亮——基态是强变形的 3/2⁻ 态，而最低的 7/2⁻ 是一个形状同质异能态，由不同 K 值的弱变形组态混合而成。这个核的自旋宇称争论了几十年，从头算能给出一个内部自洽的图像，是件不大不小的事。

路线二：多壳价空间 VS-IMSRG。 价空间 IMSRG 的思路是把整个问题切成"核芯 + 价空间"两块：核芯的复杂性被吸收进一个有效哈密顿量，价空间里则用传统的壳模型代码做精确对角化。因为价空间里是精确对角化，配对和形变天然不被当成障碍——它们只是对角化出来的结果。Miyagi、Stroberg、Holt、Shimizu 在 2020 年把这个方法从"单壳"扩展到"跨大壳"（同时包含 sd 壳和 pf 壳）[10]，直接打进了岛反转区：氖和镁同位素里，N=20 附近的基态确实被侵入组态主导，而中子确实从 sd 壳爬到了 pf 壳。这是从头算第一次在"多重壳"意义上处理这类核。代价是价空间不能无限大——大壳一多，对角化的维数又会重新爆炸。

路线三：Bogoliubov 耦合簇（BCC）。 保留耦合簇的多项式代价，同时让参考态允许粒子数不守恒。Demol、Vernik、Duguet、Tichai 用它系统扫了整条锡同位素链的电荷半径[14]。结论挺扎人：没有任何一个现有精细调节过的 χEFT 相互作用，能同时抓住全部关键特征——绝对半径对了，N=50 到 N=82 之间同位素位移的抛物线形状就不对；抛物线对了，¹³²Sn 那个"拐折"（kink）又不一定出来。换句话说，锡链变成了一个把相互作用放在火上烤的试验台。他们顺带指了条明路：要真正比出高低，得往 ¹⁰⁰Sn 方向和 ¹³⁴Sn 之外去测同位素位移，同时多体一侧得把三体项（triples）和两体电荷密度修正都补上。同一个团队在 2026 年又把锡同位素从滴线到滴线的结合能算了一遍[15]。

路线四：球形耦合簇往开壳扩展。 耦合簇传统上偏爱闭壳参考态。Marino、Bonaiti 等一组人（达姆施塔特、美因茨、橡树岭、萨克莱、鲁汶的合作）用方程运动（EOM）框架把"去掉两个核子"和"加上两个核子"的路径都做了出来，让球形耦合簇能处理子壳层附近的开壳核，并且连电偶极响应函数一起算[16][17]。

路线五：形变自洽格林函数（deformed SCGF）。 这条路最新，也最"重"。Scalesi、Duguet、Somà 在 2026 年 8 月放出的预印本，标题直接写着"标志性的 ²³⁸U"[18]。他们明说了挡在路上的两个障碍：第一，在非常大的单粒子基里做超越平均场的计算，代价高得离谱；第二，在长椭形变很大的地方，平均场能量会出现可疑的塌陷。第一个用高度优化的形变 SCGF 实现解决，第二个靠引入超越平均场之外的关联来治。结果是第一个从第一性原理算出来的铀-238。同一组人还配了一篇讲形变 SCGF 在代数图解构造二阶三阶实现的方法文章[19]。

顺便放一个对照：传统的大尺度壳模型（组态相互作用）在重变形核上早就走得很远。Otsuka、Tsunoda、Abe、Shimizu、Van Duppen 2019 年算 ¹⁵⁴Sm 和 ¹⁶⁶Er，给出的图像是 ¹⁵⁴Sm 里长椭和三轴形状共存、¹⁶⁶Er 有强烈的三轴不稳定性[35]。这不是从头算——它用的是经验壳模型相互作用。但它标出了"从头算要追上的那条线"在哪：你要么给出同样的物理图像，要么解释清楚为什么不同。

---

四、拆墙（二）：滴线那边，连续态不是修正项

前面说的都是束缚核。但核素图的边界不是束缚态结束的地方，而是最后一个还能勉强抓住中子的地方——滴线。

在滴线附近，物理换了一套规则。核子开始"漏"出去，弱束缚态和共振态主导一切，非共振的连续谱也不再是可以忽略的背景。这时候你把波函数展开在一堆束缚态基上，等于用一套圆规去画一个正在散开的墨点。

解决办法是把能量搬到复平面上。Berggren 框架的做法是：在离散的束缚态和共振态之外，沿着复能量平面上的一条围道加一串连续态基矢。这样共振的宽度、连续谱的形状、晕结构的尾巴，都能在同一套基里描述。

北京大学胡柏山、吴齐、孙中豪、徐辅仁在 2019 年把这套复能量基和 IMSRG 接起来，做出 Gamow IMSRG[20]。试金石是碳和氧：中子滴线上的 ²⁴O 共振态被复现出来，已知最重的玻罗米安核 ²²C 的晕结构清清楚楚地显示在密度分布里——而且这个晕能不能出来，关键就在连续谱的 s 波通道。他们顺手还预言了 ²²C 的低位共振激发态。

这件事的意义不在于多算了一个核，而在于：连续态不是"以后再加的小修正"，在滴线那一片它就是主角。 一套从头算方法能不能算到滴线，取决于它有没有从一开始就把连续态放进框架里。

---

五、一支误差棒的三次进化

现在回到开头那两个数：0.14–0.20 和 0.283 ± 0.071。

一支误差棒怎么画，是过去十年核物理理论里变化最快的东西之一。原因很实在：从头算有两层展开，每一层都留下一个"我少算了什么"的缺口。

```第一层：手征展开截断-------------------只算到 N3LO，再往上的阶没算。这一项误差就是你砍掉的那些项。

第二层：多体方法截断-------------------IMSRG(2)、NO2B、价空间大小、基的截断……每一项都留下一个缺口。

第三层：模型标定误差-------------------低能常数是从实验数据拟合出来的，拟合本身有误差，而且参数之间有相关性。

第四层：数值/emulator 误差-------------------用代理模型代替精确计算，代理本身也不精确。```

2015 年之前，这些缺口基本靠"看看相邻阶的变化，估个数"来填。Furnstahl、Klco、Phillips、Wesolowski 在 2015 年给出了第一套明确的贝叶斯处理：把"系数自然大小是多少"这件事写成一个先验分布，然后算出截断误差的置信区间[24]。这一步的关键贡献是：原来那个"看着像就估个数"的操作，被解释成了一种特定的自然性先验。

但 2015 年那套是"逐点"的——每个能量点上各算各的误差，点与点之间没有关联。

2019 年，Melendez、Furnstahl、Phillips、Pratola、Wesolowski 补上了这个缺失[25]。他们用高斯过程（GP）同时编码两件事：系数的自然性，以及系数之间的相关结构。这样一来，误差棒不再是每个点一个孤立的小线段，而是一整条带相关性的曲线；而且模型检验（model checking）变成了一个有力工具——你可以反过来问"数据是否与我的误差模型自洽"。他们把这套方法打包成了公开的 Python 代码 gsum。

第三件武器是特征向量延拓（eigenvector continuation，EC）emulator。

历史匹配这类分析为什么以前做不动？因为要拿几千几万组参数各算一遍多体问题，每组都是一次昂贵的计算。König、Ekström、Hebeler、Lee、Schwenk 在 2020 年给出的办法很巧妙：精确解作为参数的函数，在参数空间里其实躺在一张很低维的曲面上。那就先只算几个精确点，用它们的解张成一个子空间，然后用这个子空间去"延拓"到别的参数值[26]。在 ⁴He 的结合能和半径上验证下来，EC 比高斯过程又快又准。这一条是让"扫 10⁹ 组参数"从梦想变成工程的关键。

第四件是贝叶斯历史匹配（Bayesian history matching）。这个方法不是核物理发明的。它 2014 年出现在宇宙学里：Durham 大学的一群人要用 Galform 星系形成模型拟合观测，模型参数一大堆，跑一次模拟很贵，而且根本不知道似然函数长什么样[27]。

历史匹配对付"不知道似然"的办法很粗暴但有效：不假装知道，直接去砍。对每一组参数，先给一个宽松的先验区间；用 emulator 快速扫一遍；把"就算考虑所有不确定性也跟数据不符"的参数区域整块划掉；在被保留下来的区域里再收紧一轮。如此反复迭代，参数空间像被削铅笔一样一层层削细。

关键在这里：历史匹配不给你"唯一最优参数"，它给你一个"还没被排除的参数集合"。 这个集合的大小，本身就是认知不确定性的度量。这正是从头算需要的东西——因为低能常数本来就不是唯一确定的，而是一个带相关性的分布。

---

六、2022 与 2026：两块里程碑

有了上面四件武器，才有了 2022 年那篇 *Nature Physics*。

第一块：铅-208（2022）。

论文的作者名单是 11 个人：胡柏山、蒋伟光、宫城卓也、孙中豪、Andreas Ekström、Christian Forssén、Gaute Hagen、Jason D. Holt、Thomas Papenbrock、S. Ragnar Stroberg、Ian Vernon[28]。他们分别在橡树岭国家实验室、田纳西大学、查尔姆斯理工大学、达姆施塔特工业大学、TRIUMF、密歇根州立大学和杜伦大学。注意最后那位 Ian Vernon——他就是上面历史匹配方法那篇统计论文的作者，这是两个领域的直接嫁接。

技术上有三件东西缺一不可：

第一，三体力的存储方案。前面提过，Miyagi 等人 2022 年的新存储方案让 E3max 从 18 推到 28[9]，让 ²⁰⁸Pb 这种质量的核在算力上变得可行。

第二，emulator。他们用特征向量延拓加高斯过程，把慢的多体求解器换成快的代理[26]。

第三，历史匹配。他们探索了 10⁹ 组不同的核力参数化，通过若干轮历史匹配，最后只剩 34 组"非不可信"（non-implausible）的参数[27][28]。然后在 ⁴⁸Ca 的每核子结合能、点质子半径、第一个 2⁺ 态能量这三个观测量上做标定加权，得到后验预测分布。

结果就是那个 0.14–0.20 飞米。

这个数比 PREX-II 的 0.283 ± 0.071 飞米小[30]，与弹性质子散射、反质子原子、相干 π 介子光生这几个强子探针的结果一致，也和从双中子星并合引力波推出来的约束一致。把天体物理观测、PREX-II 和 χEFT 约束联合起来做非参数方程态分析，2021 年得到的中子皮是 0.17 ± 0.04 飞米[31]——也在两个数之间。

这篇论文真正的新东西，不是那个数，是那个数后面跟着的完整误差账单。 它第一次说明白了：一个 A=208 的核，怎么从核子-核子散射数据一路推过来，中间每一步的不确定性是多少，最后合成一支什么样的误差棒。

这条线最新的进展出现在 2026 年 6 月。Noël、Heinz、Hoferichter、Miyagi、Schwenk 做了第一个把库仑畸变效应完整纳入的从头算宇称破坏不对称度 A_PV 计算[32]，同时处理 ⁴⁸Ca 和 ²⁰⁸Pb。他们的结论很干脆：从头算偏好的 A_PV，对 ⁴⁸Ca 略小于实验值，对 ²⁰⁸Pb 略大于实验值，全球显著度 1.9σ。然后他们用理论自洽的电荷密度和弱密度，从实验 A_PV 反推出 ²⁰⁸Pb 的中子皮：0.187(25)(18) 飞米——比 PREX-II 原报告的 0.283 小了一大截。

这件事的性质值得说清楚。这不是"理论算错了要改"，也不是"实验错了"。这是理论和实验的误差棒终于精细到了能互相较劲的程度，于是过去被误差掩盖住的 1.9σ 张力浮了出来。对一门学科来说，这是成熟的标志，不是危机。

第二块：¹³²Sn 附近的 70 个核质量，用来改 r-process（2026）。

宇宙里比铁重的一半元素来自快中子俘获过程（r-process）。r-process 的路径沿着极端缺中子的核走，这些核大多做不出实验，只能靠理论质量表。长期以来用的都是纯唯象的全局质量模型——比如有限程液滴模型 FRDM。

Kuske、Miyagi、Arcones、Schwenk 在 2026 年 7 月发在 *Physical Review Letters* 上的论文，第一次把从头算质量真正接进了 r-process 网络计算[29]。作者分属达姆施塔特工业大学物理系、筑波大学计算科学中心、EMMI/GSI、马普核物理研究所和 GSI。

技术细节值得摆出来，因为它展示了这类工作现在是怎么做的：

  • 用价空间 IMSRG 算了 70 个核：Pd 122–134、Ag 123–135、Cd 124–137、In 125–139、Sn 126–140；
  • 相互作用用 NN+3N 1.8/2.0 (EM)，哈密顿量展开在 15 个大谐振子壳里，用 IMSRG++ 代码做 NO2B 截断的演化，价空间内用 KSHELL 做精确对角化；
  • 三体显式算符暂时没上，作者引用近期工作说明这对关联能的影响约 1%，而对分离能的影响更小；
  • 为了估计多体截断误差，他们用了两个不同的价空间（VS1：Z = 28–50、N = 82–126；VS2：Z = 28–50、N = 64–90）。同一个核在两个价空间里算出来的差，就是多体误差的估计——因为如果是精确计算，两个价空间应该给出一模一样的结果；
  • 还换用 ΔN2LOGO(394) 相互作用做交叉检验，分离能结果相似；
  • 质量基线是 FRDM2012 加上 AME2020 的实验值，核反应网络用 WinNet，中子俘获率用 TALYS 2.0 重算；
  • 最后跑了中子星并合模拟的 11 218 条轨迹做质量加权丰度。

结果是什么？

在 N=82 之后，Pd、Ag、Cd 这些核的中子分离能比 FRDM 低，配对关联更强。这直接改变了 r-process 的路径：N=82 的等待点被加强了，物质在那里淤积得更久，整个核合成流变慢，第二峰（A≈130）附近的丰度抬高 2 到 3 倍，而第三峰（A≈195）的位置发生更强的移动。

为什么这件事分量不轻？因为过去二十年里，从头算质量的精度一直在被质疑"到底能不能和唯象质量表比"。这篇论文第一次把这个问题换了个问法：不是"从头算质量比 FRDM 准吗"，而是"把从头算质量放进去，r-process 的预言会不会系统性地变、而且变化方向是否稳健"。他们证明会变，而且这个变化在换用 HFB24、HFBD1M 等不同质量基线、换用不同 β 衰变率集合之后依然成立。

从"理论质量表的一个备选"到"能改变天体物理结论的输入"，这是 2026 年发生的事。

---

七、2025–2026 的新边界：往铅以上、往铀

这两年最有意思的一批结果，都是往质量数更高、形变更强的方向推。

²⁶⁶Pb：超重双幻核。 Bonaiti、Hagen、Papenbrock 在 2025 年 8 月给出 ²⁰⁸Pb 和 ²⁶⁶Pb 的从头算计算[12]。²⁰⁸Pb 用来做框架验证，和实验数据对上了；²⁶⁶Pb 的结论是：它确实是双幻核，最低激发态相对基态有 2.6 MeV 的能隙，而且它就坐在中子滴线上——再多加一个中子，这个核就不束缚了。一个用来描述轻核性质定出来的相互作用，被推到了 82 个质子、184 个中子的地方，还能给出自洽的图像。

²³⁸U：元素周期表的下限。 Scalesi、Duguet、Somà 2026 年 8 月的预印本[18]，前面已经说过。这里补一句它的定位：作者在摘要里写得很直白——"核素图的绝大一部分由重变形双开壳核构成，而这一部分此前完全在射程之外"。铀-238 是第一个从第一性原理算出来的例子，它把"从头算的上端"从铅推到了锕系。

⁸⁰Zr 附近的强变形核。 Hu、Sun、Hagen、Papenbrock 2024 年在 *Physical Review C* 上用耦合簇处理了 ⁸⁰Zr 附近的强变形核[13]。这是和 IM-GCM 路线并行的另一条形变路径。

锡链的系统扫描。 前面提到的 Demol 2026 电荷半径[14]和 Vernik 2026 滴线到滴线结合能[15]，把同一个元素链从头到尾扫了一遍。这种"整条链"的工作比单点突破更有诊断价值，因为它把相互作用的系统性偏差暴露出来了——正如 Demol 他们发现的那样，没有哪个现有相互作用能同时通过全部测试。

¹⁰⁰Sn：双幻核的老标杆。 Morris、Simonis、Stroberg、Stumpf、Hagen、Holt、Jansen、Papenbrock、Roth、Schwenk 在 2018 年给出了 ¹⁰⁰Sn 的从头算[33]。¹⁰⁰Sn 是 N = Z = 50 的最重双幻核，这篇工作确认它确实是双幻的，并且算了它的四极集体性，还基于 ¹⁰⁰Sn 的三粒子-两空穴激发精确算出了 ¹⁰¹Sn 里 7/2⁺ 和 5/2⁺ 之间那个很小的劈裂。同时代的精密质量测量（Mougeot 等，*Nature Physics* 2021）给锡链提供了一批关键实验锚点[34]。

大背景文件。 NuPECC（欧洲核物理合作委员会）2024 年的长程规划在 2024 年 11 月 19 日正式发布，186 页[40]。它给的是整个欧洲核物理的优先级排序，从头算方法的推广和多体方法的可扩展性在其中占的位置不轻。

---

八、还剩三堵墙，以及一张诚实的时间表

把话收拢。按现在的能力和缺口，可以列三堵墙。

第一堵墙：双开壳 + 形变。 现在的情况是"几条路线各有试点"——IM-GCM 在镁同位素链，形变 SCGF 在铀，耦合簇在 ⁸⁰Zr 附近，BCC 在锡链。但还没有一个能像单关壳那样横扫整片区域的标准工具。每换一个质量区，往往就得换方法、换实现、重新做收敛性研究。这不是科学问题没解决，而是工程化程度不够——缺一个"哪里都能用、代价可预期、报错可解释"的通用框架。

第二堵墙：连续态。 Gamow IMSRG 证明了这条路能走[20]，SCGF 一族本来就以处理连续谱见长[21]。但大规模地把它铺到整片滴线区域，和静态关联、形变、三体力同时放进一个框架里——这件事还在路上。物理上，滴线附近恰恰是 r-process 最需要输入的地方，所以这块的缺口是有代价的。

第三堵墙：成本，具体说是三体力。 IMSRG 目前普遍用 NO2B 近似，把三体力在参考态上正规排序后只保留零体、一体、二体项。近期工作把这项近似的误差量化了：对关联能约 1%[29]，对分离能更小，但随质量数增长会不会变大，还在研究。显式三体算符的 IMSRG(3) 和 VS-IMSRG(3) 已经有实现[37]，但系统性地用到中重核上仍然很贵。另外从多参考态出发的 MR-IMSRG(2) 要保留二体和三体不可约密度，这又是另一笔账。

关于时间表，我觉得最诚实的一句话来自 Somà 那份摘要本身："长期目标是……到质量数 A~300 为止的数千个核"[1]。注意"长期"这两个字。

从 A≈100 到 A≈300，不是把现有方法调大一点就能外推的。质量数每往上走一个量级，往往就换一套方法、换一套近似、重新做一遍收敛性研究。而且目标不是"能算出一个数"，是"能给出一个带着可信误差棒的数"——后者的门槛高得多。

还有一个必须说清楚的现实：现在这套体系里的关键节点，大多是"点和线"，不是"面"。 铅-208 是一个点，铀-238 是一个点，锡链是一条线，镁同位素链是一条线。从点到面，中间差的不是算力峰值，是自动化程度、误差传播的规范化和计算成本的下降。

---

九、如果你只记三件事

第一，从头算 = 只输入核力，不输入答案。 参数在轻核上一次性定完，之后算任何核都不许再调。这个自我约束是它全部价值的来源，也是"误差棒"必须画出来的原因。

第二，卡在 A≈100 不是因为算力不够，是因为参考态选错了。 双开壳 + 形变的核，基态根本不是球形——你拿一个球形参考态出发，后面再怎么算都在描一张画错的草图。所有新方法的共同主题，就是想办法让参考态自己会变形、会配对、会混合。

第三，2022 年之后，"误差"从脚注变成了主角。 铅-208 那 34 组参数、锡链那 70 个核质量、²⁰⁸Pb 中子皮那 1.9σ 的张力——它们讲的是同一件事：理论第一次能报出一个"带完整误差账单"的数。这意味着从头算不再只是方法论上的洁癖，而是能和实验互相质证的对手。

至于那两支中子皮的误差棒什么时候能合成一支，现在没人能给出时间。但有一点已经确定：这不再是"谁算得准"的口水仗，而是两套各自带着误差评估的数在互相较劲——这本身就是进步。

---

参考文献

  1. V. Somà，"Scalable ab initio approaches"，特邀报告，European Nuclear Physics Conference 2025（EuNPC 2025），法国卡昂，2025 年 9 月 23 日，Contribution ID 371（会议摘要与演示文稿均可从 indico.in2p3.fr/event/30430 获取）。报告人署名机构为 CEA Paris-Saclay（法国）。本文第一节中"数十个同位素直到质量数 A≈100 左右，以及若干选定的重核"、"长期目标……到质量数 A≈300 为止的数千个核"、"主要挑战在于设计既能处理复杂体系（即双开壳体系）、又能随质量数温和增长的計算方案"三处表述，均直接引自该摘要。同一摘要另于 2025 年 5 月 30 日在韩国基础科学研究院（IBS）同名特邀报告中使用。
  1. E. Epelbaum, H.-W. Hammer, Ulf-G. Meißner，"Modern theory of nuclear forces"，*Reviews of Modern Physics* 81(4), 1773-1825（2009 年 10 月）。DOI: 10.1103/RevModPhys.81.1773；arXiv:0811.1338 [nucl-th]。作者署名机构为于利希研究中心核物理研究所（IKP-3）与于利希强子物理中心、波恩大学辐射与核物理研究所（HISKP）与贝特理论物理中心。手征有效场论推导核力的奠基性综述，本文第一节 χEFT 基础描述的依据。
  1. D. R. Entem, R. Machleidt，"Accurate charge-dependent nucleon-nucleon potential at fourth order of chiral perturbation theory"，*Physical Review C* 68(4), 041001（2003 年 10 月 1 日）。DOI: 10.1103/PhysRevC.68.041001；arXiv:nucl-th/0304018。作者署名机构为西班牙萨拉曼卡大学核物理组与美国爱达荷大学物理系。手征微扰论四阶（N3LO）两体核力的首篇论文，"EM"相互作用的两体部分即出自此文。
  1. A. Ekström, G. R. Jansen, K. A. Wendt, G. Hagen, T. Papenbrock, B. D. Carlsson, C. Forssén, M. Hjorth-Jensen, P. Navrátil, W. Nazarewicz（共 10 位作者），"Accurate nuclear radii and binding energies from a chiral interaction"，*Physical Review C* 91(5), 051301（2015 年）。DOI: 10.1103/PhysRevC.91.051301；勘误见 *Physical Review C* 109(5), 059901（2024）。DOI: 10.1103/PhysRevC.109.059901。即"NNLOsat"相互作用的原始文献，以轻核半径与结合能同时拟合为特色。
  1. W. G. Jiang, A. Ekström, C. Forssén, G. Hagen, G. R. Jansen, T. Papenbrock（共 6 位作者），"Accurate bulk properties of nuclei from A = 2 to ∞ from potentials with Δ isobars"，*Physical Review C* 102(5), 054301（2020 年 11 月）。DOI: 10.1103/PhysRevC.102.054301；arXiv:2006.16774 [nucl-th]。作者署名机构为田纳西大学、橡树岭国家实验室、查尔姆斯理工大学（Jiang 一人署三家）。Δ-full 手征相互作用的代表工作，即 ΔNNLOGO，把结合能与半径算到 A = 132。
  1. H. Hergert, S. K. Bogner, T. D. Morris, A. Schwenk, K. Tsukiyama（共 5 位作者），"The in-medium similarity renormalization group: A novel ab initio method for nuclei"，*Physics Reports* 621, 165-222（2016 年）。DOI: 10.1016/j.physrep.2015.12.007；arXiv:1512.06956 [nucl-th]。IMSRG 方法的标准综述。
  1. H. Hergert，"A Guided Tour of Ab Initio Nuclear Many-Body Theory"，*Frontiers in Physics* 8, 379（2020 年）。DOI: 10.3389/fphy.2020.00379；arXiv:2008.05061 [nucl-th]。本文第二节"参考态分三档""多体展开代价随质量数多项式增长"两处判断的依据；该文亦明确指出"从头算的进一步推广到形变，才能进入双开壳体系"。
  1. S. R. Stroberg, J. D. Holt, A. Schwenk, J. Simonis（共 4 位作者），"Ab Initio Limits of Atomic Nuclei"，*Physical Review Letters* 126(2), 022501（2021 年 1 月 12 日）。DOI: 10.1103/PhysRevLett.126.022501；arXiv:1905.10475 [nucl-th]。用价空间 IMSRG 预测从氦到铁的质子与中子滴线，共约 700 个同位素，并用现有实验数据标定滴线预言的理论误差。
  1. T. Miyagi, S. R. Stroberg, P. Navrátil, K. Hebeler, J. D. Holt（共 5 位作者），"Converged ab initio calculations of heavy nuclei"，*Physical Review C* 105(1), 014302（2022 年 1 月 3 日）。DOI: 10.1103/PhysRevC.105.014302；arXiv:2104.04688 [nucl-th]。作者署马普核物理研究所（Miyagi、Stroberg、Navrátil、Holt 部分署名）与达姆施塔特工业大学核物理研究所（Hebeler）。提出三体力矩阵元的新存储方案，内存降低约两个数量级，使 E3max 从 18 提升到 28，并把 ¹³²Sn 基态能量算到对 E3max 收敛在几个 MeV 以内。本文第二节与第六节的核心技术依据。
  1. T. Miyagi, S. R. Stroberg, J. D. Holt, N. Shimizu（共 4 位作者），"Ab initio multishell valence-space Hamiltonians and the island of inversion"，*Physical Review C* 102(3), 034320（2020 年 9 月 16 日）。DOI: 10.1103/PhysRevC.102.034320；arXiv:2004.12969 [nucl-th]。作者署名机构为 TRIUMF（加拿大）、华盛顿大学物理系、麦吉尔大学物理系、东京大学原子核研究中心。VS-IMSRG 扩展到跨大壳价空间，处理氖与镁同位素中 N=20 岛反转。本文第三节"路线二"的依据。
  1. E. F. Zhou, C. R. Ding, J. M. Yao, B. Bally, H. Hergert, C. F. Jiao, T. R. Rodríguez（共 7 位作者），"Ab initio nuclear shape coexistence and emergence of island of inversion around N = 20"，*Physics Letters B* 865, 139464（2025 年 4 月 9 日）。DOI: 10.1016/j.physletb.2025.139464；arXiv:2410.23113 [nucl-th]。作者署名机构为中山大学物理与天文学院（周、丁、姚、焦）、美国密歇根州立大学 FRIB 实验室（姚、Hergert）、法国 CEA Paris-Saclay 理论核物理与服务部 ESNT（Bally）、西班牙马德里康普顿斯大学物质结构、热物理与电子学系（Rodríguez）。IM-GCM（MR-IMSRG + PGCM）处理镁同位素形状共存；³³Mg 基态为强变形 K^π = 3/2⁻，最低 7/2⁻ 为形状同质异能态。本文第三节"路线一"的依据。
  1. F. Bonaiti, G. Hagen, T. Papenbrock（共 3 位作者），"Structure of the doubly magic nuclei ²⁰⁸Pb and ²⁶⁶Pb from ab initio computations"，arXiv:2508.14217 [nucl-th]（2025 年 8 月 19 日）。DOI: 10.48550/arXiv.2508.14217。（截至本文写作时尚未见期刊发表记录。）作者隶属美国橡树岭国家实验室物理部与田纳西大学。结论：²⁶⁶Pb 为双幻核，最低激发态相对基态能隙 2.6 MeV，该核位于中子滴线上。
  1. B. S. Hu, Z. H. Sun, G. Hagen, T. Papenbrock（共 4 位作者），"Ab initio computations of strongly deformed nuclei near ⁸⁰Zr"，*Physical Review C* 110(1), L011302（2024 年 7 月 24 日）。DOI: 10.1103/PhysRevC.110.L011302。耦合簇方法处理强变形核的另一条路径。
  1. P. Demol, U. Vernik, T. Duguet, A. Tichai（共 4 位作者），"Ab initio calculations of nuclear charge radii across and beyond ¹³²Sn: Putting chiral EFT nuclear interactions to the test"，*Physics Letters B* 878, 140524（2026 年）。DOI: 10.1016/j.physletb.2026.140524；arXiv:2602.22030 [nucl-th]。作者署名机构为鲁汶大学与布鲁塞尔大学跨学科天文与天体物理研究所（Demol）、鲁汶大学、达姆施塔特工业大学与 EMMI（Vernik）、鲁汶大学与 CEA Saclay IRFU（Duguet）、达姆施塔特工业大学、EMMI 与马普核物理研究所（Tichai）。Bogoliubov 耦合簇在 BCCSD 阶计算锡链电荷半径；核心结论是没有任何现有精细调节的 χEFT 相互作用能同时抓住绝对半径、N=50 到 N=82 的抛物线型同位素位移以及 ¹³²Sn 处的拐折。本文第三节"路线三"的依据。
  1. U. Vernik, P. Demol, T. Duguet, A. Tichai，"High-precision ab initio calculations of nuclear binding energies: Tin isotopes from dripline to dripline"，arXiv:2607.05086 [nucl-th]（2026 年）。DOI: 10.48550/arXiv.2607.05086。（截至本文写作时尚未见期刊发表记录。）
  1. F. Marino, F. Bonaiti, P. Demol, S. Bacca, T. Duguet, G. Hagen, G. R. Jansen, T. Papenbrock, A. Tichai（共 9 位作者），"From closed shells to open shells: Coupled-cluster calculations of atomic nuclei"，*Physical Review C* 113(4), 044301（2026 年 4 月）。DOI: 10.1103/p297-y8vq；预印本 arXiv:2512.17311 [nucl-th]。本文第三节"路线四"的依据。
  1. F. Marino, F. Bonaiti, S. Bacca, G. Hagen, G. R. Jansen（共 5 位作者），"Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory"，*Physical Review C* 112(1), 014315（2025 年 7 月）。DOI: 10.1103/PhysRevC.112.014315；arXiv:2504.11012 [nucl-th]。两粒子移除方法的实现细节。
  1. A. Scalesi, T. Duguet, V. Somà（共 3 位作者），"The iconic ²³⁸U: ab initio nuclear structure theory towards the limit of the periodic table"，arXiv:2608.28091 [nucl-th]（2026 年 8 月 28 日提交，8 月 31 日公布）。DOI: 10.48550/arXiv.2608.28091。（截至本文写作时尚未见期刊发表记录。需注意：有第三方新闻稿称此文"发表于 Physical Review Letters"并称作者隶属"Institut de Physique Nucléaire d'Orsay"，与预印本署名不符，本文以预印本署名为准。）作者署名机构为瑞典查尔姆斯理工大学物理与天文系（Scalesi）与法国 Université Paris-Saclay, CEA, IRFU（Duguet、Somà）。首个 ²³⁸U 的从头算计算；指出两大障碍为（i）超大单粒子基中超越平均场计算的代价、（ii）大长椭形变下平均场能量的可疑塌陷。本文第三节"路线五"与第七节的依据。
  1. A. Scalesi, T. Duguet, V. Somà（共 3 位作者），"Deformed self-consistent Green's function method for atomic nuclei at second and third order in the algebraic diagrammatic construction"，arXiv:2608.23700 [nucl-th]（2026 年）。DOI: 10.48550/arXiv.2608.23700。形变 SCGF 的方法学配套文章。
  1. B. S. Hu, Q. Wu, Z. H. Sun, F. R. Xu（共 4 位作者），"Ab initio Gamow in-medium similarity renormalization group with resonance and continuum"，*Physical Review C* 99(6), 061302（2019 年 6 月）。DOI: 10.1103/PhysRevC.99.061302；arXiv:1906.10539 [nucl-th]。作者署名机构为北京大学（第四作者徐辅仁为通讯作者）。在复能量 Berggren 框架下把 IMSRG 扩展到共振与连续态；复现 ²⁴O 的中子滴线共振态，给出 ²²C 的晕结构（连续谱 s 通道起决定作用）。本文第四节的依据。
  1. V. Somà，"Self-consistent Green's function theory for atomic nuclei"，*Frontiers in Physics* 8, 340（2020 年）。DOI: 10.3389/fphy.2020.00340。SCGF 方法用于有限核的标准综述，本文第八节"连续态墙"处 SCGF 定位的依据。
  1. G. Hagen, T. Papenbrock, M. Hjorth-Jensen, D. J. Dean（共 4 位作者），"Coupled-cluster computations of atomic nuclei"，*Reports on Progress in Physics* 77(9), 096302（2014 年）。DOI: 10.1088/0034-4885/77/9/096302；arXiv:1312.7872 [nucl-th]。耦合簇核结构计算的标准综述。
  1. K. Hebeler, J. D. Holt, J. Menéndez, A. Schwenk（共 4 位作者），"Nuclear forces and their impact on neutron-rich nuclei and neutron-rich matter"，*Annual Review of Nuclear and Particle Science* 65, 457-484（2015 年）。DOI: 10.1146/annurev-nucl-102313-025446；arXiv:1508.06893 [nucl-th]。作者署名机构含达姆施塔特工业大学、EMMI、TRIUMF、东京大学。核力对缺中子核与核物质影响的综述，本文第一节与第二节的辅助依据。
  1. R. J. Furnstahl, N. Klco, D. R. Phillips, S. Wesolowski（共 4 位作者），"Quantifying truncation errors in effective field theory"，*Physical Review C* 92(2), 024005（2015 年 8 月 18 日）。DOI: 10.1103/PhysRevC.92.024005；arXiv:1506.01343 [nucl-th]。作者署名机构为俄亥俄州立大学（Furnstahl、Wesolowski）与俄亥俄大学核与粒子物理研究所（Klco、Phillips）。首次把贝叶斯方法系统用于 EFT 截断误差的逐点量化。
  1. J. A. Melendez, R. J. Furnstahl, D. R. Phillips, M. T. Pratola, S. Wesolowski（共 5 位作者），"Quantifying correlated truncation errors in effective field theory"，*Physical Review C* 100(4), 044001（2019 年 10 月 10 日）。DOI: 10.1103/PhysRevC.100.044001。作者署名机构为俄亥俄州立大学物理系（Melendez、Furnstahl）、俄亥俄大学物理与天文系 / 达姆施塔特工业大学核物理研究所 / EMMI-GSI（Phillips 一人署三家）、俄亥俄州立大学统计系（Pratola）、美国索尔兹伯里大学数学与计算机科学系（Wesolowski）。用高斯过程编码 EFT 系数的自然性与相关结构，产出带相关性的截断误差置信带；配套开源代码 gsum。本文第五节的核心依据。
  1. S. König, A. Ekström, K. Hebeler, D. Lee, A. Schwenk（共 5 位作者），"Eigenvector continuation as an efficient and accurate emulator for uncertainty quantification"，*Physics Letters B* 810, 135814（2020 年 11 月 10 日）。DOI: 10.1016/j.physletb.2020.135814；arXiv:1909.08446 [nucl-th]。作者署名机构为达姆施塔特工业大学、EMMI、美国北卡罗来纳州立大学（König）、查尔姆斯理工大学（Ekström）、达姆施塔特工业大学与 EMMI（Hebeler）、密歇根州立大学（Lee）、达姆施塔特工业大学、EMMI 与马普核物理研究所（Schwenk）。以 ⁴He 的基态能量与半径为范例，证明特征向量延拓 emulator 在精度和效率上均优于高斯过程。
  1. I. Vernon, M. Goldstein, R. Bower（共 3 位作者），"Galaxy Formation: Bayesian History Matching for the Observable Universe"，*Statistical Science* 29(1), 81-90（2014 年 2 月）。DOI: 10.1214/12-STS412；arXiv:1405.4976 [stat.ME]。作者署名机构为英国杜伦大学数学科学系（Vernon、Goldstein）与杜伦大学计算宇宙学研究所与物理系（Bower）。贝叶斯历史匹配方法的原始文献，用于 Galform 星系形成模型的参数约束。本文第五节"第四件武器"的依据；该方法后被核物理（²⁰⁸Pb 计算）直接借用。
  1. B. Hu, W. Jiang, T. Miyagi, Z. Sun, A. Ekström, C. Forssén, G. Hagen, J. D. Holt, T. Papenbrock, S. R. Stroberg, I. Vernon（共 11 位作者），"Ab initio predictions link the neutron skin of ²⁰⁸Pb to nuclear forces"，*Nature Physics* 18(10), 1196-1200（2022 年 10 月；在线 2022 年 8 月 22 日）。DOI: 10.1038/s41567-022-01715-8。作者署名机构跨美国橡树岭国家实验室、田纳西大学、瑞典查尔姆斯理工大学、德国达姆施塔特工业大学、加拿大 TRIUMF、美国密歇根州立大学、英国杜伦大学等。用历史匹配探索 10⁹ 组核力参数化，收敛到 34 组非不可信参数；给出 ²⁰⁸Pb 中子皮厚度 0.14–0.20 飞米。本文开头、第一节、第五、第六节的核心文献。
  1. J. Kuske, T. Miyagi, A. Arcones, A. Schwenk（共 4 位作者），"r-Process Nucleosynthesis with Ab Initio Nuclear Masses around the N = 82 Shell Closure"，*Physical Review Letters* 137(5), 052701（2026 年 7 月 28 日）。DOI: 10.1103/k4lv-436c；arXiv:2509.19131 [astro-ph.HE]，初版 2025 年 9 月 23 日，v3 修订于 2026 年 8 月 11 日。作者署名机构为德国达姆施塔特工业大学物理系（Kuske、Arcones、Schwenk）、日本筑波大学计算科学中心（Miyagi）、EMMI / GSI 亥姆霍兹重离子研究中心（Miyagi、Arcones、Schwenk）、马普核物理研究所（Miyagi、Schwenk）。用 VS-IMSRG 计算 Pd 122–134、Ag 123–135、Cd 124–137、In 125–139、Sn 126–140 共 70 个核的质量（NN+3N 1.8/2.0 (EM)，15 个大谐振子壳，IMSRG++ 做 NO2B 演化，KSHELL 精确对角化，VS1/VS2 两价空间之差估计多体误差），接入 WinNet 核反应网络（中子俘获率用 TALYS 2.0 重算，基线质量模型 FRDM2012 + AME2020，并检验 HFB24、HFBD1M），跑了中子星并合模拟的 11 218 条轨迹。结论：N=82 等待点加强，A≈130 附近丰度抬高 2–3 倍，核合成流变慢，第三峰位移。本文第六节的核心文献。
  1. D. Adhikari 等（PREX 合作组，共 101 位作者），"Accurate Determination of the Neutron Skin Thickness of ²⁰⁸Pb through Parity-Violation in Electron Scattering"，*Physical Review Letters* 126(17), 172502（2021 年 4 月 28 日）。DOI: 10.1103/PhysRevLett.126.172502；arXiv:2102.10767 [nucl-ex]。测量宇称破坏不对称度 A_PV = 550 ± 16(stat) ± 8(syst) ppb；联合此前测量给出中子皮 R_n − R_p = 0.283 ± 0.071 飞米。注意：正式标题为 "Accurate Determination…"，不是 "Precision Determination…"，网络二手引用常写错。本文开头与第六节作为对照实验数据。
  1. E. Essick, I. Tews, P. Landry, A. Schwenk（共 4 位作者），"Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of ²⁰⁸Pb with Minimal Modeling Assumptions"，*Physical Review Letters* 127(19), 192701（2021 年 11 月 5 日）。DOI: 10.1103/PhysRevLett.127.192701；arXiv:2102.10074 [astro-ph.HE]。用基于高斯过程的非参数状态方程表示，联合重脉冲星质量、LIGO/Virgo、NICER 观测与 χEFT 约束，得到 S₀ = 33.0(+2.0/−1.8) MeV、L = 53(+14/−15) MeV、中子皮 R_skin(²⁰⁸Pb) = 0.17(+0.04/−0.04) 飞米。
  1. F. Noël, M. Heinz, M. Hoferichter, T. Miyagi, A. Schwenk（共 5 位作者），"Ab initio calculations of parity-violating electron scattering off ⁴⁸Ca and ²⁰⁸Pb"，arXiv:2606.25019 [nucl-th]（2026 年 6 月 23 日）。DOI: 10.48550/arXiv.2606.25019。（截至本文写作时尚未见期刊发表记录。第三方新闻稿称其"发表于 Physical Review Letters"，本文以 arXiv 记录为准。）作者署名机构为瑞士伯尔尼大学（Noël、Hoferichter）、美国橡树岭国家实验室（Heinz）、日本筑波大学（Miyagi）、德国达姆施塔特工业大学 / EMMI / 马普核物理研究所（Schwenk）。首个完整纳入库仑畸变效应的从头算 A_PV 计算；从头算与实验的全球张力为 1.9σ；据此重新提取的 ²⁰⁸Pb 中子皮为 R_n − R_p = 0.187(25)(18) 飞米。本文第六节"这条线的最新进展"的依据。
  1. T. D. Morris, J. Simonis, S. R. Stroberg, C. Stumpf, G. Hagen, J. D. Holt, G. R. Jansen, T. Papenbrock, R. Roth, A. Schwenk（共 10 位作者），"Structure of the Lightest Tin Isotopes"，*Physical Review Letters* 120(15), 152503（2018 年 4 月 12 日）。DOI: 10.1103/PhysRevLett.120.152503。作者署名机构为田纳西大学与橡树岭国家实验室（Morris、Hagen、Papenbrock）、达姆施塔特工业大学核物理研究所（Simonis、Stumpf、Roth）、TRIUMF 与里德学院（Stroberg）、TRIUMF（Holt）、橡树岭国家实验室国家计算科学中心（Jansen）、EMMI-GSI 与马普核物理研究所（Schwenk）。¹⁰⁰Sn 双幻性的从头算确认，¹⁰¹Sn 基于三粒子-两空穴激发的精确计算。
  1. M. Mougeot, D. Atanasov, J. Karthein, R. N. Wolf, P. Ascher, K. Blaum, K. Chrysalidis, G. Hagen, J. D. Holt, W. J. Huang, G. R. Jansen, I. Kulikov, Yu. A. Litvinov, D. Lunney, V. Manea, T. Miyagi, T. Papenbrock, L. Schweikhard, A. Schwenk, T. Steinsberger, S. R. Stroberg, Z. H. Sun, A. Welker, F. Wienholtz, S. G. Wilkins, K. Zuber（共 26 位作者），"Mass measurements of ⁹⁹⁻¹⁰¹In challenge ab initio nuclear theory of the nuclide ¹⁰⁰Sn"，*Nature Physics* 17(10), 1099-1103（2021 年 10 月；在线 2021 年 9 月 23 日）。DOI: 10.1038/s41567-021-01326-9。测量机构以德国格赖夫斯瓦尔德大学物理研究所、马普核物理研究所（Blaum 组）、达姆施塔特工业大学核物理研究所、GSI 亥姆霍兹重离子研究中心、德累斯顿工业大学为主，实验在 CERN 的 ISOLTRAP 装置完成。对 ¹⁰⁰In（¹⁰⁰Sn 的 β 衰变子核）等核做直接质量测量，精度提升约 100 倍，据此挑出 ¹⁰⁰Sn 质量两个互相矛盾的文献值中的一个。本文第七节作为锡链的实验锚点引用。（注意：常被误引的 DOI 后缀 01333-w 是另一篇，正确后缀为 01326-9。）关于 ¹⁰⁰Sn 质量争议的另一方视角，另见 *Physics Today* 新闻稿"A neighboring nuclide narrows down tin-100's mass"（2021 年）——该文说明实验倾向于 GSI 一方的 β 衰变能数值，但 ¹⁰⁰Sn 质量的最终确定仍需更高精度的直接测量。
  1. T. Otsuka, Y. Tsunoda, T. Abe, N. Shimizu, P. Van Duppen（共 5 位作者），"Underlying Structure of Collective Bands and Self-Organization in Quantum Systems"，*Physical Review Letters* 123(22), 222502（2019 年 11 月）。DOI: 10.1103/PhysRevLett.123.222502；arXiv:1907.10759 [nucl-th]。作者署名机构为东京大学物理系、理化学研究所仁科加速器科学中心、比利时鲁汶大学核与辐射物理研究所（Otsuka）、东京大学原子核研究中心（Tsunoda、Abe、Shimizu）、鲁汶大学（Van Duppen）。需注意：这是基于组态相互作用（壳模型）的计算，不是从头算；本文第三节用它标定"从头算需要追上的物理图像参考线"——¹⁵⁴Sm 的长椭/三轴共存与 ¹⁶⁶Er 的强三轴不稳定性。
  1. M. Li, G. C. McLaughlin, R. Surman（共 3 位作者），"Implications of a weakening N = 126 shell closure away from stability for r-process astrophysical conditions"，*Physics Letters B* 872, 140109（2026 年；在线 2025 年 12 月 18 日）。DOI: 10.1016/j.physletb.2025.140109；arXiv:2510.02690 [nucl-th]。作者署名机构为美国加州大学伯克利分校 N3AS 网络与物理系、圣母大学物理与天文系（Li）、北卡罗来纳州立大学（McLaughlin）。指出 N=126 壳闭合在远离 Z=82 处会减弱，而全局质量模型大多没抓住这个趋势；这是下一阶段从头算质量要啃的区域。
  1. M. Heinz, T. Miyagi, S. R. Stroberg, A. Tichai, K. Hebeler, A. Schwenk（共 6 位作者），"Improved structure of calcium isotopes from ab initio calculations"，*Physical Review C* 111(3), 034311（2025 年 3 月 10 日在线）。DOI: 10.1103/PhysRevC.111.034311。作者署名机构为马普核物理研究所（Heinz、Miyagi、Tichai、Hebeler、Schwenk）等。把 IMSRG 从标准截断 IMSRG(2) 推广到含正规排序三体算符的 IMSRG(3)-N7，用它研究 ⁴⁴,⁴⁸,⁵²Ca。两处结论值得记：一是 IMSRG(3)-N7 显著改善了 ⁴⁸Ca 第一个 2⁺ 激发能，从而改善了 N=28 壳闭合的描述；二是对电荷半径的修正并没有解决 ⁵²Ca 与 ⁴⁸Ca 之间那个反常的大半径差。文章还基于方法的尺度广延性给出了 IMSRG(2) 多体误差的估计方法。这是本文第八节"三体力这堵墙"处"显式三体算符已有实现、且多体误差开始可量化"的直接依据。（注意：VS-IMSRG 一侧的近似 triples 修正是另一篇——He & Stroberg 2024；Kuske 等 2026 年用后者做过交叉检验。）
  1. Z. H. Sun, A. Ekström, C. Forssén, G. Hagen, G. R. Jansen, T. Papenbrock（共 6 位作者），"Multiscale Physics of Atomic Nuclei from First Principles"，*Physical Review X* 15, 011028（2025 年 2 月 10 日）。DOI: 10.1103/PhysRevX.15.011028。从第一性原理连接不同尺度核物理的代表性综述。
  1. P. Arthuis, K. Hebeler, A. Schwenk（共 3 位作者），"Neutron-rich nuclei and neutron skins from chiral low-resolution interactions"，*European Physical Journal A* 62(7), 156（2026 年）。DOI: 10.1140/epja/s10050-026-01905-5；arXiv:2401.06675 [nucl-th]。缺中子核与中子皮的另一条计算路径。
  1. NuPECC（Nuclear Physics European Collaboration Committee），"NuPECC Long Range Plan 2024 for European Nuclear Physics"，arXiv:2503.15575（2025 年）。DOI: 10.48550/arXiv.2503.15575，186 页；正式发布于 2024 年 11 月 19 日（经由欧洲科学基金会主办；署名编者包括 Marek Lewitowicz、Eberhard Widmann、Gabriele-Elisabeth Körner）。欧洲核物理的官方优先级文件，本领域长期规划的背景依据。
  1. 本知识库内部交叉引用（用于衔接既有条目，非外部文献）：《核力与核结构》（20260427T162920）、《核有效场论手征微扰论》（20260518T213428）、《核物理第一性原理计算 NCSM-IMSRG-QMC-SCGF》（20260518T213343）、《壳模型与组态相互作用深度研究》（20260529T193424）、《量子多体理论与计算核物理方法》（20260702T141003）、《弱束缚核与连续谱耦合 Gamow 壳模型》（20260816T000713）。

---

*核实说明：本文引用的全部条目均经逐条检索核对作者名单、署名机构、期刊、卷、期、页码与 DOI。核对中发现并已更正的几处常见错误：(1) PREX-II 正式标题为 "Accurate Determination…" 而非 "Precision Determination…"；(2) Hu 等 2022 年摘要中"10⁹ 组参数化"在部分转载页面上因上标丢失被印成"109"，本文按 *Nature Physics* 原文取 10⁹；(3) Otsuka 等 2019 年那篇是组态相互作用（壳模型）计算，不属于从头算，本文已明确标注；(4) Scalesi 等 2026 年 ²³⁸U 一文有第三方新闻稿称其"发表于 PRL"并称作者隶属"Orsay 核物理研究所"，与预印本署名（Chalmers + CEA IRFU Saclay）不符，本文以预印本为准并标注尚未见期刊发表记录；(5) Noël 等 2026 年那篇亦有新闻稿称"发表于 PRL"，同样以 arXiv 记录为准；(6) Mougeot 等 2021 年 *Nature Physics* 一文的 DOI 后缀是 01326-9，网络上流传的 01333-w 属于另一篇，本文已更正；(7) 本文初稿曾把 Heinz 等 2025 年 *Physical Review C* 111, 034311 一文的标题误写为与 N=82 壳闭合相关，实际标题为 "Improved structure of calcium isotopes from ab initio calculations"，内容是用 IMSRG(3)-N7 研究钙同位素并量化多体误差，已按原文更正；(8) Li 等 2026 年一文的分支机构按 arXiv 预印本署名列出（Li 同时署伯克利 N3AS 与圣母大学）。*