Shake-up and shake-off spectra in the electron capture decay of atomic $^7$Be

本研究采用多构型狄拉克-福克计算来模拟原子 $^7$Be 在电子俘获衰变中的电子激发（shake-up）与电离（shake-off）光谱，揭示了尽管这些模型能够解释部分光谱特征，但材料诱导的波函数修正仍是一个挑战，并提供了一个修正后的 L/K 电子俘获比值 0.0756(20)，从而改善了对亚 MeV 级惰性中微子的约束。

要点

大局观：猎捕幽灵粒子

想象一下，宇宙是一个巨大的拼图，科学家们拥有一张描述其运作方式的图景，称为“标准模型”。但其中缺少了一些碎片。最大的谜团之一是暗物质，以及为什么物质比反物质更多。

为了寻找这些缺失的碎片，科学家们正在搜寻一种被称为**惰性中微子（sterile neutrino）**的“幽灵粒子”。这些粒子是隐形的、有质量的，它们不与普通物质发生相互作用，这使得捕捉它们变得异常困难。

BeEST实验是目前设置的最灵敏的陷阱之一，旨在捕捉这些幽灵。它使用了一种名为**铍-7（Beryllium-7, 7Be）**的放射性原子。当这种原子衰变时，它通常会吐出一个中微子并转变为锂原子。通过测量锂原子受到的微小“踢力”（反冲），科学家可以计算出中微子的质量。如果中微子很重（比如是惰性中微子），那么受到的踢力会比预期的小。

问题所在：“摇晃”效应

本文研究的是该实验中一个主要的干扰来源：电子激发（Shake-up）与电离（Shake-off）。

把原子想象成一栋房子，里面的家具（电子）按特定的房间（能层）排列。

1. 事件发生： 突然间，房子的主人（原子核）发生了变化。一个电子被捕获了，房子瞬间变成了另一种类型的房子（从铍变成了锂）。
2. 冲击： 因为房子变化得太快，家具并不会静止不动。它们会被震动。
   • 激发（Shake-up）： 一些家具被撞到了更高的架子上（激发态）。
   • 电离（Shake-off）： 一些家具被直接扔出了窗外（电离）。

过去，科学家使用粗略、陈旧的地图来预测家具会如何摇晃。这些地图就像是“卡通画”——它们没有考虑到家具之间互相碰撞（电子相关性）的情况，也没有考虑到高能物理（相对论）的影响。因为这些地图不准确，实验中的“背景噪声”非常混乱，导致很难发现幽灵粒子的信号。

本文的研究内容：高清晰度改造

本文作者决定从头开始构建一个关于这种摇晃过程的3D高清晰度模拟。

• 工具： 他们使用了一种极其先进的数学方法——多构型狄拉克-福克方法（Multiconfiguration Dirac-Fock）。你可以把它想象成一个物理引擎，能够模拟每一个电子如何与其他电子碰撞，同时兼顾爱因斯坦的相对论速度限制。
• 计算： 他们精确计算了电子被“向上摇晃”到更高架子或被“摇出房子”的概率，涵盖了“K壳层”（内层房间）和“L壳层”（外层房间）的捕获情况。
• 结果： 他们发现，这种摇晃比之前认为的更加剧烈且复杂。具体来说，当原子捕获来自外层“L壳层”的电子时，剩余电子的摇晃程度比捕获内层“K壳层”电子时要剧烈得多。

“钽（Ta）”因子：为什么模拟并不完美

论文提出了一个关键的区别：他们的完美模拟是针对在真空空间中漂浮的孤立原子进行的。然而，在实际实验中，铍原子被嵌入在**钽（Tantalum, Ta）**金属（传感器）之中。

• 类比： 想象你在模拟一个在真空中敲击鼓声的效果，但实际上你是在一个嘈杂拥挤的地铁站里敲击。金属传感器的壁面会改变电子的行为。
• 差异： 作者发现，他们完美的“真空”模拟无法完美匹配真实的“地铁站”数据。真实的峰值更宽且发生了偏移。他们怀疑金属传感器扭曲了电子波，这种现象被称为“基质效应（matrix effects）”。

主要发现：一次更精准的测量

尽管模拟结果并未与混乱的现实世界数据完全吻合，但它足以修正一个此前略有偏差的测量值。

• 旧数值： 科学家此前认为，每当原子捕获一个内层“K电子”时，它捕获外层“L电子”的次数大约是前者的7次（比例为 0.070）。
• 新数值： 利用这种更准确的摇晃模型，他们重新计算了这个比例。他们发现旧模型低估了“L壳层”的捕获量。新的、更准确的比例是 0.0756。

这为什么重要

这听起来可能只是一个微小的数字，但在寻找幽灵粒子的世界里，这意义重大。

1. 更清晰的信号： 通过了解“家具”究竟如何摇晃，科学家可以更准确地减去背景噪声。这使得“幽灵粒子”的信号更加清晰可见。
2. 避免假警报： 本文证实，电子复杂的摇晃过程并不会在科学家关注的能量范围（60–108 eV）内产生看起来像惰性中微子的假信号。这让科学家确信，如果他们在那里看到了信号，那一定是真实的。
3. 面向未来： 作者承认他们的模拟是针对孤立原子的。下一步是研究如何模拟金属传感器内部的原子，从而更接近现实情况。

