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

锕系核低能结构-梨形核集体运动与宇称态的五维版图

2026-07-22

如果核物理有一个"为什么我们还在研究这个"的终极回答，在锕系——从钍（Z=90）到锎（Z=98）的那条元素链上——那个回答就是："因为这些核不是球。"

不只是椭球。是梨。一头大、一头小、在空间中旋转时像一颗被压缩的水滴。量子力学的平均场本来应该把它们压回对称——但在大约N≈134中子数和Z≈88质子数附近的特定核区，单粒子轨道的角动量对齐得如此完美，以至于整个核集体地、自发地破缺了空间反射对称性。这就是八极形变。梨形。

这个领域在2024-2026年正在经历一轮罕见的多维度刷新——从日本北海道到英国利物浦到中国华北再到意大利的Legnaro，理论和实验以四种不同的数学语言逼近同一组数据。这次刷新最不寻常的地方在于：它不是某一维度的单点突破，而是实验（ISOLDE库仑激发）、集体模型（IBM映射、几何模型）、壳模型（PNC-CSM推转壳）、以及新物理连接（Schiff矩-EDM、贝塔衰变）——四条线同时在动，而且部分交汇了。

先从实验版图最稳固的一角开始。

利物浦大学的Peter Butler在2024年发表了ISOLDE Miniball合作组在氡（Z=86）和镭（Z=88）同位素链上多年库仑激发实验的系统性综述，发表在*Physica Scripta*（DOI:10.1088/1402-4896/ad22c9）。HIE-ISOLDE用4.23-5.08 MeV/u的二次放射性束去轰击¹²⁰Sn和⁶⁰Ni靶，24个Miniball高纯锗探测器加上分段硅粒子探测器，耦合道GOSIA分析。核心结论清晰而稳固：²²²⁻²²⁶Ra同位素的E3跃迁矩在自旋变化时大致恒定——这是旋转梨形核的确凿签名——一只梨在转动，形状不变。²²⁸Ra不再表现出同样的行为，意味着在Ra同位素链上，静态八极形变在质量226附近达到峰值、然后在此处开始软化。

但更有意思的是氡——Z差两个质子。²²²,²²⁴,²²⁶Rn用的是完全相同的实验条件——但它们不是静态梨形。它们是八极振动的：核在反射对称和不对称两个形状之间来回振荡，但基态本身不是对称破缺的。两个质子——八极形变的质变——就在这两个元素之间发生。Butler在综述中写道：为什么这些核的E3跃迁矩在实验中测得与理论比较——仍然没有哪个特定的微观模型可以一致地解释所有数据。

Butler的综述直接给出了物理的终端：电偶极矩搜索。一个静态梨形核——如²²⁵Ra——拥有一个基态加一个能量仅高出约50 keV的、自旋相同但宇称相反的态。在宇称双重态之间，极化力——任何微小的、破坏时间反演和宇称的新物理相互作用——都会将这两种态相混合，从而在实验室坐标系中产生一个可测的电偶极矩。增强因子——由于双态间距极小而产生的奇异敏感性——约100到1000倍。Butler列举了候选核的当前状态：²²⁵Ra（双态间距约50 keV）是当前最佳候选，²²³,²²⁵Ra和²²¹Rn已进入原子物理EDM实验的正式候选名单。²²⁹Pa曾被寄望拥有约200 eV间距的5/2宇称双重态——但Butler在综述中明确指出最近的工作"对这种双重态的存在提出了质疑"，并提议进行新的测量来重新评估这项关键输入。

现在往新物理的一端再推一步——从梨形核到Schiff矩——因为这个连接在2025年经历了一场理论升级。

哥伦比亚大学的Jonathan Engel在2025年的*Annual Review of Nuclear and Particle Science*（Vol. 75, pp.129-151；arXiv:2501.02744）上发表了一篇全面综述——从Schiff定理（电子对核EDM的屏蔽）讲起，走过所有主要的核结构模型（壳模型、密度泛函、集体模型），然后停在了最新的前沿：从头算方法在计算Schiff矩上的应用。这是关键一步——因为传统的Schiff矩计算——比如用Woods-Saxon或Skyrme平均场——对八极形变参数β₃的微小变化极其敏感，而β₃本身在核理论中有至少30-50%的模型依赖不确定性。从头算方法——IMSRG（介质内相似重整化群）和耦合簇——把不确定性从平均场级别缩小到两核子-三核子相互作用级别。Engel的综述发表不到半年，2024年Flambaum和Mansour在*Physical Review C* 111, 055501 (2024)（arXiv:2411.18943）发表了一个补充性的半解析方法——从E1跃迁半衰期直接提取内禀Schiff矩和内禀EDM——避免了完整的核多体计算，但以极高的不确定性（估计超过一个量级）为代价。

