Composition Of Ultra-high Energy Cosmic Rays Research Articles

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

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.

Nuclear and electromagnetic cascades induced by ultra-high-energy cosmic rays in radio galaxies: implications for Centaurus A

ABSTRACT Very high energy (VHE) γ-rays ($\gtrsim\!\! 0.1\rm ~TeV$) and neutrinos are crucial for identifying accelerators of ultra-high-energy cosmic rays (UHECRs), but this is challenging especially for UHECR nuclei. In this work, we develop a numerical code to solve the transport equation for UHECRs and their secondaries, where both nuclear and electromagnetic cascades are taken into account self-consistently, considering steady UHECR accelerators such as radio galaxies. In particular, we focus on Centaurus A, which has been proposed as one of the most promising UHECR sources in the local Universe. Motivated by observations of extended VHE γ-ray emission from its kiloparsec-scale jet by the High Energy Stereoscopic System (H.E.S.S.), we study interactions between UHECRs accelerated in the large-scale jet and various target photon fields including blazar-like beamed core emission, and present a quantitative study on VHE γ-ray signatures of UHECR nuclei, including the photodisintegration and Bethe–Heitler pair production processes. We show that VHE γ-rays from UHECR nuclei could be detected by the ground-based γ-ray telescopes given that the dominant composition of UHECRs consists of intermediate-mass (such as oxygen) nuclei.

Evaluations of uncertainties in simulations of propagation of ultrahigh-energy cosmic-ray nuclei derived from microscopic nuclear models

Photodisintegration is a main energy loss process for ultrahigh-energy cosmic-ray (UHECR) nuclei in intergalactic space. Therefore, it is crucial to understand systematic uncertainty in photodisintegration when simulating the propagation of UHECR nuclei. In this work, we calculated the cross sections using the random phase approximation (RPA) of density functional theory (DFT), a microscopic nuclear model. We calculated the E1 strength of 29 nuclei using three different density functionals. We obtained the cross sections of photonuclear reactions, including photodisintegration, with the E1 strength. Then, we implemented the cross sections in the cosmic-ray propagation code CRPropa. We found that assuming certain astrophysical parameter values, the difference between UHECR energy spectrum predictions using the RPA calculation and the default photodisintegration model in CRPropa can be more than the statistical uncertainty of the spectrum. We also found that the differences between the RPA calculations and CRPropa default in certain astrophysical parameters obtained by a combined fit of UHECR energy spectrum and composition data assuming a phenomenological model of UHECR sources can be more than the uncertainty of the data.

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.

A Method to Estimate the Primary Mass of Ultrahigh-energy Cosmic Rays from Measurements of the Depth of the Shower Maximum and the Muon Content

The origin and mass composition of ultrahigh-energy cosmic rays is one of the big open questions of modern physics. Estimating the mass of ultrahigh-energy cosmic rays always involves the comparison of data with simulations using hadronic interaction models. The state-of-the-art hadronic interaction models, however, fail to agree with each other on the expected depths of the shower maxima, and on the average number of muons produced in air showers initiated by the same primary particles. Thus the interpretation of data in terms of the primary masses of cosmic rays differs depending on which hadronic interaction model is considered as a reference. Good agreement, however, can be found in the implications from hadronic interaction models in the prediction of the change of the average depth of the shower maximum and in the relative change in the number of produced muons as a function of the nuclear mass of the primary cosmic ray. Considering these, a model-independent estimation of the primary mass of ultrahigh-energy cosmic rays can be given in relative terms. In this work we propose a method to estimate the primary mass of ultrahigh-energy cosmic rays from the depth of the shower maximum and the number of muons produced in the shower, utilizing the similarities of modern hadronic interaction models. We discuss a geometric approach that combines the depth of the shower maximum and the number of muons produced in a shower to estimate the mass of the primary particle.

Constraints on the Hosts of UHECR Accelerators

Interactions of ultra-high-energy cosmic rays (UHECRs) in the surroundings of their accelerators can naturally explain the observed spectrum and composition of UHECRs, including the abundance of protons below the ankle. Here we show that astrophysical properties of the UHECR source environment such as the temperature, size, and magnetic field can be constrained by UHECR and neutrino data. Applying this to candidate sources with a simple structure shows that starburst galaxies are consistent with these constraints, but galaxy clusters are in tension with them. For multicomponent systems like active galactic nuclei and gamma-ray bursts, the results are indicative, but a customized analysis is needed for definitive conclusions.