总结： 本文构建了一个超级精确的计算机模型，用于模拟原子在衰变时的“摇晃”过程。虽然模型显示现实世界的传感器材料会让情况变得复杂，但新的数学方法让科学家纠正了一个长期的测量误差，为他们搜寻宇宙中缺失的幽灵粒子提供了更锋利的工具。

技术摘要

技术摘要：在原子 $^{7}\text{Be}$ 的电子俘获衰变中研究振激（Shake-up）与振离（Shake-off）能谱

问题陈述
BeEST（超导隧道结中的铍电子俘获）实验旨在通过测量 $^{7}\text{Li}$ 的反冲能量谱，为亚 MeV 惰性中微子的存在设定严格的限制。目前的灵敏度受限于对电子振激（SU）和振离（SO）能谱建模的系统不确定性。这些效应源于电子俘获过程中核势能的突然变化，导致剩余的原子电子发生激发（SU）或电离（SO）。以往的模型依赖于非相对论近似，且缺乏最先进的电子相关效应，导致在描述实验数据中的谱线尾部和峰值展宽时存在误差。此外，L/K 电子俘获比是一个关键参数，对于惰性中微子的搜索以及天体物理学中的 $^{7}\text{Li}$ 产生模型均至关重要，由于上述建模局限性，该参数之前的测定存在显著的系统不确定性。

方法论
作者利用 ab initio MCDGME 代码中的多构型 Dirac-Fock（MCDF）形式体系，计算了 $^{7}\text{Be}$ 的原子结构以及生成的 $^{7}\text{Li}$ 子原子及其离子。关键方法论特征包括：

• 相对论与相关波函数： 计算采用了无对哈密顿量（no-pair Hamiltonian），包括库仑、高特（Gaunt）相互作用和延迟相互作用，以及量子电动力学（QED）辐射修正（自能和真空极化）。
• 构型相互作用： 为了考虑电子相关效应，通过单激发和双激发将基组空间扩展至 4s 轨道。
• 骤变近似（SA）： SU 和 SO 概率是基于初始 $^{7}\text{Be}$ 波函数与最终 $^{7}\text{Li}$（及电离态）波函数之间的重叠度计算得出的。该形式体系明确计算了 K 壳层和 L 壳层俘获的单次和双次振激过程（包括 SU 和 SO）。
• 解析近似： 为了便于拟合实验数据，通过使用解析函数（K-SO 使用幂律，L-SO 使用对数正态分布）对数值计算得到的 SO 能谱进行近似，使解析结果与 ab initio 结果的约化卡方值保持较低水平。

主要结果

1. 振激概率： 计算表明，与 K 壳层俘获相比，L 壳层俘获后的振激概率显著更高。具体而言，振激效应约占所有 L 俘获事件的 30%。
2. 双次振激过程： 研究计算了双次 SU 和双次 SO 事件的概率。虽然这些过程目前低于 BeEST 实验的探测阈值（对 L 俘获贡献 <2.5%，对 K 俘获贡献 <0.2%），但它们是不可忽视的，并在 136 eV 附近产生约 24% 的相对谱差异，在 225 eV 附近产生约 12% 的差异。
3. 能谱特征：
   • 振激（Shake-up）： L 俘获的 K-SU 跃迁（$1s^2 2s \to 1s^2 2s^2$）产生一个与 K-GS 峰能量相同的峰，解释了为何需要在 K-GS 拟合中引入第二个高能组分。
   • 振离（Shake-off）： K-SO 能谱延伸至比以往 Levinger 模型预测的更高能量。研究发现，L 俘获后的 K-SO 能谱是不可忽视的，在 160 至 190 eV 之间对能谱有贡献。
4. 与实验对比： 与 BeEST 第三阶段（Phase-III）数据对比时，ab initio 模拟可以解释部分但并非全部能谱特征。峰值质心和宽度的差异归因于固态效应（与钽（Ta）传感器的矩阵相互作用）以及潜在的探测器非线性。具体而言，Ta 晶格环境可能会修改波函数，从而影响靠近电离阈值的跃迁概率。
5. 修正后的 L/K 俘获比： 通过引入新的、更准确的 SU 和 SO 模型，并移除来自 Si 基底伽马相互作用的背景（通过第三阶段的符合甄别），作者修正了钽中 $^{7}\text{Be}$ 的 L/K 电子俘获比。该值从 0.070(7) 更新为 0.0756(20)。

意义与主张
论文声称，新的 ab initio 模型为分析 BeEST 能谱提供了更严谨的理论基础，降低了寻找惰性中微子的系统不确定性。修正后的 L/K 比是一个重要的成果，因为它影响了与宇宙学 $^{7}\text{Li}$ 产生相关的恒星环境中 $^{7}\text{Be}$ 衰变的计算。作者指出，尽管其计算仅限于孤立原子，但他们成功识别了特定的能谱组分（例如 K-GS 峰的高能肩部），而这些组分此前无法得到解释。他们强调，未来的工作必须纳入矩阵效应（例如通过密度泛函理论），以完全解决计算结果与测量概率之间的差异，特别是在靠近电离阈值处。研究结论认为，没有任何预测的能谱特征会在感兴趣区域（60–108 eV）内模拟出惰性中微子的信号，这增强了当前搜索限制下背景建模的有效性。