Hadronic Interactions Research Articles

On the model uncertainties for the predicted muon content of extensive air showers

Motivated by the excess of the muon content of cosmic ray induced extensive air showers (EAS), relative to EAS modeling, observed by the Pierre Auger Observatory, and by the tension between Auger data and air shower simulations on the maximal muon production depth Xmaxμ, we investigate the possibility to modify the corresponding EAS simulation results, within the Standard Model of particle physics. We start by specifying the kinematic range for secondary hadron production, which is of relevance for such predictions. We further investigate the impact on the predicted EAS muon number and on Xmaxμ of various modifications of the treatment of hadronic interactions, in the framework of the QGSJET-III model, in particular the model calibration to accelerator data, the amount of the “glue” in the pion, and the energy dependence of the pion exchange process. None of the considered modifications of the model allowed us to enhance the EAS muon content by more than 10%. On the other hand, for the maximal muon production depth, some of the studied modifications of particle production give rise up to ∼10 g/cm2 larger Xmaxμ values, which increases the difference with Auger observations.

The population of Galactic supernova remnants in the TeV range

Supernova remnants (SNRs) are likely to be significant sources of cosmic rays up to the knee of the local cosmic-ray (CR) spectrum. They produce gamma rays in the very-high-energy (VHE) ($E>0.1$ TeV) range mainly via two mechanisms: hadronic interactions of accelerated protons with the interstellar medium and leptonic interactions of accelerated electrons with soft photons. Observations with current instruments have lead to the detection of about a dozen SNRs emitting VHE gamma rays and future instruments should significantly increase this number. Yet, the details of particle acceleration at SNRs and of the mechanisms producing VHE gamma-rays at SNRs remain poorly understood. We aim to study the population of SNRs detected in the TeV range and its properties and confront it to simulated samples in order to address fundamental questions concerning particle acceleration at SNR shocks. Such questions concern the spectrum of accelerated particles, the efficiency of particle acceleration, and the gamma-ray emission being dominated by hadronic or leptonic interactions. By means of Monte Carlo methods, we simulated the population of SNRs in the gamma--ray domain and confronted our simulations to the catalogue of sources from the High Energy Stereoscopic System (H.E.S.S.) Galactic Plane Survey (HGPS). We systematically explored the parameter space defined in our model, including for example, the slope of accelerated particles alpha , the electron-to-proton ratio $K_ ep $, and the efficiency of particle acceleration $ $. In particular, we found possible sets of parameters for which $ 90$<!PCT!> of Monte Carlo realisations are found to be in agreement with the HGPS. These parameters are typically found at $ 4.2 ep $, and $0.03 0.1 $ . We are able to strongly argue against some regions of the parameter space describing the population of Galactic SNRs in the TeV range, such as $ 4.35$, or ep Our model is so far able to explain the SNR population of the HGPS. Our approach, when confronted with the results of future systematic surveys, such as the Cherenkov Telescope Array Observatory, will help remove degeneracy from the solutions and to better understand particle acceleration at SNR shocks in the Galaxy.

Observations of GRB 221009A by the MAVEN-SEP Instrument at Mars

We report on the detection of the Brightest Of All Time (the BOAT) gamma-ray burst (GRB) 221009A by the MAVEN Solar Energetic Particle (SEP) silicon (Si) detectors. The SEP instrument on board the MAVEN spacecraft at Mars is designed to measure charged particle fluxes and energies. In this work, we show that the SEP Si detectors have detected populations of secondary charged particles and low-energy primary photons from the GRB. Our analysis relies on a series of simplified Geant4 Monte Carlo simulations. We show that electromagnetic and hadronic interactions between the GRB photons and the spacecraft’s passive material and Si detectors produced the detected secondary particles.

QGSJET-III model of high energy hadronic interactions. II. Particle production and extensive air shower characteristics

The hadronization procedure of the QGSJET-III Monte Carlo (MC) generator of high energy hadronic interactions is discussed. Selected results of the model, regarding production spectra of secondary particles, are presented in comparison to experimental data and to the corresponding predictions of the QGSJET-II-04 MC generator. The model is applied to calculations of basic characteristics of extensive air showers initiated by cosmic ray interactions in the atmosphere and the results are compared to predictions of other MC generators of cosmic ray interactions. Published by the American Physical Society 2024

