跳到正文
返回核反应
至真探索 · 核反应

一次碰撞，怎么变成一条物态方程约束

2026-09-22

一、先把要问的问题说清楚

核物理里有一个听起来很朴素、做起来极难的问题：把一堆质子和中子挤在一起，它们的能量会怎么变？

这个问题有一个专门的名字，叫核物质物态方程。它不是某一种原子核的性质，而是一种"理想化材料"的性质——假设质子中子无限多、均匀铺开、不带电磁效应，那么每增加一点密度，能量会上升多少。中子星内部就是这种材料，超新星爆发的核心也是这种材料。

要描述它，最常用的是两个数。

第一个数刻画"对称核物质"，也就是质子数恰好等于中子数的那种。这个数叫不可压缩系数，符号 K₀，单位 MeV。它的意思是：你把核物质压缩一小下，它会有多"硬"、多想弹回来。通俗一点说，它就是核物质的弹簧常数。K₀ 越大，同样的压缩需要越多能量。

第二个数刻画"中子比质子多"带来的额外代价。这叫对称能，符号 S(ρ) 或者 E_sym(ρ)。核物质有一个偏好：质子中子成对的时候最舒服。一旦中子多出很多，能量就被顶上去。这个多出来的能量就是对称能。它的关键点是它随密度变化——在接近饱和密度（记作 ρ₀，大约 0.16 个核子每立方费米）附近是多少，在低于饱和密度时怎么走，在高于饱和密度时怎么走，这三段的答案可以完全不同。

一句常被引用的话：核物质在密度 ρ、同位旋不对称度 δ 下的每核子能量，可以写成

E(ρ, δ) = E(ρ, 0) + S(ρ)·δ² + 更高阶项

其中 δ = (ρ_n − ρ_p)/(ρ_n + ρ_p)。这一条在中等不对称度下相当好用，背后是一个被称为"抛物线近似"的假设。

现实中，K₀ 和 S(ρ) 都不是从第一性原理直接算出来就能信的。它们必须被测量。而测量它们的唯一办法，是把原子核撞碎。

---

二、为什么偏偏是每核子 20 到 100 MeV

重离子碰撞的能量跨度极大，从每核子几 MeV 到几千 GeV。不同的能量窗口看的是完全不同的物理。

高能端（每核子几个 GeV 到 TeV）撞出来的是夸克胶子等离子体，那是在探 QCD 相图，跟本文要讲的 EOS 参数是两条路。

低能端（每核子几 MeV）撞出来的是熔合反应，两个核粘成一个，主要看位垒。

中间这一段，也就是每核子 20 到 100 MeV，被单独称为"费米能区"。 这个名字来自一个事实：入射核子的相对速度已经和核内核子的费米速度可比。这个窗口很特殊，因为它同时具备三件事：

一是能把核物质压到 1 到 2 倍饱和密度，也就是既能探亚饱和也能碰到超饱和；二是时间尺度合适，反应持续几十到几百 fm/c，长到足以发生大量核子交换，短到不至于完全热化；三是碎片仍然带同位素信息，撞完之后飞出来的碎片，其 N/Z 比值还记着碰撞过程中"谁给了谁多少中子"。

最后这一点是整个方法的支点。因为这个窗口里，你测的不只是有多少粒子飞出去，而是飞出去的粒子里中子和质子的比例。而中子和质子的分布，恰好被对称能控制。

---

三、第一道难关：你根本不知道它撞得有多正

任何一次重离子碰撞，都有一个最基础的几何量：碰撞参数 b，也就是两个核擦肩而过时，中心之间的横向距离。b = 0 是正面对撞，b 接近两核半径之和是擦边。

问题在于，b 是不可直接测量的。探测器看到的是终态：飞出来的粒子有多少个、往哪些角度飞、各自什么电荷质量。b 只能从这些终态反推。

过去几十年最通行的做法叫"多重数切割"：统计一下带电粒子的数目，把数目最多的前百分之几叫"中心碰撞"，把数目少的叫"周边碰撞"。这个做法的问题很隐蔽：多重数和 b 的关系不是一一对应，而是一个带涨落的统计关系。同一个 b，多重数会在一个范围内波动；反过来，同一个多重数区间，对应的 b 也有一个分布。

忽略这个涨落会怎样？会系统性低估 b。

这件事在 2021 年被 J. D. Frankland、D. Gruyer 等人正式解决了。论文是《Model independent reconstruction of impact parameter distributions for intermediate energy heavy ion collisions》，发表在 *Physical Review C* 104, 034609 (2021)，作者共 20 人，署名单位以 GANIL 和法国卡昂的 LPC Caen 为主，署名合作组是 INDRA Collaboration。预印本编号 arXiv:2011.04496。

方法本身是概率反演，不是机器学习，也不是拟合。它的逻辑是这样：先假设一个条件概率 P(X | c_b)，意思是"给定碰撞参数（用归一化的 c_b 表示），某个观测量 X 会取什么值"，这个条件概率的宽度用一个 gamma 分布来描述涨落。然后用贝叶斯定理把它反过来：从测到的 X 分布，反推真实的 c_b 分布。方法本身借用的是超高能重离子碰撞里已经提出来的一套思路。

结果有两个，都值得记住。第一，用新方法重新看那些"中心碰撞"事件，推出的平均碰撞参数比老办法算出来的明显更大，而且能量越低差得越多。第二，这个方法还能给出一个以前拿不到的量：外推到 b = 0 时观测量应该在什么值。

这一篇是后面所有事情的地基。如果连"撞得多正"都算错了，那拿碰撞参数依赖去和输运模型对照就无从谈起。

---

四、第二道难关：拿什么当"温度计"

有了 b，还缺一个能读出对称能的量。

这个量叫同位旋输运比，记作 R_i。定义如下：

R(X_i) = (2·X_i − X_AA − X_BB) / (X_AA − X_BB)

其中 X 取准弹核残骸的中子质子比 N/Z，下标 A 和 B 代表两个同位素不同的反应伙伴，下标 i 表示混合体系。也就是说，你用两个不对称体系的测量值（AB 和 BA），减去两个对称参考体系（AA 和 BB）的值，再归一化。

这个构造有几层讲究。它把整体归一化因子消掉了，所以不依赖探测效率的绝对标定；它把最终结果的取值范围锁在 −1 到 +1 之间：0 表示完全均衡（两边中子质子比变得一样），正负 1 表示完全没有发生交换。最关键的是，这个比值被设计成对末态相互作用不敏感——碰撞结束之后碎片还会蒸发、还会衰变，那些过程会污染很多观测量，但对这个经过归一化的比值影响要小得多。

谱系要说准，因为这里有两个必须分清的阶段。概念上的先声是 F. Rami 等人 2000 年的工作：《Isospin tracing: A probe of nonequilibrium in central heavy-ion collisions》，*Physical Review Letters* 84(6), 1120–1123 (2000)，署名合作组是 FOPI，共 70 位作者，主要单位是斯特拉斯堡的 IReS 和克莱蒙费朗的 LPC。他们用 ⁹⁶Ru 与 ⁹⁶Zr 四种弹靶组合、在 400 A·MeV 下，用"测到的质子数"和"t/³He 产额比"两个同位旋示踪量去映射弹靶核子之间的混合程度，结论是即使在最中心的碰撞里全局均衡也没有达到。注意这个能量在费米能区之上，而且它给的是"同位旋示踪"的思想，还不是后来那个归一化的输运比公式。

把归一化输运比系统化、并明确用来探对称能密度依赖的，是 M. B. Tsang 等人 2004 年的文章《Isospin Diffusion and the Nuclear Symmetry Energy in Heavy Ion Reactions》，*Physical Review Letters* 92, 062701 (2004)，2004 年 2 月 13 日发表。他们用的是 ¹¹²Sn+¹¹²Sn 与 ¹²⁴Sn+¹²⁴Sn 作对称参考、¹¹²Sn+¹²⁴Sn 与 ¹²⁴Sn+¹¹²Sn 作不对称体系，入射能 E/A = 50 MeV，并由此判断准弹核与准靶核并未达到同位旋均衡，从而可以评估同位旋输运速率。今天 INDRA-FAZIA 用的 R_i 公式和"两个对称体系作参考"的构造，谱系在这里。

物理图像很直白：对称能如果随密度涨得很快（"硬"），那么在压缩最厉害的那一段，中子就会更强烈地被推出去，均衡得更快，R 就更靠近 0。 反过来，对称能如果涨得慢（"软"），均衡就慢，R 就离 0 远。所以测 R 随 b 的变化曲线，就等于在给对称能的密度依赖做体检。

---

五、INDRA-FAZIA：两个老家伙凑成一对

要同时满足"看得全"和"分得清"，靠一台探测器做不到。INDRA-FAZIA 的思路是把两台互补的设备拼起来。

INDRA 是 1990 年代设计的，由硅探测器、碘化铯闪烁体和电离室组成，共 564 个探测单元，排成 17 个同心环，覆盖 4π 立体角的约 90%。它对电荷的识别能到 Z 等于 92，但同位素分辨只到 A 小于等于 8——也就是说，它能数清碎片有几个质子，却分不清这个碎片是氦-3 还是氦-4 的重版本。它的长处是看得全。

