Composition Of Cosmic Rays Research Articles

Measurement of the depth of maximum of air-shower profiles with energies between 1018.5 and 1020 eV using the surface detector of the Pierre Auger Observatory and deep learning

We report an investigation of the mass composition of cosmic rays with energies from 3 to 100 EeV (1 EeV=1018 eV) using the distributions of the depth of shower maximum Xmax. The analysis relies on ∼50,000 events recorded by the surface detector of the Pierre Auger Observatory and a deep-learning-based reconstruction algorithm. Above energies of 5 EeV, the dataset offers a 10-fold increase in statistics with respect to fluorescence measurements at the Observatory. After cross-calibration using the fluorescence detector, this enables the first measurement of the evolution of the mean and the standard deviation of the Xmax distributions up to 100 EeV. Our findings are threefold: (i) The evolution of the mean logarithmic mass toward a heavier composition with increasing energy can be confirmed and is extended to 100 EeV. (ii) The evolution of the fluctuations of Xmax toward a heavier and purer composition with increasing energy can be confirmed with high statistics. We report a rather heavy composition and small fluctuations in Xmax at the highest energies. (iii) We find indications for a characteristic structure beyond a constant change in the mean logarithmic mass, featuring three breaks that are observed in proximity to the ankle, instep, and suppression features in the energy spectrum. Published by the American Physical Society 2025

Inference of the Mass Composition of Cosmic Rays with Energies from 1018.5 to 1020 eV Using the Pierre Auger Observatory and Deep Learning

We present measurements of the atmospheric depth of the shower maximum Xmax, inferred for the first time on an event-by-event level using the surface detector of the Pierre Auger Observatory. Using deep learning, we were able to extend measurements of the Xmax distributions up to energies of 100 EeV (1020 eV), not yet revealed by current measurements, providing new insights into the mass composition of cosmic rays at extreme energies. Gaining a 10-fold increase in statistics compared to the fluorescence detector data, we find evidence that the rate of change of the average Xmax with the logarithm of energy features three breaks at 6.5±0.6(stat)±1(syst) EeV, 11±2(stat)±1(syst) EeV, and 31±5(stat)±3(syst) EeV, in the vicinity to the three prominent features (ankle, instep, suppression) of the cosmic-ray flux. The energy evolution of the mean and standard deviation of the measured Xmax distributions indicates that the mass composition becomes increasingly heavier and purer, thus being incompatible with a large fraction of light nuclei between 50 and 100 EeV. Published by the American Physical Society 2025

Diffusive shock acceleration of dust grains at supernova remnants

Diffusive shock acceleration (DSA) is a prominent mechanism for energizing charged particles up to very large rigidities at astrophysical collisionless shocks. In addition to ions and electrons, it has been proposed that interstellar dust grains could also be accelerated through diffusive shock acceleration; for instance, at supernova remnants (SNRs). Considering interstellar dust grains of various sizes and compositions, we investigate the possibility of grain acceleration at young SNR shocks (throughout the free expansion and Sedov-Taylor phases) and the maximum energies reached by the accelerated grains. We investigate the potential implications for the abundance of refractory species relative to volatile elements in the cosmic-ray composition. We rely on semi-analytical descriptions of particle acceleration at strong shocks, and on self-similar solutions for the dynamics of SNR shock waves. For simplicity, type Ia thermonuclear SNRs expanding in a uniform interstellar medium are considered. We find that the acceleration of dust grains at relativistic speed is possible, up to a Lorentz factor of $ $ and a kinetic energy of $E_ k nuc 10^2$ GeV/nuc for the smaller grains of size $a $ cm. We find that the subsequent sputtering of grains can produce nuclei with a sufficient rigidity to be injected in the process of diffusive shock acceleration. Such a scenario can help naturally account for the overabundance of refractory elements in the Galactic cosmic-ray composition, provided that a fraction, $ $, of dust grains swept up by an SNR are energized through DSA.

