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

壳模型与组态相互作用

2026-05-30

1949 年，Maria Goeppert Mayer 在芝加哥大学的办公室里盯着核素图看了几个月。她注意到一个模式：当中子数或质子数恰好是 2、8、20、28、50、82、126 的时候，原子核就特别稳定、特别难被激发。她说，这不就像原子的电子壳层吗。同一年，Hans Jensen 在海德堡独立发现了同一件事。核壳模型就这么诞生了。七十七年后，壳模型和它最核心的计算方法——组态相互作用（configuration interaction，CI）——仍然是人类理解原子核结构的最主要工具。这篇文章要讲的是：这个工具在过去两年里经历了什么，它的极限在哪，以及新冒出来的几个方向会把我们带到哪里。

壳模型的操作说起来其实很朴素。先找一个"芯"——选一个填满了最低壳层的核，通常是个双幻核。芯外面的轨道叫价空间。你允许价空间里的质子和中子在这些轨道上任意分配——每一个分配方式叫一个组态。真正的核波函数，是所有符合对称性要求（总角动量 J、总同位旋 T、宇称 π）的组态的线性叠加。所以叫"组态相互作用"——不同组态之间通过剩余相互作用混合，这才是真实核的样子，而不是单一组态。

问题是组态的数量涨得比任何直觉都快。sd 壳（从氧-16 到钙-40）的 M-scheme 维度大概几百万到几亿，还勉强能对角化。pf 壳（钙-40 到镍-78）轻轻松松就上百亿。再往外走，蛮力对角化就是不可能的任务——不是难，是不存在够大的计算机。

所以壳模型的七十七年历史，本质上是在一个问题上不断做减法：怎么把指数爆炸的组态空间压到一个可计算的子空间，同时保证物理不被压丢。

先说最近两个把压缩做到极致的新工具。

Calvin W. Johnson 和 Oliver C. Gorton（圣地亚哥州立大学）在 2025 年 3 月扔出了一篇预印本，标题很平淡："A Weak Entanglement Approximation for Nuclear Structure"。里面的内容一点不平淡。他们发现：当质子数和中子数不相等的时候（N ≠ Z），原子核 CI 波函数的质子部分和中子部分之间的量子纠缠非常弱。这意味着什么？意味着你可以把质子空间和中子空间分开处理，各自对角化，只取最低能的一批本征态，再乘回来——完全不需要在联合希尔伯特空间里做对角化。

他们把这个近似叫做 PANASh（Proton And Neutron Approximate Shell-model）。验证结果令人难以相信：在 pf 壳里，PANASh 在一台笔记本电脑上算出了和全 CI 几乎一样的结果（能谱差几百 keV 以内）。需要强调的是——在 pf 壳，全 CI 的维度可达数十亿。但真正吓人的是他们拿 PANASh 试了铯-129（50 亿维）和铈-132（2.4 万亿维）。铈-132 的全 CI 维度是 2.4 × 10¹²——远超当前任何超算的处理能力——而 PANASh 给出了可验证的、物理上合理的结果。这是壳模型历史上第一次触达万亿维度的 CI 空间。

如果 PANASh 走的是"分治"路线，Johnson 的另一个 2025 年工具 TRACER（与 Keller 合作，11 月预印本）走的是"预筛选"路线。TRACER 快速计算哈密顿量在各类粒子占据数构型空间中的能量中心值（centroid），以此为基础在 BIGSTICK 对角化之前优化截断权重。计算快到可以在参数空间里做迭代优化——这在以前是不可能的。

BIGSTICK 本身是 Johnson 组开发的旗舰 CI 代码，用即时分解（on-the-fly factorization）避免存储整矩阵。2025 年 11 月版本 8.0 发布，同时上线了 RHODIUM 后处理器——专门从 BIGSTICK 波函数中提取密度矩阵、谱因子和贝塔衰变振幅。

