Hadronic Interaction Models SIBYLL Research Articles

Density of GeV muons in air showers measured with IceTop

We present a measurement of the density of GeV muons in near-vertical air showers using three years of data recorded by the IceTop array at the South Pole. Depending on the shower size, the muon densities have been measured at lateral distances between 200 m and 1000 m. From these lateral distributions, we derive the muon densities as functions of energy at reference distances of 600 m and 800 m for primary energies between 2.5 PeV and 40 PeV and between 9 PeV and 120 PeV, respectively. The muon densities are determined using, as a baseline, the hadronic interaction model Sibyll 2.1 together with various composition models. The measurements are consistent with the predicted muon densities within these baseline interaction and composition models. The measured muon densities have also been compared to simulations using the post-LHC models EPOS-LHC and QGSJet-II.04. The result of this comparison is that the post-LHC models together with any given composition model yield higher muon densities than observed. This is in contrast to the observations above 1 EeV where all model simulations yield for any mass composition lower muon densities than the measured ones. The post-LHC models in general feature higher muon densities so that the agreement with experimental data at the highest energies is improved but the muon densities are not correct in the energy range between 2.5 PeV and about 100 PeV.

Hadronic interaction model sibyll 2.3d and extensive air showers

We present a new version of the hadron interaction event generator Sibyll. While the core ideas of the model have been preserved, the new version handles the production of baryon pairs and leading particles in a new way. In addition, production of charmed hadrons is included. Updates to the model are informed by high-precision measurements of the total and inelastic cross sections with the forward detectors at the LHC that constrain the extrapolation to ultra-high energy. Minimum-bias measurements of particle spectra and multiplicities support the tuning of fragmentation parameters. This paper demonstrates the impact of these changes on air shower observables such as $X_{\rm max}$ and $N_\mu$, drawing comparisons with other contemporary cosmic ray interaction models.

The hadronic interaction model Sibyll 2.3c and muon production in extensive air-showers

One of the applications of the hadronic interaction model Sibyll is the simulation of extensive air showers of ultra-high energy cosmic rays. In recent years it has become more and more clear that simulations do not agree with measurements when it comes to observables related to muons in air showers. We discuss the processes in Sibyll that are directly related to muon production in extensive air showers and describe their relation to shower observables.

Study of themuon content of high-energy air showers with KASCADE-Grande

In this work, we report measurements on the muon content (Eth > 230 MeV) of extensive air showers (EAS) induced by cosmic rays with primary energy from 10 PeV up to 1 EeV performed with the KASCADE-Grande experiment. The measurements are confronted with SIBYLL 2.3. The results are focused on the dependence of the total muon number and the lateral density distribution of muons in EAS on the zenith angle and the total number of charged particles in the shower. We also present updated results of a detailed study of the attenuation length of shower muons, which reveal a deviation between the measured data and the predictions of the post-LHC hadronic interaction models SIBYLL 2.3, QGSJET-II-04 and EPOS-LHC.

The hadronic interaction model Sibyll – past, present and future

The hadronic interaction model Sibyll – past, present and future

The hadronic interaction model Sibyll – past, present and future

Sibyll is one of the first microscopic interaction models that was specifically developed for interpreting cosmic ray data. It combines non-perturbative concepts of simulating hadronic particle production with predictions derived from perturbative QCD calculations, focusing on forward particle production of relevance in studying cosmic ray interactions. In this contribution we briefly recall the history of Sibyll and then, in this context, describe improvements made in the different versions of the Sibyll model. The discussion focuses on the basic concepts and ideas of these improvements rather than going into detail or giving a comprehensive description of the models. We also discuss shortcomings, conceptual problems, and uncertainties in modeling hadronic interactions and make some suggestions how to address these open questions in the future.

The contribution of charm to high energy atmospheric neutrinos

We discuss the perfomance of the new version of the Sibyll hadronic interaction model Sibyll 2.xx which corrects some of the problems of Sibyll 2.1 and also generates charm hadrons and thus produces very high energy atmospheric neutrinos.

Influence of hadronic interaction models and the cosmic ray spectrum on the high-energy atmospheric muon and neutrino flux

