← 返回 PaperDaily 大模型与智能体

北航联合中科院新发现：铅靶才是测原子核质子半径的最优解

原子核就像一个挤满质子和中子的“小球宿舍”，核物理学家最关心的问题之一就是：这间宿舍到底有多大、质子住在哪一层。

龙哥读论文发布于 2026-08-20 00:20:17 阅读 2 查看原文

#AI科研 #人工智能

原论文信息如下：

论文标题:
Best Reaction Target To Determine Proton Distribution Radii of Atomic Nuclei 发表日期:
2026年08月发表单位:
School of Physics, Beihang University; Institute of Modern Physics, Chinese Academy of Sciences; Shanghai Institute of Applied Physics, Chinese Academy of Sciences; Sino-French Institute of Nuclear Engineering and Technology, Sun Yat-sen University; Medaustron; Key Laboratory of Radiation Physics and Technology of the Ministry of Education, Institute of Nuclear Science and Technology, Sichuan University; School of Nuclear Science and Technology, University of Chinese Academy of Sciences 原文链接:
https://arxiv.org/pdf/2608.14957v1.pdf

原子核就像一个挤满质子和中子的“小球宿舍”，核物理学家最关心的问题之一就是：这间宿舍到底有多大、质子住在哪一层。物质半径早就有了成熟的测量套路，可质子半径却一直是个“老大难”，尤其是对于那些寿命只有几毫秒、甚至几百纳秒的不稳定核。光靠传统电磁探针根本扫不到它们头上。最近，北航联合中科院近代物理研究所等单位，用一组“换靶”实验给出了一个又简单又反常识的答案：与其费劲做修正，不如直接换一块更重的靶——铅靶一上场，Glauber模型的麻烦自己就消失了。

核物理中的关键谜题：如何精确测量不稳定核的质子半径？

要弄懂这篇论文在解决什么问题，得先分清原子核的三把“尺子”：物质半径（核子整体分布）、中子半径和质子半径。三者分开测量，才能看出中子皮有多厚、质子分布有没有“外凸”，这些信息直接关系到中子星内部的状态方程长什么样，也关系到我们对奇异核结构的理解。可以说，原子核半径的精确测量，是整个核物理大厦的一块基石。

物质半径的测量从1980年代就起步了。Tanihata等人用高能重离子碰撞，测出相互作用截面（σ_I），再用Glauber模型反推半径。这招对不稳定核特别友好：束流打过去，数一数有多少个核发生了反应，截面就出来了。不过，相互作用截面只对核子总数敏感，质子和中子的贡献混在一起，分不开。

质子半径就得换个思路。J. Meng等人在2002年提出，用电荷改变截面（charge-changing cross section，σ_cc）来关联质子分布半径R_p。原理很直接：入射核的质子数只要变了，就算一次“电荷改变事件”。σ_cc的大小跟入射核里质子的空间分布强相关，所以可以从σ_cc里反推R_p。这个方法不需要复杂的电子束装置，只要有飞行时间谱仪加电离室，就能在次级束线上直接测，连寿命只有几百纳秒的核都能碰一碰。

但问题来了。Glauber模型假设只有弹核里的质子直接参与碰撞，中子全程当“吃瓜群众”。实际测出来的σ_cc却普遍比模型算出来的高约10%。为什么？因为碰撞中弹核可能先被撞掉一两个中子，剩下来的预碎片还处于激发态，一激动就蒸发个质子或α粒子出来，电荷数还是变了。这类带电粒子蒸发（charged-particle evaporation，CPE）通道，原始Glauber模型根本没算进去。

于是大家只能“打补丁”：先拿质子半径已知的稳定核做刻度，拟合出一个经验标度因子，再把这个标度因子搬到不稳定核上。不同课题组用的修正方案五花八门，有的用CPE修正，有的引入分离能依赖的S因子，一直没收敛成一个统一框架。更要命的是，有研究发现，氢靶和碳靶给出的同一批富中子氮同位素质子半径竟然对不上。到底谁准？没人能打包票。

这篇论文想干的事，就是把“靶核种类”这个变量一脚踢进分析里：既然标度因子这么烦人，能不能找到一种靶，让标度因子直接等于1，让Glauber模型裸算出马就行？

39组新数据揭示标度因子的靶核依赖规律

