Hybrid radio and particle detection of air showers: potential for ultra-high-energy photon identification

技术摘要：用于超高能光子鉴别的混合射电与粒子簇射探测技术

问题陈述
识别超高能（UHE）光子（能量 $> 10^{17}$ eV）是宇宙射线物理学中的一个关键挑战，因为这些粒子是定位超高能宇宙射线源并约束其起源的重要信使。然而，这一搜索过程受到主要的强子宇宙射线（主要是质子）背景的阻碍。传统的鉴别方法依赖于簇射发展的差异，例如缪子含量（质子诱导的簇射富含缪子，而光子诱导的簇射则缺乏缪子）以及最大发展深度（$X_{max}$）。

现有的搜索工作主要集中在近垂直方向的簇射，使用粒子探测器（闪烁计数器或水切连科夫坦克），并通常辅以荧光探测器或地下缪子探测器。然而，对于极高倾角的簇射（天顶角 $\theta > 60^\circ$），电磁组分和强子组分在到达地面之前已被大气吸收。因此，基于地面的粒子探测器仅能记录存留的缪子组分。虽然射电探测在处理此类倾斜簇射方面已取得进展，但区分光子诱导事件与质子诱导事件仍然困难，因为两者都会产生射电信号，且对于光子而言，其粒子信号往往过于稀疏，难以有效触发标准粒子探测器。

方法论
本研究提出了一种结合射电天线与地面闪烁计数器的混合探测策略，用于鉴别倾斜的光子诱导与质子诱导簇射。该研究利用基于 GRANDProto300 原型机（位于戈壁沙漠的一个拥有 300 根天线的射电阵列）的虚拟实验设置，并辅以放置在每根天线位置的类似 Telescope Array (TA) 的闪烁计数器。

该方法包含以下步骤：

1. 模拟： 使用 CORSIKA (v7.7550) 进行质子和光子引发的扩展大气簇射模拟，并使用 CoREAS 进行射电发射模拟。模拟覆盖了从 $10^{17}$ 到 $10^{18.5}$ eV 的能量以及从 $70^\circ$ 到 $85^\circ$ 的天顶角。分析重点在于极高倾角簇射的特定物理特性，即在地磁辐射和地磁同步辐射起重要作用的区域，且此时粒子足迹分布在较大面积上（$\sim 30$ km）。
2. 触发准则： 要求事件至少触发 4 根射电天线（在 50–200 MHz 范围内峰峰值振幅 $> 75$ $\mu$V/m），并且对于混合分析，至少触发 3 个闪烁计数器（能量沉积 $> 3$ MIPs）。
3. 鉴别观测量： 针对触发事件提取了两个关键观测量：
   • $RMS_{radio,tot}$： 所有被触发天线中射电信号轨迹的总均方根值。
   • $E_{dep,tot}$： 闪烁计数器中记录的总能量沉积。
     光子诱导的簇射预计具有更高的射电信号（由于其更发达的电磁组分），但与同能量的质子诱导簇射相比，其粒子能量沉积显著较低（由于缺乏缪子）。
4. 分类： 在 $\log_{10}(E_{dep,tot}) - \log_{10}(RMS_{radio,tot})$ 空间内训练逻辑回归分类器。为了应对剧烈下降的宇宙射线能谱以及光子的稀有性，质子事件被赋予了 10 倍的权重。分析在 16 个能量和天顶角区间内进行，以优化不同簇射几何结构的性能。
5. 优化： 优化分类阈值（$p_{\gamma,cut}$）以最小化背景污染（质子被误认为光子），同时保持可接受的信号效率。研究还通过测试带有或不带有分箱以及不同权重方案的配置，调查了重建分辨率（能量和天顶角）的影响。

主要贡献与结果

• 混合触发效率： 研究表明，虽然在所考虑的能量范围内，射电触发效率对于质子和光子都接近 100%，但条件粒子触发效率（给定射电触发的情况下）对于质子显著更高。这种差异源于质子簇射产生了能够到达地面的富含缪子的组分，而光子簇射由于粒子密度较低，往往无法触发闪烁计数器。
• 鉴别性能： 混合方法实现了两种主要类型之间的有效分离。
  • 对于能量高于 $10^{17.8}$ eV 且天顶角低于 $82^\circ$ 的情况，背景污染（质子泄漏）可以降低到 $10^{-4}$ (0.01%) 以下。
  • 在最低能量和最高天顶角处观察到的最大背景污染为 2.6%，此时粒子足迹最为稀疏。
  • 确定了一个最优分类阈值 $p_{\gamma,cut} = 0.9$，它在约 $0.12\%$ 的背景污染与随能量和角度从 0% 到 78% 变化的信号效率之间取得了平衡。
• 通量限制： 通过将此方法应用于假设的 GRANDProto300 混合阵列运行一年的情况，作者推导出了积分光子通量的 95% 置信水平上限。
  • 在 0.3 至 3 EeV 的能量范围内，预测的限制范围为 1.2 至 2.6 km$^{-2}$ sr$^{-1}$ yr$^{-1}$。
  • 这些限制与 Pierre Auger 观测站、KASCADE-Grande 和 Telescope Array 的现有约束相比具有竞争力，特别是在其他实验灵敏度有限的高倾角区域（$\theta > 70^\circ$）。
  • 将阵列扩展到完整的 GRANDProto300 规模（约为填充阵列天线数量的 3 倍）将使这些限制提高约 3 倍。

意义与主张
论文声称，混合射电与粒子探测提供了一种可扩展且具有竞争力的超高能光子鉴别方法，专门针对“未知领域”的极高倾角簇射（$\theta > 70^\circ$）。通过利用光子独特的贫缪子特性与可探测电磁组分的对比，该方法为现有的垂直簇射搜索提供了互补的探测通道。

作者谦虚地指出，由于缺乏用于验证的真实数据，本研究受限于模拟。他们强调，未来的改进方向可能包括：

• 引入额外的观测量，例如激活的闪烁计数器数量（$N_{scint}$），该指标显示出质子与光子的清晰区别。
• 利用多元分析（MVA）或机器学习技术。
• 改进重建方法（能量和天顶角），目前这些方法依赖于真实的模拟值而非重建的估计值。
• 将粒子足迹分析扩展到更大的距离，目前这受到针对低倾角簇射优化的去稀疏算法（dethinning algorithms）的限制。

最终，这项工作表明，像 GRAND 这样计划中的大规模射电天文台，如果配备了地面粒子探测器，可以显著推进对超高能光子的搜索，并为宇宙射线起源提供新的约束。