这两篇独立论文代表了两条分化路线：一条朝"精"走——从头算方法论承诺系统性可控的误差；另一条朝"快"走——解析/半解析估计法可以快速扫描多个核素但精度低。

MohanMurthy、Silwal和Winger在2024年的*Hyperfine Interactions* 245, 64里做的就是这个扫描——用Möller-Nix质量模型系统检查了所有已知核素的四极和八极形变参数，然后挑出了48个同时具有显著四极和八极形变的同位素作为EDM候选核的优选清单。他们的Top 3是²²³,²²⁵Fr和²²³Ra——都是碱土金属或碱金属，可以用激光冷却（Fr和Ra已经被冷却到微开尔文量级）——这是原子EDM实验的一个硬性技术条件。

现在把镜头转回核结构理论——从最宏观的集体模型到最微观的壳模型——看同一组数据如何在四条理论线上交汇。

北海道大学的Kosuke Nomura在2025年发表了他在映射玻色子模型（mapped IBM）上的关键升级（arXiv:2507.23282, 2025）。他的做法是从Gogny-D1M能量密度泛函的约束Hartree-Fock-Bogoliubov计算出发——在二维（β₂-β₃）势能面上扫描——然后用一组映射规则把这个平均场势能面映射到一个由s、p、d和f玻色子组成的IBM哈密顿量。所有参数由映射决定——不需要额外实验拟合。Nomura在2025年的关键创新是引入了p（偶极）玻色子——以前的sdf-IBM没有p玻色子，这意味着它无法描述电偶极（E1）跃迁和固有偶极矩——而这两者正是梨形核的签名物理量。引入p玻色子后——对²¹⁸⁻²³⁰Ra和²²⁰⁻²³²Th——负宇称yrast带和E1跃迁的描述从"质"上就比之前好——不是参数精调的结果，而是新的自由度天然地解耦了那些不能用一个纯八极场描述的集体激发。

Nomura在2025年*EPJ Web of Conferences* 329上发表了全球八极关联系统学——用sdf-IBM扫描了从镧系到锕系全象限的八极效应。最强八极关联区域一如既往地出现在质子数Z≈56、88和与它们成对的特定中子数N≈34、56、88、134——这些是单粒子能级中满足Δl=Δj=3ħ耦合条件的核子数。在锕系区——铀同位素链（²²²⁻²³²U）的八极形变在²²⁶U–²²⁸U处达到最大值，此后软化——到²³²U之后不再出现稳定八极形变。一个有趣的预言是"双八极态"——由两个f玻色子（f²）耦合而成的0⁺激发带——Nomura的模型在几个锕系核中都预言了这种态的激发能在约1.5-2.5 MeV附近——如果实验有足够的灵敏度在未来测量到它，这将是八极对涨落的首次直接观测。

第二套集体模型——来自Budaca组（INPC 2025）——不是从平均场做映射，而是直接在四极-八极集体坐标上求解含时的玻尔哈密顿量（所谓几何模型）。它的一个独特优势是能鉴别相变临界点：在A≈224-228的Ra和Th区域，Budaca发现了一个在β₂-β₃平面上从静态八极形变过渡到八极振动的平衡临界点——这类似于四极形变中著名的X(5)临界点对称性——这是第一次在实验能谱和跃迁数据中间接地"看到"八极形变中的形状相变。

第三套——Khalaf等人（2025, *International Journal of Modern Physics E*）——用一个三参数集体转动公式拟合了²²⁴,²²⁶Ra, ²²⁶,²²⁸,²³²Th, ²³⁶,²³⁸U, ²³⁸,²⁴⁰,²⁴⁴Pu的基态带和交替宇称带。他们发现了一个精准到了极端的实验事实：²³⁶U和²³⁸U的基态转动带——直到自旋24ħ——的跃迁能偏差在3 keV以内。这些核的质子数、中子数和形变都不同——但它们的基态带跃迁能近乎相同。这不是统计波动——它是"相同带"现象在正常形变锕系中的存在证据——一个在超形变核中被已知但很少在正常形变锕系中被系统化研究的效应。

