碎片会说话

一、先说清楚这篇不写什么

知识库里已经有两篇讲碎片的文章。20260603 那篇讲的是散裂——轻弹打重靶，一头扎进去砸出一堆中子，落点在加速器驱动嬗变系统和核数据上。20260825 那篇讲的是碎裂跟散裂的区别，以及全世界稀有束流装置的格局：谁走飞行中分离、谁走在线同位素分离，为什么这几年新核素像下雨一样掉出来。

两篇讲的都是"碎片怎么被造出来、用什么机器装"。

这一篇换个方向：碎片被造出来之后，我们能从它身上读出什么。

同一批碎片，上一篇文章里它是产品，这一篇里它是探针。碎片的动量分布里写着核内核子原本怎么跑，碎片的同位素比例里写着对称能，碎片产额随中子数的指数斜率里写着两个反应系统之间那点同位旋差异。你甚至可以用碎片的平行动量分布反过来去量弹核的中子皮有多厚——这是最近几年才认真提出来的一条新路。

这条路上有四件事要分开讲：碎裂到底是怎么发生的（机制），丰中子碎片自己有什么脾气（碎片性质），两个反应的产额比为什么会服从一个漂亮的指数律、又什么时候会失灵（同位旋标度），以及那套动量分布探中子皮的办法现在到底走到哪一步。

二、碎裂不是把核砸碎，是把核削皮

先纠正一个画面。

很多人第一次听到"碎裂"，脑子里出现的是一颗核被撞成几十块。不对。在束流能量为每核子几百兆电子伏时，绝大多数的碎裂事件是擦着过去的。靶核很薄，一毫米上下，弹核以接近光速飞过，两个核只在边缘重叠一小块。重叠区域里的核子互相撞，被削走；剩下的大半个核几乎没被碰到。

这个过程分两步，名字都很好懂。

第一步叫磨蚀，英文 abrasion。想象你把一个苹果贴着边切掉一片。切的时候果肉是被硬生生削走的，不是被炸开的。这一步极快，大约十的负二十二次方秒量级，快到核子之间来不及互相通知。削走的核子带走什么动量，完全看它撞上那一刻自己在核里怎么运动——是随机的。

第二步叫消融，英文 ablation。被削过的核现在缺了一块，表面能上升，形状也歪了，它带着一份激发能继续往前飞。这份激发能会通过蒸发中子、质子、α粒子慢慢放掉，一路放到冷静下来。这一段时间尺度长得多。

整个反应的总时长不到一纳秒，碎片却以几乎和原来一样的速度冲出去。 这是碎裂反应最值钱的一条性质：碎片出靶就是束流，不需要停下来、不需要化学提取、不需要等待释放，半衰期短到微秒的核都能来得及被鉴定。全世界走飞行中分离路线的装置，吃的就是这碗饭。

但这里有一个物理学家惦记了几十年的量：削掉一个核子，到底往剩下的碎片里塞进多少激发能？ 常用的大致数目是每削掉一个核子沉积约二十七兆电子伏。这个数不是随便来的，它决定后面蒸发阶段能放出多少粒子，进而决定你最后看到的那片碎片的同位素比例长什么样。而它的理论根据有两条互相竞争的路子。

一条是雨贡纽（J. Hüfner）等人在 1975 年提出的液滴路线：核被削掉一块，表面积增大，多出来的表面能就变成激发能。另一条是格拉伯（R. J. Glauber）框架下的空穴路线：核子在费米海里有固定的占据，磨蚀瞬间挖走了几个核子，费米面以下留下空穴，这些空穴的能量之和就是激发能。两条路线给出的数不一样，而且在轻碎片区差别不小。这个分歧到今天也没有完全关掉。

再说一遍那个画面：不是砸碎，是削皮，然后带着伤口蒸发。

三、动量分布：一把一九七四年的老尺子

现在讲这一篇里最基础、也最好用的一件工具。

削走的核子带走什么动量是随机的，那么留下来的碎片，动量也就跟着变模糊。这个模糊有多大？一九七四年，石溪纽约州立大学的莫里斯·戈德哈伯（A. S. Goldhaber）在一篇只有三页的文章里给出了答案，标题是《碎裂过程的统计模型》，发表在 *Physics Letters B* 53, 306–308 (1974)。

他其实给了两个模型。第一个模型假设碎裂是个"快过程"：削掉核子那一刻，你只是在随机抽取核子里原本就有的费米运动。第二个模型假设它是个"慢过程"：核先被激发到温度 T、达到平衡，再慢慢蒸发。有意思的地方在于，两个完全不同的假设，最后给出同一个形状的公式：

σ² = (p_F² / 5) · A(A_p − A) / (A_p − 1)

这里 A_p 是弹核质量数，A 是剔除核子之后、还没开始蒸发的那个"预碎片"的质量数，p_F 是核内核子的费米动量。第一个模型里，p_F 就是字面意义上的费米动量；第二个模型里，这个位置换成温度。同一个公式，两种解释——这是这个公式最容易被二手材料讲丢的地方。

公式里那个 A(A_p − A) 特别值得盯一眼。它是抛物线形的：削掉很少时它小，削掉一半时最大，削掉大半时又开始变小。物理上也说得通——你削掉一半核子的时候，剩下那半个核的动量"账单"最不平衡。

这条抛物线，实验上真的看到了。

德国重离子研究中心（GSI）的 A. Bacquias 与 V. Föhr、D. Henzlova、A. Kelić-Heil、M. V. Ricciardi、K.-H. Schmidt（共六人）在 *Physical Review C* 85, 024904 (2012) 上做了一次系统测量：用每核子 1 吉电子伏的 ¹¹²Sn 束和 ¹²⁴Sn 束分别打 ¹¹²Sn 靶和 ¹²⁴Sn 靶，把碎片的质量数、电荷数和速度都逐一测定。结果就是那条抛物线——碎片的纵动量分布宽度先涨后落，峰顶落在质量数约为弹核一半的地方。Bacquias 那篇还把磨蚀、蒸发和"热崩裂"三段的贡献接进了一个解析公式，比只用磨蚀一段要准得多。

同时要记住这把尺子在哪里会失灵。

密歇根州立大学的 D. J. Morrissey 在 *Physical Review C* 39, 460–470 (1989) 给过一套纯经验的宽度参数化。它在碎片质量接近弹核时非常好用，是很多束流线设计的默认选择。但当碎片质量小于弹核的一半时，它开始严重偏离——Bacquias 那篇的图九把这一点画得很清楚。所以凡是看到有人拿 Morrissey 公式一路外推到轻碎片，可以直接存疑。

