Sibyll$^{\bigstar}$

Felix Riehn,Anatoli Fedynitch,Ralph Engel

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

2024-04-09

Abstract:In the last decade, an increasing number of datasets have revealed a consistent discrepancy between the number of muons measured in ultra-high-energy extensive air showers (EAS) and the numbers predicted by simulations. This gap persists despite incorporating Large Hadron Collider (LHC) data into the tuning of current hadronic interaction models, leading to the phenomenon often termed the ''muon puzzle''. To gain a deeper understanding of the potential origins of this muon puzzle, we have developed Sibyll$^{\bigstar}$, a series of phenomenologically modified versions of Sibyll 2.3d. In these models, we have increased muon production by altering $\rho^0$, baryon-antibaryon pair, or kaon production in hadronic multiparticle production processes. These variants remain within bounds from provided by accelerator measurements, including those from the LHC and fixed-target experiments, notably NA49 and NA61, showing a level of consistency comparable to Sibyll 2.3d. Our findings show that these modifications can increase the muon count in EAS by up to 35%, while minimally affecting the depth of shower maximum ($X_{\rm max}$) and other shower variables. Additionally, we assess the impact of these modifications on various observables, including inclusive muon and neutrino fluxes and the multiplicities of muon bundles in deep underground and water/ice Cherenkov detectors. We aim for at least one of these model variants to offer a more accurate representation of EAS data at the highest energies, thereby enhancing the quality of Monte Carlo predictions used in training neural networks. This improvement is crucial for achieving more reliable data analyses and interpretations.

High Energy Physics - Phenomenology,High Energy Astrophysical Phenomena

What problem does this paper attempt to address?

This paper discusses the "muon puzzle," i.e. the discrepancy between the observed number of muons in ultra-high-energy air showers (EAS) and the predicted number from simulations. The study modifies the Sibyll 2.3d model to create a series of variants called Sibyll*, which increase muon production in EAS while maintaining consistency with accelerator experiment data. These modifications aim to more accurately simulate EAS data, reduce systematic uncertainties in data analysis and interpretation, and improve Monte Carlo predictions used for training neural networks.