Non-linear diffusive acceleration of heavy nuclei in supernova remnant shocks

We describe a semi-analytical approach to non-linear diffusive shock acceleration in the case in which nuclei other than protons are also accelerated. The structure of the shock is determined by the complex interplay of all nuclei, and in turn this shock structure determines the spectra of all components. The magnetic field amplification upstream is described as due to streaming instability of all nuclear species. The amplified magnetic field is then taken into account for its dynamical feedback on the shock structure as well as in terms of the induced modification of the velocity of the scattering centers that enters the particle transport equation. The spectra of accelerated particles are steep enough to be compared with observed cosmic ray spectra only if the magnetic field is sufficiently amplified and the scattering centers have high speed in the frame of the background plasma. We discuss the implications of this generalized approach on the structure of the knee in the all-particle cosmic ray spectrum, which we interpret as due to an increasingly heavier chemical composition above $10^{15}$eV. The effects of a non trivial chemical composition at the sources on the gamma ray emission from a supernova remnant when gamma rays are of hadronic origin are also discussed.

💡 Research Summary

The paper presents a semi‑analytical framework for non‑linear diffusive shock acceleration (NL‑DSA) that simultaneously treats protons and heavy nuclei (He, CNO, Fe, etc.) in supernova‑remnant (SNR) shocks. Traditional DSA models focus on protons and electrons, often adding heavy ions as a simple scaling factor, which fails to capture the observed composition‑dependent features of the Galactic cosmic‑ray spectrum, especially the “knee” around 10^15 eV. The authors address this gap by coupling three essential ingredients: (1) the dynamical feedback of all accelerated species on the shock structure, (2) magnetic‑field amplification upstream driven by the streaming instability of every nuclear component, and (3) the modification of the scattering‑center velocity (the Alfvénic drift) in the transport equation.

First, the model solves the coupled set of fluid equations (mass, momentum, and energy conservation) together with the kinetic transport equation for each species. The pressure contributed by each nuclear population enters the momentum balance, altering the compression ratio and the velocity profile across the shock. This self‑consistent treatment yields a shock structure that reflects the combined influence of all species rather than a proton‑dominated approximation.

Second, the authors compute the growth rate of resonant Alfvén waves generated by the streaming of all accelerated particles. The amplified magnetic field reduces the diffusion coefficient, thereby increasing the maximum attainable energy for each species. Importantly, the amplified field also contributes to the total pressure, further modifying the shock dynamics. The model shows that a magnetic amplification factor of order ten or more is required for the resulting particle spectra to be as steep as the observed Galactic cosmic‑ray spectrum (spectral index ≈ 2.2–2.4 rather than the canonical 2.0 of test‑particle DSA).

Third, the drift of scattering centers relative to the background plasma is explicitly included. Because the Alfvén speed in the amplified field can be a sizable fraction of the upstream flow speed, particles effectively experience a reduced compression ratio. This Alfvénic drift steepens the spectra of all species, bringing the theoretical predictions into agreement with measurements from balloon‑borne and satellite experiments.

The combined effect of strong magnetic amplification and significant Alfvénic drift leads to a natural explanation of the knee. Heavy nuclei, owing to their larger charge (Z), achieve higher rigidities for the same acceleration time, pushing their cut‑off energies to ∼ Z × 10^14 eV. As the shock evolves, the composition of the accelerated population becomes progressively heavier, and the all‑particle spectrum exhibits a gradual steepening at the knee, consistent with observations that indicate an increasingly heavy composition above 10^15 eV.

Beyond the particle spectra, the paper explores the implications for γ‑ray emission from SNRs of hadronic origin. The authors calculate the π^0‑decay γ‑ray flux assuming the derived multi‑species spectra and find that a heavy‑nuclei‑dominated population produces a harder γ‑ray spectrum at TeV energies compared with a pure‑proton scenario. This effect offers a potential diagnostic: high‑resolution γ‑ray observations (e.g., with CTA) could infer the source composition by measuring spectral curvature and cutoff features.

Numerically, the semi‑analytical method iterates between the fluid solution and the kinetic spectra until convergence. The resulting spectra reproduce the observed Galactic cosmic‑ray slope, the knee position, and predict a characteristic γ‑ray signature. The authors emphasize that the model’s predictive power hinges on two key parameters: the level of magnetic field amplification (set by the efficiency of the streaming instability) and the Alfvénic drift speed (determined by the amplified field strength and plasma density).

In summary, this work extends NL‑DSA theory to a realistic multi‑component plasma, demonstrating that the interplay between heavy nuclei, magnetic‑field growth, and scattering‑center drift can simultaneously explain the steep Galactic cosmic‑ray spectrum, the composition‑driven knee, and the high‑energy γ‑ray emission from SNRs. The framework provides a robust platform for interpreting forthcoming data from next‑generation cosmic‑ray and γ‑ray observatories, and it underscores the importance of treating heavy nuclei on an equal footing with protons in models of astrophysical particle acceleration.