Large-scale shell-model investigation of $2\nu$ECEC in $^{132}$Ba and $^{78}$Kr

本文通过对 $^{132}$Ba 和 $^{78}$Kr 中双中微子双电子俘获过程进行大规模壳模型研究，提供了基于经过验证的有效相互作用的更新后的核矩阵元和半衰期预测，旨在为未来的实验工作提供支持。

要点

以下是该论文的通俗化解释，采用了日常生活的类比。

大局观：捕捉隐形的幽灵

想象原子核是一个微小且拥挤的舞池。通常情况下，舞者（质子和中子）非常稳定，不会更换舞伴。但有时，两名舞者会决定在同一时刻交换位置。这是一个被称为**双电子捕获（Double Electron Capture）**的罕见事件。

在这种特定的“舞蹈”中，原子核内的两个质子会抓取来自原子外壳的两个电子，并转化为中子。因为这种现象极其罕见，所以需要极其漫长的时间——数万亿年——才能看到它发生一次。科学家们想要弄清楚究竟要等多久（即半衰期），因为这有助于他们理解宇宙的基本规则，例如中微子（一种微小的、幽灵般的粒子）的本质。

这篇论文的作者就像是建筑师和工程师。他们并没有建造一台新机器来捕捉这些事件；相反，他们构建了一个超级详细的计算机模拟系统，用以预测舞池的行为方式以及等待时间应该是多久。

两大主角：$^{132}\text{Ba}$ 和 $^{78}\text{Kr}$

研究人员专注于两种特定的原子核（候选者），它们是进行这种罕见舞蹈的潜在对象：

1. 钡-132 ($^{132}\text{Ba}$)：一种重原子，科学家怀疑它可能会跳这种舞，但目前还没有人亲眼目睹它在行动。我们仅根据旧有的地质线索知道它可能会发生。
2. 氪-78 ($^{78}\text{Kr}$)：一种科学家最近已确认其确实会跳这种舞的原子，但目前的测量结果仍有些模糊。

他们是如何做的：“乐高”模拟法

为了预测会发生什么，科学家们使用了一种称为**大规模壳模型（Large-Scale Shell Model）**的方法。

• 类比： 想象你要预测一个由数十亿块乐高积木组成的复杂结构如何支撑自身。你不能仅仅靠猜测；你需要准确知道每一块积木是如何与其邻居连接的。
• 工具： 科学家们使用了一套庞大的数字“乐高套装”（称为有效相互作用），它告诉计算机质子和中子如何相互作用。
  • 对于钡-132，他们使用了名为 SN100PN 的套装。
  • 对于氪-78，他们使用了 GWBXG。

升级版： 在之前对氪的研究中，他们只观察了乐高建筑的“底层”。而在这项新研究中，他们扩展了模型，将“高层”（高能级）也纳入其中。这就像是意识到，要了解一座摩天大楼如何在风中摇摆，你必须观察顶层，而不仅仅是地基。

检查蓝图：模拟成功了吗？

在信任关于这种罕见舞蹈的预测之前，科学家必须确保他们的模拟是准确的。他们通过检查相关原子的“正常”行为来做到这一点：

• 能级： 他们检查了计算机是否预测了原子正确的“振动”或能量状态。
• 形状： 他们检查了原子是呈球形还是略微扁平的蛋形（形变）。

结果： 计算机模拟与真实的实验数据几乎完美契合。这就像是建造了一个桥梁的缩放模型，并发现它的表现与真实的桥梁完全一致。这让他们对预测这种罕见舞蹈的可靠性充满了信心。

主要发现：关于“等待时间”的预测

1. 对于 钡-132（神秘候选者）

由于还没有人见过钡-132跳这种舞，科学家们提供了一个理论基准。

• 预测： 他们计算出，如果你等待大约 $7.33 \times 10^{24}$ 年（那是 7 后面跟着 24 个零！），你或许能看到它发生。
• 意义所在： 这是一个针对未来实验的“目标”。它告诉科学家：“不要指望在 100 年内看到它；你需要建造能够等待数万亿年的探测器。”他们的计算结果远高于目前科学家设定的最低限度，这意味着搜索工作仍处于开放状态。

