一片碎片的中子质子比，怎么变成对称能的一根曲线

先说清楚这篇跟之前几篇的分工

知识库里已经有好几篇跟对称能沾边的文章，各有各的活。

《中子皮厚度与对称能》是入门，《晕核与皮核》讲的是核的边缘能长多远，《一次碰撞，怎么变成一条物态方程约束》讲的是费米能区那批装置——INDRA-FAZIA 怎么凑出来的、SπRIT 往相反方向看什么、RAON 当时走到哪一步。那篇的主角是机器。同一天写成的《碎片会说话》讲的是碎裂本身——磨蚀、消融、碎片的动量分布，以及怎么用动量分布反过来量中子皮，它的主角是碎裂这件事。

这一篇不碰机器。主角是一个观测量：反应产物的同位旋分布，也就是碎片里中子和质子的比例，是怎么被读成对称能的。同一批数据，换成五种不同的读法，能读出五个不完全一样的故事。每一种读法都带着自己的假设，也带着自己的系统误差。把这些摊开看，比看任何一条最终曲线都更有用。

往后半篇换一头：丰中子核自己长什么样、怎么反应。因为所有这些探针，归根结底都建立在一件事上——中子多的核，它的边缘跟普通核不一样。这个「不一样」既是信号，也是麻烦。

一、先把两个数说清楚

核物质有一个理想化的版本：质子中子无限多、均匀铺开、不带电磁效应。它的能量密度随密度怎么变，就是核物质物态方程。里面最不确定的那一块，叫对称能。

它的来历很朴素。核物质偏好质子中子成对，一旦中子多出来，能量就被顶上去。把这一份多出来的能量写成公式，就是

```E/A(ρ, δ) = E/A(ρ, 0) + S(ρ)·δ² + O(δ⁴)

ρ = 核子总数密度δ = (N−Z)/A ，中子过剩的比例S(ρ) = 对称能，随密度变化```

奇数次幂被同位旋对称性禁掉了，δ⁴ 项很小可以忽略。所以对称能就是一个只跟密度有关的函数 S(ρ)。

真正好用的写法，是在饱和密度 ρ₀ 附近（约 0.16 个核子每立方飞米）做泰勒展开：

```S(ρ) ≈ S₀ + L·x + ½·K_sym·x² + ...

x = (ρ − ρ₀) / (3ρ₀)S₀ = 饱和密度处的对称能大小L = 斜率，管「稍微压缩一点，对称能涨多快」K_sym = 曲率```

S₀ 相对好办，三十年前从质量公式就能定出三十兆电子伏上下。难的是 L。

L 难在哪？它管的是压力——中子物质和质子物质的压力差。这个东西看不见摸不着，只能间接反推。而且它不是一个数字，是一条曲线在某个点上的导数。不同实验、不同理论给出 L 的范围能差出好几倍。

最近的几组数放在一起看更清楚。中山大学张振与上海交通大学陈立文用贝叶斯方法把 CREX 和 PREX-2 的数据合起来，得到 S(ρ₀) = 30.2 +4.1/−3.0 MeV，而 L = 15.3 +46.8/−41.5 MeV（90% 置信度）[19]。东德州农工大学的 William G. Newton 与 Gabriel Crocombe 用中子皮加纯中子物质先验，得到 J = 31.3 +4.2/−5.9 MeV、L = 40 +34/−26 MeV[18]。孙向东等人的《中国科学：物理学 力学 天文学》那篇联立重离子碰撞与中子星观测，给的是 L₀ = 64.59 +13.25/−9.38 MeV[28]。

三个数，误差棒覆盖的范围从十几到一百多兆电子伏。这就是现状：中心值差两倍以上，误差棒宽得能把对方整个吞掉。

二、中子多的核，凭什么当探针

要在实验室里碰到对称能，得找一个「中子过剩梯度」很大的地方。丰中子核正好就是。

2.1 中子皮：多出来的一层

一个铅-208 核里有 82 个质子、126 个中子。多出来的 44 个中子挤不进质子的地盘，被推到外面，形成薄薄一层。这层的厚度就是中子皮，记作 R_n − R_p。

它为什么跟对称能有关？因为把多余中子往外推的力，正是对称能提供的压力。皮越厚，说明这股压力越大，对应的 L 一般也就越大。Tsang 等 2012 年那篇 19 位作者的综述[26]把当时所有实验与理论的中子皮约束摆在一张表上比过，是这个「一般」最早的出处之一。

这个「一般」不能省。不同核、不同探针给出的趋势并不完全同步——后面第六节会看到，铅-208 和钙-48 给出的「软硬」方向就不一致。

2.2 晕：比皮更极端

皮是「多出来一层」，晕是「多出来一大团」。

锂-11 是最有名的例子。它有一个锂-9 的核心，外面挂着两个中子，飘得极远。判断标准是：价核子有一半以上的概率待在核心势场之外。这类系统通常还是「玻罗米安」结构——三个部分（核心加两个中子）合起来是束缚的，但任何两个部分单独拿出来都不束缚。

北京大学与理化学研究所的合作者在 2026 年 3 月的《Particles》上发了一篇综述，把轻丰中子核的近期进展系统过了一遍[30]。视野再宽一档的，是叶沿林、杨晓菲、樱井博仪与胡柏山四位在《自然综述·物理》上的那篇《奇特核物理》[29]，覆盖从滴线到壳演化的整个方向。里面几条值得记：硼-17 这个「弱晕核」的 s 波成分被证明是极少的；锂-11、铍-14、硼-17 里都测到了空间局域的氘核子关联；四中子体系有一个共振候选，而三中子体系里没有找到共振态；铍-16 的直接双中子衰变被刻画出来了；氦-8 的 0₂⁺ 态表现出类似「α 加两个中子对再加两个中子对」的凝聚结构。

2.3 为什么这两种东西对外面的反应特别敏感

碰撞发生的时候，两个核先接触的是表面。中子皮和晕恰好全在表面。所以核子输运、碎片形成、轻粒子的发射，第一现场都在这里。这就是为什么丰中子核诱导的反应特别适合当探针。

同时也是麻烦。表面到底有多弥散、有没有团簇、核是不是变形——这几件事都会改变表面的样子，跟中子皮的效果混在一起。2026 年 9 月复旦大学蔡宝军、方德清、马余刚三位在 arXiv 上挂了一篇 59 页的专题综述，专门讲中子皮怎么影响重离子碰撞中的粒子发射，里面把这一串「纠缠不清的因素」列得很细：形变、表面弥散、壳结构、团簇、模型依赖[16]。这是目前这个方向最完整的一份地图。

三、碰撞那一瞬间，中子往哪跑

现在说输运。两个中子过剩程度不同的核撞在一起，接触的那一小块区域里，有两个梯度同时在推核子：

第一个是密度梯度。 两个核的费米尾巴叠在一起，接触区密度高，周围低。核子会从高密度往低密度跑。

第二个是同位旋梯度。 一个核中子多、一个核中子少，接触区里中子的比例不均匀，中子会从「中子多的一侧」往「中子少的一侧」扩散。

第二个梯度靠什么驱动？靠对称能的压力。所以，如果你能测出中子到底跑了多少，就能反推这股压力有多大——也就是 L。

麻烦在于这两个梯度是同时作用的，出来的效果混在一起。而且碰撞过程还有预平衡发射、碎片形成、激发态退激这些后一步的事，它们会改写最终测到的中子质子比。所以「测到的分布」不等于「输运那一刻的分布」。

