Abstract:We present detailed microscopic simulations of high-energy cosmic-ray air showers penetrating high-altitude ice layers that can be found at the polar regions. We use a combination of the CORSIKA Monte Carlo code and the Geant4 simulation toolkit, and focus on the particle cascade that develops in the ice to describe its most prominent features. We discuss the impact of the ice layer on the total number of particles in function of depth of the air shower, and we give a general parameterization of the charge distribution in the cascade front in function of Xmax of the cosmic ray air shower, which can be used for analytical and semi-analytical calculations of the expected Askaryan radio emission of the in-ice particle cascade. We show that the core of the cosmic ray air shower dominates during the propagation in ice, therefore creating an in-ice particle cascade strongly resembling a neutrino-induced particle cascade. Finally, we present the results of microscopic simulations of the Askaryan radio emission of the in-ice particle cascade, showing that the emission is dominated by the shower core, and discuss the feasibility of detecting the plasma created by the particle cascade in the ice using RADAR echo techniques.

What problem does this paper attempt to address?

The problem that this paper aims to solve is how to simulate the particle - cascade process generated by high - energy cosmic rays in the high - altitude ice layer. Specifically, the researchers use the CORSIKA Monte Carlo code and the Geant4 simulation toolkit to simulate the development of these particle cascades in ice, with the focus on describing the main characteristics of particle cascades in ice. They explored the influence of the ice layer on the total number of particles in the air shower and gave a general parameterization of the charge distribution at the front edge of the particle cascade, which can be used for the analysis and semi - analytical calculation of the Askaryan radio emission expected to be generated by the particle cascade in ice. In addition, the study also demonstrated the possibility of detecting the plasma generated by the particle cascade in ice through the radar echo technique. ### Research Background and Objectives With the IceCube Neutrino Observatory confirming the existence of high - energy cosmic neutrinos in the TeV - PeV energy range, a new window for exploring the universe has been opened. However, at energies exceeding several PeV, even a detection volume on the cubic - kilometer scale like IceCube is difficult to collect enough events. Given the long attenuation length of radio - wave signals, the detection modules of radio neutrino observatories are usually about 1 kilometer apart and can efficiently cover a larger detection volume. Therefore, using radio detection techniques to search for interacting neutrinos may become a valuable means of expanding to the extremely high - energy region. ### Simulation Setup and Methods To study the propagation of particle cascades in ice induced by cosmic rays and their related radio - detection characteristics, the researchers developed a dedicated simulation framework, combining the CORSIKA Monte Carlo code and the Geant4 simulation toolkit. CORSIKA is used to generate proton - induced cosmic - ray air showers with different primary - particle energies and zenith angles, providing the positions and momenta of all particles in the shower. Geant4 is then used to simulate the propagation of these particles in ice, taking into account the density gradient of the Taylor Dome ice cap in Antarctica. ### Simulation Results - **Deposited Energy**: The study shows that the average energy of the electromagnetic component is lower than the critical energy (about 80 MeV) for cascade generation in ice at a distance of several meters. Therefore, in this region, the cascade will gradually disappear within a few radiation lengths. At a small radius close to the core, the average energy of the electromagnetic component is much higher than the critical energy, and the cascade will continue to develop until it reaches the maximum value. - **Longitudinal Profile**: Compared with the propagation in air, the electromagnetic part in ice is hardly affected by the air - ice interface, while the number of hadrons and muons is reduced. - **Lateral Charge Distribution**: The particle distribution within the front edge of the particle cascade is affected by multiple Coulomb scattering and mainly depends on the density of the medium. Through simulation, the researchers found that these distributions can be well described by an analytical expression that depends on the maximum depth \(X_{\text{max}}\) of the air shower. ### Application Prospects - **Askaryan Radio Emission**: Since the properties of particle cascades in ice are very similar to those of cascades induced by high - energy neutrinos, it is expected that both will produce similar Askaryan radio emissions. If not identified, this signal may become a strong background for Askaryan detectors in ice to search for high - energy neutrino interactions; but if identified, it can also serve as a natural principle verification for the Askaryan method to detect neutrino - induced particle cascades in ice. - **Radar Echo Technique**: The study also explored the possibility of detecting the plasma generated by the particle cascade in ice through the radar echo technique, which has important application value for future research.

Simulation of in-ice cosmic ray air shower induced particle cascades

Simulation of radio signals from cosmic-ray cascades in air and ice as observed by in-ice Askaryan radio detectors

In-ice Askaryan Emission from Air Showers: Implications for Radio Neutrino Detectors

Radio signatures of cosmic-ray showers with deep in-ice antennas

Estimation of aperture of the Tunka-Rex radio array for cosmic-ray air-shower measurements

Simulating radio emission from air showers with CORSIKA 8

Simulating radio emission from particle cascades with CORSIKA 8

Parametrization of the Relative Amplitude of Geomagnetic and Askaryan Radio Emission from Cosmic-Ray Air Showers using CORSIKA/CoREAS Simulations (ICRC2021)

RDSim, a fast and comprehensive simulation of radio detection of air showers

Analytic calculation of radio emission from parameterized extensive air showers, a tool to extract shower parameters

Prospects for a radio air-shower detector at South Pole

Radio wavefront of very inclined extensive air-showers: a simulation study for extended and sparse radio arrays

Simulations of cross media showers with CORSIKA 8

Estimation of $X_\mathrm{max}$ for air showers measured at IceCube with elevated radio antennas of a prototype surface station

A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers

Particle content of very inclined air showers for radio signal modeling

Towards a cosmic ray composition measurement with the IceAct telescopes at the IceCube Neutrino Observatory

Characterization of Atmosphere-Skimming Cosmic-Ray Showers in High-Altitude Experiments

Macroscopic approach to the radar echo scatter from high-energy particle cascades