Search for photons above 1018 eV by simultaneously measuring the atmospheric depth and the muon content of air showers at the Pierre Auger Observatory

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above 1017 eV. It measures extensive air showers generated by ultrahigh energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the shower maximum (Xmax) and a steeper lateral distribution function, along with a lower number of muons with respect to the bulk of hadron-induced cascades. In this work, a new analysis technique in the energy interval between 1 and 30 EeV (1 EeV=1018 eV) has been developed by combining the fluorescence detector-based measurement of Xmax with the specific features of the surface detector signal through a parameter related to the air shower muon content, derived from the universality of the air shower development. No evidence of a statistically significant signal due to photon primaries was found using data collected in about 12 years of operation. Thus, upper bounds to the integral photon flux have been set using a detailed calculation of the detector exposure, in combination with a data-driven background estimation. The derived 95% confidence level upper limits are 0.0403, 0.01113, 0.0035, 0.0023, and 0.0021 km−2 sr−1 yr−1 above 1, 2, 3, 5, and 10 EeV, respectively, leading to the most stringent upper limits on the photon flux in the EeV range. Compared with past results, the upper limits were improved by about 40% for the lowest energy threshold and by a factor 3 above 3 EeV, where no candidates were found and the expected background is negligible. The presented limits can be used to probe the assumptions on chemical composition of ultrahigh energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles. Published by the American Physical Society 2024

Mass composition of ultrahigh energy cosmic rays from distribution of their arrival directions with the Telescope Array

Mass composition of ultrahigh energy cosmic rays from distribution of their arrival directions with the Telescope Array

Isotropy of Cosmic Rays beyond 10^{20} eV Favors Their Heavy Mass Composition.

We report an estimation of the injected mass composition of ultrahigh energy cosmic rays (UHECRs) at energies higher than 10EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extragalactic magnetic fields (EGMFs), the results are consistent with a relatively heavy injected composition at E∼10 EeV that becomes lighter up to E∼100 EeV, while the composition at E>100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extragalactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant.

A Multiple Scattering-Based Technique for Isotopic Identification in Cosmic Rays

Analyzing the isotopic composition of cosmic rays (CRs) provides valuable insights into the galactic environment and helps refine existing propagation models. A particular interest is devoted to secondary-to-primary ratios of light isotopic components of CRs, the measurement of which can provide complementary information with respect to secondary-to-primary ratios like B/C. Given the complexity of the concurrent measurement of velocity and momentum required to differentiate isotopes of the same Z, a task typically accomplished using magnetic spectrometers, existing measurements of these ratios only effectively characterize the low-energy region (below 1 GeV/nucl). This study introduces a novel technique for isotopic distinction in CRs at high energies up to 100 GeV/nucl based on multiple scattering, which, combined with the proposed measurement of velocity, represent an interesting alternative to magnetic spectrometers. The performance of this technique was assessed through a dedicated simulation using the GEANT4 package, with specific emphasis on Z = 1 isotopes.

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.

Improving composition of ultrahigh energy cosmic rays with ground detector data

We show that the maximum shower depth (Xmax) distributions of ultrahigh energy cosmic rays (UHECRs), as measured by fluorescence telescopes, can be augmented by building a mapping to observables collected by surface detectors. The resulting statistical improvement of such an augmented dataset depends in a universal way on the strength of the correlation exhibited by the mapping. Building upon the publicly available data on “golden hybrid” events from the Pierre Auger Observatory, we project possible improvements in the inferred composition of UHECRs for a range of possible mappings with varying correlation strengths. Published by the American Physical Society 2024

Simple and less hadronic-model-dependent method for cosmic-ray composition determination

Simple and less hadronic-model-dependent method for cosmic-ray composition determination

Reconstruction of air shower muon lateral distribution functions using integrator and binary modes of underground muon detectors

The investigation of cosmic rays holds significant importance in the realm of particle physics, enabling us to expand our understanding beyond atomic confines. However, the origin and characteristics of ultra-high-energy cosmic rays remain elusive, making them a crucial topic of exploration in the field of astroparticle physics. Currently, our examination of these cosmic rays relies on studying the extensive air showers (EAS) generated as they interact with atmospheric nuclei during their passage through Earth’s atmosphere. Accurate comprehension of cosmic ray composition is vital in determining their source. Notably, the muon content of EAS and the atmospheric depth of the shower maximum serve as the most significant indicators of primary mass composition. In this study, we present two novel methods for reconstructing particle densities based on muon counts obtained from underground muon detectors (UMDs) at varying distances to the shower axis. Our methods were analyzed using Monte Carlo air shower simulations. To demonstrate these techniques, we utilized the muon content measurements from the UMD of the Pierre Auger cosmic ray Observatory, an array of detectors dedicated to measuring extensive air showers. Our newly developed reconstruction methods, employed with two distinct UMD data acquisition modes, showcased minimal bias and standard deviation. Furthermore, we conducted a comparative analysis of our approaches against previously established methodologies documented in existing literature.

Peculiarities of the Isotope Composition of Solar Cosmic Rays

