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

gA淬灭根源判别-Delta与两体流与100Sn

2026-09-21

# g_A 淬灭的根源判别:Δ-isobar、两体流,和一块叫 ¹⁰⁰Sn 的石头

自由中子里那个数字是 1.276。把它塞进原子核,同一个数字就缩到 0.74 到 0.77 之间[30]。缩掉的这两三成,核物理界追了五十年。

追法有两种。一种是把嫌疑犯列出来,一个个审——这是综述的干法。另一种是找一个能一锤定音的案子,逼所有嫌疑犯同时出庭对质——这是判决的干法。

这一篇干的是后者。

先说清楚分工。知识库里已经有一篇《β衰变核结构:g_A 淬灭五十年谜题走到哪了》(2026-08-22),那是全景地图:三派学说、⁹⁰Zr、电荷交换、0νββ、天体物理、中国团队的位置,都铺开了。还有一篇《重离子敲出与谱因子淬灭》(2026-07-08)讲的是另一码事——敲出反应里那个"谱因子"的淬灭,机制不同,别混。本篇不重复前两篇的内容。本篇只回答一个问题:

那少掉的两三成,到底该记在谁头上?

而这个问题的答案,出乎意料地不在理论那边。它卡在一块石头的质量上——一块叫 ¹⁰⁰Sn 的石头,两个实验室量出来的"体重"差了 0.6 MeV,而这 0.6 MeV 直接决定"是两体流干的,还是核力之外还有别的东西"。

一、先把"淬灭"这两个字钉死

1.1 一个数字,和它的两个版本

β 衰变是核子里把一个质子换成中子(或者反过来),同时吐出一个电子和一个中微子。这件事的快慢,取决于两样东西:一个是"相空间"(能吐出多少能量、多少种配置),一个是"跃迁矩阵元"(初态和末态有多像)。

后一样东西里藏着一个数字,叫 g_A,轴矢耦合常数。自由中子测出来是 1.276。但核物理学家早就发现,如果拿这个 1.276 去算核内的 β 衰变,算出来总是比实测快。快多少?大约快 1/0.75 倍,也就是要用 0.75 乘一下才对得上。

这个 0.75 不是理论算出来的,是先被数据逼出来的。1973 年,D. H. Wilkinson 系统整理了镜像核与超允许 β 衰变,指出核介质里的轴向耦合常数需要重整化[2]。再往前,1973 年 M. Ericson、A. Figureau 和 C. Thévenet 提出了物理解释:核里到处都是π介子场,而轴向流正好耦合到π介子上,于是核内的 g_A 被π介子场"稀释"了[3]。这个思路很快被 M. Rho 接手,他用 Barshay-Brown-Rho 的π-核光学势模型,把淬灭和π介子凝聚的 Lorentz-Lorenz 修正直接挂上了钩[4]。两年后,E. Oset 和 M. Rho 把这个图像直接落到 GT 矩阵元上[5]。

所以"淬灭"(quenching)这个词从一开始就带着两层含义:一层是记账层面的——数据要求你打七五折;一层是物理层面的——谁吃掉了那 25%。

1.2 为什么这个数字这么要命

两个理由。

第一个理由在 0νββ(无中微子双β衰变)。这件事的衰变率正比于 g_A 的四次方[31]。g_A 从 1.276 缩到 1.0,衰变率就掉到原来的 0.38 倍;缩到 0.7,掉到 0.09 倍。这直接决定了:一台探测器跑了几年之后,你到底是能不能从"没看见"推出中微子质量的上限。而 g_A^eff 取 1.276 还是取 1.0,是两套不能直接比较的约定——这是文献里 0νββ 矩阵元数值"散开"的人为来源之一[31]。从头算的现代做法是:用裸的 g_A,另外把两体流显式算进去,这样两套记账终于可以对上[32][33]。

第二个理由在天体物理。超新星塌缩时的电子俘获率、进而前身星的中心电子丰度、铁核质量,都依赖 GT(Gamow-Teller,自旋-同位旋)强度,而壳模型计算里这个强度是被 0.75 平方过的。

1.3 怎么分量:求和规则是唯一的裁判

要判断"少了不少",得先知道"本来该有多少"。

这件事有一个模型无关的答案,叫 Ikeda 求和规则:S(β⁻) − S(β⁺) = 3(N − Z),其中 S 是 GT 强度对全部能量的积分[7]。这个式子只依赖中子和质子的数目差,跟核结构细节无关,1963 年给出,至今没人推翻。

实测呢?Ejiri 2025 年的综述性分析指出:用电荷交换反应测出来的 GT 总强度和自旋偶极(spin-dipole)总强度,只有基于核子的求和规则极限的一半左右[12]。也就是说,故事不是"低能区少了一点",而是"整个自旋-同位旋通道被砍掉了一半"。

"砍掉一半"这件事,就是全部争论的起点。

二、四个嫌疑人,一个卧底

要谈"根源判别",先把候选名单摆全。核物理界五十年里提出过四条路线,我把它们按"要不要引入核子以外的新自由度"排一排。

嫌疑人一号:Δ-isobar(Δ(1232) 共振)。 核子不是基本粒子,它有一个自旋 3/2、同位旋 3/2 的激发态,叫 Δ(1232),质量比核子重约 300 MeV。核内的核子可以被"踢"成 Δ,形成一个 Δ-空穴(Δh)激发。这个过程会分掉一部分自旋-同位旋强度,而因为 Δ 激发需要的能量很高(约 300 MeV),分走的那部分在低能区就"看不见"了——于是低能跃迁看起来被淬灭了。

嫌疑人二号:两体流(two-body currents)。 弱相互作用在核里不只是"一个个核子各自衰变"(那叫单体流),还有"两个核子通过交换介子一起参与"的过程——两个核子的轴矢流之间夹一个π介子,或者有一个短程接触项。这些叫两体流,或者老派叫介子交换流(MEC)。它们对跃迁振幅的修正是负的(削减),大小约 20%–30%。

嫌疑人三号:多体关联。 这个说法是:不是算符错了,是波函数错了。GT 强度并没有消失,只是被强关联推到了很高的激发能区,低能区看到的那一条只剩一点点。你用简单波函数算,自然会高估低能强度。

卧底:基本重整化(FR)。 这一派说:前面三个嫌疑人,本质上都是"核效应对算符和波函数的修正"。但如果三个都算完了、账还是不平,那说明g_A 本身在核介质里就被真空改变了,这不是"核效应",是 QCD 层面的东西。这一派的旗手是 Mannque Rho。

看清楚这四条,本篇要说的核心判断就可以先亮出来:

嫌疑人一号和二号,很可能不是两个人,是同一个人的两个名字。

三、嫌疑人一号:Δ-isobar —— 它一直躲在两体流的口袋里

3.1 历史:Δ-hole 曾经是主流解释

Δ 进入淬灭这个故事的时间很早。1982 年,E. Oset、H. Toki 和 W. Weise 在《Physics Reports》上发了一篇 100 页的长综述《Pionic modes of excitation in nuclei》[6],系统建立了"核子 + Δ(1232) 通过介子交换相互作用"的多体图像。这篇综述最后一部分专门讲什么?讲的就是 Δ(1232) 在低能 GT 跃迁和同位旋矢量磁自旋跃迁淬灭里的角色。结论是:在长波长极限下,Δ-空穴的屏蔽机制可以解释 GT 强度的淬灭。同期 I. S. Towner 在《Physics Reports》上把自旋矩阵元的淬灭系统整理成一本 115 页的手册[1],Δ-空穴和介子交换流都是里面的标准条目。

这一派给出的定量工具是 Landau-Migdal 参数。在有效粒子-空穴相互作用里,Δ-空穴通道的强度由一个参数 g'_{NΔ} 控制,而这个参数原则上可以从测量反推。

2001 年,A. Arima、W. Bentz 和两位同名的 T. Suzuki 干了这件事[8]。他们拿 ⁹⁰Zr(p,n)⁹⁰Nb 反应的 GT 求和规则淬灭数据,反推出 Landau-Migdal 参数 g'_{NΔ},再和 π+ρ 交换模型(显式考虑反对称化和短程关联)算出来的值对比。结论是:π+ρ 模型倾向给出比实验更强的淬灭,但如果在参数里用夸克模型值 f_Δ/f_π,两者在实验误差内一致。

注意这篇论文的作者单位很有意思,值得记一笔:Arima 当时的署名是东京参议院(House of Councillors, Tokyo,他那时已从东京大学退休转任参议员),Bentz 是东海大学(Tokai University,平塚),两位 Toshio Suzuki 分别署名日本大学和福井大学 + 理化学研究所和光。一篇论文里两个同名作者、四个不同单位,引用时最容易出错。

3.2 现在:Ejiri 给了 Δ 一个具体数字

Δ 路线最新的、也是最直白的一份工作,是 2025 年 10 月 20 日大阪大学核物理研究中心(RCNP)的 Hiroyasu Ejiri(江尻宏康)独著论文[12]。