Testing hadronic-model predictions of depth of maximum of air-shower profiles and ground-particle signals using hybrid data of the Pierre Auger Observatory

We test the predictions of hadronic interaction models regarding the depth of maximum of air-shower profiles, Xmax, and ground-particle signals in water-Cherenkov detectors at 1000 m from the shower core, S(1000), using the data from the fluorescence and surface detectors of the Pierre Auger Observatory. The test consists of fitting the measured two-dimensional (S(1000), Xmax) distributions using templates for simulated air showers produced with hadronic interaction models pos-, et--04, 2.3d and leaving the scales of predicted Xmax and the signals from hadronic component at ground as free-fit parameters. The method relies on the assumption that the mass composition remains the same at all zenith angles, while the longitudinal shower development and attenuation of ground signal depend on the mass composition in a correlated way. The analysis was applied to 2239 events detected by both the fluorescence and surface detectors of the Pierre Auger Observatory with energies between 1018.5 eV to 1019.0 eV and zenith angles below 60°. We found, that within the assumptions of the method, the best description of the data is achieved if the predictions of the hadronic interaction models are shifted to deeper Xmax values and larger hadronic signals at all zenith angles. Given the magnitude of the shifts and the data sample size, the statistical significance of the improvement of data description using the modifications considered in the paper is larger than 5σ even for any linear combination of experimental systematic uncertainties. Published by the American Physical Society 2024

Cosmic-ray diffusion in two local filamentary clouds

Context. Hadronic interactions between cosmic rays (CRs) and interstellar gas have been probed in γ rays across the Galaxy. A fairly uniform CR distribution is observed up to a few hundred parsecs from the Sun, except in the Eridu cloud, which shows an unexplained 30–50% deficit in GeV to TeV CR flux. Aims. To explore the origin of this deficit, we studied the Reticulum cloud, which shares notable traits with Eridu: a comparable distance in the low-density region of the Local Valley and a filamentary structure of atomic hydrogen extending along a bundle of ordered magnetic-field lines that are steeply inclined to the Galactic plane. Methods. We measured the γ-ray emissivity per gas nucleon in the Reticulum cloud in the 0.16–63 GeV energy band using 14 years of Fermi-LAT data. We also derived interstellar properties that are important for CR propagation in both the Eridu and Reticulum clouds, at the same parsec scale. Results. The γ-ray emissivity in the Reticulum cloud is fully consistent with the average spectrum measured in the solar neighbourhood, but this emissivity, and therefore the CR flux, is 1.57 ± 0.09 times larger than in Eridu across the whole energy band. The difference cannot be attributed to uncertainties in gas mass. Nevertheless, we find that the two clouds are similar in many respects: both have magnetic-field strengths of a few micro-Gauss in the plane of the sky; both are in approximate equilibrium between magnetic and thermal pressures; they have similar turbulent velocities and sonic Mach numbers; and both show magnetic-field regularity with a dispersion in orientation lower than 10°–15° over large zones. The gas in Reticulum is colder and denser than in Eridu, but we find similar parallel diffusion coefficients around a few times 1028 cm2 s−1 in both clouds if CRs above 1 GV in rigidity diffuse on resonant, self-excited Alfvén waves that are damped by ion-neutral interactions. Conclusions. The loss of CRs in Eridu remains unexplained, but these two clouds provide important test cases to further study how magnetic turbulence, line tangling, and ion-neutral damping regulate CR diffusion in the dominant gas phase of the interstellar medium.

Energy spectra of elemental groups of cosmic rays with the KASCADE experiment data and machine learning