太平洋另一边，清水则孝（Noritaka Shimizu，东京大学核子研究中心）的 KSHELL 代码走了类似的即时生成路子，加上宇都野穣（Yutaka Utsuno，JAEA）和角田佑介（Yusuke Tsunoda）2019 年开发的厚重启块 Lanczos 算法。KSHELL 的标志性战绩：镍-56 在 pf 壳的全 CI 计算——10.87 亿维矩阵，8192 个 CPU 核，145 秒。2025 年 Patel、Srivastava、Shimizu 和 Utsuno 用这套工具劈开了铟-99 到铟-131 整整 33 个奇 A 核——计算和实验符合极好，从能级、电磁矩到谱因子一次性全部拿下。

但这些代码都必须选择一个有限价空间。没有一种有效相互作用是普适的——它在某个壳层里管用，出去就不灵了。过去五年，一个全新的范式在成形：不要唯象地调参，从手征有效场论（chiral EFT）的核子-核子相互作用出发，通过重整化群流动直接导出价空间的等效哈密顿量。一个自由参数都没有。

这个方向的核心引擎是 VS-IMSRG（价空间介质相似重整化群）。2024 年，Alexander Tichai（当时在达姆施塔特工业大学）、Kapás、Miyagi 等人实现了 VS-IMSRG 和 DMRG（密度矩阵重整化群）的结合——叫 VS-DMRG。用矩阵乘积态（MPS）描述波函数，而不是 CI 的线性展开。结果：同样的精度下，希尔伯特空间需求减了两个数量级。对 N=50 同中子素的能谱和 E2 跃迁做了系统计算，发表在 Physics Letters B 上。

但一个令人不安的真相也浮出了水面。Stroberg、Henderson、Hackman、Hagen 和 Holt 在 2022 年（PRC 105, 034333）就系统确认：VS-IMSRG 在 sd 壳里系统地低估了 E2 跃迁强度，而且这个低估幅度与同位旋基本无关。问题不出在相互作用——出在截断。VS-IMSRG 丢掉的多粒子-多空穴激发对 E2 有系统性贡献。这个"缺失的 E2 强度"问题是当前从头算壳模型最棘手的开放性难题。

说了这么多计算方法，壳模型到底算出了什么？

2025 年 1 月 Nature Physics 上，Alexandra Gade（密歇根州立大学/FRIB）领导的合作组用束流伽马谱学测了铬-62 的 0₂⁺ 态。铬-62 就在 N=40 反转岛的心脏边上。壳模型告诉他们：基态是 4 粒子-4 空穴闯入组态主导的，0₂⁺ 来自 2 粒子-2 空穴组态——这是教科书级的形状共存。同一个月，ISOLDE（CERN）的 CRIS 实验重新测了铬-61 的基态自旋：不是之前大家以为的 5/2⁻，而是 1/2⁻。LNPS 壳模型计算（Lenzi-Nowacki-Poves-Sieja，其相互作用 2010 年首次发表在 PRC 82, 054301）完美预言了这个自旋，并指出铬-61 是三相形变——过渡核，刚好站在反转岛的"西岸登陆点"。

2025 年 12 月，Ceulemans 等（PRL 135, 252502）在 ISOLDE 用 (d,p) 反应找到了镍-69 中 νd₅/₂ 轨道在 2.56 MeV 的主导碎片。这个 d₅/₂ 轨道就是驱动 N=40 反转岛形变的关键角色——它和 g₉/₂ 之间的能隙在大质量富中子区变得极窄，一丁点剩余相互作用就足以把核子推进闯入组态、引发形变。LNPS 的有效性在这一带反复被验证——Nowacki、Obertelli 和 Poves 2021 年在 PPNP 上的长篇综述将 N=40 反转岛描述为"壳模型预言能力的最佳展示之一"。