还有一处小地方：费米动量本身随核的质量变化并不剧烈，常用的参数化形式是 p_F(A) = 281 · (1 − A 的负零点五六八次方) 兆电子伏每光速，这个式子来自 Giacomelli 等人 2004 年的工作，基于 Moniz 等人 1971 年的电子散射数据。

为什么这把老尺子对本篇这么重要？因为宽度和费米动量挂钩，而中子和质子的费米动量是不一样的。 只要你能把碎片按中子数分开测，那条抛物线就不再是一条，而是一族——差的那一点，就是中子和质子分布差异的指纹。

四、丰中子碎片：它们怎么活下来的

现在换一个角度看同一件事：不看动量，看同位素。

一片碎片从磨蚀出来的时候，中子数和质子数大致还保留着弹核的比例。但接下来它要蒸发。蒸发先放谁？在丰中子核里，中子的分离能低，主要放中子。这意味着蒸发阶段会系统性地把碎片往质子多的方向推。

于是实验室里看到的碎片，已经是"热碎片"冷却之后的"冷碎片"，两者不是一回事。这个区别在论文里经常被一笔带过，但它决定了一件事：你从冷碎片反推出来的对称能系数，可能不是碎片形成那一刻的真实值。

德克萨斯 A&M 大学回旋加速器研究所的 D. V. Shetty、A. S. Botvina、S. J. Yennello、G. A. Souliotis 等人做过一项针对性的研究：比较 ⁵⁸Fe + ⁵⁸Ni、⁵⁸Fe + ⁵⁸Fe 与作为参照的 ⁵⁸Ni + ⁵⁸Ni 三个系统在每核子 30、40、47 兆电子伏下的碎片产额。结论很干脆：初级碎片分布会被后续衰变强烈改造，从碎片产额反推出的对称能明显低于正常核在饱和密度下的值，说明这些碎片是在降低了的密度下、并且高度激发的时候形成的。这项工作的预印本是 arXiv:nucl-ex/0409019。

那么丰中子碎片到底"丰"到什么程度？

过去两年，全世界三台机器分别交出了答卷，而且答案的口味完全不同。

日本理化学研究所的放射性同位素束流工厂（RIBF）2026 年在 *Progress of Theoretical and Experimental Physics*（PTEP）上报告了第一次用 ²⁰⁸Pb 初级束流做新核素搜索的结果。束流能量每核子 345 兆电子伏，打铍靶，用 BigRIPS 逐事例鉴别。他们一共找到二十二个新核素，方向要分开数：十七个丰中子的，包括 ¹⁷⁹Ho、¹⁸¹,¹⁸²Er、¹⁸⁴⁻¹⁸⁶Tm、¹⁸⁸,¹⁸⁹Yb、¹⁹¹Lu、¹⁹³,¹⁹⁴Hf、¹⁹⁵,¹⁹⁶Ta、¹⁹⁸,¹⁹⁹W、²⁰⁰,²⁰¹Re；还有五个丰质子的，¹¹⁸Ce、¹²⁰Pr、¹⁵²,¹⁵³Hf、¹⁵⁴Ta，十七加五正好二十二。作者是 N. Fukuda、S. Michimasa、H. Suzuki、Y. Shimizu、H. Takeda、Y. Togano、M. Yoshimoto、Y. Yanagisawa、H. Sakurai 等人，通讯由 S. Michimasa 执笔。

为什么要用铅而不用铀？ 这是这篇实验最聪明的一手。用 ²³⁸U 束虽然产额高，但铀自身会裂变，裂变产物会污染整个碎片谱。而 ²⁰⁸Pb 几乎不裂变，在同样的质量区里给出干净得多的谱。这就是为什么同样是在 N 等于 126 附近找核，铅束能一路推进到 ²⁰¹Re。

美国稀有同位素束流装置（FRIB）走的是另一条路。O. B. Tarasov、B. M. Sherrill、A. C. Dombos、K. Fukushima、A. Gade、K. Haak、M. Hausmann、D. Kahl、D. Kaloyanov、E. Kwan、H. K. Matthews、P. N. Ostroumov、M. Portillo、I. Richardson、M. K. Smith、S. Watters（全部来自 FRIB 与密歇根州立大学物理与天文系）在 *Physical Review C* 112, 034604 (2025) 上报告：用 ⁸²Se 束，每核子 228 兆电子伏，打碳靶，用先进的稀有同位素分离器（ARIS）鉴别，首次看到 ⁶³Sc、⁶⁵Ti、⁶⁶Ti、⁶⁸V 四个新同位素，另外有一个事例与 ⁶¹Ca 相符。这篇最值钱的地方不在发现本身，而在它提出的新系统学：作者用一套叫 ΔBE 的经验方法描述截面随中子数的走向，比传统的 Qg 方法在丰中子区更灵敏——因为截面本身携带核结合能的信息。用 Hartree-Fock-Bogoliubov 质量表 HFB-22 配合磨蚀-消融计算，能相当好地重现测到的截面。

中国这边，惠州强流重离子加速器装置（HIAF）在 2026 年 7 月 21 日开始试运行，不到一个月就交出了第一个物理结果。实验用增强器（BRing）提供的 ²⁰⁹Bi 束，能量是每核子 836.66 兆电子伏，打石墨靶，产物经高能放射性束流线（HIRIBL）筛选后注入高精度环形谱仪（SRing），在等时性光学模式下工作。团队累计观测到 十个 ¹⁵³Hf 离子，估算产生截面在皮巴量级，对应的回转时间是 1159.348 纳秒。结论是 ¹⁵³Hf 是一个束缚或弱束缚核素。这项工作由中国科学院近代物理研究所重离子科学与技术全国重点实验室牵头，联合中国科学院大学、德国 GSI 亥姆霍兹重离子研究中心、德国科隆大学、东华理工大学和先进能源科学与技术广东省实验室完成，发表在 *Science Bulletin* 上。几乎同时，日本 RIBF 也独立报告了 ¹⁵³Hf 的观测——两项独立测量互相印证。这从侧面说明，在核素图这个位置上，人类离边界有多近了。

五、同位旋标度：两个反应，一条指数

前面讲的是单个反应的碎片。现在要做一件更巧的事：把两个只差同位旋的反应放在一起，比它们的产额。

