Two-neutrino $ββ$ decay to excited states at next-to-leading order

该研究利用核壳模型和手征有效场论，计算了多种双贝塔衰变实验核素跃迁至第一激发态的核矩阵元及次领头阶修正，分析了核形变与组态结构对衰变半衰期的影响，并指出其预测值与当前实验限制基本一致。

要点

这篇论文就像是一群物理学家在**“原子核宇宙”里进行的一次精密的“双生子探险”**。

为了让你轻松理解，我们可以把这篇论文的核心内容想象成在研究一种极其罕见的**“原子核变身术”**。

1. 核心故事：原子核的“双重变身”

想象一下，原子核里住着两个调皮的中子兄弟。在一种极其罕见的情况下（叫做双贝塔衰变），这两个中子兄弟决定同时变身成两个质子兄弟。

• 变身过程：为了保持平衡，他们必须吐出两个电子和两个反中微子（就像变身时留下的“烟雾”）。
• 速度极慢：这个过程慢得令人发指，一个原子核要完成一次变身，平均需要等待1000 亿亿年（比宇宙年龄还长得多）。
• 为什么要研究它？
  • 科学家想通过观察这种变身，来寻找一种更神秘的“无中微子”变身（即没有吐出反中微子）。如果找到了，就能证明中微子就是它自己的反物质，这将彻底改变我们对宇宙起源（为什么物质多于反物质）的理解。
  • 为了找到那个“无中微子”的奇迹，我们必须先极其精确地算出普通“双中微子”变身的概率。

2. 这次探险的新发现：不仅看“地面”，还要看“二楼”

以前的研究主要关注原子核从“地面层”（基态）直接变到“地面层”的情况。但这篇论文把目光投向了**“二楼”**（激发态）。

• 比喻：以前我们只研究人从一楼跳回一楼；现在我们要研究人从一楼跳到二楼，再跳回一楼的过程。
• 挑战：跳上二楼的过程更复杂，中间涉及很多看不见的“中间状态”（就像跳高时的起跳点）。

3. 他们用了什么工具？（核壳模型与“调音师”）

为了计算这种变身有多难，作者们使用了一种叫做**“核壳模型”**的超级计算器。

• 比喻：这就像是用乐高积木搭建原子核的模型。但是，积木的拼接规则（哈密顿量）有很多种，不同的规则搭出来的形状略有不同。
• 问题：以前的计算发现，用“裸奔”的规则（不加修饰）算出来的变身概率太快了，和实验对不上。
• 解决方案：他们像**“调音师”**一样，给计算加了“阻尼”（淬火因子），或者换了一种更高级的“有效规则”，让计算结果能贴合已知的实验数据。

4. 两个关键的新发现

A. “次级效应”：虽然微小，但有时很关键

作者们引入了**“次领头阶（NLO）”**的修正。

• 比喻：如果你开车从 A 到 B，通常只算距离（主项）。但这次他们连**“风阻”、“轮胎摩擦”和“路面颠簸”**（次级项）都算进去了。
• 结果：在大多数情况下，这些额外的因素只让结果改变了5% 以内，就像你开车时忽略风阻，误差也不大。
• 例外：但是，如果主项（距离）因为某种原因互相抵消变得很小（就像两个力方向相反，合力几乎为零），那么这些“风阻”和“颠簸”就会变得非常重要，甚至主导结果。这就好比在走钢丝，稍微一点风就能把你吹倒。

B. “形状差异”：原子核长得越像，变身越容易

这是论文最精彩的部分。作者发现，原子核的**“长相”（形变）**对变身概率影响巨大。

• 比喻：想象两个原子核，一个是圆滚滚的球，另一个是橄榄球。
  • 如果变身前后的原子核形状非常相似（都是橄榄球），变身就比较容易，概率大（半衰期短）。
  • 如果形状差异很大（从球变成橄榄球），变身就很困难，概率小（半衰期长）。
• 三轴性：作者还特别研究了原子核是不是“歪”的（三轴形变）。就像鸡蛋如果歪着长，和正着长，变身难度完全不同。
• 结论：他们发现，不同计算模型预测的巨大差异，很大程度上是因为它们对原子核“长相”的预测不同。

5. 最终结论：我们能观测到吗？

• 对于 76Ge（锗）：计算出的变身时间非常长，但最乐观的预测已经接近目前实验能探测到的极限了。这意味着，如果我们再努力一点，或者实验精度再提高一点，有可能在不久的将来看到它。
• 对于 82Se（硒）：计算结果与最近的一个实验迹象吻合，这让人很兴奋。
• 对于其他元素：比如 136Xe（氙），预测的变身时间比目前实验能达到的极限还要长100 倍。这意味着在短期内，我们很难在这些元素里看到这种“跳二楼”的变身。

