Reassessing the gallium anomaly using self-consistent electron wave functions

本文通过使用由狄拉克-库仑方程导出的自洽电子波函数来计算中微子俘获截面，重新评估了长期存在的镓异常，从而更新了该差异的全局显著性及其关于惰性中微子的解释。

要点

想象一下，你正试图计算有多少滴雨水落入了一个特定的水桶中。你有一个非常精确的数学公式，可以根据降雨量的大小和水桶的大小，精准地预测应该有多少滴雨水落下。然而，每次你在现实生活中实际计数时，都会发现落下的雨滴比你的公式预测的要少。这种“丢失的雨水”就是物理学家所说的“镓异常”（Gallium Anomaly）。

三十多年来，使用镓-71（一种金属类型）作为“水桶”来捕捉中微子（来自太阳或放射源的微小、幽灵般的粒子）的实验，一致发现捕捉到的中微子比预期的要少。这种预测与现实之间的差距已经不断扩大，以至于它现在被视为一个重大的物理学谜题。

Cadeddu及其同事的这篇论文，就像是一群顶尖机械师决定从头开始重建那个预测公式的引擎，看看是否在计算过程中出了错。

旧方法 vs. 新方法

旧引擎（近似法）：
此前，科学家们使用一个“草稿版”的数学版本来计算中微子如何与镓原子发生相互作用。他们将电子（围绕原子核运动的微小粒子）视为简单的、平滑的波，且在原子内部几乎不发生变化。这就像是通过观察一张平坦的地图来估算一条崎岖道路的形状。他们假设电子波在微小的原子核内部任何地方都是相同的。

新引擎（精确解）：
在这项新研究中，作者决定不再使用那张平坦的地图。相反，他们使用了高精度的GPS，通过求解描述电子实际行为的精确方程（称为狄拉克-库仑方程，Dirac-Coulomb equation）。

• 类比： 想象原子核是一个拥挤的舞池。旧方法假设舞池里的每个人都静止地站在一个完美的圆圈里。而新方法则实际统计了每一位舞者，考虑了他们在拥挤空间中如何互相碰撞以及如何移动。他们利用专门的计算机程序，求解了既包括束缚在原子上的电子、也包括飞出的电子的精确“波形”。

“平均化”技巧

这篇论文中的另一个关键变化是他们处理原子核大小的方式。

• 旧方法： 他们选择了原子核中心的一个单一点（就像通过将温度计插入房间正中间来测量房间温度），并假设这个点代表了整个原子核。
• 新方法： 他们意识到原子核是有体积的，因此他们在整个原子核的体积内对电子的行为进行了“平均”。这就像是在测量房间内每个位置的温度并找到真正的平均值，而不是仅仅根据中心点进行猜测。

他们发现了什么？

当他们运行新的、更精确的计算时：

1. 预测发生了变化： 他们新的、更准确的公式预测捕捉到的中微子会比旧公式预测的要少。
2. 差距扩大了： 因为新的预测值更低，所以“预期看到的内容”与“实验实际看到的内容”之间的差异变得更大了。
3. 结果： “缺失的中微子”现在成为了一个 5.5-sigma 的问题。在科学界，“sigma”是衡量置信度的一种度量。5-sigma 的结果是发现性的金标准，意味着这种差异仅仅是随机巧合的可能性小于 350 万分之一。

“惰性中微子”假说

物理学家们有一个解释这些丢失雨水的热门理论：惰性中微子（Sterile Neutrinos）。

• 隐喻： 想象中微子是向水桶飞行的鸟群。该理论认为，其中一些鸟是“隐形的”（惰性的）。它们完全不会与水桶发生相互作用；它们只是直接从旁边飞过。如果这些隐形的鸟确实存在，就能解释为什么水桶比预期的要空。

作者更新了数学模型，以观察这个“隐形鸟”理论是否仍然成立。

• 好消息： 数学仍然允许这些隐形的鸟存在。镓实验的数据仍然强烈指向它们的存在。
• 坏消息： 其他实验（观察反应堆中微子、太阳中微子以及质量测量）已经设置了“围栏”，表明这些隐形的鸟不能以镓数据所暗示的那种方式飞行。镓数据希望这些鸟是非常活跃的，但其他的围栏却说它们必须非常害羞。

结论

作者并没有解决这个谜团；他们实际上让这个谜团变得更加神秘了。通过使用更好的数学模型和更精确的电子模型，他们证实了“缺失的中微子”是一个真实且稳健的问题，而不是一个计算错误。

他们得出结论，虽然“惰性中微子”仍然是头号嫌疑对象，但它目前正处于与其他实验证据的僵持状态。谜团依然悬而未决，作者建议可能需要一种使用不同类型探测器的实验，才能最终抓获这个“罪魁祸首”。

技术摘要

技术摘要：利用自洽电子波函数重新评估镓异常