第四套理论线进入微观——南京航空航天大学贺晓涛的课题组和广西师范大学章骏的合作——在2024年*原子核物理评论*（NPR）上发表的推转壳模型粒子数守恒方法（PNC-CSM）研究。它处理的是同一个物理——八极形变——但用的是完全不兼容的数学：不把核当成一系列四极和八极形变的玻色子或几何波函数——而是当成几百个核子在变了形的Nilsson势场中运动、费米面附近的核子被两两配对（对关联）、转动让科里奥利力逐渐瓦解配对、费米面附近单个核子的角动量轨道被转动拉到同一条轴上——然后整个系统的能量最小值从一种形变蹦到另一种形变。

贺晓涛组在Nilsson势中显式引入八极关联项——在U和Pu同位素链上做PNC-CSM计算（²³⁶,²³⁸U，²³⁸,²⁴⁰Pu）——他们得到了一系列实验可对标的结论。第一：转动惯量在低频率区与实验值一致——低频区的转动惯量比反射对称核偏高约10-20%——因为八极形变削弱了质子对关联。第二：八极形变越大——对关联削弱越显著。第三：奇A核²³⁷U和²³⁹Pu的s=−i内禀带转动惯量比相邻偶偶核高50-60%——这个巨大增幅不是来自八极——而是来自泡利堵塞：非成对的中子占了一条离费米面很近的、具有高角动量的单粒子轨道（这是内禀带的定义）——它"挡"住了其他核子的配对相位空间——等效于有效对间隙的剧烈下降。第四：博戈留波夫准粒子的阻挡机制——不是单个轨道——而是推转壳模型自动包含了所有靠近费米面的轨道在对关联破缺过程中重新排列的集体效应。转动导致对关联衰减、八极形变也导致对关联衰减——这两个效应在低频区是同方向的——推转壳模型能把它们同时量化。

有一个与此高度互补的模型——在同一组物理问题上从几何-IBM的路走过来——Efimov、Koval、Izosimov和Usmanov在2025年*Nuclear Physics A*上用IBM的有效转动惯量来研究重核和超重核的回弯现象。他们发现²²⁰Th、²⁴²Pu和²⁴⁴Pu展示了回弯行为——并且在IBM框架下回弯的来源可以追溯到高自旋下玻色子系统的集体激发——等效于微观语言中i₁₃/₂中子轨道和h₁₁/₂质子轨道的转动顺排。Alhendi等人2024年的工作进一步用一个加了对衰减修正的指数模型成功复现了²³⁰Th, ²³²Th, ²⁴⁰Pu, ²⁴²Pu等核的转动惯量回弯。

现在环顾全局——五个维度各自送来了什么样的一幅图景。

实验维度（A）：ISOLDE Miniball给了稳固的Ra/Rn——明确了静态vs振动的边界。但更重的锕系——Th和U——所需的放射性束——统计量极低——目前还没有类似的高质量库仑激发数据。EDM搜索需要核结构的Schiff矩输入——而前端的核实验（明确确定²²⁹Pa的宇称双重态是否真实存在）尚未完成。这是一个急迫且清晰的实验空白。

理论维度（B）：四套框架并行而互补——微观（PNC-CSM）、平均场映射（IBM）、几何（Bohr哈密顿量）和从头算（IMSRG/耦合簇、用于Schiff矩）——但它们在同一个八极形变参数（β₃）上的误差棒仍然大到涵盖彼此的中心值——而且不同框架之间缺乏系统性的联合基准比较。

天体物理维度（C）：极其薄弱。八极形变理论上会影响裂变位垒的高度和形状——进而影响r过程中裂变循环产物的质量分布——但还没有天体物理观测可以单独约束八极形变的强度。中子星方面：八极形变会影响中子星内壳层的核物质不均一性（所谓的"核面食"——nuclear pasta）——但这里的理论连接极其间接。目前的天体物理限制是空白——不是"还没有精确"——而是"根本就没有"。

