跳到正文
返回核结构
至真探索 · 核结构

原子核电磁矩-从成绩单到试金石

2026-08-16

原子核电磁矩：当成绩单变成试金石

原子的核，转起来像个微小的陀螺。这个陀螺带多少磁，叫磁矩；它的电荷分布有多"扁"，叫电四极矩。几十年来，这两个数在物理学家手里是张成绩单——每测一个新核，就拿去跟壳模型对一下。对上了，说明波函数编得没错；对不上，就调一调那个"有效g因子"。这个因子的物理来源——介子交换、核心极化——理论上有解释，但几十年里没人能直接从第一性原理把它算出来，它始终是个要按数据定标的参数。

2024年，这张成绩单被人撕了。

CERN的ISOLDE装置把镍的五个奇质量同位素——⁵⁹Ni、⁶¹Ni、⁶³Ni、⁶⁵Ni、⁶⁷Ni——逐个送进激光束，扫出它们的超精细结构，反推出每个核的磁矩和四极矩。镍不是第一次被人测，但这次测得很细。⁵⁹Ni的磁矩是头一回测到，−0.7439(15)核磁子，括号里是误差，小到可以放心；⁶⁵Ni的磁矩0.7575(7)，比旧值0.69(6)的误差缩小了近两个量级（0.06压到0.0007）；⁶⁵Ni和⁵⁹Ni的四极矩也都是第一次见天日。测量细节相当讲究：1.4 GeV的质子打碳化铀靶，打出镍原子，用激光在离子源里电离，质量分离，再聚成一束，与一束同向激光重合——这叫共线激光谱，是CERN-ISOLDE的招牌菜（COLLAPS装置）。论文的第一作者P. Müller在德国达姆施塔特工业大学，理论部分的T. Miyagi和A. Schwenk也在达姆施塔特，实验组长W. Nörtershäuser坐镇，作者名单里还有比利时的G. Neyens、英国的B. Cheal、中国的X.F. Yang（北京大学）。

测完，他们把数据交给了当时最前沿的从头计算方法。方法的名字叫VS-IMSRG，说人话就是：从核力出发，一步步把多体问题"重正化"成壳模型的形式，中间不塞任何经验参数。结果很难看。理论算出来的磁矩系统性偏低，五个核全偏。把理论与实验的均方偏差量化一下——大了3到5倍，具体看你用哪套手征核力。这就是你知识库里那句"偏差×3-5"的出处，一字不虚。

然后他们做了一件关键的事：把"双体流"显式加了进去。

双体流是什么？想象两个核子挨得很近，中间交换着π介子。你转动这个双核子系统的磁矩，不光是单个核子自己转出来的——两个核子"联手"那部分，也就是交换介子的那个过程，也在转。老壳模型把这种贡献揉进有效g因子，谁也不知道准不准。VS-IMSRG则从手征有效场论里把这个贡献直接、显式地算出来。加进去之后，偏差立刻塌缩：最好的组合（ΔN2LOGO核力）把均方偏差从0.101压到0.022核磁子平方，压缩了约5倍；另外两套核力也压了约3倍。这篇论文的摘要原话是："双体流贡献显著改善与实验数据的一致度，均方偏差取决于所用相互作用降低3到5倍。"镍链的结果表明，在中重质量区、整条同位素链的系统性比较里，含双体流的从头计算第一次整体压过了调好参数的老壳模型——轻核区量子蒙特卡洛早就能做到这一点，镍链的意义在于把战场推到了中重质量区。

