核力与核结构
1935年，一个27岁的日本物理学家说：两个核子之间有一种力，由一种还没被发现的粒子传递。粒子的质量大约是电子的200倍。他把这个想法写成论文，投了出去。
两年后，那个粒子被找到了。汤川秀树拿到了1949年的诺贝尔奖。核力研究从此有了一个起点。
但起点不是终点。汤川的预言只说对了核力的长程部分——大约2飞米以外的吸引力，确实来自π介子的交换。2飞米以内，核力变得更复杂。而核力到底从哪来，为什么长这样，能不能从更基本的理论推导出来——这些问题又花了九十年，到现在还没完全解决。
这篇文章要讲的就是这段九十年的追问：从汤川的π介子到手征有效场论，从猜参数到系统性展开，从两体力到三体力，从壳模型到密度泛函，从稳定核到奇特核，从美国的FRIB到中国的HIAF。
---
核力有一个让人头疼的根本矛盾：它的母理论是量子色动力学（QCD），但QCD在核力起作用的能区——大约1 GeV以下——微扰论完全失效。就是说，你不能用费曼图一个一个算，因为每一阶的贡献都一样大，展开式不收敛。
这在之前的QCD禁闭研究里已经说过了。质子质量的91%来自夸克的结合能和胶子，不是夸克本身。质子半径约0.84飞米。而核力的力程大约1到2飞米——跟质子的大小差不多。两个质子凑到这么近，内部的夸克和胶子开始互相感受对方的存在，但又不至于完全搅成一锅粥。核力就是在这么一个尴尬的距离上运作的。
所以从QCD直接算核力，相当于要你在不打开盒子的情况下描述两个盒子里的东西怎么互相推拉。格点QCD可以做一些计算——把时空切成网格，在超级计算机上数值求解——但目前只能在比物理π介子质量更重的条件下算两核子散射，离日常使用的核力还差得远。
那怎么办？退一步。用有效理论。
---
有效理论的思想很简单：如果你只关心低能现象，就不需要知道高能细节。就像你计算苹果落地不需要知道苹果的原子结构，你算核力也不需要知道夸克在10飞米以下干什么。
Weinberg在1990和1991年做了这件事。他注意到QCD有一个近似的对称性——手征对称性——在低能区被自发打破了。打破的后果是产生了π介子，作为这种打破的"残余振动"——物理学家叫它Goldstone玻色子。π介子很轻（约140 MeV），正好是核力长程部分的传递者。
Weinberg说：既然π介子是QCD对称性破缺的产物，那我就可以用对称性来约束核力的形式。具体来说，把核力按（Q/Λ_χ）的幂次展开，Q是典型动量（约π介子质量），Λ_χ ≈ 1 GeV是手征破缺能标。低阶的项由对称性完全决定，高阶的项需要拟合实验数据，但贡献越来越小。
这个纲领把核力的推导从"猜参数"变成了"系统性展开"。你不再需要猜核力长什么样——对称性告诉你它必须长什么样，你只需要确定少数几个低能常数。
手征有效场论（χEFT）的展开阶次是这样排的：LO（领头阶）只有一个π交换图加两个接触相互作用，6个参数。NLO多了双π交换和更多接触项。N2LO第一次出现三体力——这是关键转折。N3LO是当前最常用的精度，Entem和Machleidt在2003年构建了N3LO NN势，24个参数。2024年，他们在Progress in Particle and Nuclear Physics上发表了最新综述，报告了N4LO核力的进展——参数更多，但系统性更好，与NN散射数据的拟合χ²/datum接近1。
但手征EFT有一个根本性的困扰：收敛性问题。每一阶的修正应该是上一阶的（Q/Λ_χ）²倍——大约1/3到1/4。但实际计算中发现，某些可观测量在N3LO和N4LO之间的改进不够大，让人怀疑展开是否真的在快速收敛。2026年1月，Epelbaum和Machleidt在Springer出版了新的专著章节，系统梳理了χEFT的当前状态，承认收敛性问题仍然存在但并非致命——只是需要更仔细地处理截止参数的依赖性。
---
核力还有一个让人不安的事实：两体力不够用。
如果你用最好的两体核力去算核物质——无限大的均匀核子系统——你会得到一个结果：每个核子的结合能大约是16 MeV，平衡密度大约是0.16 fm⁻³。但两体力给出的结果偏松——结合不够紧，密度偏低。这叫"核物质饱和性的欠束缚问题"。
三体力修好了这个缺口。最经典的三体力图像来自Fujita和Miyazawa在1957年的工作：两个核子交换一个π介子，其中一个核子在中间被激发成Δ共振态，然后又退激发回来。这个"2π交换+中间Δ"的过程自然产生了三体力。
在手征EFT里，三体力在N2LO才出现，包含三种类型：2π交换图（叫Fujita-Miyazawa型的推广）、1π交换-接触图、纯接触图。三种加在一起，才能正确描述核物质的饱和点和轻核的结合能。
2026年3月，一篇发表在MDPI上的论文系统研究了三体力对轻到中等质量晕核的影响，发现三体力对晕核半径和滴线位置有决定性作用——没有三体力，你甚至无法正确预言哪个核是最后一个稳定核。
三体力不只是理论上的修修补补。它关系到一些根本性的问题：中子星的状态方程——中子星内部有多密实，取决于三体力怎么压——以及r过程核合成路径中不稳定核的结合能。这些是核天体物理最核心的输入量。
---
知道了核力，下一步是算核结构。这里有三条主要路线，各有各的地盘。
壳模型是最老的。Mayer和Jensen在1949年发现，原子核里的核子跟原子里的电子一样，也按能级一层一层填充——2、8、20、28、50、82、126这些"幻数"对应特别稳定的核。1963年他们拿了诺贝尔奖。壳模型的核心假设是：每个核子在一个平均场里独立运动，剩余相互作用当微扰处理。
壳模型算得准，但算不了大的核。原因是计算量随核子数指数增长——你追踪的组态数在中等质量核（比如锡-132）里就已经超过10^15，再大的核就是天文数字。2025年，东京大学的MCSM（蒙特卡洛壳模型）团队在钙同位素链上做了前所未有的计算，但也只能到A≈60左右。
核密度泛函理论（DFT）是另一条路。它不追踪每个核子的量子态，而是直接算原子核的能量——作为核子密度分布的函数。想法跟凝聚态物理里的DFT一样，但核DFT的泛函形式不是从基本原理推导的，而是参数化的。主流的泛函有三种：Skyrme力（零程、简单、快）、Gogny力（有限程、稍慢但更好）、相对论平均场（从量子场论出发、自旋-轨道劈裂自然出现）。
核DFT能算整个核素图上的所有核——从最轻到最重、从滴线到稳定线——但精度有限。结合能的偏差大约1-2 MeV，对超重核和滴线核的外推更不可靠。2025年2月，一篇发表在INSPIRE-HEP上的论文尝试把从头算的结果注入核DFT的泛函参数——用精确的小核计算来约束泛函形式，再用约束好的泛函算大核。这是两条路线的桥梁。2025年8月，Nature Communications上发表了一篇用机器学习改进无轨道密度泛函的文章——神经网络学会了处理形变核的量子效应，解决了这个方向上一个长期存在的教科书难题。
最野心勃勃的路线是从头算（ab initio）。它直接从手征EFT核力出发，不做任何唯象拟合，用数值方法求解多体薛定谔方程。方法有好几种：量子蒙特卡洛（QMC）、耦合簇（CC）、介质中相似重整化群（IM-SRG）、无核心壳模型（NCSM）。它们的共同点是：从核力出发，不引入额外参数。
从头算的边界在快速推进。2024年4月，一篇发表在Physics Letters B上的Letter报告了首次对开壳钙、镍、锡同位素链的从头算Bogoliubov耦合簇计算。这意味着从头算已经从中等质量核（A≈40-48的钙）推进到了重核（A≈132的锡）。这在五年前被认为是不可能的。2025年3月，上海科技大学和中科院物理所合作报告了利用神经网络量子蒙特卡洛方法求解核多体问题的进展——FeynmanNet架构把领头阶无π有效场论核力的多体计算推到了新精度。
但从头算仍然有它的天花板。锡-132有132个核子，希尔伯特空间的维度仍然是天文数字。耦合簇方法需要截断到CCSD或CCSDT级别，截断带来的系统误差很难估计。IM-SRG方法对3体力以上贡献的处理仍然是近似的。
三条路线正在靠近。壳模型用从头算的核力做有效相互作用，DFT用从头算的结果约束泛函，从头算用DFT的结果做初始猜测。这不是三条独立的小路，而是同一座山的三个坡面，正在汇合。
---
理论在进步，实验也没闲着。事实上，实验正在以前所未有的速度产出新数据。
2022年5月，美国密歇根州立大学的FRIB正式开始运行。它是目前世界上最强的放射性核束装置，400 MeV/u的重离子直线加速器，理论上能产生已知约5000种同位素的80%。2024年2月，FRIB宣布发现了5个前所未见的重元素同位素：铥-182、铥-183、镱-186、镱-187、镥-190。这些同位素在地球上不存在——它们只在超新星爆发或中子星碰撞的极端条件下短暂出现。2026年3月，日本RIKEN又发现了8个新同位素，包括铬-72、锰-74、铁-77等，进一步拓展了核素图的边界。
FRIB对奇特核的研究尤其重要。奇特核是远离稳定线的不稳定核素，它们的行为经常违背从稳定核外推的预期。
最戏剧性的例子是晕核。1985年，Tanihata等人在伯克利实验室发现锂-11的截面异常大——远超核物质密度外推的预期。后来才明白：锂-11由一个锂-9的核心和两个远离核心的中子组成，这两个中子的轨道半径是核心的几倍，形成一圈"晕"。锂-11看起来像一个房子大的棉花糖——中心一小团硬芯，外面一大圈蓬松的云。
晕核改变了核物理学家对原子核的基本认识。传统图像里，原子核是一个密实的球，所有核子挤在一起。晕核告诉你：核的边缘可以是模糊的、松散的，最后一个核子甚至可以游离到核心之外很远。
另一个让人意外的现象是新幻数的出现。传统的幻数（2、8、20、28、50、82、126）在稳定核里坚如磐石，但到了滴线附近，有些幻数消失了，新的幻数出现了。比如钙-54（N=34）被提议为新幻数的候选——在丰中子区，N=34的能隙可能重新打开。2025年，多个实验组在RIKEN和FRIB上对钙-54进行了谱学测量，结果正在分析中。如果确认，这就意味着核"壳层"不是一成不变的——它会随同位素链条上的位置而演化。
更根本的问题是：滴线在哪里？滴线是核素图上最后一个束缚核的边界——超过这条线，多加一个中子或质子，核就不再束缚。实验上确定滴线非常困难，因为你需要产生并探测极其短命的核素。氧的滴线是当前最热的争议之一：传统预期N=16（氧-24），但有些实验和计算暗示可能是N=28（氧-44）。这个差别意味着8个中子的不确定性——对核结构理论来说，这是生死攸关的问题。
---
中国在这场竞赛中有一张重要的牌：HIAF。
2025年10月28日，中国科学院近代物理研究所宣布，强流重离子加速器装置（HIAF）在广东惠州完成了全线束流贯通调试。这是目前国际规模最大的全离子加速器。2026年4月9日，HIAF的高精度环形谱仪实验终端通过了工艺测试，完成了首个终端实验。HIAF的束流强度和能量精度都将达到世界领先水平，能够精确测量远离稳定线的核素质量和寿命。
与HIAF同址建设的还有CiADS（加速器驱动嬗变系统），将是世界上首个兆瓦级加速器驱动次临界系统。虽然CiADS的主要目标是核废料嬗变，但它的质子加速器也能产生放射性核束，为核结构研究提供补充。
中国的核结构理论也在发力。北京大学核物理与核技术全国重点实验室在协变密度泛函理论方面有深厚积累，中山大学核理论组在连续谱与形变的耦合方面做了原创性工作。中科院近物所在核质量测量方面是世界领先的——他们在储存环上测量的核质量精度达到了10 keV级别。
---
碳-12的霍伊尔态是核天体物理里最精巧的一个故事。
恒星内部的碳是怎么来的？三个氦-4核（α粒子）碰巧凑到一起，形成碳-12。但这个过程在标准条件下几乎不可能——两个α粒子先合成铍-8，但铍-8的寿命只有10^-16秒，第三个α粒子必须在这段时间内撞上去。概率极低。
1954年，Fred Hoyle指出，碳-12必须有一个特殊的激发态，能量恰好位于三个α粒子的总能量附近，这样才能共振增强三α反应率。他预言这个态在7.68 MeV。实验一找，果然在7.65 MeV——这就是霍伊尔态。如果这个态不存在，恒星内部几乎产不出碳，也就没有碳基生命。你正在读这篇文章，得感谢这个态。
2025年11月，一篇发表在Physics Letters B上的论文报告了霍伊尔态γ衰变概率的新测量。他们用了一种新方法——三重符合探测散射α粒子、反冲碳-12核和γ射线——得到了比以前更精确的辐射衰变分支比。这个数字直接决定三α过程在恒星内部的反应率，进而影响大质量恒星的演化和超新星爆发的条件。
2023年，Nature Communications上发表了一篇用从头算框架研究碳-12核态几何的工作，给出了霍伊尔态的模型无关密度图。结果显示，霍伊尔态的几何确实倾向于三个α粒子形成松散的等边三角形——跟Hoyle当年的直觉一致，但现在有了第一性原理的支撑。
---
回到核力本身。我们走到哪了？
手征EFT给出了一个系统性框架，从QCD的对称性出发推导核力。N4LO的NN势已经达到了接近唯象势的精度，三体力也已经被纳入。但收敛性问题还没解决，截止参数的依赖性仍然是一个方法论上的裂缝。
从头算把核力塞进多体方程里直接求解，已经推进到了锡-132。但计算量仍然限制了精度和范围，三体力以上贡献的处理仍不完善。
核DFT能覆盖整个核素图，精度约1-2 MeV，外推到超重核和滴线仍然不可靠。机器学习正在被引入来改善泛函，但这还处于早期。
奇特核不断带来惊喜——晕核、新幻数、滴线的移动——每一次都在挑战从稳定核外推的预期。FRIB和HIAF将在未来十年产出海量的新数据。
核力和核结构的研究本质上是在回答一个简单的问题：物质为什么长这样？为什么碳-12有一个恰好让生命存在的激发态？为什么核素图在滴线处突然终止？为什么中子星可以那么密实而不坍缩？每一个答案都牵扯到核力的具体形式——它怎么吸引、怎么排斥、怎么在三个核子之间产生额外的力。
90年前汤川给出的只是核力的一个轮廓。现在我们有了更精细的画笔，也有了更清晰的画布。但画面还远未完成。
---
关键时间线
| 时间 | 事件 | 来源 |
|---|---|---|
| 1935 | 汤川秀树预言π介子传递核力 | Proc. Phys. Math. Soc. Japan |
| 1947 | π介子在宇宙线中被发现 | Bristol团队 |
| 1949 | Mayer和Jensen提出核壳层模型（1963诺贝尔奖） | PRL |
| 1953 | Bohr和Mottelson提出集体模型（1975诺贝尔奖） | Phys. Rev. |
| 1957 | Fujita-Miyazawa三体力 | Prog. Theor. Phys. |
| 1985 | Tanihata等发现晕核Li-11 | PRL |
| 1990-91 | Weinberg提出手征EFT纲领 | PLB 251, NPB 363 |
| 2003 | Entem-Machleidt N3LO核势 | PRC 68 |
| 2018 | CEPC概念设计报告 | IHEP |
| 2022.05 | FRIB开始运行 | MSU |
| 2024.02 | FRIB发现5个新同位素 | APS Physics |
| 2024.04 | 首次从头算开壳Ca/Ni/Sn链（BCC） | PLB |
| 2024.05 | 手征EFT核力最新综述（N4LO进展） | Prog. Part. Nucl. Phys. |
| 2025.02 | 核DFT：从头算约束泛函 | INSPIRE-HEP |
| 2025.08 | 机器学习无轨道DFT解决形变核难题 | Nature Commun. |
| 2025.10 | HIAF全线束流贯通调试 | 中科院近物所 |
| 2025.11 | 霍伊尔态γ衰变概率新测量 | PLB |
| 2026.01 | Epelbaum-Machleidt手征EFT新专著 | Springer |
| 2026.03 | RIKEN发现8个新同位素 | FRIB最新同位素列表 |
| 2026.03 | 三体力对晕核影响系统研究 | MDPI |
| 2026.04 | HIAF高精度环形谱仪首终端实验 | 中科院近物所 |