所有的读法，本质上都是在想办法绕开这层污染。

四、五种读法

4.1 同位旋扩散比：最老牌的做法

第一种做法很聪明：不直接测「跑了多少」，而是测「跑完之后还剩多少差异」。

具体是用三个反应系统。设 a = 锡-124（中子多），b = 锡-112（中子少）。做 a+a、b+b、a+b 三种碰撞，能量和碰撞参数都一样。然后定义

``R_i = ( 2·X^{a+b} − X^{a+a} − X^{b+b} ) / ( X^{a+a} − X^{b+b} )``

X 是某个对同位旋敏感的观测量，比如弹核残余的平均同位旋不对称度 δ。

这个式子的构造很有讲究。a+a 和 b+b 是对称碰撞，两边中子比例一样，不会发生同位旋扩散。它们的作用是当标尺，把预平衡发射、二次衰变这些跟同位旋无关的效应约掉。剩下 a+b 相对这两把标尺的偏差，才归给扩散。

读法：R_i = ±1 表示完全没扩散；R_i ≈ 0 表示完全混匀。

这个量最早的思想来自 FOPI 合作组 2000 年的一篇论文，70 位作者，用 ⁹⁶Ru + ⁹⁶Zr 在每核子 400 兆电子伏下做「同位旋示踪」[2]。真正把它接到对称能上的是 Tsang 等人 2004 年在《物理评论快报》上的工作，21 位作者，用 ¹²⁴Sn + ¹¹²Sn 在每核子 50 兆电子伏、碰撞参数 6 飞米[1]。

Tsang 那篇算下来，扩散主要发生在碰撞后约 30 到 80 飞米每光速那段时间，对应的中心密度从 1.2ρ₀ 掉到 0.3ρ₀。用自由空间核子核子截面拟合，得到 S(ρ) ≈ 31.6·(ρ/ρ₀)^1.05；换成介质中的截面，指数变成约 0.7[1]。

R_i 这条线后来被 INDRA-FAZIA 接手。这个合作组从 2022 年第一次出结果开始，一路做下来：PRC 106, 024603 (2022) 是装置首次结果，47 位作者[11]；PRC 108, 054611 (2023) 转向准弹破裂通道，42 位作者[12]；2025 年 4 月的 PRC 111, 044601 给出「模型无关」的镍镍体系同位旋扩散测量，42 位作者[13]；同年 9 月的 PLB 868, 139815 给出对称能约束，34 位作者[14]。到 2026 年 7 月，S. Mallik 等四人在《Universe》上把这条线做了一次完整综述，把同位旋输运比和同位旋扩散流两条路都摆出来了[15]。

R_i 的优点很清楚：它把模型无关的部分做到了最大——标尺系统自己约掉了污染。缺点也清楚：它只告诉你扩散的总量，不直接给出「在哪个密度上」。那段密度是模型算出来的，不是量出来的。

4.2 同位旋标度：从产额比里读数

第二种做法换了个角度：不看某一个碎片，看很多碎片的产额比值。

两个体系（中子多的和中子少的前身核）碎裂之后，同一个碎片 (N, Z) 在两边的产额之比，往往满足一条指数关系：

``Y₂(N, Z) / Y₁(N, Z) = C · exp( α·N + β·Z )``

α 和 β 就是「同位旋标度参数」。这条规律的漂亮之处在于，α 大致正比于两个前身核的同位旋不对称度平方之差：

``α ≈ 4 · (C_sym / T) · Δ[(Z_s / A_s)²]``

C_sym 是对称能系数，T 是温度，Z_s/A_s 是前身核的质子分数。所以只要测出产额、量出温度，就能反推 C_sym。

这个规律的原始论文是 Tsang 等 2001 年在《物理评论快报》上的《核反应中的同位素标度》，6 位作者[3]。同年《物理评论快报》上还有一篇讲同位旋分馏的，20 位作者，讲的是碎裂过程中中子和质子往不同碎片里分配[4]。

这条路 2025 年有了新数据。印度可变能量回旋加速器中心（VECC，加尔各答）的 Payel Karmakar 等 18 位作者，用 K500 超导回旋加速器的束流，把 ¹⁴N 和 ²⁰Ne 分别打到 ¹¹²Sn、¹¹⁶Sn、¹²⁴Sn 上，束流能量每核子 18 到 30 兆电子伏[5]。在所有体系里都清楚看到了同位旋标度行为。温度用双同位素比法测出来，最后给出 C_sym 约从 24 掉到 17 兆电子伏，对应的激发能每核子 2.1 到 2.8 兆电子伏。

他们同时发现一个反例：在激发能每核子 1.8 兆电子伏处，C_sym 的值跟不上下降趋势。作者自己的解释是，那个能量下多重碎裂可能正在消失，机制换了，所以这个点不能跟其它点放在一条线上看[5]。这种「作者主动标出体系失效点」的做法，读起来很舒服。

4.3 等压产额比：把 a_sym/T 直接读出来

第三种做法是中国人做得比较多的一条路，代表是河南师范大学马春旺组。

思路是这样的：取同一个质量数 A 的一组同量异位素（isobar，就是 A 相同、Z 不同的那些核），看它们的产额比怎么随 Z 变化。用修正的 Fisher 模型推下来，产额比的对数与一堆能量系数成简单线性关系：

``ln R ~ 线性组合( (μ_n − μ_p)/T, a_c/T, a_sym/T, a_p/T )``

其中 a_sym 就是 Weizsäcker 半经验质量公式里的对称能项，a_c 是库仑项，a_p 是对能项，T 是温度。也就是说，只要量一组同量异位素的产额，就能把这几个系数一起解出来。

马春旺等人 2011 年在《物理评论 C》上把这个框架立起来，4 位作者，用的是 ⁴⁸Ca 和 ⁶⁴Ni 弹核碎裂的碎片[6]。署名单位包括河南师范大学物理系与中科院上海应用物理研究所核物理研究室。

2013 年他们又往前推了一步，8 位作者，用统计磨蚀-烧蚀模型算 ⁶⁰Ni + ¹²C 在每核子 80 和 140 兆电子伏下的碎片产额[7]。这篇里几个结论值得记：

• 入射能量对 a_sym/T 的影响很小。
• 中子过剩 I = N − Z 相同的碎片，a_sym/T 随 A 增大而上升；而同量异位素之间，a_sym/T 随 A 增大而下降。两个趋势方向相反，很容易写混。
• 前身碎片的 a_sym/T 比最终碎片小不少——说明这条路上二次衰变的污染是实的，不是可忽略的。
• 碰撞参数小的碎片 a_sym/T 比碰撞参数小的大，作者解释为温度更高时 a_sym/T 下降。
• 参考 IYR 的选法会影响结果，对中子更富的碎片影响尤其大[7]——这一条是方法本身的软肋，作者自己写出来了。

这条路 2026 年有了机器学习的版本，后面 4.5 节讲。

4.4 轻团簇比：对「杂音」最不敏感的一种

第四种做法不看重碎片，看最轻的那几个：氚（t，一个质子两个中子）和氦-3（³He，两个质子一个中子）。它们的产额比 t/³He 也是一个同位旋指标。

