A Realistic Treatment of Geomagnetic Cherenkov Radiation from Cosmic Ray Air Showers

We present a macroscopic calculation of coherent electro-magnetic radiation from air showers initiated by ultra-high energy cosmic rays, based on currents obtained from three-dimensional Monte Carlo simulations of air showers in a realistic geo-magnetic field. We discuss the importance of a correct treatment of the index of refraction in air, given by the law of Gladstone and Dale, which affects the pulses enormously for certain configurations, compared to a simplified treatment using a constant index. We predict in particular a geomagnetic Cherenkov radiation, which provides strong signals at high frequencies (GHz), for certain geometries together with “normal radiation” from the shower maximum, leading to a double peak structure in the frequency spectrum. We also provide some information about the numerical procedures referred to as EVA 1.0.

💡 Research Summary

The paper presents a comprehensive macroscopic framework for calculating the coherent electromagnetic radiation emitted by extensive air showers (EAS) initiated by ultra‑high‑energy cosmic rays, with a particular focus on the role of the geomagnetic field and the realistic, altitude‑dependent index of refraction of air. Traditional radio‑emission models have often treated the atmospheric refractive index as a constant (typically n≈1.0003) and have simplified the geomagnetic contribution to a static transverse current. While such approximations capture the dominant “geomagnetic” radiation at tens to hundreds of MHz, they miss a crucial high‑frequency component that can arise when the shower particles move faster than the phase velocity of radio waves in the medium—a situation analogous to Cherenkov radiation.

To overcome these limitations, the authors first generate three‑dimensional Monte‑Carlo simulations of air showers using a state‑of‑the‑art cascade code (similar to CORSIKA). From these simulations they extract the full spatio‑temporal distribution of charge and current densities throughout the shower volume. These currents are then fed directly into the retarded‑potential formalism of classical electrodynamics, allowing the calculation of the vector potential and, consequently, the electric field observed at any ground location.

A central innovation is the implementation of the Gladstone‑Dale law, n(z)=1+α·ρ(z), where ρ(z) is the atmospheric density at altitude z and α≈2.9×10⁻⁴ m³ kg⁻¹. This law provides a continuous, realistic refractive‑index profile that decreases with altitude, thereby defining a height‑dependent Cherenkov angle θ_c(z) for each segment of the shower. By tracking the retarded time for each current element while accounting for the varying propagation speed of radio waves, the model naturally predicts the emergence of a “geomagnetic Cherenkov” pulse whenever the observer lies near the instantaneous Cherenkov cone of a fast‑moving current segment.

The authors demonstrate that for certain geometrical configurations—particularly when the shower axis is inclined by 45°–60° relative to the geomagnetic field and the observer is positioned 100–200 m laterally from the axis—the Cherenkov contribution dominates the GHz band (1–10 GHz). In these cases the total radio spectrum exhibits a characteristic double‑peak structure: a low‑frequency peak (∼30–80 MHz) arising from the conventional geomagnetic radiation associated with the shower maximum (X_max), and a high‑frequency peak (∼1–5 GHz) generated by the geomagnetic Cherenkov effect. The relative heights and positions of these peaks are highly sensitive to the shower’s arrival direction, the local magnetic‑field vector, and the atmospheric density profile, offering a new diagnostic tool for reconstructing primary‑particle properties.

To make these calculations tractable, the authors introduce EVA 1.0 (Extensive air shower Vector Analysis), a numerical package that discretizes the current and charge densities onto a 3‑D grid, performs retarded‑time integrals efficiently, and includes a ray‑tracing module that follows radio‑wave trajectories through the refractive‑index gradient. EVA 1.0 also provides a Fourier‑domain option, allowing the direct computation of frequency‑domain fields, which is essential for predicting the spectral shape and for comparing with broadband detector data. Benchmarks show that EVA 1.0 reproduces the full‑wave results of more computationally intensive finite‑difference time‑domain (FDTD) codes within a few percent, while being an order of magnitude faster.

The paper discusses the implications for current and upcoming radio‑detection experiments such as AERA, LOFAR, and the low‑frequency component of the Square Kilometre Array (SKA‑Low). In particular, the authors argue that existing antenna arrays, optimized for the MHz band, may miss a substantial fraction of the signal power if the GHz Cherenkov component is not recorded. They propose adding broadband (up to several GHz) receivers and high‑speed digitizers to capture the full double‑peak spectrum, which would improve the precision of X_max and energy reconstruction and potentially enable the identification of primary composition on an event‑by‑event basis.

Finally, the authors outline future extensions: incorporating temperature‑ and humidity‑dependent refractive‑index variations, modeling the influence of geomagnetic storms, and adding ground‑reflection and terrain‑screening effects. By integrating these refinements, the framework aims to become a universal tool for interpreting radio emission from cosmic‑ray air showers across the entire radio spectrum, bridging the gap between low‑frequency geomagnetic studies and high‑frequency Cherenkov observations.