Ultra-high Energy Cosmic Rays Research Articles

Quantum Gravity Phenomenology Induced in the Propagation of UHECR, a Kinematical Solution in Finsler and Generalized Finsler Spacetime

It is well-known that the universe is opaque to the propagation of Ultra-High-Energy Cosmic Rays (UHECRs) since these particles dissipate energy during their propagation interacting with the background fields present in the universe, mainly with the Cosmic Microwave Background (CMB) in the so-called GZK cut-off phenomenon. Some experimental evidence seems to hint at the possibility of a dilation of the GZK predicted opacity sphere. It is well-known that kinematical perturbations caused by supposed quantum gravity (QG) effects can modify the foreseen GZK opacity horizon. The introduction of Lorentz Invariance Violation can indeed reduce, and in some cases making negligible, the CMB-UHECRs interaction probability. In this work, we explore the effects induced by modified kinematics in the UHECR lightest component phenomenology from the QG perspective. We explore the possibility of a geometrical description of the massive fermions interaction with the supposed quantum structure of spacetime in order to introduce a Lorentz covariance modification. The kinematics are amended, modifying the dispersion relations of free particles in the context of a covariance-preserving framework. This spacetime description requires a more general geometry than the usual Riemannian one, indicating, for instance, the Finsler construction and the related generalized Finsler spacetime as ideal candidates. Finally we investigate the correlation between the magnitude of Lorentz covariance modification and the attenuation length of the photopion production process related to the GZK cut-off, demonstrating that the predicted opacity horizon can be dilated even in the context of a theory that does not require any privileged reference frame.

Acceleration of ultrahigh-energy cosmic rays in the early afterglows of gamma-ray bursts: Concurrence of jet dynamics and wave-particle interactions

The origin of ultra-high energy cosmic rays $($UHECRs$)$ remains a mystery. It has been suggested that UHECRs can be produced by the stochastic acceleration in relativistic jets of gamma-ray bursts $($GRBs$)$ at the early afterglow phase. Here, we develop a time-dependent model for proton energization by cascading compressible waves in GRB jets considering the concurrent effect of the jet's dynamics and the mutual interactions between turbulent waves and particles. Considering fast mode of magnetosonic wave as the dominant particle scatterer and assuming interstellar medium $($ISM$)$ for the circumburst environment, our numerical results suggest that protons can be accelerated up to $\textrm{10}^{\textrm{19}}\,$eV during the early afterglow. An estimation shows ultra-high energy nuclei can easily survive photodisintegration in the external shocks in most cases, thus allowing the acceleration of $\textrm{10}^{\textrm{20}}\,$eV cosmic-ray nuclei in the proposed frame. The spectral slope can be as hard as $\textrm{d}\textit{N}/\textrm{d}\textit{E} \propto \textit{E}^\textrm{0}$, which is consistent with the requirement for the interpretation of intermediate-mass composition of UHECR as measured by the Pierre Auger Observatory.

Ultra-High-energy Cosmic Rays from beyond the Greisen–Zatsepin–Kuz’min Horizon

Abstract Ultra-high-energy (UHE) cosmic rays (CRs) of energies ∼(1018–1020) eV, accelerated in violent astrophysical environments, interact with cosmic background radiation fields via photo-hadronic processes, leading to strong attenuation. Typically, the Universe would become “opaque” to UHE CRs after several tens of megaparsecs, setting the boundary of the Greisen–Zatsepin–Kuz’min (GZK) horizon. In this work, we investigate the contribution of sources beyond the conventional GZK horizon to the UHE CR flux observed on Earth, when photospallation of the heavy nuclear CRs is taken into account. We demonstrate that this contribution is substantial, despite the strong attenuation of UHE CRs. A significant consequence is the emergence of an isotropic background component in the observed flux of UHE CRs, coexisting with the anisotropic foreground component that is associated with nearby sources. Multi-particle CR horizons, which evolve over redshift, are determined by the CR nuclear composition. Thus, they are dependent on the source populations and source evolutionary histories.

Extensive Air Shower Reconstruction using the timing information from the RF-system of the Astroneu array