这个量的好处非常具体。陈、柯、李、雍几位在早期的工作里算过：软对称能和硬对称能给出的 t/³He 比相差约 17%（在 ⁵²Ca + ⁴⁸Ca 每核子 80 兆电子伏的中心碰撞里）。而在同一个计算里，把不可压缩系数从 380 兆电子伏改成 201 兆电子伏、或者把核子核子截面从自由空间值换成介质中值，t/³He 比的变化不到 4%[43]。

换句话说，它对对称能敏感，对同位旋标量那部分几乎不敏感。这正是你要的：一个干净的同位旋探针，杂音越小越好。

还有一个动力学细节很好玩：t/³He 比随团簇动能的变化方向，软硬两种对称能恰好相反——软对称能下它随动能上升，硬对称能下它随动能下降。所以真正有信息的是高动能那一端。

这条线后来被雍高等人推到更高密度，用 ¹³²Sn + ¹²⁴Sn 在每核子 400 兆电子伏下算 t-³He 的相对流与微分流[44]。他们的结论是：在这个反应里，对称能对 t-³He 相对流的影响比对 π⁻/π⁺ 比的影响还强。

4.5 机器学习与贝叶斯：2026 年的新面孔

前面四种读法都是「先算一个观测量，再跟模型比」。机器学习换了个玩法：直接把整个分布交给网络，让它学「什么样的分布对应什么样的对称能」。

2026 年 2 月，河南师范大学的彭丹与马春旺在《物理评论 C》上发了一篇，2 位作者，用贝叶斯神经网络做同位旋标度分析[8]。数据是质子轰击钙、镍、氙同位素靶产生的散裂碎片。他们在不同体系和不同碎片组之间，看到同位旋标度参数 α 与源不对称度差 Δ[(Z_s/A_s)²] 有很强的线性关系——这本身支持了「碎片形成时确实处于统计平衡」这个假设。然后用这个关系反推 Z = 4 到 26 的碎片的 C_sym。

结论里最有价值的是限制条件：C_sym 对碎片的选择、温度、反应体系都很敏感，而且重靶（比如氙）上这个方法的表现明显变差[8]。作者没有把方法吹成万能，这点是诚实的。

同年 6 月，西南大学杨祖星、樊晓华与扬州大学陈鹏辉、理化学研究所西村俊二（4 位作者）在《物理评论 C》上发表了另一个方向的机器学习工作[9]。他们用同位旋相关的 Boltzmann-Uehling-Uhlenbeck 输运模型模拟 ²³⁸U + ²³⁸U 在每核子 500 兆电子伏下的极中心碰撞，然后训练一个「卷积取向滤波器」网络来判断两个变形铀核的相对取向。只用末态粒子的各向异性流，验证集准确率约 80%。这个工作的意义在于：变形核碰撞里，反应取向本身是个混淆变量，网络能在一定程度上把它分出来，这是评估对称能效应可靠性的前提。

更早的一个代表性结果是湖州师范学院王永佳组的深度卷积网络工作，用 UrQMD 模拟的质子与中子横动量-快度二维分布训练。他们的结论很具体：逐事件用质子谱分类软硬对称能，准确率只有约 60%（涨落太大）；改成事件求和的质子谱，准确率升到 98%。回归任务里，预测 L 的平均绝对误差在质子谱上是约 20.4 兆电子伏，中子谱上是约 14.8 兆电子伏[10]。

这条线要提醒一句：网络的输出精度不等于物理精度。它只是把「模型内部区分的难易」翻译成了数字。如果训练用的输运模型本身有偏，网络会把这个偏差学得很扎实。

五、探针的另一头：丰中子核自己怎么反应

前面所有读法都建立在丰中子核的反应上。那丰中子核自己有什么特别的反应特征？

5.1 反应截面：最直接的招牌

晕核的经典招牌就是反应截面偏大。1980 年代中期，谷畑勇夫等人用高能放射性束测相互作用截面，发现好几个丰中子核的物质半径远超 R ∝ A^(1/3) 的液滴公式预期。锂-11 最夸张。

这个数字到今天还在用。2026 年 3 月那篇《Particles》综述里，把这类截面测量列为晕核判据的第一条[30]。

5.2 电荷拾取：丰中子核的指数式增长

这是一个 2025 年的新结果，而且反差极大。

北京航空航天大学孙保华与谷畑勇夫领导的团队，用德国重离子研究中心（GSI）的装置，把 24 种 p 壳层原子核（含大量丰中子核）以 86% 光速轰击氢靶和碳靶，测「电荷拾取」总截面——就是一个中子翻转成质子的过程，它跟 β 衰变直接相关。

结果是：截面随弹核中子数增加呈指数式增长。从稳定的碳-12 到丰中子的碳-19，截面涨了约 50 倍，而基于稳定核建立的经验公式只预测 2.5 倍[34]。

作者的解释是，丰中子核可以打开更多「中子到质子」的超允许跃迁通道。这个结果不只是核反应的事——它直接连到长期困扰核物理界的 β 衰变强度猝灭问题上。

这篇发表在《物理评论 X》15(3), 031004 (2025)，37 位作者，第一作者张寂潮（北航博士毕业生，现为大阪大学博士后），共同通讯作者孙保华与谷畑勇夫，合作单位包括 GSI、加拿大 TRIUMF、大阪大学等十所机构。

5.3 转移与敲出：读单粒子结构

单核子转移反应（敲出一个质子、加进一个中子之类）是读单粒子结构的老工具。它给出三个东西：谱因子、自旋宇称、有效单粒子能。

2026 年 5 月，北京师范大学刘红娜在《Frontiers in Physics》上发了一篇小综述，专门讲用敲出反应追踪幻数的迁移[31]。主线是两个：N = 20 和 N = 28 的幻数在丰中子核里失效；N = 32 和 34 出现了新的壳隙。她也点出了观测量的链条：核质量、电荷半径、2₁⁺ 能量、B(E2) 跃迁概率、直接反应的谱因子——这几个互补的量合起来才能给出「壳闭不闭」的一致图像。

同一个逻辑在《中国科学通报》2026 年的一篇综述里被写得更多细节，署名北京大学[33]。里面两条具体数据：氘核打铍-11 生成铍-12，测到 d 波侵入组态占 57%，跟铍-11 里 s 波主导（约 80%）形成鲜明对比；氘核打碳-16 生成碳-17，给出了 s₁/₂ 与 d₅/₂ 轨道简并的证据，支持碳同位素里 N = 14 壳闭的弱化。同一条思路上还有理论侧的新提案：中西勇纪、田中纯贵、民井淳与远藤慎平四人 2026 年 4 月提出用低能背角 (d,p) 反应，靠库仑位垒把反应局限在核外，专门去探近阈 s 波成分的渐近尾巴[39]。他们还建了一套大接收度带电粒子探测器阵列（北大的 LACPU），能同时测 (d,p)、(d,³He)、(d,t) 多个道。

顺带一提，同样的方法在稳定核上也还在产出细活。佛罗里达州立大学 B. Kelly 等 10 位作者 2026 年 9 月在 arXiv 上挂了一篇，用超大型分裂磁谱仪测 ⁵⁰Ti 的 (d,p) 与 (d,pγ) 反应，一共认出 82 个激发态，并对 2p₃/₂、2p₁/₂、1f₅/₂、1g₉/₂、2d₅/₂ 五个中子轨道算了空位求和规则[40]（已被《物理评论 C》接收）。

5.4 融合：最丰中子的氧同位素

2026 年 3 月，FRIB 完成了第一次用氧-21 束流做的融合测量。氧-21 半衰期 3.4 秒，是这条同位素链上被测到的最丰中子的成员。

