[Paper Review] Cosmic ray escape from supernova remnants

S. Gabici|arXiv (Cornell University)|Aug 24, 2011

Astrophysics and Cosmic Phenomena4 references3 citations

TL;DR

This paper reviews cosmic ray (CR) escape mechanisms from supernova remnants (SNRs), emphasizing that CRs escape gradually during the Sedov phase as shock velocity decreases, with higher-energy particles escaping first. The key contribution is that CR spectra escaping SNRs are expected to be steep (index ~2.4), which reconciles SNR acceleration with the observed Galactic CR spectrum up to the knee energy (~3 PeV), especially when Alfven drift and magnetic field amplification are included in models.

ABSTRACT

Galactic cosmic rays are believed to be accelerated at supernova remnants via diffusive shock acceleration. Though this mechanism gives fairly robust predictions for the spectrum of particles accelerated at the shock, the spectrum of the cosmic rays which are eventually injected in the interstellar medium is more uncertain and depends on the details of the process of particle escape from the shock. Knowing the spectral shape of these escaping particles is of crucial importance in order to assess the validity of the supernova remnant paradigm for cosmic ray origin. Moreover, after escaping from a supernova remnant, cosmic rays interact with the surrounding ambient gas and produce gamma rays in the vicinity of the remnant itself. The detection of this radiation can be used as an indirect proof of the fact that the supernova remnant was indeed accelerating cosmic rays in the past.

Motivation & Objective

• To understand the spectral shape of cosmic rays escaping supernova remnants, which determines their contribution to the Galactic cosmic ray spectrum.
• To assess whether supernova remnants can be the dominant source of Galactic cosmic rays, particularly up to the knee energy (~3 PeV).
• To evaluate the role of particle escape mechanisms—especially diffusion and shock velocity evolution—in shaping the observed CR spectrum.
• To explore how gamma-ray emission from CR interactions with ambient gas (e.g., molecular clouds) can serve as indirect evidence of past CR acceleration in SNRs.
• To examine the impact of Alfven drift and magnetic field amplification on CR spectrum steepening, resolving discrepancies with observed spectra.

Proposed method

• Uses a Hillas-like confinement criterion (diffusion length ≤ SNR radius) to determine when particles escape during the Sedov phase.
• Applies scaling laws for SNR radius and shock velocity during the Sedov phase: R_sh ∝ t^{2/5}, u_sh ∝ t^{-3/5}.
• Estimates particle diffusion length l_d ∝ D/u_sh, with D ∝ E/B (Bohm diffusion), to model escape timescales.
• Derives energy-dependent escape times by equating l_d ≈ R_sh, yielding E_max ∝ t^{-1/5} for escaping particles.
• Incorporates Alfven drift by modifying the effective compression ratio felt by CRs, leading to steeper spectra than predicted by standard shock acceleration.
• Models CR diffusion and gamma-ray production from hadronic interactions in surrounding gas, focusing on TeV energy range detectable by Cherenkov telescopes.

Experimental results

Research questions

• RQ1 How do the spectral properties of cosmic rays escaping supernova remnants compare to the observed Galactic cosmic ray spectrum?
• RQ2 What physical mechanisms—especially diffusion and shock evolution—determine the energy-dependent escape of cosmic rays from SNRs?
• RQ3 Can the steepening of the cosmic ray spectrum due to Alfven drift explain the observed knee in the cosmic ray spectrum at ~3 PeV?
• RQ4 To what extent can gamma-ray emission from molecular clouds near SNRs serve as indirect evidence of past cosmic ray acceleration?
• RQ5 How do magnetic field amplification and non-resonant instabilities affect the predicted CR escape spectrum and detectability?

Key findings

• Cosmic rays escape gradually during the Sedov phase, with the highest-energy particles escaping first as the shock slows down.
• The maximum energy of confined particles decreases slowly over time, scaling as E_max ∝ t^{-1/5}, due to increasing diffusion length relative to SNR size.
• The spectrum of escaping cosmic rays is predicted to be steep, with a power-law index of approximately 2.4, which matches the observed Galactic cosmic ray spectrum near the knee.
• Incorporating Alfven drift—where scattering centers move differently upstream and downstream—can steepen the CR spectrum beyond the standard E^{-2} prediction, resolving a key discrepancy.
• Gamma-ray emission from hadronic interactions in nearby molecular clouds can serve as indirect evidence of past CR acceleration, especially if detected within ~100–200 pc and within ~10^4 years of the SNR's explosion.
• The excess of CRs near SNRs lasts for ~10^4 years, depending on local diffusion coefficients, and is most detectable in regions with enhanced CR streaming or magnetic field amplification.