Uncertainties in the Estimation of Air Shower Observables from Monte Carlo Simulation of Radio Emission

Abstract

The detection of extensive air showers (EAS) induced by cosmic rays via radio signals has undergone significant advancements in the last two decades. Numerous ultra-high energy cosmic ray experiments routinely capture radio pulses in the MHz to GHz frequency range emitted by EAS. The Monte Carlo simulation of these radio pulses is crucial to enable an accurate reconstruction of the primary cosmic ray energy and to infer the composition of the primary particles. In this work, a comprehensive comparison of the predicted electric field in EAS simulated with CoREAS and ZHAireS was conducted to estimate the systematic uncertainties arising from the use of different simulation packages in the determination of two key shower observables namely, the electromagnetic energy of the EAS and the depth of maximum development ($X_{\max }$). For this comparison, input parameters and settings as similar as possible were used in both simulations, along with the same realistic atmospheric refractive index depending on altitude, which is crucial for the prediction of radio emission properties of EAS. In addition, simulated EAS with very similar values of depth of maximum development were selected. Good agreement was found between CoREAS and ZHAireS, with discrepancies in the dominant electric field components generally remaining below 10\% across the frequency range of a few MHz to hundreds of MHz, relevant for most radio detection experiments, translating into uncertainties in the determination of energy below $5\%$ and $\simeq 10g/c{\mathrm{m}}^2$ in $X_{\max }$. Our work underscores the need for further studies to clarify their origin and impact on $X_{\max }$ inference in composition analyses.

24 pages, 11 figures, version accepted in Astroparticle Physics

References in corpus (13)

• Radio detection of Cosmic-Ray Air Showers and High-Energy Neutrinos
• A method for high precision reconstruction of air shower Xmax using two-dimensional radio intensity profiles
• Ultra-High-Energy Cosmic Rays: The Intersection of the Cosmic and Energy Frontiers
• Polarized radio emission from extensive air showers measured with LOFAR
• Measuring a Cherenkov ring in the radio emission from air showers at 110-190 MHz with LOFAR
• Depth of shower maximum and mass composition of cosmic rays from 50 PeV to 2 EeV measured with the LOFAR radio telescope
• Radio Measurements of the Depth of Air-Shower Maximum at the Pierre Auger Observatory
• Demonstrating Agreement between Radio and Fluorescence Measurements of the Depth of Maximum of Extensive Air Showers at the Pierre Auger Observatory
• The SLAC T-510 experiment for radio emission from particle showers: detailed simulation study and interpretation
• Density and magnetic intensity dependence of radio pulses induced by energetic air showers
• Simulating radio emission from particle cascades with CORSIKA 8
• A LOFAR-style reconstruction of cosmic-ray air showers with SKA-Low
• Characterization of Atmosphere-Skimming Cosmic-Ray Showers in High-Altitude Experiments