镍链不是孤例。铋的五个丰中子同位素²¹⁴⁻²¹⁸Bi，2025年也被ISOLDE扫了一遍，用另一种技术——源内共振电离激光谱（RILIS），结果发表在Physical Review C 112, 034304上，署名93人。第一兼通讯作者A. Barzakh来自俄罗斯的PNPI（经CERN协议），作者名单里能看到欧洲一大串实验室，还有中科院近代物理研究所的Z. Liu。磁矩本身没有大新闻：加和关系基本对得上——拿已知的单粒子磁矩拼出奇奇核的磁矩，跟测量值相符。真正的发现在电荷半径。铋同位素链的电荷半径在N≈132、134处出现拐点——这正是钋（Z=84）到砹（Z=85）过渡时八极形变突然冒出来的区域。铋链被钉在了"Z>82八极形变区"（钋、氡、镭一带）的南边界上。这个边界为什么值钱？因为八极形变的核，原子核内禀电偶极矩（EDM）会被放大——而EDM是寻找CP破坏新物理的最锋利探针之一（²²⁵Ra现在把上限压到|d|<1.4×10⁻²³ e·cm）。圈定八极形变区的边界，等于给未来的EDM实验画靶。这也是"铋测量"真正意义上的"新物理"分量——不在于磁矩本身，而在于它给新物理搜索选址。

如果你以为故事到此圆满，那就错了。双体流救了磁矩，但至少还有三块硬骨头。

第一块，⁴⁸Ca的10.23 MeV强磁偶极（M1）跃迁，三种测量互相打架。电子散射给B(M1)≈3.9–4.0核磁子平方；质子散射给3.45–4.1；光子核反应（γ,n）却给6.8±0.5。从头计算倾向大的那个值——但耦合簇算出的范围从7.0到10.2核磁子平方，比实验分歧本身还宽，说明理论自己也没收敛。可问题还更深：领头阶的双体流里，pion-in-flight项和seagull项几乎完全对消，剩下什么得算到更高阶才知道。谁对谁错，现在没人敢拍板。这是个还没打完的官司。

第二块，电四极矩没被救。磁矩被双体流修好了，四极矩却还欠着账。VS-IMSRG在sd壳的E2强度普遍欠算约2倍——¹⁷O的基态四极矩实验−2.558 e·fm²，理论只算到−1.40。镍链里⁶¹Ni的实验四极矩16.3 e·fm²，理论只有6.2–8.3；⁵⁹Ni更惨，实验+8.8，理论在−0.80到+0.98之间晃，连正负号都定不下来。原因大概是模型空间外面的多粒子-多空穴激发，这部分VS-IMSRG还没吃进去。磁矩和四极矩，一个被双体流救活，一个还在等更高阶的算法——两条腿，明显不同步。

第三块，双体流本身还是半成品。"3到5倍取决于相互作用"这句话反过来读就是：不同核力给的双体流大小不一样，还没收敛。轻核区更扎眼：用Norfolk势做量子蒙特卡洛算A≤10的磁结构，A=9的核（⁹Be、⁹B）里双体流贡献高达约33%，而且手征展开到N³LO时贡献反常地又变大——幂次计数到底收不收敛，是悬案。也就是说，双体流这个"解药"，剂量还没配准。

这背后其实有一条更大的暗线。双体流不只救磁矩，还救β衰变。2019年Gysbers等人的Nature Physics论文给出了一个主流解释：实验上看到的轴矢耦合常数在核里变小（约0.75倍），主要不是基本常数变了，而是双体弱流加核子关联的功劳。这个解释仍有人争论——Δ共振的贡献、更高阶弱流、在无中微子双β衰变里的规范匹配都还没完全收口——但它把"淬灭"从一个神秘现象变成了可计算的效应。磁矩和β衰变，两个看似不搭界的难题，指向同一个根源——电磁算符和弱算符都含双核子项，这是规范不变性逼出来的。正因如此，物理学家现在对"淬灭因子"格外警惕：既然M1跃迁和Gamow-Teller跃迁的双体流贡献性质完全不同，那对两者用同一个淬灭因子就站不住脚。而这个"统一淬灭因子"恰恰是无中微子双β衰变（0νββ）矩阵元计算里的常用招数。一个核矩问题，就这样捅到了粒子物理最前沿的腰眼上——你算0νββ的衰变率，得先搞清楚淬灭因子，而淬灭因子的根基，正在被磁矩数据动摇。

