Two-component $γ$-Ray Structure from the CR Sources Within Dense Clouds

Recent observations have revealed that several cosmic ray (CR) sources themselves exhibit pronounced double power-law features in their radiation spectra. Combined with the phenomenon of two-component structure in the observed CR energy spectrum supported by multi-messenger data, this raises a fundamental question: can the two-component structure of the cosmic ray energy spectrum and the double power-law feature of the gamma-ray radiation energy spectrum from supernova remnants be understood within a unified picture? In this study, we propose a two-component model that incorporates the re-acceleration of background ``sea" CR particles by astrophysical sources to systematically explain the formation of double power-law spectra within those sources. Our model successfully reproduces the gamma-ray observations of multiple CR sources. The results support that double power-law structures may be a generic feature of Galactic CR sources within crushed clouds. This work offers a new theoretical perspective on the origin and propagation of cosmic rays, and its predictions may be further tested with future observations of a larger sample of CR sources.

💡 Research Summary

In this paper the authors address a striking observational puzzle: many Galactic cosmic‑ray (CR) sources, especially supernova remnants (SNRs) interacting with dense molecular clouds, display gamma‑ray spectra that consist of two distinct power‑law segments. At the same time, the all‑particle CR spectrum measured near Earth shows a hardening above a few hundred GeV, a feature that has been interpreted either as a change in the source injection spectrum or as a transition in the propagation regime. The authors propose a unified “two‑component” framework that naturally produces both phenomena.

The core idea is that the high‑energy gamma‑ray emission is dominated by freshly injected, thermal particles that are accelerated at the SNR shock (the “thermal‑injection component”). These particles generate gamma rays through proton‑proton (pp) collisions (for hadrons) or inverse‑Compton scattering (for electrons). The low‑energy component, however, originates from the pre‑existing Galactic CR “sea”. When the SNR shock runs into a crushed cloud, the background CRs experience diffusive shock re‑acceleration and subsequent adiabatic compression in the radiative shell. This process lifts the low‑energy part of the spectrum, creates a spectral break at a momentum p_br set by ion‑neutral damping, and hardens the spectrum if the background CR slope is steeper than the shock‑induced index.

To quantify the background CR distribution, the authors adopt a spatially‑dependent propagation (SDP) model. In this model the diffusion coefficient Dxx(r,z,R) is anti‑correlated with the density of CR sources: regions rich in SNRs have a reduced diffusion coefficient. The diffusion coefficient is written as D0 F(r,z) βη(R/R0)δ0 F(r,z), where F(r,z) encodes the source distribution, δ0≈0.63 and D0≈4.8×10^28 cm^2 s⁻¹ are fitted to reproduce the B/C ratio, the proton spectrum, and other secondary‑to‑primary data using the GALPROP code. The resulting background spectrum f_GCR(p) serves as the upstream boundary condition for the re‑acceleration calculation.

Re‑acceleration is treated analytically following non‑linear shock acceleration theory. The steady‑state transport equation upstream and downstream of the shock yields a solution

f_acc(p)=α p^{−α}∫_{p_m}^{p}p’^{α−1}f_GCR(p’)dp',

with α=3v_sh/(v_sh−u_d), where v_sh is the shock speed and u_d the downstream flow speed. The minimum momentum p_m is set by the low‑energy cutoff of the Galactic sea, while the maximum momentum p_max is limited by the acceleration time versus the age of the remnant. For slow shocks propagating in dense media, ion‑neutral damping steepens the spectrum above a break momentum p_br, producing the observed low‑energy turnover.

After re‑acceleration, the gas in the radiative shell is compressed. The compression factor s=n_m/n_d (n_m is the post‑cooling density, n_d the immediate downstream density) typically lies between 10 and 30 for the sources considered. The compressed spectrum is

f_ad(p)=s^{2/3} f_acc(s^{−1/3}p),

which amplifies the particle density by s^{2/3} and consequently boosts the pion‑decay gamma‑ray emission.

The total gamma‑ray output is obtained by adding the thermal‑injection component (a power law N(E)∝E^{−α1} with α1≈1.5–2.5) and the re‑accelerated‑compressed component. For hadrons the authors compute the pp → π^0 → γγ channel; for electrons they include inverse‑Compton scattering on the interstellar radiation field and bremsstrahlung. Model parameters (distance, age, ambient density n0, break momentum p_br, compression factor s, etc.) are fitted to the spectral energy distributions of four well‑studied SNRs: Cas A, Tycho, W44, and IC 443. Using Markov Chain Monte Carlo (MCMC) techniques, they achieve excellent agreement with data from Fermi‑LAT, LHAASO, MAGIC, and other instruments across the GeV–PeV range.

Key results include:

• The low‑energy gamma‑ray slope is reproduced by the re‑accelerated‑compressed sea CRs, while the high‑energy tail is dominated by freshly accelerated particles.
• For W44 and IC 443, which are deeply embedded in dense clouds, the compression factor s≈20–30 is required, explaining the pronounced low‑energy bump that single‑component models cannot fit.
• The model simultaneously respects the measured B/C ratio and the local proton spectrum, confirming that the adopted SDP propagation is compatible with standard CR observables.
• The double‑power‑law feature emerges naturally from the superposition of two physically distinct populations, without invoking ad‑hoc spectral breaks in the source injection spectrum.

The authors argue that this two‑component picture provides a generic explanation for the widespread occurrence of double power‑law gamma‑ray spectra in SNRs interacting with crushed clouds. It also offers a possible link to the hardening observed in the all‑particle CR spectrum: the re‑accelerated sea component may contribute a harder sub‑population at Earth, especially if local sources reside in regions of reduced diffusion.

In the discussion, the paper emphasizes that the re‑acceleration and compression processes are expected whenever a shock encounters a dense, partially ionized medium, making the mechanism broadly applicable to other Galactic accelerators (e.g., pulsar wind nebulae interacting with filaments). Future high‑resolution gamma‑ray observations of a larger sample of SNRs, combined with detailed molecular‑cloud mapping, will allow direct measurement of the compression factor and the efficiency of re‑acceleration, providing stringent tests of the model.

In summary, the study presents a coherent theoretical framework that unifies the double‑power‑law gamma‑ray spectra of individual SNRs with the spectral hardening of the Galactic CR spectrum. By incorporating background CR re‑acceleration and adiabatic compression within crushed clouds, the authors demonstrate that these phenomena are natural outcomes of standard diffusive shock acceleration operating in realistic interstellar environments. The model’s success in reproducing multi‑instrument data for several archetypal remnants suggests that double‑power‑law structures may indeed be a generic hallmark of Galactic CR sources.