Proof of principle for template synthesis approach for the radio emission from vertical extensive air showers

Mitja Desmet,Stijn Buitink,Tim Huege,David Butler,Ralph Engel,Olaf Scholten

DOI: https://doi.org/10.1016/j.astropartphys.2023.102923

2024-01-09

Abstract:The radio detection technique of cosmic ray air showers has gained renewed interest in the last two decades. While the radio experiments are very cost-effective to deploy, the Monte-Carlo simulations required to analyse the data are computationally expensive. Here we present a proof of concept for a novel way to synthesise the radio emission from extensive air showers in simulations. It is a hybrid approach which uses a single microscopic Monte-Carlo simulation, called the origin shower, to generate the radio emission from a target shower with a different longitudinal evolution, primary particle type and energy. The method employs semi-analytical relations which only depend on the shower parameters to transform the radio signals in the simulated antennas. We apply this method to vertical air showers with energies ranging from $10^{17}$ eV to $10^{19}$ eV and compare the results with CoREAS simulations in two frequency bands, namely the broad [20, 500] MHz band and a more narrow one at [30, 80] MHz. We gauge the synthesis quality using the maximal amplitude and energy fluence contained in the signal. We observe that the quality depends primarily on the difference in $X_{\text{max}}$ between the origin and target shower. After applying a linear bias correction, we find that for a shift in $X_{\text{max}}$ of less than 150 $\text{g}/\text{cm}^2$ , template synthesis has a bias of less than 2% and a scatter up to 6%, both in amplitude, on the broad frequency range. On the restricted [30, 80] MHz range the bias is similar, but the spread on amplitude drops down to 3%. These fluctuations are on the same level as the intrinsic scatter we observe in Monte-Carlo ensembles. We therefore surmise the observed scatter in amplitude to originate from intrinsic shower fluctuations we do not explicitly account for in template synthesis.

High Energy Astrophysical Phenomena

What problem does this paper attempt to address?

### Problems the paper attempts to solve

The paper aims to solve the problem of high computational cost when using traditional Monte Carlo simulation methods to analyze the radio emission data of extensive air showers (EAS). Specifically, the current state - of - the - art Monte Carlo - based simulation codes (such as CoREAS and ZHAireS) are very time - consuming when calculating the radio emission of air showers, which limits the number of simulations that can be generated. To overcome this limitation, the authors propose a new method - template synthesis.

### Main objectives of the template synthesis method

1. **Reduce computational cost**: Generate radio emission signals of target showers with different longitudinal evolutions, primary particle types and energies by using a single microscopic Monte Carlo simulation (called the origin shower).
2. **Improve simulation efficiency**: Use semi - analytic relations, which only depend on shower parameters, to transform radio signals in simulated antennas.
3. **Verify the effectiveness of the method**: Evaluate the performance of the template synthesis method in two frequency ranges (20 - 500 MHz and 30 - 80 MHz), especially the deviation and scatter in terms of maximum amplitude and energy flux, by comparing with CoREAS simulation results.

### Method overview

1. **Slice simulation**: Divide the atmosphere into layers of constant atmospheric depth, consider radio emission in each layer separately, and then reconstruct the complete signal.
2. **Parameterize amplitude spectra**: Fit the amplitude spectra of each slice using parameterized functions that depend on the number of particles in the slice and the maximum depth \(X_{\text{max}}\) of the shower.
3. **Extract spectral functions**: Obtain spectral functions by fitting the relationship between spectral coefficients and \(X_{\text{max}}\), which describe the spectral parameter values of a specific antenna in a given slice.
4. **Construct templates**: Generate templates using a single microscopic simulation (the origin shower), including normalized amplitude spectra and phase spectra.
5. **Synthesize radio emission**: Synthesize radio emission signals of each slice according to the longitudinal evolution and \(X_{\text{max}}\) of the target shower, using templates and spectral functions.

### Main results

- **Deviation and scatter**: For cases where the difference in \(X_{\text{max}}\) is less than 150 g/cm², the template synthesis method has a deviation of less than 2% and a scatter of less than 6% in the wide frequency band (20 - 500 MHz); in the narrow frequency band (30 - 80 MHz), the deviation is similar, but the scatter drops to 3%.
- **Intrinsic fluctuations**: The observed scatter mainly originates from intrinsic shower fluctuations, which are not explicitly considered in template synthesis.

### Conclusion

The template synthesis method can effectively generate radio emission signals with different longitudinal evolutions, primary particle types and energies while reducing computational cost. The method shows good performance in both wide and narrow frequency bands, and its scatter level is comparable to the intrinsic fluctuations in the Monte Carlo simulation ensemble. Future work will further generalize this method to be applicable to showers of arbitrary geometries.