总结

这篇论文就像是在告诉我们要**“知己知彼”**：

1. 我们不仅要看原子核变身的“主菜”（主要过程），还要小心那些不起眼的“配菜”（次级修正），因为它们在某些情况下会翻盘。
2. 原子核的**“身材”**（形变）决定了它变身有多难。
3. 虽然理论预测还有很大的不确定性（就像用不同尺子量出来的长度不一样），但我们的计算正在逼近实验的边界。特别是对于锗（Ge）和硒（Se），未来的实验极有可能验证这些理论，甚至揭开中微子神秘面纱的一角。

简单来说，这是一次高精度的理论“排雷”行动，告诉实验物理学家：“别在那些太远的地方浪费时间，把重点放在锗和硒上，那里最有希望发现新大陆！”

技术摘要

这是一篇关于双中微子双贝塔衰变（$2\nu\beta\beta$）到激发态的理论物理研究论文。作者利用核壳模型（NSM）和手征有效场论（$\chi$EFT），对多种原子核的衰变半衰期进行了次领头阶（NLO）精度的计算。

以下是该论文的详细技术总结：

1. 研究背景与问题 (Problem)

• 物理过程：$2\nu\beta\beta$ 衰变是两个中子同时转化为两个质子，并发射两个电子和两个反中微子的过程。它是已知最慢的衰变过程之一，半衰期可超过 $10^{21}$ 年。
• 科学意义：该过程与无中微子双贝塔衰变（$0\nu\beta\beta$）密切相关。理解 $2\nu\beta\beta$ 的核矩阵元（NMEs）对于提高 $0\nu\beta\beta$ 的预测精度至关重要，后者是验证中微子是否为马约拉纳粒子及探索宇宙物质 - 反物质不对称性的关键。
• 现有挑战：
  • 大多数多体方法（如 NSM、QRPA、IBM）在使用裸算符时往往会高估 $2\nu\beta\beta$ 的 NMEs，通常需要通过唯象“淬灭”（quenching）因子来修正。
  • 现有的计算主要集中在基态到基态（$0^+_{gs} \to 0^+_{gs}$）的跃迁，而针对基态到第一激发态（$0^+_{gs} \to 0^+_2$） 的衰变研究较少，且缺乏高精度的 NLO 修正。
  • 核形变（特别是三轴形变）对 NMEs 的影响尚需深入分析。

2. 方法论 (Methodology)

作者采用**核壳模型（NSM）作为主要计算框架，并结合手征有效场论（$\chi$EFT）**引入高阶修正。

• 核模型与哈密顿量：
  • 针对 $^{48}\text{Ca}$、$^{76}\text{Ge}$、$^{82}\text{Se}$、$^{124}\text{Sn}$、$^{130}\text{Te}$ 和 $^{136}\text{Xe}$ 等核素，使用了多种经过验证的壳模型哈密顿量（如 KB3G, GXPF1A, GCN2850, JUN45, RG, GCN5082, QX 等），以确保能准确描述初态和末态的能谱。
  • 使用了 ANTOINE 和 NATHAN 代码进行计算。
• 算符处理：
  • 比较了两种 Gamow-Teller (GT) 算符：
    1. 裸 GT 算符：配合唯象淬灭因子（$q_{\beta\beta}$ 或 $q_\beta$）以重现实验测得的基态到基态半衰期。
    2. 有效 GT 算符：在壳模型价空间内重整化，理论上已包含部分多体关联，通常不需要额外淬灭（或淬灭较小）。
• 高阶修正（NLO）：
  • 轻子能量展开：对轻子能量分母进行泰勒展开，保留前三项，引入高阶相空间因子（PSF）和次领头阶 NMEs。
  • $\chi$EFT 修正：首次在这些衰变中评估了 NLO 长程项，包括弱磁矩（weak magnetism）和单π介子交换图（one-pion-exchange）。计算中保留了能量分母展开的三项。
• 核结构分析：
  • 利用形状不变量（Shape invariants）分析初态和末态的形变参数（$\beta_2, \gamma$），特别是**三轴形变（triaxiality）**的影响。
  • 分析了态的**同位旋（seniority）**结构对 NMEs 的影响。

3. 主要结果 (Key Results)

A. 半衰期预测

• 总体量级：预测的 $0^+_{gs} \to 0^+_2$ 半衰期通常比当前的实验下限长至少两个数量级，意味着近期在这些核素中探测到该过程非常困难。
• 例外情况：
  • $^{76}\text{Ge}$：预测的半衰期下限接近当前的实验限制（约 $10^{24}$ 年量级），但不同哈密顿量导致的预测范围跨度达两个数量级（$2-180 \times 10^{24}$ 年）。
  • $^{82}\text{Se}$：预测范围（$0.4-190 \times 10^{23}$ 年）与最近的实验迹象（Indication）一致。
