[论文解读] The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 34th International Cosmic Ray Conference (ICRC 2015)

IceCube Collaboration, M. G. Aartsen+892 more|arXiv (Cornell University)|Jan 1, 2015

Astrophysics and Cosmic Phenomena参考文献 2被引用 7

一句话总结

本文结合皮埃尔·让宇宙射线观测站（Pierre Auger Observatory）与望远镜阵列（Telescope Array）的超高能宇宙射线（UHECRs）数据，以及冰立方中微子观测站（IceCube Neutrino Observatory）的高能中微子候选事件，开展联合分析。通过堆叠法与交叉相关法，搜索UHECR到达方向与中微子事例之间的相关性，发现级联事例的预试验p值为2.7×10⁻⁴——在考虑宇宙射线各向异性后，其显著性约为2.4σ，虽尚未达到统计显著性水平。

ABSTRACT

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

研究动机与目标

• 研究皮埃尔·让宇宙射线观测站与望远镜阵列探测的超高能宇宙射线（UHECRs）与冰立方中微子观测站高能中微子候选事件之间潜在的相关性。
• 评估观测到的UHECR到达方向各向异性（如望远镜阵列热点）是否可能与来自相同源的中微子通量相关联。
• 评估堆叠分析对点状中微子源的探测灵敏度，结合能量与角度分辨率模型，这些源与UHECR源相关。
• 确定在不同能量阈值与磁偏转假设下，实现5σ发现潜力所需的每源中微子通量。
• 比较不同角度距离假设（D = 3°, 6°, 9°）与事件类型（轨迹与级联）的结果，以评估一致性与显著性。

提出的方法

• 采用无参量似然分析方法，基于已知的UHECR到达方向对中微子信号事例进行堆叠，将信号概率密度函数（PDF）建模为角度距离与能谱的函数。
• 使用似然函数 lnL(nsν,γ) = ∑ᵢ ln[nsν/Nν × Si(γ,Ei) + (1−nsν/Nν) × Bi]，其中 nsν 为总信号事例数，γ 为谱指数，Bi 表示大气背景。
• 通过 σj = √(σ²_M_D + σ²_exp) 建模UHECR的角度不确定性，其中 σ_exp = 0.9°（奥古尔）或 1.5°（TA），σ_M_D = D × 100 EeV / ECR 表示磁偏转。
• 设定能量阈值 Eth = 85 EeV 以限制UHECR源的数量并提高灵敏度，基于实现5σ发现所需通量的模拟结果。
• 计算预试验与事后试验p值以评估显著性，对不同角度距离D进行多重检验校正。
• 将结果与各向同性中微子通量及各向同性UHECR方向进行比较，以隔离宇宙射线各向异性对显著性的影响。

实验结果

研究问题

• RQ1 皮埃尔·让宇宙射线观测站与望远镜阵列探测的UHECR到达方向与冰立方高能中微子事例之间是否存在统计显著的相关性？
• RQ2 该潜在相关性的显著性如何依赖于UHECR磁偏转的假设（即 σ_M_D = D × 100 EeV / ECR）？
• RQ3 在堆叠分析中，实现5σ发现潜力所需的每源中微子通量是多少？该值如何随能量阈值 Eth 变化？
• RQ4 观测到的UHECR到达方向各向异性（如望远镜阵列热点）是否会影响潜在中微子相关性的显著性？
• RQ5 冰立方数据集中，高能轨迹与高能级联事例的结果有何差异？

主要发现

• 对高能级联事例的堆叠分析在 D = 6° 时得到预试验p值为2.7×10⁻⁴，事后试验p值为8×10⁻⁴，表明可能存在有趣信号。
• 在考虑UHECR到达方向各向异性（如TA热点）后，级联相关性的显著性上升至约2.4σ。
• 对于相同的角度距离 D = 3°，最小事后试验p值为25%，拟合得到的超额事例数约为123个，谱指数 γ = −3.24。
• 在 D = 6° 的分析中，事后试验p值大于50%，表明无显著超额。
• 实现5σ发现潜力所需的每源中微子通量在 Eth = 85 EeV 时最小，此时灵敏度与源数量达到最佳平衡。
• 交叉相关分析推断的角度尺度（≈22°）与磁偏转尺度（~10°）及实验不确定性（级联约15°，宇宙射线约1°）一致，支持结果的内部一致性。