⚛️ nuclear theory

Dark Matter Inelastic Scattering with Nuclei for Direct Detection

本文通过使用相对论构型相互作用密度泛函理论，研究了暗物质直接探测中 WIMP-原子核散射的核响应，特别关注了 $^{129}$ Xe 和 $^{131}$ Xe 同位素中的非弹性通道，旨在证明非弹性贡献可能显著超过弹性贡献，并对探测结果的解释产生实质性影响。

原作者： Shao-Feng Ge, Oleg Titov, Yakun Wang

发布于 2026-09-16

📖 1 分钟阅读🧠 深度阅读

CC BY 4.0

原作者： Shao-Feng Ge, Oleg Titov, Yakun Wang

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

几十年来，物理学家一直在寻找构成我们宇宙大部分物质的不可见物质。他们称之为暗物质，尽管他们无法直接看到它，但他们相信它无处不在，在恒星和行星之间的空间中漂浮。主流理论认为，这种暗物质是由沉重、运动缓慢的粒子组成的，它们极少与普通物质发生相互作用。为了捕捉它们的踪迹，科学家们在深层地下建造了巨大的探测器，里面充满了数吨液氙。这个想法很简单：如果一个暗物质粒子撞击一个氙原子，它应该会给原子一个微小的踢力，产生一道闪光，从而被灵敏的传感器记录下来。多年来，研究人员一直关注最常见的碰撞类型，即暗物质粒子撞击原子并弹开，使原子保持在正常的、平静的状态。这被称为弹性散射。

然而，还有另一种可能性，直到现在才被很大程度上忽视了。就像台球可以撞击另一颗球并使其旋转或振动一样，暗物质粒子可能会以足够的能量撞击氙原子，使其进入更高的激发态。这被称为非弹性散射。在这种情况下，原子不仅会发生反冲，还会保留一部分能量，等待稍后释放，以回到稳定状态。由于激发原子所需的能量是特定的，因此这一过程只有在暗物质粒子运动足够快且足够重的情况下才会发生。长期以来，科学家们一直认为这些被“摇晃”的原子非常罕见，或者难以计算。但一项新的研究表明，忽略这第二种类型的碰撞可能是一个严重的错误，可能会掩盖研究人员正在寻找的正信号。

来自上海交通大学和北航的研究团队对暗物质如何与氙原子相互作用进行了全新的审视，特别关注了这些激发态。他们意识到，要理解碰撞过程中原子内部发生了什么，你需要一张非常详细的原子内部结构图。之前的研究使用了一种称为核壳层模型的计算方法，该方法将原子视为一个具有固定房间数量（粒子居住其中）的建筑。虽然这种方法很有用，但对于像氙这样的重原子来说，使用起来变得非常困难，因为可能的排列组合数量庞大到如果不进行简化处理就无法处理。新团队使用了一种更现代、更强大的方法，称为相对论配置相互作用密度泛函理论。这种方法将原子内的每一个粒子都视为活跃的参与者，从而能够更完整、更准确地描绘出原子在受到撞击时的行为。

他们的计算表明，非弹性通道并非微不足道的注脚；对于某些类型的暗物质相互作用，它的重要性可以与弹性通道并驾齐驱，在某些特定情况下，甚至可以强上千倍。研究人员发现，当暗物质撞击两种具有核自旋的特定类型的氙原子时，原子可以被激发到更高的能量水平。在某些情景下，这些激发原子的信号非常强烈，以至于完全主导了总信号，使标准的弹性碰撞显得微不足道。此外，他们还发现，不仅仅是第一激发态很重要。如果暗物质粒子足够重，它可以将原子踢入更高的能量级，而这些高能级状态会对探测到的总事件数做出显著贡献。事实上，对于某些相互作用，第二激发态的信号强度几乎与第一激发态的信号强度相当。

