No Arabic abstract
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.
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.
In the past few years, gamma-ray astronomy has entered a golden age. At TeV energies, only a handful of sources were known a decade ago, but the current generation of ground-based imaging atmospheric Cherenkov telescopes has increased this number to more than one hundred. At GeV energies, the Fermi Gamma-ray Space Telescope has increased the number of known sources by nearly an order of magnitude in its first 2 years of operation. The recent detection and unprecedented morphological studies of gamma-ray emission from shell-type supernova remnants is of great interest, as these analyses are directly linked to the long standing issue of the origin of the cosmic-rays. However, these detections still do not constitute a conclusive proof that supernova remnants accelerate the bulk of Galactic cosmic-rays, mainly due to the difficulty of disentangling the hadronic and leptonic contributions to the observed gamma-ray emission. In this talk, I will review the most relevant cosmic ray related results of gamma ray astronomy concerning supernova remnants.
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.
In the last few years several experiments have shown that the cosmic ray spectrum below the knee is not a perfect power-law. In particular, the proton and helium spectra show a spectral hardening by ~ 0.1-0.2 in spectral index at particle energies of ~ 200-300 GeV/nucleon. Moreover, the helium spectrum is found to be harder than that of protons by ~ 0.1 and some evidence for a similar hardening was also found in the spectra of heavier elements. Here we consider the possibility that the hardening may be the result of a dispersion in the slope of the spectrum of cosmic rays accelerated at supernova remnant shocks. Such a dispersion is indeed expected within the framework of non-linear theories of diffusive shock acceleration, which predict steeper (harder) particle spectra for larger (smaller) cosmic ray acceleration efficiencies.
We present results from an optical study of two young Balmer-dominated remnants of SNIa in the Large Magellanic Cloud, 0509$-$67.5 and 0519$-$69.0, in an attempt to search for signatures of efficient cosmic-ray (CR) acceleration. We combine proper motion measurements from HST with corresponding optical spectroscopic measurements of the H$alpha$ line at multiple rim positions from VLT/FORS2 and SALT/RSS and compare our results to published Balmer shock models. Analysis of the optical spectra result in broad H$alpha$ widths in the range of 1800-4000 km s$^{-1}$ for twelve separate Balmer-dominated filaments that show no evidence for forbidden line emission, the corresponding shock speeds from proper motion measurements span a range of 1700-8500 km s$^{-1}$. Our measured values of shock speeds and broad H$alpha$ widths in 0509$-$67.5 and 0519$-$69.0 are fit well with a Balmer shock model that does not include effects of efficient CR acceleration. We determine an upper limit of 7%/$chi$ (95% confidence) on the CR acceleration efficiency for our ensemble of data points, where $chi$ is the ionization fraction of the pre-shock gas. The upper limits on the individual remnants are 6%/$chi$ (0509$-$67.5) and 11%/$chi$ (0519$-$69.0). These upper limits are below the integrated CR acceleration efficiency in the Tycho supernova remnant, where the shocks predominantly show little H$alpha$ emission, indicating that Balmer-dominated shocks are not efficient CR accelerators.