四极矩还有个更俗气的麻烦：要从超精细结构里反推四极矩，你得知道原子内的电场梯度（EFG），而EFG经常算不准。按Stone为IAEA编的最新汇编（INDC(NDS)-0833），17个元素的参考四极矩标准到现在都是缺的，V、Cr、Zn这些常用元素的EFG误差还挂着9%以上。新的破解路径是把放射性原子做成分子——比如RaF——分子里的电场梯度大得多，超精细劈裂一眼就能读出四极矩。2025年有人在CRIS上测出²²³RaF分子的超精细结构，从中读出短寿命放射性分子里的第一个电四极矩——不过这只是会议报告，还没正式发表。分子谱学成了四极矩测量的新门。

说到数据，得提一句核矩的"家谱"。1989年Raghavan编了第一本核矩表，2005年牛津的N.J. Stone接过棒，编了那本被引用上千次的《Table of nuclear magnetic dipole and electric quadrupole moments》，现在推荐值挂在IAEA的网站上（nds.iaea.org/nuclearmoments）。但请注意一个怪现象：NUBASE和AME——这两个核物理的"户口本"——只记自旋宇称，不记矩。矩的评价散落在ENSDF和Stone的汇编里，没有统一版本管理。所以核矩数据的权威性，其实比核质量数据弱一个档次。而且评价永远追着实验跑：⁶¹Ni的参考磁矩刚被镍论文重新审了一遍，移动了约0.0023核磁子——按新参考值自己的误差算，相当于六七个σ——说明老参考值本身在漂。做核矩的人，脚下是沙。

还有一条很少人注意的线：磁矩在强磁场里的天体物理作用。俄罗斯的V.N. Kondratyev提出，超新星和中子星并合里磁场强到10¹⁴高斯以上，核的磁矩被场钉住，Zeeman作用改变单粒子壳结构，进而改写核统计平衡的丰度——反幻数核产量被指数放大，幻数核被压制。这个机制预言⁴⁴Ti产量大增，指向INTEGRAL卫星在仙后座A和SN1987A残骸里看到的⁴⁴Ti过量（"钛问题"）。但要说明，这是Kondratyev一个人的非主流提案——⁴⁴Ti过量还有别的标准解释，学界并未普遍接受磁化这一说。另有研究指出，在10¹⁴⁻¹⁶高斯的磁化等离子体里，⁵⁴Fe的电子俘获率被Landau量子化压掉两个量级（Luo、Kusakabe和Kajino，2022年发表于EPJ Web of Conferences），直接影响超新星铁芯的质量。也就是说，磁矩不只是地上的成绩单，还是天上核合成的隐形推手。天体物理里还有几个核的矩格外金贵：⁷Be的磁矩−1.3995(5)核磁子（Nörtershäuser等人2009年用激光谱测的），是约束⁷Be(p,γ)⁸B反应率——也就是太阳硼中微子来源——的硬输入；²⁶Al的磁矩+2.804(4)、四极矩+0.265(32)b，关系到银河系里那颗最亮的1.809 MeV伽马线怎么读。可惜，⁷Be和⁴⁴Ti激发态的四极矩至今没人直接测过——空白本身，就是机会。

那中国在这张图上站哪？三条线，一根理论两根实验。

理论这条线最亮。孟杰-李剑学派的协变密度泛函理论（CDFT），把"时间反演破缺场"堂堂正正做进磁矩计算——这是相对论框架下别人没系统做过的。他们算双幻核±1的四个核（²⁰⁷Pb、²⁰⁹Pb、²⁰⁷Tl、²⁰⁹Bi）的磁矩，相对偏差做到6.1%，优于非相对论方法的13.2%；³³Mg的磁矩算到−0.9134核磁子，复现实验−0.7456(5)的物理趋势；还加入了单π介子交换流修正。这套工作是全球相对论框架下最系统的磁矩计算之一，综述发表在Frontiers of Physics 13, 132109 (2018)。

实验第一条，中国原子能科学研究院（CIAE）的β-NMR。注意一个常见误植：中国β-NMR的大本营是CIAE，不是北师大——北师大是理论组（彭婧做磁转动/反磁转动）。CIAE跟日本大阪大学的Minamisono/Matsuta组合作，在北京的串列加速器上建起中国第一台β-NMR/β-NQR谱仪，专挑质子滴线附近的核测。¹⁷F的四极矩0.076(4)b，由此判定它是质子皮而不是质子晕；²⁸P的磁矩0.309(9)是首次测定；¹²B的磁矩经Knight位移修正后精到0.99993(48)核磁子。这些成果不新，但独到——国际上没有直接对标。