束流是 40Ar 主束打 10 毫米厚碳靶产生的，经 ARIS 谱仪筛选、气体阻止室热化、EBIT 电荷增育，最后由 ReA3 直线加速器加速到每核子 3.2 兆电子伏，强度约每秒 800 到 1000 个离子。靶探测器是印第安纳大学的 MuSIC，一个横向场、带 Frisch 栅的多采样电离室，用甲烷气同时当靶和探测介质[36]。

这项工作的价值在于把氧同位素链的融合激发函数补齐到 ²¹O。这个团队此前已经在 GANIL 测过 ¹⁹O、²⁰O，在圣母大学测过 ¹⁷,¹⁸O，三个实验室的数据能拼成一条完整的链。

5.5 多核子转移：壳效应是把双刃剑

要合成很丰中子的重核，多核子转移是一条热门路线。但壳效应在这里的作用是双向的。

中科院近代物理研究所的戴等 8 位作者 2026 年 8 月在《物理快报 B》上发表了一篇，系统算了 ¹³⁶Xe 打 ²⁰⁸Pb、²⁰⁴Hg、²⁰⁸Hg 三个靶的反应率[37]。结论是：对靠近入射道的少核子转移产物，壳效应增强产额；但对需要转移大量核子的丰中子产物，壳效应强烈抑制生成。抑制在双幻核铅-208 靶上最厉害。

所以他们的建议是：别拿幻数核当碰撞对象，改用非幻数靶，并具体推荐了「²³⁸U 打 ²⁰⁴Hg」这个体系，论证它能显著提高合成铱-203、锇-202 等丰中子核素的概率。

5.6 结构侧的顺带收获

这两个月还有几条值得一提的结构结果。

美国 FRIB 用 ⁸²Se 束打碳靶，在 ARIS 上首次产生并鉴别了 ⁶³Sc、⁶⁵Ti、⁶⁶Ti、⁶⁸V 四个新核素，另有一个事件与 ⁶¹Ca 相符[35]。第一作者是 O. B. Tarasov，16 位作者。

HIAF 在 2026 年 7 月 21 日通过工艺验收、进入调试阶段后，首个物理成果是成功产生并鉴别了缺中子核素铪-153，一共观测到 10 个离子，测量表明它是束缚或弱束缚核素，与理论预测一致。8 月以短通讯形式发表在《科学通报》英文版，33 位作者[41]。

理论侧，埃及的 W. M. Seif 与 A. Adel 用 Skyrme 泛函做了一套系统计算，把 N = 50、82、126 以及尚未观测到的 184 四个中子壳的基态性质、α 衰变半衰期放在一起比较，找出了明显的「壳驱动」规律：中子皮厚度随同位旋不对称度增大而增大，但在更重的壳上幅度反而变小[32]。他们对 N = 184 同位素的合成区间给出了 Z = 112 到 120 的估计。

六、读数有多靠谱

这一节是全篇最该记住的部分。

6.1 关于中子皮，不同探针给出的数差多少

拿铅-208 来说。

PREX 合作组用宇称破缺电子散射，把 PREX-1 与 PREX-2 合并，给出 0.283 ± 0.071 飞米。

九州大学若狭智嗣、田上慎吾、八寻正道三位，用反应截面的 Love-Franey t 矩阵折叠模型，重新分析 534.1、549、806 兆电子伏下质子打铅-208 的散射数据，给出 0.324 ± 0.047(实验) ± 0.009(理论) 飞米[22]。这个值比 PREX 还厚，而且误差更小。顺带他们给钙-48 的值是 0.144 ± 0.075 飞米。

质子弹性散射本身是老路。京都大学（现东北大学）钱广等 11 位作者 2010 年在 295 兆电子伏下测了 ²⁰⁴Pb、²⁰⁶Pb、²⁰⁸Pb 的中子密度分布[21]。

2025 年在钙-48 上又出了一条新结果。西南大学杨高华、匡勇与理化学研究所杨祖星、李志攀（4 位作者）用相对论冲量近似加深度神经网络，从质子弹性散射数据里提取钙-48 的中子皮，得到 0.199(17) 飞米，发表在《物理快报 B》862, 139301[23]。

这里有一个必须注意的版本陷阱：这篇的 arXiv 早期版本给的值是 0.211(11) 飞米，正式发表版改成了 0.199(17)。引数值必须认正式版。

还有一条容易被忽略的路：镜像核的电荷半径差。北京师范大学安荣、孙帅、曹李刚、张丰收（4 位作者）2023 年在《核科学与技术》上，用镜像对核的电荷半径差来约束对称能[24]；同一个团队 2024 年又用钙-48 与铅-208 的中子皮做了一个新的对称能量化[25]。

6.2 「合成」这件事本身正在被当成一门学问

把上面那些数放在一起求平均，是个陷阱。

2026 年 2 月，Adrian Azizi、Carlos A. Bertulani、Carlos Davila（东德州农工大学与达姆施塔特工业大学，3 位作者）在 arXiv 上挂了一篇，提出用层次贝叶斯框架来合成不同探针给出的中子皮[20]。核心想法是：不要假设所有测量都是同一个真值的有偏估计，而是给每种方法单独设一个偏差参数和内在离散度。他们拿锡同位素练手，推断从 ¹⁰⁰Sn 到 ¹⁴⁰Sn 的中子皮趋势，发现稳定核附近不确定度最小，往两端走越来越大。最终的后验把对称能斜率 L 压得很紧，因为中子皮主要对亚饱和密度的对称能压力敏感[20]。这篇已被《物理评论 C》接收，正式出版在 114(1), 015805 (2026)。

简单平均为什么错？因为它隐含假设「所有测量都是同一件事的无偏估计，而且报出来的误差已经包含了所有来源」。这两条都不成立。

6.3 ρ 的 2/3 次方：一条基准线

东德州农工大学的李宝安 2026 年 1 月在《欧洲物理杂志特刊》上发了一篇 60 页的综述[17]。

他讨论的是西门子在 1970 年代观察到的一个规律：在饱和密度附近，对称能大致按 (ρ/ρ₀)^(2/3) 走。这不是新发现——这篇综述的价值在于划边界。李宝安的结论是三条：

第一，这个 2/3 标度是稳健的，可以当基准用到 2ρ₀ 为止，但超过这个密度就失效。

第二，在 2ρ₀ 以下，实验和理论的结论大体一致，但曲率 K_sym 和高阶参数的不确定度仍然很大。

第三，超过 2ρ₀ 之后，不确定度暴涨，原因包括自旋-同位旋相关的张量力、三体力，以及由此产生的核子短程关联——这几样都还没被约束住[17]。

所以「对称能在 2ρ₀ 以下基本清楚了」这句话，只说对了一半：形状大体有共识，斜率还在打架。

6.4 一条不能忘的底层限制

回过头看第一节那三个 L 的数（15.3、40、64.59），误差棒都很大。这不是因为谁算错了。

根本原因是：所有地面实验测的都是有限核、有限温度、有限同位旋不对称度的系统，而你想要的是无限核物质在零温下的那条曲线。 从前者到后者，中间隔着输运模型、统计模型、能量密度泛函——每一层都要引入假设。R_i 依赖模型算出的密度-时间关系；同位旋标度依赖「碎片形成时处于统计平衡」；IYR 依赖 Fisher 模型的推广形式。

