14–24 Jul 2025
CICG - International Conference Centre - Geneva, Switzerland
Europe/Zurich timezone

Fast and accurate simulation of the radio emission from cosmic ray air showers using template synthesis

17 Jul 2025, 13:20
15m
Room 8

Room 8

Talk Cosmic-Ray Indirect CRI

Speaker

Mitja Desmet

Description

In order to interpret the radio data from extensive air showers detectors, we rely on accurate simulations. The state-of-the-art simulation frameworks use Monte-Carlo techniques which pose computational challenges. This is a limiting factor for the next generation of radio arrays, for example the upcoming Square Kilometer Array (SKA), which will have orders of magnitude more antennas than the current arrays. Therefore, we developed template synthesis, a novel, fast and accurate forward model to calculate the radio emission, which achieves the same accuracy as full microscopic simulations in only a fraction of the time. This speeds up the simulation-heavy reconstruction techniques used in Auger and LOFAR by orders of magnitude. Moreover, it can revolutionize the field by allowing efficient production of simulations for more advanced techniques such as interferometry and reconstruction of the longitudinal shape of the shower. Since template synthesis is a fully differentiable model, it can even be used for machine learning or information field theory applications

Our method synthesises the radio emission using a microscopically (Monte-Carlo) simulated origin shower and a set of semi-analytical scaling relations. These relations were extracted from a large set of CoREAS simulations, thus capturing the behaviour of the microscopic models. We divide the atmosphere into slices of constant atmospheric depth and synthesise the emission from each slice separately. In this process we correct for the shower age inside each slice, which we found to be one of the crucial parameters determining the radio emission. The computation time of the synthesis process is negligible compared to the runtime of the microscopic simulation. Crucially, when comparing the synthesised traces to CoREAS simulations, the amplitudes are typically within 5% of each other.

In this contribution we present the complete framework, which can be used for showers of all geometries without any modifications, addressing an important limitation in the previous iterations. It can also account for different experimental conditions, such as the atmosphere and magnetic field, as well as the influence of the air shower geometry. We show that the scaling relations are indeed universal and can be used in many different conditions, without losing its accuracy.

Since the template synthesis framework is now ready to be used in analyses, we are happy to make it available to the community as a Python package. It will be made public during the conference, complete with the necessary parameters to use it.

Author

Co-authors

Keito Watanabe Stijn Buitink (Vrije Universiteit Brussel (VUB)) Tim Huege

Presentation materials