FAZIA（Four Pi A and Z Identification Array）2015 年投入运行。它用的是叠层望远镜结构，一层一层堆了多种探测器，配合数字化电子学和脉冲形状分析，同位素分辨能力延伸到 Z 约 25。对一片从准弹核上剥下来的中等质量碎片，它能把质量数和中子数都定下来。它的长处是分得清。

耦合方式是：拆掉 INDRA 最前面 5 个环，换上 12 个 FAZIA 模块，让 FAZIA 占据最前角——因为最重的、最像准弹核的那些碎片都往小角度飞。

参数记清楚：FAZIA 采用 wall 配置，覆盖实验室系极角约 1.7° 到 13.7°，共 192 个望远镜；INDRA 覆盖互补的 14° 到 176°。合起来的几何接收度约 80%——注意这个数字比 INDRA 单独工作时的 90% 要低，也比 FAZIA 单独工作的约 70% 要高，是耦合后的折中。

装置装在法国 GANIL 的 D5 实验厅。

---

六、第一次实战：2019 年的 E789 实验

INDRA-FAZIA 耦合完成后，第一批数据是 2019 年的 E789 实验。它选了四个体系：

  • ⁵⁸Ni + ⁵⁸Ni（对称参考）
  • ⁶⁴Ni + ⁶⁴Ni（对称参考）
  • ⁵⁸Ni + ⁶⁴Ni（不对称）
  • ⁶⁴Ni + ⁵⁸Ni（不对称，前后对调）

并在两个能量下各测一遍：每核子 32 MeV 和 52 MeV。

为什么挑镍？因为 ⁵⁸Ni 的质子数和中子数都是 28，几乎完全对称；⁶⁴Ni 是中子富集的稀有同位素，多出 8 个中子。把同一对核前后对调打一遍，可以很大程度上抵消几何效应和探测器不均匀性——这就是 R_i 公式里 AB 和 BA 两项存在的理由。

这一批数据陆续产出了一串文章，路径很清晰：

  1. **C. Ciampi 等，*Physical Review C* 106, 024603 (2022)**。这是 INDRA-FAZIA 的第一批结果，同时给出 32 和 52 A·MeV 两个能量。
  2. **C. Ciampi 等，*Physical Review C* 108, 054611 (2023)**。把分析细化到准弹蒸发通道。
  3. **L. Baldesi 等，*Journal of Physics: Conference Series* 2586, 012039 (2023)。会议报告，聚焦半周边和周边碰撞。它给出一个符合直觉但很关键的结论：32 MeV/nucleon 下的同位旋均衡比 52 MeV/nucleon 更强**。原因是能量低时两个核待在一起的时间更长，交换更充分。
  4. **C. Ciampi 等，*Physical Review C* 111, 044601 (2025)**。这是本篇的核心方法论文。
  5. **C. Ciampi 等，*Physics Letters B* 868, 139815 (2025)**。这是本篇的核心结果论文。

---

七、2025 年的两篇：方法篇与结果篇

先把作者和单位说准，因为这里最容易出错。

《Model-independent measurement of isospin diffusion in Ni-Ni systems at intermediate energy》，*Physical Review C* 111, 044601 (2025)，2025 年 4 月发表，DOI 10.1103/PhysRevC.111.044601，预印本 arXiv:2412.13648。第一作者 C. Ciampi，通讯单位是 GANIL（法国国家重离子加速器实验室，CEA-DRF 与 CNRS-IN2P3 联合单位，在卡昂）。署名合作组为 INDRA-FAZIA Collaboration，共 34 位作者。合作者名单里能识别出 J. D. Frankland、D. Gruyer、N. Le Neindre、S. Mallik、R. Bougault、A. Chbihi、L. Baldesi、S. Barlini、E. Bonnet、B. Borderie、A. Camaiani、G. Casini、D. Dell'Aquila、J. A. Dueñas、Q. Fable、F. Gramegna、B. Hong、S. Kim 等人。值得一提的是名单里有韩国作者（B. Hong、S. Kim），说明这是真正的国际合作。

这一篇做的事情，是把前面两条线第一次拼在一起：用 INDRA 数据集做模型无关的碰撞参数重建，用 FAZIA 数据集测准弹残骸的 N/Z，然后把结果画成"同位旋输运比随碰撞参数变化"的曲线。结论是：随着碰撞越来越正，两个不对称反应给出的准弹核 N/Z 越来越接近——这就是同位旋均衡正在发生的直接证据。文章还专门说明，以前两篇文章（Ciampi 2022、Ciampi 2023）里用 AMD+GEMINI++ 模拟提取碰撞参数，与这次的模型无关办法在重叠的 b 区间内结果一致，从而验证了老办法在这个区间里的可用性。

《Constraining the nuclear symmetry energy with Fermi-energy heavy ion collisions》，*Physics Letters B* 868, 139815 (2025)，2025 年 8 月 10 日发表，DOI 10.1016/j.physletb.2025.139815，法国国家开放档案编号 hal-04851007。作者是 C. Ciampi（GANIL）、S. Mallik（加尔各答 VECC 与孟买 HBNI）、F. Gulminelli（LPC Caen）、D. Gruyer（LPC Caen）、J. D. Frankland（GANIL） 等共 34 位作者。注意 Mallik 的署名机构：孟买 Homi Bhabha National Institute 与加尔各答 Variable Energy Cyclotron Centre，这两处常被写漏。

这一篇做的事情是正面对照：第一次把 BUU 输运计算与 INDRA-FAZIA 数据推出的、带碰撞参数依赖的同位旋输运比放在一起比较。用的泛函来自两条不同路线——一条是从头算（ab initio）结果拟合得到的，一条是唯象方法给出的。

结果是一句话：在 ρ/ρ₀ = 1.01 处，对称能 S = 29.0 ± 0.7 MeV。

同时，文章做了一件很有分量的事：研究重子密度与同位旋流密度随时间的演化，去确定这次实验到底显著地探到了哪个密度区间。这正是它敢把结果标在一个具体密度上的原因——不是随便挑一个 ρ₀ 附近的值，而是从时间演化里算出敏感窗口。

这里必须点出一个容易被误传的细节：这次实验实际探测的密度是 ρ/ρ₀ = 1.01，也就是略高于饱和密度，只比 ρ₀ 低 1%。

它不是"在 0.5 倍饱和密度深处的严格亚饱和测量"。之所以常被归入"亚饱和方向"，是因为与 SπRIT 的约 2ρ₀ 相比，它整体落在饱和密度附近及偏下的一段，属于低密度侧的约束。把 S = 29.0 ± 0.7 MeV 说成"亚饱和深处的对称能"是过度解读，把 ρ/ρ₀ = 1.01 说成"0.5ρ₀"更是错的。

---

八、第二次实战：2022 年的 E818 实验

E789 之后，INDRA-FAZIA 做了第二个 campaign：E818，执行时间是 2022 年 5 月 7 日到 5 月 17 日。

体系换成两个：⁵⁸Ni + ⁵⁸Ni 和 ³⁶Ar + ⁵⁸Ni，入射能 74 MeV/nucleon。换体系的目的很明确：用两个质量不同的弹核，去制造两种大小不同的准弹核，从而研究准弹核的"气化"和"类气化"事件，并把介质内效应（in-medium effect）的测量扩展到更重的团簇上。

题外但很有意思的一点：E818 之前，INDRA 换了全新的电子学和电源系统，因为要替掉 30 年前的老电子学和电缆。这笔钱来自诺曼底大区的专项拨款，改造花了两年。

人员信息记准确：E818 的两位发言人是 Rémi Bougault 和 Giovanni Casini；数据负责人是 John D. Frankland 和 Diego Gruyer；GANIL 方面的科学协调人是 Nicolas Le Neindre。实验数据集的 DOI 是 10.26143/GANIL-2022-E818_20，数据体量 10.83 TB，托管在 IN2P3 计算中心的 CC-IN2P3。

E818 的科学动机来自更早的一批分析：用 INDRA 的重离子数据，从 ²H、³H、³He、⁴He、⁶He 的化学平衡常数里提取核物质与星体物质（含轻团簇）的物态方程，并用贝叶斯方法给团簇的内部配分函数加一个介质内修正。那次分析曾发现数据支持一个普适的标量 sigma 介子耦合修正，且这个修正比早先用德克萨斯 A&M 数据（未含介质内效应）得到的结果更大——意味着团簇结合能受影响更小、融化密度更大、超新星物质里团簇的贡献更多。

顺带记一笔：FAZIA 的数据集里，有明确的中国作者署名（Sun Baohua），以及韩国作者（Hong Byungsik、Kim Sunji、Kweon Min Jung、Nam Seonho），说明这个合作组的成员构成相当国际化。

---

九、另一条路：SπRIT 往相反方向看

如果说 INDRA-FAZIA 是从饱和密度附近往下看，那 SπRIT 就是从饱和密度的两倍往上走。

SπRIT 的全名是 SAMURAI Pion-Reconstruction and Ion-Tracker。它是一台时间投影室（TPC），装在日本理化学研究所 RIKEN 的放射性同位素束流工厂 RIBF 的 SAMURAI 谱仪磁铁里。它也可以当"活性靶"用。

技术参数：漂移长度 50.5 cm，pad 平面 86.4 cm × 134.4 cm，共 12096 个读出 pad，用 GET（Generic Electronics for TPCs）电子学读出。