这篇论文的论证链条非常清楚:

第一,电荷交换反应测出的 τσ(同位旋-自旋)总强度和自旋偶极总强度,相对核子求和规则极限减半。

第二,用有效 NN 和 NΔ τσ 相互作用的准粒子 RPA 分析,说明这个减半来自非核子的 Δ-isobar 共振效应。

第三,把效应参数化成一句极简的话:引入有效耦合

g^Δ_τσ / g_τσ ≈ 0.7

也就是说,Δ-isobar 把 τσ 跃迁矩阵元砍到自由核子值的约 0.7 倍。

第四,把这个 0.7 用到弱作用、电磁作用和核相互作用三类矩阵元的 τσ 分量上。

0.7 这个数字,请记住它。后面会看到它不是孤立的。

顺便记一个引用陷阱:这篇论文在 arXiv 上(2502.06122,2025-02-10 提交)的题目是《Delta-isobar resonance effects on nuclear matrix elements for neutrinoless ββ decays and astro-neutrino inverse-β decays》,而正式发表的题目改成了《Δ-isobar resonance effects studied by τσ summed strengths and nuclear matrix elements for β and ββ decays》[12]。你从 arXiv v1 抄题目,就会抄到一个不存在的标题。

再记一个:这篇论文的 DOI 是 10.1103/gtwh-dqmm,不是 10.1103/PhysRevC.112.045505。APS 从 2025 年起对新文章改用随机式 DOI,按"卷.页码"拼出来的 DOI 会直接 404。这类错误在二手材料里已经泛滥。

3.3 关键的一步:在 χEFT 里,Δ 被"打包"进了两个常数

现在到了这一节真正要讲的东西。

现代核力理论的主流框架叫手征有效场论(χEFT)。它的做法是把核子、π介子当作"活性自由度",而把 Δ 当作"不活性"的——通过"积掉"(integrate out)这个自由度,Δ 的影响就转移到了低能常数(LEC)里。

具体说:χEFT 的 NN 相互作用里有两个低能常数 c₃ 和 c₄,它们描述π-N 相互作用的特定结构。而这两个常数在数值上不小——为什么?因为它们被 Δ 共振主导。 物理上,这就是"Δ 的尾巴留在了 c₃、c₄ 里"。

那么问题来了:χEFT 的轴矢两体流,在次次领头阶(N²LO)有哪些项?主要就是 c₃、c₄ 的π交换项,加上一个短程接触项(系数记作 c_D,它和三核子力的接触项 c_E 是同一个常数,来自同一条手征拉氏量)。

结论已经很清楚了:

当你在两体流里画出那个 c₃、c₄ 的π交换图时,你算的就已经包含了 Δ 的一部分贡献。 "Δ-isobar 机制"和"两体流机制"在 χEFT 的记账法里不是两个独立的项,而是同一笔账的两种记法。

这不是我的推断,从反面也验证了:如果你显式地把 Δ 当成活性自由度放进 χEFT(所谓 Δ-full χEFT,或者叫 SSE 展开),你会得到同样好的物理,而且收敛更快——因为 Δ 与核子的质量差只有 300 MeV,不算大,硬把它积掉会让展开的收敛变慢。Pascalutsa、Vanderhaeghen 和 S. N. Yang 在《Physics Reports》上的 108 页 Δ(1232) 综述[17]里把这件事讲得很透:Δ 是核子最低的一级激发,任何想描述核子轴向性质的框架都绕不开它。

3.4 最硬的一条证据:氚的 GT 矩阵元

判"Δ 到底包不包含在两体流里",有一个极干净的实验案子:氚(³H)的 β 衰变。

氚是三个核子,从头算(第一性原理)算得动,实验精度也高。A. Baroni 等 11 人 2018 年在《Physical Review C》上的工作[16],做了这样一件事:用显式包含 Δ 中间态的手征 NN + 3N 相互作用构造三核子波函数,把轴矢流算到 N³LO(同样带 Δ)。

结果:预测值只比实验值高 1% 到 4%。

这句话的分量在于:三核子这么小的系统,"多体关联"没什么文章可做(就三个核子),如果淬灭还有二三十个百分点没解释,那一定是算符的问题。而现在,把 Δ 显式放进算符后,剩下的偏差只有个位数百分比。反过来说——Δ 的贡献,确实就在那个算符里。

反过来看,如果你用不带 Δ 的框架、只留到 N²LO,想凑出同样的结果,你就得让 c₃、c₄ 去干额外的活。也就是说:不带 Δ 的框架里,c₃、c₄ 的"有效值"会比带 Δ 时更大。 这就是"两体流解释了一部分淬灭"这句话的微观来源。

我还想补一条 2007 年的旁证:M. Procura、B. U. Musch、T. R. Hemmert 和 W. Weise 用格点 QCD 的夸克质量外推算核子 g_A 时发现,标准的 πN 手征微扰论到 p⁴ 阶根本无法把格点数据内插到物理点,而显式加入 Δ(1232) 的版本就成功了[28]。

一句话总结嫌疑人一号的现状:它不是一个独立的嫌疑人,而是"两体流"这个嫌疑人的曾用名。 但它在历史上留下了一个可检验的数字:约 0.7[12]。

四、嫌疑人二号:两体流 —— 从两个核子一路算到 ¹⁰⁰Sn

4.1 起点:2011 年那篇三页纸

两体流路线进入现代阶段的标志,是 J. Menéndez、D. Gazit 和 A. Schwenk 2011 年在 PRL 上的三页论文[14]。这篇论文的标题直接把两件事并排放在一起:《Chiral Two-Body Currents in Nuclei: Gamow-Teller Transitions and Neutrinoless Double-Beta Decay》。

它的贡献是:用 χEFT 推导的轴矢两体流,去修正轻核和中等质量核的 GT 跃迁,同时指出这件事对 0νββ 矩阵元意味着什么。当时的问题是:从头算方法算出的 GT 跃迁系统性偏快——注意,这时候"多体关联"已经被算进去了,还是偏快。多出来的那部分,就只能记在算符上。

4.2 2019 年:一次号称"解决五十年谜题"的计算

真正把这件事推到大标题上的,是 P. Gysbers 等 11 人 2019 年发表在《Nature Physics》上的工作[13]。这篇论文的标题用的是很重的词:《Discrepancy between experimental and theoretical β-decay rates resolved from first principles》(实验与理论 β 衰变率的差异,从头算解决了)。

它的作者阵容横跨北美和欧洲,值得逐单位记一遍,因为二手材料常把单位搞错:

  • TRIUMF(加拿大温哥华):P. Gysbers、J. D. Holt、P. Navrátil、S. R. Stroberg
  • 不列颠哥伦比亚大学:P. Gysbers
  • 橡树岭国家实验室(ORNL)物理部:G. Hagen、G. R. Jansen、T. D. Morris、T. Papenbrock
  • 田纳西大学:G. Hagen、T. D. Morris、T. Papenbrock
  • ORNL 国家计算科学中心:G. R. Jansen
  • ORNL 计算科学与工程部:T. D. Morris
  • 劳伦斯利弗莫尔国家实验室(LLNL)核与化学科学部:S. Quaglioni、K. A. Wendt
  • 达姆施塔特工业大学核物理研究所:A. Schwenk
  • GSI 亥姆霍兹重离子研究中心 EMMI:A. Schwenk
  • 马普核物理研究所(海德堡):A. Schwenk
  • 里德学院:S. R. Stroberg
  • 华盛顿大学:S. R. Stroberg

通讯作者是 G. Hagen。方法是耦合簇(CC)算多体波函数、VS-IMSRG 算有效算符、手征 EFT 提供相互作用和两体流。

结果:对一直到 ¹⁰⁰Sn 的核,把两体流算进去之后,理论与实验的 β 衰变率终于对上了。 论文摘要里的原话是:这个淬灭"在很大程度上来自弱力与两个核子的耦合,以及核内的强关联"[13]。文中给出的 ¹⁰⁰Sn 淬灭因子范围是 0.73 到 0.85。

这个 0.73–0.85,请也记住。

4.3 2022 年:把两体流搬到重核

两体流路线有一个天然的软肋:轻核算得动,重核怎么办?

E. M. Ney、J. Engel 和 N. Schunck 2022 年在《Physical Review C》上回答了这个问题[15],《Two-body weak currents in heavy nuclei》。三人的单位是北卡罗来纳大学教堂山分校物理与天文系(J. Engel、N. Schunck 系里)和LLNL 核与化学科学部。他们用 Skyrme 密度泛函的框架,把两体弱流做成可以在重核里推广的形式。

4.4 应用:0νββ 的记账方式被改了

两体流路线的实际影响,最集中地体现在 0νββ 上。