问题陈述
本文探讨了“镓异常”（gallium anomaly）问题，这是一个在利用 $^{71}\text{Ga}$ 作为靶材的中微子俘获实验中观察到的持续差异。使用 $^{51}\text{Cr}$ 和 $^{37}\text{Ar}$ 放射性源（GALLEX、SAGE 以及最近的 BEST）的实验一致测量到，中微子诱导的反 $\beta$ 衰变（IBD）速率系统性地低于理论预测值。这一亏损目前的显著性已超过 $4\sigma$，挑战了理论截面计算的稳健性，并引发了关于超越标准模型物理（如惰性中微子）的假设。作者指出，以往的理论评估依赖于电子波函数的领先阶（LO）近似，且对核矩阵元与轻子矩阵元的处理存在不一致性。

方法论
作者通过放弃传统的 LO 近似，转而采用自洽数值方法，对 $^{71}\text{Ga}$ 的理论中微子俘获截面进行了全面的修正。关键方法组成部分包括：

1. 精确波函数解： 作者不再使用近似波函数，而是通过数值求解狄拉克-库仑方程（Dirac-Coulomb equation），以获得结合态（电子俘获，EC）和连续态（IBD）的电子态。这是利用基于 RADIAL 库的定制软件程序包，通过狄拉克-哈特里-福克-斯莱特（DHFS）方程实现的。
2. 自洽势场： 计算采用了包含核势、电子势和交换势的 DHFS 势场。核势采用两参数费米分布进行建模，并纳入了对 $^{71}\text{Ge}$ 核电荷半径（$R_{ch} = 4.032 \pm 0.002$ fm）的首次直接测量结果。交换势参数（$C_{ex}$）取平均值 $1.25$，以平衡 Slater（$3/2$）与 Kohn-Sham（$1$）近似之间的模型差异。
3. 平均化程序： 作者超越了在单一点（例如 $r_0 = 0$ 或 $r_0 = R_{box}$）评估波函数的惯例，提出对电子波函数和费米函数在核电荷密度（$\rho_{ch}(r)$）上进行平均。该方法同时应用于 IBD 截面计算以及通过细致平衡原理用于提取核矩阵元的电子俘获（EC）速率计算。
4. 更新后的输入参数： 分析整合了最新的实验输入，包括：
   • $^{71}\text{Ge}$ EC 半衰期的四个独立测量值的加权平均值（$t_{1/2} = 11.465 \pm 0.003$ d）。
   • 更新后的 EC 过程 Q 值（$232.47 \pm 0.09$ keV）。
   • 精炼后的电子俘获概率比（$P_L/P_K$ 和 $P_M/P_K$）及结合能。
   • 源自 $(p,n)$ 和 $(^3\text{He}, ^3\text{H})$ 散射数据的激发态伽莫夫-泰勒（Gamow-Teller）强度。

主要贡献与结果

• 修正后的截面： 由于采用了平均化程序，新的计算得出的基态截面比之前的估计值低约 2.8%。在计入激发态贡献后，$^{51}\text{Cr}$ 和 $^{37}\text{Ar}$ 源的总截面比最近文献（文献 [27]）中报道的数值小约 2%。
• 异常的显著性： 通过将更新后的理论预测与 GALLEX、SAGE 及 BEST 的实验数据进行对比，作者重新评估了该异常的显著性：
  • 使用 $(p,n)$ 散射数据进行激发态贡献计算时，差异达到 $5.5\sigma$（$R = 0.835^{+0.030}_{-0.048}$）。
  • 使用 $(^3\text{He}, ^3\text{H})$ 散射数据时，差异为 $4.7\sigma$（$R = 0.811^{+0.037}_{-0.035}$）。
  • 使用仅基态截面的保守估计，其异常显著性为 $3.8\sigma$。
• 自洽性： 作者强调，通过在用于 IBD 的费米函数和用于 EC 的核处电子密度计算中使用相同的软件和底层假设，系统性偏差在比例中被抵消，从而使结果更加稳健。

意义与解释
论文声称，镓异常的显著性现在比之前报告的更高，达到了 $5.5\sigma$。作者重新审视了将此异常解释为短基线主动-惰性中微子振荡（3+1 模型）的可能性。其结果证实，镓数据（包括精确的 BEST 测量）需要较大的主动-惰性中微子混合角（在 $2\sigma$ 置信度下，对于 $\Delta m^2_{41} \gtrsim 0.98 \text{ eV}^2$，有 $0.16 \lesssim \sin^2 2\vartheta_{ee} \lesssim 0.54$）。

然而，作者指出，这种所需的较大混合角与反应堆反中微子数据、太阳中微子限制以及 KATRIN 中微子质量测量结果存在张力。因此，本文结论认为，尽管该异常的统计显著性有所增加，但惰性中微子的解释仍与其它实验数据存在冲突。作者表示，解决镓异常的问题仍然是一个开放性问题，可能需要使用不同探测器的全新源实验来解决这一分歧。