实验第二条，北京大学的杨晓菲课题组，从2018年从零建激光核谱。2022年，他们主导了第一个中国领头的ISOLDE核矩测量——Sc同位素链和N=28同中子素，⁴⁹Sc的四极矩精度提升约10倍，实验数据同时拷问了大规模壳模型和VS-IMSRG。2025年，课题组拿出PLASEN谱仪：标称约100 MHz分辨率、1:200总效率，课题组认为已达国际最高水准——第三方独立比对还没人做过——已经在中国原子能院的BRIF装置上在线跑通⁹²,⁹⁵Rb。杨晓菲本人也出现在镍论文的作者名单里——中国学者已经深度参与ISOLDE的核矩测量。

但中国的短板同样硬。ab initio核矩计算——就是镍和铋论文里那种VS-IMSRG、耦合簇——本土基本空白，得靠海外华人（比如TRIUMF的B.S. Hu）撑场面。HIAF建成了：2025年10月27日首次调试出束，2026年7月21日通过工艺验收转入试运行，但电磁矩实验能力还停在"规划+预研"，一篇核矩论文都还没出。重丰中子区（N>60、N>90）的矩测量，中国几乎没有公开发表的工作。这是追赶的地图，也是机会的地图。

现在把镜头拉远，看这个领域整体的重心在怎么移。

二十年前，物理学家手里是有效g因子加有效电荷，一调就平。今天的图景完全不同：数据不再被拿去"拟合"，而是被拿去"拷问算符"。你的电磁算符里到底有没有双体流？算到第几阶才收敛？四极矩为什么还不收敛？这些问题才是镍、铋、锑、铟、铬、锕一长串同位素链真正在问的问题。ISOLDE、FRIB、RIKEN、TRIUMF这四台装置，2024–2026年一口气贡献了镍（磁矩+四极矩）、铋（电荷半径+电磁矩+八极边界）、铬-61（N=40反转岛西界，自旋修订为1/2，μ=+0.541(6)）、铟（逼近¹⁰⁰Sn，双重幻核画像）、钍-229（离子阱里μ=+0.365(3)、Q=3.11×10⁻²⁸ m²，为核光钟铺路）、锕（TRIUMF）等一大批结果。每一条同位素链，都是一次对从头计算的公开考试。

而考试还在出题。双体流的高阶项何时收敛？⁴⁸Ca的实验矛盾谁来仲裁？重核区（Z>50）的从头验证何时到来？四极矩的符号问题怎么解决？0νββ的淬灭因子到底怎么处理？每一个开放问题后面，都站着一个具体的测量和一个具体的算法。

回到开头那个陀螺。一个磁场强度，一个电荷形状，两个小数字，牵出来的是整张核结构地图——从双体流算符到八极形变区边界，从g_A淬灭到0νββ，从太阳中微子到银河的钛。这张地图还在重画。而画图的笔，一半在CERN的激光束里，一半在超级计算机的算法里，剩下那一小截，正在从中国自己的装置里长出来。

---

