Search for Annual Modulation in Combined ANAIS 112 and COSINE 100 Data using Bayesian Model Comparison

通过对约三年 ANAIS-112 与 COSINE-100 联合数据的贝叶斯模型比较，该研究未发现暗物质诱导的年调制信号的证据，因为贝叶斯因子强烈支持恒定背景模型而非正弦模型。

要点

以下是用通俗语言和日常类比对该论文的解读。

大谜团：什么是暗物质？

想象我们的银河系是一片巨大而隐形的海洋。我们知道这片海洋存在，是因为在其中航行的“船只”（恒星）表现出某种行为，暗示它们正被某种看不见的东西推动。科学家将这种看不见的物质称为暗物质。

几十年来，主流理论认为这种暗物质是由微小的、幽灵般的粒子构成的，这些粒子被称为WIMPs（弱相互作用大质量粒子）。这些粒子非常“害羞”，极少与普通物质发生碰撞，但如果真的发生了碰撞，它们可能会留下一道微小的“划痕”或信号。

争议："年度调制”信号

1998 年，一个名为DAMA/LIBRA的实验声称发现了这些“划痕”。他们表示观察到了一种特定的模式：“划痕”的数量每年像心跳一样上下波动一次。

类比：想象你站在海滩上。你注意到每年六月左右，海浪会稍微变大，而到了十二月左右，海浪会变小。于是你得出结论：“海洋在呼吸！”这就是 DAMA 所声称的：地球穿行于暗物质的“风”中，随着我们绕太阳公转，我们在六月撞上更多的暗物质，而在十二月撞上更少。

怀疑论者：COSINE-100 和 ANAIS-112

其他科学家对此持怀疑态度。他们建造了自己的探测器，即位于韩国的COSINE-100和位于西班牙的ANAIS-112，使用了与 DAMA 完全相同的材料，看看是否能听到同样的“海洋呼吸”声。

很长一段时间里，这些新实验表示：“我们什么也没听到。”他们查看数据后，没有发现明显的上下波动模式。

新研究：用新工具重新审视

在这篇论文中，两位研究人员（Om Godse 和 Shantanu Desai）决定重新审视 COSINE-100 和 ANAIS-112 的合并数据（约三年的数据）。

他们没有仅仅查看数字，而是使用了一种名为贝叶斯模型比较的统计工具。

类比：想象你是一名试图破案的侦探。你有两个理论：

1. 理论 A（余弦模型）：海浪是有节奏的。存在真实的模式（即“海洋在呼吸”）。
2. 理论 B（恒定模型）：海浪只是随机噪声。没有模式；那只是一条平坦、枯燥的线。

通常，你可能只是查看数据并猜测哪个理论更吻合。但这篇论文使用了一个“数学天平”（称为贝叶斯因子）来权衡证据。它问的是：“数据来自有节奏的模式而非随机噪声的可能性，究竟大了多少？”

结果：天平没有倾斜

研究人员使用不同的“规则”（称为先验）来测试这个天平，看看答案是否会因数学设置方式的不同而改变。他们在两个不同的能量范围内查看了数据（就像用两种不同颜色的滤镜观察海洋）。

裁决：
结果很明确：天平几乎没有移动。

• “数学天平”给出的分数低于1.15。
• 在统计学界，如此低的分数就像在嘈杂的房间里听到微弱的耳语。它几乎不值一提。

这意味着：COSINE-100 和 ANAIS-112 的数据不支持存在有节奏的“海洋呼吸”（年度调制）这一观点。数据看起来更像是随机噪声（恒定值），而非来自暗物质的信号。

结论

该论文得出结论：当你将这两个实验的数据合并并通过这种严格的数学测试时，没有证据支持 DAMA 声称看到的年度调制信号。

研究人员还向公众分享了他们的“配方”（计算机代码），以便任何人都可以检查他们的工作。这证实了早期的发现：“海洋呼吸”可能只是光线造成的错觉，而非真实现象。

技术摘要

技术摘要：ANAIS-112 与 COSINE-100 联合数据的贝叶斯模型比较