S. J. Novario 等 9 人 2021 年的 PRL[32],用耦合簇 + 手征 EFT 算 ⁴⁸Ca 的 0νββ 矩阵元,其中淬灭因子不再人为给,而是从 ¹⁰⁰Sn 区域两体流计算得到的 Ikeda 求和规则结果推出来。

A. Belley 等 12 人 2024 年的 PRL[33],把这个做法推到 ⁷⁶Ge——这是当前最受关注的 0νββ 候选核之一,作者阵容里有 TRIUMF、中山大学(J. M. Yao)、萨克雷 IRFU、北卡罗来纳大学。

这一路线的精神是:不要再用一个拍脑袋的 g_A^eff,而是用裸 g_A 加显式两体流。 这样,同一个理论框架既能算 β 衰变,也能算 0νββ,两套数字可以直接比较。当然,从头算方法本身也有能力边界——它之所以能从 A=2 一路推到锡区,靠的是手征相互作用与多体方法的组合,而不是无条件的普适性[34]。

五、嫌疑人三号:多体关联 —— "你以为你算的是 ¹⁰⁰Sn,其实不是"

5.1 反方最响的一炮

如果两体流是主角,那 2025 年 9 月出现的一篇预印本给了它一记闷棍。

H. Zhou(周昊)、L.-J. Wang(王龙军)和 Y. Sun(孙扬)的《Many-body correlations as the origin of Gamow-Teller quenching in nuclear β-decay》[28]——题目本身就是结论:多体关联才是淬灭的根源。

三位作者的单位是:西南大学物理科学与技术学院(H. Zhou、L.-J. Wang,王龙军为通讯作者)和上海交通大学物理与天文学院(Y. Sun)。

这篇论文的论证方式很有针对性。它说:最近的研究"过度强调算符重整化(包括手征两体流),而多体关联的作用——特别是在重开壳核里——被低估了"。然后它做了一件平衡的事:用大组态壳模型同时装进手征两体弱流和完整的核关联,让两个机制在同一套波函数里公平竞争。

案子选的是 ⁷⁶Ge。结论有三层:

第一,强核关联把大量 GT 强度推到了高激发能区,于是低能强度被大幅压低——这就是从外面看到的"淬灭"。

第二,这些关联的具体来源是形变、跨壳关联、以及密集高激发态之间的混合。

第三,也是最狠的一句:手征两体流的贡献只有 5% 到 15%,取决于所用的耦合常数。

一个说两体流是主角(Gysbers[13]),一个说两体流是龙套(Zhou 等[28])。这是目前这个领域最刺眼的对冲。

(注:这篇预印本 2025-09-29 挂出 v1,2026-04-01 更新到 v2。截至 2026-09-21,未在 CrossRef 查到期刊版。)

5.2 2026 年的一个折中方案

2026 年 5 月 20 日出现了第三种声音,做法是把两边都算进去再加起来。

论文题目本身就是一份和解协议:《The quenching of the axial-vector coupling constant g_A in β-decay: joint effects from chiral two-body currents and many-body correlations》(β 衰变中 g_A 的淬灭:手征两体流与多体关联的联合效应)[25]。

作者七人,单位分布很广,逐一记:

  • B.-L. Wang(王滨磊):北京师范大学物理学与天文学学院,射线束技术教育部重点实验室
  • W.-L. Lv(吕万利):兰州大学核科学与技术学院 + 稀有同位素前沿科学中心
  • L.-G. Cao(曹李刚):北京师范大学 + 北京市科学技术研究院辐射技术研究所
  • Y.-F. Niu(牛一斐):兰州大学核科学与技术学院 + 稀有同位素前沿科学中心
  • G. Colò:米兰大学物理系 + INFN 米兰分部
  • H. Sagawa:理化学研究所仁科中心 + 会津大学 + 中国科学院理论物理研究所
  • F.-S. Zhang(张丰收):北京师范大学 + 北京市科学技术研究院 + 兰州重离子加速器国家实验室理论核物理中心

方法:微观自洽 RPA + 粒子-振动耦合(PVC),两体流来自手征 EFT,在 ⁵⁶Ni、¹⁰⁰Sn、¹³²Sn 三个双幻核上做,用多种 Skyrme 能量密度泛函。

结果:把 PVC 带来的多体关联和手征两体流合起来,不需要任何额外调节,就能定量复现实验数据;提取出的淬灭因子落在 q ≈ 0.73 到 0.80。

这个 0.73–0.80,请第三次记住。

5.3 还有一条"有效算子"路线

还有一派人走的是第四条路:既不改 g_A,也不改 g_A,而是问"壳模型里那个简化的 GT 算子,和真实弱流算子差多少"。

代表是 L. Coraggio、N. Itaco、G. De Gregorio、A. Gargano、Z. H. Cheng(程志辉)、Y. Z. Ma、F. R. Xu(许甫荣,北京大学)、M. Viviani 八人 2024 年在《Physical Review C》上的工作[24]。他们用 Q-box 折叠图这种多体微扰方法,从手征 EFT 的有效场论出发,把壳模型 GT 算子本身重新归一化。那不勒斯 INFN 是主力。

这条路线的好处是:它算出来的"有效算子"可以直接塞进大规模的壳模型计算里,比如 0νββ 的核矩阵元。

5.4 三个嫌疑人的数字排在一起

到这里可以做个表了。三个嫌疑人给出的"淬灭因子",用的是完全不同的方法、完全不同的核:

路线代表工作核给出的淬灭量
Δ-isobar(Δ-hole)Ejiri 2025 [12]电荷交换反应综合g^Δ/g ≈ 0.7
两体流 + 强关联(从头算)Gysbers 等 2019 [13]轻核到 ¹⁰⁰Snq ≈ 0.73–0.85
两体流 + 多体关联(RPA+PVC)Wang 等 2026 [25]⁵⁶Ni、¹⁰⁰Sn、¹³²Snq ≈ 0.73–0.80
多体关联为主(大组态壳模型)Zhou 等 2025/26 [28]⁷⁶Ge两体流仅 5–15%
有效算子重整化Coraggio 等 2024 [24]壳模型有效算子视模型空间而定
Δ 显式、三核子Baroni 等 2018 [16]³H只剩 1–4% 偏差

注意第一、二、三行的数字:0.7、0.73–0.85、0.73–0.80。 三条完全不同的路子,指向同一个区间。这不像是三个独立的机制,倒像是同一个机制被三把不同的尺子量了三遍。

这正好呼应第三节的结论:Δ-isobar 与两体流,在 χEFT 的记账法里是同一笔账。 Δ 的贡献被积掉、转写进了 c₃、c₄;你在两体流里用 c₃、c₄ 画π交换图,就已经在算它了。反过来,Ejiri 那种"显式 NΔ 相互作用"的做法,得到的也是同一个约 0.7。

有一个独立佐证很支持这个看法:⁹⁰Zr 的电荷交换测量给出的淬灭因子是 Q = 0.88 ± 0.06[10]。这是 K. Yako 等 20 人用 ⁹⁰Zr(n,p) 反应在 293 MeV 下测的,能量窗口覆盖 0 到 70 MeV,用多极分解把 GT 成分从连续谱里抠出来。0.88 比 0.75 轻得多——为什么?因为能量窗口不同。 如果你把高激发区也算进来,被"挪走"的强度就有一部分回来了。这说明"淬灭"这个数字本身带窗口依赖,不是一个绝对的常数[9][10][11]。

(顺带纠一个作者顺序的错:那篇用 Skyrme RPA 提取 GT 与自旋偶极淬灭因子的 PRC 100, 054324 (2019)[11],作者顺序是 L.-G. Cao、S.-S. Zhang、H. Sagawa,不是常被写成"Sagawa 打头"。三个人的署名单位分别是华北电力大学数理学院、北京航空航天大学物理与核能工程学院、理化学研究所仁科中心 + 会津大学。同一个曹李刚,2026 年在[25]里的署名单位已经是北京师范大学了——引用时必须按论文发表时的单位写。)

六、卧底:如果前三个都算完了,账还是不平

现在轮到第四条路线。

6.1 Rho 的"基本重整化"

Mannque Rho 的立场,用一句话说是:你们说的那个"淬灭",有一半根本不是淬灭。

他的逻辑是这样的。核物理文献里说的"淬灭的 g_A",其实是指"用单个核子的轴向流、配上简单波函数去算核跃迁时,需要的那个经验缩放因子"。但真正该问的问题是:在核介质里,单个核子的轴向流本身,还是自由空间里那个吗?

Rho 的答案是:不是。原因是 QCD 的尺度对称性破缺(scale anomaly)。在核介质里,真空结构变了,迹反常(trace anomaly)通过胶子能量动量张量的反常维度 β',改变真空的性质,进而改变在介质中"跑动"的 g_A。他指出这一切有一个红外不动点(IR fixed point)在背后支撑——在无质量夸克极限下,这个不动点出现在 Nambu-Goldstone 模式里。