装置的设计论文有两代：

  • **R. Shane 等，*Nuclear Instruments and Methods in Physics Research A*（2015）**，提出与设计阶段。
  • **J. Barney 等，*Review of Scientific Instruments* 92, 063302 (2021)，DOI 10.1063/5.0041191，预印本 arXiv:2005.10806。这是完整的装置论文，作者 39 人，署名机构包括密歇根州立大学 NSCL、RIKEN 仁科中心、京都大学、韩国高丽大学、德州 A&M、西密歇根大学、波兰科学院核物理研究所、意大利 Gran Sasso 国家实验室、利物浦大学、清华大学的 Z. G. Xiao。注意这里容易搞错：有人会把它记成发表在 Nuclear Instruments and Methods A 上，实际是 Review of Scientific Instruments。**（arXiv 检索页面里那一行"Nuclear Inst. and Methods A 986 (2021) 164713"属于另一篇不相关的论文，是列表串行导致的误读。）

它测什么？测π⁻ 与 π⁺ 的谱比。反应体系是 ¹³²Sn+¹²⁴Sn、¹¹²Sn+¹²⁴Sn 和 ¹⁰⁸Sn+¹¹²Sn，入射能 E/A = 270 MeV。这个能量下碰撞会压缩到约 2 倍饱和密度，而 π 介子主要产生在高密区，π⁻/π⁺ 的比值因此对超饱和密度处的对称能敏感。

两篇关键结果，方向相反：

**第一篇：Estee 等，《Probing the Symmetry Energy with the Spectral Pion Ratio》，*Physical Review Letters* 126, 162701 (2021)，预印本 arXiv:2103.06861。第一作者 J. Estee，作者列表很长（W. G. Lynch、C. Y. Tsang、M. B. Tsang、J. Barney、G. Jhang、T. Isobe、M. Kurata-Nishimura、B. Hong、Y. Leifels、J. Łukasik、A. B. McIntosh、H. Sakurai、S. J. Yennello、Z. G. Xiao、Y. Zhang 等，署名 SπRIT collaboration；理论合作者 M. D. Cozma 单列）。结果给出 42 < L < 117 MeV**。

**第二篇：Jhang 等，《Symmetry energy investigation with pion production from Sn+Sn systems》，*Physics Letters B* 813, 136016 (2021)，预印本 arXiv:2012.06976。这一篇的结论要老实说：把数据与七个主流输运模型对照，模型之间定性趋势一致，但定量的差异太大，无法可靠提取对称能约束**。文章自己指出，这或许能解释此前从 Au+Au 的 π 数据得出的相互矛盾的对称能结论。

这两篇放在一起，恰好说明了这个领域最真实的困难：同一批高质量数据，喂给不同的输运模型，会得到不同的对称能。

后面有两条理论线沿着这批数据往下走，而且结论分得很开：

其一：Gao-Chan Yong（中科院近代物理研究所，兰州；中国科学院大学核科学与技术学院）。 论文《Symmetry energy extracted from the SπRIT pion data in Sn+Sn systems》，*Physical Review C* 104, 014613 (2021)，DOI 10.1103/PhysRevC.104.014613，预印本 arXiv:2104.12098。他用改进的 impIBUU 模型（包含了核子-核子短程关联），在 270 MeV/nucleon 下算 Sn+Sn 的 π⁻、π⁺ 产额与比值，发现当采用软的对称能、L(ρ₀) = 66.75 ± 24.75 MeV 时能与 SπRIT 数据吻合得很好。文章还特别强调，初始核子动量分布的高动量尾（HMT）强烈影响阈值附近的 π 产生。

其二：Liu Hengjin、Zhang Ban、Guo Yuan-Qing、Chen Si-Min、Feng Zhaoqing（华南理工大学物理与光电学院）。 论文《Shedding light on the pion production in heavy-ion collisions and applications to the neutron star matter properties》，*Journal of Physics G: Nuclear and Particle Physics*，2026 年 8 月 27 日接收，DOI 10.1088/1361-6471/ae9f77。他们用兰州量子分子动力学（LQMD）输运模型，第一次把共振衰变以及 π 介子再吸收过程中的能量守恒、以及非弹核子-核子和核子-共振碰撞都考虑进去，分析了 ¹³²Sn+¹²⁴Sn 与 ¹⁰⁸Sn+¹¹²Sn 在 270 MeV/nucleon 下的 π⁻/π⁺ 动能谱，按 0.2ρ₀–0.8ρ₀ 和 1.2ρ₀–1.8ρ₀ 两个密度区分别追踪同位旋扩散。从 SπRIT 数据提取的结果是 L(ρ₀) = 42 ± 25 MeV（1σ 标准误差分析），并由此给出最大质量 2 倍太阳质量、半径 11 到 13 公里的中子星。

66.75 对 42，两个数差了 25 MeV。这不是谁算错了，这就是当前模型依赖的真实体量。

---

十、RAON：得说清楚它现在在哪儿

原始材料把 INDRA-FAZIA、SπRIT 和韩国的 RAON 并列为"主力平台"。前两个成立，第三个需要更正。

RAON 的全名是 Rare isotope Accelerator complex for ON-line experiments，由韩国基础科学研究院 IBS 旗下的稀有同位素科学研究所 IRIS 建设，2011 年立项。

它的超导直线加速器分低能段 SCL3 和高能段 SCL2。

SCL3 已建成并完成调试。 调试过程是：先让氩束通过 22 个四分之一波长谐振腔（QWR）模块加速到 2.5 MeV/u，再通过 32 个半波长谐振腔（HWR）模块加速到 18 MeV/u，然后把氩束送到 KoBRA（Korea Broad Acceptance Recoil spectrometer and Apparatus）做装置调试。ISOL 系统用质子束打 SiC 靶，产生并鉴别了 Na、Al 等同位素，随后把 ²⁵Na 加速到 16 MeV/u 送到 KoBRA；用 LaC₂ 靶则产生并鉴别了 Cs、Ba 的同位素。共线激光谱学系统 CLaSsy 用 ²¹Na 和 ²²Na 束流完成调试。NDPS（核数据生产系统）和 MR-TOF（多次反射飞行时间质量测量系统）仍在调试中。2024 年 3 月召开了首次面向韩国国内用户的 PAC 会议，2024 年夏天开始用户服务。

SCL2 仍在研制中，设计目标是加速重离子到 200 MeV/u（²³⁸U⁷⁹⁺，电荷态在剥离后分布在 +77 到 +81，平均 +79）。

RAON 的第一篇物理结果发表于 2026 年 6 月 8 日，发表在 *Results in Physics* 上。内容是：用 RAON 加速的 ⁴⁰Ar 束流打氢靶，在库仑位垒附近的低能区测量 ⁴⁰Ar + p 弹性散射。实验把束流能量分别设为 4.4、5.9、8.3 MeV/u，用 IBS 稀有核研究中心（CENS）自研的 ELARK 探测器测量散射质子的角度和能量。研究团队由 CENS 主任 Han In-sik 领导，通讯作者是 An Sung-hoon。结论是：现有的全局光学势无法正确复现这次测到的低能区散射数据，团队据此导出了一个适用于该能区的光学势。

这个事实必须写清楚：RAON 目前做的还是库仑位垒附近的低能核反应数据，不是费米能区（20 到 100 MeV/nucleon）的重离子碰撞，更没有做对称能约束。 它的 KoBRA、NDPS、CLaSsy、MR-TOF 分别面向核反应、核数据、激光谱学和精密质量测量。要等 SCL2 建成并达到设计指标，它才有可能进入费米能区这一战场。把它与已经产出对称能约束的 INDRA-FAZIA 并列，是提前兑现了一张还没到期的支票。

---

十一、K₀ 到底是多少？——一个查不到出处的数字

现在处理最硬的一处。

原始材料给出的依据是"K₀ = 249 ± 29 MeV"。经过多轮定向检索——中文"不可压缩系数 K₀ 249"、英文 "K0 = 249"、"249 ± 29"、"K_0 = 249"、arXiv 全文检索、INSPIRE 检索、会议论文集口径——没有找到任何文献使用这个数值口径。

能查到的、可引用的 K₀ 谱系是这样的：

方法/来源K₀ 取值说明
巨单极共振（GMR）主流拟合240 ± 20 MeV长期被认为的"公认值"
Stone, Stone, Moszkowski, *Phys. Rev. C* 89, 044316 (2014)250 < K₀ < 315 MeV放开表面/体积系数比 c = K_surf/K_vol 后得到更宽区间
相对论点耦合 DD-PCX 泛函213.03 ± 3.54 MeV用 ²⁰⁸Pb 同位旋偶极与同位旋标量单极激发联合优化
FOPI 椭圆流（IQMD）190 ± 30 MeVLe Fèvre 等
FOPI 椭圆流（UrQMD）220 ± 40 MeVWang 等
Huth 等采用的合并流约束200 ± 25 MeV见下文 Nature 2022
Cozma，dcQMD 输运模型230 ⁺⁹₋₁₁ MeV集体流与阻止本领联合
Xie 与 Bao-An Li，贝叶斯192 ⁺¹²₋₁₆ MeV高能重离子框架
低能核反应融合阻挫（BUU）240–260 MeV用冷热熔合反应概率导出
Sun, Liu, Li, Xiao, Zhang，*中国科学·物理学 力学 天文学* 56(7) (2026)250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV集体流 + π 产生 + 中子星观测的贝叶斯联合