这些发现的影响是深远的，例如对于像 LUX-ZEPLIN (LZ) 探测器这样最近报告了一个看起来像暗物质碰撞事件的实验。该事件的能量相当高，约为 248 千电子伏特。之前的分析假设这是一个弹性碰撞，并试图使用旧的计算方法来解释它。然而，当新团队将他们更准确的模型应用于这个特定事件时，情况发生了戏剧性的变化。他们的计算显示，通过标准弹性碰撞观察到如此高能事件的概率实际上比之前认为的要低得多。相反，新模型表明，在该能量下的事件率主要由非弹性通道主导，但即便如此，预测的事件数量也远低于旧模型所暗示的数量。这意味着，LZ 合作组观察到的这个单次事件可能根本不是暗物质的迹象，或者至少不是他们正在测试的那种暗物质。

这项研究还强调了如何解读探测器数据的一个更深层次的问题。信号的形状——即事件数量随能量反冲变化的规律——取决于使用哪种计算方法，其表现形式截然不同。旧方法预测事件数量会在 LZ 探测器观察到信号的能量附近出现峰值。而新的、更详细的方法则预测在同一区域会出现凹陷，即信号实际上处于最低点。这种差异并非微小的误差，而是二十五倍的差距。如果新的计算是正确的，那么 LZ 看到的信号就很难被解释为暗物质探测，因为来自暗物质的预期背景远低于之前的假设。

最终，这项工作作为一个至关重要的提醒：用于解释实验数据的工具与数据本身同样重要。通过使用一种更复杂的建模方式来模拟原子核，研究人员展示了支配暗物质如何与物质相互作用的规则比此前认为的更加复杂和多样。他们证明了，忽略原子被摇晃进入激发态的可能性，可能会导致科学家错过故事的主要部分，或者相反，将一个随机的波动误认为是发现。随着探测器变得越来越灵敏并能够观察到这些细微细节，对同样精确的理论模型的需求变得至关重要。寻找暗物质的过程不仅仅是建造更大的液氙罐，也是关于如何以唯有最先进理论才能提供的清晰度，去理解那些原子内部复杂的粒子之舞。

技术摘要：暗物质与原子核的非弹性散射及其在直接探测中的应用

问题陈述
传统的弱相互作用大质量粒子（WIMP）直接探测主要关注弹性散射，即靶核保持在其基态。然而，对于足够重的 WIMP，动能转移可以将靶核激发到低能激发态，从而导致非弹性散射。尽管先前的研究已经承认了非弹性通道的潜力，特别是针对氙同位素中的自旋相关（SD）相互作用，但现有的理论框架面临两个主要局限：

模型空间限制： 大多数计算依赖于核壳模型（NSM），由于组合爆炸，这要求对像氙（ $^{129,131}$ Xe）这类重核进行配置空间的截断。这些截断引发了对核响应函数收敛性和准确性的质疑。
能级覆盖有限： 先前的 NSM 研究主要集中在向第一激发态的跃迁。尽管对于重暗物质（DM）而言，高阶激发态的贡献在运动学上是可达的，但这些状态对总事件率的贡献尚未得到系统的探索。

方法论
作者通过采用一种结合了非相对论有效场论（NR EFT）与相对论构型相互作用密度泛函（ReCD）理论的统一理论框架来解决这些局限性。

运动学与形式体系： 本研究建立了非弹性散射的运动学约束，推导了激发天然氙中特定核能级所需的最小暗物质质量。研究将 NR EFT 形式体系推广到包含非弹性跃迁，利用有效算符的多极展开。这种方法允许根据角动量、宇称和时间反演对称性导出严格的选择定则，从而确定哪些算符对弹性与非弹性通道有贡献。
核结构计算（ReCD）： 作者没有使用 NSM，而是采用了 ReCD 理论。该方法将所有核子视为大型单粒子空间中的活跃粒子，避免了 NSM 中常见的严重模型空间截断。它通过构型混合（不同准粒子构型的混合）和角动量投影（以恢复由平均场近似打破的旋转对称性）来引入多体相关性。
对比： 作者明确对比了其 ReCD 结果与现有的 NSM 计算（特别是用于参考文献 [20] 和 [22] 的计算），以量化核结构建模对预测事件率的影响。