要回答这个问题，光靠零散的旧数据远远不够。研究团队在兰州重离子研究装置（HIRFL）的RIBLL2次级束线上，用280 MeV/nucleon的¹⁸O初级束流轰击铍靶，产生了从⁹Be到¹⁸O的一批p壳核素，然后让这些次级束依次穿过氢、碳、天然银和天然铅四种反应靶，一次性测得了39个全新的σ_cc数据——注意，是同一套实验装置、同一批核素、同一能量标定下拿到的四靶对比数据，这在此前是从未有过的。实验布局和粒子鉴别如图1所示。

图1：（a）包含RIBLL2和外部靶装置（ETF）的实验设施示意图；（b）反应靶前的粒子鉴别谱；（c）标记入射离子¹⁶N后，反应靶后方的电荷数Z鉴别谱。红、蓝、灰直方图分别对应C、Pb和空靶数据。Pb和空靶的Z=7峰高已按C靶归一化，阴影区域为N_out选择窗口。

束流在靶中心附近的能量统一压低到约240 MeV/nucleon。测量的核素覆盖⁹¹⁰Be、¹⁰⁻¹³B、¹²⁻¹⁶C、¹⁴⁻¹⁷N和¹⁶⁻¹⁸O，共18个核素，每个核素在四种靶上各测一个σ_cc。氢靶的截面用聚乙烯靶数据减去碳靶数据得到。所有39个数据点都是首次报道。表1列出了完整测量结果。

表1：本工作在氢（σ_cc,H^exp）、碳（σ_cc,C^exp）、天然银（σ_cc,Ag^exp）和天然铅（σ_cc,Pb^exp）靶上测得的电荷改变截面。注1为同一装置已发表的数据（来自参考文献[27]）。

新数据到手后，作者分别用零程和有限程两种光学极限Glauber模型做了理论计算，两种计算都统一低估实验值6%到10%左右。如果把实验值和理论值的比值σ_cc^exp/σ_cc^cal拎出来，就能看到“标度因子”到底在怎么变化。

结果非常干净：标度因子不再是一个“神秘经验常数”，而是随着靶核和弹核的性质呈现单调的系统性变化。作者巧妙地构造了一个新物理量A_T^(1/3)×S，其中A_T是靶核质量数，S = S_n(A,Z) + Min(S_p(A−1,Z), S_α(A−1,Z))。这里S_n(A,Z)是弹核的单中子分离能，S_p(A−1,Z)和S_α(A−1,Z)分别是弹核去掉一个中子后的“预碎片”的单质子分离能和单α分离能。Min函数取两者中较小的那个，对应预碎片退激时择优选择质子发射还是α发射。

把标度因子对A_T^(1/3)×S画出来，有意思的画面出现了：碳、氮、氧同位素的数据从氢靶的约1.7，沿着一条平滑的曲线单调降至铅靶的约1.0，几乎线性收敛。铍和硼同位素也落在同一条趋势带上（见图2的小插图）。换言之，靶核越重，Glauber模型的直接计算就越接近实验值；到铅靶时，经验修正项基本可以扔进垃圾桶了。

图2：实验电荷改变截面与理论值之比σ_cc^exp/σ_cc^cal作为A_T^(1/3)×S的函数，S以MeV为单位。面板（a）中实心方块为本文最新数据，插图为Be和B同位素；（b）中开放符号为文献数据，按两个Z区间分组。不同靶以不同颜色区分。红、蓝、黑虚线代表弹核¹⁴C在不同靶上(σ_evap+σ_cc^cal)/σ_cc^cal的计算比值，对应的激发能分布分别由GS、AMD和IQMD模型给出。

高Z靶为何能免去经验修正？物理机制深度解析

看起来“靶越重越省事”，但重靶的电磁相互作用明明更强，为什么σ_cc反而更接近理论值？这就得把间接过程一笔一笔拆开算清楚。实验测到的σ_cc^exp，理论上可以分解成三部分贡献：直接带电粒子移除截面σ_cc^cal，纯中子移除后预碎片带电粒子蒸发产生的贡献σ_evap，以及电磁 dissociation（EMD）中质子发射的贡献σ_EMD^p。完整的表达式如图下方公式所示。

其中，σ_evap对应弹核与靶核碰撞时被“擦掉”i个中子形成预碎片，预碎片再蒸发带电粒子的贡献，P_i^evap就是对应的带电粒子蒸发概率，这个概率强烈依赖预碎片的激发能分布；N(E_γ)是靶核产生的虚光子谱，σ_γ^p是光核质子发射截面，两者的卷积对应EMD的质子发射贡献。