The Astroneu cosmic ray telescope is a distributed hybrid array consisting of both scintillator counters and RF antenna detectors used for the detection of extensive air showers (EAS). The array is deployed at the Hellenic Open University campus, on the outskirts of the urban area of Patras in Greece. In the present development phase, the Astroneu telescope includes two stations consisting of 3 scintillation detectors modules (SDM) and one RF antenna while a third station includes 3 particle detectors and 4 RF antennas (3SDM-4RF). In each station, the RF-detectors are operating receiving a common trigger upon a 3-fold coincidence between the particle detectors of the station. In this study we present recent results from the 3SDM-4RF autonomous station related to the estimation of the direction of the incoming cosmic air shower using only the timing information from the 4 RF detectors. The directions of the reconstructed showers using the RF timing are in agreement with the corresponding results using the SDMs timing as well as with the simulation predictions. This verifies that the RF signal emitted from EAS originating form Ultra High Energy Cosmic Rays (UHECR), can be detected even in areas with strong electromagnetic background.

Ultrahigh energy cosmic ray deflection by the intergalactic magnetic field

The origin and composition of ultrahigh energy cosmic rays (UHECRs) remains a mystery. The common lore is that UHECRs are deflected from their primary directions by the Galactic and extragalactic magnetic fields. Here we describe an extragalactic contribution to the deflection of UHECRs that does not depend on the strength and orientation of the initial seed field. Using the IllustrisTNG simulations, we show that outflow-driven magnetic bubbles created by feedback processes during galaxy formation deflect approximately half of all ${10}^{20}\text{ }\text{ }\mathrm{eV}$ protons by 1\ifmmode^\circ\else\textdegree\fi{} or more, and up to 20\ifmmode^\circ\else\textdegree\fi{}--30\ifmmode^\circ\else\textdegree\fi{}. This implies that the deflection in the intergalactic medium must be taken into account in order to identify the sources of UHECRs.

Cosmic rays, neutrinos, and GeV-TeV gamma rays from starburst galaxy NGC 4945

The detection of high-energy astrophysical neutrinos and ultra-high-energy cosmic rays (UHECRs) provides a new way to explore sources of cosmic rays. One of the highest energy neutrino events detected by IceCube, tagged as IC35, is close to the UHECR anisotropy region detected by Pierre Auger Observatory. The nearby starburst (SB) galaxy, NGC 4945, is close to this anisotropic region and inside the mean angular error of the IC35 event. Considering the hypernovae contribution located in the SB region of NGC 4945, which can accelerate protons up to $\ensuremath{\sim}{10}^{17}\text{ }\text{ }\mathrm{eV}$ and inject them into the interstellar medium, we investigate the origin of this event around this starburst galaxy. We show that the interaction of these protons with the SB region's gas density could explain Fermi-LAT gamma-ray and radio observations if the magnetic field's strength in the SB region is the order of $\ensuremath{\sim}\mathrm{mG}$. Our estimated PeV neutrino events, in ten years, for this source is approximately 0.01 ($4\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}$) if a proton spectral index of 2.4 (2.7) is considered, which would demonstrate that IC35 is not produced in the central region of this SB galaxy. Additionally, we consider the superwind region of NGC 4945 and show that protons can hardly be accelerated in it up to UHECRs.

A Simulation Study of Ultra-relativistic Jets. II. Structures and Dynamics of FR-II Jets

Abstract We study the structures of ultra-relativistic jets injected into the intracluster medium and the associated flow dynamics, such as shocks, velocity shear, and turbulence, through three-dimensional relativistic hydrodynamic (RHD) simulations. To that end, we have developed a high-order accurate RHD code, equipped with a weighted essentially non-oscillatory scheme and a realistic equation of state. Using the code, we explore a set of jet models with the parameters relevant to FR-II radio galaxies. We confirm that the overall jet morphology is primarily determined by the jet power, and the jet-to-background density and pressure ratios play secondary roles. Jets with higher powers propagate faster, resulting in more elongated structures, while those with lower powers produce more extended cocoons. Shear interfaces in the jet are dynamically unstable, and hence, chaotic structures with shocks and turbulence develop. We find that the fraction of the jet-injected energy dissipated through shocks and turbulence is greater in less powerful jets, although the actual amount of the dissipated energy is larger in more powerful jets. In lower power jets, the backflow is dominant in the energy dissipation owing to the broad cocoon filled with shocks and turbulence. In higher power jets, by contrast, both the backflow and jet-spine flow are important for the energy dissipation. Our results imply that different mechanisms, such as diffusive shock acceleration, shear acceleration, and stochastic turbulent acceleration, may be involved in the production of ultra-high energy cosmic rays in FR-II radio galaxies.

