Combined KM3NeT/ARCA and ANTARES search for compact neutrino sources

技术摘要：KM3NeT/ARCA 与 ANTARES 的联合紧凑中微子源搜索

问题与目标
本研究旨在解决来自点状天体源的高能中微子搜索问题。主要目标是展示结合两个不同中微子探测器数据的潜力：ANTARES 望远镜（2008–2022 年运行于法国图龙附近）和 KM3NeT/ARCA 探测器（目前正在西西里岛附近建设中）。这两种仪器都为南天和银河系中心提供了最佳的可视性。本分析旨在通过建立一个基于分箱似然法（binned likelihood approach）的联合搜索框架，以提高相对于独立分析的灵敏度和发现潜力。

方法论
本分析利用了一个包含四个不同组成部分的组合数据集：

1. ANTARES 径迹（Track）数据集： 4541 天的活时间（livetime），分为不重叠的径迹和簇射（shower）数据集。
2. ANTARES 簇射数据集： 包含在上述 4541 天的分组中。
3. ARCA19 径迹： 来自 KM3NeT/ARCA 的径迹类事件，包含 19 条检测线（48.4 天）。
4. ARCA21 径迹： 来自 KM3NeT/ARCA 的径迹类事件，包含 21 条检测线（287.4 天）。

定义： “径迹”（Tracks）是指涉及缪子产生的事件，而“簇射”（Showers）是由所有中微子风味的中性流相互作用以及电子中微子的带电流相互作用所诱发的。

统计框架
本分析采用分箱似然法。将观测到的事件直方图 ($N$) 与基于参考通量 ($\Phi_0 = 10^{-8} \text{ GeV cm}^{-2} \text{ s}^{-1}$) 的可缩放信号估计值 ($S$) 以及背景估计值 ($B$) 进行比较。似然函数定义为：
$$\ln L = \sum_{i}^{N_{bins}} N_i \ln(B_i + \zeta S_i) - (B_i + \zeta S_i)$$
其中 $\zeta$ 是信号强度参数。测试统计量 $\lambda$ 是最佳拟合信号强度 ($\hat{\zeta}$) 与仅背景假设 ($\zeta=0$) 之间的对数似然比。

• 发现潜力（Discovery Potential）： 定义为实现 50% 伪实验中达到 $5\sigma$ 显著性所需的信号强度 ($\zeta_{5\sigma}$)。
• 灵敏度（Sensitivity）： 定义为使用 Neyman 方法计算出的 $\zeta$ 的中值 90% 置信水平（CL）上限。

背景建模与系统不确定性
为了解决与背景建模相关的系统不确定性，针对不同赤纬带进行了四种独立的分析案例：

• 案例 1： ARCA 背景完全由数据驱动；ANTARES 生成过程由数据驱动，评估过程由蒙特卡洛（MC）驱动。
• 案例 2： 与案例 1 类似，但 ANTARES 生成背景数据经过随机化处理（均匀方位分布）。
• 案例 3： 所有数据集在生成和评估阶段均使用 DEFT 算法进行背景估计。
• 案例 4： ARCA 由数据驱动；ANTARES 完全基于蒙特卡洛（MC）模型。

这些四个案例的结果分布定义了系统不确定性带。直方图根据角距离 ($\log_{10} \alpha$)、重建能量 ($\log_{10} E$) 以及对于 ANTARES 径迹而言的角误差估计 ($\log_{10} \beta$) 进行构建。

关键结果
本分析假设能量谱为 $E^{-2}$ 且风味比例为 1:1:1。

• 性能提升： 虽然 ANTARES 数据在联合分析结果中起到了主要驱动作用，但与 ARCA 数据的结合相比，联合分析的性能比 ANTARES 独立分析提高了约 10%。
• 候选源 MG3 J225517+2409： 该源产生了最低的预试验（pre-trial）p 值，为 $(1.2 \pm 0.4) \cdot 10^{-4}$，对应于 $(3.69 \pm 0.09)\sigma$ 的显著性。
  • 最佳拟合通量：$E^2\phi = (1.39 \pm 0.02) \cdot 10^{-8} \text{ GeV cm}^{-2} \text{ s}^{-1}$。
  • 试验后（post-trial）概率：通过二项公式计算为 $0.013 \pm 0.004$ ($(2.26 \pm 0.13)\sigma$)。基于模拟的方法得出的试验后 p 值为 $0.018 \pm 0.006$ ($(2.12 \pm 0.14)\sigma$)。
  • 考虑了从 1255 个天体中进行初始选择及随后目录缩减的简化试验后估计表明，观察到此类候选源的概率约为 $\approx 0.1$。
• 其他候选源： 本分析还评估了其他已知源，包括 3C403 ($p \approx 1.5 \cdot 10^{-3}$) 和 TXS 0506+056，并为 106 个源的目录提供了 90% CL 上限。
• 天空图： 在 MG3 J225517+2409 附近的星图显示，有一个 ARCA21 事件的误差框包含了该源的位置，其信号似然度比仅背景场景高出 6.3 倍。

意义与主张
本文声称证明了 ANTARES 与 KM3NeT/ARCA 之间进行联合分析的可行性和潜力。研究表明，虽然目前的结果主要由成熟的 ANTARES 数据集驱动，但整合 ARCA 数据可以带来可衡量的灵敏度提升。作者指出，随着 KM3NeT/ARCA 规模和曝光量的持续增长，预计在几年内其性能将达到 ANTARES 的水平，这使得联合数据分析在未来的发现中变得愈发重要。本工作提供了一个通过多种背景建模策略处理系统不确定性的框架，并为广泛的天体物理候选源设定了上限。