σ_cc^cal本身就是Glauber模型计算的那一项。后两项都是模型一开始没算进去的“额外收入”，也就是导致标度因子大于1的根源。那为什么铅靶上这两项变弱了呢？关键在于纯中子移除这个前序步骤。

纯中子移除主要发生在碰撞参数比较大的“擦边球”事件里。一个简单的几何模型就能说明问题：

其中b_max ≈ 1.2(A_P^(1/3) + A_T^(1/3))是最大碰撞参数，Δb是发生中子移除的反应重叠区宽度。靶核越大，b_max越大，而Δb近似不变，于是中子移除截面在总反应道中的占比σ_{-xn}/σ_cc就会系统性下降。换句话说，大靶把“擦边球”的相对面积摊薄了，留给CPE的“原料”自然就少了。

这个几何图像靠谱吗？论文用微观输运模型做了交叉验证。他们把IQMD（同位旋依赖的量子分子动力学）和GEMINI衰变模型耦合起来，模拟¹⁴C与碳靶、铅靶的碰撞，并把电荷改变道和CPE道的碰撞参数分布分别提取出来。结果见图3（a）：虽然铅靶的总电荷改变截面明显大于碳靶，但CPE道的碰撞参数分布宽度（FWHM）在两种靶上几乎一致，都是1.65 fm左右，直接印证了Δb近似不变的前提。代入公式2画出的中子移除份额随靶核大小的变化曲线，正好解释了为什么低Z靶的蒸发贡献更强。

图3：（a）IQMD+GEMINI计算中，¹⁴C与碳（红）和铅（蓝）碰撞时电荷改变截面随碰撞参数的变化，柱宽0.157 fm，空心和实心直方图分别表示总电荷改变截面和CPE贡献；（b）几何模型中中子移除份额的靶核大小依赖性。

那么EMD呢？毕竟铅靶的电荷高达82，库仑场强得吓人。这里有一个关键区分：EMD对相互作用截面σ_I影响巨大，尤其是中子晕核，E1光子一照，中子就飞了。但在σ_cc这个观测量上，EMD需要敲出的是质子，而不是中子。质子在核里通常被束缚得更深，240 MeV/nucleon能区的等效光子能量又不够高，加上重靶-重核之间的库仑位垒对带电粒子蒸发有压制作用，导致EMD对σ_cc的贡献微乎其微。论文用等效光子法估算，¹⁸O在铅靶上约300 MeV/nucleon时EMD对σ_cc的贡献低于1%；¹⁰⁷Ag在铅靶上所有非相干激发的总贡献也只有约0.2%。

这两股力量一推一拉，就把高Z靶变成了“免修正窗口”：一方面CPE的原料（中子移除事件）相对减少；另一方面EMD对σ_cc几乎不构成干扰。结果就是，铅靶上实验值与裸Glauber模型计算的比值几乎正好等于1，经验标度因子直接下岗。

从p壳到sd壳：普适性验证与未来展望

一个新规律如果只在p壳核素上成立，说服力终归有限。作者特意翻遍了文献，把sd壳核素（²⁰Ne、²⁸Si、³⁶⁴⁰Ar）的30个数据点，以及更多p壳核素的42个数据点全部拉进来，反应能量覆盖267到400 MeV/nucleon，靶核从氢一直排到铅。把这些独立数据集跟本工作的结果摆在一起，所有数据都落在同一条标度曲线上——没有例外。这说明，标度因子对A_T^(1/3)×S的依赖不是某个同位素链的巧合，而是跨壳层、跨能量、跨靶核体系的普适行为。

当然，论文也没有把话说满。⁹Be和¹⁰B两个核素被明确排除在标度关系之外，因为它们的单中子移除后伴随碎裂，分离能语言已经失效。此外，目前验证的核素都集中在p壳和sd壳的稳定区附近，更靠近滴线的极端同位旋核素是否仍然满足这条曲线，还是个开放问题。作者在总结中特别呼吁，未来应该对靠近滴线的核素开展更多σ_cc测量，把这个标度行为的适用范围彻底画清楚。

从实用角度看，这篇工作的直接贡献是给实验方案设计“指了一条明路”：以后想测不稳定核的质子半径，优先选铅靶，而不是传统的碳靶或氢靶。虽然铅靶上的事例率会低一些，但换来的好处是模型不确定性大幅下降，不再需要依赖那些纠缠不清的经验标度因子。对追求系统精度的核结构研究来说，这笔账相当划算。