Bi-Spectral Infrared Algorithm for Cloud Coverage over Oceans by the JEM-EUSO Mission Program.

The need to monitor specific areas for different applications requires high spatial and temporal resolution. This need has led to the proliferation of ad hoc systems on board nanosatellites, drones, etc. These systems require low cost, low power consumption, and low weight. The work we present follows this trend. Specifically, this article evaluates a method to determine the cloud map from the images provided by a simple bi-spectral infrared camera within the framework of JEM-EUSO (The Joint Experiment Missions-Extrem Universe Space Observatory). This program involves different experiments whose aim is determining properties of Ultra-High Energy Cosmic Ray (UHECR) via the detection of atmospheric fluorescence light. Since some of those projects use UV instruments on board space platforms, they require knowledge of the cloudiness state in the FoV of the instrument. For that reason, some systems will include an infrared (IR) camera. This study presents a test to generate a binary cloudiness mask (CM) over the ocean, employing bi-spectral IR data. The database is created from Moderate-Resolution Imaging Spectroradiometer (MODIS) data (bands 31 and 32). The CM is based on a split-window algorithm. It uses an estimation of the brightness temperature calculated from a statistical study of an IR images database along with an ancillary sea surface temperature. This statistical procedure to obtain the estimate of the brightness temperature is one of the novel contributions of this work. The difference between the measured and estimation of the brightness temperature determines whether a pixel is cover or clear. That classification requires defining several thresholds which depend on the scenarios. The procedure for determining those thresholds is also novel. Then, the results of the algorithm are compared with the MODIS CM. The agreement is above 90%. The performance of the proposed CM is similar to that of other studies. The validation also shows that cloud edges concentrate the vast majority of discrepancies with the MODIS CM. The relatively high accuracy of the algorithm is a relevant result for the JEM-EUSO program. Further work will combine the proposed algorithm with complementary studies in the framework of JEM-EUSO to reinforce the CM above the cloud edges.

Magnetic diffusion and interaction effects on ultrahigh energy cosmic rays: Protons and nuclei

The flux of ultrahigh energy cosmic rays reaching the Earth is affected by the interactions with the cosmic radiation backgrounds as well as with the magnetic fields that are present along their trajectories. We combine the SimProp cosmic ray propagation code with a routine that allows to account for the average effects of a turbulent magnetic field on the direction of propagation of the particles. We compute in this way the modification of the spectrum which is due to the magnetic horizon effect, both for primary nuclei as well as for the secondary nuclei resulting from the photo-disintegration of the primary ones. We also provide analytic parameterizations of the attenuation effects, as a function of the magnetic field parameters and of the density of cosmic ray sources, which make it possible to obtain the expected spectra in the presence of the magnetic fields from the spectra that would be obtained in the absence of magnetic fields. The discrete nature of the distribution of sources with finite density also affects the spectrum of cosmic rays at the highest energies where the flux is suppressed due to the interactions with the radiation backgrounds, and parameterizations of these effects are obtained.

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.

Energetics of ultrahigh-energy cosmic-ray nuclei

Energetics of the ultrahigh-energy cosmic rays (UHECRs) generated in the universe is crucial for pinning down their candidate sources. Using the recent Auger data on UHECR spectra, we calculate the UHECR energy generation rate density for different species of nuclei at the injection, considering intermediate and heavy nuclei as well as protons, through scanning over source parameters on the spectral index, maximum energy and redshift evolution. We find the resulting UHECR energy generation rate density to be $\approx(0.2-2)\times$10$^{44}$~erg~Mpc$^{-3}$~yr$^{-1}$ at $10^{19.5}$~eV with a nontrivial dependence on the spectral index. Nuclei other than protons and irons favor hard spectral indices at the injection, and the required value of energy budget is smaller for intermediate nuclei. Given significant uncertainties in hadronic interaction models and astrophysical models for the Galactic-extragalactic transition, our results can be regarded as conservative. The composition data on $X_{\rm max}$ give additional constraints, but the results are consistent within the model uncertainties.