Peculiarities of the Isotope Composition of Solar Cosmic Rays

Results from recent analysis of KASCADE-Grande data

KASCADE and its extension array of KASCADE-Grande were devoted to measure individual air showers of cosmic rays in the primary energy range of 100 TeV to 1 EeV. The experiment has substantially contributed to investigate the energy spectrum and mass composition of cosmic rays in the transition region from galactic to extragalactic origin of cosmic rays as well as to quantify the characteristics of hadronic interaction models in the air shower development through validity tests using the multi-detector information from KASCADE-Grande. Although the data accumulation was completed in 2013, data analysis is still continuing. Recently, we investigated the reliability of the new hadronic interactions models of the SIBYLL version 2.3d only with the energy spectra from the KASCADE-Grande data. The evolution of the muon content of high energy air showers in the atmosphere is studied as well, using EPOS-LHC, SIBYLL 2.3, QGSJET-II-04 and SIBYLL 2.3c. In this talk, recent results from KASCADE-Grande and the update of the KASCADE Cosmic Ray Data Centre (KCDC) will be discussed.

Probing hadronic interaction models with the hybrid data of the Pierre Auger Observatory

Presently large systematic uncertainties remain in the description of hadronic interactions at ultra-high energies and a fully consistent description of air-shower experimental data is yet to be reached. The amount of data collected by the Pierre Auger Observatory using simultaneously the fluorescence and surface detectors in the energy range 10^{18.5}18.5 -10^{19.0}19.0 eV has provided opportunity to perform a multi-parameter test of model predictions. We apply a global method to simultaneously fit the mass composition of cosmic rays and adjustments to the simulated depth of shower maximum, and hadronic signal at ground level. The best description of the hybrid data is obtained for a deeper scale of simulated depth of shower maximum than predicted by hadronic interaction models tuned to the LHC data. Consequently, the deficit of the simulated hadronic signal at ground level, dominated by muons, is alleviated with respect to the unmodified hadronic interaction models. Because of the size of the adjustments to simulated depth of shower maximum and hadronic signal and the large number of events in the sample, the statistical significance of these assumed adjustments is large, greater than 5\sigma_{stat}σstat, even for the combination of the systematic experimental shifts within 1\sigma_{sys}σsys that are the most favorable for the models.

The composition of cosmic rays according to the data on EAS cores

The main purpose of this work is to find the causes of the break of the cosmic ray spectrum at an energy of 3 PeV, which is called the knee. The solution of the problem is associated with the determination of the nuclear composition of cosmic rays in the knee area. The conclusions of this work are based on the analysis of the characteristics of the EAS cores obtained using X-ray emulsion chambers. According to these data, a number of anomalous effects are observed in the knee region, such as scaling violation in the spectra of secondary hadrons, an excess of muons in EAS with gamma families and others. At the same energies equivalent to 1-100 PeV the laboratory system colliders show scaling behavior. So analysis of the data on the EAS cores suggests that the knee in their spectrum is formed by a component of cosmic rays of a non-nuclear nature, possibly consisting of stable (quasi-stable) particles of hypothetical strange quark matter, which named strangelets. This is the only model of the knee compatible with the magnetic rigidity of the nuclear spectra break R=100 TV. In fact, stranglets are stable heavy quasi-nuclei with a positive electric charge of Z=30-1000, so the mechanism of their acceleration coincides with the nuclear one. The break of the cosmic ray spectrum can be associated with a significantly larger mass of strangelets compared to nuclei.

Highlights of the results from the GRAPES-3 experiment

The GRAPES-3 experiment is a unique, extensive air shower experiment consisting of 400 scintillator detectors spread over 25000 m^22 and a 560 m^22 muon telescope. The experiment located at Ooty, India, has been collecting data for the past two decades. The unique capabilities of GRAPES-3 have allowed the study of cosmic rays over energies from a few TeV to tens of PeV and beyond. The measurement of the directional flux of muons (E_\muμ≥1 GeV) by the large muon telescope permits an excellent gamma-hadron separation, which then becomes a powerful tool in the study of multi-TeV gamma-ray sources and the composition of primary cosmic rays. However, the high precision measurements also enable studies of transient atmospheric and interplanetary phenomena such as those produced by thunderstorms and geomagnetic storms. This paper presents some exciting new and recent results, including updates on various ongoing analyses.

Machine learning approach to the background reduction in singly charged cosmic-ray isotope measurements with AMS-02