再往重核区走。2024 年 9 月，Dao 和 Nowacki（斯特拉斯堡大学，IPHC）在 arXiv 上投了一篇标志性论文：用离散非正交壳模型加变分后投影（DNO-SM(VAP)），首次对超重核锘-254 做了完整的壳模型描述。锘-254 的质子数 102，中子数 152——这是超重核区的门户。Dao 和 Nowacki 用了多达 32 个非正交 Slater 行列式，完美复现了基态转动带、K=8⁻ 长寿命同核异能态、新测得的 0₂⁺ 态和一个 (K=4) 4⁺ 态。壳模型正式进入了超重核区。他们的相互作用底子是 Kuo-Herling 势——从 Hamada-Johnston 核子-核子势导出，加上唯象的单极修正。

超重核区还有一个更戏剧性的理论进展来自大塚孝治（Takaharu Otsuka，东京大学/RIKEN）的蒙特卡罗壳模型（MCSM）。MCSM 的本质是用随机抽样加变分优化避免全维度对角化。加上其变体准粒子真空壳模型（QVSM），这套方法能处理 10³³ 维的希尔伯特空间——是传统 CI 极限的几十个数量级以上。2025 年，大塚团队（角田佑介、清水则孝、Duguet、宇都野穣、阿部喬）在 EPJA 上得出了一个颠覆性结论：在重变形核（A > 140）中，三轴形变不是例外，是常态。铒-166 的"伽马带"——教科书教了几十年的伽马振动——经 MCSM 确认是三轴转动。真正的伽马振动态在更高的能区（约 3 MeV 以上），带头是 0₃⁺。连一直被视为轴对称典范的钐-154，MCSM 也给出了有限的三轴度。大塚团队进一步提出了"PT-plot"（相位指定 T 图），可视化波函数在三轴空间的节面结构——这对识别振动态和三轴转动态是质的飞跃。

壳模型给出的不止是谱学。它的核矩阵元直接进入中微子物理的核心方程。2025 年，Castillo、Jokiniemi 和 Menéndez（巴塞罗那大学）在 PLB 860, 139181 上完成了对 ⁷⁶Ge、¹⁰⁰Mo 和 ¹³⁶Xe 无中微子双贝塔衰变矩阵元的 N²LO（手征展开的次次领头阶）完整计算——同时用了壳模型和 pnQRPA 两种独立方法。结果揭示了一个严重的理论问题：N²LO 修正对壳模型是 -(5-10)%，对 QRPA 却是 +(10-15)%，连符号都不一致。Kota 和 Sahu 同一年的统计壳模型计算（arXiv:2507.02707）给出的矩阵元更是比传统壳模型结果小一倍左右。这个理论分歧的大小，直接决定了吨级无中微子双贝塔衰变实验（如 LEGEND、nEXO）的灵敏度标定能否可靠——目前不能。

连续谱耦合是传统壳模型最大的盲区——它用的是谐振子基，天然束缚，无法描述粒子发射。Gamow 壳模型（GSM）在这个盲区上打了二十多年的楔子。Nicolas Michel（中科院近物所/兰州）、Witek Nazarewicz（MSU/FRIB）和 Marek Płoszajczak（GANIL）是 GSM 的奠基人和主要推手。GSM 用复合能量的 Berggren 基统一处理束缚态、共振态和散射连续谱。2024-2025 年 GSM 的进展集中在两个方向：一是从头算 GSM（NCGSM）直接算到了四中子体系，与 2022 年 Nature 上的四中子共振实验结果接近；二是 2025 年徐志诚（复旦大学）等人在 GSM 框架内实现了从头算价空间有效算符，成功复现了滴线核 ²²Si 与其镜像核 ²²O 之间的贝塔衰变同位旋不对称——其中 s₁/₂ 连续谱的贡献是关键。

机器学习正在从另外一个方向渗透进来。Roth 组（达姆施塔特工业大学）在 2024-2025 年连续甩出几篇高影响力文章。Wolfgruber、Knöll 和 Roth 在 PRC 110, 014327 (2024) 里训练了一个 ANN——让它从轻核（²H、³H、⁴He）的 NCSM 计算结果中"学会"截断空间如何向无穷空间外推，然后把学到的模式直接用到更重的 p 壳核（⁶Li、⁹Be、¹²C、¹⁶O）——不经过任何重新训练。这个"训一次、到处用"的性质极为重要，因为它绕开了 NCSM 的维度爆炸——外推误差被量化了出来，不再是一个含糊的"系统误差"。

