[论文解读] High-Dimensional Bayesian Likelihood Normalisation for CRESST's Background Model

G. Angloher, S. Banik+60 more|arXiv (Cornell University)|Jul 19, 2023

Dark Matter and Cosmic PhenomenaPhysics and Astronomy被引用 3

一句话总结

本文提出了一种用于CRESST电磁背景模型的高维贝叶斯似然归一化方法，无需假设长期平衡即可改善能谱拟合。通过贝叶斯推断将Geant4模拟的能谱与实验数据进行拟合，该方法在[1 keV, 40 keV]范围内解释了高达82.7%的背景，相比之前的方法提高了18.6%，并通过完整协方差处理增强了鲁棒性，减少了对先验假设的依赖。

ABSTRACT

Using CaWO$_4$ crystals as cryogenic calorimeters, the CRESST experiment searches for nuclear recoils caused by the scattering of potential Dark Matter particles. A reliable identification of a potential signal crucially depends on an accurate background model. In this work we introduce an improved normalisation method for CRESST's model of the electromagnetic backgrounds. Spectral templates, based on Geant4 simulations, are normalised via a Bayesian likelihood fit to experimental background data. Contrary to our previous work, no assumption of partial secular equilibrium is required, which results in a more robust and versatile applicability. Furthermore, considering the correlation between all background components allows us to explain 82.7% of the experimental background within [1 keV, 40 keV], an improvement of 18.6% compared to our previous method.

研究动机与目标

• 开发一种在低能暗物质搜索背景下，更稳健且更精确的CRESST电磁背景模型归一化方法。
• 消除背景归一化中对显式长期平衡假设的依赖，从而提升模型的灵活性与可靠性。
• 改善模拟与实验背景数据在多个同位素组分中的能谱一致性。
• 通过正确处理背景组分之间的相关性，实现更全面且可扩展的背景模型。

提出的方法

• 采用贝叶斯似然框架，将Geant4模拟的能谱模板归一化至实验背景数据。
• 该方法使用一个似然模型，通过完整协方差矩阵考虑所有背景组分之间的相关性。
• 为污染水平和干扰参数分配先验，并对干扰参数进行边缘化处理，以提高鲁棒性。
• 通过最大后验估计（MAP）实现归一化，从而支持系统不确定性传播。
• 通过似然函数让数据约束相对活度，避免显式依赖长期平衡假设。
• 使用马尔可夫链蒙特卡洛（MCMC）采样实现方法，以探索后验分布并评估不确定性。

实验结果

研究问题

• RQ1 是否可以通过贝叶斯似然归一化方法，在不假设长期平衡的前提下，改善模拟与实验背景数据之间的能谱一致性？
• RQ2 在背景组分之间引入完整协方差后，对归一化方法的准确性和鲁棒性有何影响？
• RQ3 与以往方法相比，新方法在[1 keV, 40 keV]能量范围内能多大程度上解释观测到的实验背景？
• RQ4 不同背景组分的推导污染活度如何比较？这些结果对残留放射性提供了哪些新见解？

主要发现

• 新提出的贝叶斯似然归一化方法在[1 keV, 40 keV]能量范围内，最多可解释82.7%的实验背景，相比之前使用的高斯拟合方法，最大提升达18.6%。
• 该方法减少了对长期平衡假设的依赖，使背景模型更具鲁棒性和通用性，适用于多种探测器模块。
• 对于内部放射性背景，该方法得出了具有小不确定度的一致活度估计，例如238U为1110 ± 10 µBq/kg，并显著改善了对低活度组分（如214Bi和210Po）的约束。
• 该方法成功处理了相似污染物之间的相关性（如238U和235U衰变链中的污染物），从而得到更一致且更符合物理实际的活度估计。
• 对于近外部和额外外部放射性背景，该方法提供了更紧致的约束，许多组分被抑制至亚µBq/kg水平，表明高背景组分得到有效抑制。
• 该方法在能谱一致性方面表现更优，特别是在1–10 keV范围，该区域电磁与核 recoil 辨识性能下降，因此显著提升了对潜在暗物质信号的探测灵敏度。