这个思路的关键一步,写在 Y.-L. Ma(马永亮)和 M. Rho 2020 年的 PRL 上[21]:《Quenched g_A in Nuclei and Emergent Scale Symmetry in Baryonic Matter》。触发点正是 RIKEN 对 ¹⁰⁰Sn 的测量结果。 两位作者当时的署名单位是吉林大学理论物理中心与物理学院、国科大杭州高等研究院基础物理与数学科学学院、国科大亚太理论物理中心(ICTP-AP,北京),以及巴黎萨克雷大学理论物理研究所(IPhT,CEA/CNRS)。

Rho 2021 年在《Symmetry》上又把这件事展开了一遍[22],把"淬灭"这个词专门加了引号,并明确说:这不是壳模型意义上的淬灭,而是 QCD 真空的改变对 g_A 的"基本"重整化。

6.2 2024 年,把话说到最重

Rho 和 L.-Q. Shao(邵龙奇)2024 年 10 月的预印本[23],把论证压到了最紧的一段。两位作者的署名是:

  • Mannque Rho:巴黎萨克雷大学、CNRS、CEA,理论物理研究所 IPhT(Gif-sur-Yvette)
  • Long-Qi Shao:国科大杭州高等研究院基础物理与数学科学学院 + 中国科学院理论物理研究所 + 中国科学院大学

论文的论证结构只有三步:

第一步,把 GSI 和 RIKEN 的 ¹⁰⁰Sn 数据放在一起:B(GT) 分别是 9.1(+3.0/−2.6) 和 4.4(+0.9/−0.7)。

第二步,用"极端单粒子壳模型"(ESPM,对双幻核 ¹⁰⁰Sn 适用)换算成淬灭因子。因为 ¹⁰⁰Sn 的 GT 跃迁几乎是一条跃迁(π0g₉/₂ 上完全填满的质子,掉进完全空的 ν0g₇/₂),所以

q = sqrt( B(GT)_实测 / 17.78 )

这里 17.78 是 ESPM 给出的裸值。

第三步,得到:q^GSI ≈ 0.73,q^RIKEN ≈ 0.50。

然后指出:如果这个 ~30%–40% 的差别不是被误差棒淹没的假象,那它不能被两体流和多体关联解释掉——因为那些机制的定量贡献我们已经量过,就是前面表里的 0.73 到 0.85。差出来的那一截,指向"基本重整化"。

论文最后点出这件事的杀伤力在哪:0νββ 的衰变率正比于 g_A 的四次方,如果 g_A 真的被"基本重整化"到 0.5 附近,0νββ 的预期信号要弱一个数量级以上[23]。

这就是卧底的力量:它不需要解释 ¼ 的淬灭,它只需要证明剩下的那一截没人能认领。

七、¹⁰⁰Sn:一块石头,两条裂缝

现在,全部重量都压在 ¹⁰⁰Sn 这个案子上。

7.1 为什么是它

¹⁰⁰Sn 是最重的、还能稳定存在的 N=Z 双幻核(50 个质子,50 个中子)。它的特殊性来自"极端单粒子壳模型"的适用性:

  • 质子的 π0g₉/₂ 壳完全填满;
  • 中子的 ν0g₇/₂ 壳完全空着;
  • 于是 β⁺ 衰变几乎只有一条通道:一个 π0g₉/₂ 质子翻成 ν0g₇/₂ 中子。

结果就是"超允许 GT 跃迁":GT 强度不需要分散到几十条跃迁上去,它集中在一条上。这意味着,你只要测出半衰期和衰变能,就能直接读出总 GT 强度。这是核物理里少有的干净案子。ESPM 给的裸值是 B(GT) = 17.78[20]。

这一双幻性质在近年得到独立支持:对铟(Z=49)同位素的精密激光谱学测出,电荷半径与四极矩随中子数呈抛物线趋势,在 N=50 与 N=82 两个闭壳处都明显下沉[27]。也正因为如此,它成了检验淬灭的"理想实验室"——它把"淬灭到底多大"这个问题,压缩成了一个可以精测的数。

7.2 GSI 2012:259 个核

2012 年,C. B. Hinke 等 62 人的国际合作组,在 GSI 亥姆霍兹重离子研究中心做成了这件事[19]。他们的办法是:用 1 GeV/A 的 ¹²⁴Xe 初级束流打铍靶,弹核碎裂产生 ¹⁰⁰Sn,用碎片分离器挑出来,再植入高度分段的硅条探测器里做衰变谱学,外面包一圈锗探测器阵列看γ射线。

数字:

  • 成功产出并观测 259 个 ¹⁰⁰Sn,比此前任何实验多一个数量级(此前约 70 个衰变的量级);
  • 半衰期 T₁/₂ = 1.16 ± 0.20 s;
  • β 端点能量 3.29 ± 0.20 MeV;
  • 得到的 GT 强度 B(GT) = 9.1。

论文摘要用的是很重的措辞:"迄今为止在允许核 β 衰变中测得的最大 GT 强度",确立了这条跃迁的"超允许"性质[19]。

误差棒要小心:期刊版正文/补充材料里广泛引用的是 +3.0/−2.6(Lubos 等[20]和 Rho & Shao[23]都按这个引);而 Hinke 的博士论文里给的是 +4.8/−2.3,中心值同为 9.1。两个误差棒来自不同的分析场景。引用时必须注明你用的是哪一个版本。

7.3 RIKEN 2019:2500 个核,结果反而更小

2019 年,D. Lubos 等 52 人的 EURICA 合作组在理化学研究所仁科加速器研究中心重做了这件事[20]。他们的探测器是 WAS3ABi + EURICA,检出约 2500 个 ¹⁰⁰Sn——是 GSI 那次约十倍。

数字:

  • 半衰期 T₁/₂ = 1.18 ± 0.08 s(误差比此前平均值小一半);
  • β 端点能量 Q_β = 3.91 ± 0.15 MeV;
  • 得到的 GT 强度 B(GT) = 4.4(+0.9/−0.7)。

作者在理研的年度报告里特意写了这么一句:这次的数据"第一次使得区分不同模型计算成为可能",并且"在所有理论中,手征有效场论似乎最能够描述数据"[20]。

于是,问题变成了这样:

GSI 2012 [19]RIKEN 2019 [20]
统计量259 个核~2500 个核
半衰期1.16 ± 0.20 s1.18 ± 0.08 s
β 端点能量3.29 ± 0.20 MeV3.91 ± 0.15 MeV
B(GT)9.1(+3.0/−2.6)4.4(+0.9/−0.7)
换算淬灭因子 q≈ 0.73≈ 0.50

半衰期几乎一致(都在 1.17 s 附近)[19][20]。真正的分歧全在 β 端点能量:一个是 3.29 MeV,一个是 3.91 MeV,差了 0.62 MeV。 两个误差一合并,差别接近 2σ。Lubos 等自己也承认:新值比此前最好的结果更大,差值"约为两边误差之和的两倍",并说"这个不一致可能来自两次测量中都未被计入的系统误差"[20]。

7.4 为什么 0.62 MeV 会翻倍地改变结论

关键在 B(GT) 和 f t 值的关系。GT 强度反比于比较半衰期 ft:

B(GT) ∝ 1 / (f · t₁/₂)

其中 f 是相空间积分,它强烈依赖衰变能。端点能量多 0.62 MeV,相空间 f 就变大一截,于是从同一个半衰期推出来的 B(GT) 就变小一截。 这就是为什么两个组的半衰期几乎一样,B(GT) 却差了一倍。

再往前一步:q² = B(GT)_实测 / 17.78,所以

  • 若 B(GT) = 9.1 → q ≈ 0.72
  • 若 B(GT) = 4.4 → q ≈ 0.50

"根源判别"的胜负,就这样被压缩成了一个 0.62 MeV 的端点能量差上。 而这是一个纯实验问题。

八、裁决:两次独立的质量测量,都指向同一个方向

β 端点能量难测,因为它是一个连续谱的端点,还牵扯到 ¹⁰⁰Sn 的质量(而 ¹⁰⁰Sn 太难直接称重)。但有一条绕路可以走:去称它的邻居。

¹⁰⁰Sn 的质量,可以间接地从 ¹⁰⁰In 的质量加上 β 衰变能推出来。¹⁰⁰In 是可以精测的。 于是问题转化为:拿 GSI 的 3.29 MeV 去推,和拿 RIKEN 的 3.91 MeV 去推,会得到两个不同的 ¹⁰⁰Sn 质量。哪个更合理,可以拿核结构的趋势去检验。

这种"用结构性判据挑选实验值"的裁决,做过两次。两次都指向同一个方向。

8.1 第一次:2021 年 ISOLTRAP 称 ¹⁰⁰In

M. Mougeot 等 26 人 2021 年在《Nature Physics》上的工作[18],用 CERN-ISOLDE 的 ISOLTRAP 质量谱仪精确测量 ¹⁰⁰In 的质量。面对的困难是:产额极低、污染严重。他们用的是当时最灵敏的多反射飞行时间(MR-TOF)技术。

