No Arabic abstract
Supernova remnants are known to accelerate cosmic rays on account of their non-thermal emission of radio waves, X-rays, and gamma rays. Although there are many models for the acceleration of cosmic rays in Supernova remnants, the escape of cosmic rays from these sources is yet understudied. We use our time-dependent acceleration code RATPaC to study the acceleration of cosmic rays and their escape in post-adiabatic Supernova remnants and calculate the subsequent gamma-ray emission from inverse-Compton scattering and Pion decay. We performed spherically symmetric 1-D simulations in which we simultaneously solve the transport equations for CRs, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in a volume large enough to keep all CRs in the simulation. The transport equations for cosmic-rays and magnetic turbulence are coupled via the cosmic-ray gradient and the spatial diffusion coefficient of the cosmic rays, while the cosmic-ray feedback onto the shock structure can be ignored. Our simulations span 100kyrs, thus covering the remnants evolution until the beginning of the post-adiabatic phase. At later stages of the evolution cosmic rays over a wide range of energy can reside outside of the remnant, creating spectra that are softer than predicted by standard DSA and feature breaks in the 10-100 GeV-range. The total spectrum of cosmic rays released into the interstellar medium has a spectral index of s~2.4 above roughly 10 GeV which is close to that required by Galactic propagation models. We further find the gamma-ray luminosity to peak around an age of 4,000 years for inverse-Compton-dominated high-energy emission. Remnants expanding in low-density media emit generally more inverse-Compton radiation matching the fact that the brightest known supernova remnants - RCW86, Vela Jr, HESSJ1731-347 and RXJ1713.7-3946 - are all expanding in low density environments.
Supernova remnants (SNRs) are believed to accelerate particles up to high energies through the mechanism of diffusive shock acceleration (DSA). Except for direct plasma simulations, all modeling efforts must rely on a given form of the diffusion coefficient, a key parameter that embodies the interactions of energetic charged particles with the magnetic turbulence. The so-called Bohm limit is commonly employed. In this paper we revisit the question of acceleration at perpendicular shocks, by employing a realistic model of perpendicular diffusion. Our coefficient reduces to a power-law in momentum for low momenta (of index $alpha$), but becomes independent of the particle momentum at high momenta (reaching a constant value $kappa_{infty}$ above some characteristic momentum $p_{rm c}$). We first provide simple analytical expressions of the maximum momentum that can be reached at a given time with this coefficient. Then we perform time-dependent numerical simulations to investigate the shape of the particle distribution that can be obtained when the particle pressure back-reacts on the flow. We observe that, for a given index $alpha$ and injection level, the shock modifications are similar for different possible values of $p_{rm c}$, whereas the particle spectra differ markedly. Of particular interest, low values of $p_{rm c}$ tend to remove the concavity once thought to be typical of non-linear DSA, and result in steep spectra, as required by recent high-energy observations of Galactic SNRs.
The spectrum of cosmic ray protons and electrons released by supernova remnants throughout their evolution is poorly known, because of the difficulty in accounting for particle escape and confinement in the downstream of a shock front, where both adiabatic and radiative losses are present. Here we calculate the spectrum of cosmic ray protons released during the evolution of supernovae of different types, accounting for the escape from upstream and for adiabatic losses of particles advected downstream of the shock and liberated at later times. The same calculation is carried out for electrons. The magnetic field in the post-shock region is calculated by using an analytic treatment of the magnetic field amplification due to non--resonant and resonant streaming instability and their saturation. We find that when the field is the result of the growth of the cosmic-ray--driven non--resonant instability alone, the spectrum of electrons and protons released by a supernova remnant are indeed different, but such a difference becomes appreciable only at energies $gtrsim 100-1000$ GeV, while observations of the electron spectrum require such a difference to be present at energies as low as $sim 10$ GeV. An effect at such low energies requires substantial magnetic field amplification in the late stages of the supernova remnant evolution (shock velocity $ll 1000$ km/s), perhaps not due to streaming instability but hydrodynamical processes. We comment on the feasibility of such conditions and speculate on the possibility that the difference in spectral shape between electrons and protons may reflect either some unknown acceleration effect, or additional energy losses in cocoons around the sources.
We discuss recent observations of high energy cosmic ray positrons and electrons in the context of hadronic interactions in supernova remnants, the suspected accelerators of galactic cosmic rays. Diffusive shock acceleration can harden the energy spectrum of secondary positrons relative to that of the primary protons (and electrons) and thus explain the rise in the positron fraction observed by PAMELA above 10 GeV. We normalize the hadronic interaction rate by holding pion decay to be responsible for the gamma-rays detected by HESS from some SNRs. By simulating the spatial and temporal distribution of SNRs in the Galaxy according to their known statistics, we are able to then fit the electron (plus positron) energy spectrum measured by Fermi. It appears that IceCube has good prospects for detecting the hadronic neutrino fluxes expected from nearby SNRs.
The emission mechanism for hard $gamma$-ray spectra from supernova remnants (SNRs) is still a matter of debate. Recent multi-wavelength observations of TeV source HESS J1912+101 show that it is associated with an SNR with an age of $sim 100$ kyrs, making it unlikely produce the TeV $gamma$-ray emission via leptonic processes. We analyzed Fermi observations of it and found an extended source with a hard spectrum. HESS J1912+101 may represent a peculiar stage of SNR evolution that dominates the acceleration of TeV cosmic rays. By fitting the multi-wavelength spectra of 13 SNRs with hard GeV $gamma$-ray spectra with simple emission models with a density ratio of GeV electrons to protons of $sim 10^{-2}$, we obtain reasonable mean densities and magnetic fields with a total energy of $sim 10^{50}$ ergs for relativistic ions in each SNR. Among these sources, only two of them, namely SN 1006 and RCW 86, favor a leptonic origin for the $gamma$-ray emission. The magnetic field energy is found to be comparable to that of the accelerated relativistic ions and their ratio has a tendency of increase with the age of SNRs. These results suggest that TeV cosmic rays mainly originate from SNRs with hard $gamma$-ray spectra.
Despite significant efforts over the last decade, the origin of the cosmic ray positron excess has still not been unambiguously established. A popular class of candidates are pulsars or pulsar wind nebulae but these cannot account for the observed hard spectrum of cosmic ray antiprotons. We revisit the alternative possibility that the observed high-energy positrons are secondaries created by spallation in supernova remnants during the diffusive shock acceleration of the primary cosmic rays, which are further accelerated by the same shocks. The resulting source spectrum of positrons at high energies is then naturally harder than that of the primaries, as is the spectrum of other secondaries such as antiprotons. We present the first comprehensive investigation of the full parameter space of this model -- both the source parameters as well as those governing galactic transport. Various parameterisations of the cross-sections for the production of positrons and antiprotons are considered, and the uncertainty in the model parameters discussed. We obtain an excellent fit to the recent precision measurements by AMS-02 of cosmic ray protons, helium, positrons and antiprotons, as well as of various primary and secondary nuclei. The only notable deviation is an excess of antiprotons around ~10 GeV. This model thus provides an economical explanation of the spectra of all secondary species -- from a single well-motivated population of sources.