问题陈述
二十多年来，DAMA/LIBRA 实验一直在钠碘（NaI）探测器的剩余事件率中报告具有统计显著性的年调制信号，这一特征与弱相互作用大质量粒子（WIMP）暗物质相互作用一致。然而，这一主张尚未被其他使用相同靶材料的直接探测实验所证实，特别是韩国的 COSINE-100 和西班牙的 ANAIS-112。虽然此前对 COSINE-100 和 ANAIS-112 联合数据的前三年（称为 C25）进行的联合分析利用频率学派回归未发现调制证据，但本研究旨在通过贝叶斯框架补充这些发现。所解决的核心问题是：这两个实验的联合剩余数据是否提供了统计证据，支持优于常数背景模型的正弦调制模型（预期来自暗物质），且独立于先前工作中使用的具体回归技术。

方法论
作者对两个假设进行了贝叶斯模型比较：

• 模型 $M_2$（余弦信号）： 剩余数据遵循正弦函数 $R(t) = S_m \cos(\omega(t - t_0))$，其中 $S_m$ 为振幅，$\omega$ 为角频率，$t_0$ 为相位。
• 模型 $M_1$（常数背景）： 剩余数据由常数 $A$ 描述。

分析利用了 COSINE-100 和 ANAIS-112 约三年的联合数据，总曝光量为 485 千克·年。数据被划分为 15 天的时间间隔，涵盖两个能量范围：1–6 keV 和 2–6 keV。

关键的方法步骤包括：

1. 似然构建： 假设观测到的剩余值服从高斯似然，将误差视为高斯分布。
2. 先验选择： 为确保稳健性，作者针对自由参数（$S_m$、$\omega$、$t_0$）测试了四组不同的先验：
   • 振幅（$S_m$）： 均匀先验，范围从最小到最大观测剩余值，以及单独从 0 到最大绝对剩余值。
   • 频率（$\omega$）和相位（$t_0$）： 均匀先验（分别覆盖整个数据集持续时间和 0–365 天），以及以 DAMA/LIBRA 报告的最佳拟合值为中心的正态先验（$T \approx 1$ 年，$t_0 \approx 145$ 天）。
3. 证据计算： 使用 Dynesty 中实现的嵌套采样算法计算两个模型的贝叶斯证据（边缘似然）。比较指标为贝叶斯因子的自然对数 $\ln(B_{21})$，其中 $B_{21}$ 是余弦模型证据与常数模型证据的比率。

主要贡献

• 独立验证： 这项工作使用贝叶斯模型比较方法，而非原始联合分析中使用的最小二乘法和马尔可夫链蒙特卡洛（MCMC）回归，对 C25 结果提供了独立验证。
• 先验稳健性： 通过测试多种先验配置（包括那些基于 DAMA/LIBRA 最佳拟合值 informed 的配置），该研究证明了关于不存在信号的结论对先验假设的变化具有稳健性。
• 开放科学： 作者已公开其分析代码，以确保贝叶斯证据计算的透明度和可重复性。

结果
分析得出了所有四个先验组在两个能量区间内的贝叶斯因子自然对数（$\ln(B_{21})$）。结果总结如下：

• 对于所有测试的先验组合，$\ln(B_{21})$ 的值均小于 1.15。
• 根据解释贝叶斯因子的 Jeffreys 标度，$\ln(B_{21}) < 1.15$（对应 $B_{21} \leq \sqrt{10}$）的值表明，余弦模型相对于常数模型的证据“几乎不值一提”。
• 具体而言，在 1–6 keV 和 2–6 keV 区间内，计算值范围约为 -0.13 至 1.36，大多数值聚集在 1.15 阈值附近或低于该阈值。
• 振幅、频率和相位的边缘化后验轮廓显示，对于 DAMA 偏好的参数，没有表现出比常数速率的零假设更强的偏好。

意义与主张
本文结论认为，贝叶斯模型比较没有提供证据支持联合 ANAIS-112 和 COSINE-100 数据包含与暗物质相互作用一致的年调制信号的假设。结果与常数背景假设一致。

作者谦逊地将他们的贡献表述为对 C25 分析的补充。他们指出，他们的发现加强了以下结论：这些实验前三年的联合数据不支持 DAMA 的年调制主张。该研究并未声称能彻底排除一般的暗物质相互作用，但具体发现，无论关于信号相位和频率的先验假设如何，在该特定数据集中，正弦信号相对于常数背景均无统计偏好。