成果:¹⁰⁰In 质量的精度提高了 90 倍。

然后是裁决部分。他们用两个文献中的 β 衰变能,分别算出 ¹⁰⁰Sn 的质量,再算一个核结构量——双中子经验壳隙 Δ2n,也就是 N=50 壳隙,把它和附近 Z 更小的同类核素的趋势比较。

结果:

  • 用 GSI 的值算出来的 ¹⁰⁰Sn 质量,得到的 Δ2n 跟随轻核素设定的趋势;
  • 用 RIKEN 的值算出来的,在 Z=50 处出现一个陡峭的跳变,用核结构模型很难解释。

Mougeot 等因此判定:他们的新测量支持 GSI 的值[18]。

这件事的物理意义:如果 ¹⁰⁰Sn 真的"特别不束缚"(RIKEN 的值意味着它更不束缚),那 Z=50 的壳隙就会反常。而 Z=50 是核物理里最结实的幻数之一,它在锡同位素链上被反复验证过(例如激光谱学在 N=82 处看到的电荷半径不连续)。为了迁就一个 β 端点能量,把 Z=50 的壳隙推翻,代价太大。

8.2 第二次:2026 年 FRIB 称 ¹⁰¹Sn

如果只有一次裁决,可以说是一家之言。2026 年出现了第二次。

C. M. Ireland 等 26 人的工作,发表在《Physical Review C》113(3), L021302 (2026)[26]——注意这篇在 arXiv 上的题目《Mass of ¹⁰¹Sn and Bayesian extrapolations to the proton drip line》(arXiv:2510.11815)和期刊题目基本一致,但期刊版要引用随机式 DOI 10.1103/vz5k-yf93。

他们的做法是:

第一,在密歇根州立大学 FRIB(稀有同位素束流装置)做 ¹⁰¹Sn 的首次彭宁阱质量测量,得到质量过剩 −59 889.89(96) keV,精度比此前提高 300 倍(约 0.1 keV 量级)。

第二,用七种基于核密度泛函(DFT)的全球质量模型,加上统计机器学习,搭一个贝叶斯模型组合(BMC)框架,先验证它在 48 ≤ Z ≤ 52 的链上都能在 1σ 内复现实验质量。

第三,用这个框架外推,给出 ¹⁰⁰Sn 的质量过剩预测为 −58 090(800) keV。

结论一句话:这个预测"在 1σ 内偏向 GSI 测出的 β 延迟 Q 值所推出的 ¹⁰⁰Sn 质量"[26]。

(论文里还有一句话值得单独记:他们指出,AME2020 里采用的 ¹⁰⁰Sn 质量,就是从 GSI 的 β 衰变能间接算出来的[26]。也就是说,国际质量评估体系早就在这两组数据里选了 GSI。)

8.3 天平现在偏向哪边

三次独立的、方向一致的信号:

  • 2021 年 ISOLTRAP 精测 ¹⁰⁰In(精度提升 90 倍)→ 核结构趋势支持 GSI[18]
  • 2026 年 FRIB 精测 ¹⁰¹Sn(精度提升 300 倍)+ 贝叶斯模型组合外推 → 1σ 内偏向 GSI[26]
  • AME2020 原子质量评估采用的 ¹⁰⁰Sn 质量源自 GSI 的 β 衰变能[26]

如果 GSI 的端点能量(3.29 MeV)是对的,那么:

  • B(GT) = 9.1;
  • q = 0.73;
  • 而两体流 + 多体关联给的是 0.73–0.85[13] 和 0.73–0.80[25];
  • Δ-isobar 那条独立路线给的是 0.7[12]。

三者重合。账平了。

反过来,如果 RIKEN 的 3.91 MeV 是对的,那么 q = 0.50,比所有已知机制能提供的淬灭少了约 30%–40%。那一截没人认领,卧底就该出场了[23]。

本篇的判断:目前的证据链偏向 GSI,也就是偏向"两体流 + 多体关联已经够了"这一侧。 但要说清楚——这不是判决,是举证责任的方向变了。RIKEN 那次统计量大约十倍于 GSI[19][20],在统计上更强;而它的弱点集中在端点能量这个系统性问题上。一个统计更强的测量在一个特定的系统量上偏离,比一个统计更弱的测量整体偏离,更常见,也更容易被后续修正。要真正结案,需要的是一次独立的、高精度的 ¹⁰⁰Sn β 端点能量测量——这件事目前还没做。

至于 Rho & Shao 那篇点出的"30%–40% 的差别是不是被误差棒淹没的假象"——他们自己也在文里承认:他们"不引用误差棒,因为在所引结果里注入了理论论证,其可靠性在我们的框架里难以量化"[23]。这句话是对同行的批评,也是对自己论证的限定。

九、这件事对谁最要命

9.1 0νββ:g_A 的四次方

最直接的后果在 0νββ。衰变率 ∝ (g_A^eff)⁴[29]。q 取 0.73 还是 0.50:

  • q = 0.73 → (0.73)⁴ ≈ 0.28
  • q = 0.50 → (0.50)⁴ ≈ 0.0625

两者相差 4.5 倍。这意味着:同一台探测器、同样的本底、同样的曝光量,如果 g_A 真的是 0.5,那么能达到的中微子质量灵敏度会差一倍以上——因为信号率掉到 1/4.5,你要么加四五倍曝光,要么接受差一倍以上的质量上限。

值得注意的是,现代从头算的做法已经避免了这个"人为的散开":用裸 g_A 加显式两体流[13][32][33]。Belley 等 2024 年的 ⁷⁶Ge 工作[33]同时给出统计误差和系统误差(模型空间、相互作用、算符截断),这是让 0νββ 矩阵元"可比较"的正确方向。

9.2 天体物理:超新星和中子星

超新星塌缩时的电子俘获率直接决定前身星的中心电子丰度和铁核质量。壳模型弱率用的 GT 强度是被 (0.75)² 标度过的。如果真实淬灭更浅(q 更大),电子俘获会更快;如果更深,会更慢。这个差别会沿着一整条恒星演化链传下去。

9.3 逆 β 衰变与中微子探测

Ejiri 那篇工作的另一个落点是逆 β 衰变(中微子-核反应)[12]。这是超新星中微子探测和"astrо-中微子"研究的基础截面。他把 0.7 这个因子同时应用到弱、电磁、核相互作用三类矩阵元的 τσ 分量上,理由就是它们共享同一个 τσ 相关性。

十、还没解决的,和中国团队在做什么

10.1 三件没解决的事

第一,¹⁰⁰Sn 的 β 端点能量还没有独立复测。 这是目前最硬的一块缺口。GSI 的 3.29(20) 和 RIKEN 的 3.91(15) 之间,需要一个第三方的高精度测量来判决。而它很难——¹⁰⁰Sn 的产额极低,端点又靠近连续谱的高端。

第二,两体流到底占多少,仍是开放的。 Gysbers 等说它"很大程度上"解释淬灭[13];Wang 等 2026 说它和多体关联"合起来"是 0.73–0.80[25];Zhou 等说它单独只有 5%–15%[28]。三个说法并不必然矛盾——关键在于"多体关联"在不同框架里的定义不同。在 VS-IMSRG 里,被吸收进有效算符的关联,和在大组态壳模型里被显式展开的关联,不是同一批。这个"定义差"是当前文献混乱的主要来源,也是我建议读者读这些论文时的第一道过滤。

第三,Δ 与两体流的记账重叠还没有一个统一的定量说法。 本篇的判断是"它们是同一笔账"(见第三节),但严格说,这需要有人把同一批图用两种记账法各算一遍、给出严格的重复计数估计。目前还没有这样的工作。

10.2 中国团队的位置

有意思的是,在"根源判别"这场辩论里,中国团队已经从"参与"走到了"主导某一侧":

  • 西南大学 + 上海交大(H. Zhou、L.-J. Wang、Y. Sun)站在"多体关联是主因"这一侧,是反方最响的一炮[28]。
  • 北京师范大学 + 兰州大学 + 米兰大学 + 理研/会津大学 + 中科院理论所的多单位组合(Wang、Lv、Cao、Niu、Colò、Sagawa、Zhang)提出"联合效应"的折中方案,把 q 钉在 0.73–0.80[25]。
  • 北京大学许甫荣参与的壳模型有效算子重整化[24],提供了把理论接到大尺度壳模型计算上的桥梁。
  • 马永亮 + Rho 站在"基本重整化"这一侧[21],这条路线从 2020 年一直走到 2024 年的 IR 不动点论证[23]。
  • 北京大学(罗迪雯一作,李智焕/华辉通讯,华东师大程奕源共同通讯) 在 ¹³²Sn 南邻 ¹²⁸Ag 上的工作,是这一区域结构信息的重要补充——提醒一下这篇的一个常见错误:¹²⁸Ag 的中子数是 81,不是 82,而且该工作近物所并未参与。

