Theoretical analysis and predictions for the double electron capture of $^{124}$Xe

本文通过改进核结构与原子结构计算，对$^{124}$Xe中的双中微子双电子捕获过程进行了全面的理论分析，由此得到了精炼的核矩阵元，预测了各种衰变道的特定俘获份额（显著地，KK道为74%，累积的KL$_1$-KO$_1$道为24%），并为液氙实验中的本底建模提供了更新后的原子弛豫能。

要点

想象一下，原子是一个繁忙的小型城市。在这个城市内部，原子核就像是市议会，挤满了质子和中子。通常情况下，这些居民非常稳定，但有时，他们会决定重新排列组合，以便让自己感觉更舒适。

这篇论文是关于发生在名为 氙-124 (Xenon-124) 的城市中一次非常罕见且特定的“重组”事件。在这种事件中，市议会决定从外部社区抓取两个属于自己的居民（电子）并把它们拉入原子核内。当这种情况发生时，这座城市会转变为一座名为 � tellurium-124 (獌-124/Tellurium-124) 的新城市，并吐出两个微小的、隐形的信使——中微子。

科学家们称之为双电子捕获（具体指两中微子版本，即 $2\nu$ECEC）。这就像是在泳池里进行“双重浸入”，只不过浸入的不是水，而是亚原子粒子。

以下是研究人员的工作内容，通过简单的语言进行了解释：

1. 构建更好的蓝图（理论部分）

在过去，科学家试图预测这种“双重浸入”发生的频率，但他们的蓝图有些粗糙。他们忽略了关于电子如何移动以及原子核如何反应的一些细节。

本文作者决定构建一个更加精确的蓝图。

• “泰勒展开”类比： 想象你在描述一辆车的路径。一个简单的描述可能只是说“它向前行驶”。一个更好的描述会加上“它在加速”。而最好的描述会加上“它在加速，然后轻微转向，然后减速”。作者使用了名为“泰勒展开”的数学工具来增加这些额外的细节层级（高达能量的四次方）。这使得他们能够观察到之前模型所遗漏的衰变过程中的“转向与减速”。
• “新比例”： 由于增加了这些额外细节，他们发现了比较该过程不同部分的新方法（称为 $\xi$ 比例）。你可以将这些视为赛车场上的新检查点，有助于科学家在日后更好地理解这场比赛。

2. 观察周边社区（原子部分）

为了计算这种事件发生的可能性，你需要准确知道电子居住在哪里。

• “泡利阻塞”隐喻： 想象一个拥挤的电梯。如果电梯已经满了，你就不能直接把别人挤进去；你必须等有人出来。在原子核内部，“最内层”的位置就像是一个满员的电梯。作者意识到，被捕获的电子不能随心所欲地进入任何地方；它们会被已经存在的其他电子所阻挡。他们考虑了这种“拥挤”规则，这改变了计算结果。
• 扩大搜索范围： 之前的研究只关注紧邻原子核的两个社区（称为 K 和 L1 层）。作者说：“让我们观察所有的社区，甚至是那些更远的社区（直到 O 层）。”他们发现，虽然外围社区被捕获的可能性较低，但它们仍然对总事件有所贡献。

3. 模拟市议会（原子核部分）

模拟原子核是最困难的部分，因为它是粒子组成的混乱人群。作者使用了两种不同的“模拟引擎”来预测原子核的行为：

• 引擎 A (ISM)： 这就像是一个详细的、逐房间进行的模拟。他们使用不同的规则（称为“哈密顿量”）运行模拟，以观察结果是否成立。他们发现，当他们包含了原子核采取的所有可能的“中间”步骤时，预测的事件“强度”比以前更简单的模型要低。
• 引擎 B (pn-QRPA)： 这是另一种类型的模拟。他们调整了这个引擎的设置，直到它与我们已有的现实世界数据相匹配。他们发现，使用这种引擎进行更细致的计算后，得到的“强度”值比之前的尝试要小得多。