We report the reconstruction of the mass component spectra of cosmic rays (protons, helium, carbon, silicon and iron) and their mean mass composition, at energies from 1.4 to 100 PeV. The results are derived from the archival data of the extensive air shower experiment KASCADE. We use a novel machine learning technique developed specifically for this reconstruction, and post-LHC hadronic interaction models: QGSJet-II.04, EPOS-LHC and Sibyll 2.3c. We have found an excess of the proton component and a deficit of intermediate and heavy nuclei components compared to the original KASCADE results. The spectra of protons and helium show a knee-like behavior at ∼ 4.4 PeV and ∼ 11 PeV, with significances 5.2σ and 3.9σ, respectively. The spectrum of the iron component has a hint (2.4σ) of a hardening at ∼ 4.5 PeV, which can be interpreted as a counterpart of a hardening in the proton spectrum at 166 TeV, recently reported by the GRAPES-3 experiment. The systematic uncertainties of our analysis were found to be smaller than those of the original KASCADE, as well as those of IceTop and TALE experiments, over the most part of the energy range studied. We also estimated separately the uncertainty related to the difference between the three mentioned hadronic interaction models. We also compute a mean logarithm mass of CR flux as a function of energy. It is in agreement with the results of IceTop, TALE and LHAASO within the uncertainties.

Studies of ultra-high-energy cosmic rays by the Yakutsk air shower array, based on the Cherenkov light detection technique

The paper presents the results of long-term, continuous observations of the Cherenkov light at the Yakutsk array. The Cherenkov light measurements allow one to reconstruct the longitudinal development of the shower, which is associated with both the model of hadron interactions and the mass composition of cosmic rays. The paper also provides a review of the results that are important for the physics of air showers and the astrophysics of cosmic rays, obtained over a long period of time at the Yakutsk array.

Double scattering in deuteron electrodisintegration

We demonstrate that at sufficiently high energies when the eikonal regime is established for hadronic interactions, the double scattering subprocess can be clearly identified and isolated in quasi-elastic deuteron electro-disintegration processes. Comparing theoretical calculations with the recent high precision experimental data we present a “proof of principle” that these processes can be used to study advanced issues related to hadron formation in QCD. In this case, the double scattering represents as a fermi-scale “detector” which probes products of high Q2 scattering from the bound nucleon through their rescattering from the spectator nucleon in the deuteron.

Axionlike-particle contributions to the μ→e conversion

We study the μ→e conversion process in nuclear targets arising in models of axionlike particles (ALPs) with hadronic and charged lepton flavor violating (CLFV) interactions. Contributions to this process generally fall into two categories: spin-independent (SI) and spin-dependent (SD). While the SI contribution can be generated by a dipole operator through purely leptonic ALP interactions, the SD contribution can also be present through ALP-quark interactions at tree-level. It is naively anticipated that the SI contribution would be dominant due to its coherent enhancement. In this paper, we show that is not generically the case; in particular, for naturally sized ALP couplings to quarks of order ∼mq/fa, the SD interaction induced by ALP-π0 mixing turns out to be the leading contribution to μ→e conversion. Intuitively, this stems from the suppressed dipole contribution by the QED one-loop factor which counters the effect of SI coherent enhancement. Our study highlights the importance of μ→e conversion searches in exploring the parameter space of generic ALP models, and demonstrates the competitiveness of these searches in probing the CLFV ALP parameter space in the heavy mass range of ma≳mμ. Published by the American Physical Society 2024

Muon measurements at the Pierre Auger Observatory

The Pierre Auger Observatory is the world’s largest detector for observation of ultra-high-energy cosmic rays (UHECRs) (above the energy of 10^{17}1017 eV). It consists of a Fluorescence Detector (FD) and an array of particle detectors known as the Surface Detector (SD). Observations of extensive air showers by the Observatory can be used to probe hadronic interactions at high energy, in a kinematic and energy region inaccessible to experiments at man-made accelerators and to measure the muon component of the shower. Air showers induced by different primaries have different muon contents. With increasing mass of the primary cosmic ray particle, it is expected that the muon content in the corresponding air showers should also increase. Therefore, the determination of the muon component in the air shower is crucial to infer the mass of the primary particle. This is a key ingredient in the searches conducted to pinpoint the sources of UHECRs. Recent results obtained from the Pierre Auger Observatory and other experiments indicate that all the shower simulations underestimate the number of muons in the showers compared to the data. This is the so-called muon deficit. In this paper we briefly review the muon measurements, and present in more detail recent results on fluctuations in the muon number. These results provide new insights into the origin of the muon deficit in air shower simulations and constrain the models of hadronic interactions at ultrahigh energies. With the current design of the surface detectors it is also difficult to reliably separate the contributions of muons to the SD signal from the contributions of photons, electrons, and positrons. Therefore, we also present a new method to extract the muon component of the signal time traces recorded by each SD station using recurrent neural networks. The combination of such algorithms, with the future data collected by the upgraded Pierre Auger Observatory, will be a major step forward, as we are likely to achieve an unprecedented resolution in mass estimation on an event-by-event basis.