短板仍在实验。 ¹⁰⁰Sn 那个关键端点能量,目前只能靠 GSI 和 RIKEN 的历史数据,加上 ISOLTRAP 和 FRIB 的间接裁决。中国自己的 HIAF 和 RIBLL2 在电荷交换测量上还在建设,短期内替不上。

10.3 怎么读这个领域的论文

最后给读者三条实用的过滤规则:

规则一:先问"这个淬灭因子是对哪个能量窗口说的"。 ⁹⁰Zr 的 0.88[10] 和 ¹⁰⁰Sn 的 0.73 不是同一件事。窗口不同,数字不可直接比较。

规则二:先问"多体关联"指的是什么。 是显式展开的组态空间里的关联,还是被折叠进有效算符的关联?这两个数不能相加,否则重复计数。

规则三:看到 g_A^eff 这个符号,先问它是不是把两体流藏进去了。 如果一篇论文用 g_A^eff < 1.276 又同时说"不含两体流",那它已经把两体流的效应(其中包含 Δ)混进有效值里了。这是最常见的重复计数来源。

十一、结语

回到开头那个问题:少掉的两三成,记在谁头上?

本篇的回答是:

第一,"Δ-isobar"和"两体流"在物理上不是两个嫌疑人。 在 χEFT 里,Δ 被积掉、转写进了 c₃ 和 c₄;两体流里的π交换项用的正是这两个常数。三条独立路线——Ejiri 的 NΔ 相互作用给出 0.7[12]、Gysbers 的从头算给出 0.73–0.85[13]、Wang 等的 RPA+PVC 给出 0.73–0.80[25]——落进同一个区间,这本身就是"同一笔账"的最强旁证。Procura 等的格点外推工作[28]和 Baroni 等的三核子计算[16]从两个方向支持这一点。

第二,"多体关联"和"两体流"的分工还没定,而且可能永远定不下来,因为它们在不同框架里的边界不同。但它对实际应用的影响有限——Wang 等的"联合效应"工作[25]说明:不管你怎么分摊,合起来都是 0.73–0.80。

第三,也是最关键的:这场"根源判别"的胜负手,不在理论,在一块石头的质量上。 ¹⁰⁰Sn 的 β 端点能量,3.29 MeV[19] 还是 3.91 MeV[20],决定了 q 是 0.73 还是 0.50;而 0.73 能被现有机制吃干净,0.50 不能——那多出来的 30%–40% 就得记在"基本重整化"头上[21][23]。

而目前的天平偏向 3.29 MeV:ISOLTRAP 对 ¹⁰⁰In 的 90 倍精度提升支持它[18],FRIB 对 ¹⁰¹Sn 的 300 倍精度提升加贝叶斯外推支持它[26],AME2020 采用的也是它[26]。

所以,本篇的判断是:如果证据链就这样走下去,"g_A 淬灭需要一个 QCD 层面的新解释"这个说法会失去它的主要实验支柱。 但那不是定论——只差一次对 ¹⁰⁰Sn β 端点能量的独立精测。 在那次测量做出来之前,这个五十年的谜题的最后一句话,还写在实验记录本上,而不是写在论文摘要里。

---