新物理维度（E）：EDM-Sffiff矩——这是八个维度里唯一和锕系/镧系低能核结构有直接定量联系的超越标准模型方向。静态梨形核——²²⁵Ra——把宇称双重态的间距从典型核的几MeV压缩到了约50 keV——这意味着EDM灵敏度被放大了两个到三个量级。但反过来——也需要核结构理论把Schiff矩的不确定性从当前的"一个量级以上"压缩到——最好低于50%——才能让EDM实验导出有意义的CPV参数约束。这一点Engel在2025年的综述中非常明确：从头算（IMSRG/耦合簇）是目前唯一可能提供这种量级精度缩减的路径——它还需要几年的时间来覆盖足够多种锕系核。

中国维度（D）：两个在方法上正交的原创贡献——贺晓涛/章骏的PNC-CSM（壳模型、推转、对关联）和王媛媛/张振华的RAT-PRM（粒子-转子、非轴对称、宇称双重带）——各自覆盖了锕系低能核结构不同但互补的方面。贺晓涛组做的是偶偶核和奇A核的基态带和交替宇称带的转动惯量、对关联和八极——依赖Nilsson势；张振华组做的是特定的奇质量核（¹⁴⁵Ba）的宇称双重带的精密谱学——依赖协变密度泛函势。两条线没有会合——贺晓涛不用PRM、张振华不用PNC-CSM——这是机会也是裂缝。如果能将PNC-CSM的对关联-转动定量结果与RAT-PRM的谱学-旋称劈裂分析统一在一个框架中——裂口将变成桥梁。

这些多线叙事在2024-2026年的不寻常之处——它们不是随机分布的独立发现——而是围绕一个共同的核心物理题目的收敛。那个题目是：在量子壳层和集体运动的不稳定平衡下——一个核的形状能维持多久、在多大角动量下它崩盘、它的空间反演对称性被打破后的精确签名是什么——以及这种破缺是否可以通过一个放大镜（梨形核→Schiff矩→EDM）看到超越标准模型的新物理。

目前——最紧的链接在实验端（Butler/ISOLDE）和新物理端（Engel/Schiff矩）之间——核实验能确定哪个候选核拥有最小的宇称双态间距（即最强的EDM增强），核理论能告诉你在给定间距下Schiff矩的值是多少，然后原子物理EDM实验用冷原子阱去测。但这条链上有三个弱点：第一——²²⁹Pa的双态间距是否真的约200 eV尚未被确认（Butler 2024对这一点提出了质疑）；第二——Schiff矩的从头算精度目前还不足以将EDM实验的不确定度压缩到有意义约束CPV参数的量级；第三——从头算方法目前只能覆盖最轻的锕系——到Ra和Th——而较重的核（Pa, U, Pu）——多数最具八极形变潜力的候选核——还在计算能力的范围之外。

这是三道锁——每道锁锁住一段信息流——但每道锁的钥匙已经在打了。

---