Air showers and hadronic interactions with CORSIKA 8

The CORSIKA 8 project is a collaborative effort aiming to develop a versatile C++ framework for the simulation of extensive air showers, intended to eventually succeed the long-standing FORTRAN version. I present an overview of its current capabilities, focusing on aspects concerning the hadronic and muonic shower components. In particular, I demonstrate the “cascade lineage” feature and its application to quantify the importance of certain phase-space regions in hadronic interactions for muon production. Additionally, I show first results using Pythia 8.3, which as of late is usable as interaction model in cosmic-ray applications and is currently being integrated into CORSIKA 8.

Sub-TeV hadronic interaction model differences and their impact on air showers

In the sub-TeV regime, the most widely used hadronic interaction models disagree significantly in their predictions of particle spectra from cosmic ray induced air showers. We investigate the nature and impact of model uncertainties, focussing on air shower primaries with energies around the transition between high and low energy hadronic interaction models, where the dissimilarities are largest and which constitute the bulk of the interactions in air showers.

Insights into Neutron Star Equation of State by Machine Learning

Due to its powerful capability and high efficiency in big data analysis, machine learning has been applied in various fields. We construct a neural network platform to constrain the behaviors of the equation of state of nuclear matter with respect to the properties of nuclear matter at saturation density and the properties of neutron stars. It is found that the neural network is able to give reasonable predictions of parameter space and provide new hints into the constraints of hadron interactions. As a specific example, we take the relativistic mean field approximation in a widely accepted Walecka-type model to illustrate the feasibility and efficiency of the platform. The results show that the neural network can indeed estimate the parameters of the model at a certain precision such that both the properties of nuclear matter around saturation density and global properties of neutron stars can be saturated. The optimization of the present modularly designed neural network and extension to other effective models is straightforward.

Sibyll[formula omitted]

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★, a series of phenomenologically modified versions of Sibyll 2.3d. In these models, we have increased muon production by altering ρ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 (Xmax) 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.

Examining the influence of hadronic interactions on the directed flow of identified particles in RHIC Beam Energy Scan energies using UrQMD model

The directed flow of identified particles can serve as a sensitive tool for investigating the interactions during initial and final states in heavy ion collisions. This study examines the rapidity distribution ([Formula: see text]), rapidity-odd directed flow ([Formula: see text]) and its slope ([Formula: see text]) for [Formula: see text], [Formula: see text], p, and [Formula: see text] in Au+Au collisions at different collision centralities and beam energies ([Formula: see text], 11.5, 14.5, 19.6, 27, and 39[Formula: see text]GeV) using the UrQMD model. We investigate the impact of late-stage hadronic interactions on charge-dependent [Formula: see text] and its slope by modifying the duration of the hadronic cascade lifetime ([Formula: see text]). The energy dependence of [Formula: see text] for p ([Formula: see text]) exhibits distinct pattern compared to [Formula: see text] and [Formula: see text]. Notably, we observe a change in the sign reversal position of proton [Formula: see text] at different beam energies with varying [Formula: see text] in central and mid-central collisions. Moreover, the difference in [Formula: see text] between positively and negatively charged hadrons ([Formula: see text]) demonstrates a stark centrality dependence for different particle species. The deuteron displays a significant increase in [Formula: see text] with increasing [Formula: see text] compared to p and n. This investigation underscores the importance of considering the temporal evolution and duration of the hadronic phase when interpreting the sign reversal, charge splitting of [Formula: see text] and light nuclei formation at lower RHIC energies.

LeHaMoC: A versatile time-dependent lepto-hadronic modeling code for high-energy astrophysical sources