其中最接近 249 ± 29 的，是最后一行的 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV。中心值差 1.56 MeV，但误差结构完全不同（−18.79/+16.51 是非对称区间，不是 ±29）。

结论有两条，都得说：

第一，249 ± 29 MeV 这个口径在当前文献里不可引用。它不是 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ 的等价写法（误差差了一倍多），也不是任何其他已知结果的转写。

第二，如果确实需要一个"最新、综合、可引用"的 K₀ 锚点，应该用 Cozma 的 230 ⁺⁹₋₁₁ MeV（纯地面重离子数据）或 李昂组的 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV（地面加天体的贝叶斯联合），并注明出处与置信区间口径。

还有一处需要点出的细节：K₀ 本身与对称能是两件事。 K₀ 描述的是对称核物质（δ = 0）的曲率，S(ρ) 描述的是同位旋不对称那部分的密度依赖。INDRA-FAZIA 系列文章给的是 S 的约束，不是 K₀ 的约束。"用 INDRA-FAZIA 测了 K₀ = 249"这种说法在逻辑上就是串了线。K₀ 的现代约束主要来自巨单极共振、集体流和停止本领。

---

十二、把地面上的数据和天上的数据拼起来

费米能区的碰撞和对称能之所以值得做，最终目的是约束中子星。这几年最有力的进展，是把两边的数据放进同一套贝叶斯推断里。

**第一篇：Huth 等，《Constraining neutron-star matter with microscopic and macroscopic collisions》，*Nature* 606(7913), 276–280 (2022)**，DOI 10.1038/s41586-022-04750-w，预印本 arXiv:2107.06229，2022 年 6 月 9 日发表。

作者 11 人：S. Huth、P. T. H. Pang、I. Tews、T. Dietrich、A. Le Fèvre、A. Schwenk、W. Trautmann、K. Agarwal、M. Bulla、M. W. Coughlin、C. Van Den Broeck。署名单位包括德国马普引力物理研究所（波茨坦 AEI-Golm）、GSI 亥姆霍兹重离子研究中心等。

做法是把三类输入用贝叶斯推断揉在一起：天体多信使观测（中子星质量半径、引力波潮汐形变度）、相对论能区金核碰撞数据、以及微观核理论计算。结论有两个：第一，加入重离子碰撞数据之后，致密物质的压强相对此前分析上升了；第二，这个压强上升把中子星半径往更大的方向推，与 NICER 的观测相符。文中给出的一个代表性约束是 K = 200 ± 25 MeV。

这篇文章的历史地位在于：它是第一次把重离子碰撞和天体观测放进同一个贝叶斯框架，并证明两者高度自洽，还能互相补盲——重离子管中间密度，引力波和 X 射线管高密度。

**第二篇：C. Y. Tsang、M. B. Tsang、W. G. Lynch、R. Kumar、C. J. Horowitz，《Determination of the equation of state from nuclear experiments and neutron star observations》，*Nature Astronomy* 8(3), 328–336 (2024)**，DOI 10.1038/s41550-023-02161-z，2024 年 1 月 5 日在线发表，3 月刊出。作者五人，含印第安纳大学、密歇根州立大学 NSCL 等。

这篇的特点是"数据中心式"：用12 项核实验约束 + 3 项天文观测，在贝叶斯框架下合成一套完整的核物态方程，并给出对核子自由度、中子星组分和以中微子辐射冷却的推论。

第三篇（国内）：Sun Xiangdong、Liu Yangyang、Li Ang、Xiao Zhigang、Zhang Yingxun，《Astronomical and terrestrial studies on the equation of state of dense matter》，*Scientia Sinica Physica, Mechanica & Astronomica*（《中国科学：物理学 力学 天文学》）56(7)，2026 年，DOI 10.1360/SSPMA-2025-0201。

署名单位为厦门大学天文学系、中山大学等。作者阵容值得记一下：李昂（厦门大学）与张英逊（中国原子能科学研究院）为通讯作者，肖志刚（当时在湖州师范学院/中山大学方向）也在列，第一作者是 Sun Xiangdong。

做法是先用地面上核实验提取的核子相互作用信息，结合核多体理论与非平衡输运模型，再引入高精度中子星观测，做自洽联合分析。纳入的数据包括重离子碰撞的集体流与 π 介子产生数据，以及中子星的质量、半径、潮汐形变度。用贝叶斯推断得到的 1σ 区间是：

  • K₀ = 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV
  • **m*/m = 0.6 ⁺⁰·⁰¹₋₀·₀₁**（核子有效质量比）
  • S₀ = 28.46 ⁺⁴·³⁶₋₂·₂₄ MeV（饱和密度处对称能）
  • L₀ = 64.59 ⁺¹³·²⁵₋₉·³⁸ MeV（对称能斜率）

这组数字是目前能看到的最新一版"地面加天上"的中国团队结果，也是本文建议替代 249 ± 29 的锚点之一。

---

十三、中子星那一侧的刻度在变精

交叉约束的另一端是天文观测本身。这一侧的进展直接决定了联合分析的紧度。

最新的一把尺子是 PSR J0437−4715。

**Miller 等，《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。

作者 15 人：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. Morsink、P. S. Ray、R. A. Remillard、Z. Wadiasingh、M. T. Wolff。第一作者在马里兰大学，合作单位包括NASA 戈达德太空飞行中心、伊利诺伊大学香槟分校、Haverford College 等。

这一篇的技术要点很值得说：此前 NuSTAR 观测已确认这颗脉冲星有一个被调制的非热成分，但早先对 NICER 数据的分析没有把它建模进去。这次加上之后，贝叶斯证据显著提高。团队采用的头部模型是"三个等温圆形热斑 + 调制非热成分"，推出在独立测得的质量 M = 1.418 ± 0.044 倍太阳质量下，半径的 68% 置信区间为 11.9 到 15.5 km，与同质量量级的 PSR J0030+0451 此前结果一致。

中文界还有一条平行工作值得一提：Qi Liqiang、Zhang Juan、Xu Weiwei、Zheng Shijie、Ge Mingyu、Li Ang、Zhang Shuang-Nan、Feng Hua、Lu Fangjun，《Multiwavelength Analysis of PSR J0437−4715 with Pulse Profile Modeling》，预印本 arXiv:2605.05429。

署名单位是中国科学院高能物理研究所粒子天体物理重点实验室（第一作者 Qi Liqiang，通讯作者 Ge Mingyu）、中国科学院国家天文台、北京师范大学、中国科学院大学，以及厦门大学天文学系（Li Ang）。他们用哈勃望远镜的远紫外、ROSAT 的软 X 射线和 XMM-Newton 的 X 射线数据做多波段联合建模，配合射电偏振位置角给出的热斑几何先验，得到 M = 1.38 ± 0.03 倍太阳质量、赤道周向半径 R = 13.25 ⁺⁰·³⁴₋₀·₃₅ km（68% 置信区间）。文中指出，相比只用 NICER 的拟合，多波段数据把半径后验分布往更大的方向推。

两个半径不一样（一个给区间 11.9–15.5 km，一个给 13.25 km 的点估计），这不是矛盾——它们的建模假设、数据组合和功率区间口径都不同。读这类结果必须带上建模条件。

---

十四、剩下的模型依赖有多大：一篇 2026 年的综述怎么说

2026 年 5 月出现了一篇专门总结这条技术路线的综述，它的判断很值得引用。

S. Mallik、F. Gulminelli、C. Ciampi、D. Gruyer，《Probing the density dependence of nuclear symmetry energy through isospin transport in heavy-ion reactions》。预印本 arXiv:2605.06138，2026 年 5 月 7 日；正式版发表在 Universe 2026, 12, 202（对应 Universe 期刊"第 10 周年·从中能到超高能核反应中的强相互作用物质"专题）。

作者署名机构要写准：C. Ciampi 在 GANIL（法国国家重离子加速器实验室，CEA-DRF 与 CNRS-IN2P3 联合单位）；S. Mallik 在加尔各答 Variable Energy Cyclotron Centre，同时属孟买 Homi Bhabha National Institute；F. Gulminelli 与 D. Gruyer 在 LPC Caen（诺曼底大学、ENSICAEN、UNICAEN、CNRS-IN2P3 联合实验室）。法国国家开放档案编号 hal-05627602。

它明确的立场是：要从重离子反应里拿到稳健的对称能约束，观测量必须满足两个条件——对末态相互作用尽量不敏感，并且被输运模型可靠地预言、模型依赖有限。它把 INDRA-FAZIA 的同位旋输运比数据与 BUU 计算、连同从头算与唯象两类现代核泛函一起梳理，给出对称能的置信区域，并特别强调把敏感密度区定位准确这件事的价值。

它同时坦白了一件事，也是这个领域现在最诚实的自我评价：目前的约束仍然建立在单一 BUU 输运框架之上。 换句话说，S = 29.0 ± 0.7 MeV 这个极其漂亮的误差棒，只统计了实验误差与泛函散布，没有包含输运模型本身的不确定性。