• 不确定性来源：主要的不确定性来源于核哈密顿量的选择，而非算符形式（裸算符与有效算符的结果通常较为接近）。

B. 高阶修正贡献

• 轻子能量展开（$\epsilon_{\text{Taylor}}$）：通常贡献在百分之几（~1-12%）。但在 $^{82}\text{Se}$ 和 $^{136}\text{Xe}$ 中，由于领头阶 NME 发生显著抵消，这些项变得更为重要。
• NLO 项（$\epsilon_{\text{NLO}}$）：
  • 在大多数情况下，NLO 贡献小于 5%。
  • 显著增强效应：当领头阶 GT NME 因不同组分的抵消而变得非常小（甚至接近零）时，NLO 项（特别是单π介子交换项）的贡献会显著增加，甚至主导总结果。例如在 $^{124}\text{Sn}$ 和 $^{136}\text{Xe}$ 中，NLO 贡献可达 10-30%，甚至在特定相互作用下超过领头阶项。

C. 核结构与形变的影响

• 形变差异与 NME：发现初态和末态之间的**形变差异（$\delta_{\text{def}}$）**与 NME 大小呈线性关系。形变越相似，NME 越大；形变差异越大，NME 越小。
• 三轴形变的关键作用：忽略三轴形变（$\gamma$）会破坏 NME 与形变差异的线性相关性。例如在 $^{76}\text{Ge}$ 中，不同哈密顿量预测的 $\gamma$ 值不同，导致 NME 差异巨大。
• 同位旋结构（Seniority）：在 $^{82}\text{Se}$ 中，JJ45 哈密顿量预测的 NME 较大，是因为在该哈密顿量下，初末态不同同位旋分量间的典型抵消效应未发生。这表明除了形变，态的同位旋结构也是决定 NME 的关键因素。
• 哈密顿量选择：对于 $^{76}\text{Ge}$，能最好重现实验谱四极矩（$Q_s$）的 RG 相互作用预测了最长的半衰期；对于 $^{82}\text{Se}$，JUN45 和 JJ4BB 与实验数据符合较好，但预测的半衰期差异巨大。

4. 关键贡献 (Key Contributions)

1. 首次 NLO 计算：首次在手征有效场论框架下，对 $2\nu\beta\beta$ 到激发态的衰变进行了次领头阶（NLO）计算，量化了弱磁矩和π介子交换的贡献。
2. 激发态系统研究：系统计算了多种实验候选核素（$^{48}\text{Ca}, ^{76}\text{Ge}, ^{82}\text{Se}, ^{124}\text{Sn}, ^{130}\text{Te}, ^{136}\text{Xe}$）的 $0^+_{gs} \to 0^+_2$ 衰变半衰期。
3. 形变与结构的关联分析：深入揭示了核形变（特别是三轴形变）和同位旋结构对 NMEs 的调控机制，解释了不同壳模型哈密顿量预测结果差异巨大的物理根源。
4. 实验对比：将理论预测与最新的实验数据（特别是 $^{82}\text{Se}$ 的实验迹象和 $^{76}\text{Ge}$ 的极限）进行了详细对比，指出了理论的不确定性范围。

5. 意义与展望 (Significance)

• 理论精度提升：通过引入 NLO 修正和更精细的核结构分析，提高了对 $2\nu\beta\beta$ 衰变理论预测的可靠性，减少了模型依赖性。
• 实验指导：明确了 $^{76}\text{Ge}$ 和 $^{82}\text{Se}$ 是探测 $0^+_{gs} \to 0^+_2$ 衰变最有希望的候选者。未来的实验若能确认 $^{82}\text{Se}$ 的信号或在 $^{76}\text{Ge}$ 中观测到该过程，将是对当前核结构模型和 NME 计算的重要验证。
• 对 $0\nu\beta\beta$ 的启示：由于 $2\nu\beta\beta$ 和 $0\nu\beta\beta$ 的 NMEs 存在相关性，更精确的 $2\nu\beta\beta$ 计算有助于约束 $0\nu\beta\beta$ 的核矩阵元，从而更准确地推断中微子质量。
• 未来工作：作者计划在未来工作中包含 NLO 短程项（尽管耦合常数未知）以及显式的两核子流（two-nucleon currents），并进一步探索核结构与 NME 的深层联系。

总结：该论文通过高精度的壳模型计算和有效场论修正，系统研究了双贝塔衰变到激发态的物理机制，揭示了核形变和同位旋结构在决定衰变率中的核心作用，并为未来的实验探测提供了重要的理论基准。