[论文解读] A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers

A. Corstanje, S. Buitink+18 more|arXiv (Cornell University)|Jan 1, 2023

Astrophysics and Cosmic Phenomena参考文献 28被引用 3

一句话总结

本文提出一种针对宇宙射线簇射无线电脉冲的高精度插值方法，采用傅里叶级数与三次样条插值，在极坐标网格上重建任意地面位置的完整脉冲时间序列。该方法仅需约200个天线即可实现精确模拟，显著降低计算成本，同时保持亚纳秒级时间分辨率与亚百分之一幅度精度，使其适用于SKA时代簇射研究。

ABSTRACT

Analysis of radio signals from cosmic-ray induced air showers has been shown to be a reliable method to extract shower parameters such as primary energy and depth of shower maximum.The required detailed air shower simulations take 1 to 3 days of CPU time per shower for a few hundred antennas.With nearly 60,000 antennas envisioned to be used for air shower studies at the Square Kilometre Array (SKA), simulating all of these would come at unreasonable costs.We present an interpolation algorithm to reconstruct the full pulse time series at any position in the radio footprint, from a set of antennas simulated on a polar grid.Relying on Fourier series representations and cubic splines, it significantly improves on existing linear methods.We show that simulating about 200 antennas is sufficient for high-precision analysis in the SKA era, including e.g. interferometry which relies on accurate pulse shapes and timings.We therefore propose the interpolation algorithm and its implementation as a useful extension of radio simulation codes, to limit computational effort while retaining accuracy.

研究动机与目标

• 为减少模拟宇宙射线簇射无线电信号的计算负担，特别是针对SKA-Low等大规模阵列。
• 解决在SKA-Low设想的60,000个天线全部模拟时因CPU成本过高而不切实际的问题。
• 开发一种插值方法，以在无线电覆盖范围内保持高精度的脉冲波形与时间信息。
• 通过最少数量的模拟天线，实现对干涉测量与参数重建（如能量、Xmax）的精确模拟。

提出的方法

• 该方法使用模拟天线的极坐标网格来表示簇射无线电信号的覆盖范围。
• 采用傅里叶级数建模总无线电能通量在覆盖范围内的角度依赖性。
• 利用三次样条插值，从模拟天线数据中重建任意位置的电场完整时间序列。
• 该算法估计可靠的截止频率以进行滤波，提升高频成分的重建精度。
• 支持非均匀的径向天线间距，允许在核心附近实现更密集采样，而在远距离区域采用较稀疏分布。
• 该方法通过高密度模拟覆盖范围进行验证，并与先前的径向基函数（RBF）插值方法进行比较。

实验结果

研究问题

• RQ1 能否通过高精度插值方法，仅基于稀疏模拟天线数据，在极低计算成本下重建完整的脉冲时间序列？
• RQ2 与全模拟基准相比，该方法在脉冲波形、幅度与时间信息上的保真度如何？
• RQ3 为实现SKA时代应用所需的亚纳秒级时间分辨率与亚百分之一幅度精度，最少需要多少个天线？
• RQ4 在误差与频率响应方面，该方法与现有插值技术（如径向基函数RBFs）相比表现如何？
• RQ5 非均匀天线间距能否有效提升精度，同时最小化所需模拟数量？

主要发现

• 插值脉冲与模拟脉冲之间的归一化互相关最高达0.9998，表明脉冲波形几乎完全一致。
• 在所有极化方式与核心距离下，脉冲幅度与能量通量均实现了亚百分之一精度的重建。
• 即使在远距离（最高达270 m）仍保持高精度，多数位置的互相关值高于0.99。
• 与先前的RBF方法相比，插值误差显著降低，最大误差低于峰值通量的1%，而此前方法的误差约为3%。
• 仅使用208个战略性布置的天线——尤其在核心200 m以内实现更密集采样——即可获得与全模拟相当的高精度结果。
• 通过精确估计截止频率进行滤波，该方法能够可靠重建高频成分，包括切伦科夫环结构。