主要贡献与结果

非弹性与弹性速率的重新评估：
使用 ReCD，作者计算了所有 15 个 NR EFT 算符的非弹性与弹性事件率之比（ $R_{inel}/R_{el}$ ）。
- 对于算符 $O_1, O_3, O_8, O_{11}, O_{12}, O_{15}$ ，非弹性贡献可以忽略不计（ $\lesssim 1\%$ ）。
- 对于算符 $O_4, O_5, O_6, O_9, O_{10}$ ，非弹性贡献显著（ $\sim 10\%$ ）。
- 对于算符 $O_7, O_{13}, O_{14}$ ，非弹性通道可以主导或与弹性通道相当。具体而言，对于 $O_{13}$ ，非弹性速率与弹性速率相当。对于 $O_7$ 和 $O_{14}$ ，非弹性贡献可以超过弹性速率 $10^2$ 到 $10^3$ 倍。
高阶激发态的影响：
本文定量地证明了（除了第一激发态之外的）高阶激发态并不总是次要的。
- 对于 $m_\chi = 10$ TeV， $^{131}$ Xe 的第二激发态对某些算符有显著贡献。
- 对于算符 $O_9$ ，第二激发态的贡献约为弹性速率的 $70\%$ ，超过了第一激发态（ $\sim 50\%$ ）的贡献。包括这些高阶态可以使 $^{131}$ Xe 的总预测事件率增加 2 倍。
与核壳模型（NSM）的差异：
ReCD 计算揭示了与 NSM 结果的显著偏差，特别是在高动量转移（ $|q| > 100$ MeV）情况下。
- 响应函数： 通过 ReCD 计算的核响应函数（ $W_T$ ）与 NSM 结果存在实质差异。例如，对于算符 $\Sigma''$ ，ReCD 的结果大约是 NSM 结果的五倍。
- 振荡行为： NSM 结果在特定的动量转移处表现出剧烈的凹陷，而 ReCD 结果则显示出更平滑的振荡。这些差异源于底层核密度矩阵的变化，特别是来自特定单粒子轨道（如 $1d_{3/2}$ 与 $2s_{1/2}$ ）贡献的相对振幅。
对 LZ 事件解释的影响：
作者将他们的发现应用于近期 LUX-ZEPLIN (LZ) 观测到的一个具有约 $E_R \approx 248$ keV 互作用能的单事件。
- NSM 预测： 基于 NSM 的分析预测在 200–250 keV 附近存在微分事件率的局部极大值，这可以容纳观测到的事件。
- ReCD 预测： ReCD 计算预测在 $E_R \approx 250$ keV 处存在一个局部极小值，在该能量下，事件率大约是 NSM 预测值的 25 倍之小。此外，ReCD 预测在低反冲能量（ $\sim 50$ keV）处存在全局极大值，这与未观测到低能事件的事实相矛盾。
- 结论： 本文指出，将 LZ 事件解释为暗物质信号高度依赖于所使用的核模型；在 ReCD 输入下，该信号假设变得显著不再理想。

意义
本文声称，对于精确预测直接探测（特别是针对重 WIMP 和自旋相关相互作用），需要一个可靠的核多体结构描述，这必须超越 NSM 受限的构型空间。通过利用 ReCD 理论，作者证明了：

非弹性散射是特定相互作用算符的总信号中的关键组成部分，其贡献可能主导事件率。
在分析中必须包含高阶激发态，以避免低估事件率。
核模型的选择（ReCD 与 NSM）会导致预测的事件率和谱形状出现数量级的差异，从而显著改变对实验数据（如近期 LZ 事件）的解释。

这项工作为解释当前及未来的基于氙的直接探测实验提供了更稳健的理论基础，强调了在寻找暗物质的过程中进行先进核结构计算的必要性。

类似论文