这也是为什么最近几年贝叶斯方法、模型无关重建、多信使联合约束这些做法会一起涌上来。它们对付的是同一件事：把「模型依赖」这个黑箱的边界画清楚。顺带说一句，把这批输运模型横向摆开做系统比对，本身就是一件大工程——2024 年那篇 136 位作者的综述[27]干的就是这件事；而把地面核实验与中子星观测放进同一个贝叶斯框架的代表工作，是 Tsang 等五人 2024 年在《自然·天文学》上的那篇[42]。

七、装置侧只讲三句

这一节跟《一次碰撞，怎么变成一条物态方程约束》划清边界，只交代现在能提供什么，不重复那篇讲过的装置细节。

FRIB 从 2022 年 5 月开始运行，服务着一个约 1800 人的用户群体[38]。它的价值主要在于把丰中子核的产生能力整体抬高了一档——前面 5.6 节那四个新核素、5.4 节那次氧-21 融合测量，都是这个能力的直接产物。

HIAF 2025 年 10 月首次出束，2026 年 7 月 21 日通过工艺验收进入调试。调试期间氧束和铋束的强度分别超过此前国际纪录的 3.0 倍与 7.5 倍。首个物理成果是铪-153[41]。

RAON 在费米能区这条线上不是主力——这一点在旧文里已经展开，这里不重复。

八、我的判断

实验这一侧，最扎实的是反应截面和电荷拾取这类「单观测量」测量。 碳-12 到碳-19 截面涨 50 倍这种结果，几乎没有模型依赖的余地，数字摆在那里。它们不直接给对称能，但把探针本身的尺寸标定了，是下面所有推论的起点。

同位旋分布这一侧，最成熟的是 R_i，最危险也是 R_i。 它的构造确实聪明，污染被标尺系统约掉了；但它给的是一个积分的、被模型翻译过的结果。你得到的不是一个「在 ρ = 0.8ρ₀ 处的对称能」，而是「在模型算出的那段密度历史里的平均效果」。2025 年 Ciampi 等人把它做到「模型无关测量」，那是相对于早先同一模型内比较而言的进步，不是说这个量本身脱离模型了。

最新的那条线——机器学习——现在最需要的不是更准的网络，而是更诚实的误差棒。 2026 年那两篇（彭丹与马春旺的贝叶斯网络、杨祖星等人的取向分类网络）都做了一件对的事：主动写出方法的失效区间。彭丹那篇明说氙靶上不行；杨祖星那篇把反应取向这个混淆变量单独拎出来处理。这两个「限制」比它们的中心结果更有价值。

理论上，最值得盯的是那本 59 页的综述[16]。 它做的事情是把这个领域所有能想到的混淆因素列成一张清单——形变、表面弥散、壳结构、团簇、模型依赖。任何一篇说「我测到了对称能」的文章，都该拿这张清单过一遍，看哪些因素没有排除。这份清单本身就是这个领域当前最诚实的自我评估。

如果只让我挑一个数字盯着看，我挑 L。 S₀ 的共识已经够好了（三十上下，误差几个兆电子伏），K_sym 现在谁都说不准（所以先别急）。L 是唯一一个「误差棒大到荒谬、但所有方法都在往中间挤」的量：15.3、40、64.59——三个数摆在那里，谁也没法宣布胜利。等哪天这三条线在同一个误差棒里重叠，那时候再说「我们知道对称能了」才不算吹牛。

九、最后一句话

碎片有多胖，是它在碰撞的那一瞬间被中子过剩的梯度推着走的结果；把这个推的过程反过来读，就得到核物质有多硬。中间隔着五层翻译，每一层都会走音。

参考文献

正文中的引用标注 [N] 与本列表第 N 条一一对应。

1. M. B. Tsang, T. X. Liu, L. Shi, P. Danielewicz, C. K. Gelbke, X. D. Liu, W. G. Lynch, W. P. Tan, G. Verde, A. Wagner, H. S. Xu 等（共 21 位作者），"Isospin Diffusion and the Nuclear Symmetry Energy in Heavy Ion Reactions"，*Physical Review Letters* 92(6), 062701 (2004)，在线 2004-02-10。DOI: 10.1103/PhysRevLett.92.062701。第一作者单位为密歇根州立大学国家超导回旋加速器实验室（NSCL），合作方含阿贡国家实验室、得克萨斯 A&M 大学商业分校等。反应体系 ¹²⁴Sn + ¹¹²Sn，每核子 50 兆电子伏，碰撞参数 6 飞米。

1. F. Rami, Y. Leifels, B. de Schauenburg, A. Gobbi, B. Hong 等（FOPI 合作组，共 70 位作者），"Isospin Tracing: A Probe of Nonequilibrium in Central Heavy-Ion Collisions"，*Physical Review Letters* 84(6), 1120–1123 (2000)，在线 2000-02-07。DOI: 10.1103/PhysRevLett.84.1120。实验为 ⁹⁶Ru + ⁹⁶Zr 在每核子 400 兆电子伏，在德国重离子研究中心（GSI）完成；这是同位旋示踪思想的原始论文，与后来 Tsang 等人定义的 R_i 公式是两段不同的谱系。

1. M. B. Tsang, W. A. Friedman, C. K. Gelbke, W. G. Lynch, G. Verde, H. S. Xu（共 6 位作者），"Isotopic Scaling in Nuclear Reactions"，*Physical Review Letters* 86(22), 5023–5026 (2001)，在线 2001-05-28。DOI: 10.1103/PhysRevLett.86.5023。同位旋标度（isoscaling）规律的原始论文。

1. H. S. Xu, M. B. Tsang, T. X. Liu, X. D. Liu, W. G. Lynch 等（共 20 位作者），"Isospin Fractionation in Nuclear Multifragmentation"，*Physical Review Letters* 85(4), 716–719 (2000)，在线 2000-07-24。DOI: 10.1103/PhysRevLett.85.716。

1. Payel Karmakar, Samir Kundu, T. K. Rana, S. Mallik, S. Manna, Ratnesh Pandey, K. Banerjee, Pratap Roy, S. Mukhopadhyay, Arijit Sen, T. K. Ghosh, P. Pant, G. Mukherjee, A. Sultana, D. Mondal, D. Pandit, R. Shil, K. Rani（共 18 位作者），"Symmetry energy and isoscaling behavior of fragments emitted in ¹⁴N,²⁰Ne + ¹¹²,¹¹⁶,¹²⁴Sn at 18–30 MeV/nucleon"，*Physical Review C* 112(2), 024614 (2025)，在线 2025-08-22。DOI: 10.1103/xtx3-mh6c（APS 2025 年起改用随机式 DOI）。实验在印度可变能量回旋加速器中心（VECC）K500 超导回旋加速器上完成。给出 C_sym 约 24–17 兆电子伏（对应激发能每核子 2.1–2.8 兆电子伏）。

1. Chun-Wang Ma, Fang Wang, Yu-Gang Ma, Chan Jin（共 4 位作者），"Isobaric yield ratios in heavy-ion reactions, and symmetry energy of neutron-rich nuclei at intermediate energies"，*Physical Review C* 83(6), 064620 (2011)。DOI: 10.1103/PhysRevC.83.064620。arXiv:1107.0131。署名单位含河南师范大学物理系与中科院上海应用物理研究所。用修正 Fisher 模型分析 ⁴⁸Ca 与 ⁶⁴Ni 弹核碎裂碎片。