参考文献

  1. P. Müller, S. Kaufmann, T. Miyagi 等（TU Darmstadt/MPI-K/GSI/CERN/KU Leuven/北京大学等27人）, *Electromagnetic moments of the odd-mass nickel isotopes ⁵⁹⁻⁶⁷Ni*, Physics Letters B 854 (2024) 138737; arXiv:2405.13668
  2. A. Barzakh 等（PNPI/JAEA/CERN/KU Leuven/IMP兰州等93人）, *Charge radii and electromagnetic moments of ²¹⁴⁻²¹⁸Bi: Exploring the "southern" border of the Z>82 octupole-deformation region*, Physical Review C 112 (2025) 034304; DOI 10.1103/sltw-g4d6
  3. T. Miyagi, X. Cao, R. Seutin, S. Bacca, R.F. Garcia Ruiz, K. Hebeler, J.D. Holt, A. Schwenk, *Impact of Two-Body Currents on Magnetic Dipole Moments of Nuclei*, Physical Review Letters 132 (2024) 232503; arXiv:2311.14383
  4. B. Acharya, B.S. Hu, S. Bacca, G. Hagen, P. Navrátil, T. Papenbrock, *Magnetic Dipole Transition in ⁴⁸Ca*, Physical Review Letters 132 (2024) 232504; arXiv:2311.11438
  5. C. Brase, T. Miyagi, J. Menéndez, A. Schwenk, *Two-body currents at finite momentum transfer and applications to M1 transitions*, Physical Review C 113 (2026) 014317; arXiv:2504.08711
  6. P. Gysbers, G. Hagen, J.D. Holt 等, *Discrepancy between experimental and theoretical β-decay rates resolved from first principles*, Nature Physics 15 (2019) 428–431; arXiv:1903.00047
  7. S. Lechner 等, *Electromagnetic moments of the antimony isotopes ¹¹²⁻¹³³Sb*, Physics Letters B 847 (2023) 138278; arXiv:2311.01110
  8. S.R. Stroberg, J. Henderson, G. Hackman, P. Ruotsalainen, G. Hagen, J.D. Holt, *Systematics of E2 strength in the sd shell with the valence-space in-medium similarity renormalization group*, Physical Review C 105 (2022) 034333; arXiv:2203.15342
  9. R. Seutin, O.J. Hernandez, T. Miyagi, S. Bacca, K. Hebeler, S. König, A. Schwenk, *Magnetic dipole operator from chiral effective field theory for many-body expansion methods*, Physical Review C 108 (2023) 054005; arXiv:2308.00136
  10. L. Lalanne 等, *⁶¹Cr as a doorway to the N=40 island of inversion*, Physical Review C 112 (2025) L031301; arXiv:2409.07324
  11. J. Karthein 等, *Electromagnetic properties of indium isotopes illuminate the doubly magic character of ¹⁰⁰Sn*, Nature Physics 20 (2024) 1719; arXiv:2310.15093
  12. R. Chambers-Wall, A. Gnech, G.B. King, S. Pastore, M. Piarulli, R. Schiavilla, R.B. Wiringa, *Magnetic structure of A≤10 nuclei using the Norfolk nuclear models with QMC*, Physical Review C (2024); arXiv:2407.04744
  13. V.N. Kondratyev, *Nucleosynthesis at Strong Magnetization and the Titanium Problem*, Physics of Particles and Nuclei 50 (2019) 576–580
  14. S. Bai, Á. Koszorús, B.S. Hu, X.F. Yang 等, *Electromagnetic moments of scandium isotopes and N=28 isotones in the distinctive 0f7/2 orbit*, Physics Letters B 829 (2022) 137064; arXiv:2203.09111
  15. J. Li, J.X. Wei, J.N. Hu, P. Ring, J. Meng, *Relativistic description of magnetic moments in nuclei with doubly closed shells plus or minus one nucleon*, Physical Review C 88 (2013) 064307; arXiv:1307.2044
  16. H.R. Hu, X.F. Yang 等, *Development and characterization of a high-resolution, high-sensitivity collinear resonance ionization spectroscopy system*, Science Bulletin 70 (2025) 2721
  17. G. Zitzer, J. Tiedau, Ch.E. Düllmann, M.V. Okhapkin, E. Peik, *Laser spectroscopy on the hyperfine structure and isotope shift of sympathetically cooled ²²⁹Th³⁺ ions*, Physical Review A 111 (2025) L050802; arXiv:2504.00974
  18. N.J. Stone, *Table of nuclear magnetic dipole and electric quadrupole moments*, Atomic Data and Nuclear Data Tables 90(1) (2005) 75–176
  19. Y. Luo, M. Kusakabe, T. Kajino, *Impacts of Strongly Magnetized Degenerate Plasma on the Electron-capture Rates*, EPJ Web of Conferences 260 (2022) 07004
  20. R. Li, A. Teigelhöfer 等, *Collinear laser spectroscopy on neutron-rich actinium isotopes*, Physical Review C 111 (2025) 054320; arXiv:2503.13732