Search for neutrino emission from blazar $γ$-ray flares accounting for possible neutrino time delays

技术摘要：考虑可能存在的中微子时间延迟对耀变体 $\gamma$ 射线耀发中中微子发射的搜索

问题陈述
尽管 IceCube 中微子天文台已经证实了高能天体中微子（$E_\nu \gtrsim 1$ TeV）的存在，但其弥散通量所贡献的具体源头在很大程度上仍未被识别，尽管已有证据将特定事件与耀变体（例如 TXS 0506+056）和活动星系（例如 NGC 1068）联系起来起来。以往针对耀变体与 IceCube 中微子相关性的搜索通常假设中微子到达时间与主要的电磁（特别是 $\gamma$ 射线）耀发在时间上是重合的。然而，理论考量表明这种重合可能并不成立。如果中微子是通过 $p\gamma$ 相互作用产生的，那么喷流框架内的质子加速时间（$t'_{acc}$）和质子能量损失时间（$t'_{p\gamma}$）可能比负责电磁耀发的电子时间尺度长出几个数量级。因此，高能中微子（$\gtrsim 100$ TeV）可能会在相关的 $\gamma$ 射线耀发峰值之后数月甚至数年才到达。本文旨在解决缺乏考虑此类潜在时间滞后情况的搜索问题。

方法论
作者在 IceCat-1 公开目录中对高能中微子警报与两个不同的耀变体源样本进行了交叉匹配分析，扫描了广泛的喷流框架时间延迟（$t'_{delay}$) 范围内的相关性。

1. 数据集：

   • 中微子： 来自 IceCat-1 目录的 348 个警报（$E_\nu \gtrsim 100$ TeV），排除了与 IceTop 信号重合的警报。
   • 源列表：
     • H21+： 一个包含 2 94 个来自 MOJAVE 数据库（Homan et al. 2021）的 AGN 样本，具有已知的多普勒因子（$D$）、红移（$z$）和 Fermi-LAT 光变曲线。
     • R24+： 一个包含来自 CGRaBS 目录（Rodrigues et al. 2024b）的 295 个耀变体的样本，其多普勒因子通过能谱建模推导得出，并与 Fermi-LAT 光变曲线进行了交叉匹配。
   • 光变曲线： 来自 Fermi-LAT 光变曲线库（LCR）的每月分箱光子通量光变曲线。

2. 统计框架：

   • 空间关联： 如果源落在扩大的角不确定性区域内（包含系统不确定性 $\Delta = 0.78^\circ$），则将中微子与耀变体关联起来。
   • 权重方案： 对于每个警报，根据警报的显著性、针对大角不确定性的惩罚因子以及时间因子 $G_{\nu a}$ 计算权重 $w_\nu(t'_{delay})$。
   • 时间因子： 该因子将中微子到达时间（$t_\nu$）与耀变体在到达之前的全局 $\gamma$ 射线光变曲线最大值（$t^{\gamma \max}_a$）联系起来。该因子在匹配预期观测者框架延迟（$t_{delay} = t'_{delay}(1+z)/D$）的箱体大小范围内分配最大权重，并对偏差应用高斯惩罚。
   • 检验统计量 (TS)： 全局 TS 定义为所有警报个体权重的总和。
   • 显著性检验： 将真实数据的 TS 与通过对 49,999 个模拟数据集（其中中微子的赤经和到达时间被随机化）生成的 TS 分布进行比较，计算预试验 $p$ 值，随后进行试验校正（通过扫描 $t'_{delay}$）。

主要贡献

• 新颖的搜索策略： 本研究代表了首次尝试专门考虑由缓慢质子加速或能量损失引起的通用喷流框架时间延迟（$t'_{delay}$），来搜索 IceCube 中微子与耀变体 $\gamma$ 射线耀发之间关联的研究。
• 双重源样本： 研究利用了两个独立的耀变体目录，这些目录采用了不同的多普勒因子确定方法（亮度温度论证对比 SED 建模），以测试潜在相关性对喷流运动学不确定性的稳健性。
• 系统延迟扫描： 分析扫描了从 1 个月到 $\sim 3.5 \times 10^4$ 天的 $t'_{delay}$，覆盖了中微子产生可能滞后于电磁辐射的理论范围。

结果

• 统计显著性： 搜索结果显示，在 IceCat-1 中微子与 H21+ 或 R24+ 源列表中的 AGN 之间不存在统计显著的相关性证据。
  • 对于 H21+ 列表，最小预试验 $p$ 值为 0.037（在 $t'_{delay} \approx 7.7 \times 10^3$ d 处）。
  • 对于 R24+ 列表，最小预试验 $p$ 值为 0.026（在 $t'_{delay} \approx 1.9 \times 10^3$ d 处）。
  • 经过试验校正后，这些最小值分别对应于 0.14 (1.06$\sigma$) 和 0.11 (1.21$\sigma$)，与零假设一致。
• 顶级候选对象： 最显著的关联是由警报 IC120523B（与 3C 454.3 的 2010 年耀发相关）和 IC150926A（与 3C 279 的耀发相关）驱动的。虽然这些特定事件显示出与预期延迟的时间对齐，但不足以建立整个数据集的全局相关性。
• 多普勒因子的差异： 两个源列表的最优延迟时间存在显著差异（达 4 倍），这主要是由于 H21 和 R24 方法在估计 3C 454.3 和 3C 279 的多普勒因子时存在差异所致。

意义与主张
作者谨慎地得出结论，即使考虑到时间延迟，其结果也不支持 IceCat-1 中微子与所研究的特定耀变体族群之间存在显著相关性。他们提出了几种可能并非互斥的原因，解释了这些零结果：

1. 低效率： 耀变体可能是低效的 $\sim 100$ TeV 中微子源，仅贡献了 IceCube 通量的极小部分（约 0.5%）。
2. 能量尺度： 如果耀变体是高效的中微子源，其发射峰值可能位于更高的能量水平（$\gtrsim 10$ PeV），超出了当前 IceCube 警报样本的主要灵敏度范围。
3. 模型假设： 假设存在通用的喷流框架时间延迟（$t'_{delay}$）可能是无效的，因为不同源的加速机制和磁场可能存在差异。
4. 解耦： 质子加速和能量耗散在空间或时间上可能与电子过程解耦，这意味着中微子的产生与观测到的 $\gamma$ 射线耀发之间不存在因果联系。

本文指出，未来的改进可能需要更大的数据集（例如即将到来的 IceCat-2 和 CAZ 目录）以及包含射电光变曲线，因为已知射电光变曲线会滞后于 $\gamma$ 射线耀发，并且可能更好地与中微子产生的时间尺度相关联。