而这恰恰是 SπRIT 那边最痛的地方——七个模型的散布大到无法提取约束。两边合起来看，就是当前这个领域的真实处境：实验精度已经不是瓶颈，理论框架的差异才是。

---

十五、前沿还在往哪走

观测侧：

INDRA-FAZIA 会继续从 E818 数据里挖。它瞄准的是准弹核的气化与类气化事件、更重团簇的介质内效应、以及用动力学温度、激发态温度和同位素温度三种"温度计"同时测温。这些工作的目标不是再报一个 S 值，而是把探测密度区间往更低处推——因为越靠低密度，对称能的行为越不确定，而低密度那一段恰好是中子星壳层和超新星物质最关心的。

装置侧：

FAZIA 的同位素分辨能力（Z 到约 25）和 INDRA 的角覆盖（90%）这个组合目前没有对手。下一步的竞争点在于能不能用放射性束流做同样的测量——用中子极端富集的束流可以放大同位旋梯度，让 R_i 的信号更明显。这需要 FAZIA 适应更高强度、更复杂的束流条件。

实验室侧：

  • RAON：等 SCL2 建成。这是唯一一条明确指向"放射性束 + 费米能区"的新路线。
  • HIAF（中国惠州）：2026 年 7 月 21 日通过工艺验收、进入试运行。它是世界上第一台把超导直线加速器与快循环同步加速器集成在一起的重离子研究装置。调试中氧束达到每脉冲 2.5×10¹¹ 个粒子、最高能量 4.25 GeV，铋束达到每脉冲 3.2×10¹⁰ 个粒子，均创同类装置记录。调试实验里还观测到两个新同位素，短寿命核素 Au-202 的质量精度比原世界记录提高一倍。HIAF 的六个实验终端面向"核存在极限、核天体物理过程、重离子交叉应用"三大目标——目前它还不是 EOS 约束的主力，但它的高流强特性使它未来具备做费米能区联动测量的条件。
  • FRIB、RIBF：继续提供中子富集束流。

方法与理论侧：

这一侧的变化可能最大。

一是机器学习进入这条链条。湖州师范学院与中科院近物所联合团队（Yongjia Wang、Zepeng Gao、Hongliang Lü、Qingfeng Li）用 LightGBM 做逐事件的对称能反演，论文《Decoding the nuclear symmetry energy event-by-event in heavy-ion collisions with machine learning》发表在 *Physics Letters B* 835, 137508 (2022)。传统做法只用观测量在所有事件上的期望值，而实验和模拟其实都能给出逐事件的观测量——这一篇把逐事件信息用起来了，从 UrQMD 测试数据里提取的对称能斜率参数平均散布约 30 MeV，并对模型参数变化稳健。

顺带说，这个组的另两项工作也值得记：Yongjia Wang、Qingfeng Li 与 GSI 的 Yvonne Leifels、Arnaud Le Fèvre 合作，用 UrQMD 加 FOPI 质子椭圆流数据，把对称能斜率约束到 L₀ = 43 ± 20 MeV（95% 置信水平），发表在 *Physics Letters B* 802, 135249 (2020)；以及 Quantum Li（Qingfeng Li）早期关于重建判据对对称势敏感探针影响的分析（arXiv:0808.0624），署名 FIAS 法兰克福高等研究院与湖州师范学院。

二是用神经网络重建碰撞参数。这条路要和研究第三节讲的"模型无关重建"严格区分开——它们是两件事。

  • Frankland 等 2021 年的模型无关重建，是概率反演方法：条件概率 + gamma 分布涨落核 + 贝叶斯定理，不含任何机器学习。它已经用在真实 INDRA-FAZIA 数据上，是 PRC 111, 044601 的地基。
  • Valente 与 Lopez（LPC Caen）的工作是在模拟上训练神经网络（HIPSE 与 ELIE 模拟生成训练集），达到亚飞米级别的碰撞参数精度，再配合贝叶斯推断去估计介质中的核子-核子截面和达到的最大密度。相关结果见 *Il Nuovo Cimento C* 48(2), 1–7 (2025)，另有 EuNPC 2025 会议海报。

把这两者说成同一件事，或者把神经网络那一条说成"INDRA-FAZIA 用的方法"，都是错的。

三是多模型对照。既然单一 BUU 框架的残余模型依赖无法自证，那就必须做跨代码的横向比对——这正是 SπRIT 那篇 PLB 813, 136016 用七个模型做的事情，也是 Sorensen 等人那份 132 位以上作者的白皮书《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，预印本 arXiv:2301.13253）花大力气做的事。

---

十六、一页总结

把这篇的脉络压成几句话：

问题的形状。 核物质的物态方程有两个关键数：K₀（对称核物质的硬度）和 S(ρ)（中子富集带来的额外能量，随密度变化）。费米能区（每核子 20 到 100 MeV）的重离子碰撞，是目前唯一能在实验室里同时压出密度变化、又保住同位素信息的手段。

方法的链条。 一次碰撞要变成一条约束，必须走完五步：一，用模型无关的概率反演从多重数反推碰撞参数 b（Frankland 等，*Phys. Rev. C* 104, 034609 (2021)）；二，用同位旋输运比 R_i 把对称参考体系和不称体系对照起来（谱系到 Rami 2000 与 Tsang 2004）；三，把 R_i 随 b 的变化与 BUU 输运计算对照；四，从重子密度与同位旋流密度的时间演化里定出真正被探到的密度区间；五，给出对称能随密度的置信区域。

最新的结果。 Ciampi 等（INDRA-FAZIA 合作组，34 位作者）用 2019 年 E789 的 ⁵⁸·⁶⁴Ni + ⁵⁸·⁶⁴Ni 数据，在 *Phys. Rev. C* 111, 044601 (2025) 完成模型无关测量，在 *Phys. Lett. B* 868, 139815 (2025) 给出 S = 29.0 ± 0.7 MeV @ ρ/ρ₀ = 1.01。注意敏感密度是接近饱和密度的 1.01ρ₀，不是深亚饱和。

另一个方向。 SπRIT（RIKEN RIBF）测 π⁻/π⁺ 谱比，探的是约 2ρ₀ 的超饱和密度。Estee 等给出 42 < L < 117 MeV（*Phys. Rev. Lett.* 126, 162701）；Jhang 等用七个输运模型对照后判定无法可靠提取约束（*Phys. Lett. B* 813, 136016）。后续理论解读分化：Gao-Chan Yong 的 impIBUU 得 L = 66.75 ± 24.75 MeV（*Phys. Rev. C* 104, 014613），华南理工 Feng Zhaoqing 组的 LQMD 得 L = 42 ± 25 MeV（*J. Phys. G*，2026）。

三处必须更正的说法。

  1. K₀ = 249 ± 29 MeV 查无出处，不可引用。可用的替代锚点是 Cozma 的 230 ⁺⁹₋₁₁ MeV，或李昂组 2026 年贝叶斯联合的 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV。另外，INDRA-FAZIA 系列约束的是 S(ρ) 不是 K₀，不能串线。
  2. RAON 目前不是费米能区 EOS 研究的主力平台。SCL3 只到 18 MeV/u，SCL2（设计 200 MeV/u）仍在研制。它的第一篇物理成果是 2026 年 6 月发表的 ⁴⁰Ar + p 低能弹性散射。
  3. SπRIT 与 INDRA-FAZIA 探的是相反方向（超饱和 vs 饱和密度附近偏下），不能当成"三台机器做同一件事"；机器学习碰撞参数重建（Valente 与 Lopez）与模型无关概率反演（Frankland 等）也不是同一件事。

最真实的困难。 综述作者自己点破：现有约束仍建立在单一 BUU 框架上，所以那个 ±0.7 MeV 的误差棒没有包含输运模型的不确定性。而 SπRIT 那边七个模型的散布大到无法提取约束。把两边放在一起看，这个领域现在卡在的地方不是探测器，是理论框架之间还没有收敛。

---

参考文献

