PRL里程碑突破：第一性原理终结核物理数十年“淬灭之谜”

远离稳定线的丰中子原子核的β衰变半衰期，是决定宇宙中快中子俘获过程（r-过程）核合成走向及重元素丰度分布的核心天体物理输入量。长期以来，中重质量区丰中子魔数核的微观计算受限于多体方法的多重动力学关联和手征弱电流算符的复杂性，传统理论高度依赖唯象模型的参数调校。2026年5月，发表于《物理评论快报》的标志性成果——《Ab Initio Calculations of β-Decay Half-Lives for N=50 Neutron-Rich Nuclei》，首次实现了针对N=50丰中子同量异位素β衰变半衰期的全自洽、不依赖唯象参数的第一性原理（ab initio）微观计算。

一、 引言：极端丰中子核与宇宙元素起源的纽带

宇宙中比铁更重的元素（如金、铂、铀等）约有一半是通过快中子俘获过程创生的。这一极端的天体物理过程通常发生在双中子星合并或超新星爆发等剧烈环境中。在r-过程的路径上，当物质流到达中子魔数（如N=50, 82, 126）时，由于中子分离能的骤降，核合成会在此处发生停滞，这些魔数核也因此被称为“等待点”（Waiting Points）核。

在N=50的等待点核中，⁷⁸Ni（镍-78）及其附近的极丰中子核素是连接轻重元素合成的关键纽带。这些核素的β衰变半衰期直接决定了物质流克服等待点向更重核区推进的速度。然而，由于这些极端非平衡态原子核在实验室中的合成产额极低、寿命极短，其实验测量极其困难。尽管近年来诸如日本理化学研究所（RIKEN）等前沿放射性核束装置（RIBF）取得了一系列突破，但仍有大量更靠近中子滴落线的关键核素寿命属于实验“无人区”。因此，构建具有强大预测能力的微观核结构理论显得至关重要。

二、 传统理论的困境与第一性原理的兴起

在过去数十年中，针对中重质量区原子核β衰变寿命的理论预测主要依赖于两类方法：

• 唯象模型：如准粒子随机相位近似（QRPA）或密度泛函理论（DFT）。这类方法计算效率高，但严重依赖于对已知核素数据的拟合，在外推到未知的极端丰中子区时，预测结果往往存在数个数量级的偏差。
• 传统壳模型：虽然能提供较好的谱学描述，但其有效相互作用和弱流算符同样包含大量经验经验参数。

最令人诟病的理论问题在于所谓的“轴矢量耦合常数淬灭（g_A Quenching）”。在传统唯象计算中，为了使理论计算的伽莫夫-泰勒（Gamow-Teller, GT）跃迁强度与实验值相符，研究者必须人为地将自由核子的轴矢量耦合常数g_A≈1.27乘以一个小于1的淬灭因子（通常在0.6到0.8之间）。这种缺乏微观物理基础的经验调校，极大地削弱了模型在未知核区的预测可信度。

第一性原理核结构理论的宗旨则是从现实的核子-核子相互作用出发，不引入任何针对特定原子核的微观调校，直接求解多体薛定谔方程。长期以来，由于中重核区的多体关联过于复杂，第一性原理计算一直被局限在A≤40的轻核区。如何将这一理论外推到A~80且具有强关联特征的N=50区域，是当代核物理理论的最前沿挑战。

三、 计算流水线：手征有效场论与介质中相似重整化群的融合

Zhen Li、Takayuki Miyagi与Achim Schwenk教授等人组成的国际团队，通过集成当今计算核物理的最先进技术，构建了一条高度自洽的第一性原理计算流水线：

1. 手征有效场论（Chiral EFT）的微观基石

研究的起点是基于量子色动力学（QCD）手征对称性破缺的手征有效场论。团队采用了包含次次领头阶或更高阶的两体核力（NN）与三体核力（3N）的现实核力模型。更重要的是，在描述弱流相互作用时，他们不仅使用了传统的单体电流（1B Current），还严格根据手征膨胀的自洽性，引入了两体弱电流（2B Currents / Meson-Exchange Currents）。