2. 对于 氪-78（已确认的候选者）

科学家已经观察到氪-78跳这种舞，但测量值各不相同。

• 预测： 这个更详细的新模拟预测的等待时间为 $8.78 \times 10^{22}$ 年。
• 改进之处： 在他们之前的研究中（使用较小的乐高套装），他们预测的时间略有不同。通过在模型中加入“高层”，他们的新预测与近期实际观测到的结果更加接近。这就像是从模糊的照片升级到了高清图像；画面现在变得更清晰、更准确了。

“音量旋钮”（轴向耦合）

模拟中一个棘手的部分是，计算机并不知道宇宙中每一个微小的力。为了解决这个问题，科学家使用了一个被称为有效轴向耦合常数（$g_{\text{eff}}^A$）的“音量旋钮”。

• 类比： 想象你在录制一首歌，但你的麦克风漏掉了一些高音。你调大音量（旋钮）来补偿麦克风漏掉的部分。
• 科学家测试了不同的“音量设置”，以观察它如何改变预测的等待时间。即使在不同的设置下，他们的结果仍然与我们目前所知的知识保持一致。

结论：他们学到了什么？

论文的结论是：

1. 模拟非常可靠： 他们的计算机模型非常擅长描述这些原子的行为。
2. 钡-132： 他们提供了目前最好的理论推测，说明我们要等待多久才能看到它衰变。这有助于实验人员了解他们的探测器需要具备多高的灵敏度。
3. 氪-78： 通过使用更大、更复杂的模型，他们改进了预测，使其比以前更符合现实世界的数据。

简而言之，这些科学家构建了一张更好的原子舞池地图。他们还没有捕捉到舞者（对于钡而言），但他们对何时何地寻找有了更清晰的认识；而对于氪，他们的地图现在比旧版要精确得多。

技术摘要

技术摘要：关于 $^{132}$Ba 和 $^{78}$Kr 中 2νECEC 的大规模壳模型研究

问题陈述
对二阶弱相互作用过程（特别是两中微子双电子捕获，2νECEC）的研究，对于探测基本中微子性质和测试核结构模型至关重要。虽然两中微子双 $\beta$ 衰变（$2\nu\beta\beta$）已得到充分研究，但由于 2νECEC 模式的半衰期通常较长且 Q 值较低，其研究相对较少。可靠的核矩阵元（NME）和半衰期理论预测，对于指导未来的实验搜索至关重要，特别是针对像 $^{132}$Ba 这样目前仅存在半衰期下限的候选核素，以及最近观察到 $^{78}$Kr 的情况。现有的理论估算表现出显著的模型依赖性，因此需要稳健且经过独立验证的方法来约束这些稀有衰变模式。

方法论
作者采用大规模核壳模型（NSM）计算来研究 $^{132}$Ba 和 $^{78}$Kr 中的 2νECEC 过程。

• 有效相互作用： 研究对 $^{132}$Ba 使用了 SN100PN 有效相互作用，对 $^{78}$Kr 使用了 GWBXG 相互作用。
• 模型空间：
  • 对于 $^{132}$Ba（母核），模型空间包括质子和中子的 $0g_{7/2}1d_{2s}0h_{11/2}$ 轨道，并以 $^{100}$Sn 为惰性核芯。应用了截断以确保计算可行性，将构型限制在特定的粒子-空穴极限内。
  • 对于 $^{78}$Kr，本研究通过扩展中子模型空间扩展了前人的工作。早期的计算受限于 $N=50$ 壳层闭合，而当前的工作包含了跨越 $N=50$ 进入 $0g_{7/2}$、$1d_{5/2}$、$1d_{3/2}$ 和 $2s_{1/2}$ 轨道的单粒子-单空穴（1p-1h）激发。