这件事的起点是二〇〇〇年。密歇根州立大学国家超导回旋加速器实验室（NSCL）的 H. S. Xu、M. B. Tsang、T. X. Liu、X. D. Liu、W. G. Lynch、W. P. Tan、C. K. Gelbke 等人（后面还有 L. Beaulieu、B. Davin、Y. Larochelle、T. Lefort、R. T. De Souza、R. Yanez、V. E. Viola、R. J. Charity、L. G. Sobotka 等一长串），在 *Physical Review Letters* 85(4), 716–719 (2000) 上报告了 ¹¹²Sn、¹²⁴Sn 四种组合在每核子 50 兆电子伏下的测量。他们把同位素、同中子素、同量异位素的产额比全都摆上桌，得出的结论很有画面感：气相的中子比液相多。

也就是说，一个被撞热、开始膨胀的核，会自发地把自己分成两部分——一部分是还结合成核的"液滴"，一部分是跑出来的自由核子组成的"气"。而中子更愿意待在气相里。这就是同位旋分馏。

一年之后，同一个团队的 M. B. Tsang、C. K. Gelbke、X. D. Liu、W. G. Lynch、W. P. Tan 等人在 *Physical Review C* 64, 054615 (2001) 把这套东西提炼成了一个公式。同一年他们还写了另一篇《核反应中同位旋标度的条件》，*Physical Review C* 64, 041603 (2001)。同一个团队、同一年、同一卷、只差页码——这两篇是这一行里最容易串号的一对，引的时候务必看页码。

公式本身干净得不像话。设反应 1 和反应 2 只在整体同位旋上不同，某个碎片（中子数 N、质子数 Z）在这两个反应里的产额比记为 R₂₁，那么

R₂₁(N, Z) = C · exp(αN + βZ)

C 是归一化常数，α 和 β 叫同位旋标度参数。整个规律就叫同位旋标度，英文 isoscaling。

漂亮的不是这个指数形式本身，而是 α 和 β 的物理含义。在统计框架里，α 正比于对称能系数除以温度，再乘上两个源的同位旋不对称度之差。换句话说，你量一个斜率，就有机会反推对称能。

这条关系很快就被两个不同角度的工作撑住。德国 GSI 的 A. S. Botvina、俄罗斯赫洛平镭研究所的 O. V. Lozhkin 与 GSI 的 W. Trautmann 在 *Physical Review C* 65, 044610 (2002) 上，用轻离子诱发反应的数据做了统计解释（预印本 nucl-th/0112049）。日本东北大学的 A. Ono（小野章） 与美国同行 P. Danielewicz、W. A. Friedman、W. G. Lynch、M. B. Tsang 则在 *Physical Review C* 68, 051601 (2003)（预印本 nucl-th/0305038）中指出，在动力学核碰撞里，同位旋分馏和同位旋标度其实是同一件事的两面。同组一年后又写了续篇，*Physical Review C* 70, 041604 (2004)。

然后是实验侧的数字。

最近、也最干净的一组数据来自印度加尔各答可变能量回旋加速器中心（VECC）的 K500 超导回旋加速器。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 等十八位作者，在 *Physical Review C* 112(2), 024614 (2025) 上报告了 ¹⁴N 与 ²⁰Ne 束打 ¹¹²Sn、¹¹⁶Sn、¹²⁴Sn 在每核子 18 到 30 兆电子伏的结果：所有体系里同位旋标度都清楚地成立，α 和 β 都提了出来。他们用双同位素比法另测温度，再反算对称能系数，得到 C_sym 大约从 24 兆电子伏降到 17 兆电子伏，对应激发能每核子 2.1 到 2.8 兆电子伏。文章还专门讨论了每核子 1.8 兆电子伏那个点上的不一致——作者自己判断，那个能量下多重碎裂已经减弱甚至消失，方法不适用了。

这里必须插一句老实话：同位旋标度不是一台直接读出对称能的机器。 α 里同时装着对称能和温度，温度要用另一套方法（双同位素比）单独测，误差会叠上去；再乘以不对称度之差，每一步都是误差来源。把 α 说成"对称能读数"，是过度简化。

六、另一条算盘：等压产额比

说完同位旋标度，得说一套容易跟它混起来、但其实是两回事的东西：等压产额比，英文 isobaric yield ratio，简称 IYR。

差别在比较的对象。同位旋标度比的是两个不同反应；等压产额比比的是同一个反应内部不同的同量异位素。 一个是横向跨系统，一个是纵向跨同量异位素。中文里这两个概念都常被简化成"同位旋标度"，一定要留神。

IYR 这条路的方法基础是修正的费雪模型。德克萨斯 A&M 大学的 M. Huang、Z. Chen、S. Kowalski、Y. G. Ma（马余刚）、R. Wada、T. Keutgen、K. Hagel、M. Barbui、A. Bonasera、C. Bottosso、T. Materna、J. B. Natowitz、L. Qin、M. R. D. Rodrigues、P. K. Sahu、J. Wang 等人在 *Physical Review C* 81, 044620 (2010) 上分析了多重碎裂的碎片产额，把对称能系数与温度之比 a_sym/T 作为碎片质量的函数提了出来，发现它从质量数 9 时的约 5 涨到质量数 37 时的约 16。他们还用反对称化分子动力学配合 GEMINI 衰变程序做了对照——从初级同量异位素算出来的值是 4 到 5，几乎不随质量变。这个反差说明了一件要紧的事：从最终碎片观测到的对称能系数，可能跟碎片形成那一刻的真实值差别很大。

一年后，马春旺（河南师范大学）、王芳、马余刚、金晨（上海应用物理研究所）在 *Physical Review C* 83, 064620 (2011) 上更进一步，用 ⁴⁸Ca 和 ⁶⁴Ni 弹核碎裂的碎片，把 IYR 同时跟对称能项、库仑能项、对能项和中子质子化学势差拉上了线性关系（预印本 arXiv:1107.0131）。

这里有一对兄弟文献，混淆率极高，请盯紧页码：

PRC 81, 044620 (2010) 是 Huang 等人做的，主导单位是德克萨斯 A&M 大学回旋加速器研究所；PRC 83, 064620 (2011) 才是马春旺等人做的，主导单位是河南师范大学。两篇都用修正费雪模型、都研究碎片产额比、都挂在 *Physical Review C* 上，只是差一个字的页码——044620 对 064620。并排放在一起看，差别一目了然；分开看，串号是迟早的事。

还有一篇同卷的也容易撞上：Sanjeev Kumar、Y. G. Ma、G. Q. Zhang、C. L. Zhou 的《用亚饱和密度下的多重碎裂探测对称能的密度依赖》，*Physical Review C* 84, 044620 (2011)。它和上面那篇 2011 年的论文同卷、页码只差六万——一个 83 卷 064620，一个 84 卷 044620。同一卷里 044620 和 064620 同时存在，这是这一行里最阴的一个坑。

七、同位旋标度什么时候会失灵——失灵本身也带信息