2025 年 Knöll 等（PRC 111, 064304）把框架推广到锂同位素链，Knöll、Agel 等（PRC 112, 034306）进一步训练了"转码器"神经网络——利用能级、半径、电磁矩之间的内在关联，从一个已收敛的观测量预测另一个没收敛的观测量。2025 年 Wolfgruber、Gesser 等（PRC 112, 014306）用同样的方法算了硼同位素链的电荷半径——提示硼-8 可能存在质子晕。

机器学习在壳模型中的另一个应用方向是直接预测核数据。2024 年 NST 上的一篇论文用 LightGBM（轻量梯度提升机）预测了 642 个核的 2⁺ 激发能，平均偏差比壳模型小 18 倍，比贝叶斯神经网络好 30%。2025 年 Communications Physics 发表了 ML 轨道自由 DFT 的工作——用核岭回归在无需单粒子轨道信息的情况下捕获了核壳效应和形变，这是 DFT 历史上第一次做到这一点。

中国贡献方面。CI-RHF（组态相互作用相对论 Hartree-Fock）模型是近几年一个重要的理论发展。它用一个统一的密度依赖拉氏量导出不同核区的有效哈密顿量，不需要跨壳调参。2025 年在 CPC 正式发表的结果显示：在 Ne 同位素链上，从 sd 壳的 ¹⁸⁻²⁸Ne 一直算到跨壳的 ³⁰Ne 和 ³²Ne（其中 ³⁰Ne 的极低 2⁺ 能量标志了 N=20 幻数在丰中子区的消失），全部和实验对得上。核心机制是芯极化效应被自洽地包含在 CI 基矢的构造中。

中科院近物所（兰州）的陈宁、李建国等在 2024 年 PRC 上针对 ⁷⁸Ni 以外的核区做了大规模壳模型计算（LSSM），特别是 N=51 和 N=52 同中子素的单粒子演化和集体性。对 ⁸⁴Ge 的 B(E2) 预测是目前所有理论中与实验趋势最接近的。同一时期，Guo Li 等在中国科学院预印本平台（ChinaXiv）提交了将 FCIQMC 引入核结构的首次系统研究——FCIQMC 从量子化学借来的核心优势是：它能处理强关联体系，而这正是传统 CI 的致命弱点。在 sdpf 壳 Mg 同位素上的测试证明了这个方向的可行性。

还有两条不能遗漏的新线索。一是量子信息对壳模型的诊断。Savage 等 2025 年 3 月在 PRC 111, 034317 上把 BIGSTICK 生成的 p 壳（12 量子比特）和 sd 壳（24 量子比特）波函数映射到量子比特寄存器上，计算了多体纠缠度（n-tangle）和"魔法"（magic，非稳定子度）。核心发现：sd 壳波函数有显著的量子纠缠——主要由质子-中子组态混合驱动——也有显著的量子魔法（超过稳定子态的模拟能力）。这有两层含义：第一，这些核在经典计算机上的高效模拟确实非常困难；第二，一旦量子计算机的规模和纠错能力达标，这些核反而是天然的量子模拟对象。

二是密度矩阵重整化群（DMRG）的入侵。Tichai 等 2024 年那篇 VS-DMRG 文章（PLB 855, 138841）证明了矩阵乘积态的 ansatz 在核结构中不仅可行，而且在压缩效率上远远优于 CI 的 Slater 行列式展开——同样的能谱精度只需不到百分之一的基矢数量。这暗示了一个方向：未来大规模壳模型计算可能不再以 Slater 行列式为基矢，而是以矩阵乘积态为基矢。