Studying the isotopic composition of single-charge cosmic rays (CRs) provides essential data to investigate the CR propagation processes in our Galaxy. While current measurements are rare above 4 GeV/nucleon, the Alpha Magnetic Spectrometer (AMS-02) is able to measure the isotopic fluxes up to 10 GeV/n by combining the momentum measured by the silicon tracker with the precise measurements of the velocity provided by its Ring Imaging Cherenkov Detector (RICH). The correct measurement of the particles’ velocity is essential for identifying isotopes through their mass. This is particularly challenging for single-charge particles due to the low number of photons they produce in their Cherenkov rings, which makes the reconstruction easily disrupted by noise. Hence, identifying the sources and cleaning the sample from the background is essential for ensuring the quality of the rings. In this paper, we propose a novel approach to track the events whose mass is misidentified due to interactions inside the AMS-02 detector. Based on the actual location of these interactions, we propose a novel strategy to mitigate the background effectively and with high efficiency, which includes using cut-based selection criteria and a multivariate estimator based on the signals detected by the RICH.

Tachyonic Neutrinos: From the Cosmic Rays to Extragalactic Supernovae

The possibility of neutrinos moving at faster-than-light speeds can be modeled using terms in the Lagrangian that violate Lorentz symmetry, but the question of whether m2<0 and v>c tachyonic neutrinos exist is an empirical question. It remains unresolved despite evidence from cosmic ray and other data that the electron neutrino has an effective mν2≈−0.25eV2, which would require that one or more mass states is also tachyonic. In 2013, the 3+3 model of the neutrino masses, which includes one tachyonic mass state, was proposed based on supernova SN 1987A neutrino data. Here, we update evidence for tachyonic electron neutrinos and the 3+3 model and discuss one test which could prove conclusive. The update of earlier evidence includes many new elements, including new data which make the earlier cosmic ray evidence more robust, new results on cosmic ray composition, the ankle of the spectrum, leptonic cosmic ray data, and the statistical significance of finding the three large neutrino masses stipulated in the 3+3 model. Barring a galactic supernova, which occur only around twice a century, a decisive test of the 3+3 model could involve observing an extragalactic supernova neutrino burst, that is, a cluster of neutrinos in a specific time window well beyond what chance would predict. Even though existing searches for such bursts have yielded only upper limits on the extragalactic supernova frequency within a certain distance, it is shown that the choice of a one-day window for possible neutrino clusters in time might be far more sensitive. A search using a one-day time window could be conducted using existing data, and if a signal is found it would confirm the 3+3 model. Of course, the absence of any day-long neutrino burst would not disprove the model, since it could mean only that the nearest supernova during the period when detectors were active was simply too far to be detected. Finally, apart from testing the 3+3 model, an alternative type of search is suggested using existing hadronic cosmic ray data (from the IceCube Collaboration) that might verify the tachyonic neutrino hypothesis.

Interpretation of fluxes of cosmic rays’ nuclei and electrons in the nonclassical diffusion model

We discuss new scenario of the formation of the observed spectra of electrons/positrons and nuclei of ultrahigh-energy cosmic rays in the framework of the nonclassical (superdiffusion) model of particles propagation in a highly inhomogeneous interstellar and intergalactic media. It is shown that the proposed scenario provides a consistent description of the experimental data of precision satellite measurements, ground-based hybrid EAS arrays and Cherenkov telescopes on the spectra of leptons and nuclei, as well as the mass composition of cosmic rays in the range of super-high and ultrahigh energies.

Effects of large-scale magnetic fields on the observed composition of ultrahigh-energy cosmic rays

Ultra high-energy (UHE) cosmic rays (CRs) from distant sources interact with intergalactic radiation fields, leading to their spallation and attenuation. They are also deflected in intergalactic magnetic fields (IGMFs), particularly those associated with Mpc-scale structures. These deflections extend the propagation times of CR particles, forming a magnetic horizon for each CR species. The cumulative cooling and interactions of a CR ensemble also modifies their spectral shape and composition observed on Earth. We construct a transport formulation to calculate the observed UHE CR spectral composition for 4 classes of source population. The effects on CR propagation brought about by IGMFs are modeled as scattering processes during transport, by centers associated with cosmic filaments. Our calculations demonstrate that IGMFs can have a marked effect on observed UHE CRs, and that source population models are degenerate with IGMF properties. Interpretation of observations, including the endorsement or rejection of any particular source classes, thus needs careful consideration of the structural properties and evolution of IGMFs. Future observations providing tighter constraints on IGMF properties will significantly improve confidence in assessing UHE CR sources and their intrinsic CR production properties.