到这一步要问一个更细致的问题：那条漂亮的指数律，是不是到处都是？

不是。而且在什么地方不成立，比在什么地方成立更有意思。

德克萨斯 A&M 的 G. A. Souliotis 等人用 ²³⁴U、²³⁹U 和 ²⁴⁸Cm、²⁴⁴Cm 的自发裂变数据做过一项分析（预印本 arXiv:nucl-ex/0306009），发现同位旋标度在几个中子数位置明显被破坏：N 大约 62 附近、N 大约 80 附近、可能还有 N 大约 86 到 88。

这三个位置不是随便出现的。N 等于 82 是一个已知的中子壳闭合。N 等于 62 到 64 这一带，是通常认为有形变壳效应的地方。换言之，同位旋标度在哪里失灵，就等于在哪里标出了壳结构的位置。因为在这个框架里，碎片最终产物的基态结合能（也就包含基态壳修正）本来是被消掉的，留下的只是整个激发系统的性质——一旦有强壳效应插进来，那个"消掉"的操作就不再成立。

有意思的是，同在 N 等于 50 那个更硬的球壳附近，同位旋标度没有表现出超出统计涨落的不规则。同一篇工作给出的解释是：那个位置上的不规则可能来自互补碎片那一侧的壳结构，而不是这一侧。这种"从碎片 A 的行为推断碎片 B 的壳"，听起来绕，但正是裂变两分性带来的机会。

所以同位旋标度这把尺子有两面。它成立的那部分告诉你碰撞区域的温度和对称能；它失灵的那部分告诉你壳结构在哪里。 一份数据，两条收获。

八、用碎片动量分布量中子皮

现在讲这一篇里最新的一条路，也是用户最关心的那一块。

先把问题摆出来。中子皮的厚度，也就是中子分布半径与质子分布半径之差，是核物理里最想要、又最难量的数之一。它同时牵着核结构、对称能的密度依赖、和中子星的半径。目前唯一被公认"模型无关"的测量手段，是宇称破坏电子散射——但那需要极高精度的极化电子束和极长的运行时间，全世界能做的机器一只手数得过来。

所以，能不能在碎裂反应里找到便宜的替代方案？

河南师范大学物理学院的马春旺团队给出了一个具体方案。论文《弹核碎裂反应中碎片平行动量分布作为中子皮厚度的一个可能探针》，发表在 *Nuclear Science and Techniques*（英文期刊，简称 NST）35, 文章编号 99 (2024)，DOI 是 10.1007/s41365-024-01455-9。作者是 Chun-Wang Ma、Yi-Jie Duan、Ya-Fei Guo、Chun-Yuan Qiao、Yu-Ting Wang、Jie Pu、Kai-Xuan Cheng、Hui-Ling Wei，署名单位为河南师范大学物理学院（新乡）与河南省科学院核科学与技术研究所（郑州）。

思路是这样的：在兰州量子分子动力学（LQMD）模型里模拟 每核子 140 兆电子伏的 ⁴⁸Ca 打 ⁹Be 的弹核碎裂。初始中子密度用费米型分布给，通过改变这个分布来改变 ⁴⁸Ca 的中子皮厚度。然后去看剩余碎片的平行动量分布长什么样。

关键发现是：这个分布是不对称的。 如果用左右两侧宽度不同的组合高斯函数去拟合，就得到两个宽度参数，左边叫 Γ_L，右边叫 Γ_R。而在边缘碰撞里，Γ_L 与 ⁴⁸Ca 中子皮的厚度呈现出敏感的关联。文章拿丰中子硫同位素做例子，展示了这条关联。由于平行动量分布在实验上是容易测的（只要你有速度分辨），作者提出把 Γ_L 当作测量弹核中子皮厚度的一个可能探针。

同一团队还有两条相关路线，值得一起看。

第二条：构型信息熵。 论文《用构型信息熵确定中子皮厚度》，*Nuclear Science and Techniques* 33(1), 60–66 (2022)，文章编号 6，DOI 10.1007/s41365-022-00997-0，预印本 arXiv:2201.01442。作者是 Chun-Wang Ma、Yi-Pu Liu、Hui-Ling Wei、Jie Pu、Kai-Xuan Cheng、Yu-Ting Wang。他们用修正的统计磨蚀-消融模型计算了 每核子 350 兆电子伏的 ⁴⁰⁻⁶⁰Ca 打 ⁹Be 的碎裂反应，从核密度、同位素分布、质量分布和电荷分布里各算出一个信息熵，再去看它跟弹核中子皮厚度的关系。结论是：用同位素、质量和电荷分布算出来的那几个信息熵，与中子皮厚度都呈线性关系。

第三条：类弹碎片的截面。 复旦大学的 丁梦琪（Meng-Qi Ding）、方德清、马余刚在 *Nuclear Science and Techniques* 35(12), 文章编号 211 (2024)，DOI 10.1007/s41365-024-01584-1 上发了一篇综述《中子皮及其在重离子碰撞中的效应》，系统梳理了这件事。里面提到的一条是：用改进的 ImQMD-L 模型算 每核子 200 兆电子伏的 ¹²⁴Sn 与 ¹³²Sn 打 ¹²⁴Sn，质量数大于 100 的类弹碎片截面与中子皮厚度、对称能斜率呈正相关；为了压掉模型带来的系统误差，还可以用两个系统之间这个截面之差来进一步提纯信号。同一篇综述还整理出另一条：同位旋标度参数 α 与 ⁵⁰Ca 的中子皮厚度呈线性负相关，而类弹碎片的平均中子质子比也与中子皮厚度呈负线性关系。

这篇综述里还列了一长串与中子皮相关的观测量：中子移除截面、中子质子产额比、氚与氦三的产额比、中子质子动量差、同量异位素标度参数、光子产生、镜像核电荷改变截面差，等等。中子皮越厚，质子移除截面越小、中子移除截面越大，因为中子分布更往外扩，更容易在边缘碰撞里被削掉。

但这里必须踩一脚刹车。

上面这四条——Γ_L、构型信息熵、类弹碎片截面、同位旋标度参数 α——全都是理论模拟提出来的探针，目前都还没有被实验真正验证过。 它们的共同前提是：你先得相信那个输运模型或者统计模型对中子皮如何影响碎裂产物的描述是对的。而这恰恰是需要检验的那一环。用模型算出来的关联去测一个量，再拿这个量去约束模型——逻辑上是绕回来的。

所以正确的表述是：这是一条有前途、但还没有落地的新路。 它值得投入，因为它便宜、能覆盖不稳定核，而宇称破坏电子散射既贵又只能打少数几个稳定靶。但它现在给出的不是"中子皮的测量值"，而是"一个可能的中子皮探针"。