参考文献

  1. I. S. Towner, *Quenching of spin matrix elements in nuclei*, Physics Reports 155(5), 263–377 (1987). DOI 10.1016/0370-1573(87)90138-4。【已 DOI 直查;单作者】
  1. D. H. Wilkinson, *Renormalization of the Axial-Vector Coupling Constant in Nuclear β Decay*, Physical Review C 7(3), 930–936 (1973). DOI 10.1103/PhysRevC.7.930。【已 DOI 直查;单作者】
  1. M. Ericson, A. Figureau, C. Thévenet, *Pionic field and renormalization of the axial coupling constant in nuclei*, Physics Letters B 45(1), 19–22 (1973). DOI 10.1016/0370-2693(73)90242-6。【已 DOI 直查;3 作者】
  1. M. Rho, *Quenching of axial-vector coupling constant in β-decay and pion-nucleus optical potential*, Nuclear Physics A 231(3), 493–503 (1974). DOI 10.1016/0375-9474(74)90512-0。【已 DOI 直查;单作者;署名单位为当时的纽约州立大学石溪分校】
  1. E. Oset, M. Rho, *Axial Currents in Nuclei: The Gamow-Teller Matrix Element*, Physical Review Letters 42(1), 47–50 (1979). DOI 10.1103/PhysRevLett.42.47。【已 DOI 直查;2 作者】
  1. E. Oset, H. Toki, W. Weise, *Pionic modes of excitation in nuclei*, Physics Reports 83(4), 281–380 (1982). DOI 10.1016/0370-1573(82)90123-5。【已 DOI 直查;3 作者;署名单位为萨拉曼卡大学、密歇根州立大学回旋加速器实验室、雷根斯堡大学理论物理所】
  1. K. Ikeda, S. Fujii, J. I. Fujita, *The (p,n) reactions and beta decays*, Physics Letters 3, 271–272 (1963). DOI 10.1016/0031-9163(63)90255-5。【已 DOI 直查;3 作者】
  1. A. Arima, W. Bentz, T. Suzuki, T. Suzuki, *Delta–hole interaction in nuclei and the Gamow–Teller strength in ⁹⁰Nb*, Physics Letters B 499(1–2), 104–108 (2001). DOI 10.1016/S0370-2693(00)01385-X;arXiv:nucl-th/0008029。【已 DOI 直查;4 作者(含两位同名 T. Suzuki);署名单位为东京参议院、东海大学、日本大学、福井大学 + 理研和光】
  1. M. Ichimura, H. Sakai, T. Wakasa, *Spin–isospin responses via (p,n) and (n,p) reactions*, Progress in Particle and Nuclear Physics 56(2), 446–531 (2006). DOI 10.1016/j.ppnp.2005.09.001。【已 DOI 直查;3 作者】
  1. K. Yako 等 20 人, *Determination of the Gamow–Teller quenching factor from charge exchange reactions on ⁹⁰Zr*, Physics Letters B 615(3–4), 193–199 (2005). DOI 10.1016/j.physletb.2005.04.032;arXiv:nucl-ex/0411011。结论 Q = 0.88 ± 0.06。【已 DOI 直查;署名单位含东京大学、大阪大学 RCNP、国际基督教大学、LANL、东北大学 CYRIC、俄亥俄大学、理研、九州大学】
  1. L.-G. Cao(曹李刚), S.-S. Zhang(张时声), H. Sagawa, *Quenching factor of Gamow-Teller and spin dipole giant resonances*, Physical Review C 100(5), 054324 (2019). DOI 10.1103/PhysRevC.100.054324;arXiv:1911.09795。注意作者顺序为 Cao 打头。 署名单位:华北电力大学数理学院、北京航空航天大学物理与核能工程学院、理化学研究所仁科中心、会津大学。【已 DOI 直查】
  1. H. Ejiri(江尻宏康,独著), *Δ-isobar resonance effects studied by τσ summed strengths and nuclear matrix elements for β and ββ decays*, Physical Review C 112(4), 045505 (2025-10-20). DOI 10.1103/gtwh-dqmm(随机式,非按卷期页拼);arXiv:2502.06122。arXiv v1 题名为 *Delta-isobar resonance effects on nuclear matrix elements for neutrinoless ββ decays and astro-neutrino inverse-β decays*,与期刊题名不同。署名单位:大阪大学核物理研究中心(RCNP)。【已 DOI 直查(随机式 DOI)+ 官网目次】
  1. P. Gysbers, G. Hagen, J. D. Holt, G. R. Jansen, T. D. Morris, P. Navrátil, T. Papenbrock, S. Quaglioni, A. Schwenk, S. R. Stroberg, K. A. Wendt(11 人), *Discrepancy between experimental and theoretical β-decay rates resolved from first principles*, Nature Physics 15(5), 428–431 (2019). DOI 10.1038/s41567-019-0450-7;arXiv:1903.00047。通讯作者 G. Hagen。署名单位:TRIUMF、不列颠哥伦比亚大学、橡树岭国家实验室(物理部 / 国家计算科学中心 / 计算科学与工程部)、田纳西大学、劳伦斯利弗莫尔国家实验室(核与化学科学部)、达姆施塔特工业大学、GSI-EMMI、马普核物理研究所、里德学院、华盛顿大学。【已 DOI 直查 + 官网作者单位页】
  1. J. Menéndez, D. Gazit, A. Schwenk, *Chiral Two-Body Currents in Nuclei: Gamow-Teller Transitions and Neutrinoless Double-Beta Decay*, Physical Review Letters 107(6), 062501 (2011). DOI 10.1103/PhysRevLett.107.062501。【已 DOI 直查;3 作者】
  1. E. M. Ney, J. Engel, N. Schunck, *Two-body weak currents in heavy nuclei*, Physical Review C 105(3), 034349 (2022). DOI 10.1103/PhysRevC.105.034349。署名单位:北卡罗来纳大学教堂山分校物理与天文系、劳伦斯利弗莫尔国家实验室核与化学科学部。【已 DOI 直查(含机构字段)】
  1. A. Baroni, R. Schiavilla, L. E. Marcucci, L. Girlanda, A. Kievsky, A. Lovato, S. Pastore, M. Piarulli, S. C. Pieper, M. Viviani, R. B. Wiringa(11 人), *Local chiral interactions, the tritium Gamow-Teller matrix element, and the three-nucleon contact term*, Physical Review C 98(4), 044003 (2018). DOI 10.1103/PhysRevC.98.044003;arXiv:1806.10245。结论:N³LO 含 Δ 的轴矢流对氚 GT 矩阵元只高估 1%–4%。【已 DOI 直查;注意 DOI 尾号是 044003,不是 054003】
  1. V. Pascalutsa, M. Vanderhaeghen, S. N. Yang, *Electromagnetic excitation of the Δ(1232)-resonance*, Physics Reports 437(5–6), 125–232 (2007). DOI 10.1016/j.physrep.2006.09.006;arXiv:hep-ph/0609004。署名单位:ECT(意大利 Trento)、杰斐逊实验室理论中心、威廉与玛丽学院、台湾大学(中国台湾)。【已 DOI 直查】
  1. M. Mougeot 等 26 人, *Mass measurements of ⁹⁹⁻¹⁰¹In challenge ab initio nuclear theory of the nuclide ¹⁰⁰Sn*, Nature Physics 17(10), 1099–1103 (2021). DOI 10.1038/s41567-021-01326-9(尾号 01326-9,不是 01333-w)。装置:ISOLTRAP @ CERN-ISOLDE。结论:¹⁰⁰In 质量精度提升 90 倍,核结构趋势支持 GSI 的 β 衰变能。【已 DOI 直查】
  1. C. B. Hinke 等 62 人, *Superallowed Gamow–Teller decay of the doubly magic nucleus ¹⁰⁰Sn*, Nature 486(7403), 341–345 (2012). DOI 10.1038/nature11116。数字:259 个 ¹⁰⁰Sn;T₁/₂ = 1.16 ± 0.20 s;β 端点能量 3.29 ± 0.20 MeV;B(GT) = 9.1。误差棒:期刊版引用 +3.0/−2.6;Hinke 博士论文为 +4.8/−2.3。 署名单位以慕尼黑工业大学物理系 E12 与 GSI 为主。【已 DOI 直查 + GSI 机构库记录 + 博士论文】
  1. D. Lubos, J. Park, T. Faestermann, R. Gernhäuser, R. Krücken 等 52 人(EURICA 合作组), *Improved Value for the Gamow-Teller Strength of the ¹⁰⁰Sn Beta Decay*, Physical Review Letters 122(22), 222502 (2019). DOI 10.1103/PhysRevLett.122.222502。数字:约 2500 个 ¹⁰⁰Sn;T₁/₂ = 1.18 ± 0.08 s;Q_β = 3.91 ± 0.15 MeV;B(GT) = 4.4(+0.9/−0.7)。署名单位含慕尼黑工业大学、理化学研究所仁科中心、隆德大学、TRIUMF、GANIL、东京大学。【已 DOI 直查 + RIKEN 年度报告】
  1. Y.-L. Ma(马永亮), M. Rho, *Quenched g_A in Nuclei and Emergent Scale Symmetry in Baryonic Matter*, Physical Review Letters 125(14), 142501 (2020). DOI 10.1103/PhysRevLett.125.142501;arXiv:2002.03310。署名单位:吉林大学理论物理中心与物理学院、国科大杭州高等研究院、国科大 ICTP-AP(北京)、巴黎萨克雷大学理论物理研究所(IPhT)。【已 DOI 直查】
  1. M. Rho(独著), *Multifarious Roles of Hidden Chiral-Scale Symmetry: "Quenching" g_A in Nuclei*, Symmetry 13(8), 1388 (2021). DOI 10.3390/sym13081388。署名单位:巴黎萨克雷大学 / CNRS / CEA 理论物理研究所(IPhT),91191 Gif-sur-Yvette。【已 DOI 直查;单作者】
  1. M. Rho, L.-Q. Shao(邵龙奇), *Probing for an IR-fixed Point in QCD by Superallowed Gamow-Teller Transitions in Doubly Magic Nuclei*, arXiv:2410.04991(2024-10-07 提交,v2 更新至 2026-09-21;注明 "Revised and prepared for invited paper")。截至 2026-09-21 未在 CrossRef 查到期刊版。 给出 q^GSI ≈ 0.73、q^RIKEN ≈ 0.50。署名单位:巴黎萨克雷大学 / CNRS / CEA 理论物理研究所(IPhT);国科大杭州高等研究院、中国科学院理论物理研究所、中国科学院大学。【只有 arXiv】
  1. L. Coraggio, N. Itaco, G. De Gregorio, A. Gargano, Z. H. Cheng, Y. Z. Ma, F. R. Xu(许甫荣), M. Viviani(8 人), *The renormalization of the shell-model Gamow-Teller operator starting from effective field theory for nuclear systems*, Physical Review C 109(1), 014301 (2024). DOI 10.1103/PhysRevC.109.014301。署名单位以 INFN 那不勒斯为主,含北京大学。【已 DOI 直查】
  1. B.-L. Wang(王滨磊), W.-L. Lv(吕万利), L.-G. Cao(曹李刚), Y.-F. Niu(牛一斐), G. Colò, H. Sagawa, F.-S. Zhang(张丰收)(7 人), *The quenching of the axial-vector coupling constant g_A in β-decay: joint effects from chiral two-body currents and many-body correlations*, arXiv:2605.21218(2026-05-20 提交)。截至 2026-09-21 未在 CrossRef 查到期刊版。 结论 q ≈ 0.73–0.80。署名单位:北京师范大学、兰州大学(含稀有同位素前沿科学中心)、北京市科学技术研究院辐射技术研究所、米兰大学 + INFN 米兰、理化学研究所仁科中心、会津大学、中国科学院理论物理研究所、兰州重离子加速器国家实验室。【只有 arXiv(附 HTML 版作者单位)】
  1. C. M. Ireland, G. Bollen, S. E. Campbell, X. Chen 等 26 人, *Mass of ¹⁰¹Sn and Bayesian extrapolations to the proton drip line*, Physical Review C 113(3), L021302 (2026). DOI 10.1103/vz5k-yf93(随机式,非按卷期页拼);arXiv:2510.11815。要点:FRIB 首次彭宁阱测 ¹⁰¹Sn,质量过剩 −59 889.89(96) keV,精度提升 300 倍;贝叶斯模型组合预测 ¹⁰⁰Sn 质量过剩 −58 090(800) keV,1σ 内偏向 GSI 的 β 延迟 Q 值。署名单位以密歇根州立大学 FRIB / LEBIT 为主,含阿贡国家实验室、圣母大学、埃尔朗根-纽伦堡大学(Reinhard)等。【已 DOI 直查(随机式 DOI)+ FRIB/FAU 机构库记录】
  1. J. Karthein, C. M. Ricketts, R. F. Garcia Ruiz 等 27 人, *Electromagnetic properties of indium isotopes illuminate the doubly magic character of ¹⁰⁰Sn*, Nature Physics 20(11), 1719–1725 (2024). DOI 10.1038/s41567-024-02612-y;arXiv:2310.15093。装置:CERN-ISOLDE 激光谱学。署名单位含麻省理工学院、曼彻斯特大学、CERN。【已 DOI 直查】
  1. M. Procura, B. U. Musch, T. R. Hemmert, W. Weise, *Chiral extrapolation of g_A with explicit Δ(1232) degrees of freedom*, Physical Review D 75(1), 014503 (2007). DOI 10.1103/PhysRevD.75.014503;arXiv:hep-lat/0610105。署名单位:慕尼黑工业大学理论物理系。结论:不含显式 Δ 的标准 πN 手征微扰论到 p⁴ 阶无法内插格点数据,加入 Δ 后成功。【已 DOI 直查】
  1. H. Zhou(周昊), L.-J. Wang(王龙军,通讯), Y. Sun(孙扬), *Many-body correlations as the origin of Gamow-Teller quenching in nuclear β-decay*, arXiv:2509.24542(v1 2025-09-29,v2 2026-04-01)。截至 2026-09-21 未在 CrossRef 查到期刊版。 结论:手征两体流仅贡献 5%–15%,多体关联是主因。署名单位:西南大学物理科学与技术学院(重庆)、上海交通大学物理与天文学院。【只有 arXiv(附 HTML 版作者单位)】
  1. J. T. Suhonen, *Value of the Axial-Vector Coupling Strength in β and ββ Decays: A Review*, Frontiers in Physics 5, 55 (2017). DOI 10.3389/fphy.2017.00055;arXiv:1712.01565。署名单位:于韦斯屈莱大学。【已 DOI 直查;单作者】
  1. J. Engel, J. Menéndez, *Status and future of nuclear matrix elements for neutrinoless double-beta decay: a review*, Reports on Progress in Physics 80(4), 046301 (2017). DOI 10.1088/1361-6633/aa5bc5;arXiv:1610.06548。【已 DOI 直查;2 作者】
  1. S. J. Novario, P. Gysbers, J. Engel, G. Hagen 等 9 人, *Coupled-Cluster Calculations of Neutrinoless Double-β Decay in ⁴⁸Ca*, Physical Review Letters 126(18), 182502 (2021). DOI 10.1103/PhysRevLett.126.182502。要点:淬灭因子由 ¹⁰⁰Sn 区域两体流的 Ikeda 求和规则结果推出,非人为选取。【已 DOI 直查;9 作者】
  1. A. Belley, J. M. Yao(姚建明), B. Bally, J. Pitcher, J. Engel 等 12 人, *Ab Initio Uncertainty Quantification of Neutrinoless Double-Beta Decay in ⁷⁶Ge*, Physical Review Letters 132(18), 182502 (2024). DOI 10.1103/PhysRevLett.132.182502;arXiv:2308.15634。署名单位含 TRIUMF、不列颠哥伦比亚大学、中山大学(珠海)、萨克雷 IRFU、北卡罗来纳大学。【已 DOI 直查】
  1. S. R. Stroberg, J. D. Holt, A. Schwenk, J. Simonis(4 人), *Ab Initio Limits of Atomic Nuclei*, Physical Review Letters 126(2), 022501 (2021). DOI 10.1103/PhysRevLett.126.022501。【已 DOI 直查】