以下 38 条是本文全部引用的汇总索引。正文中改为随文写出作者、期刊、卷期与 DOI 的完整信息，便于直接核对；本表按主题分组排列（INDRA-FAZIA 主线在前，SπRIT 主线随后，K₀ 谱系、中子星交叉约束、装置与国内团队工作依次在后），不按正文出现顺序编号。

  1. C. Ciampi, J. D. Frankland, D. Gruyer, N. Le Neindre, S. Mallik, R. Bougault, A. Chbihi, L. Baldesi, S. Barlini, E. Bonnet, B. Borderie, A. Camaiani, G. Casini, I. Dekhissi, D. Dell'Aquila, J. A. Dueñas, Q. Fable, F. Gramegna, C. Gouyet, M. Henri, B. Hong, S. Kim 等（INDRA-FAZIA Collaboration，共 34 位作者），"Model-independent measurement of isospin diffusion in Ni-Ni systems at intermediate energy"，*Physical Review C* 111, 044601 (2025)。DOI: 10.1103/PhysRevC.111.044601；arXiv:2412.13648。第一作者单位：GANIL（法国卡昂，CEA-DRF / CNRS-IN2P3）。
  1. C. Ciampi, S. Mallik, F. Gulminelli, D. Gruyer, J. D. Frankland 等（共 34 位作者），"Constraining the nuclear symmetry energy with Fermi-energy heavy ion collisions"，*Physics Letters B* 868, 139815 (2025)。DOI: 10.1016/j.physletb.2025.139815。作者单位：C. Ciampi（GANIL）；S. Mallik（Variable Energy Cyclotron Centre, Calcutta 与 Homi Bhabha National Institute, Mumbai）；F. Gulminelli、D. Gruyer（LPC Caen）；J. D. Frankland（GANIL）。核心结果：S = 29.0 ± 0.7 MeV @ ρ/ρ₀ = 1.01。
  1. J. D. Frankland, D. Gruyer, E. Bonnet, B. Borderie, R. Bougault, A. Chbihi, J. E. Ducret, D. Durand, Q. Fable, M. Henri, J. Lemarié, N. Le Neindre, I. Lombardo, O. Lopez, L. Manduci, M. Pârlog, J. Quicray, G. Verde, E. Vient, M. Vigilante（INDRA Collaboration，共 20 位作者），"Model independent reconstruction of impact parameter distributions for intermediate energy heavy ion collisions"，*Physical Review C* 104, 034609 (2021)。DOI: 10.1103/PhysRevC.104.034609；arXiv:2011.04496。第一作者单位：GANIL。
  1. C. Ciampi 等（INDRA-FAZIA Collaboration），"First results from the INDRA-FAZIA apparatus on isospin diffusion in ⁵⁸·⁶⁴Ni + ⁵⁸·⁶⁴Ni collisions at 32 and 52 MeV/nucleon"，*Physical Review C* 106, 024603 (2022)。arXiv:2205.10171。
  1. C. Ciampi 等（INDRA-FAZIA Collaboration），*Physical Review C* 108, 054611 (2023)。INDRA-FAZIA 同位旋扩散分析的后续细化（准弹蒸发通道）。
  1. L. Baldesi 等（INDRA-FAZIA Collaboration），"Nuclear symmetry energy at work in heavy ion reactions"，*Journal of Physics: Conference Series* 2586, 012039 (2023)。DOI: 10.1088/1742-6596/2586/1/012039。结论：32 MeV/nucleon 的同位旋均衡强于 52 MeV/nucleon。
  1. J. D. Frankland, R. Bougault, G. Casini, D. Gruyer, N. Le Neindre 等，E818 实验（"Extending our knowledge of Warm Dense Nuclear Matter in the low density region"），GANIL 数据集，DOI: 10.26143/GANIL-2022-E818_20。执行时间 2022-05-07 至 2022-05-17；体系 ⁵⁸Ni+⁵⁸Ni 与 ³⁶Ar+⁵⁸Ni @ 74 MeV/nucleon；发言人 R. Bougault、G. Casini；数据负责人 J. D. Frankland、D. Gruyer。
  1. F. Rami, Y. Leifels, B. de Schauenburg, A. Gobbi, B. Hong 等（FOPI Collaboration，共 70 位作者），"Isospin tracing: A probe of nonequilibrium in central heavy-ion collisions"，*Physical Review Letters* 84(6), 1120–1123 (2000)。DOI: 10.1103/PhysRevLett.84.1120；arXiv:nucl-ex/9908014。体系：⁹⁶Ru + ⁹⁶Zr 四种弹靶组合 @ 400 A·MeV。主要单位：IReS（斯特拉斯堡）与 LPC（克莱蒙费朗）。注意：这是同位旋示踪概念的先声，不是归一化输运比公式本身；其能量也高于费米能区。
  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, W. A. Friedman, L. Beaulieu, B. Davin, R. T. de Souza, Y. Larochelle, T. Lefort, R. Yanez, V. E. Viola, R. J. Charity, L. G. Sobotka，"Isospin Diffusion and the Nuclear Symmetry Energy in Heavy Ion Reactions"，*Physical Review Letters* 92, 062701 (2004)。DOI: 10.1103/PhysRevLett.92.062701。第一作者单位：密歇根州立大学 NSCL。
  1. S. Mallik, F. Gulminelli, C. Ciampi, D. Gruyer，"Probing the density dependence of nuclear symmetry energy through isospin transport in heavy-ion reactions"，*Universe* 2026, 12, 202。arXiv:2605.06138（2026-05-07）；HAL: hal-05627602。作者单位：C. Ciampi（GANIL）；S. Mallik（VECC Calcutta 与 HBNI Mumbai）；F. Gulminelli、D. Gruyer（LPC Caen）。
  1. A. Sorensen, K. Agarwal, K. W. Brown, Z. Chajęcki, P. Danielewicz 等（132 位以上作者），"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；arXiv:2301.13253。
  1. S. Huth, P. T. H. Pang, I. Tews, T. Dietrich, A. Le Fèvre, A. Schwenk, W. Trautmann, K. Agarwal, M. Bulla, M. W. Coughlin, C. Van Den Broeck，"Constraining neutron-star matter with microscopic and macroscopic collisions"，*Nature* 606(7913), 276–280 (2022)。DOI: 10.1038/s41586-022-04750-w；arXiv:2107.06229。
  1. C. Y. Tsang, M. B. Tsang, W. G. Lynch, R. Kumar, C. J. Horowitz，"Determination of the equation of state from nuclear experiments and neutron star observations"，*Nature Astronomy* 8(3), 328–336 (2024)。DOI: 10.1038/s41550-023-02161-z。2024-01-05 在线。
  1. X. Sun, Y. Liu, A. Li, Z. Xiao, Y. Zhang，"Astronomical and terrestrial studies on the equation of state of dense matter"，*Scientia Sinica Physica, Mechanica & Astronomica*（《中国科学：物理学 力学 天文学》）56(7), 2026。DOI: 10.1360/SSPMA-2025-0201。通讯作者：李昂（厦门大学）、张英逊（中国原子能科学研究院）。结果：K₀ = 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV，m*/m = 0.6 ⁺⁰·⁰¹₋₀·₀₁，S₀ = 28.46 ⁺⁴·³⁶₋₂·₂₄ MeV，L₀ = 64.59 ⁺¹³·²⁵₋₉·₃₈ MeV。
  1. J. R. Stone, N. J. Stone, S. A. Moszkowski，"Incompressibility in finite nuclei and nuclear matter"，*Physical Review C* 89, 044316 (2014)。DOI: 10.1103/PhysRevC.89.044316；arXiv:1404.0744。作者单位：牛津大学物理系、田纳西大学、加州大学洛杉矶分校。
  1. M. D. Cozma，"Equation of state of nuclear matter from collective flows and stopping in intermediate-energy heavy-ion collisions"，arXiv:2407.16411。dcQMD 输运模型结果：K₀ = 230 ⁺⁹₋₁₁ MeV，L = 63 ⁺¹⁰₋₁₃ MeV，m*/m = 0.695 ⁺⁰·⁰¹⁴₋₀·₀₁₈。
  1. M. B. Tsang, J. R. Stone, F. Camera, P. Danielewicz, S. Gandolfi 等，"Constraints on the symmetry energy and neutron skins from experiments and theory"，*Physical Review C* 86, 015803 (2012)。DOI: 10.1103/PhysRevC.86.015803；arXiv:1204.0466。
  1. J. Estee, W. G. Lynch, C. Y. Tsang, J. Barney, G. Jhang, M. B. Tsang, R. Wang, M. Kaneko, J. W. Lee, T. Isobe, M. Kurata-Nishimura, T. Murakami, D. S. Ahn, H. Baba, Z. Chajecki, M. Famiano, B. Hong, Y. J. Kim, Y. Leifels, J. Łukasik, A. B. McIntosh, H. Otsu, P. Pawłowski, K. Pelczar, H. Sakurai, R. Shane, H. Simon, H. Suzuki, S. Tangwancharoen, S. J. Yennello, Z. G. Xiao, Y. Zhang 等（SπRIT collaboration）与 M. D. Cozma，"Probing the Symmetry Energy with the Spectral Pion Ratio"，*Physical Review Letters* 126, 162701 (2021)。arXiv:2103.06861。结果：42 < L < 117 MeV。
  1. G. Jhang, J. Estee, J. Barney, G. Cerizza, M. Kaneko, J. W. Lee, W. G. Lynch, T. Isobe, M. Kurata-Nishimura, T. Murakami, C. Y. Tsang, M. B. Tsang, R. Wang, H. Baba, Z. Chajecki, M. Famiano, B. Hong, T. Kobayashi, P. Lasko, J. Łukasik, N. Nakatsuka, R. Olsen, H. Otsu, P. Pawłowski, K. Pelczar, H. Sakurai, C. Santamaria, H. Setiawan, A. Taketani, J. R. Winkelbauer, Z. Xiao, S. J. Yennello, J. Yurkon（SπRIT collaboration），"Symmetry energy investigation with pion production from Sn+Sn systems"，*Physics Letters B* 813, 136016 (2021)。arXiv:2012.06976。结论：七个输运模型分散过大，无法可靠提取对称能约束。
  1. Gao-Chan Yong，"Symmetry energy extracted from the SπRIT pion data in Sn+Sn systems"，*Physical Review C* 104, 014613 (2021)。DOI: 10.1103/PhysRevC.104.014613；arXiv:2104.12098。作者单位：中国科学院近代物理研究所（兰州）与中国科学院大学核科学与技术学院。结果：impIBUU，L(ρ₀) = 66.75 ± 24.75 MeV。
  1. Hengjin Liu, Ban Zhang, Yuan-Qing Guo, Si-Min Chen, Zhaoqing Feng，"Shedding light on the pion production in heavy-ion collisions and applications to the neutron star matter properties"，*Journal of Physics G: Nuclear and Particle Physics*，2026 年 8 月 27 日接收。DOI: 10.1088/1361-6471/ae9f77。作者单位：华南理工大学物理与光电学院（通讯作者 Z. Feng）。结果：LQMD 模型，L(ρ₀) = 42 ± 25 MeV。
  1. J. Barney, J. Estee, W. G. Lynch, T. Isobe, G. Jhang, M. Kurata-Nishimura, A. B. McIntosh, T. Murakami, R. Shane, S. Tangwancharoen, M. B. Tsang, G. Cerizza, M. Kaneko, J. W. Lee, C. Y. Tsang, R. Wang, C. Anderson, H. Baba, Z. Chajecki, M. Famiano, R. Hodges-Showalter, B. Hong, T. Kobayashi, P. Lasko, J. Łukasik, N. Nakatsuka, R. Olsen, H. Otsu, P. Pawłowski, K. Pelczar, H. Sakurai, C. Santamaria, H. Setiawan, A. Taketani, J. R. Winkelbauer, Z. Xiao, S. J. Yennello, J. Yurkon, Y. Zhang，"The SπRIT time projection chamber"，*Review of Scientific Instruments* 92, 063302 (2021)。DOI: 10.1063/5.0041191；arXiv:2005.10806。署名机构：密歇根州立大学 NSCL、RIKEN 仁科中心、京都大学、高丽大学、德州 A&M 大学回旋加速器研究所、西密歇根大学、波兰科学院核物理研究所、意大利 Gran Sasso 国家实验室、利物浦大学、清华大学。
  1. R. Shane, A. B. McIntosh, T. Isobe, W. G. Lynch, H. Baba, J. Barney, Z. Chajecki, M. Chartier, J. Estee, M. Famiano, B. Hong, K. Ieki, G. Jhang, R. Lemmon, F. Lu, T. Murakami, N. Nakatsuka, M. Nishimura, R. Olsen, W. Powell, H. Sakurai, A. Taketani, S. Tangwancharoen, M. B. Tsang, T. Usukura, R. Wang, S. J. Yennello, J. Yurkon，"SπRIT: A time-projection chamber for symmetry-energy studies"，*Nuclear Instruments and Methods in Physics Research A* 784, 513–517 (2015)。DOI: 10.1016/j.nima.2015.01.026；arXiv:1409.6343。署名单位含密歇根州立大学 NSCL、德州 A&M 回旋加速器研究所、RIKEN 仁科中心、利物浦大学、西密歇根大学、高丽大学、立教大学、STFC Daresbury 实验室、中国科学院上海应用物理研究所、京都大学。文中明确说明设计目标是超饱和密度（约 2ρ₀）的对称能。
  1. 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。结果：M = 1.418 ± 0.044 M☉，R 的 68% 置信区间 11.9–15.5 km。
  1. Liqiang Qi, Juan Zhang, Weiwei Xu, Shijie Zheng, Mingyu Ge, Ang Li, Shuang-Nan Zhang, Hua Feng, Fangjun Lu，"Multiwavelength Analysis of PSR J0437−4715 with Pulse Profile Modeling"，arXiv:2605.05429（2026）。作者单位：中国科学院高能物理研究所粒子天体物理重点实验室、中国科学院国家天文台、北京师范大学、中国科学院大学、厦门大学天文学系。结果：M = 1.38 ± 0.03 M☉，R = 13.25 ⁺⁰·³⁴₋₀·₃₅ km。
  1. Yongjia Wang, Zepeng Gao, Hongliang Lü, Qingfeng Li，"Decoding the nuclear symmetry energy event-by-event in heavy-ion collisions with machine learning"，*Physics Letters B* 835, 137508 (2022)。作者单位：湖州师范学院理学院、中山大学中法核工程与技术学院、华为海思研究部、中国科学院近代物理研究所。
  1. Yongjia Wang, Qingfeng Li, Yvonne Leifels, Arnaud Le Fèvre，"Study of the nuclear symmetry energy from the rapidity-dependent elliptic flow in heavy-ion collisions around 1 GeV/nucleon regime"，*Physics Letters B* 802, 135249 (2020)。DOI: 10.1016/j.physletb.2020.135249。结果：L₀ = 43 ± 20 MeV（95% 置信水平）。
  1. S. Valente, O. Lopez，"Toward a better characterization of the nuclear equation of state with a neural network for impact parameter determination"，*Il Nuovo Cimento C* 48(2), 1–7 (2025)。作者单位：LPC Caen（诺曼底大学、ENSICAEN、UNICAEN、CNRS-IN2P3）。另有 EuNPC 2025 会议报告。
  1. L. Baldesi 等（INDRA-FAZIA Collaboration），"The FAZIA apparatus: status and perspectives"，WPCF 2024 会议报告。内容含 INDRA-FAZIA 耦合方式（INDRA 前 5 环换成 12 个 FAZIA 模块）、E789 与 E818 实验清单、同位旋输运比定义（并明确溯源到 F. Rami 等 PRL 84, 1120 (2000)）。
  1. P. Danielewicz, R. Lacey, W. G. Lynch，"Determination of the equation of state of dense matter"，*Science* 298, 1592–1596 (2002)。DOI: 10.1126/science.1078070；arXiv:nucl-th/0208016。
  1. Bao-An Li, Lie-Wen Chen, Che Ming Ko，"Recent Progress and New Challenges in Isospin Physics with Heavy-Ion Reactions"，*Physics Reports* 464, 113–281 (2008)。DOI: 10.1016/j.physrep.2008.04.005；arXiv:0804.3580。
  1. J. M. Lattimer, M. Prakash，"Neutron Star Structure and the Equation of State"，*The Astrophysical Journal* 550, 426 (2001)。DOI: 10.1086/319702；arXiv:astro-ph/0002232。
  1. 中国科学院高能物理研究所/自然科学基金委相关公开报道，"China's High Intensity Heavy-ion Accelerator Facility Completes Construction, Begins Trial Operations"，中国科学院英文网站 2026 年 7 月 17 日发布（HIAF 2026-07-21 通过工艺验收并进入试运行；氧束 2.5×10¹¹ ppp、最高能量 4.25 GeV；铋束 3.2×10¹⁰ ppp）。
  1. Seung-Woo Hong（Institute for Rare Isotope Science, IBS），"Status of the RAON Facility"，IBS 会议报告。内容含 SCL3 调试（22 个 QWR 至 2.5 MeV/u，32 个 HWR 至 18 MeV/u）、²⁵Na 加速至 16 MeV/u 送 KoBRA、CLaSsy 用 ²¹Na/²²Na 调试、SCL2 研制中（设计 200 MeV/u）。
  1. IBS 与 IRIS 联合研究团队（CENS 主任 Han In-sik；通讯作者 An Sung-hoon），⁴⁰Ar + p 弹性散射测量结果，*Results in Physics*，2026 年 6 月 8 日发表。RAON 首篇物理成果；束流能量 4.4、5.9、8.3 MeV/u；使用 KoBRA 与 ELARK 探测器。
  1. S. Barlini, E. Bonnet, B. Borderie, R. Bougault, A. M. Buccola, A. Camaiani, G. Casini, A. Chbihi, C. Ciampi, J. Duenas, D. Fabris, J. Frankland, C. Frosin, F. Gramegna, D. Gruyer, F. Gulminelli, M. Henri, Hong Byungsik, Kim Sunji, N. Le Neindre, I. Lombardo, O. Lopez, G. Mantovani, T. Marchi, Nam Seonho, P. Ottanelli, H. Pais, M. Parlog, G. Pasquali, S. Piantelli, G. Poggi, C. Providencia, J. Quicray, A. Stefanini, Sun Baohua, A.-M. Sanchez-Benitez, S. Upadhyaya, S. Valdre, G. Verde, E. Vient（INDRA-FAZIA Collaboration），E818 实验科学方案（轻团簇物态方程与介质内修正），FAZIA 合作组官方发布页。
  1. Esra Yüksel, Tomohiro Oishi, Nils Paar，"Nuclear Equation of State in the Relativistic Point-Coupling Model Constrained by Excitations in Finite Nuclei"，*Universe* 7(3), 71 (2021)。DOI: 10.3390/universe7030071。作者单位：伊尔迪兹技术大学（伊斯坦布尔）、新潟大学、萨格勒布大学。文中给出 DD-PCX 泛函在纳入 ²⁰⁸Pb 的同位旋标量巨单极共振能量后，K₀ = 213.03 ± 3.54 MeV、J = 31.12 ± 0.32 MeV、L = 46.32 ± 1.68 MeV。
  1. M. Veselsky, J. Klimo, Yu-Gang Ma, G. A. Souliotis，"Constraining the equation of state of nuclear matter from fusion hindrance in reactions leading to the production of superheavy elements"，*Physical Review C* 94(6), 064608 (2016)。DOI: 10.1103/PhysRevC.94.064608；arXiv:1604.01963。第一作者单位：斯洛伐克科学院物理研究所（布拉迪斯拉发）；合作者含中国科学院上海应用物理研究所的 Yu-Gang Ma。结论：由实验熔合概率导出 K₀ = 240–260 MeV，同时排除过硬的对称能密度依赖。

