Energy estimation of cosmic rays by air shower radio signals

技术摘要：基于空气簇射无线电信号的宇宙射线能量估计

问题陈述
确定宇宙射线的一级能量对于理解其起源、加速和传播机制至关重要。虽然由宇宙射线引发的广延大气簇射（EAS）会破坏一级粒子的直接信息，但次级粒子和无线电发射提供了重建路径。先前的实验（如 AERA、Tunka-Rex、LOPES）利用无线电信号的横向分布函数（LDF），将特定距离处信号强度与一级能量进行关联。然而，这些方法通常需要针对每个站点的几何结构识别一个最优距离，并且可能对簇射核心位置较为敏感。位于伊朗塞姆南的 SURA 实验是一个由四个天线组成的无线电阵列，该实验已成功重建了到达方向和核心位置，但需要一种鲁棒的方法来估计一级能量，以最大限度地减少对核心位置变化的敏感性。

方法论
作者提出了一种基于将实验无线电信号与使用 CoREAS 扩展的 CORSIKA 进行高保真模拟进行对比的重建方法。该方法论的核心步骤如下：

1. 模拟设置：

   • SURA 配置： 模拟旨在匹配 SURA 实验的具体参数，包括其四个对数周期偶极天线（LPDA）、局部磁场分量（$B_x = 39.53 \, \mu T$, $B_z = 28.1 \, \mu T$）以及海拔高度（1130 m）。
   • 密集参考阵列： 模拟了一个包含 12,342 个天线、间距为 1 米的虚拟密集阵列，以生成无线电信号强度的高分辨率参考。该阵列作为比较的基准。
   • 事件生成： 模拟了一级能量范围在 $10^{17}$ eV 至 $7 \times 10^{17}$ eV 之间的质子诱导簇射。一级到达方向设定为天顶角 $35^\circ$ 和方位角 $40^\circ$（基于 SURA 的峰值检测率），并进行了针对 $60^\circ$ 和 $65^\circ$ 天顶角的额外模拟，以测试倾斜效应。

2. 数据校准：

   • 研究采用了两阶段校准过程：实时 FPGA 滤波（高通、带通、陷波滤波器）和后处理（到达方向重建、频率响应应用、小波去噪以及信号上采样至 0.2 ns 分辨率）。

3. 比例因子 ($C_{ij}$) 方法：

   • 定义了一个比例因子 $C_{ij}$，它是模拟信号的平均峰值强度（来自密集阵列 $\bar{A}_{sim}$）与 SURA 配置信号的峰值强度（$\bar{A}_{exp}$）之比：
     $$C_{ij} = \frac{\bar{A}_{sim}}{\bar{A}_{exp}}$$
   • 该因子是通过最小化比较实验（模拟 SURA）和密集阵列信号在整个天线阵列上的 $\chi^2$ 函数而得出的。
   • 该方法假设，虽然绝对信号强度取决于核心位置，但稀疏的 SURA 阵列与密集参考阵列之间的比例（比例因子）表现出与一级能量的特定相关性。

主要贡献

• 核心位置无关性： 论文证明了比例因子 $C_{ij}$ 在很大程度上独立于宇宙射线簇射的核心位置。即使对于远离阵列中心的簇射核心位置（在不同测试中高达 $\pm 30$ m 和 $\pm 125$ m），对于给定的能量，比例因子仍保持有效恒定。
• 能量相关性： 作者建立了比例因子与一级能量之间的明确反相关关系。随着一级能量的增加，比例因子减小。这种关系通过一个反一次函数进行建模：
  $$C_{ij} = y_0 + \frac{a}{E}$$
  其中 $E$ 是以 $10^{17}$ eV 为单位归一化的一级能量，$y_0$ 和 $a$ 是依赖于到达方向的拟合常数。
• 对到达方向的鲁棒性： 该方法在倾斜簇射（天顶角为 $60^\circ$ 和 $65^\circ$）上进行了测试。在这些角度下，$C_{ij}$ 与能量之间的反相关关系依然成立，尽管拟合参数（$a$ 和 $y_0$）随天顶角的变化略有不同。

结果

• 重建精度： 使用导出的比例因子和拟合的反函数，重建了模拟宇宙射线的一级能量。该方法实现了约 11% 的最大重建误差（在 $35^\circ$ 天顶角下，于 $5 \times 10^{17}$ eV 处观察到）。
• 误差分析： 计算得出相对误差 ($\Delta E / E$) 在某些能量水平下低至 2%（例如 $2 \times 10^{17}$ eV），并在测试的能量谱范围内普遍保持在 1%–11% 之间。
• 倾斜簇射： 对于更倾斜的簇射（$60^\circ$ 和 $65^\circ$ 天顶角），该方法仍然有效，尽管精度有所变化。在 $65^\circ$ 角度下，最低精度出现在 $1 \times 10^{17}$ eV 和 $5 \times 10^{17}$ eV；而在 $60^\circ$ 角度下，最大误差出现在 $7 \times 10^{17}$ eV。

意义与主张
论文声称，该技术通过降低对簇射核心位置的敏感性，提供了优于以往方法的显著优势，而核心位置是无线电重建中常见的误差来源。通过利用密集参考阵列来归一化稀疏的实验数据，该方法提供了一个可靠的一级能量估计器，而无需精确确定每个特定事件距离核心的最优距离。

作者得出结论，该方法对垂直和倾斜的簇射均有效，并且可以在实际的 SURA 实验数据收集完成后应用于其数据。他们指出，未来的工作将涉及将此技术应用于真实的 SURA 实验数据，并扩展模拟以包含具有更远核心位置和更多样化到达方向的簇射，从而进一步验证该方法的鲁棒性。该研究并不声称取代现有方法，而是提供了一种专门针对 SURA 阵列的几何结构和约束条件定制的补充方法。