---

---

核实与纠错记录

  1. Ejiri 2025 的 DOI 不能按卷期页拼。 期刊信息是 PRC 112(4), 045505 (2025),但真实 DOI 是 10.1103/gtwh-dqmm;按 10.1103/PhysRevC.112.045505 查询会直接 404(本次实测确认)。这是 APS 自 2025 年起对新文章改用随机式 DOI 的结果,2025 年前的老文章不受影响。[12]
  1. Ejiri 2025 的 arXiv v1 题目与期刊题名不同。 arXiv:2502.06122 v1 题为 *Delta-isobar resonance effects on nuclear matrix elements for neutrinoless ββ decays and astro-neutrino inverse-β decays*;期刊题名为 *Δ-isobar resonance effects studied by τσ summed strengths and nuclear matrix elements for β and ββ decays*。抄 v1 题目会引用到一个不存在的标题。[12]
  1. PRC 100, 054324 的作者顺序常被写错。 正确顺序是 L.-G. Cao、S.-S. Zhang、H. Sagawa(CrossRef 作者数组、APS 官网作者行、ADS 条目三处一致);二手材料常写成 "Sagawa, Cao, Zhang"。[11]
  1. 同一作者的单位会变,引用须按发表时。 曹李刚在 PRC 100, 054324 (2019) 的署名单位是华北电力大学数理学院;在 arXiv:2605.21218 (2026) 的署名单位已是北京师范大学与北京市科学技术研究院辐射技术研究所。同一人被写成两个单位不是矛盾,是时间差。[11][25]
  1. Hinke 2012 的 B(GT) 误差棒有两个版本。 Hinke 博士论文给出 B(GT) = 9.1 (+4.8/−2.3);后续文献(Lubos 等 2019、Rho & Shao 2024)普遍引用期刊版 9.1 (+3.0/−2.6)。中心值一致,误差棒不同,引用时须注明来源。若看到"9.1 正4.8负2.3"被当成 Nature 正文值引用,那就是把学位论文的数值挂到了期刊论文上。[19][20][23]
  1. Hinke 2012 的产核数是 259 个,不是 2500 个。 259 是 GSI 那次的数字(比此前实验多一个数量级,此前约 70 个衰变的量级);约 2500 个是 RIKEN/Lubos 2019 的数字。两个数字经常被互换。[19][20]
  1. ¹⁰⁰Sn 的两组半衰期几乎相同,分歧全在端点能量。 T₁/₂:GSI 1.16 ± 0.20 s;RIKEN 1.18 ± 0.08 s。β 端点能量:GSI 3.29 ± 0.20 MeV;RIKEN 3.91 ± 0.15 MeV。若只记"两组测出的半衰期差很多",就会把真正的问题(0.62 MeV 的端点差)漏掉。[19][20]
  1. Mougeot 2021 的 DOI 尾号是 01326-9,不是 01333-w。 标题讲的是 ⁹⁹⁻¹⁰¹In 的质量测量挑战 ¹⁰⁰Sn 的从头算理论。[18]
  1. Ireland 2026 的 DOI 也是随机式。 期刊信息 PRC 113(3), L021302 (2026),真实 DOI 10.1103/vz5k-yf93。同时注意:该工作标题虽以 ¹⁰¹Sn 为主,但核心结论之一是关于 ¹⁰⁰Sn 质量偏向 GSI 的 Q 值。[26]
  1. arXiv 页面没有 journal-ref 字段,不等于未发表。 本次通过 CrossRef 按题名直查确认:Rho & Shao (arXiv:2410.04991)、Zhou/Wang/Sun (arXiv:2509.24542)、Wang 等 (arXiv:2605.21218) 三篇确实尚无期刊版(而非漏挂);而 Baroni 等 (arXiv:1806.10245)、Karthein 等 (arXiv:2310.15093)、Yako 等 (arXiv:nucl-ex/0411011)、Arima 等 (arXiv:nucl-th/0008029) 则均已正式发表。两向都要查。[16][23][25][26][28]
  1. **"ECT"这个缩写后面常带一个星号(ECT*)**,在 markdown 渲染里会被当成斜体标记吃掉,导致整段文字排版错乱。本篇统一写作 "ECT(意大利 Trento)"。[17]
  1. ¹⁰⁰Sn 的 β 端点能量在 GSI 一方是"端点能量"、RIKEN 一方也写作"Q_β",两者在本篇按同一物理量比较(两篇论文均自称测的是 β 谱端点)。但请读者注意:Q_EC 与 β⁺ 端点能量差 2m_e c² ≈ 1.022 MeV,文献中三者常被混用而不加区分。比较不同论文的数字前,先确认它说的是哪一个。[19][20]
  1. ¹²⁸Ag 的中子数是 81,不是 82(它是最重的双幻核 ¹³²Sn 的"南邻")。该工作的主导单位是北京大学(一作罗迪雯,通讯李智焕、华辉,华东师范大学程奕源为共同通讯),中科院近代物理研究所未参与。这一点在二手综述里常被写错。
  1. PRC 111(3), 034311 (2025)(Heinz 等,6 人)是关于钙同位素从头算结构的工作(IMSRG(3)-N7 研究 ⁴⁴,⁴⁸,⁵²Ca),与 N=82 壳闭合无关。若见到它被引用在"锡同位素 N=82"的语境里,那是张冠李戴。
  1. Oset-Toki-Weise 1982 那篇 Physics Reports 的 DOI 是 10.1016/0370-1573(82)90123-5。 同期同卷同刊还有一个 …90134-X,是 Ghosh 与 Deb 关于"密度、密度泛函与电子流体"的完全无关的论文(Physics Reports 92, 1–44, 1982)。按印象填 DOI 极易串到那一篇上。[6]
  1. Arima 等 2001 那篇 PLB 的作者单位非常分散且含两位同名作者:Arima, A.(东京参议院);Bentz, W.(东海大学平塚);Suzuki, T.(日本大学);Suzuki, T.(福井大学 + 理研和光)。引用时只写 "T. Suzuki" 无法区分。[8]
  1. Baroni 等 2018 的 DOI 是 10.1103/PhysRevC.98.044003(PRC 98, 044003),不是 054003。本人在首轮检索时按印象写了 054003,CrossRef 直接 404,改为 044003 后命中。[16]
  1. "两体流贡献 5%–15%"这个数字,只出自 Zhou/Wang/Sun 预印本,且限定于 ⁷⁶Ge 的大组态壳模型框架。它不能被推广成"两体流在任何核里都只有 5%–15%"。[28]
  1. 本篇与知识库既有文章的分工。 既有《β衰变核结构:g_A 淬灭五十年谜题走到哪了》(20260822T131703)是全景综述;本篇是"根源判别",核心增量是:Δ-isobar 与两体流在 χEFT 记账法里的同一性论证(§3.3)、Ejiri 2025 的 0.7 与两体流的 0.73–0.85 的数值重合(§5.4)、以及 2021 年 ISOLTRAP 与 2026 年 FRIB 两次独立质量测量对 ¹⁰⁰Sn 分歧的裁决(§8)。既有文章称"两条测量还在打架",本篇给出了新的裁决证据与方向。[12][13][18][26][25]