1. Chun-Wang Ma, Heng-Li Song, Jie Pu, Tong-Lin Zhang, Sha Zhang, Shan-Shan Wang, Xin-Li Zhao, Li Chen（共 8 位作者），"Symmetry energy from neutron-rich fragments in heavy-ion collisions, and its dependence on incident energy, and impact parameters"，*Chinese Physics C* 37(2), 024102 (2013)。DOI: 10.1088/1674-1137/37/2/024102。用统计磨蚀-烧蚀模型计算 ⁶⁰Ni + ¹²C 在每核子 80 与 140 兆电子伏下的碎片，提取 a_sym/T、b_s/T、b_v/T。

1. Dan Peng, Chun Wang Ma（共 2 位作者），"Probing nuclear symmetry energy via isoscaling analysis of spallation fragments with a Bayesian neural network"，*Physical Review C* 113(2), 024605 (2026)，在线 2026-02-05。DOI: 10.1103/37mt-m922。

1. Zu-Xing Yang, Xiao-Hua Fan, Peng-Hui Chen, Shunji Nishimura（共 4 位作者），"Determining the orientation dependence of observables in heavy-ion collisions via deep neural networks"，*Physical Review C* 113(6), 064604 (2026)，在线 2026-06-10。DOI: 10.1103/6dtk-1tw2。署名单位：西南大学物理科学与技术学院、扬州大学物理科学与技术学院、理化学研究所仁科加速器科学中心。

1. Yongjia Wang, Fupeng Li, Qingfeng Li, Hongliang Lü, Kai Zhou，"Finding signatures of the nuclear symmetry energy in heavy-ion collisions with deep learning"，arXiv:2107.11012 [nucl-th]。署名单位含湖州师范学院、浙江工业大学、中科院近代物理研究所、华为海思、法兰克福高等研究院。

1. C. Ciampi, S. Piantelli, G. Casini, G. Pasquali 等（INDRA-FAZIA 合作组，共 47 位作者），"First results from the INDRA-FAZIA apparatus on isospin diffusion in ⁵⁸,⁶⁴Ni + ⁵⁸,⁶⁴Ni reactions at 32 and 52 MeV/nucleon"，*Physical Review C* 106(2), 024603 (2022)，在线 2022-08-08。DOI: 10.1103/PhysRevC.106.024603。

1. C. Ciampi, S. Piantelli, G. Casini, A. Ono 等（共 42 位作者），"Quasiprojectile breakup and isospin equilibration at Fermi energies: Potential indication of longer projectile-target contact times"，*Physical Review C* 108(5), 054611 (2023)，在线 2023-11-28。DOI: 10.1103/PhysRevC.108.054611。注意与第 11 条是两篇不同论文，作者数与侧重点均不同。

1. C. Ciampi, J. D. Frankland, D. Gruyer, N. Le Neindre 等（共 42 位作者），"Model-independent measurement of isospin diffusion in Ni-Ni systems at intermediate energy"，*Physical Review C* 111(4), 044601 (2025)，在线 2025-04-01。DOI: 10.1103/PhysRevC.111.044601。

1. C. Ciampi, S. Mallik, F. Gulminelli, D. Gruyer 等（共 34 位作者），"Constraining the nuclear symmetry energy with Fermi-energy heavy ion collisions"，*Physics Letters B* 868, 139815 (2025 年 9 月号)。DOI: 10.1016/j.physletb.2025.139815。

1. S. Mallik, F. Gulminelli, C. Ciampi, D. Gruyer（共 4 位作者），"Probing the Density Dependence of Nuclear Symmetry Energy Through Isospin Transport in Heavy-Ion Reactions"，*Universe* 12(7), 202 (2026)，在线 2026-07-06。DOI: 10.3390/universe12070202。arXiv:2605.06138。署名单位：加尔各答可变能量回旋加速器中心（VECC）、孟买 Homi Bhabha 国家研究院、卡昂诺曼底大学 LPC Caen 与法国国家科学研究中心、法兰西大学研究院、GANIL。

1. Bao-Jun Cai, De-Qing Fang, Yu-Gang Ma（共 3 位作者），"Neutron Skin Effects on Particle Emission in Heavy-Ion Collisions: A Topic Review with Astrophysical and Nuclear Structure Connections"，arXiv:2609.15856 [nucl-th] (2026-09-14)。59 页、60 图、469 条引用。署名单位：复旦大学现代物理研究所核物理与离子束应用教育部重点实验室、上海理论核物理研究中心（国家自然科学基金委与复旦大学共建）；马余刚另兼华东师范大学物理学院。这是一篇专题综述，不是新的实验测量。

1. Bao-An Li（1 位作者），"Beyond ρ^(2/3) scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy"，*The European Physical Journal Special Topics* (2026)，在线 2026-01-20。DOI: 10.1140/epjs/s11734-026-02124-1。arXiv:2510.05508（v1 2025-10-07，v2 2026-01-06）。作者单位为东德州农工大学。

1. William G. Newton, Gabriel Crocombe（共 2 位作者），"Nuclear symmetry energy from neutron skins and pure neutron matter in a Bayesian framework"，*Physical Review C* 103(6), 064323 (2021)。DOI: 10.1103/PhysRevC.103.064323。arXiv:2008.00042。署名单位为得克萨斯 A&M 大学商业分校物理与天文系。用非信息先验得 J = 31.3 +4.2/−5.9 MeV、L = 40 +34/−26 MeV；用纯中子物质先验得 J = 31.9 ± 1.3 MeV、L = 37 +9/−8 MeV。

1. Zhen Zhang, Lie-Wen Chen（共 2 位作者），"Bayesian inference of the symmetry energy and the neutron skin in ⁴⁸Ca and ²⁰⁸Pb from CREX and PREX-2"，*Physical Review C* 108(2), 024317 (2023)，在线 2023-08-21。DOI: 10.1103/PhysRevC.108.024317。arXiv:2207.03328。署名单位：中山大学中法核工程与技术学院（张振）、上海交通大学物理与天文学院与粒子物理与宇宙学重点实验室（陈立文）。合并 CREX 与 PREX-2 得 S(ρ₀) = 30.2 +4.1/−3.0 MeV、L = 15.3 +46.8/−41.5 MeV（90% 置信度）。

1. Adrian Azizi, Carlos A. Bertulani, Carlos Davila（共 3 位作者），"Hierarchical Bayesian analysis of neutron-skin thicknesses and implications for the symmetry-energy slope"，*Physical Review C* 114(1), 015805 (2026)。DOI: 10.1103/w977-dnxq。arXiv:2602.04794（v2 2026-02-04）。署名单位：东德州农工大学物理系；Bertulani 另兼达姆施塔特工业大学核物理研究所。

1. J. Zenihiro, H. Sakaguchi, T. Murakami, M. Yosoi, Y. Yasuda, S. Terashima, Y. Iwao, H. Takeda, M. Itoh, H. P. Yoshida, M. Uchida（共 11 位作者），"Neutron density distributions of ²⁰⁴,²⁰⁶,²⁰⁸Pb deduced via proton elastic scattering at E_p = 295 MeV"，*Physical Review C* 82(4), 044611 (2010)。DOI: 10.1103/PhysRevC.82.044611。