The recent observations of muon charge ratio up to about 10 TeV and of atmospheric neutrinos up to energies of about 400 TeV has triggered a renewed interest into the high-energy interaction models and cosmic ray primary composition. A reviewed calculation of lepton spectra produced in cosmic ray induced extensive air showers is carried out with a primary cosmic ray spectrum that fits the latest direct measurements below the knee. In order to achieve this, we used a full Monte Carlo method to derive the inclusive differential spectra (yields) of muons, muon neutrinos and electron neutrinos at the surface for energies between 80 GeV and hundreds of PeV. Using these results the differential flux and the flavor ratios of leptons were calculated. The air shower simulator CORSIKA 6.990 was used for showering and propagation of the secondary particles through the atmosphere, employing the established high energy hadronic interaction models SIBYLL 2.1, QGSJet-01 and QGSJet-II-03. We show that the performance of the interaction models allows makes it possible to predict the spectra within experimental uncertainties, while SIBYLL generally yields a higher flux at the surface than the QGSJet models. The calculation of the flavor and charge ratios has lead to inconsistent results, mainly influenced by the different representations of the K/π ratio within the models. The influence of the knee of cosmic rays is reflected in the secondary spectra at energies between 100 and 200 TeV. Furthermore, we could quantify systematic uncertainties of atmospheric muon- and neutrino fluxes, associated to the models of the primary cosmic ray spectrum and the interaction models. For most recent parametrizations of the cosmic ray primary spectrum, atmospheric muons can be determined with an uncertainty smaller than +15/-13% of the average flux. Uncertainties of the muon and electron neutrino fluxes can be calculated within an average error of +32/-22% and +25/-19%, respectively. See the published paper [1].

Influence of hadronic interaction models and the cosmic ray spectrum on the high-energy atmospheric muon and neutrino flux

The recent observations of muon charge ratio up to about 10 TeV and of atmospheric neutrinos up to energies of about 400 TeV has triggered a renewed interest into the high-energy interaction models and cosmic ray primary composition. A reviewed calculation of lepton spectra produced in cosmic ray induced extensive air showers is carried out with a primary cosmic ray spectrum that fits the latest direct measurements below the knee. In order to achieve this, we used a full Monte Carlo method to derive the inclusive differential spectra (yields) of muons, muon neutrinos and electron neutrinos at the surface for energies between 80 GeV and hundreds of PeV. Using these results the differential flux and the flavor ratios of leptons were calculated. The air shower simulator CORSIKA 6.990 was used for showering and propagation of the secondary particles through the atmosphere, employing the established high energy hadronic interaction models SIBYLL 2.1, QGSJet-01 and QGSJet-II-03. We show that the performance of the interaction models allows makes it possible to predict the spectra within experimental uncertainties, while SIBYLL generally yields a higher flux at the surface than the QGSJet models. The calculation of the flavor and charge ratios has lead to inconsistent results, mainly influenced by the different representations of the K/π ratio within the models. The influence of the knee of cosmic rays is reflected in the secondary spectra at energies between 100 and 200 TeV. Furthermore, we could quantify systematic uncertainties of atmospheric muon- and neutrino fluxes, associated to the models of the primary cosmic ray spectrum and the interaction models. For most recent parametrizations of the cosmic ray primary spectrum, atmospheric muons can be determined with an uncertainty smaller than +15/-13% of the average flux. Uncertainties of the muon and electron neutrino fluxes can be calculated within an average error of +32/-22% and +25/-19%, respectively. See the published paper [1].

Influence of hadronic interaction models and the cosmic ray spectrum on the high energy atmospheric muon and neutrino flux

The recent observations of muon charge ratio up to about 10 TeV and of atmospheric neutrinos up to energies of about 400 TeV have triggered a renewed interest into the high-energy interaction models and cosmic-ray primary composition. A reviewed calculation of lepton spectra produced in cosmic-ray-induced extensive air showers is carried out with a primary cosmic-ray spectrum that fits the latest direct measurements below the knee. In order to achieve this, we used a full Monte Carlo method to derive the inclusive differential spectra (yields) of muons, muon neutrinos and electron neutrinos at the surface for energies between 80 GeV and hundreds of PeV. Using these results, the differential flux and the flavor ratios of leptons were calculated. The air shower simulator CORSIKA 6.990 was used for showering and propagation of the secondary particles through the atmosphere, employing the established high-energy hadronic interaction models SIBYLL 2.1, QGSJet-01 and QGSJet-II-03. We show that the performance of the interaction models makes it possible to predict the spectra within experimental uncertainties, while SIBYLL generally yields a higher flux at the surface than the QGSJet models. The calculation of the flavor and charge ratios has lead to inconsistent results, mainly influenced by the different representations of the $K/\ensuremath{\pi}$ ratio within the models. The influence of the knee of cosmic rays is reflected in the secondary spectra at energies between 100 and 200 TeV. Furthermore, we could quantify systematic uncertainties of atmospheric muon and neutrino fluxes, associated to the models of the primary cosmic-ray spectrum and the interaction models. For most recent parametrizations of the cosmic-ray primary spectrum, atmospheric muons can be determined with an uncertainty smaller than $\genfrac{}{}{0}{}{+15}{\ensuremath{-}13}%$ of the average flux. Uncertainties of the muon and electron neutrino fluxes can be calculated within an average error of $\genfrac{}{}{0}{}{+32}{\ensuremath{-}22}%$ and $\genfrac{}{}{0}{}{+25}{\ensuremath{-}19}%$, respectively.