• 计算工具： 使用 NUSHELLX 代码获取单体跃迁密度（OBTDs），使用 KSHELL 代码计算能谱和跃迁概率。
• 验证： 通过将计算出的光谱性质（低能能谱和约化四极矩跃迁概率 $B(E2)$）与实验数据进行对比，首先验证了有效相互作用的可靠性。
• NME 计算： 使用标准的 Gamow-Teller (GT) 单体算符计算 GT NME。分析了 NME 相对于中间核（$^{132}$Cs 和 $^{78}$Br）中 $1^+$ 态能量的累积贡献，以确保收敛。应用了有效轴矢量耦合常数（$g^A_{eff}$）以解释缺失的多体相关性和二体电流。

主要贡献

1. 首次实现 $^{132}$Ba 的壳模型 NME： 本工作展示了首次针对 $^{132}$Ba 中 2νECEC 过程进行的壳模型 NME 计算。
2. 改进了 $^{78}$Kr 的模型空间： 研究通过将中子模型空间扩展到 $N=50$ 壳层闭合之外，提供了更新且改进的 $^{78}$Kr 壳模型估算，解决了此前构型受限研究中的局限性。
3. 光谱验证： 作者证明了所采用的有效相互作用成功重现了参与衰变链的核素的实验光谱性质（能级和 $B(E2)$ 强度），为 NME 计算建立了可靠的基础。
4. 敏感性分析： 针对 $^{132}$Cs 中最低 $1^+$ 态未确认能量的情况，对 $^{132}$Ba 的计算进行了敏感性分析。

结果

• 光谱： 计算结果显示与 $^{132}$Ba、$^{132}$Cs、$^{132}$Xe、$^{78}$Kr、$^{78}$Br 和 $^{78}$Se 的实验数据高度吻合。对于 $^{132}$Ba，结果表明其具有中度长椭球形形变（$\beta_2 \approx +0.18$），而 $^{78}$Kr 则表现出具有形状混合特性的扁椭球型内在形状。
• $^{132}$Ba NME 与半衰期：
  • 累积 NME 在中间核约 6.3 MeV 处饱和，值为 0.0654。
  • 使用 $g^A_{eff} = 0.94$，预测的半衰期为 $7.33 \times 10^{24}$ 年。这比目前的实验下限（$> 2.2 \times 10^{21}$ 年）高出约三个数量级。
  • 只要 $^{132}$Cs 的最低 $1^+$ 态保持在 100 keV 以下，该计算就是稳健的；只有当该态位于 500 keV 以上时，才会出现 NME 的显著抑制。
• $^{78}$Kr NME 与半衰期：
  • 累积 NME 在饱和后达到最终值 0.2311，相比之前受限模型空间得到的 0.1997 有所提高。
  • 使用 $g^A_{eff} = 0.8$，预测的半衰期为 $8.78 \times 10^{22}$ 年。
  • 该结果与最近测得的实验值 $[1.9^{+1.3}_{-0.7} \pm 0.3] \times 10^{22}$ 年基本一致。
  • 研究还估算了竞争性的正电子发射通道（$2\nu\beta^+\beta^+$ 和 $2\nu EC\beta^+$）的半衰期，分别为 $4.68 \times 10^{26}$ 年和 $1.19 \times 10^{23}$ 年。

意义与主张
作者声称，这项工作为 $^{132}$Ba 提供了一个基准理论预测，以协助未来的实验努力来约束这一稀有衰变模式。对于 $^{78}$Kr，该研究提供了相对于早期更受限模型空间的更新且改进的壳模型估算。结果被认为在典型的理论不确定性范围内是合理的，并承认使用有效轴矢量耦合部分吸收了未被标准单体算符显式包含的缺失多体效应（如二体电流）。论文强调，虽然光谱的一致性验证了波函数，但衰变算符仍需通过 $g^A_{eff}$ 进行重整化。作者得出结论，其研究结果支持了用于评估 2νECEC NMEs 的底层核波函数的可靠性，尽管他们也指出，未来通过更大的模型空间和改进的同位旋对称性相互作用可以进一步完善这些预测。