Extreme Universe Space Observatory on a Super Pressure Balloon 1 calibration: from the laboratory to the desert

The Extreme Universe Space Observatory on a Super Pressure Balloon 1 (EUSO-SPB1) instrument was launched out of Wanaka, New Zealand, by NASA in April, 2017 as a mission of opportunity. The detector was developed as part of the Joint Experimental Missions for the Extreme Universe Space Observatory (JEM-EUSO) program toward a space-based ultra-high energy cosmic ray (UHECR) telescope with the main objective to make the first observation of UHECRs via the fluorescence technique from suborbital space. The EUSO-SPB1 instrument is a refractive telescope consisting of two 1m$^2$ Fresnel lenses with a high-speed UV camera at the focal plane. The camera has 2304 individual pixels capable of single photoelectron counting with a time resolution of 2.5$\mu$s. A detailed performance study including calibration was done on ground. We separately evaluated the properties of the Photo Detector Module (PDM) and the optical system in the laboratory. An end-to-end test of the instrument was performed during a field campaign in the West Desert in Utah, USA at the Telescope Array (TA) site in September 2016. The campaign lasted for 8 nights. In this article we present the results of the preflight laboratory and field tests. Based on the tests performed in the field, it was determined that EUSO-SPB1 has a field of view of 11.1$^\circ$ and an absolute photo-detection efficiency of 10%. We also measured the light flux necessary to obtain a 50% trigger efficiency using laser beams. These measurements were crucial for us to perform an accurate post flight event rate calculation to validate our cosmic ray search. Laser beams were also used to estimated the reconstruction angular resolution. Finally, we performed a flat field measurement in flight configuration at the launch site prior to the launch providing a uniformity of the focal surface better than 6%.

Surface detectors of the TAx4 experiment

Telescope Array (TA) is the largest ultrahigh energy cosmic-ray (UHECR) observatory in the Northern Hemisphere. It explores the origin of UHECRs by measuring their energy spectrum, arrival-direction distribution, and mass composition using a surface detector (SD) array covering approximately 700 km2 and fluorescence detector (FD) stations. TA has found evidence for a cluster of cosmic rays with energies greater than 57 EeV. In order to confirm this evidence with more data, it is necessary to increase the data collection rate. We have begun building an expansion of TA that we call TAx4. In this paper, we explain the motivation, design, technical features, and expected performance of the TAx4 SD. We also present TAx4’s current status and examples of the data that have already been collected.

Detecting ultra-high-energy cosmic ray anisotropies through harmonic cross-correlations

We propose an observable for ultra-high-energy cosmic ray (UHECR) physics: the harmonic-space cross-correlation power spectrum between the arrival directions of UHECRs and the large-scale cosmic structure mapped by galaxies. This cross-correlation has not yet been considered in the literature, and it permits a direct theoretical modelling of the main astrophysical components. We describe the expected form of the cross-correlation and show how, if the distribution of UHECR sources traces the large-scale cosmic structure, it could be easier to detect with current data than the UHECR auto-correlation. Moreover, the cross-correlation is more sensitive to UHECR anisotropies on smaller angular scales, it is more robust to systematic uncertainties, and it could be used to determine the redshift distribution of UHECR sources, making it a valuable tool for determining their origins and properties.

Application of Neural Networks to Classification of Data of the TUS Orbital Telescope

We employ neural networks for classification of data of the TUS fluorescence telescope, the world’s first orbital detector of ultra-high energy cosmic rays. We focus on two particular types of signals in the TUS data: track-like flashes produced by cosmic ray hits of the photodetector and flashes that originated from distant lightnings. We demonstrate that even simple neural networks combined with certain conventional methods of data analysis can be highly effective in tasks of classification of data of fluorescence telescopes.

Constraining light dark matter upscattered by ultrahigh-energy cosmic rays