4. 结果：他们发现了什么？

通过结合更好的蓝图、详细的社区地图和这两个模拟引擎，他们做出了几项预测：

• 主要事件 (KK 通道)： 他们预测，大约 74% 的情况下，两个电子会被从最靠近的社区（K 层）捕获。这与之前实验中使用的 72.4% 略有不同，是一个微小但重要的调整。
• “次优”事件： 他们预测，大约 19% 的情况下，一个电子来自最近的社区，而另一个来自紧随其后的社区 (KL1)。
• “累积”预测： 如果把所有频率稍低的事件（从 KL1 到 KO1）相加，它们约占总事件的 24%。这大约是主要事件的三分之一。
• “弛豫”能量： 当电子被捕获时，新的城市（獌/Tellurium）处于激发态并需要冷静下来。它通过释放能量（如 X 射线）来实现这一点。作者计算了每种捕获类型释放的精确能量。这就像是为科学家提供了一个可以用来在探测器中寻找的特定“指纹”。

为什么这很重要？

这篇论文并不声称它能治愈疾病或为城市提供动力。相反，它扮演着探索者精炼地图的角色。

使用液态氙的大型实验（例如那些寻找暗物质的实验）一直在密切监测这种特定的“双重浸入”事件。然而，这种事件看起来可能像“背景噪声”，从而干扰数据。通过提供关于这种事件发生的频率、释放多少能量以及电子来自哪些“社区”的更准确地图，作者帮助实验人员将真实的信号与背景噪声区分开来。

简而言之，他们将一张模糊、低分辨率的罕见原子事件照片，转化为了高清晰度、3D 模型的图像，帮助科学家明确知道在探测器中应该寻找什么。

技术摘要

技术摘要：$^{124}$Xe 双电子俘获过程的理论分析与预测

问题陈述
双电子俘获（ECEC）过程，特别是两中微子模式（$2\nu$ECEC），是一种对于探测核结构和中微子性质具有重要意义的稀有核跃迁。虽然 $2\nu\beta^-\beta^-$ 衰变已得到广泛研究，但由于跃迁概率较低，ECEC 模式受到的关注较少。然而，$^{124}$Xe 是一个独特的案例，因为它在能够发生 ECEC 跃迁的同位素中拥有最大的 Q 值，使其成为大规模液氙实验（例如寻找暗物质和无中微子双 $\beta$ 衰变的实验）的主要目标。

以往对 $^{124}$Xe $2\nu$ECEC 半衰期的理论计算在核物理和原子模型方面都存在局限性。具体而言，早期的工作通常依赖于核矩阵元（NMEs）的闭合近似（closure approximation），仅考虑 K 和 L$_1$ 壳层的俘获，并使用了不够严谨的原子屏蔽和泡利阻塞处理。随着 $^{124}$Xe 半衰期实验测量值的获得，需要一个更严谨且完整的理论描述来解释这些数据并预测未观测到的衰变通道。

方法论
作者提供了一个综合性的理论框架，改进了针对 $^{124}$Xe 的形式体系及具体计算：

1. 改进的形式体系： 使用泰勒展开法导出衰变率，保留至电子能量的四次方项。该方法引入了高阶修正，并定义了新的 NME 比率 $\xi_{31}^{2\nu\text{ECEC}}$ 和 $\xi_{51}^{2\nu\text{ECEC}}$，这两个比率取决于中间 $1^+$ 态的能量分母。
2. 核结构计算：
   • 相互作用壳模型（ISM）： 计算在 $jj55$模型空间（轨道包括$0g_{7/2}, 1d_{5/2}, 1d_{3/2}, 2s_{1/2}, 0h_{11/2}$）内进行。不同于以往研究假设单态主导（SSD）的方法，本研究通过强度函数法实现了对中间 $1^+$ 态的全求和。研究采用了两种有效哈密顿量（GCN5082 和 SVD），其淬灭因子（$q_H$）经过校准以重现 $^{136}$Xe 的 $2\nu\beta^-\beta^-$ 数据。计算中包含了允许多达四个核子从低能轨道激发的截断。
   • 质子-中子准粒子随机相位近似（pn-QRPA）： 计算利用了同位旋恢复。质子-中子相互作用强度（$g_{pp}^{T=0}$）经过调整以重现 $^{124}$Xe 的实验半衰期。研究使用了现实的核-核势（Argonne V18 和 CD-Bonn），并考虑了 BCS 真空重叠。
