Simulating radio emission from particle cascades with CORSIKA 8

TL;DR

This paper introduces the radio emission simulation module in CORSIKA 8, a flexible framework for particle cascade modeling, and validates its performance against existing codes and formalism methods.

Contribution

The paper presents the design, implementation, and validation of a new radio emission module in CORSIKA 8, enabling customizable and multi-formalism simulations of particle showers in complex media.

Findings

CORSIKA 8 radio module matches results from CORSIKA 7 and ZHAireS.

Simulation details like step size affect radio emission accuracy.

Endpoints and ZHS formalisms produce consistent results.

Abstract

CORSIKA 8 is a new framework for simulations of particle cascades in air and dense media implemented in modern C++17, based on past experience with existing codes, in particular CORSIKA 7. The flexible and modular structure of the project allows the development of independent modules that can produce a fully customizable particle shower simulation. The radio module in particular is designed to treat the electric field calculation and its propagation through complex media to each observer location in an autonomous and flexible way. It already allows for the simultaneous simulation of the radio emission calculated with two independent time-domain formalisms, the "Endpoint formalism" as previously implemented in CoREAS and the "ZHS" algorithm as ported from ZHAireS. The design acts as the baseline interface for current and future development for the simulation of radio emission from particle showers in standard and complex scenarios, such as cross-media showers penetrating from air into ice. In this work, we present the design and implementation of the radio module in CORSIKA 8, along with validation studies and a direct comparison of the radio emission from air showers simulated with CORSIKA 8, CORSIKA 7 and ZHAireS. We also present the impact of simulation details such as the step size of simulated particle tracks on radio-emission simulations and perform a direct comparison of the "Endpoints" and "ZHS" formalisms for the same underlying air showers. Finally, we present an in-depth comparison of CORSIKA 8 and CORSIKA 7 for optimum simulation settings and discuss the relevance of observed differences in light of reconstruction efforts for the energy and mass of cosmic rays.