A Hierarchical Shock Model of Ultra-High-Energy Cosmic Rays

We propose that a hierarchical shock model$\unicode{x2014}$including supernova remnant shocks, galactic wind termination shocks, and accretion shocks around cosmic filaments and galaxy clusters$\unicode{x2014}$can naturally explain the cosmic ray spectrum from ~1 GeV up to ~200 EeV. While this framework applies to the entire cosmic ray spectrum, in this work, we focus on its implications for ultra-high-energy cosmic rays (UHECRs). We perform a hydrodynamic cosmological simulation to investigate the power processed at shocks around clusters and filaments. The downstream flux from nearby shocks around the local filament accounts for the softer, lower-energy extragalactic component around the ankle, and the upstream escaping flux from nearby clusters accounts for the transition to a hard spectral component at the highest energies. This interpretation is in agreement with UHECR observations. We suggest that a combination of early-Universe galactic outflows, cosmic ray streaming instabilities, and a small-scale turbulent dynamo can increase magnetic fields enough to attain the required rigidities. Our simulation suggests that the available volume-averaged power density of accretion shocks exceeds the required UHECR luminosity density by three orders of magnitude. We show that microgauss magnetic fields at these shocks could explain both the origin of UHECRs and potentially contribute to the diffuse radio synchrotron background below 10 GHz. The shock-accelerated electrons produce a hard radio background without overproducing diffuse inverse Compton emission. These results motivate further observational tests with upcoming facilities to help distinguish accretion shocks from other UHECR sources.

💡 Research Summary

The paper proposes a unified “hierarchical shock” framework to explain the origin of ultra‑high‑energy cosmic rays (UHECRs) across the entire cosmic‑ray spectrum, from ∼1 GeV up to ∼200 EeV. The hierarchy consists of three classes of non‑relativistic, collisionless shocks: (i) supernova‑remnant (SNR) shocks inside galaxies, (ii) termination shocks of galactic winds, and (iii) large‑scale accretion shocks that form around galaxy clusters and the filamentary cosmic web. In the low‑energy regime (∼1 GeV–10 EeV) particles are first accelerated by SNRs, then re‑accelerated at wind termination shocks, and finally injected into the filamentary shocks that dominate the extragalactic component below the “ankle”. In the highest‑energy regime (∼100 EeV–200 EeV) the quasi‑spherical accretion shocks of massive clusters provide the necessary shock speed (u₁≈1000 km s⁻¹, with rare patches up to 5000 km s⁻¹) and long acceleration times (∼Gyr) to push heavy nuclei to rigidities of order 10 EV, yielding energies up to a few × 10² EeV.

To quantify the energetics, the authors run a radiation‑hydrodynamic AMR simulation with Enzo, including a Haardt‑Madau UV background and non‑equilibrium chemistry for nine species. They identify shocks with Mach number M>5 and pre‑shock overdensity δ<10³, thereby selecting strong, low‑density shocks that are efficient particle accelerators. The simulation volume‑averaged power density processed by these large‑scale structure shocks is P_LSSS≈10⁴⁰ erg s⁻¹ Mpc⁻³, three orders of magnitude larger than the required UHECR luminosity density L_UHECR≈3×10³⁶ erg s⁻¹ Mpc⁻³. For individual halos, the shock radius is taken as r_sh≈2 r_vir, with typical kinetic power ∼10⁴⁴ erg s⁻¹ for a Virgo‑mass cluster (M₁₄≈1). Filamentary shocks contribute roughly 10 % of this power and have typical upstream speeds of ∼400 km s⁻¹.

A central challenge is the generation of μG‑level magnetic fields in the low‑density outskirts where the accretion shocks reside. The authors argue that a combination of early‑Universe galactic outflows (seeding seed fields), cosmic‑ray‑driven streaming instabilities (amplifying them locally), and a small‑scale turbulent dynamo driven by pressure gradients can raise the field to the required strength over several Gyr. Such fields enable diffusive shock acceleration (DSA) to reach the observed maximum rigidities, especially for heavy nuclei (Z≈26), consistent with Auger composition measurements that indicate a trend toward heavier elements at the highest energies.

The paper also connects the shock model to the diffuse radio synchrotron background observed between 22 MHz and 10 GHz. Shock‑accelerated electrons with a power‑law index p≈2.2 produce synchrotron emission with spectral index α≈0.6, matching the measured isotropic radio background. Because the electrons radiate in μG fields, the associated inverse‑Compton X‑ray and γ‑ray emission remains below current background limits, satisfying an important observational constraint.

Observational predictions are outlined: (1) The transition from the softer extragalactic component (below the ankle) to the hard UHECR component should be dominated by nearby filament shocks (downstream flux) and nearby cluster shocks (upstream escaping flux), leading to anisotropies aligned with the local filament and Virgo‑like clusters. (2) Cosmogenic neutrinos from cluster accretion shocks should have a redshift distribution distinct from AGN or GRB models, potentially detectable by upcoming facilities such as GRAND. (3) The hard radio background from shock‑accelerated electrons should be observable with low‑frequency arrays (e.g., SKA‑Low), providing a direct probe of the magnetic field strength in large‑scale structure shocks. (4) Heavy‑element dominance at the highest energies and a gradual hardening of the spectrum around 100 EeV are natural outcomes of the hierarchical model.

In summary, the authors demonstrate that the power available in large‑scale accretion shocks far exceeds the energy budget required for UHECRs, that plausible magnetic‑field amplification mechanisms can provide the necessary rigidity, and that the same shocks can simultaneously account for the observed diffuse radio background. The hierarchical shock model thus offers a self‑consistent, testable explanation for the full cosmic‑ray spectrum, linking low‑energy Galactic processes to the most energetic particles in the Universe. Future multi‑messenger observations—UHECR composition and anisotropy, high‑energy neutrinos, and low‑frequency radio surveys—will be decisive in confirming or refuting this scenario.