判断这条路的进度，有一个很简单的指标：看有没有实验数据点落在 Γ_L 对 δnp 的那张图上。截至现在，没有。

九、中国这边：两台机器和一整套方法

中国在这个方向上的位置，可以从"方法"和"装置"两头看。

方法这一头，河南师范大学马春旺团队是最活跃的一支。除了上面提到的三篇，他们还做过一件对整个领域都有用的事：把机器学习引进来做碎片产额预测。Chun-Wang Ma、Xiao-Bao Wei、Xi-Xi Chen、Dan Peng、Yu-Ting Wang、Jie Pu、Kai-Xuan Cheng、Ya-Fei Guo、Hui-Ling Wei（署名单位是河南师范大学物理学院粒子与核物理研究所）在 *Chinese Physics C* 46(7), 118–128 (2022) 上，用贝叶斯神经网络训练了 6393 个碎片、53 个已测碎裂反应的数据，做出了一个直接模型和一个用 FRACS 参数化引导的物理引导模型，能覆盖从每核子 40 兆电子伏到 1 吉电子伏的入射能、从 ⁴⁰Ar 到 ²⁰⁸Pb 的弹核、以及多种靶核。最新的一步是 Dan Peng 与 Chun-Wang Ma 在 *Physical Review C* 113(2), 024605 (2026) 上用同类方法处理散裂碎片的同位旋标度分析——他们发现 α 与两个源的不对称度平方差 Δ[(Z_s/A_s)²] 之间有很强的线性相关，并据此在温度指定的条件下提取了电荷数 4 到 26 的碎片对应的对称能系数。文章自己也点出了局限：靶更重（比如氙）时方法开始吃力。

装置这一头，有两台机器在动。

RIBLL2 在 2025 年完成了"全线路开通"。 这是兰州重离子加速器冷却储存环（HIRFL-CSR）的第二条放射性束流线，此前只用到中间焦平面，束线长度约 26 米，主要做 p 壳和 sd 壳的核。新的实验平台建在最后一个焦平面 F4，把整条 55 米的束线都用上了。肖晓东（Xiao-Dong Xu）、郑勇、孙志宇、宋玉楠、孙保华、寺岛（S. Terashima）、王长江、郭革、李广顺、魏秀林、许俊尧、张继超、曹勇、高丙水、韩家兴、刘锦荣、卢辰贵、金树亚、翁（H. J. Ong）、齐浩天、秦云、孙亚洲、谷畑勇夫（I. Tanihata）等三十余位作者，联合中科院近代物理研究所、中国科学院大学和北京航空航天大学，在 *Science Bulletin* 70, 1026–1029 (2025) 上报告了这件事情，用的是 每核子 400 兆电子伏的 ⁴⁰Ar 束和每核子 350 兆电子伏的 ⁷⁸Kr 束做测试，验证了粒子鉴别能力。注意这是装置论文，不是物理结果论文——它的价值在于确认这条束线从此能干中重质量区的活。

HIAF 则刚刚交出了它的第一个物理结果，就是前面说过的 ¹⁵³Hf。它用的是弹核碎裂产生放射性核素、HIRIBL 分离传输、SRing 等时性质谱鉴别的完整链条。十亿分之一都无法描述的极低截面下，一共十个离子，这就是它交出的第一份答卷。近物所在报道里说，随着性能提升，未来实验亮度有望提高约两个数量级。

对碎裂物理来说，这两件事意味着什么？RIBLL2 的 55 米束线和 HIAF 的高流强，恰好是这条路上最缺的两样东西：更好的分离纯度和更高的统计量。所有前面提到的探针——Γ_L、信息熵、截面系统学——都卡在统计量上。中子皮信号是一点一点抠出来的，没有统计量就没有信号。

十、顺便说一个争议：Zipf 定律

讲碎裂，绕不开一个漂亮得让人不安的规律。

一九九九年，当时在中科院上海原子核研究所（现在的中科院上海应用物理研究所）的马余刚在 *Physical Review Letters* 83, 3617–3620 (1999) 上提出：在核的解体过程中，把一次事件里的碎片按大小排序，第 n 名的平均电荷数或质量数与 n 成反比。这正是语言学里的 Zipf 定律，讲的是文章里第 n 常见的词，出现频率差不多是 1/n。马余刚发现，这个规律恰好出现在核的液气相变温度上。这就给液气相变提供了一个全新的判据。同年他还在 *European Physical Journal A* 6, 367 (1999) 上做了有限尺度分析。

这个想法很有吸引力，也确实在实验数据里看到过。但有一篇反驳值得知道。

X. Campi 与 H. Krivine（预印本 arXiv:cond-mat/0508541）指出：Zipf 定律其实只是碎片大小分布为幂律（指数约等于 2）的一个数学推论。谁都知道这两个指数之间有关系：Zipf 指数 λ 等于 1 除以（τ 减 1）。而幂律碎片分布本身并不能可靠地指示液气相变——它也可以在别的地方出现。所以把 Zipf 定律当成相变信号，逻辑上是把结论当成了前提。

这个争议到今天没有彻底关掉。 提到 Zipf 定律的时候，把反驳一起讲，才是诚实的做法。

十一、一页总结

回到最初那个问题：碎片能告诉我们什么。

机制上，碎裂是"削皮加蒸发"两步，不是砸碎。削走核子时带走的动量是随机的，于是在剩下碎片的动量分布里留下一道抛物线形的模糊——戈德哈伯 1974 年的公式管到今天，但它在轻碎片区需要 Bacquias 那一类把蒸发和热崩裂接进去的改进，Morrissey 的经验公式则在碎片小于弹核一半时严重失效。

碎片性质上，你看到的碎片早已不是刚形成的样子。蒸发把丰中子碎片往质子多的方向推，所以从冷碎片反推的对称能会偏低。丰中子极端碎片的生产，现在由 RIBF 的铅束（二十二个新核素，其中十七个丰中子）、FRIB 的硒束（四个新核素加一个 ⁶¹Ca 候选）、和 HIAF 的铋束（十个 ¹⁵³Hf 离子）三头并进，而截面系统学本身就在提供结合能的信息。

同位旋标度上，两个只差同位旋的反应，产额比服从一条指数律，斜率里装着对称能和温度的组合。链子是 Xu 等 2000 年的同位旋分馏、Tsang 等 2001 年的公式、Botvina 等 2002 年和 Ono 等 2003 年的理论解释。最新实验值来自 VECC，C_sym 大约 24 到 17 兆电子伏。姊妹方法等压产额比走的是另一条算盘，两篇同卷兄弟文献请认准页码。

