No Arabic abstract
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.
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.
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.
Balmer emission may be a powerful diagnostic tool to test the paradigm of cosmic ray (CR) acceleration in young supernova remnant (SNR) shocks. The width of the broad Balmer line is a direct indicator of the downstream plasma temperature. In case of efficient particle acceleration an appreciable fraction of the total kinetic energy of the plasma is channeled into CRs, therefore the downstream temperature decreases and so does the broad Balmer line width. This width also depends on the level of thermal equilibration between ions and neutral hydrogen atoms in the downstream. Since in general in young SNR shocks only a few charge exchange (CE) reactions occur before ionization, equilibration between ions and neutrals is not reached, and a kinetic description of the neutrals is required in order to properly compute Balmer emission. We provide a method for the calculation of Balmer emission using a self-consistent description of the shock structure in the presence of neutrals and CRs. We use a recently developed semi-analytical approach, where neutral particles, ionized plasma, accelerated particles and magnetic fields are all coupled together through the mass, momentum and energy flux conservation equations. The distribution of neutrals is obtained from the full Boltzmann equation in velocity space, coupled to Maxwellian ions through ionization and CE processes. The computation is also improved with respect to previous work thanks to a better approximation for the atomic interaction rates. We find that for shock speeds >2500km/s the distribution of broad neutrals never approaches a Maxwellian and its moments differ from those of the ionized component. These differences reflect into a smaller FWHM than predicted in previous calculations, where thermalization was assumed. The method presented here provides a realistic estimate of particle acceleration efficiency in Balmer dominated shocks.
It is widely believe that galactic cosmic rays are originated in supernova remnants (SNRs) where they are accelerated by diffusive shock acceleration process at supernova blast waves driven by expanding SNRs. In recent theoretical developments of the diffusive shock acceleration theory in SNRs, protons are expected to accelerate in SNRs at least up to the knee energy. If SNRs are true generator of cosmic rays, they should accelerate not only protons but also heavier nuclei with right proportion and the maximum energy of heavier nuclei should be atomic mass (Z) times that of protons. In this work we investigate the implications of acceleration of heavier nuclei in SNRs on energetic gamma rays those are produced in hadronic interaction of cosmic rays with ambient matter. Our findings suggest that the energy conversion efficiency has to be nearly double for the mixed cosmic ray composition instead of pure protons to explain the observation and secondly the gamma ray flux above few tens of TeV would be significantly higher if cosmic rays particles can attain energies Z times of the knee energy in lieu of 200 TeV, as suggested earlier for non-amplified magnetic fields. The two stated maximum energy paradigm will be discriminated in future by the upcoming gamma ray experiments like Cherenkov Telescope array (CTA).
Context: Cosmic rays are thought to be accelerated at supernova remnant (SNR) shocks, but conclusive evidence is lacking. Aims: New data from ground-based gamma-ray telescopes and the Large Area Telescope on the Fermi Gamma-ray Space Telescope are used to test this hypothesis. A simple model for gamma-ray production efficiency is compared with measured gamma-ray luminosities of SNRs, and the GeV to TeV fluxes ratios of SNRs are examined for correlations with SNR ages. Methods: The supernova explosion is modeled as an expanding spherical shell of material that sweeps up matter from the surrounding interstellar medium (ISM). The accumulated kinetic energy of the shell, which provides the energy available for nonthermal particle acceleration, changes when matter is swept up from the ISM and the SNR shell decelerates. A fraction of this energy is assumed to be converted into the energy of cosmic-ray electrons or protons. Three different particle radiation processes---nuclear pion-production interactions, nonthermal electron bremsstrahlung, and Compton scattering---are considered. Results: The efficiencies for gamma-ray production by these three processes are compared with gamma-ray luminosities of SNRs. Our results suggest that SNRs become less gamma-ray luminous at >~ 10^4 yr, and are consistent with the hypothesis that supernova remnants accelerate cosmic rays with an efficiency of ~10% for the dissipation of kinetic energy into nonthermal cosmic rays. Weak evidence for an increasing GeV to TeV flux ratio with SNR age is found.