回到最根本的问题：壳模型到底能走多远。乐观的回答是上面写的一切——PANASh 把可以触及的维度推高了三四个数量级，MCSM 推高了二十个数量级，ML 外推绕开了维度瓶颈，从头算 VS-IMSRG 正在消灭自由参数，DMRG 提供了一种比 CI 高效百倍的压缩方式。所有这些都在同时向前走。

诚实的回答是这样的：我们至今不能从第一性原理出发，对任何一个中重核的全希尔伯特空间做精确对角化。今天每一个"壳模型成功"都依赖某种数学上不严格但物理上合理的截断。质子-中子弱纠缠假设在 PANASh 里表现神勇——但它对 N=Z 核区不成立，因为那里的质子壳和中子壳完全对称，纠缠天然大。LNPS 在 N=40 反转岛一统天下——但往 N=50 方向每推一步都需要"改一改"单粒子能。从头算 VS-IMSRG 已经基本消灭了价空间内的自由参数——却在 E2 跃迁强度上困了整整三年，因为那个截断把关键的多粒子-多空穴关联丢了。

这些不是挫折。它们是同一个深层次物理事实的不同侧面：原子核恰好卡在一个最尴尬的尺度上。它不够小，不能像轻核那样用量子力学的"蛮力"算到底。又不够大，不能用统计力学——涨落太大，平均不了一切。它是中等尺度的量子多体系统，有丰富的关联，有涌现的集体行为，有形状共存和量子相变。这正是最难处理的那一类系统——也是壳模型在这个尺度上最有价值的原因。

也许最重要的一个转变正在发生——壳模型的终极目标已经不是"在对角化意义上解决问题"，而是"在物理意义上建立理解"。上面所有这些进展——纠缠度量、魔法度量、转码器神经网络、PT-plot 可视化——它们的共同特征是：让计算机不再只是算出一个数字，而是帮你看到波函数的结构。看到质子怎么跟中子耦合，看到组态怎么混合，看到三轴度怎么从张量力的一个特殊分量中长出来。

1949 年 Mayer 和 Jensen 看到的是壳。2025 年我们看到壳里那些壳之间的东西。