Ultra-High-Energy Cosmic Rays from a Population of Non-identical Sources

Astrophysical candidates for the sources of ultra-high-energy cosmic rays (UHECRs) exhibit a large diversity in terms of their properties relevant for the acceleration of charged particles, such as luminosity, Lorentz factor, size and magnetic field. Yet, fits of the observed UHECR spectrum and composition often assume identical sources. Here we investigate a population of sources with a power-law distribution of maximum energies. We show that the allowed source-to-source variance of the maximum energy must be small to describe the UHECR data. Even in the most extreme scenario, with a very sharp cutoff of individual source spectra and negative redshift evolution of the accelerators, the maximum energies of 90% of sources must be identical within a factor of three – in contrast to the variance expected for astrophysical sources. However, the overall population variance can be large when maximum rigidities are distributed as a broken power law, with a steep decline above the break and with hard source spectra. In this scenario, most of the observed UHECR flux is produced by sources near the break.

UHECR Signatures and Sources

We discuss recent results on the clustering, composition and distribution of Ultra-High Energy Cosmic Rays (UHECR) in the sky; from the energy of several tens of EeV in the dipole anisotropy, up to the highest energy of a few narrow clusters, those of Hot Spots. Following the early UHECR composition records deviations from proton we noted that the UHECR events above 40 EeV can be made not just by any light or heavy nuclei, but mainly by the lightest ones as He,D, Li,Be. The remarkable Virgo absence and the few localized nearby extragalactic sources, such as CenA, NGC 253 and M82, are naturally understood: lightest UHECR nuclei cannot reach us from the Virgo distance of twenty Mpc, due to their nuclei fragility above a few Mpc distances. Their deflection and smearing in wide hot spots is better tuned to the lighter nuclei than to the preferred proton or heavy nuclei candidate courier. We note that these lightest nuclei still suffer of a partial photodistruction even from such close sources. Therefore, their distruption in fragments, within few tens EeV multiplet chain of events, have been expected and later on observed by Auger collaboration, nearly a decade ago. These multiplet presences, strongly correlate with the same CenA, NGC253 sources. The statistical weight of such correlation is reminded. We conclude that the same role of NGC 253 clustering at lower energies could also feed the Auger dipole anisotropy at lower energy ranges. Such lower energy anisotropy could be fed and integrated by nearest Vela, Crab, LMC and Cas A contributes. In our present UHECR model, based on lightest nuclei in local volumes of a few Mpcs, closest AGN, Star-Burst or very close SNR are superimposing their signals, frozen in different epochs, distances and directions, feeding small and wide anisotropy. Possible tests to confirm, or untangle the current model from alternative ones, are suggested and updated.

Observational Constraints on Cosmic-Ray Escape from Ultrahigh-energy Accelerators

Interactions of ultra-high energy cosmic rays (UHECRs) accelerated in specific astrophysical environments have been shown to shape the energy production rate of nuclei differently from that of the secondary neutrons escaping from the confinement zone. Here, we aim at testing a generic scenario of in-source interactions through phenomenological modeling of the flux and composition of UHECRs. We fit a model in which nucleons and nuclei follow different particle energy distributions to the all-particle energy spectrum and proton spectrum below the ankle energy and distributions of maximum shower depths above this energy, as inferred at the Pierre Auger Observatory. We obtain that the data can be reproduced using a spatial distribution of sources that follows the density of extragalactic matter on both local and large scales, providing hence a realistic set of constraints for the emission mechanisms in cosmic accelerators, for their energetics, and for the abundances of elements at escape from their environments. While the quasi monoelemental increase of the cosmic-ray mass number observed on Earth from ≃2 EeV up to the highest energies calls for nuclei accelerated with a hard spectral index, the inferred flux of protons down to ≃0.6 EeV is shown to require for this population a spectral index significantly softer than that generally obtained up to now. We demonstrate that modeling UHECR data across the ankle substantiates the conjecture of in-source interactions in a robust statistical framework, although pushing the mechanism to the extreme.

What can be learnt from UHECR anisotropies observations