1. Tomotsugu Wakasa, Shingo Tagami, Masanobu Yahiro（共 3 位作者），"Skin values of ²⁰⁸Pb and ⁴⁸Ca determined from reaction cross sections"，*Results in Physics* 43, 106101 (2022 年 12 月号)。DOI: 10.1016/j.rinp.2022.106101。arXiv:2211.05312（v2 于 2025-06-28 更新以补期刊信息）。署名单位为九州大学物理系。给出铅-208 中子皮 0.324 ± 0.047(实验) ± 0.009(理论) 飞米、钙-48 中子皮 0.144 ± 0.075 飞米。

1. G. H. Yang, Y. Kuang, Z. X. Yang, Z. P. Li（共 4 位作者），"Extracting neutron skin from elastic proton-nucleus scattering with deep neural network"，*Physics Letters B* 862, 139301 (2025 年 3 月号)。DOI: 10.1016/j.physletb.2025.139301。arXiv:2311.11676。署名单位：西南大学物理科学与技术学院；杨祖星另兼理化学研究所仁科加速器科学中心。正式版给出钙-48 中子皮 0.199(17) 飞米，与 arXiv 早期版的 0.211(11) 飞米不同。

1. Rong An, Shuai Sun, Li-Gang Cao, Feng-Shou Zhang（共 4 位作者），"Constraining nuclear symmetry energy with the charge radii of mirror-pair nuclei"，*Nuclear Science and Techniques* 34(8), 119 (2023 年 8 月)。DOI: 10.1007/s41365-023-01269-1。arXiv:2303.14667。

1. Rong An 等（共 4 位作者），"New quantification of symmetry energy from neutron skin thicknesses of ⁴⁸Ca and ²⁰⁸Pb"，*Nuclear Science and Techniques* 35(10), 182 (2024)。DOI: 10.1007/s41365-024-01551-w。

1. M. B. Tsang, J. R. Stone, F. Camera, P. Danielewicz, S. Gandolfi 等（共 19 位作者），"Constraints on the symmetry energy and neutron skins from experiments and theory"，*Physical Review C* 86(1), 015803 (2012)。DOI: 10.1103/PhysRevC.86.015803。

1. A. Sorensen, K. Agarwal, K. W. Brown, Z. Chajęcki, P. Danielewicz 等（共 136 位作者），"Dense nuclear matter equation of state from heavy-ion collisions"，*Progress in Particle and Nuclear Physics* 134, 104080 (2024 年 1 月号)。DOI: 10.1016/j.ppnp.2023.104080。

1. 孙向东、刘阳阳、李昂、肖志刚、张宇（共 5 位作者），《致密物质状态方程的天文与地面研究》，*中国科学：物理学 力学 天文学* 56(7), 272010 (2026)，在线 2026-02-03。DOI: 10.1360/SSPMA-2025-0201。注意刊名为中文刊的英文题名 *Scientia Sinica Physica, Mechanica & Astronomica*，与 Science China Physics, Mechanics & Astronomy 是不同的刊。给出 K₀ = 250.56 +16.51/−18.79 MeV、m*/m = 0.60 ± 0.01、S₀ = 28.46 +4.36/−2.24 MeV、L₀ = 64.59 +13.25/−9.38 MeV。

1. Y. L. Ye, X. F. Yang, H. Sakurai, B. S. Hu（共 4 位作者），"Physics of exotic nuclei"，*Nature Reviews Physics* 7(1), 21–37 (2025 年 1 月号)；DOI 编号为 10.1038/s42254-024-00782-5（线上 2024 年），另有更正声明同卷 7(1), 62。署名单位含北京大学、东京大学、理化学研究所仁科中心。

1. 林等（共 2 位作者），"Halos and Multineutron Correlations in Light Neutron-Rich Nuclei"，*Particles* 9(1), 27 (2026)，在线 2026-03-16。DOI: 10.3390/particles9010027。署名单位：北京大学核物理与技术国家重点实验室、理化学研究所仁科中心。属「晕核发现四十周年国际研讨会」专辑。

1. Hongna Liu（刘红娜，1 位作者），"Exploring the migration of magic numbers in neutron-rich systems via knockout reactions"，*Frontiers in Physics* 14, 1778872 (2026)，在线 2026-05-12。DOI: 10.3389/fphy.2026.1778872。署名单位：北京师范大学物理与天文学院、射线束技术教育部重点实验室。

1. W. M. Seif, A. Adel（共 2 位作者），"Shell-driven similarities in structure and decay characteristics across N = 50, 82, 126, and 184 isotones"，*Physical Review C* 113(3), 034304 (2026)，在线 2026-03-04。DOI: 10.1103/yft9-22r3。署名单位为开罗大学物理系。

1. 北京大学团队（综述，2026），《单核子转移反应研究轻丰中子核壳演化的实验进展》，*科学通报*（中文版）2026 年。综述内容含 ²H(¹¹Be,p)¹²Be 的 d 波侵入组态约 57%、²H(¹⁶C,p)¹⁷C 的 s₁/₂–d₅/₂ 轨道简并证据，以及北京大学大接收度带电粒子探测器阵列（LACPU）。

1. J.-C. Zhang, B.-H. Sun, I. Tanihata, S. Terashima 等（共 37 位作者），"Charge Pickup Reaction Cross Section for Neutron-Rich p-Shell Isotopes at 900A MeV"，*Physical Review X* 15(3), 031004 (2025)，在线 2025-07-02。DOI: 10.1103/k9tj-jq8s。第一作者张寂潮（北京航空航天大学博士毕业生，现为大阪大学博士后），共同通讯作者孙保华与谷畑勇夫；合作单位含德国 GSI、加拿大 TRIUMF、日本大阪大学等十所机构。报道碳-12 到碳-19 的电荷拾取截面增长约 50 倍，而基于稳定核的经验预期只有 2.5 倍。

1. O. B. Tarasov 等（共 16 位作者），"Discovery of new isotopes in the fragmentation of ⁸²Se and insights into their production"，*Physical Review C* 112(3), 034604 (2025)。DOI: 10.1103/573p-7fjp（APS 随机式 DOI，按卷期页拼装 DOI 会返回 404）。实验在美国稀有同位素束流装置（FRIB）的 ARIS 分离器上完成，⁸²Se 束能量每核子 228 兆电子伏，首次产生并鉴别 ⁶³Sc、⁶⁵Ti、⁶⁶Ti、⁶⁸V。

1. A. Leon 等（共 21 位作者），"Neutron dynamics in the fusion of neutron-rich O ions with ¹²C"，arXiv:2609.03914 [nucl-ex] (2026-09-03)。实验于 2026 年 3 月在美国 FRIB 完成，使用 MuSIC 活性靶探测器与 ReA3 再加速束流，²¹O 强度约每秒 800–1000 个离子、每核子 3.2 兆电子伏。

1. 戴等（共 8 位作者），"Role of shell effects on producing neutron-rich nuclei in the N = 126 region via multinucleon transfer reactions"，*Physics Letters B* 879, 140614 (2026 年 8 月号)。DOI: 10.1016/j.physletb.2026.140614。署名单位为中科院近代物理研究所等。

1. A. Gade, A. O. Macchiavelli（共 2 位作者），"Toward the Limits of Nuclear Existence: First Science with FRIB"，*Annual Review of Nuclear and Particle Science* 76(1), 245–268 (2026)。DOI: 10.1146/annurev-nucl-121423-100941。

1. Yuki Nakanishi, Junki Tanaka, Atsushi Tamii, Shimpei Endo（共 4 位作者），"Probing Near-Threshold s-Wave Components in Heavy Nuclei via Coulomb-Assisted Neutron Transfer"，arXiv:2604.27605 [nucl-th] (2026-04-30)。