参考文献

  1. P.A. Butler (Univ. of Liverpool), "Studies of reflection asymmetry in heavy nuclei," *Physica Scripta* 99, 032501 (2024). DOI:10.1088/1402-4896/ad22c9. — ISOLDE Miniball综述：²²²⁻²²⁶Ra静态梨形（E3矩=常数），Rn八极振动；²²⁵Ra宇称双重态~50 keV→EDM 100-1000×增强；²²⁹Pa双重态存疑。Butler，利物浦大学。
  1. K. Nomura (Hokkaido Univ.), "Dipole excitations in quadrupole-octupole deformed nuclei from the mapped interacting boson model," arXiv:2507.23282 (2025). — spdf-IBM引入p玻色子，首次定量描述E1跃迁和固有偶极矩；应用于²¹⁸⁻²³⁰Ra和²²⁰⁻²³²Th。Nomura，北海道大学。
  1. K. Nomura, "Global study of octupole collectivity with the mapped sdf-IBM," *EPJ Web of Conferences* 329, 02010 (2025). — 全球八极关联系统学：N≈134八极峰、铀²²⁶U-²²⁸U峰→²³²U以上消失；"双八极态"激发态预言。同前。
  1. J. Engel (Columbia Univ.), "Nuclear Schiff Moments and CP Violation," *Annual Review of Nuclear and Particle Science* 75, 129-151 (2025). arXiv:2501.02744. — 从头算Schiff矩（IMSRG、耦合簇）综述；呼吁压缩核结构的不确定性至<50%以服务EDM实验。Engel，哥伦比亚大学。
  1. V.V. Flambaum, H. Mansour, "Enhanced nuclear Schiff and electric dipole moments in nuclei with an octupole deformation," *Physical Review C* 111, 055501 (2024). arXiv:2411.18943. — E1半衰期→内禀Schiff矩/EDM半解析提取；静态及振动八极效应；不确定度~一个量级。Flambaum（新南威尔士大学）, Mansour（开罗大学）。
  1. U. MohanMurthy, U. Silwal, J.A. Winger, "A survey of nuclear quadrupole deformation in order to estimate the nuclear MQM and its relative contribution to the atomic EDM," *Hyperfine Interactions* 245, 64 (2024). arXiv:2401.11055. — Möller-Nix扫描→48候选核；²²³,²²⁵Fr+²²³Ra最优（激光冷却可行）。MohanMurthy（MIT/现多机构）、Silwal（密西西比州立大学）、Winger（密西西比州立大学）。
  1. 王媛媛、张振华、李杨（华北电力大学数理学院），"梨形原子核¹⁴⁵Ba中的宇称双重带，" *中国科学：物理学 力学 天文学* 56(2), 2026（2025年网络出版）。DOI:10.1360/SSPMA-2025-0015. — 反射不对称三轴粒子-转子模型(RAT-PRM)；两套基于不同内禀组态的宇称双重带（νh₉/₂⊕νi₁₃/₂）；旋称劈裂定量分析。
  1. 章骏（广西师范大学）、贺晓涛（南京航空航天大学），"转动、堵塞效应和八极形变对U和Pu同位素对关联的影响，" *原子核物理评论* 41(1), 2024. DOI:10.11804/NuclPhysRev.41.2023CNPC40. — PNC-CSM推转壳+八极：对关联减弱(~10-20%)；奇A核转动惯量+50-60%（泡利堵塞）。章骏、贺晓涛。
  1. R. Budaca et al., "Quadrupole-octupole geometric model for alternating parity and parity doublet bands," INPC 2025, Daejeon, Korea, May 2025. — 四极-八极几何模型；A≈224-228静态→振动临界点（类似于X(5)）。Budaca，罗马尼亚科学院Horia Hulubei国家物理与核工程研究所。
  1. Y. Yang et al., "Pear shape and tetrahedral shape competition in actinide nuclei," *Chinese Physics C* 49, 074001 (2025). — Woods-Saxon平均场：N=134/136四面体更优；²²⁴,²²⁶Ra, ²²⁸Th pear-tetrahedral共存。
  1. A.M. Khalaf et al., "Nuclear structures, moments of inertia and energy staggering for some actinide isotopes using a suggested collective model," *International Journal of Modern Physics E* 34, 2550058 (2025). — CRF3三参数转动公式；²³⁶U和²³⁸U基态带跃迁能直到24ħ偏差<3 keV——"相同带"。Khalaf，Al-Azhar大学。
  1. H.A. Alhendi et al., "Improved exponential model with pairing attenuation for backbending of moment of inertia," 2024 (revised). — 对衰减修正指数模型复现²³⁰,²³²Th, ²⁴⁰,²⁴²Pu回弯。
  1. A.D. Efimov, P.N. Koval, I.N. Izosimov, P.N. Usmanov, "Moments of inertia in IBM," *Nuclear Physics A* (2025). — IBM有效转动惯量；²²⁰Th, ²⁴²,²⁴⁴Pu回弯；高自旋超重核CDFT+IBM比较。
  1. L. Bonneau et al., "Octupole deformation in two-quasiparticle states of even–even deformed nuclei," *Physica Scripta* (2024). — Skyrme-HFB：²³⁴U基态反射对称；2qp态八极形变条件（阻塞孔态+相邻粒子态八极耦合）。
  1. P. Ring et al., "Octupole shape coexistence in ²⁸⁶No," PCM 2025, March 2025. — 三维格点CDFT：²⁸⁶No基态非轴八极形变+四面体同质异能态共存(间距0.12 MeV)。Ring，慕尼黑工业大学。