Muon Acceleration in Cosmic-ray Sources

Many models of ultra-high energy cosmic-ray production involve acceleration in linear accelerators located in Gamma-Ray Bursts magnetars, or other sources. These source models require very high accelerating gradients, $10^{13}$ keV/cm, with the minimum gradient set by the length of the source. At gradients above 1.6 keV/cm, muons produced by hadronic interactions undergo significant acceleration before they decay. This acceleration hardens the neutrino energy spectrum and greatly increases the high-energy neutrino flux. We rule out many models of linear acceleration, setting strong constraints on plasma wakefield accelerators and on models for sources like Gamma Ray Bursts and magnetars.

💡 Research Summary

The paper conducts a thorough examination of linear‑accelerator scenarios for ultra‑high‑energy cosmic‑ray (UHECR) production, focusing on environments such as gamma‑ray bursts (GRBs), magnetars, and plasma wakefield accelerators (PWFA). The authors begin by reviewing the conventional requirement that these sources must sustain accelerating gradients of order 10¹³ keV cm⁻¹ to reach the observed cosmic‑ray energies, a condition traditionally set by the source size. They then introduce a previously underappreciated effect: muons produced in hadronic interactions can themselves be accelerated by the same electric fields before they decay.

A key quantitative threshold emerges from their analysis: once the gradient exceeds roughly 1.6 keV cm⁻¹, the acceleration length available to a muon (set by its relativistically dilated lifetime) becomes sufficient for the muon to gain a non‑negligible amount of energy. The energy gain ΔE ≈ G · ℓ (where G is the gradient and ℓ the effective acceleration distance) can be comparable to or larger than the muon’s initial energy, especially for the very high Lorentz factors typical of UHECR environments. Consequently, the neutrinos produced in muon decay inherit this extra energy, hardening the neutrino spectrum from the canonical E⁻² shape toward a flatter E⁻¹·⁵ or even harder distribution.

The authors perform a two‑dimensional parameter scan over G and ℓ, identifying the region G·ℓ ≳ 10⁶ keV as the regime where the high‑energy (>PeV) neutrino flux would be amplified by an order of magnitude or more relative to standard expectations. They compare this prediction with current limits from IceCube, ANTARES, and KM3NeT, finding that such an enhanced flux is not observed. Therefore, any linear‑accelerator model that inevitably produces gradients above 1.6 keV cm⁻¹ over a sufficient distance is in tension with existing neutrino data.

Applying these constraints to specific source classes, the paper shows that while PWFA can theoretically reach gradients of 10¹⁴ keV cm⁻¹, realistic acceleration lengths are limited to a few meters, yielding G·ℓ far below the critical 10⁶ keV threshold. GRB internal shock models can achieve gradients of ~10⁶ keV cm⁻¹, but the required acceleration zones would need to extend over hundreds of kilometers, which is inconsistent with observed variability timescales and jet dynamics. Magnetar models, despite their extreme magnetic fields, also suffer from limited acceleration path lengths, preventing significant muon energization.

A detailed kinetic treatment is presented, incorporating muon energy loss mechanisms (synchrotron radiation, ionization, and hadronic interactions) alongside the electric‑field acceleration term. The resulting transport equation yields an analytic expression for the steady‑state muon energy distribution, which is then folded with the decay kinematics to produce the modified neutrino spectrum. The authors demonstrate that the spectral hardening is robust across a wide range of source parameters, provided the G·ℓ condition is met.

In the discussion, the paper emphasizes that the muon‑acceleration effect provides a powerful, model‑independent diagnostic for any linear accelerator hypothesis. The absence of the predicted high‑energy neutrino excess effectively rules out a large swath of parameter space for GRB, magnetar, and PWFA scenarios that rely on ultra‑high gradients over macroscopic distances. The authors also outline how next‑generation neutrino observatories (IceCube‑Gen2, KM3NeT‑ARCA, Baikal‑GVD) will sharpen these limits, potentially closing the remaining viable region or uncovering new physics.

In conclusion, the study establishes a stringent, physics‑driven bound on linear acceleration mechanisms in cosmic‑ray sources: gradients above ~1.6 keV cm⁻¹, when combined with realistic source sizes, inevitably lead to muon acceleration that would overproduce high‑energy neutrinos. Since such an overproduction is not observed, most existing linear‑accelerator models for UHECRs must be revised or abandoned, and any future proposals must explicitly account for the muon‑acceleration constraint.