Context.In recent years, evidence for an anisotropic distribution of ultra-high-energy cosmic rays (UHECRs) has been claimed, notably a dipole modulation in right ascension has been reported by the Auger collaboration above the 5σsignificance threshold.Aims.We investigate the implications of the current data regarding large-scale anisotropies, including higher order multipoles, and we examine to what extent they can be used to shed some light on the origin of UHECRs and constrain the astrophysical and/or physical parameters of the source scenarios. We investigate the possibility of observing an associated anisotropy of the UHECR composition and discuss the potential benefit of a good determination of the composition and of the separation of the different nuclear components. We also discuss the interest and relevance of observing the UHECR sky with larger exposure future observatories.Methods.We simulated realistic UHECR sky maps for a wide range of astrophysical scenarios satisfying the current observational constraints, taking into account the energy losses and the photo-dissociation of the UHE protons and nuclei, as well as their deflexions by intervening magnetic fields. We investigated scenarios in which the UHECR source distribution follows that of the galaxies in the Universe (with possible biases), varying the UHECR source composition and spectrum, as well as the source density and the magnetic field models. For each of them, we simulated 300 realizations of independent datasets corresponding to various assumptions for the statistics and sky coverage, and we applied similar analyses as those used by the Auger collaboration for the search of large-scale anisotropies.Results.We find the following. First, reproducing the amplitude of the first-order (dipole) anisotropy observed in the Auger data, as well as its evolution as a function of energy, is relatively easy within our general assumptions. Second, this general agreement can be obtained with different sets of assumptions on the astrophysical and physical parameters, and thus it cannot be used, at the present stage, to derive strong constraints on the UHECR source scenarios or draw model-independent constraints on the various parameters individually. Third, the actual direction of the dipole modulation reconstructed from the Auger data, in the energy bin where the signal is most significant, appears highly unnatural in essentially all scenarios investigated, and this calls for their main assumptions to be reconsidered, either regarding the source distribution itself or the assumed magnetic field configuration, especially in the Galaxy. Fourth, the energy evolution of the reconstructed dipole direction contains potentially important information, which may become constraining for specific source models when larger statistics is collected. Fifth, for such high-statistics datasets, most of our investigated scenarios predict a significant quadrupolar modulation, especially if the light component of UHECRs can be extracted from the all-particle dataset. Sixth, except for protons, the energy range in which the GZK horizon strongly reduces is a key target for anisotropy searches for each given nuclear species. Seventh, although a difference in the average composition of the UHECRs in regions having a different count rate is naturally expected in our models, it is unlikely that the composition anisotropy recently reported by Auger can be explained by this effect, unless the reported amplitude is a strong positive statistical fluctuation of an intrinsically weaker signal.

Explaining the UHECR spectrum, composition and large-scale anisotropies with radio galaxies

Radio galaxies are promising candidates as the sources of ultrahigh energy cosmic rays (UHECRs). In this work, we examine if the stringent constraints imposed by the dipole and quadropole anisotropies as well as the UHECR spectrum and composition allow that radio galaxies are the dominant extragalactic cosmic ray sources. In order to calculate the UHECR flux emitted by individual radio galaxies, we constrain their properties using information from the radio-CR correlation and a dynamical evolution model. In addition to the UHECR flux from individual, local sources, we include the diffuse flux emitted by the bulk of non-local radio galaxies based on their radio luminosity distribution. Analyzing the source parameters within a range around their expected properties, we finally determine the configurations of local sources describing well the UHECR spectrum, composition and large-scale anisotropies. We obtain a good description of all data even in the case that we include only a small number of local sources. In particular, we find that scenarios where few sources like Fornax A and Virgo A dominate the flux above the ankle, while low-luminosity radio galaxies contribute an isotropic background dominating below the ankle, provide a good fit to the data.

Mass Composition of UHECRs from Xmax Distributions Recorded by the Pierre Auger and Telescope Array Observatories