Context. Recent associations of high-energy neutrinos with active galactic nuclei (AGN) have revived the interest in leptohadronic models of radiation from astrophysical sources. The rapid increase in the amount of acquired multi-messenger data will require fast numerical models that may be applied to large source samples. Aims. We develop a time-dependent leptohadronic code, LeHaMoC, that offers several notable benefits compared to other existing codes, such as versatility and speed. Methods. LeHaMoC solves the Fokker-Planck equations of photons and relativistic particles (i.e. electrons, positrons, protons, and neutrinos) produced in a homogeneous magnetized source that may also be expanding. The code utilizes a fully implicit difference scheme that allows fast computation of steady-state and dynamically evolving physical problems. Results. We first present test cases where we compare the numerical results obtained with LeHaMoC against exact analytical solutions and numerical results computed with ATHEvA, a well-tested code of similar philosophy but a different numerical implementation. We find a good agreement (within 10–30%) with the numerical results obtained with ATHEvA without evidence of systematic differences. We then demonstrate the capabilities of the code through illustrative examples. First, we fit the spectral energy distribution from a jetted AGN in the context of a synchrotron-self Compton model and a proton-synchrotron model using Bayesian inference. Second, we compute the high-energy neutrino signal and the electromagnetic cascade induced by hadronic interactions in the corona of NGC 1068. Conclusions. LeHaMoC is easily customized to model a variety of high-energy astrophysical sources and has the potential to become a widely utilized tool in multi-messenger astrophysics.

QGSJET-III model of high energy hadronic interactions: The formalism

The physics content of the QGSJET-III Monte Carlo generator of high energy hadronic collisions is described. In particular, a phenomenological implementation of higher twist corrections to hard parton scattering processes is discussed in some detail. Additionally addressed is the treatment of the so-called “color fluctuation” effects related to a decomposition of hadron wave functions into a number of Fock states characterized by different spatial sizes and different parton densities. Selected model results regarding the energy dependence of the total, elastic, and diffractive proton-proton cross sections are presented. Published by the American Physical Society 2024

Constraints on the proton fraction of cosmic rays at the highest energies and the consequences for cosmogenic neutrinos and photons

Over the last decade, observations have shown that the mean mass of ultra-high-energy cosmic rays (UHECRs) increases progressively toward the highest energies. However, the precise composition is still unknown and several theoretical studies hint at the existence of a subdominant proton component up to the highest energies. Motivated by the exciting prospect of performing charged-particle astronomy with ultra-high-energy (UHE) protons we quantify the level of UHE-proton flux that is compatible with present multimessenger observations and the associated fluxes of neutral messengers produced in the interactions of the protons. We study this scenario with numerical simulations of two independent populations of extragalactic sources and perform a fit to the combined UHECR energy spectrum and composition observables, constrained by diffuse gamma-ray and neutrino observations. We find that up to of order 10% of the cosmic rays at the highest energies can be UHE protons, although the result depends critically on the selected hadronic interaction model for the air showers. Depending on the maximum proton energy (E max p) and the redshift evolution of sources, the associated flux of cosmogenic neutrinos and UHE gamma rays can significantly exceed the multimessenger signal of the mixed-mass cosmic rays. Moreover, if E max p is above the GZK limit, we predict a large flux of UHE neutrinos above EeV energies that is absent in alternate scenarios for the origin of UHECRs. We present the implications and opportunities afforded by these UHE proton, neutrino and photon fluxes for future multimessenger observations.

ALICE luminosity determination for Pb–Pb collisions at √s NN = 5.02 TeV

Luminosity determination within the ALICE experiment is based on the measurement, in van der Meer scans, of the cross sections for visible processes involving one or more detectors (visible cross sections). In 2015 and 2018, the Large Hadron Collider provided Pb–Pb collisions at a centre-of-mass energy per nucleon pair of √s NN = 5.02 TeV. Two visible cross sections, associated with particle detection in the Zero Degree Calorimeter (ZDC) and in the V0 detector, were measured in a van der Meer scan.This article describes the experimental set-up and the analysis procedure, and presents the measurement results. The analysis involves a comprehensive study of beam-related effects and an improved fitting procedure, compared to previous ALICE studies, for the extraction of the visible cross section. The resulting uncertainty of both the ZDC-based and the V0-based luminosity measurement for the full sample is 2.5%. The inelastic cross section for hadronic interactions in Pb–Pb collisions at √s NN = 5.02 TeV, obtained by efficiency correction of the V0-based visible cross section, was measured to be 7.67 ± 0.25 b, in agreement with predictions using the Glauber model.