Light halo dark matter (DM) particles upscattered by high-energy cosmic rays (CRs) can be energetic, and become detectable by conventional direct detection experiments. The current constraints derived from space-based direct CR measurements can reach O(10−31) cm2 for a constant DM-nucleon scattering cross section. We show that if the CR energy spectrum follows a power law of type ∼E−3, the derived constraints on the scattering cross section will be highly insensitive to DM particle mass. This suggests that ultrahigh-energy CRs (UHECRs) indirectly measured by ground-based detectors can be used to place constraints on ultralight DM particles, as E−3 is a very good approximation of the UHECR energy spectrum up to energy ∼1020 eV. Using the recent UHECR flux data, we show that the current constraints derived from space-based CR measurements can in principle be extended to ultralight DM particles far below eV scale.

“Super GZK” Particles in a Classic Kramers’ Diffusion-over-a-barrier Model. I. The Case of Protons

Abstract In interactions of ultra-high-energy cosmic ray (UHECR) protons with cosmic microwave background photons, we focus in this work on photopion production reactions and the effects of the measured, broad, energy-loss distributions in these reactions on the evolution of the protons’ density functions in energy space. We rely on a Fokker–Planck transport equation in energy space whose transport coefficients are calculated using laboratory measurements. We also derive a Fokker–Planck potential that accounts for both systematic (drift) and stochastic (dispersive) energy losses due to photopion production reactions. Our results show that dispersive energy losses have significant effects on estimating the protons’ horizon distance and their energy spectrum, as well as to elucidate a broadness in the GZK cutoff. We use the derived Fokker–Planck potential to assign a characteristic probability for a proton to clear the potential barrier as a function of energy. This estimate of probability can be used to assist observations in distinguishing between protons and heavy ions as charged particles. Our model is able to account for the so-called super GZK particles as a classic diffusion-over-a-barrier manifestation of the transport of UHECR protons in energy space in quantifying the extent and broadness of the GZK cutoff.

Small shower CORSIKA simulations

Extensive Air Showers (EAS) induced by cosmic ray particles of very low energies, owing to the significantly steep cosmic ray energy spectrum, dominate the secondary particle flux measured by single detectors and small shower arrays. Such arrays connected in extended networks can be used to determine potentially interesting spatial correlations between showers, which may shed new light on the nature of ultra high-energy cosmic rays. The quantitative interpretation of showers recorded by small local arrays requires a methodology that differs from that used by ordinary large EAS arrays operating in the "knee" region and above. We present "small EAS generator," a semi-analytical method for integrating cosmic ray spectra over energies of interest and summing over the mass spectra of primary nuclei in arbitrary detector configurations. Furthermore, we provide results on the EAS electron and muon fluxes and particle density spectra.

Hunting super-heavy dark matter with ultra-high energy photons

At any epoch, particle physics must be open to completely unexpected discoveries, and that is reason enough to extend the reach of searches for ultra-high energy (UHE) photons. The observation of a population of photons with energies E≳100EeV would for example imply the existence of either a completely new physical phenomena, or particle acceleration mechanisms heretofore never seen or imagined. But as we outline in this Letter of Interest, there are also good arguments for super-heavy dark matter (SHDM) in a parameter range such that it could be discovered via its decays to, in particular, UHE photons. Only ultra-high energy cosmic ray observatories have capabilities to detect UHE photons. We first investigate how current and future observations can probe and constrain SHDM models in important directions, and then outline some of the scenarios that motivate such searches. We also discuss connections between constraints on SHDM and on the parameter values of cosmological models.

Probing UHECR production in Centaurus A using secondary neutrinos and gamma-rays

In this paper the production of neutrinos and photons by ultra high energy cosmic rays (UHECR) interacting with the extragalactic background radiation is studied. Centaurus A is assumed as the prime source of UHECR and the possibility to identify this source by detecting the secondary neutrinos and photons produced in the propagation of UHECR is investigated. Fifteen astrophysical models regarding three extragalactic magnetic fields (EGMF) and five composition abundances are simulated. The flux and arrival direction of neutrinos and photons are investigated. It is shown that the detection of a signal from Cen A with statistical significance is achievable by current observatories in a few years and by proposed experiments in the near future. The dependence of the results on the models is also presented.