In this paper we infer the mass composition of the ultra high energy cosmic rays (UHECRs) from measurements of Xmax distributions recorded at the Pierre Auger (2014) and Telescope Array (TA) (2016) Observatories, by fitting them with all possible combinations of Monte Carlo (MC) templates from a large set of primary species (p, He, C, N, O, Ne, Si and Fe), as predicted by EPOS-LHC, QGSJETII-04 and Sibyll 2.1 hadronic interaction models. We use the individual fractions of nuclei reconstructed from one experiment in each energy interval to build equivalent MC Xmax distributions, which we compare with the experimental Xmax distributions of the other experiment, applying different statistical tests of compatibility. The results obtained from both experiments confirm that the mass composition of the UHECRs is dominated (≳70%) by protons and He nuclei in the energy range investigated lgE(eV) = [17.8–19.3] (Auger) and lgE(eV) = [18.2–19.0] (TA). The indirect comparisons between the Xmax distributions recorded by the two experiments show that the degree of compatibility of the two datasets is good, even excellent in some high energy intervals, especially above the ankle (lgE(eV)∼18.7). However, our study reveals that, at low energies, further effort in data analysis is required in order to harmonize the results of the two experiments.

Spectra and composition of ultrahigh-energy cosmic rays and the measurement of the proton-air cross section

The shape of the longitudinal development of the showers generated in the atmosphere by very high energy cosmic ray particles encodes information about the mass composition of the flux, and about the properties of hadronic interactions that control the shower development. Studies of the energy dependence of the average and width of the depth of maximum distribution of showers with $E \gtrsim 10^{17.3}$ eV measured by the Pierre Auger Observatory, suggest, on the basis of a comparison with current models, that the composition of the cosmic ray flux undergoes a very important evolution, first becoming lighter and then rapidly heavier. These conclusions, if confirmed, would have profound and very surprising implications for our understanding of the high energy astrophysical sources. Studies of the shape of the depth of maximum distribution in the same energy range have been used by Auger and by the Telescope Array Collaboration to measure the interaction length of protons in air, a quantity that allows to estimate the $pp$ cross sections for values of $\sqrt{s}$ well above the LHC range. In this paper we argue that it is desirable to combine the studies of the cosmic ray composition with those aimed at the measurement of the $p$--air cross section. The latter allow to obtain estimates for the fraction of protons in the flux that can be of great help in decoding the composition and its energy dependence. Studies that consider multiple parameters to characterize the depth of maximum distributions also offer the possibility to perform more sensitive tests of the validity of the models used to describe high energy showers.

Modeling the spectrum and composition of ultrahigh-energy cosmic rays with two populations of extragalactic sources

We fit the ultrahigh-energy cosmic-ray (UHECR, Egtrsim 0.1 EeV) spectrum and composition data from the Pierre Auger Observatory at energies Egtrsim 5cdot 10^{18} eV, i.e., beyond the ankle using two populations of astrophysical sources. One population, accelerating dominantly protons (^1H), extends up to the highest observed energies with maximum energy close to the GZK cutoff and injection spectral index near the Fermi acceleration model; while another population accelerates light-to-heavy nuclei (^4He, ^{14}N, ^{28}Si, ^{56}Fe) with a relatively low rigidity cutoff and hard injection spectrum. A significant improvement in the combined fit is noted as we go from a one-population to two-population model. For the latter, we constrain the maximum allowed proton fraction at the highest-energy bin within 3.5sigma statistical significance. In the single-population model, low-luminosity gamma-ray bursts turn out to match the best-fit evolution parameter. In the two-population model, the active galactic nuclei is consistent with the best-fit redshift evolution parameter of the pure proton-emitting sources, while the tidal disruption events could be responsible for emitting heavier nuclei. We also compute expected cosmogenic neutrino flux in such a hybrid source population scenario and discuss possibilities to detect these neutrinos by upcoming detectors to shed light on the sources of UHECRs.

Identifying nearby sources of ultra-high-energy cosmic rays with deep learning

We present a method to analyse arrival directions of ultra-high-energy cosmic rays (UHECRs) using a classifier defined by a deep convolutional neural network trained on a HEALPix grid. To illustrate a high effectiveness of the method, we employ it to estimate prospects of detecting a large-scale anisotropy of UHECRs induced by a nearby source with an (orbital) detector having a uniform exposure of the celestial sphere and compare the results with our earlier calculations based on the angular power spectrum. A minimal model for extragalactic cosmic rays and neutrinos by Kachelrieß, Kalashev, Ostapchenko and Semikoz (2017) is assumed for definiteness and nearby active galactic nuclei Centaurus A, M82, NGC 253, M87 and Fornax A are considered as possible sources of UHECRs. We demonstrate that the proposed method drastically improves sensitivity of an experiment by decreasing the minimal required amount of detected UHECRs or the minimal detectable fraction of from-source events several times compared to the approach based on the angular power spectrum. We also test robustness of the neural networks against different models of the large-scale Galactic magnetic fields and variations of the mass composition of UHECRs, and consider situations when there are two nearby sources or the dominating source is not known a priori. In all cases, the neural networks demonstrate good performance unless the test models strongly deviate from those used for training. The method can be readily applied to the analysis of data of the Telescope Array, the Pierre Auger Observatory and other cosmic ray experiments.