2. 价空间介质中相似重整化群（VS-IMSRG）

面对A~80的多体系统，全空间的第一性原理计算（如无核心壳模型）其希尔伯特空间维度会呈指数级爆炸，在计算上是不可行的。该工作采用了价空间介质中相似重整化群（VS-IMSRG）方法。该方法通过一个连续的幺正变换，系统地将高能态（核心激发）的效应重整化到选定的价空间中（例如针对 ⁷⁸Ni周边核素选择的特定价空间），从而实现多体哈密顿量与弱流算符的自洽解耦。这使得原本无法求解的巨型多体问题，转化为了可在传统壳模型框架下对角化的微观有效哈密顿量问题。

3. 跃迁流算符的精细处理

在衰变动力学方面，计算不仅严格处理了占主导地位的允许伽莫夫-泰勒（Gamow-Teller, GT）跃迁，还首次在第一性原理框架下系统评估了一阶禁戒跃迁（First-Forbidden, FF）对总半衰期的贡献。一阶禁戒跃迁在原子核远离稳定线、衰变能巨大时尤为关键，它涉及宇称的改变和更高阶的算符结构，计算难度极大。

四、 核心发现：破解“淬灭之谜”与两体电流的微观干涉

该项研究取得了两个在核物理学界引发广泛关注的核心突破：

1. 彻底解决轴矢量耦合常数的“淬灭之谜”

当研究团队在微观计算中同时包含三体核力（3N）与弱两体电流（2B Currents）时，奇迹出现了：计算出的β衰变GT跃迁强度在没有引入任何唯象淬灭因子的情况下，自发地与已有实验数据达成了高度吻合。

从微观物理机制上看，两体电流（如π介子交换电流、Δ激发态等）与单体弱流算符在多体矩阵元的求和中产生了系统的干涉效应。这一发现确凿地证明了：过去唯象模型中所谓的g_A淬灭，本质上并不是基础弱耦合常数的改变，而是由于过去的多体模型粗糙地忽略了三体核力带来的多体关联以及复杂的两体弱电流机制。 第一性原理计算成功“还原”了这一丢失的物理。

2. 两体电流的非平庸相消干涉

研究进一步揭示，在N=50的丰中子同量异位素中，两体电流的加入实际上对弱流跃迁矩阵元起到了压低作用（即增加了理论计算出的半衰期）。这种微观上的相消干涉解释了为什么某些极丰中子核的寿命比早期单纯使用单体电流预测的要长。由于第一性原理不包含外推盲区，这一机制的阐明为预测更深处的N=82以及N=126等完全不可及区域的核素寿命确立了坚实的理论范式。

五、 天体物理效应：重塑r-过程核合成路径与千新星模拟

⁷⁸Ni及其同量异位素的半衰期数据是r-过程网络计算的关键“阀门”。这篇论文提供的不依赖唯象参数的半衰期计算结果，对天体物理学界具有直接的冲击力：

• 重塑元素丰度曲线的第一峰（First Peak）与过渡区：精确的第一性原理半衰期输入，改变了数值模拟中物质流在A~80质量区的驻留时间，从而直接修正了计算出的太阳系及贫金属星中的元素丰度观测对比图，特别是精细调整了从硒（Se）、溴（Br）到氪（Kr）等元素的丰度比例。
• 千新星（Kilonova）的光变曲线预测：在双中子星合并事件（如著名的 GW170817 引力波事件）中，抛出物的放射性衰变热是驱动千新星光学/红外辐射的唯一能量来源。在中子俘获停止后的早期阶段（秒至分钟量级），N=50等待点核的衰变热贡献至关重要。该工作提供的更长的、更精确的微观半衰期，使得天体物理学家能够构建更精细的辐射转移模型，从而利用千新星的颜色和亮度演化更准确地推算合并事件中抛出的总物质质量。