Gamma-ray emission expected from Kepler's SNR

Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to investigate the properties of Kepler's SNR and, in particular, to predict the gamma-ray spectrum expected from this SNR. Observations of the nonthermal radio and X-ray emission spectra as well as theoretical constraints for the total supernova (SN) explosion energy E_sn are used to constrain the astronomical and particle acceleration parameters of the system. Under the assumption that Kepler's SN is a type Ia SN we determine for any given explosion energy E_sn and source distance d the mass density of the ambient interstellar medium (ISM) from a fit to the observed SNR size and expansion speed. This makes it possible to make predictions for the expected gamma-ray flux. Exploring the expected distance range we find that for a typical explosion energy E_sn=10^51 erg the expected energy flux of TeV gamma-rays varies from 2x10^{-11} to 10^{-13} erg/(cm^2 s) when the distance changes from d=3.4 kpc to 7 kpc. In all cases the gamma-ray emission is dominated by \pi^0-decay gamma-rays due to nuclear CRs. Therefore Kepler's SNR represents a very promising target for instruments like H.E.S.S., CANGAROO and GLAST. A non-detection of gamma-rays would mean that the actual source distance is larger than 7 kpc.Comment: 6 pages, 4 figures. Accepted for publication in Astronomy and Astrophysics, minor typos correcte

New evidence for strong nonthermal effects in Tycho's supernova remnant

For the case of Tycho's supernova remnant (SNR) we present the relation between the blast wave and contact discontinuity radii calculated within the nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. It is demonstrated that these radii are confirmed by recently published Chandra measurements which show that the observed contact discontinuity radius is so close to the shock radius that it can only be explained by efficient CR acceleration which in turn makes the medium more compressible. Together with the recently determined new value $E_{sn}=1.2\times 10^{51}$ erg of the SN explosion energy this also confirms our previous conclusion that a TeV gamma-ray flux of $(2-5)\times 10^{-13}$ erg/(cm$^2$s) is to be expected from Tycho's SNR. Chandra measurements and the HEGRA upper limit of the TeV gamma-ray flux together limit the source distance $d$ to $3.3\leq d\leq 4$ kpc.Comment: 5 pages, 4 figures. Accepted for publication in Astrophysics and Space Science, Proc. of "The Multi-Messenger Approach to High-Energy Gamma-ray Sources (Third Workshop on the Nature of Unidentified High-Energy Sources)", Barcelona, July 4-7, 200

TeV Gamma Rays Expected from Supernova Remnants in Different Uniform Interstellar Media

Calculations of the expected TeV $\gamma$-ray emission, produced by accelerated cosmic rays (CRs) in nuclear collisions, from supernova remnants evolving in a uniform interstellar medium (ISM) are presented. The aim is to study the sensitivity of $\gamma$-ray production to a physical parameter set. Apart from its general proportionality to N_H, it is shown that the $\gamma$-ray production essentially depends upon the ratio of the CR diffusion coefficient $\kappa$ to a critical value $\kappa_{crit}=10(B_0/5 \mu{G})(N_H/0.3 {cm}^{-3})^{-1/3}\kappa_B$, where B_0 and N_H are the magnetic field and hydrogen number density of the ISM, and $\kappa_B$ denotes the Bohm diffusion coefficient. If $\kappa$ is of the same order or lower than $\kappa_{crit}$, then the peak TeV $\gamma$-ray flux in the Sedov evolutionary phase, normalized to a distance of 1 kpc, is about 10^{-10}(N_H/0.3 {cm}^{-3}) photons cm^{-2} s^{-1}. For a CR diffusion coefficient that is significantly larger than $\kappa_{crit}$, the CR cutoff energy is less than 10 TeV and the expected $\gamma$-ray flux at 1 TeV is considerably below the presently detectable level of 10^{-12} photons cm^{-2} s^{-1}. The same is of course true for a SNR in the rarified, so-called hot ISM.Comment: 9 pages, 2 figures, to appear in Astroparticle Physic

Cosmic Ray Acceleration by Supernova Shocks

We analyse the results of recent measurements of nonthermal emission from individual supernova remnants (SNRs) and their correspondence to the nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. It is shown that the theory fits these data in a satisfactory way and provides the strong evidences for the efficient CR production in SNRs accompanied by significant magnetic field amplification. Magnetic field amplification leads to considerable increase of CR maximum energy so that the spectrum of CRs accelerated in SNRs is consistent with the requirements for the formation of Galactic CR spectrum up to the energy ~10^17 eV.Comment: 30 pages, 8 figures, Solicited talk at 36th COSPAR Scientific Assembly; Beijing, Chin