---

核实纠错表

以下逐条列出原始材料中需要更正或必须注意说法的地方。

  1. K₀ = 249 ± 29 MeV 查无出处。 经中文与英文多轮定向检索、arXiv 全文、INSPIRE 与会议论文集口径核对，未发现任何文献使用该数值口径。最接近的是 Sun/Liu/Li/Xiao/Zhang 在《中国科学：物理学 力学 天文学》56(7) (2026) 给出的 K₀ = 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV——中心值差 1.56 MeV，但误差结构完全不同（非对称区间，不是 ±29）。建议不引用 249 ± 29，改用 Cozma 的 230 ⁺⁹₋₁₁ MeV 或李昂组的 250.56 ⁺¹⁶·⁵¹₋₁₈·₇₉ MeV。
  1. "用 INDRA-FAZIA 测 K₀ = 249"在逻辑上串了线。 INDRA-FAZIA 系列论文（PRC 111, 044601；PLB 868, 139815）约束的是对称能 S(ρ)，不是不可压缩系数 K₀。K₀ 描述对称核物质（δ = 0）的曲率，其现代约束主要来自巨单极共振、集体流与核阻止本领。
  1. RAON 被误列为"主力平台"。 事实是：SCL3 已建成并调试，最高只到 18 MeV/u（Ar 束）；SCL2 设计目标 200 MeV/u（U⁷⁹⁺）仍在研制中。RAON 的第一篇物理成果（2026 年 6 月 8 日，*Results in Physics*）是 ⁴⁰Ar + p 在库仑位垒附近的低能弹性散射（4.4 / 5.9 / 8.3 MeV/u），不是费米能区重离子碰撞，也没有做对称能约束。RAON 目前不能与 INDRA-FAZIA 并列。
  1. SπRIT 的探测方向被忽略。 SπRIT 测 π⁻/π⁺ 谱比，探的是约 2ρ₀ 的超饱和密度；INDRA-FAZIA 探的是ρ/ρ₀ ≈ 1.01，即饱和密度附近略偏下。两者方向相反、互补而非重复，不能写成"三台机器做同一件事"。
  1. INDRA-FAZIA 的敏感密度被夸大。 PLB 868, 139815 给出的结果标注在 ρ/ρ₀ = 1.01，即只比饱和密度低 1%。说成"深亚饱和（如 0.5ρ₀）的对称能"是过度解读。
  1. "ML 碰撞参数重建"与"模型无关重建"是两件不同的事。
  • Frankland 等，*Phys. Rev. C* 104, 034609 (2021)：概率反演方法（条件概率 + gamma 分布涨落核 + 贝叶斯定理），不含机器学习，已用于真实 INDRA-FAZIA 数据，是 PRC 111, 044601 的地基。
  • Valente 与 Lopez（LPC Caen）：神经网络在 HIPSE / ELIE 模拟上训练，达到亚飞米精度，配合贝叶斯推断估计介质中核子-核子截面与最大密度。见 *Il Nuovo Cimento C* 48(2), 1–7 (2025) 及 EuNPC 2025 海报。

