核图谱系统学
# 原子核的地图正在被重画
一张画了四千多个格子的地图，正在被一格一格地重画。
地图上的每个格子，是一种原子核。横着数质子，竖着数中子，从只有一个质子的氢，一路铺到一百多个质子的超重元素。这张地图叫核素图，是核物理最底层的坐标系。一百年前，人们开始往里填格子，以为把它填完只是时间问题。填到边界就填不动了——边界之外，多一个质子或多一个中子，核当场散架。这条边界，核物理叫它滴线。
地图上有几个特别显眼的位置，叫幻数：2、8、20、28、50、82、126。凑到这些数字的核格外结实，像地图上的堡垒。质子中子都凑满的，叫双幻核，是最结实的堡垒。比如铅-208，82个质子126个中子，稳到被全世界拿来当质量的基准。
堡垒是地图上最可靠的地标。造地图的人相信，幻数写在哪，堡垒就在哪。这个信念撑了半个多世纪。
直到一张照片把地图撕开一个口子。
``氧 Z=8 滴线: 24O (N=16) 25O 26O 27O 28O 全散架氟 Z=9 滴线: 31F (N=22) 多一个质子，边界跳了6个中子氖 Z=10 滴线: 34Ne (N=24)钠 Z=11 39Na 出现 (N=28) 28 这个幻数也不稳``
2023年，日本理化学研究所的加速器里，一个叫氧-28的原子核被造了出来，然后当场散架。这件事本身不稀奇，稀奇的是它本该是堡垒。氧-28有8个质子20个中子，20是幻数。在双幻核的候选名单里，它是最后一个还没被实验看到的。按壳模型的说法，它应该稳得像块石头。结果它是个非束缚的共振态，衰变能量只有0.46兆电子伏，裂成氧-24加四个中子，存在时间以仄秒计——十的负二十一次方秒。
带队的是东京工业大学的近藤洋介。论文根据实验测出的敲出截面认为，氧-28没有闭合壳——N=20这个幻数，在氧这里不顶用了。这个解读还有理论家在争论，但数据本身很硬。
这不是孤例。整条氧链都在说同一件事：中子滴线停在氧-24，那里的N是16，不是20。氧-25到氧-28全部非束缚。可只要多一个质子，变成氟，滴线一下子伸出去六个中子，到氟-31，N=22。差一个质子，边界跳了六个中子。这张地图的边缘不是平滑的曲线，是一道锯齿。
钠那边也出了怪事。2019年，理研的数据里闪过一个钠-39的痕迹，只有一个事件。2022年，正式确认，这次看到九个。钠-39有28个中子，比之前已知最丰中子的钠同位素多出两个。28这个幻数，在钠这里也摇摇欲坠。
这些不是零星的怪事。它们指向同一个东西：幻数会流动。
幻数不是刻在石头上的。它来自核子之间某种力的平衡，而这个平衡随质子中子的比例改变。东京大学的大塚孝治团队2005年指出，核子间的张量力会让自旋和轨道耦合的强度跟着质子数变化，幻数随之移动。三体力也插了一脚——氧的滴线能停在N=16，靠的是三体力的排斥性贡献，把N=20那层轨道的能量抬上去，让这个幻数不再牢靠。
钙-52的故事把这一点演得淋漓尽致。2013年，欧洲核子中心ISOLDE的潘宁阱测出钙-52的质量，发现N=32的壳隙大约4兆电子伏，和经典幻数一个量级。紧跟着，钙-54的激发态谱显示N=34也是幻数——大塚团队2001年预言，2013年理研实验证实。两个新幻数，都站住了。
然后，多一个质子。钙-54变成钪-55，N=34的壳隙崩到1.59兆电子伏，在所有N=34同中子素里最低。再添一个质子到钛-56，第一激发态只有1.127兆电子伏，幻数性荡然无存。一个质子，摧毁一个幻数。这不是比喻，是测量的结果。
幻数也会在某些方向变强。2024年，兰州近代物理研究所合成了锇-160和钨-156，发现N=82的壳闭合在缺中子一侧反而增强了，指向一个还没出生、但理论上该很稳的双幻核，铅-164。
物理学家给这种变动起了个名字，叫幻数演化。它意味着"魔法"不是核的固有属性，而是某个区域的局部性质。钙附近是魔法，钪和钛附近不是。这逼着理论家回答一个问题：到底是什么在控制壳层的升降？
答案是力。具体说，是张量力和三体力的组合。壳模型、从头算、密度泛函理论，三家给出的答案并不一致。有的相互作用能复现钙-54的闭合壳，却错过钪-55；有的正好相反。理论在同一个区域给出互斥的预言，这在核物理里相当罕见。
堡垒出了问题，那地图本身呢？
核素图的骨架是质量。每个格子离边界多远，都由质量决定——中子分离能是质量的差，滴线的位置由分离能决定，r过程核合成的路径由质量的高低决定。没有质量，整张地图就是一堆空框。
所以质量的精确测量，是核图谱系统学的地基。地基上站着几样工具。潘宁阱精度最高，最优情形下相对不确定度能到十的负十二次方，但只能测寿命够长的核。飞行时间质谱快，精度差一截。储存环等时质谱卡在中间，专治短寿命核。
中国在这块地基上，位置很特殊。
兰州近代物理研究所的CSRe储存环，发展出一种叫Bρ-定义等时质谱的技术，单离子测量精度能做到约150千电子伏。2023年，这套方法的核心论文发表在《自然·物理》上。用这套技术，兰州团队测了滴线附近的硅-23、磷-26、硫-27/28、氩-31，把质子滴线的位置钉得更牢，还发现氩-31可能是新的双质子晕核。同一条环，还在质子数等于中子数的核上测了一批质量——这些核直接检验三体力，测出来的数值，现有模型普遍对不上。
更特别的是，全球核素图的质量基准——AME原子质量评价——五人核心编辑团队里，有两位是中国学者：近代物理研究所的王猛，和黄文嘉。2021年发布的AME2020，三篇论文合计一千多页，发表在《中国物理C》上。王猛团队自己产数据，自己进评价表，自己再拿数据检验理论——产数据、评数据、用数据，一条闭环。世界上能玩这条完整链条的，数得出几个国家。
闭环的另一端是理论。质量模型的竞争是一场精确度的竞赛。目前领跑的是一个叫WS4的宏观-微观模型，广西师范大学王宁等人2014年发表，对已知质量的均方根偏差约298千电子伏，是所有全局模型里最准的。对手们各有短长：Skyrme型密度泛函的偏差从1.4到7兆电子伏不等；协变密度泛函DRHBc，用2021年后新测量的质量检验，偏差约1.35兆电子伏，比所有Skyrme都好，仍不如WS4；从头算从核力出发，不拟合全局质量，却最接近"第一性原理"的承诺。
差距有多重要？关键在滴线以外。
对已知质量的核，模型和模型之间偏差约1兆电子伏。越往滴线走，偏差越大，能大到20兆电子伏以上。不同的模型对滴线的位置给出不同的答案。一个叫RCHB的模型说，束缚核总共有9035个。另一个叫DRHBc的模型预言偶偶束缚核2583个，按比例折算总数大约一万个。两个模型连"地图上到底有多少格"都对不上，差出几千个。DRHBc还说滴线外面有"半岛"——一些孤立存在、周围全是散架的核。
质量不准，后果直接落在天体物理上。2017年，两颗中子星并合，引力波探测器和全世界几十台望远镜同时捕捉到。千新星的光谱里，丹麦哥本哈根大学的沃森团队鉴定出锶——人类第一次在一个并合事件里，直接看到r过程造出的元素。中子星并合从此成为r过程的一个确证来源。
要算出r过程造了多少金、多少银、多少铀，你得知道路径上一百多个核的质量。圣母大学的Mumpower团队2016年做过定量评估：现在最好的全局质量模型，误差约400到650千电子伏，太大，做不出可靠的丰度预测。需要压到100千电子伏以下，模型之间才能区分。WS4的298千电子伏，还差三倍。2022年，Zhu等人把七种质量模型塞进r过程计算，得到的丰度差异高达两到三个数量级。换句话说：核质量的不确定性，比天体物理环境本身的不确定性还大。
地图的另一端，是超重元素。
周期表的尽头，人们期待一个"稳定岛"——在Z=114、120或126附近，中子数凑到172或184时，超重核的寿命突然变长的区域。这个岛的准确位置，各家预言不一样。宏观-微观模型说中心在Z=114、N=184。协变密度泛函说在Z=120。Skyrme密度泛函也说Z=120。N=184相对稳固，质子幻数到底落在114、120还是126，悬而未决。2024年的综述《超重元素的追寻》把共识收拢了一下：那个岛存在，但它是"增强稳定"的岛——寿命从几秒延长到几分钟、几小时——不是真正稳定的超重核。
实验这边，进展慢得让人焦虑。2024年，伯克利用钛-50打钚-244，合成了两个𫟷-116原子——22天，两个原子。这是第一次用非"魔数"束流合成超重元素，验证了通往120号元素的新路线。但更早的尝试里，德国GSI的TASCA装置用钛-50打锫-249找119号元素、打锎-249找120号元素，连续照射四个月，什么都没看到，只给出排除上限：65飞靶和200飞靶。元素119和120，至今没有出生。
还有一个安静但深刻的战场：无中微子双β衰变。这是寻找中微子是否等于它自己的反粒子的关键实验。要从实验测到的半衰期反推中微子质量，必须知道一个中间量，叫核矩阵元。你可以把它想成一把扳机——同一个扳机，不同方法算出来的力度，差好几倍。北京大学的赵鹏巍团队2024年提出一个新框架，声称在最低阶就能严格算出这个量，不用依赖那些不确定的短程修正。如果对了，这是把扳机校准的一条新路。
回头说中国。除了CSRe和AME，中国还有一张深地王牌：锦屏地下实验室里的JUNA加速器。2400米深的岩层下，宇宙线被挡到几乎为零，这台400千伏、10毫安束流的加速器，是目前世界上最强的深地核天体物理装置。2022年，JUNA团队在《自然》上发表氟-19(p,γ)氖-20的测量，首次进入第一代恒星的伽莫夫窗口，在225千电子伏处发现一个新共振，反应率比旧推荐值高到7.4倍——这是中国核物理平台的第一篇《自然》。领衔作者是北京师范大学的张立勇，通讯作者包括何建军、美国圣母大学的韦舍尔、以及中国原子能科学研究院的柳卫平。锦屏之下，那颗号称"圣杯"的反应——碳-12(α,γ)氧-16，恒星演化绕不开的关键一步——也在被一步步逼近。
接下来十年，地图会怎么变？
美国的FRIB从2022年出束到现在，新同位素一个接一个被发现，铀束功率正往设计的400千瓦爬。日本的理研，过去几年发现的新核素占全球一半以上。德国的FAIR要上Super-FRS。中国的HIAF已经通过工艺验收，这是超导直线加速器加同步加速器加储存环的组合，直指119号和120号元素。四台机器，各自盯着一块地图边缘。
理论那边也有新动静。2025年，橡树岭国家实验室牵头的从头算工作预言，钙的中子滴线可能一直伸到中子数50附近——钙-70也许还束缚着。如果对，钙将是滴线反常外伸的第一个中重元素，N=50会在钙这里重新变成幻数。
但机器再强，也绕不开那个最根本的问题：当幻数本身会流动，我们以为的地图，到底还有多少是对的？
氧-28不是双幻核。钙-54是幻数核，多一个质子变钪-55，幻数又突然不是了。N=34在一个质子之间崩解。滴线的位置，模型和模型之间差着几千个核。稳定岛的中心，三家理论三个答案。核矩阵元，方法和方法之间差着几倍。这张地图的每一根线条，都在被测量和理论来回拉扯。
也许这正是这个领域最诚实的现状：核图谱系统学不是一张正在被誊抄的成品图，而是一张正在被重新测绘的草图。堡垒会移动，边界会跳变，幻数会生灭。而每一次移动，都逼着理论交出更精确的核力。关于核力的知识，就是这样一格一格被逼出来的。
*参考文献*
[1] Y. Kondo 等（东京工业大学/理研，SAMURAI21合作组），"First observation of 28O", *Nature* 620, 965-970 (2023). DOI:10.1038/s41586-023-06352-6. arXiv:2308.06433.
[2] D. S. Ahn 等（理研仁科中心），"Location of the Neutron Dripline at Fluorine and Neon", *Phys. Rev. Lett.* 123, 212501 (2019). DOI:10.1103/PhysRevLett.123.212501.
[3] D. S. Ahn 等（理研仁科中心），"Discovery of 39Na", *Phys. Rev. Lett.* 129, 212502 (2022). DOI:10.1103/PhysRevLett.129.212502.
[4] F. Wienholtz 等（欧洲核子中心ISOLTRAP），"Masses of exotic calcium isotopes pin down nuclear forces", *Nature* 498, 346-349 (2013).
[5] D. Steppenbeck 等（理研RIBF），"Evidence for a new nuclear 'magic number' from the level structure of 54Ca", *Nature* 502, 207-210 (2013). DOI:10.1038/nature12522.
[6] T. Otsuka, T. Suzuki, R. Fujimoto, H. Grawe, Y. Akaishi（东京大学/理研/GSI），"Evolution of Nuclear Shells due to the Tensor Force", *Phys. Rev. Lett.* 95, 232502 (2005).
[7] T. Otsuka, A. Gade, O. Sorlin, T. Suzuki, Y. Utsuno（东京大学/理研/密歇根州立大学等），"Evolution of shell structure in exotic nuclei", *Rev. Mod. Phys.* 92, 015002 (2020). arXiv:1805.06501.
[8] M. Wang, Y. H. Zhang, X. Zhou 等（中国科学院近代物理研究所），"Mass Measurement of Upper fp-Shell N=Z-2 and N=Z-1 Nuclei and the Importance of Three-Nucleon Force along the N=Z Line", *Phys. Rev. Lett.* 130, 192501 (2023).
[9] X. Zhou 等（中国科学院近代物理研究所），"Brho-defined isochronous mass spectrometry", *Nature Physics* 19, 1091-1097 (2023).
[10] Y. Yu 等（中国科学院近代物理研究所），"Nuclear Structure of Dripline Nuclei Elucidated through Precision Mass Measurements of 23Si, 26P, 27,28S, and 31Ar", *Phys. Rev. Lett.* 133, 222501 (2024). arXiv:2410.17701.
[11] W. J. Huang, M. Wang, F. G. Kondev, G. Audi, S. Naimi， "The AME2020 atomic mass evaluation (I)", *Chinese Physics C* 45, 030002 (2021). DOI:10.1088/1674-1137/abddb0.
[12] M. Wang, W. J. Huang, F. G. Kondev, G. Audi, S. Naimi，"The AME2020 atomic mass evaluation (II)", *Chinese Physics C* 45, 030003 (2021). DOI:10.1088/1674-1137/abddaf.
[13] F. G. Kondev, M. Wang, W. J. Huang, S. Naimi, G. Audi，"The NUBASE2020 evaluation of nuclear physics properties", *Chinese Physics C* 45, 030001 (2021). DOI:10.1088/1674-1137/abddae.
[14] N. Wang, M. Liu, X. Z. Wu, J. Meng（广西师范大学/中国原子能科学研究院），"Surface diffuseness correction in global mass formula", *Physics Letters B* 734, 215 (2014)（WS4）.
[15] B. S. Hu, A. Ekstrom, C. Forssen, G. Hagen, W. G. Jiang, T. Miyagi, T. Papenbrock（橡树岭国家实验室/查尔姆斯理工/田纳西大学），"The neutron dripline in calcium isotopes from a chiral interaction", *Phys. Rev. Lett.* (2025). arXiv:2512.11723.
[16] O. R. Smits, C. E. Duellmann, P. Indelicato, W. Nazarewicz, P. Schwerdtfeger（GSI/美因茨大学/梅西大学/索邦大学/密歇根州立大学），"The quest for superheavy elements and the limit of the periodic table", *Nature Reviews Physics* 6, 86-98 (2024). DOI:10.1038/s42254-023-00668-y.
[17] J. M. Gates 等（劳伦斯伯克利国家实验室），"Toward the Discovery of New Elements: Production of Livermorium (Z=116) with 50Ti", *Phys. Rev. Lett.* 133, 172502 (2024). arXiv:2407.16079.
[18] J. Khuyagbaatar 等（GSI TASCA），"Search for elements 119 and 120", *Phys. Rev. C* 102, 064602 (2020). DOI:10.1103/PhysRevC.102.064602.
[19] H. B. Yang 等（中国科学院近代物理研究所），"Discovery of New Isotopes 160Os and 156W: Revealing Enhanced Stability of the N=82 Shell Closure on the Neutron-Deficient Side", *Phys. Rev. Lett.* 132, 072502 (2024), Editors' Suggestion.
[20] M. R. Mumpower, R. Surman, G. C. McLaughlin, A. Aprahamian（圣母大学/洛斯阿拉莫斯/北卡州立），"The impact of individual nuclear properties on r-process nucleosynthesis", *Prog. Part. Nucl. Phys.* 86, 86-126 (2016). arXiv:1508.07352.
[21] D. Watson 等（哥本哈根大学），"Identification of strontium in the merger of two neutron stars", *Nature* 574, 497-500 (2019). DOI:10.1038/s41586-019-1676-3.
[22] L. Y. Zhang, J. J. He, R. J. deBoer, M. Wiescher 等（JUNA合作组，北师大/中国原子能科学研究院/圣母大学等），"Measurement of 19F(p,gamma)20Ne reaction suggests CNO break-out in first stars", *Nature* 610, 656-660 (2022). DOI:10.1038/s41586-022-05230-x.
[23] S. Wang, P. W. Zhao, J. Meng（北京大学/山东大学），"Relativistic model-free prediction for neutrinoless double beta decay at leading order", *Physics Letters B* 855, 138782 (2024).
[24] Y. Zhu, J. Barnes, K. A. Lund, T. M. Sprouse, N. Vassh, G. C. McLaughlin, M. R. Mumpower, R. Surman，"Modeling nuclear mass uncertainties for r-process nucleosynthesis", *EPJ Web of Conferences* 260, 03004 (2022). DOI:10.1051/epjconf/202226003004.
[25] E. Leistenschneider 等（NSCL/MSU LEBIT、TRIUMF TITAN），"Precision Mass Measurements of Neutron-Rich Scandium Isotopes Refine the Evolution of N=32 and N=34 Shell Closures", *Phys. Rev. Lett.* 126, 042501 (2021). arXiv:2006.01302.