Electrostatic Potentials in Supernova Remnant Shocks

Recent advances in the understanding of the properties of supernova remnant shocks have been precipitated by the Chandra and XMM X-ray Observatories, and the HESS Atmospheric Cerenkov Telescope in the TeV band. A critical problem for this field is the understanding of the relative degree of dissipative heating/energization of electrons and ions in the shock layer. This impacts the interpretation of X-ray observations, and moreover influences the efficiency of injection into the acceleration process, which in turn feeds back into the thermal shock layer energetics and dynamics. This paper outlines the first stages of our exploration of the role of charge separation potentials in non-relativistic electron-ion shocks where the inertial gyro-scales are widely disparate, using results from a Monte Carlo simulation. Charge density spatial profiles were obtained in the linear regime, sampling the inertial scales for both ions and electrons, for different magnetic field obliquities. These were readily integrated to acquire electric field profiles in the absence of self-consistent, spatial readjustments between the electrons and the ions. It was found that while diffusion plays little role in modulating the linear field structure in highly oblique and perpendicular shocks, in quasi-parallel shocks, where charge separations induced by gyrations are small, and shock-layer electric fields are predominantly generated on diffusive scales.Comment: 7 pages, 2 embedded figures, Accepted for publication in Astrophysics and Space Science, as part of the HEDLA 2006 conference proceeding

Diffusive propagation of cosmic rays from supernova remnants in the Galaxy. II: anisotropy

We investigate the effects of stochasticity in the spatial and temporal distribution of supernova remnants on the anisotropy of cosmic rays observed at Earth. The calculations are carried out for different choices of the diffusion coefficient D(E) for propagation in the Galaxy. The propagation and spallation of nuclei are taken into account. At high energies we assume that $D(E)\sim(E/Z)^{\delta}$, with $\delta=1/3$ and $\delta=0.6$ being the reference scenarios. The large scale distribution of supernova remnants in the Galaxy is modeled following the distribution of pulsars with and without accounting for the spiral structure of the Galaxy. Our calculations allow us to determine the contribution to anisotropy resulting from both the large scale distribution of SNRs in the Galaxy and the random distribution of the nearest remnants. The naive expectation that the anisotropy amplitude scales as D(E) is shown to be an oversimplification which does not reflect in the predicted anisotropy for any realistic distribution of the sources. The fluctuations in the anisotropy pattern are dominated by nearby sources, so that predicting or explaining the observed anisotropy amplitude and phase becomes close to impossible. We find however that the very weak energy dependence of the anisotropy amplitude below $10^{5}$ GeV and the rise at higher energies, can best be explained if the diffusion coefficient is $D(E)\sim E^{1/3}$. Faster diffusion, for instance with $\delta=0.6$, leads in general to an exceedingly large anisotropy amplitude. The spiral structure introduces interesting trends in the energy dependence of the anisotropy pattern, which qualitatively reflect the trend seen in the data. For large values of the halo size we find that the anisotropy becomes dominated by the large scale regular structure of the source distribution, leading indeed to a monotonic increase of $\delta_A$ with energy.Comment: 21 Pages, to appear in JCA

Cosmic-ray acceleration in supernova remnants: non-linear theory revised

A rapidly growing amount of evidences, mostly coming from the recent gamma-ray observations of Galactic supernova remnants (SNRs), is seriously challenging our understanding of how particles are accelerated at fast shocks. The cosmic-ray (CR) spectra required to account for the observed phenomenology are in fact as steep as $E^{-2.2}--E^{-2.4}$, i.e., steeper than the test-particle prediction of first-order Fermi acceleration, and significantly steeper than what expected in a more refined non-linear theory of diffusive shock acceleration. By accounting for the dynamical back-reaction of the non-thermal particles, such a theory in fact predicts that the more efficient the particle acceleration, the flatter the CR spectrum. In this work we put forward a self-consistent scenario in which the account for the magnetic field amplification induced by CR streaming produces the conditions for reversing such a trend, allowing --- at the same time --- for rather steep spectra and CR acceleration efficiencies (about 20%) consistent with the hypothesis that SNRs are the sources of Galactic CRs. In particular, we quantitatively work out the details of instantaneous and cumulative CR spectra during the evolution of a typical SNR, also stressing the implications of the observed levels of magnetization on both the expected maximum energy and the predicted CR acceleration efficiency. The latter naturally turns out to saturate around 10-30%, almost independently of the fraction of particles injected into the acceleration process as long as this fraction is larger than about $10^{-4}$.Comment: 24 pages, 5 figures, accepted for publication in JCA

