No Arabic abstract
Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code {sc arepo} for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.
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.
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 middle-aged supernova remnant (SNR) W44 has recently attracted attention because of its relevance regarding the origin of Galactic cosmic-rays. The gamma-ray missions AGILE and Fermi have established, for the first time for a SNR, the spectral continuum below 200 MeV which can be attributed to neutral pion emission. Confirming the hadronic origin of the gamma-ray emission near 100 MeV is then of the greatest importance. Our paper is focused on a global re-assessment of all available data and models of particle acceleration in W44, with the goal of determining on a firm ground the hadronic and leptonic contributions to the overall spectrum. We also present new gamma-ray and CO NANTEN2 data on W44, and compare them with recently published AGILE and Fermi data. Our analysis strengthens previous studies and observations of the W44 complex environment and provides new information for a more detailed modeling. In particular, we determine that the average gas density of the regions emitting 100 MeV - 10 GeV gamma-rays is relatively high (n= 250 - 300 cm^-3). The hadronic interpretation of the gamma-ray spectrum of W44 is viable, and supported by strong evidence. It implies a relatively large value for the average magnetic field (B > 10^2 microG) in the SNR surroundings, sign of field amplification by shock-driven turbulence. Our new analysis establishes that the spectral index of the proton energy distribution function is p1 = 2.2 +/- 0.1 at low energies and p2 = 3.2 +/- 0.1 at high energies. We critically discuss hadronic versus leptonic-only models of emission taking into account simultaneously radio and gamma-ray data. We find that the leptonic models are disfavored by the combination of radio and gamma-ray data. Having determined the hadronic nature of the gamma-ray emission on firm ground, a number of theoretical challenges remains to be addressed.
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 have long been regarded as sources of the Galactic cosmic rays up to petaelectronvolts, but convincing evidence is still lacking. In this work we explore the common origin of the subtle features of the cosmic ray spectra, such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical calculation shows that those features of cosmic ray spectra can be well reproduced in a scenario with e$^+$e$^-$ pair production by interactions between high energy cosmic rays and background photons in an environment similar to the young supernova remnant. The success of such a coherent explanation serves in turn as an evidence that at least a portion of cosmic rays might be accelerated at young supernova remnants.