Ultra High Energy Particles Propagation and the Transition from Galactic to Extra-Galactic Cosmic Rays

We discuss the basic features of the propagation of Ultra High Energy Cosmic Rays in astrophysical backgrounds, comparing two alternative computation schemes to compute the expected fluxes. We also discuss the issue of the transition among galactic and extra-galactic cosmic rays using theoretical results on fluxes to compare different models.

💡 Research Summary

The paper provides a comprehensive review of ultra‑high‑energy cosmic‑ray (UHECR) propagation and the transition from Galactic to extragalactic origins. It begins by outlining the fundamental energy‑loss processes that shape the observable spectrum: electron‑positron pair production on the cosmic microwave background (CMB) dominates at ≲10¹⁸ eV, photodisintegration of nuclei on the infrared‑optical background becomes important in the 10¹⁸–10¹⁹ eV range, and photopion production (the Greisen‑Zatsepin‑Kuzmin, or GZK, effect) sharply suppresses the flux above ≈5×10¹⁹ eV.

Two computational frameworks are compared. The first is a continuous‑energy‑loss (CEL) kinetic‑equation approach, which treats the average loss rate as a deterministic function of energy and redshift. This method yields rapid estimates of the propagated spectrum but neglects stochastic fluctuations in interaction depth and the detailed evolution of nuclear composition. The second framework employs Monte Carlo simulations (e.g., SimProp, CRPropa) that follow each particle’s trajectory, randomly sampling interaction probabilities, secondary production, and cosmological expansion. Monte Carlo results reproduce the CEL spectrum at low energies but diverge near the GZK cutoff, where they predict a more pronounced steepening and a rapid shift toward heavier nuclei due to photodisintegration.

The authors then examine the Galactic‑to‑extragalactic transition using two representative models. In the “dip” scenario, the extragalactic component is proton‑dominated; the ankle at ~10¹⁸.⁵ eV is interpreted as a spectral dip caused by pair production on the CMB. This model yields a relatively sharp transition and matches the observed average depth of shower maximum (Xmax) and its fluctuations if the composition remains light. In the “mixed‑composition” scenario, sources emit a spectrum of nuclei (He, CNO, Fe) with a rigidity‑dependent cutoff. The transition is smoother, extending from ~10¹⁸ to 10¹⁹ eV, and the Xmax evolution reflects an increasing average mass with energy. By confronting both models with current data from the Pierre Auger Observatory and the Telescope Array, the paper finds that the mixed‑composition framework provides a marginally better fit to the combined spectrum‑and‑composition measurements, although systematic uncertainties (energy scale, hadronic interaction models) prevent a decisive discrimination.

Sensitivity studies highlight the impact of source evolution (parameterized as (1+z)ᵐ) and the maximum acceleration energy (Emax). Strong positive evolution enhances the high‑energy tail, sharpening the GZK feature, while higher Emax pushes the transition to larger energies and allows heavier nuclei to survive to Earth.

In the concluding section, the authors argue that future progress hinges on high‑statistics, composition‑sensitive observations. Upgraded facilities such as AugerPrime, TA×4, and planned space‑based detectors (e.g., POEMMA) will deliver the necessary precision to disentangle stochastic propagation effects from source characteristics. Integrating detailed Monte Carlo propagation with multi‑messenger data (neutrinos, gamma rays) is identified as the most promising pathway to resolve the long‑standing question of where Galactic cosmic rays end and extragalactic UHECRs begin.