参考文献

  1. M.G. Mayer, Phys. Rev. 75, 1969 (1949); O. Haxel, J.H.D. Jensen, H.E. Suess, Phys. Rev. 75, 1766 (1949) — 核壳模型的独立发现
  2. C.W. Johnson, O.C. Gorton et al., arXiv:2504.00273 (2025) — PANASh: 质子-中子弱纠缠近似壳模型
  3. C.W. Johnson, K.S. Keller, arXiv:2511.03161 (2025) — TRACER: 能量中心值与截断优化
  4. C.W. Johnson, W.E. Ormand, K.S. McElvain, arXiv:1801.08432 (updated 2025) — BIGSTICK v8.0 CI 壳模型代码
  5. RHODIUM post-processor, arXiv:2510.23787 (2025) — BIGSTICK 波函数分析工具
  6. C.W. Johnson, W.E. Ormand, P.G. Krastev, arXiv:1303.0905 (2013) — BIGSTICK 即时分解算法基础
  7. N. Shimizu, arXiv:1310.5431 — KSHELL: 大规模并行壳模型代码
  8. N. Shimizu, T. Mizusaki, Y. Utsuno, Y. Tsunoda, arXiv:1902.02064 (2019) — 厚重启块 Lanczos 方法
  9. Patel, Srivastava, N. Shimizu, Y. Utsuno et al. (2025) — 铟-99 至铟-131 大规模壳模型计算
  10. T. Otsuka, Y. Tsunoda, N. Shimizu, Y. Utsuno et al., Eur. Phys. J. A 61, 126 (2025) — MCSM/QVSM 重变形核三轴性系统研究
  11. Y. Tsunoda, T. Otsuka, N. Shimizu, T. Duguet, Y. Utsuno, T. Abe, arXiv:2507.20275 (2025) — 重变形核的振动模式
  12. A. Tichai, R. Kapás, T. Miyagi, V. Werner, Ö. Legeza, A. Schwenk, G. Zarand, Phys. Lett. B 855, 138841 (2024) — VS-DMRG: N=50 同中子素谱学
  13. S.R. Stroberg, J. Henderson, G. Hackman, P. Ruotsalainen, G. Hagen, J.D. Holt, Phys. Rev. C 105, 034333 (2022) — VS-IMSRG E2 跃迁强度系统性低估
  14. S.M. Lenzi, F. Nowacki, A. Poves, K. Sieja, Phys. Rev. C 82, 054301 (2010) — LNPS 有效相互作用
  15. F. Nowacki, A. Obertelli, A. Poves, Prog. Part. Nucl. Phys. 120, 103866 (2021) — N=40 反转岛全面综述
  16. D.D. Dao, F. Nowacki, arXiv:2409.08210 (2024) — DNO-SM(VAP): ²⁵⁴No 首个完整壳模型描述
  17. A. Gade et al., Nat. Phys. 21, 37 (2025) — ⁶²Cr 形状共存
  18. CRIS/ISOLDE Collaboration, Phys. Rev. C (2025, Editor's Suggestion) — ⁶¹Cr 基态自旋重测定
  19. Ceulemans et al., Phys. Rev. Lett. 135, 252502 (2025) — ⁶⁹Ni 中 νd₅/₂ 轨道的直接证据
  20. Castillo, Jokiniemi, Menéndez, Phys. Lett. B 860, 139181 (2025) — 无中微子双贝塔衰变 N²LO 完整矩阵元
  21. V.K.B. Kota, R. Sahu, arXiv:2507.02707 (2025) — 统计壳模型 0νββ 矩阵元
  22. M.J. Savage et al., Phys. Rev. C 111, 034317 (2025) — 原子核中的量子魔法与多体纠缠
  23. T. Wolfgruber, M. Knöll, R. Roth, Phys. Rev. C 110, 014327 (2024) — ANN 增强 NCSM 外推
  24. M. Knöll et al., Phys. Rev. C 111, 064304 (2025) — Li 同位素 ANN 外推
  25. M. Knöll, Agel, Wolfgruber, Maris, Roth, Phys. Rev. C 112, 034306 (2025) — ML 关联外推电磁观测量
  26. T. Wolfgruber, Gesser, M. Knöll, Maris, Roth, Phys. Rev. C 112, 014306 (2025) — 硼同位素半径高精度从头算
  27. 陈宁, 李建国 et al., Phys. Rev. C (2024) — ⁷⁸Ni 以上核区大规模壳模型计算
  28. Guo Li et al., ChinaXiv (2025); Nucl. Sci. Tech. 36, 212 (2025) — FCIQMC 引入核结构计算
  29. CI-RHF Collaboration, Chin. Phys. C 49, 064104 (2025) — 组态相互作用相对论 Hartree-Fock 模型
  30. N. Michel, W. Nazarewicz, M. Płoszajczak et al. — Gamow 壳模型系列工作（2024-2025 进展报告及 GSM 教程）
  31. N. Michel, M. Płoszajczak, "Gamow Shell Model", Lecture Notes in Physics 972, Springer (2021) — GSM 标准教科书
  32. 徐志诚 et al., INPC 2025 报告 — 从头算 GSM 有效算符与连续谱 β 衰变
  33. K.D. Launey, T. Dytrych et al., arXiv:2510.15171 (2025) — 对称性适配从头算核反应综述
  34. Sharypov, Mazur, Shirokov, arXiv:2504.06594 (2025) — ML 外推 NCSM 至无穷基
  35. NST (2024) — LightGBM 预测 642 核 2⁺ 激发能
  36. B.A. Brown, W.D.M. Rae, Nucl. Data Sheets 120, 115 (2014) — NuShellX@MSU 壳模型代码
  37. 周先荣 et al. (华东师范大学), 2025 — Gamow 壳模型拓展至超核体系