Particle Acceleration at Relativistic Shocks

I review the current status of Fermi acceleration theory at relativistic shocks. I first discuss the relativistic shock jump conditions, then describe the non-relativistic Fermi mechanism and the differences introduced by relativistic flows. I present numerical calculations of the accelerated particle spectrum, and examine the maximum energy attainable by this process. I briefly consider the minimum energy for Fermi acceleration, and a possible electron pre-acceleration mechanism.Comment: 17 pages, 4 figures. To appear in "Relativistic Flows in Astrophysics", A.W. Guthmann, M. Georganopoulos, A. Marcowith and K. Manolokou, eds., Lecture Notes in Pysics, Springer Verla

Particle Acceleration in Cosmic Sites - Astrophysics Issues in our Understanding of Cosmic Rays

Laboratory experiments to explore plasma conditions and stimulated particle acceleration can illuminate aspects of the cosmic particle acceleration process. Here we discuss the cosmic-ray candidate source object variety, and what has been learned about their particle-acceleration characteristics. We identify open issues as discussed among astrophysicists. -- The cosmic ray differential intensity spectrum is a rather smooth power-law spectrum, with two kinks at the "knee" (~10^15 eV) and at the "ankle" (~3 10^18 eV). It is unclear if these kinks are related to boundaries between different dominating sources, or rather related to characteristics of cosmic-ray propagation. We believe that Galactic sources dominate up to 10^17 eV or even above, and the extragalactic origin of cosmic rays at highest energies merges rather smoothly with Galactic contributions throughout the 10^15--10^18 eV range. Pulsars and supernova remnants are among the prime candidates for Galactic cosmic-ray production, while nuclei of active galaxies are considered best candidates to produce ultrahigh-energy cosmic rays of extragalactic origin. Acceleration processes are related to shocks from violent ejections of matter from energetic sources such as supernova explosions or matter accretion onto black holes. Details of such acceleration are difficult, as relativistic particles modify the structure of the shock, and simple approximations or perturbation calculations are unsatisfactory. This is where laboratory plasma experiments are expected to contribute, to enlighten the non-linear processes which occur under such conditions.Comment: accepted for publication in EPJD, topical issue on Fundamental physics and ultra-high laser fields. From review talk at "Extreme Light Infrastructure" workshop, Sep 2008. Version-2 May 2009: adjust some wordings and references at EPJD proofs stag

Diffusive propagation of cosmic rays from supernova remnants in the Galaxy. I: spectrum and chemical composition

In this paper we investigate the effect of stochasticity in the spatial and temporal distribution of supernova remnants on the spectrum and chemical composition of cosmic rays observed at Earth. The calculations are carried out for different choices of the diffusion coefficient D(E) experienced by cosmic rays during propagation in the Galaxy. In particular, at high energies we assume that D(E)\sim E^{\delta}, with $\delta=1/3$ and $\delta=0.6$ being the reference scenarios. The large scale distribution of supernova remnants in the Galaxy is modeled following the distribution of pulsars, with and without accounting for the spiral structure of the Galaxy. We find that the stochastic fluctuations induced by the spatial and temporal distribution of supernovae, together with the effect of spallation of nuclei, lead to mild but sensible violations of the simple, leaky-box-inspired rule that the spectrum observed at Earth is $N(E)\propto E^{-\alpha}$ with $\alpha=\gamma+\delta$, where $\gamma$ is the slope of the cosmic ray injection spectrum at the sources. Spallation of nuclei, even with the small rates appropriate for He, may account for slight differences in spectral slopes between different nuclei, providing a possible explanation for the recent CREAM observations. For $\delta=1/3$ we find that the slope of the proton and helium spectra are $\sim 2.67$ and $\sim 2.6$ respectively at energies above 1 TeV (to be compared with the measured values of $2.66\pm 0.02$ and $2.58\pm 0.02$). For $\delta=0.6$ the hardening of the He spectra is not observed. We also comment on the effect of time dependence of the escape of cosmic rays from supernova remnants, and of a possible clustering of the sources in superbubbles. In a second paper we will discuss the implications of these different scenarios for the anisotropy of cosmic rays.Comment: 28 pages, To appear in JCA