中子皮探针上，马春旺团队的平行动量分布方案（Γ_L）是最新的一条，构型信息熵、类弹碎片截面、同位旋标度参数 α 各成一路。四条都是理论模拟的产物，都还没有实验验证。

最后一句实话：碎裂反应是这个领域最老的工具之一，二十世纪七十年代就在用了。但"怎么读它"这件事一直在更新——从读产额，到读比值，到读熵，到读动量分布的非对称。这条路的尽头，是把一台用于生产的机器，变成一台用于测量的机器。

参考文献

本列表按正文讨论的先后顺序编号。正文各处论断均已随文写出作者、期刊、卷期页与 DOI，本列表作为汇总索引，供逐条复核。

A. S. Goldhaber，"Statistical models of fragmentation processes"，*Physics Letters B* 53, 306–308 (1974)。DOI: 10.1016/0370-2693(74)90388-8。作者单位：石溪纽约州立大学。
A. Bacquias, V. Föhr, D. Henzlova, A. Kelić-Heil, M. V. Ricciardi, K.-H. Schmidt（共 6 位作者），"Dispersion of longitudinal momentum distributions induced in fragmentation reactions"，*Physical Review C* 85, 024904 (2012)。DOI: 10.1103/PhysRevC.85.024904；预印本 arXiv:1111.1862。第一作者单位：GSI 亥姆霍兹重离子研究中心；合作单位含斯特拉斯堡大学、于韦斯屈莱大学。
D. J. Morrissey，"Systematics of momentum distributions from reactions with relativistic ions"，*Physical Review C* 39, 460–470 (1989)。DOI: 10.1103/PhysRevC.39.460。
J. Hüfner, K. Schäfer, B. Schürmann，"Abrasion-ablation in reactions between relativistic heavy ions"，*Physical Review C* 12, 1888–1898 (1975)。DOI: 10.1103/PhysRevC.12.1888。作者单位：海德堡马克斯·普朗克核物理研究所。
D. V. Shetty, A. S. Botvina, S. J. Yennello, G. A. Souliotis, E. Bell, A. Keksis，"Effect of neutron composition and excitation energy on the primary fragment yield distribution in multifragmentation reaction"，预印本 arXiv:nucl-ex/0409019。作者单位：德克萨斯 A&M 大学回旋加速器研究所；A. S. Botvina 另属俄罗斯科学院核研究所。
N. Fukuda, S. Michimasa, H. Suzuki, Y. Shimizu, H. Takeda, Y. Togano, M. Yoshimoto, Y. Yanagisawa, K. Kusaka, M. Ohtake, H. Sato, T. Sumikama, K. Yoshida, S. Motomura, H. Baba 等，"First observation of twenty-two exotic isotopes using a ²⁰⁸Pb primary beam at RIBF"，*Progress of Theoretical and Experimental Physics* 2026（PTEP 编辑推荐）。作者单位：日本理化学研究所仁科加速器科学研究中心。通讯执笔：S. Michimasa。
O. B. Tarasov, B. M. Sherrill, A. C. Dombos, K. Fukushima, A. Gade, K. Haak, M. Hausmann, D. Kahl, D. Kaloyanov, E. Kwan, H. K. Matthews, P. N. Ostroumov, M. Portillo, I. Richardson, M. K. Smith, S. Watters，"Discovery of new isotopes in the fragmentation of ⁸²Se and insights into their production"，*Physical Review C* 112, 034604 (2025)。DOI: 10.1103/573p-7fjp。作者单位：美国密歇根州立大学稀有同位素束流装置（FRIB）与物理与天文系。
中科院近代物理研究所，"Discovery of isotope hafnium-153 near the proton drip line via isochronous mass spectrometry: extending the nuclear landscape with HIAF"，*Science Bulletin* (2026)。DOI: 10.1016/j.scib.2026.06.056。牵头单位：中国科学院近代物理研究所重离子科学与技术全国重点实验室；合作单位：中国科学院大学、德国 GSI 亥姆霍兹重离子研究中心、德国科隆大学、东华理工大学、先进能源科学与技术广东省实验室。
H. S. Xu, M. B. Tsang, T. X. Liu, X. D. Liu, W. G. Lynch, W. P. Tan, A. Vander Molen, G. Verde, A. Wagner, H. F. Xi, C. K. Gelbke, L. Beaulieu, B. Davin, Y. Larochelle, T. Lefort, R. T. De Souza, R. Yanez, V. E. Viola, R. J. Charity, L. G. Sobotka，"Isospin fractionation in nuclear multifragmentation"，*Physical Review Letters* 85(4), 716–719 (2000)。DOI: 10.1103/PhysRevLett.85.716；预印本 nucl-ex/9910019。第一作者单位：密歇根州立大学国家超导回旋加速器实验室。
M. B. Tsang, C. K. Gelbke, X. D. Liu, W. G. Lynch, W. P. Tan 等，"Isoscaling in statistical models"，*Physical Review C* 64, 054615 (2001)。DOI: 10.1103/PhysRevC.64.054615。作者单位：密歇根州立大学。
M. B. Tsang, W. A. Friedman, C. K. Gelbke, W. G. Lynch, G. Verde 等，"Conditions for isoscaling in nuclear reactions"，*Physical Review C* 64, 041603 (2001)。DOI: 10.1103/PhysRevC.64.041603；预印本 nucl-ex/0106009。
A. S. Botvina, O. V. Lozhkin, W. Trautmann，"Isoscaling in light ion induced reactions and its statistical interpretation"，*Physical Review C* 65, 044610 (2002)。DOI: 10.1103/PhysRevC.65.044610；预印本 nucl-th/0112049。作者单位：德国 GSI 亥姆霍兹重离子研究中心；赫洛平镭研究所。
A. Ono, P. Danielewicz, W. A. Friedman, W. G. Lynch, M. B. Tsang，"Isospin fractionation and isoscaling in dynamical nuclear collisions"，*Physical Review C* 68, 051601 (2003)。DOI: 10.1103/PhysRevC.68.051601；预印本 nucl-th/0305038。第一作者单位：日本东北大学。
A. Ono, P. Danielewicz, W. A. Friedman, W. G. Lynch, M. B. Tsang，"Symmetry energy for fragmentation in dynamical nuclear collisions"，*Physical Review C* 70, 041604 (2004)。DOI: 10.1103/PhysRevC.70.041604；预印本 nucl-th/0409027。
P. Karmakar, S. Kundu, T. K. Rana, S. Mallik, S. Manna, R. Pandey, K. Banerjee, P. Roy, S. Mukhopadhyay, A. Sen, T. K. Ghosh, P. Pant, G. Mukherjee, A. Sultana, D. Mondal, D. Pandit, R. Shil, K. Rani，"Symmetry energy and isoscaling behavior of fragments emitted in ¹⁴N,²⁰Ne + ¹¹²,¹¹⁶,¹²⁴Sn at 18–30 MeV/nucleon"，*Physical Review C* 112(2), 024614 (2025)。DOI: 10.1103/xtx3-mh6c。作者单位：印度加尔各答可变能量回旋加速器中心（VECC）K500 超导回旋加速器。
M. Huang, Z. Chen, S. Kowalski, Y. G. Ma, R. Wada, T. Keutgen, K. Hagel, M. Barbui, A. Bonasera, C. Bottosso, T. Materna, J. B. Natowitz, L. Qin, M. R. D. Rodrigues, P. K. Sahu, J. Wang，"Isobaric yield ratios and the symmetry energy in heavy-ion reactions near the Fermi energy"，*Physical Review C* 81, 044620 (2010)。DOI: 10.1103/PhysRevC.81.044620。主导单位：德克萨斯 A&M 大学回旋加速器研究所；合作单位含中科院近代物理研究所、中科院上海应用物理研究所、西里西亚大学、鲁汶天主教大学、意大利南方国家实验室。
C.-W. Ma, F. Wang, Y.-G. Ma, C. Jin，"Isobaric yield ratios in heavy-ion reactions, and symmetry energy of neutron-rich nuclei at intermediate energies"，*Physical Review C* 83, 064620 (2011)。DOI: 10.1103/PhysRevC.83.064620；预印本 arXiv:1107.0131。作者单位：河南师范大学（新乡）、中科院上海应用物理研究所。
S. Kumar, Y. G. Ma, G. Q. Zhang, C. L. Zhou，"Probing the density dependence of the symmetry energy via multifragmentation at subsaturation densities"，*Physical Review C* 84, 044620 (2011)。DOI: 10.1103/PhysRevC.84.044620；预印本 arXiv:1111.4163。作者单位：中科院上海应用物理研究所。
G. A. Souliotis 等，"Isoscaling studies of fission: a sensitive probe into the dynamics of scission"，预印本 arXiv:nucl-ex/0306009。作者单位：德克萨斯 A&M 大学回旋加速器研究所。
G. A. Souliotis 等，"Heavy-residue isoscaling as a probe of the symmetry energy of hot fragments"，预印本 arXiv:nucl-ex/0511008。作者单位：德克萨斯 A&M 大学回旋加速器研究所，使用 MARS 反冲分离器与超导螺线管束线。
C.-W. Ma, Y.-J. Duan, Y.-F. Guo, C.-Y. Qiao, Y.-T. Wang, J. Pu, K.-X. Cheng, H.-L. Wei，"A possible probe to neutron-skin thickness by fragment parallel momentum distribution in projectile fragmentation reactions"，*Nuclear Science and Techniques* 35, 文章编号 99 (2024)。DOI: 10.1007/s41365-024-01455-9。作者单位：河南师范大学物理学院（新乡）、河南省科学院核科学与技术研究所（郑州）。
C.-W. Ma, Y.-P. Liu, H.-L. Wei, J. Pu, K.-X. Cheng, Y.-T. Wang，"Determination of neutron-skin thickness using configurational information entropy"，*Nuclear Science and Techniques* 33(1), 60–66 (2022)，文章编号 6。DOI: 10.1007/s41365-022-00997-0；预印本 arXiv:2201.01442。作者单位：河南师范大学物理学院、河南省光电传感集成应用重点实验室。
M.-Q. Ding, D.-Q. Fang, Y.-G. Ma，"Neutron skin and its effects in heavy-ion collisions"，*Nuclear Science and Techniques* 35(12), 文章编号 211 (2024)。DOI: 10.1007/s41365-024-01584-1。作者单位：复旦大学。通讯作者：方德清、马余刚。
C.-W. Ma, X.-B. Wei, X.-X. Chen, D. Peng, Y.-T. Wang, J. Pu, K.-X. Cheng, Y.-F. Guo, H.-L. Wei，"Precise machine learning models for fragment production in projectile fragmentation reactions using Bayesian neural networks"，*Chinese Physics C* 46(7), 118–128 (2022)。DOI: 10.1088/1674-1137/ac5efb。作者单位：河南师范大学物理学院粒子与核物理研究所。
D. Peng, C.-W. Ma，"Probing nuclear symmetry energy via isoscaling analysis of spallation fragments with a Bayesian neural network"，*Physical Review C* 113(2), 024605 (2026)。DOI: 10.1103/37mt-m922。2026 年 2 月 5 日发表。
X.-D. Xu, Y. Zheng, Z.-Y. Sun, Y.-N. Song, B.-H. Sun, S. Terashima, C.-J. Wang, G. Guo, G.-S. Li, X.-L. Wei, J.-Y. Xu, J.-C. Zhang, Y. Cao, B.-S. Gao, J.-X. Han, J.-R. Liu, C.-G. Lu, S.-Y. Jin, H. J. Ong, H.-T. Qi, Y. Qin, Y.-Z. Sun, I. Tanihata, L.-P. Wan, K.-L. Wang, S.-T. Wang, X.-X. Wang, T.-Y. Wu, X.-T. Wu, M.-X. Zhang, W.-W. Zhang, X.-B. Zhang, X.-H. Zhang, J.-W. Zhao, Z.-C. Zhou, L.-H. Zhu，"Full realization of the RIBLL2 separator at the HIRFL-CSR facility"，*Science Bulletin* 70, 1026–1029 (2025)。DOI: 10.1016/j.scib.2025.01.020；预印本 arXiv:2505.00053。作者单位：中科院近代物理研究所、中国科学院大学、北京航空航天大学。
Y. G. Ma，"Application of information theory in nuclear liquid gas phase transition"，*Physical Review Letters* 83, 3617–3620 (1999)。DOI: 10.1103/PhysRevLett.83.3617；预印本 nucl-th/0102019。作者单位：CCAST 世界实验室（北京）、中科院上海原子核研究所。
X. Campi, H. Krivine，"Zipf's law in multifragmentation"，预印本 arXiv:cond-mat/0508541。作者单位：巴黎第十一大学（CNRS UMR 8626）。
C.-W. Ma, H.-L. Wei, X.-Q. Liu, J. Su, H. Zheng, W.-P. Lin, Y.-X. Zhang，"Nuclear fragments in projectile fragmentation reactions"，*Progress in Particle and Nuclear Physics* 121, 103911 (2021)。DOI: 10.1016/j.ppnp.2021.103911。署名单位含四川大学、河南师范大学、中国原子能科学研究院、陕西师范大学。
T. Kubo 等（RIKEN RIBF New Isotope Collaboration），"Neutron dripline search for fluorine, neon and sodium and the discovery of ³⁹Na at RIKEN RIBF"，*Nuclear Physics A* 1062, 123154 (2025)。DOI: 10.1016/j.nuclphysa.2025.123154。作者单位：日本理化学研究所仁科加速器科学研究中心。