3. 原子结构与相空间因子（PSFs）：
   • 使用 Dirac-Hartree-Fock-Slater (DHFS) 自洽框架生成电子波函数。
   • 改进之处包括：严格处理原子屏蔽、弥散核表面修正、现实的核电荷密度以及电子交换相关效应。
   • 至关重要的是，计算考虑了最内层核子的泡利阻塞，并将电子俘获求和扩展到了所有可用的 s 波电子（直至 O$_1$ 子壳层），而非像以往研究那样仅限于 K 和 L$_1$ 壳层。
   • 计算了最终 $^{124}$Te 原子的原子弛豫能，以辅助背景建模。

关键结果

• 核矩阵元（NMEs）：
  • ISM 计算得到的 NME 对截断和有效哈密顿量的变化具有鲁棒性。由于采用了全求和而非 SSD，GCN5082 哈密顿量的 NME 相比截断方法降低了约 25%。
  • 当受到实验半衰期约束时，pn-QRPA 的结果显示其 NME 值显著低于之前的 pn-QRPA 计算（此前计算通常采用不同的 $g_{pp}$ 调整策略）。新的 $^{124}$Xe NME 被发现与 $^{128,130}$Te 的 $2\nu\beta^-\beta^-$ 衰变相当，反映了相似的配对特征。
• 半衰期预测：
  • 使用 ISM，计算出的 $^{124}$Xe 总半衰期在 $(1.1 \pm 0.2_{\text{stat}} \pm 0.1_{\text{sys}}) \times 10^{22}$ 年的实验值误差两倍以内。
  • 引入泰勒展开项和精炼后的 PSF 有助于使理论预测更接近实验数据。
• 俘获份额与通道：
  • 预测的主导 KK 通道的俘获份额为 74.1%，这与之前实验分析中使用的 72.4% 不同。这种差异源于本理论模型排他性地考虑了 s 波电子，而实验信号模型则包含了所有电子。
  • 次高概率的通道是 KL$_1$，其份额约为 19%。
  • KL$_1$ 到 KO$_1$ 通道的累积俘获份额预计约为 24%。
• 原子弛豫能：
  • 计算了各种空穴构型（如 KK 为 64.62 keV，KL$_1$ 为 37.05 keV）的精确弛豫能。这些数值旨在用于液氙实验的背景建模。

意义与主张
本文声称提供了第一个完整的 $^{124}$Xe $2\nu$ECEC 理论描述，该描述同时改进了核结构处理（通过全求和与同位旋恢复）和原子物理处理（通过 DHFS、泡利阻塞及扩展轨道俘获）。

作者断言，其 ISM 结果与实验数据一致，验证了所计算 NME 及新预测俘获份额的可靠性。其主要贡献在于预测了 KL$_1$-KO$_1$ 通道（代表约 24% 的衰变）应在 31.93–37.05 keV 的能量区间内被实验观测到。此外，该工作表明，部分半衰期之比（例如 $T_{KK}^{1/2}/T_{KL_1}^{1/2}$）可以作为对 $\xi_{31}^{2\nu\text{ECEC}}$ 参数的实验约束，从而为探测低能态与高能态之间的相互作用提供新途径。研究强调，以往的半衰期理论估计由于使用了闭合近似和不完整的轨道求和，往往与最精确的 PSF 不兼容。