Particle dynamics and GZK limit in relativity with a preferred frame

The theory termed ‘special relativity (SR) with a preferred frame’ incorporates the preferred frame into special relativity while retaining the fundamental space-time symmetry which, in the standard SR, manifests itself as Lorentz invariance. In this paper, the theory is extended to particle dynamics based on the concept of the ‘modified spacetime symmetry’ which allows to develop the theory along the lines of the standard relativity theory. The modified relativistic dynamics is applied to describing the interactions of the Ultra High Energy Cosmic Rays (UHECR) with universal diffuse background radiation and, in particular, to calculating the Greisen–Zatsepin–Kuzmin (GZK) threshold for the energy suppression due to pion photoproduction by UHECR protons. The result, that the GZK cutoff energy depends on the distance (cosmological redshift z) to the source of particles, may contribute to interpretation of the data on the mass composition of UHECR.

Systematic parameter space study for the UHECR origin from GRBs in models with multiple internal shocks

ABSTRACT We scrutinize the paradigm that conventional long-duration gamma-ray bursts (GRBs) are the dominant source of the ultrahigh energy cosmic rays (UHECRs) within the internal shock scenario by describing UHECR spectrum and composition and by studying the predicted (source and cosmogenic) neutrino fluxes. Since it has been demonstrated that the stacking searches for astrophysical GRB neutrinos strongly constrain the parameter space in single-zone models, we focus on the dynamics of multiple collisions for which different messengers are expected to come from different regions of the same object. We propose a model that can describe both stochastic and deterministic engines, which we study in a systematic way. We find that GRBs can indeed describe the UHECRs for a wide range of different model assumptions with comparable quality albeit with the previously known problematic energy requirements; the heavy mass fraction at injection is found to be larger than 70 per cent ($95 {{\ \rm per\ cent}}$ CL). We demonstrate that the post-dicted (from UHECR data) neutrino fluxes from sources and UHECR propagation are indeed below the current sensitivities but will be reached by the next generation of experiments. We finally critically review the required source energetics with the specific examples found in this study.

Evidence for a Supergalactic Structure of Magnetic Deflection Multiplets of Ultra-high-energy Cosmic Rays

Abstract Evidence for a large-scale supergalactic cosmic-ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high-energy cosmic-ray (UHECR) energies above 1019 eV using 7 years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy–position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, because the supergalactic plane is a fit to the average matter density of the local large-scale structure, UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy–angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of the TA (the Hotspot and Coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data posttrial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2σ. The 10 years of data posttrial significance is 4.1σ. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

Galactic magnetic field bias on inferences from UHECR data

A consequence of Liouville's theorem indicates that the recently observed large scale anisotropy in the arrival direction of Ultra-High-Energy Cosmic Rays (UHECRs) cannot be produced from an isotropic extragalactic flux, thus some anisotropy already needs to be present outside our Galaxy. But in this case, the observed spectrum and composition of UHECRs differ from the one outside of the Milky Way, due to the suppression or the amplification of the UHECR flux from certain directions by the Galactic magnetic field. In this work, we investigate this effect for the case of a dipole and a quadrupole anisotropy, respectively, for the widely-used JF12 and PT11 magnetic field models. We investigate bounds on the maximal amplitude of the observed anisotropy and the maximal charge number of UHECRs. Furthermore, the flux modification is discussed in the light of the recent Auger results about the dipolar anisotropy and the mass composition. We find that this effect yields a modification of the observed flux of up to ∼ 25% for the investigated magnetic field model and the observed dipole, in particular for a heavy chemical composition of UHECRs as suggested by Auger data interpreted assuming the `Sibyll2.3c' or the `EPOS-LHC' model of hadronic interactions in air showers.