核实纠错表

以下每一条都是本轮核查中发现的、在中文与英文二手材料里反复出现或极易出现的问题。

#	常见说法	实际情况	依据
1	戈德哈伯 1974 年那篇在 PRL 上	在 Physics Letters B 53, 306–308 (1974)，DOI 10.1016/0370-2693(74)90388-8	期刊官网与 INSPIRE 一致
2	戈德哈伯只提出"费米运动"一个模型	他提了两个：快过程（费米运动）与慢过程（热平衡）。两者给出同一个宽度公式，只是参数含义不同	原文标题为 "Statistical models"（复数）
3	同位旋标度就是"同位旋相关的产额规律"	它特指两个只差同位旋的反应之间的产额比服从指数律；等压产额比（IYR）是同一反应内同量异位素之比，两者不可混称	Tsang 等 PRC 64, 054615 (2001) 与 Ma 等 PRC 83, 064620 (2011) 定义不同
4	Tsang 2001 年只有一篇同位旋标度论文	同年同团队有两篇同卷文章：PRC 64, 054615（"Isoscaling in statistical models"）与 PRC 64, 041603（"Conditions for isoscaling"）	页码不同，务必对页
5	IYR 那篇 PRC 81, 044620 是马春旺组做的	那是 Huang 等（德克萨斯 A&M 主导）；*马春旺等的才是 PRC* 83, 064620**	两篇作者与单位完全不同
6	PRC 84, 044620 与前一条是同一篇	不是。PRC 84, 044620 是 Kumar、Ma、Zhang、Zhou 的多重碎裂论文，与 PRC 83, 064620 同卷不同篇	卷号 83 与 84 只差一位，页码 064620 与 044620 也极像
7	马春旺的单位就是河南师范大学	早期（2021 年 PPNP 综述）署名含四川大学、中国原子能科学研究院、陕西师范大学；2022 年后的工作出现河南省科学院核科学与技术研究所（郑州）与河南省光电传感集成应用重点实验室	各篇署名页
8	NST 33 (2022) 那篇的卷期页只有一种写法	数据库里同时存在 "Vol. 33, no. 1, pp. 60–66" 与 "Volume 33, article number 6" 两种记法。NST 官方用文章编号	Springer 官网目次
9	Tarasov 2025 那篇 DOI 可写成 10.1103/PhysRevC.112.034604	2025 年后 APS 启用随机式 DOI，真值是 10.1103/573p-7fjp；按老式拼装会返回 404	APS 官网
10	FRIB 用 ⁸²Se 打了五个新同位素（含 ⁶⁷Cr 或 Cr-71）	论文摘要与 FRIB 官方新核素列表都只认四个（⁶³Sc、⁶⁵,⁶⁶Ti、⁶⁸V），另有一个事例与 ⁶¹Ca 相符。机构新闻稿里出现的"第五个"未进论文	APS 摘要在不同平台还有 1.07 与 1.83 pµA 两种束流强度，以 APS 官网的 1.07 pµA 为准
11	RIBF 那次是"发现了二十二个新同位素"	要分成十七个丰中子 + 五个丰质子，方向不能省。丰质子的五个是 ¹¹⁸Ce、¹²⁰Pr、¹⁵²,¹⁵³Hf、¹⁵⁴Ta	PTEP 2026 原文
12	²⁰⁸Pb 束与 ²³⁸U 束的碎裂实验可以互换着讲	用铅束的关键优势是铅几乎不裂变，没有裂变产物污染，因此能给出更干净的谱。RIBF 的稀土新核素（Cs-152 等七个）用的是铀束的飞行中裂变，两条线不同	JPS Hot Topics 与 PTEP 原文
13	HIAF 的首个物理成果是 2025 年	2025 年 10 月是完成束流调试；HIAF 2026 年 7 月 21 日开始试运行，首个物理结果（¹⁵³Hf）2026 年 7 月底至 8 月初公布	近物所中英文报道
14	RIBLL2 那篇是物理结果论文	它是装置论文，用 ⁴⁰Ar（400 MeV/u）与 ⁷⁸Kr（350 MeV/u）验证束线与探测器性能	Sci. Bull. 70, 1026–1029 (2025) 摘要
15	Morrissey 的经验宽度公式可以外推到轻碎片	它在碎片质量小于弹核一半时严重失效，不能外推	Bacquias 等 PRC 85, 024904 (2012) 图 9
16	同位旋标度参数 α 可以直接读出对称能	α 里同时装着对称能与温度，温度须用双同位素比法另测，误差叠加；且正比于两源不对称度之差	Tsang 等 (2001)、Karmakar 等 (2025)
17	Zipf 定律是液气相变的可靠证据	Campi 与 Krivine 指出 λ = 1/(τ−1)，它只是幂律碎片分布的数学推论，而幂律本身不必然指示相变	arXiv:cond-mat/0508541
18	用碎片动量分布测中子皮已经实现	现有四条探针（Γ_L、构型信息熵、类弹碎片截面、α）全部是理论模拟提出，尚无实验验证	来源均为模型计算论文
19	复旦那篇中子皮综述与马春旺组 2024 那篇是同一篇	两篇*同为 NST* 第 35 卷：马春旺组的文章编号 99，复旦丁梦琪、方德清、马余刚的综述文章编号 211**	Springer 官网目次
20	印度那篇 isoscaling 论文的 S. Mallik 与 INDRA-FAZIA 的 S. Mallik 是两个人	是同一人，署名单位为加尔各答 VECC 与孟买 Homi Bhabha National Institute	两篇的署名页
21	碎片产额直接等于碎片刚形成时的分布	蒸发阶段会系统性地把丰中子碎片推向质子多的方向，从冷碎片反推的对称能会低于真实值	Shetty 等 arXiv:nucl-ex/0409019；Huang 等 PRC 81, 044620 (2010)
22	碎裂就是"把核撞成很多块"	束流能量数百兆电子伏时，绝大多数事件是边缘擦碰，分"磨蚀"（削走）与"消融"（蒸发）两步，不是整核炸开	Goldhaber (1974)、Hüfner 等 (1975)