1. B. Kelly, M. Spieker, L. T. Baby, S. Baker, A. L. Conley, D. Houlihan, K. W. Kemper, E. Litvinova, N. Tsoneva, A. Volya（共 10 位作者），"Detailed experimental study of excited states in ⁵⁰Ti via the (d,p) and (d,pγ) reactions"，arXiv:2609.09412 [nucl-ex] (2026-09-08)，已被《物理评论 C》接收。署名单位：佛罗里达州立大学物理系；利特维诺娃另属西密歇根大学；措涅娃属罗马尼亚 ELI-NP 与 IFIN-HH。

1. 中科院近代物理研究所等（共 33 位作者，第一作者李），"Discovery of isotope ¹⁵³Hf near the proton drip line via isochronous mass spectrometry"，*Science Bulletin* 71(16), 4099–4102 (2026)。DOI: 10.1016/j.scib.2026.06.056。依托强流重离子加速器装置（HIAF）的谱仪环 SRing，用等时性质谱技术完成；HIAF 于 2026 年 7 月 21 日通过工艺验收，这是其调试阶段的首个物理成果。

1. C. Y. Tsang, M. B. Tsang, W. G. Lynch, R. Kumar, C. J. Horowitz（共 5 位作者），"Determination of the equation of state from nuclear experiments and neutron star observations"，*Nature Astronomy* 8(3), 328–336 (2024)，在线 2024-01-05。DOI: 10.1038/s41550-023-02161-z。arXiv:2310.11588。署名单位含密歇根州立大学、印第安纳大学。

1. Lie-Wen Chen, Che Ming Ko, Bao-An Li（共 3 位作者），"Light cluster production in intermediate energy heavy-ion collisions induced by neutron-rich nuclei"，*Nuclear Physics A* 729, 809–834 (2003)。DOI: 10.1016/j.nuclphysa.2003.09.010。arXiv:nucl-th/0306032。软硬对称能给出的 t/³He 比相差约 17%、而改不可压缩系数与核子核子截面影响不到 4% 这组数字出自此文。

1. Gao-Chan Yong, Bao-An Li, Lie-Wen Chen, Xun-Chao Zhang（共 4 位作者），"Triton-³He relative and differential flows as probes of the nuclear symmetry energy at supra-saturation densities"，*Physical Review C* 80(4), 044608 (2009)，在线 2009-10-15。DOI: 10.1103/PhysRevC.80.044608。arXiv:0906.0939。反应体系 ¹³²Sn + ¹²⁴Sn，每核子 400 兆电子伏。

核实与纠错记录

下表逐条列出本轮核查中发现的、容易被二手材料或本知识库旧文记述错误的地方。

序号	常见说法	实际情况	依据
1	若狭等人的中子皮论文是 2025 年发表	正式发表于 *Results in Physics* 43, 106101，2022 年 12 月号；arXiv:2211.05312 的 v2 是 2025 年才更新以补上期刊信息，容易误读成年份	CrossRef 直查 DOI 10.1016/j.rinp.2022.106101，印刷期 2022-12
2	FRIB 首次观测 ⁶³Sc、⁶⁵,⁶⁶Ti、⁶⁸V 那篇的第一作者是 Shimizu	第一作者是 O. B. Tarasov，共 16 位作者	CrossRef 直查 DOI 10.1103/573p-7fjp
3	该文的 DOI 可按「期刊缩写加卷期文章号」拼出来	拼出来的 10.1103/PhysRevC.112.034604 返回 404；真 DOI 是随机式的 10.1103/573p-7fjp	CrossRef 直查，构造 DOI 报 404
4	用神经网络从质子弹性散射提取的钙-48 中子皮是 0.211(11) 飞米	那是 arXiv 早期版本；正式发表版是 0.199(17) 飞米	PLB 862, 139301 (2025) 正式版；arXiv:2311.11676 早期稿
5	INDRA-FAZIA 的同位旋扩散结果只有一篇	至少有四篇定位不同：PRC 106, 024603 (2022) 首次结果（47 作者）、PRC 108, 054611 (2023) 准弹破裂通道（42 作者）、PRC 111, 044601 (2025) 模型无关镍镍测量（42 作者）、PLB 868, 139815 (2025) 对称能约束（34 作者）；作者数各不相同，混引会对调团队	CrossRef 逐条直查
6	R_i 这个量是 Tsang 2004 提出的	公式的思想先声来自 Rami 等人 2000 年（FOPI 合作组，70 作者，⁹⁶Ru+⁹⁶Zr 每核子 400 兆电子伏）；Tsang 2004 把它用在锡锡体系并接上对称能约束，是两段谱系	PRL 84, 1120 (2000)；PRL 92, 062701 (2004)
7	t/³He 比对核物质性质「普遍敏感」	它对对称能敏感（软硬差约 17%），但对不可压缩系数（380 改 201 兆电子伏）与介质中核子核子截面只差不到 4%——这恰恰是它作为干净同位旋探针的优点	轻团簇产生系列计算（arXiv:nucl-th/0306032）；arXiv:1310.3967
8	李宝安 2026 年那篇综述提出了 ρ 的 2/3 次方标度	该标度由 Siemens 在 1970 年代提出；这篇综述做的是划边界——确认它到 2ρ₀ 有效、更高密度失效，并指出曲率 K_sym 与高阶参数仍无约束	EPJ Special Topics (2026)，DOI 10.1140/epjs/s11734-026-02124-1
9	arXiv:2609.15856 是一篇新的实验测量	是复旦大学蔡宝军、方德清、马余刚三人的专题综述，59 页、60 图、469 条引用，用来梳理中子皮影响粒子发射的全部混淆因素	arXiv:2609.15856 v1 (2026-09-14)，作者 3 人
10	不同探针给出的铅-208 中子皮应当彼此接近	差距是实的：PREX 合并值 0.283 ± 0.071 飞米；反应截面路线给 0.324 ± 0.047(实验) ± 0.009(理论) 飞米；同一路线给钙-48 只有 0.144 ± 0.075 飞米。铅与钙给出的「软硬」方向并不同步，所以才有「单独看 CREX 与 PREX-2 在 68.3% 置信度上强冲突、到 90% 才相容」这回事	数值分别见本表依据的第 19、22 条文献
11	《Physics of exotic nuclei》是 2024 年发表	DOI 编号含 2024，但卷期是 *Nature Reviews Physics* 7(1), 21–37（2025 年 1 月号）；另有一篇更正声明在同卷 7(1), 62	CrossRef 直查 DOI 10.1038/s42254-024-00782-5
12	旧记忆把 Ciampi 的两篇都记成 34 位作者	PLB 868, 139815 是 34 作者，PRC 111, 044601 是 42 作者——「均 34 作者」是错的	CrossRef 逐条直查
13	《中国科学》那篇刊名可写作 Science China Physics, Mechanics & Astronomy	该刊英文题名是 *Scientia Sinica Physica, Mechanica & Astronomica*，与 Science China 系列是两本不同的刊	CrossRef 直查 DOI 10.1360/SSPMA-2025-0201
14	「碎片越丰中子，越能给出强的对称能约束」	方向要分开：等压产额比方法在更丰中子的碎片上对参考 IYR 的选取更敏感，作者自己指出这是方法软肋；贝叶斯同位旋标度方法在重靶（氙）上表现明显变差	PRC 37, 024102 (2013)；PRC 113, 024605 (2026)