龙迷三问

下面是龙哥对于大家可能的一些问题的解答：

这篇论文到底在解决什么问题？北京航空航天大学联合中科院近代物理研究所等机构利用兰州重离子加速器HIRFL，系统测量了18个p壳核素在氢、碳、银、铅靶上的39组电荷改变截面。

这篇工作最值得看的点是什么？实验测量了39个新的电荷改变截面数据，覆盖18个p壳核素在4种靶上的数据，与已有文献数据一致性好，并揭示了标度因子对靶核质量和分离能的普适依赖关系。

这篇工作的边界或风险在哪里？优点：实验数据量大且系统性强，首次在统一框架下比较多种靶核的电荷改变截面；揭示了标度因子对靶核质量和分离能的普适依赖关系；提出高Z靶可直接提取质子半径，简化了实验分析流程。缺点：仅覆盖p壳核素，对更重核素的普适性有待验证；EMD贡献的估计依赖于理论模型；对滴线附近核素的适用性尚不清楚。

如果你还有哪些想要了解的，欢迎在评论区留言或者讨论~

龙哥点评

论文创新性分数：★★★★☆

通过系统测量18个p壳核素在氢、碳、银、铅四种靶上的电荷改变截面，揭示标度因子对靶核质量和分离能的普适依赖关系，证明高Z靶（Pb）可直接提取质子半径而无需经验修正。

实验合理度：★★★☆☆

现有材料未完整覆盖数据划分、基线公平性和统计显著性，因此按中性评价处理。

学术研究价值：★★★★☆

通过系统测量18个p壳核素在氢、碳、银、铅四种靶上的电荷改变截面，揭示标度因子对靶核质量和分离能的普适依赖关系，证明高Z靶（Pb）可直接提取质子半径而无需经验修正；更关键的是问题定义是否可复用到同类任务。

稳定性：★★★☆☆

现有材料未提供充分的极端条件、重复运行或扰动测试，稳定性暂按中性评价。

适应性以及泛化能力：★★★☆☆

现有材料未完整展示跨数据集、跨场景或分布外实验，泛化能力仍需进一步验证。

硬件需求及成本：★★★☆☆

现有材料缺少完整训练资源、参数量、显存和推理时延信息，成本暂按中性评价。

复现难度：★★★☆☆

现有材料未确认完整代码、配置、数据处理脚本和权重是否齐备，复现难度暂按中性评价。

产品化成熟度：★★★☆☆

论文验证以研究实验为主，真实部署中的时延、成本、维护和异常场景仍需补充验证。

可能的问题：仅覆盖p壳核素，对更重核素的普适性有待验证；EMD贡献的估计依赖于理论模型；对滴线附近核素的适用性尚不清楚。

主要参考文献

[1] I. Tanihata et al., Phys. Rev. Lett. 55, 2676 (1985).

[2] J. Meng, S.-G. Zhou, and I. Tanihata, Phys. Lett. B 532, 209 (2002).

[3] T. Yamaguchi et al., Phys. Rev. Lett. 107, 032502 (2011).

[4] J. Zhang et al., Sci. Bull. 69, 1647 (2024).

[5] J. Zhao et al., Phys. Lett. B 858, 139082 (2024).

[6] M. Tanaka et al., Phys. Rev. C 106, 014617 (2022).

[7] W. R. Webber et al., Astrophys. J. 508, 949 (1998).

[8] 论文原文：https://arxiv.org/pdf/2608.14957v1.pdf

*本文仅代表个人理解及观点，不构成任何论文审核或者项目落地推荐意见，具体以相关组织评审结果为准。欢迎就论文内容交流探讨，理性发言哦～想了解更多原文细节的小伙伴，可以点击"阅读原文"，查看更多原论文细节哦！

质子在哪儿？铅靶告诉你！想第一时间跟进更多核物理、重离子物理、大模型+科研等前沿解读，欢迎加入龙哥读论文粉丝群，扫描下方二维码或者添加龙哥助手微信号加群：kangjinlonghelper。一定要备注：研究方向+地点+学校/公司+昵称（如核物理+兰州+中科院近物所+小核），根据格式备注，可更快被通过且邀请进群。『龙哥读论文』微信群目前包含：图像处理、大模型及智能体、自动驾驶及机器人、AI医疗及AI金融5个群，核物理群正在筹建中，快来占座！
wechat_helper