幻数变迁
1949年，Mayer和Jensen用壳模型解释了原子核的幻数——2、8、20、28、50、82、126。这些数字标记着核子填满一个壳层的时刻，壳与壳之间有能隙，像楼梯间的平台。幻数核特别稳定，球形，不易变形。这个图景对稳定线附近的核非常好使。
然后加速器造出了远离稳定线的核。幻数开始动了。
有的消失了，有的冒了出来。同一个数字，在核图的一个角落是铁律，挪到另一个角落就不管用了。这不是壳模型错了，是壳模型背后的参数在变——核力各分量随质子和中子填充方式的改变，把壳间距重新排了一遍。
旧幻数的崩塌
最出名的例子是N=20的"反转岛"。按标准壳模型，中子数20应该闭合sd壳层，核保持球形。但在丰中子的Ne、Na、Mg同位素里，³²Mg的2⁺激发态出奇地低，B(E2)出奇地大——这不是球形核该有的样子。Warburton等人在1990年代的壳模型计算指出，这些核的基态已经被跨壳激发组态（侵入态）占领，正常sd壳组态反而成了激发态。反转岛的南界一路推进：从³²Mg到²⁹F，2023年终于到了²⁸O。
Kondo等人在RIKEN的RIBF上，用SAMURAI磁谱仪和NeuLAND中子阵列，通过²⁹F的单质子移除反应产生了²⁸O，看到它通过四中子衰变消失[1]。谱因子测量给出C²S=0.48，远大于sd闭合壳层预言的0.13——²⁸O的中子根本没老老实实待在sd壳层里，pf壳层的侵入组态占了主导。N=20壳隙从Z=14处的7 MeV缩小到Z=8处只剩2 MeV。本该是双幻核的²⁸O，连幻数的影子都没有。
N=28也没好到哪里。从双幻核⁴⁸Ca只移走两个质子，⁴⁶Ar的E(2⁺)就掉了2255 keV。⁴⁴S出现形变与球形的共存，⁴²Si变成了扁椭形变转子。N=20和N=28两个反转岛在Mg同位素链上甚至合流了。
新幻数的诞生
旧的去，新的来。最干净的新幻数故事发生在Ca同位素链。
N=32的新幻数由Wienholtz等人在CERN ISOLDE上用ISOLTRAP质量测量首次确认[2]：⁵²Ca的双中子分离能在N=32处出现明显的拐折。Steppenbeck等人在RIKEN用⁵⁵Sc的质子敲出反应测到了⁵⁴Ca的2⁺激发态，确认N=34也是一个新壳隙[3]。Ca同位素链的E(2⁺)在N=32和N=34处都出现了峰值——双幻特征。
但N=32幻数有多"牢"？Koszorus和杨晓菲等人在ISOLDE上用共线激光谱测了⁵²K的电荷半径，结果出乎意料[4]：K同位素的电荷半径在N=32处没有表现出幻核应有的"收拢"特征。仅仅少一个质子，N=32的幻性就从Ca到了K就不灵了。这让理论家头疼——DFT和ab initio的计算都在这个点上卡了壳。
2025年，杨晓菲课题组又测了丰中子Sc同位素的电荷半径[5]。Sc（Z=21）比Ca多一个质子，电荷半径的演化趋势却跟Ca截然不同。0f₇/₂轨道的seniority对称性在Sc同位素中展现出奇偶震荡，但Ca那套成功的Fayans-DFT理论放到Sc上就不好使了。一个质子改变了整条同位素链的面貌。
张量力：幻数变迁的幕后推手
为什么幻数会变？Otsuka在2010年的PRL里给出了核心答案[6]：核力中的张量力分量通过单极相互作用改变单粒子能。简单说，π0f₇/₂质子轨道和ν0f₅/₂中子轨道之间有张量吸引——质子填进f₇/₂，就把中子的f₅/₂往下拉。所以从Ca到Ni，随着π0f₇/₂被填满，ν0f₅/₂越来越低，N=34壳隙越来越窄，最终在⁶²Ni完全消失。
2026年，Kumar、Miyagi和清水则范用VS-IMSRG从头算框架给出了定量分解[7]：N=34壳隙的消失中，NN张量力贡献约83%，3N张量力贡献约17%。3N力虽然占比小，但不可忽略——它的主要来源是双π交换中的c₃项。这个结论在两种手征EFT相互作用下都成立，不依赖参数选择。
Otsuka、Gade、Sorlin、Suzuki和Utsuno在2020年的Reviews of Modern Physics长文综述里把这些机制系统化了[8]：壳演化是由中心力、张量力、二体自旋-轨道力和三体力的单极分量共同驱动的。每个分量的贡献可以从有效单粒子能（ESPE）的位移中逐一读出。
幻数的起源：从核力到壳结构
幻数不是天上掉下来的。它从核力的哪个分量里长出来？2025年，中山大学尧江明团队给出了最直接的回答。
Ding、Wang、尧江明、Hergert、梁昊梓和Bogner在PRL上发表了"From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure"[9]。他们用VS-IMSRG演化核力，追踪从高分辨率到低分辨率的壳结构变化，发现了一个普适的转化过程：高分辨尺度下自旋-轨道劈裂小，自旋对称性强；随着尺度降低，SO劈裂显著增大，赝自旋对称性浮现——传统幻数在这个转化中自然涌现。分解各力的贡献，三体力是增强SO劈裂的主导因素，其中双π交换的c₃项占3N总贡献的80%–90%。用相对论Hartree-Fock计算在⁹⁰Zr中验证，得到一致的图像。
APS Physics在2026年2月以"Exposing Nuclear Magic"为题专文评论了这项工作[10]：幻数不是壳模型的人为拟合参数，而是核力各分量在粗粒化过程中的自然产物。
双幻核的坚守与侵蚀
⁷⁸Ni（Z=28, N=50）是核图上最丰中子的双幻候选者。Taniuchi等人在2019年的Nature上给出了⁷⁸Ni的第一张谱学画像[11]：2⁺态能量支持双幻性，但旁边不远处就有一条形变带在等着——形状共存。71位作者的跨国团队，用RIKEN RIBF的RIBF装置完成了这个实验。
但N=50幻数能撑到多远？Hu、Sun、Hagen、Jansen和Papenbrock在2024年用耦合团簇从头算从⁷⁸Ni一路算到了⁷⁰Ca[12]：N=50幻数向⁷⁰Ca方向逐渐被侵蚀，⁷⁶Fe、⁷⁴Cr、⁷²Ti的基态趋向形变，⁷⁴Cr被确认为N=50反转岛内第一个同中子异位素。⁷⁰Ca的势能面变得平坦，没有明确定义的极小值——这个核到底什么形状，目前的计算方法还给不出可靠答案。
同位旋对称的反转岛
2025年最大的惊喜来自N=Z核。
Ha、Recchia、Lenzi、Iwasaki等人在Nature Communications上报道了⁸⁴Mo和⁸⁶Mo第一个2⁺态的寿命测量[13]。⁸⁴Mo（N=Z=42）的B(E2)高达1740 e²fm⁴，⁸⁶Mo（N=Z+2=44）只有707 e²fm⁴——两者之间出现了急剧的结构跳变。壳模型DNO-SM计算揭示，⁸⁴Mo的基态由8粒子-8空穴侵入组态主导，质子和中子同步跨越壳隙——这不是此前只在丰中子侧看到的反转岛，而是质子和中子一起动的"同位旋对称反转岛"。此前所有反转岛都在中子富集侧（N=8, 20, 28, 40），⁸⁴Mo是第一个在质子富集侧被确认的。
关键机制：二核子力算不出正确的B(E2)，必须加三核子力。3N力缩小了N=Z=40的谐振子壳隙（有利于侵入组态），同时恢复了N=50的自旋-轨道壳隙（让⁸⁶Mo回归近球形）。这是3N力在变形核中发挥关键作用的第一个明确案例。
Z=14：新质子幻数
2025年7月，Xing（邢玉明）、Luo（罗宇飞）、Zhang（张玉虎）等40位作者在PRL上发表了²²Si的质量测量[14]。在兰州重离子加速器冷却储存环上，用改进的Bρ-defined等时性质量谱法，首次测到了质子滴线核²²Si的质量。单质子分离能Sₚ=+1412(114) keV，双质子分离能S₂ₚ=+229(54) keV——²²Si是束缚的，Si元素的质子滴线位置由此确定。
更重要的是，质量差异分析表明²²Si展现出与其镜像伙伴²²O类似的双幻特征。Gamow壳模型计算证实，²²Si中价质子的平均占据数与²²O中价中子几乎相同——Z=14是一个新质子幻数。但镜像对称性有轻微破缺：²²Si的质子分布比²²O的中子分布更延展，因为非束缚的π2s₁/₂轨道有微小贡献。
全局视野
叶沿林、杨晓菲、樱井博仪和胡柏山在2024年的Nature Reviews Physics综述中描绘了奇特核物理的全景[15]：过去三十年，核版图从3000多种扩展到远未到达边界，Halo结构、集团结构、壳演化、新集体模式、与连续态的强耦合——这些都是弱束缚开放量子系统的特征。壳演化与幻数变迁是其中最核心的现象之一，因为它直接反映了核力各分量在改变环境下的行为。
刘洪娜在2026年的Frontiers in Physics小型综述中聚焦了N=20/N=28的消失和N=32/N=34的涌现[16]，通过敲出反应的光谱因子追踪壳隙的演化。她提出了几个关键的未来实验：⁴⁸S的在束γ谱学将判别N=32亚壳在低Z方向的命运；RIKEN的HYPATIA升级装置和中国的HIAF将把γ谱学推到⁶⁰Ca甚至N=126的丰中子区域。
Koszorus和杨晓菲的⁵²K电荷半径已经警告过我们：N=32的幻性高度局域化，在Ca最强，到K就被挑战。N=34也类似——Z < 20时快速减弱。幻数不是一堵墙，更像一道帘，风一吹就变形，风大了就吹开。但帘子的位置不是随机的，它由核力的张量力、三体力、自旋-轨道力精确地锚定。把这些力的参数搞对了，帘子在哪里、怎么变形，原则上都能算出来。
现在算对了大部分，但还不是全部。⁷⁰Ca到底什么形状？⁴⁸S是不是双幻核？N=126在丰中子侧还能撑住吗？这些问题等待FRIB、HIAF和HYPATIA来回答。
---
参考文献
[1] Y. Kondo *et al.*, "First observation of ²⁸O," *Nature* 620, 965–970 (2023). arXiv:2306.08348
[2] F. Wienholtz *et al.*, "Masses of exotic calcium isotopes pin down nuclear forces," *Nature* 498, 346–349 (2013).
[3] D. Steppenbeck *et al.*, "Evidence for a new nuclear 'magic number' from the level structure of ⁵⁴Ca," *Nature* 502, 207–210 (2013).
[4] Á. Koszorus, X. F. Yang *et al.*, "Charge radii of exotic potassium isotopes challenge nuclear theory and the magic character of N=32," *Nature Physics* 17, 439–443 (2021).
[5] S. W. Bai *et al.*, "Charge Radii of Neutron-Rich Scandium Isotopes and the Seniority Symmetry in the 0f₇/₂ Shell," *Phys. Rev. Lett.* (2025).
[6] T. Otsuka, T. Suzuki, M. Honma, Y. Utsuno, N. Tsunoda, K. Tsukiyama, and M. Hjorth-Jensen, "Novel Features of Nuclear Forces and Shell Evolution in Exotic Nuclei," *Phys. Rev. Lett.* 104, 012501 (2010).
[7] A. Kumar, T. Miyagi, and N. Shimizu, "NN and 3N Tensor Force in the N=34 Shell Evolution: An Ab Initio Perspective," arXiv:2605.16822 (2026). University of Tsukuba.
[8] T. Otsuka, A. Gade, O. Sorlin, T. Suzuki, and Y. Utsuno, "Evolution of shell structure in exotic nuclei," *Rev. Mod. Phys.* 92, 015002 (2020).
[9] C. R. Ding, C. C. Wang, J. M. Yao, H. Hergert, H. Z. Liang, and S. Bogner, "From spin to pseudospin symmetry: The origin of magic numbers in nuclear structure," *Phys. Rev. Lett.* (2025). arXiv:2504.09148. 中山大学/FRIB/东京大学.
[10] "Exposing Nuclear Magic," *Physics* 19, 11 (2026). APS Viewpoint.
[11] R. Taniuchi *et al.*, "⁷⁸Ni revealed as a doubly magic stronghold against nuclear deformation," *Nature* 569, 53–58 (2019).
[12] B. S. Hu, Z. H. Sun, G. Hagen, G. R. Jansen, and T. Papenbrock, "Ab initio computations from ⁷⁸Ni towards ⁷⁰Ca along neutron number N=50," *Phys. Lett. B* 860, 139010 (2025). arXiv:2408.07856. Oak Ridge National Laboratory.
[13] J. Ha, F. Recchia, S. M. Lenzi, H. Iwasaki *et al.*, "Abrupt structural transition in exotic molybdenum isotopes unveils an isospin-symmetric island of inversion," *Nature Communications* 16, 10631 (2025).
[14] Y. M. Xing, Y. F. Luo, Y. H. Zhang, M. Wang, X. H. Zhou, J. G. Li *et al.*, "Z=14 Magicity Revealed by the Mass of the Proton Dripline Nucleus ²²Si," *Phys. Rev. Lett.* 135, 012501 (2025). 中国科学院近代物理研究所/北京大学/INFN/GSI/MPIK等.
[15] Y. L. Ye, X. F. Yang, H. Sakurai, and B. S. Hu, "Physics of exotic nuclei," *Nature Reviews Physics* 7, 21–37 (2025).
[16] H. N. Liu, "Exploring the migration of magic numbers in neutron-rich systems via knockout reactions," *Frontiers in Physics* 14, 1778872 (2026). 北京师范大学.