把后者说成"INDRA-FAZIA 用的方法"是错的。

  1. R_i 的谱系溯源需要分两段写。 第一段是概念先声：**F. Rami 等（FOPI 合作组，70 位作者），*Phys. Rev. Lett.* 84(6), 1120–1123 (2000)，用 ⁹⁶Ru+⁹⁶Zr 在 400 A·MeV 下做"同位旋示踪"，用的是质子数与 t/³He 产额比——能量高于费米能区，且尚未给出归一化的输运比公式。第二段才是本体：M. B. Tsang 等，*Phys. Rev. Lett.* 92, 062701 (2004)，用两个对称体系作参考、把归一化输运比系统化地用于对称能密度依赖。把 R_i 直接说成"Rami 2000 提出"是不准确的，把它的能量说成费米能区也是错的。**
  1. SπRIT 装置论文的期刊被记错。 完整装置论文是 **J. Barney 等，*Review of Scientific Instruments* 92, 063302 (2021)**，DOI 10.1063/5.0041191。常被误记为 *Nuclear Instruments and Methods in Physics Research A*。arXiv 检索页面里"Nuclear Inst. and Methods A 986 (2021) 164713"属于另一篇不相关的论文，是列表串行造成的误读。
  1. Jhang 等 PLB 813, 136016 (2021) 的结论方向被反转。 这篇的结论不是给出了对称能约束，而是七个输运模型分散过大、无法可靠提取约束。把 L 的数值（42 < L < 117 MeV）归到这一篇是错的——那是 **Estee 等，*Phys. Rev. Lett.* 126, 162701 (2021)** 的结果。
  1. PLB 868, 139815 的作者数与机构。 共 34 位作者；S. Mallik 的署名机构是加尔各答 Variable Energy Cyclotron Centre 与孟买 Homi Bhabha National Institute（不是 GANIL）；F. Gulminelli 与 D. Gruyer 在 LPC Caen；C. Ciampi 与 J. D. Frankland 在 GANIL。这四处常被混淆。
  1. E818 不是 2019 年，也不在 E789 之后立刻做。 E818 执行时间是 2022 年 5 月 7 日至 5 月 17 日，体系为 ⁵⁸Ni + ⁵⁸Ni 与 ³⁶Ar + ⁵⁸Ni，能量 74 MeV/nucleon。发言人 R. Bougault、G. Casini；数据负责人 J. D. Frankland、D. Gruyer。此次实验前 INDRA 更换了全新电子学与电源（诺曼底大区专项拨款，耗时两年）。
  1. FAZIA 的角度与望远镜数。 FAZIA wall 配置覆盖实验室系 1.7° 到 13.7°，共 192 个望远镜；INDRA 覆盖 14° 到 176°。耦合总接收度约 80%（INDRA 单独 90%，FAZIA 单独约 70%）。把耦合后的接收度写成 90% 是错的。
  1. INDRA-FAZIA 的名称与位置。 装置在法国 GANIL 的 D5 实验厅；FAZIA 是 Four Pi A and Z Identification Array 的缩写，2015 年投入运行；INDRA 是 1990 年代设计，564 个探测单元、17 个环、覆盖 4π 的约 90%、同位素分辨只到 A ≤ 8。两者的分辨能力差别是耦合的全部理由。
  1. "²³²Sn"式的笔误要避免。 SπRIT 的反应体系是 ¹³²Sn+¹²⁴Sn、¹¹²Sn+¹²⁴Sn、¹⁰⁸Sn+¹¹²Sn，入射能 E/A = 270 MeV。SπRIT 可以当活性靶用，漂移长度 50.5 cm，pad 平面 86.4 cm × 134.4 cm，12096 个 pad。
  1. K₀ 与"核物质的弹簧常数"这个比喻要用对。 K₀ 的定义是 K₀ = 9ρ₀²·d²(E/A)/dρ² 在 ρ = ρ₀ 处取值，描述的是对称核物质的曲率。不要把它说成"对称能的硬度"。
  1. 两篇 2026 年的中国团队工作容易混淆。 华南理工 Feng Zhaoqing 组的 LQMD 结果（*J. Phys. G*，2026，DOI 10.1088/1361-6471/ae9f77）给 L = 42 ± 25 MeV；中科院近物所 Gao-Chan Yong 的 impIBUU 结果（*Phys. Rev. C* 104, 014613 (2021)）给 L = 66.75 ± 24.75 MeV。同源数据、不同模型、结果相差约 25 MeV——这正是当前模型依赖的体量，不能合并成一句话。
  1. 综述的正式发表信息。 Mallik/Gulminelli/Ciampi/Gruyer 的综述不只是 arXiv 预印本，已正式发表：***Universe* 2026, 12, 202**。引用时应给正式卷期。arXiv 编号为 2605.06138，HAL 编号 hal-05627602。
  1. 李昂组论文的期刊名要写对。 是 *Scientia Sinica Physica, Mechanica & Astronomica*（《中国科学：物理学 力学 天文学》），不是 *Science China Physics, Mechanics & Astronomy*——这是两家不同期刊。DOI 为 10.1360/SSPMA-2025-0201，卷期为 56(7)，2026 年。
  1. NICER PSR J0437−4715 的两个半径不要混用。 Miller 等（*ApJL* 1000(2), L48, 2026）给 68% 区间 11.9–15.5 km（@ M = 1.418 ± 0.044 M☉）；Qi 等（arXiv:2605.05429）给 13.25 ⁺⁰·³⁴₋₀·₃₅ km（@ M = 1.38 ± 0.03 M☉）。两者建模假设、数据组合与区间口径不同，不是矛盾，但绝不能互换引用。另需注意 Miller 等论文的预印本摘要曾给 11.8–15.1 km，正式版为 11.9–15.5 km，应引正式版。
  1. "七个输运模型"与"两个理论解读"的关系要分层。 Jhang 等的七模型对照是合作组内部的系统比对（结论：无法提取约束）；Estee 等给出 42 < L < 117 MeV 是基于单一模型的提取；Yong 与 Feng 组的后续工作是各自独立模型的再解读。三层不能压成一句"重离子 π 数据给出 L = 42–117 MeV"。