A nonlinear kinetic theory, combining cosmic-ray (CR) acceleration in supernova remnants (SNRs) with their gas dynamics, is used to re-examine the nonthermal properties of the remnant of SN 1987A for an extended evolutionary period of 5-50 yr. This spherically symmetric model is approximately applied to the different features of the SNR which consist of (i) a blue supergiant wind and bubble, and (ii) of the swept-up red supergiant (RSG) wind structures in the form of an H II region, an equatorial ring (ER), and an hourglass region. The RSG wind involves a mass loss rate that decreases significantly with elevation above and below the equatorial plane. The model adapts recent three-dimensional hydrodynamical simulations by Potter et al. in 2014 that use a significantly smaller ionized mass of the ER than assumed in the earlier studies by the present authors. The SNR shock recently swept up the ER, which is the densest region in the immediate circumstellar environment. Therefore, the expected gamma-ray energy flux density at TeV energies in the current epoch has already reached its maximal value of $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. This flux should decrease by a factor of about two over the next 10 years.

We analyze the results of recent measurements of Galactic cosmic ray (GCRs) energy spectra and the spectra of nonthermal emission from supernova remnants (SNRs) in order to determine their consistency with GCR origin in SNRs. It is shown that the measured primary and secondary CR nuclei energy spectra as well as the observed positron-to-electron ratio are consistent with the origin of GCRs up to the energy 10^17 eV in SNRs. Existing SNR emission data provide evidences for efficient CR production in SNRs accompanied by significant magnetic field amplification. In some cases the nature of the detected gamma-ray emission is difficult to determine because key SNR parameters are not known or poorly constrained.

We present the energy spectrum of antiproton cosmic ray (CR) component calculated on the basis of the nonlinear kinetic model of CR production in supernova remnants (SNR). The model includes reacceleration of already existing in interstellar medium antiprotons as well as creation of antiprotons in nuclear collisions of accelerated protons with gas nuclei and their subsequent acceleration by SNR shock. It is shown that antiprotons production in SNRs produces considerable effect in their resultant energy spectrum making it essentially flatter above 10 GeV so that the spectrum at TeV-energies increases by a factor of five. Calculated antiproton spectrum is well consistent with the PAMELA data, which correspond to energies below 100 GeV. As a consistency check we have also calculated within the same model the energy spectra of secondary nuclei and show that the measured boron-to-carbon ratio is consistent with the significant SNR contribution.

Recent observations of the diffuse Galactic \gr emission (DGE) by the {\it Fermi} Large Area Telescope ({\it Fermi}-LAT) have shown significant deviations, above a few GeV until about 100 GeV, from DGE models that use the GALPROP code for the propagation of cosmic ray (CR) particles outside their sources in the Galaxy and their interaction with the target distributions of the interstellar gas and radiation fields. The surplus of radiation observed is most pronounced in the inner Galaxy, where the concentration of CR sources is strongest. The present study investigates this "{\it Fermi}-LAT Galactic Plane Surplus" by estimating the \gr emission from the sources themselves, which is disregarded in the above DGE models. It is shown that indeed the expected hard spectrum of CRs, still confined in their sources (SCRs), can explain this surplus. The method is based on earlier studies regarding the so-called EGRET GeV excess which by now is generally interpreted as an instrumental effect. The contribution from SCRs is predicted to increasingly exceed the DGE models also above 100 GeV, up to \gr energies of about ten TeV, where the corresponding surplus exceeds the hadronic part of the DGE by about one order of magnitude. Above such energies the emission surplus should decrease again with energy due to the finite life-time of the assumed supernova remnant sources. Observations of the DGE in the inner Galaxy at 15 TeV with the Milagro \gr detector and, at TeV energies, with the ARGO-YBJ detector are interpreted to provide confirmation of a significant SCR contribution to the DGE.

The nature of the recently detected HE and VHE gamma-ray emission of Tycho's supernova remnant (SNR) is studied. A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is employed to investigate the properties of Tycho's SNR and their correspondence to the existing experimental data, taking into account that the ambient interstellar medium (ISM) is expected to be clumpy. It is demonstrated that the overall steep gamma-ray spectrum observed can be interpreted as the superposition of two spectra produced by the CR proton component in two different ISM phases: The first gamma-ray component, extending up to about $10^{14}$ eV, originates in the diluted warm ISM, whereas the second component, extending up to 100 GeV, comes from numerous dense, small-scale clouds embedded in this warm ISM. Given the consistency between acceleration theory and the observed properties of the nonthermal emission of Tycho's SNR, a very efficient production of nuclear CRs in Tycho's SNR is established. The excess of the GeV gamma-ray emission due to the clouds' contribution above the level expected in the case of a purely homogeneous ISM, is inevitably expected in the case of type Ia SNe.

The properties of the Galactic supernova remnant (SNR) SN 1006 are theoretically re-analyzed in the light of the recent H.E.S.S. results. Nonlinear kinetic theory is used to determine the momentum spectrum of cosmic rays (CRs) in space and time in the supernova remnant SN 1006. The physical parameters of the model - proton injection rate, electron-to-proton ratio and downstream magnetic field strength - are determined through a fit of the result to the observed spatially-integrated synchrotron emission properties. The only remaining unknown astronomical parameter, the circumstellar gas number density, is determined by a normalization of the amplitude of the gamma-ray flux to the observed amplitude. The bipolar morphology of both nonthermal X-ray and gamma-ray emissions is explained by the preferential injection of suprathermal nuclei and subsequent magnetic field amplification in the quasi-parallel regions of the outer supernova shock. The above parameters provide an improved fit to all existing nonthermal emission data, including the TeV emission spectrum recently detected by H.E.S.S., with the circumstellar hydrogen gas number density $N_H\approx 0.06 cm^{-3}$ close to values derived from observations of thermal X-rays. The hadronic and leptonic gamma-ray emissions are of comparable strength. The overall energy of accelerated CRs at the present epoch is of the order of 5 % of the total hydrodynamic explosion energy, and is predicted to rise with time by a factor of $\approx 2$. The relevance of CR escape from the SNR for the spectrum of the gamma-ray emission is demonstrated. The sum of the results suggests that SN 1006 is a CR source with a high efficiency of nuclear CR production, as required for the Galactic CR sources, both in flux as well as in cutoff energy.

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants is employed to re-examine the nonthermal properties of the remnant of SN 1987A for an extended evolutionary period of 5--100 yr. It is shown that an efficient production of nuclear CRs leads to a strong modification of the outer supernova remnant shock and to a large downstream magnetic field $B_\mathrm{d}\approx 20$ mG. The shock modification and the strong field are required to yield the steep radio emission spectrum observed, as well as to considerable synchrotron cooling of high energy electrons which diminishes their X-ray synchrotron flux. These features are also consistent with the existing X-ray observations. The expected \gr energy flux at TeV-energies at the current epoch is nearly $\epsilon_{\gamma}F_{\gamma}\approx 4\times 10^{-13}$ erg cm$^2$s$^{-1}$ under reasonable assumptions about the overall magnetic field topology and the turbulent perturbations of this field. The general nonthermal strength of the source is expected to increase roughly by a factor of two over the next 15 to 20 yrs; thereafter it should decrease with time in a secular form.

The properties of the - presumably - youngest Galactic supernova remnant (SNR) G1.9+0.3 are investigated within the framework of nonlinear kinetic theory of cosmic ray acceleration in SNRs. The observed angular size and expansion speed as well as the radio and X-ray emission measurements are used to determine relevant physical parameters of this SNR. Under the assumption that SNR G1.9+0.3 is the result of a Type Ia supernova near the Galactic center (at the distance d=8.5 kpc) the nonthermal properties are calculated. In particular, the expected TeV gamma-ray spectral energy density is predicted to be as low as $\epsilon_{\gamma}F_{\gamma} \approx 5\times 10^{-15}$ erg cm$^{-2}$ s$^{-1}$, strongly dependent ($F_{\gamma}\propto d^{-11}$) upon the source distance d.

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is employed to investigate the properties of SNR RX J1713.7-3946. Observations of the non-thermal radio and X-ray emission spectra as well as the H.E.S.S. measurements of the very high energy gamma-ray emission are used to constrain the astronomical and CR acceleration parameters of the system. It is argued that RX J1713.7-3946 is a core collapse supernova (SN) of type II/Ib with a massive progenitor, has an age of ~1600 yr and is at a distance of ~1 kpc. It is in addition assumed that the CR injection/acceleration takes place uniformly across the shock surface for this kind of core collapse SNR. The theory gives a consistent description for all the existing observational data, including the non-detection of thermal X-rays and the spatial correlation of the X-ray and gamma-ray emission in the remnant. Specifically it is shown that an efficient production of nuclear CRs, leading to strong shock modification and a large downstream magnetic field strength B_d ~140 mkG can reproduce in detail the observed synchrotron emission from radio to X-ray frequencies together with the gamma-ray spectral characteristics as observed by the H.E.S.S. telescopes. The calculations are consistent with RX J1713.7-3946 being an efficient source of nuclear cosmic rays.

Aims. The properties of the Galactic supernova remnant (SNR) RX J0852.0-4622 are theoretically analysed. Methods. An explicitly time-dependent, nonlinear kinetic model of cosmic ray (CR) acceleration in SNRs is used to describe the properties of SNR RX J0852.0-4622, the accelerated CRs and the nonthermal emission. The source is assumed to be at a distance of ~1 kpc in the wind bubble of a massive progenitor star. An estimate of the thermal X-ray flux in such a configuration is given. Results. We find that the overall synchrotron spectrum of RX J0852.0-4622 as well as the filamentary structures in hard X-rays lead to an amplified magnetic field B > 100 muG in the SNR interior. This implies that the leptonic very high energy (VHE) gamma-ray emission is suppressed, and that the VHE gamma-rays are hadronically dominated. The energy spectrum of protons produced over the life-time of the remnant until now may well reach ''knee'' energies. The derived gamma-ray morphology is consistent with the H.E.S.S. measurements. The amount of energy in energetic particles corresponds to about 35% of the hydrodynamic explosion energy. A remaining uncertainty concerns the thermal X-ray flux at 1 keV. A rough estimate, possibly not quite appropriate for the assumed wind bubble configuration, results in it being larger than the nonthermal flux at this energy. Conclusions. It is concluded that this SNR expanding into the wind bubble of a massive star in a dense gas environment can be a hadronic gamma-ray source that is consistent with all existing multi-wavelength constraints, except possibly the thermal X-ray emission.

The properties of the Galactic supernova remnant SN 1006 are theoretically reanalysed. Nonlinear kinetic theory is used to determine the acceleration efficiency of cosmic rays (CRs) in the supernova remnant SN 1006. The known range of astronomical parameters and the existing measurements of nonthermal emission are examined in order to define the values of the relevant physical parameters which determine the CR acceleration efficiency. It is shown that the parameter values -- proton injection rate, electron to proton ratio and downstream magnetic field strength -- are determined with the appropriate accuracy. In particular also the observed azimuthal variations in the gamma-ray morphology agree with the theoretical expectation. These parameter values, together with the reduction of the gamma-ray flux relative to a spherically symmetric acceleration geometry, allow a good fit to the existing data, including the recently detected TeV emission by H.E.S.S. SN 1006 represents the first example where a high efficiency of nuclear CR production, required for the Galactic CR sources, is consistently established.

The composition of the overall spectrum of cosmic rays (CRs) is studied under the assumption that ultra high energy CRs above the energy 10^{17} eV are produced at the shock created by the expanding cocoons around active galactic nuclei (AGNs). It is shown that the expected CR composition is characterised by two peaks in the energy dependence of the mean CR atomic number <A(\epsilon)>. The first one at the energy \epsilon \approx 10^{17} eV corresponds to the very end of the Galactic CR component, produced in supernova remnants (SNRs). It is followed by a sharp decrease of <A(\epsilon)> within the energy interval from 10^{17} to 10^{18} eV. This is a signature of the transition from Galactic to extragalactic CRs. The second peak, with <ln A> \approx 2, at energy \epsilon\approx 10^{19} eV, expected at the beginning of the GZK cutoff, is the signature of the CR production by the nonrelativistic cocoon shocks. The calculated CR composition is consistent with the existing data. The alternative scenario, which suggests reacceleration increasing the energy of CRs produced in SNRs by a factor of 30, is also examined.

Aims. The hadronic vs. leptonic origin of the gamma-ray emission from the Supernova Remnant RX J1713.7-3946 is discussed both in the light of new observations and from a theoretical point of view. Methods. The existing good spatial correlation of the gamma-ray emission and the nonthermal X-ray emission is analyzed theoretically. In addition, the recently published new H.E.S.S. observations define the energy spectrum more precisely, in particular at the high and low energy ends of the instrument's dynamical range. There now exist much more constraining X-ray observations from Suzaku that extend substantially beyond 10 keV. These new data are compared with the authors' previous theoretical predictions, both for dominant hadronic and for simple inverse Compton models. Results. Apart from the well-known MHD correlation between magnetic field strength and plasma density variations, emphasized by the wind-bubble-structure of the remnant, it is argued that the regions of magnetic field amplification also are correlated with enhanced densities of accelerated nuclear particles and the associated streaming instabilities. Therefore a correlation of nonthermal X-ray and \gr emission is not only possible but even to be expected for a hadronic emission scenario. A leptonic origin of the gamma-ray emission would require an implausibly uniform strength of the magnetic field. The observational and theoretical inferences about substantial field amplification in this remnant agree very well with the recent X-ray and \gr observations. Conclusions. All this argues strongly for the dominance of hadronic \grs in the \gr emission spectrum and a fortiori for an overwhelming contribution of nuclear cosmic ray particles to the nonthermal energy in this remnant.

We theoretically and phenomenologically investigate the question whether the gamma-ray emission from the remnants of the type Ia supernovae SN 1006, Tycho's SN and Kepler's SN can be the result of electron acceleration alone. The observed synchrotron spectra of the three remnants are used to determine the average momentum distribution of nonthermal electrons as a function of the assumed magnetic field strength. Then the inverse Compton emission spectrum in the Cosmic Microwave Background photon field is calculated and compared with the existing upper limits for the very high energy gamma-ray flux from these sources. It is shown that the expected interstellar magnetic fields substantially overpredict even these gamma-ray upper limits. Only rather strongly amplified magnetic fields could be compatible with such low gamma-ray fluxes. However this would require a strong component of accelerated nuclear particles whose energy density substantially exceeds that of the synchrotron electrons, compatible with existing theoretical acceleration models for nuclear particles and electrons. Even though the quantitative arguments are simplistic, they appear to eliminate simplistic phenomenological claims in favor of a inverse Compton gamma-ray scenario for these sources.

Cosmic ray (CR) acceleration at the shock created by the expanding cocoons around active galactic nuclei (AGNs) is studied. It is shown that above the energy $10^{18}$ eV the overall energy spectrum of CRs, produced during the AGN evolution and released in the intergalactic space, has the form $N\propto \epsilon^{-\gamma}$, with $\gamma\approx 2.6$, which extends up to $\epsilon_{max}\sim 10^{20}$ eV. It is concluded that cocoons shocks have to be considered as a main source of extragalactic CRs, which together with Galactic supernova remnants provide the observed CR spectrum.

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.

The consequences of a newly suggested value for the SN explosion energy 1.2x10^{51} erg are explored for the case of Tycho's supernova remnant (SNR). A nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs is employed to investigate the properties of Tycho's SNR and their correspondence to the existing experimental data. It is demonstrated that the large mean ratio between the radii of the contact discontinuity and the forward shock is consistent with the very effective acceleration of nuclear energetic particles at the forward shock. It is also argued that consistency of the value E_{sn} = 1.2x10^{51} erg with the gas dynamics, acceleration theory, and the existing gamma-ray measurements requires the source distance to be greater than 3.3 kpc. The corresponding ambient gas number density is lower than 0.4 cm^{-3}. Since the expected gamma-ray flux strongly depends on the source distance, F_{\gamma}\propto d^{-7}, a future experimental determination of the actual gamma-ray flux from Tycho's SNR will make it possible to determine the values of the source distance and of the mean ambient gas density. A simple inverse Compton model without a dominant population of nuclear CRs is not compatible with the present upper limit for the gamma-ray emission for any reasonable ambient interstellar B-field. Given the consistency between acceleration theory and overall, as well as internal, gas dynamics, a future gamma-ray detection would make the case for nuclear particle acceleration in Tycho's SNR incontrovertible in our view.

Explicitly time-dependent, nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) has been employed to investigate the properties of SNR RX J1713.7-3946. Observations of the nonthermal radio and X-ray emission spectra as well as earlier H.E.S.S. measurements of the very high energy $\gamma$-ray emission were used to constrain the astronomical and the particle acceleration parameters of the system. The model assumes that the object was a core collapse supernova (SN) with a massive progenitor, has an age of $\approx 1600$ yr and is at a distance of $\approx 1$ kpc. It is shown that an efficient production of nuclear CRs, leading to strong shock modification and a large downstream magnetic field strength $B_{\mathrm{d}}\sim100$ $\mu$G, can reproduce the observed synchrotron emission from radio to X-ray frequencies together with the \gr spectral characteristics as observed by the H.E.S.S. telescopes. Small-scale filamentary structures observed in nonthermal X-rays provide empirical confirmation for this field amplification scenario which leads to a strong depression of the inverse Compton and Bremsstrahlung fluxes. The results are compared with the latest H.E.S.S. observations.

Explicitly time-dependent, nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) has been used to investigate the properties of the very large SNR RX J0852.0-4622. The available observations do not clearly distinguish between a ``nearby'' (at $\sim 200$ pc) and a ``distant'' (at $\sim 1$ kpc) source scenario. Therefore two correspondingly different models were analyzed. While the 200 pc solution can not be a priory excluded, the 1 kpc solution turns out to be clearly preferable for physical reasons. It requires a core collapse supernova (SN) with a massive progenitor in a molecular cloud $\sim 4000$ yrs ago. The overall synchrotron spectrum and the filamentary structures in hard X-rays both consistently lead to an amplified magnetic field $B > 100 \mu$G in the SNR interior. This implies a suppression of the leptonic TeV $\gamma$-ray emission to about 1 percent of the flux measured by the H.E.S.S. telescope system which therefore must be hadronic, consistent with the theoretical solution. Up to the present the 1 kpc solution has already converted $\sim 10$ percent of the explosion energy into nonthermal energy, as expected for a Galactic CR source. Also the derived $\gamma$-ray morphology is consistent with the H.E.S.S. measurements. For the ``nearby'' solution the leptonic and hadronic $\gamma$ ray fluxes are in the ratio 1:10 which means that this case is also hadronically dominated. However, the magnetic field strength, consistent with the overall synchrotron spectrum, differs significantly from that derived from the X-ray filaments. Finally, the total mechanical energy released amounts to only $1.8 \times 10^{50}$ erg, uncomfortably low even for a core collapse event.

Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants is employed to calculate CR spectra. The magnetic field in SNRs is assumed to be significantly amplified by the efficiently accelerating nuclear CR component. It is shown that the calculated CR spectra agree in a satisfactory way with the existing measurements up to the energy $10^{17}$ eV. The power law spectrum of protons extends up to the energy $3\times 10^{15}$ eV with a subsequent exponential cutoff. It gives a natural explanation for the observed knee in the Galactic CR spectrum. The maximum energy of the accelerated nuclei is proportional to their charge number $Z$. Therefore the break in the Galactic CR spectrum is the result of the contribution of progressively heavier species in the overall CR spectrum so that at $10^{17}$ eV the CR spectrum is dominated by iron group nuclei. It is shown that this component plus a suitably chosen extragalactic CR component can give a consistent description for the entire Galactic CR spectrum.

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.

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants is employed to investigate the properties of the remnant SN 1987A. It is shown that a large downstream magnetic field ~10 mG is required to fit the existing observational data. Such a strong field together with the strong shock modification due to CR backreaction provides the steep and concave radioemission spectrum and considerable synchrotron cooling of high energy electrons which diminish their X-ray synchrotron flux below the observed Chandra flux which has to be considered as an upper limit for nonthermal X-ray emission. The expected gamma-ray energy flux at TeV-energies at the current epoch is 2x10^{-13} erg/(cm^2 s).

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.

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is employed to investigate the properties of SNR RX J1713.7-3946. Observations of the nonthermal radio and X-ray emission spectra as well as the H.E.S.S. measurements of the very high energy gamma-ray emission are used to constrain the astronomical and the particle acceleration parameters of the system. Under the assumptions that RX J1713.7-3946 was a core collapse supernova (SN) of type II/Ib with a massive progenitor, has an age of \approx 1600 yr and is at a distance of \approx 1 kpc, the theory gives indeed a consistent description for all the existing observational data. Specifically it is shown that an efficient production of nuclear CRs, leading to strong shock modification, and a large downstream magnetic field strength B_d ~ 100 mkG can reproduce in detail the observed synchrotron emission from radio to X-ray frequencies together with the gamma-ray spectral characteristics as observed by the H.E.S.S. telescopes. Small-scale filamentary structures observed in nonthermal X-rays provide empirical confirmation for the field amplification scenario which leads to a strong depression of the inverse Compton and Bremsstrahlung fluxes. Going beyond that and using a semi-empirical relation for young SNRs between the resulting CR pressure and the amplified magnetic field energy upstream of the outer SN shock as well as a moderate upper bound for the mechanical explosion energy, it is possible to also demonstrate the actual need for a considerable shock modification in RX J1713.7-3946. It is consistent with RX J1713.7-3946 being an efficient source of nuclear cosmic rays.

We use nonlinear kinetic theory to study the remnant dynamics and the particle acceleration as well as the properties of the nonthermal emission from the supernova remnant SN 1006. The known range of astronomical parameters is examined to determine whether it encompasses the existing synchrotron emission data. Given the present-day spatial extent and expansion rate of the object, it is shown that the hadronic gamma-ray flux is very sensitive to the ambient gas density N_H and that the existing H.E.S.S. upper limit requires N_H < 0.1 cm^{-3}. The strength of the amplified magnetic field downstream of the shock is about 150 \mu G.

It is shown that amplification of the magnetic field in supernova remnants (SNRs) occurs in all six objects where morphological measurements are presently available in the hard X-ray continuum at several keV. For the three archetypical objects (SN 1006, Cas A and Tycho's SNR) to which nonlinear time-dependent acceleration theory has been successfully applied up to now, the global theoretical and the local observational field strengths agree very well, suggesting in addition that all young SNRs exhibit the amplification effect as a result of very efficient acceleration of nuclear cosmic rays (CRs) at the outer shock. Since this appears to be empirically the case, we may reverse the argument and consider field amplification as a measure of nuclear CR acceleration and it has indeed been argued that acceleration in the amplified fields allows the CR spectrum from SNRs to reach the knee in the spectrum or, in special objects, even beyond. The above results are furthermore used to investigate the time evolution of field amplification in young SNRs. Although the uncertainties in the data do not allow precise conclusions regarding this point, they rather clearly show that the ratio of the magnetic field energy density and the kinetic energy density of gas flow into the shock is of the order of a few percent if the shock speed is high enough V_s > 10^3 km/s, and this ratio remains nearly constant during the SNR evolution. The escape of the highest energy nuclear particles from their sources becomes progressively important with age, reducing also the cutoff in the \pi^0 -decay gamma-ray emission spectrum with time after the end of the sweep-up phase. Simultaneously the leptonic gamma-ray channels will gain in relative importance with increasing age of the sources.

The time-dependent nonlinear kinetic theory for cosmic ray (CR) acceleration in supernova remnants (SNRs) is applied studying the properties of the synchrotron emission from SNRs, in particular, the surface brightness-diameter ($\Sigma-D$) relation. Detailed numerical calculations are performed for the expected range of the relevant physical parameters, namely the ambient density and the supernova explosion energy. The magnetic field in SNRs is assumed to be significantly amplified by the efficiently accelerating nuclear CR component. Due to the growing number of accelerated CRs the expected SNR luminosity increases during the free expansion phase, reaches a peak value at the beginning of the Sedov phase and then decreases again, since in this stage the overall CR number remains nearly constant, whereas the effective magnetic field diminishes with time. The theoretically predicted brightness-diameter relation in the radio range in the Sedov phase is close to $\Sigma_\mathrm{R}\propto D^{-17/4}$. It fits the observational data in a very satisfactory way. The observed spread of $\Sigma_\mathrm{R}$ at a given SNR size $D$ is the result of the spread of supernova explosion energies and interstellar medium densities.

The contribution of Source Cosmic Rays (SCRs), accelerated and still confined in Supernova Remnants (SNRs), to the diffuse high energy \gr emission above 1 GeV from the Galactic disk is studied. The \grs produced by the SCRs have a much harder spectrum than those generated by the Galactic Cosmic Rays (GCRs). Extending a previous paper, a simple SNR population synthesis is considered and the Inverse Compton emission from the SCR electrons is evaluated in greater detail. Then the combined spectrum of \gr emission from the Galactic Supernova Remnant population is calculated and this emission at low Galactic latitudes is compared with the diffuse \gr emission observed by the EGRET and ground based instruments. The average contribution of SCRs is comparable to the GCR contribution already at GeV energies, due to Supernovae of types II and Ib exploding into the wind bubbles of quite massive progenitor stars, and becomes dominant at \gr energies above 100 GeV. At TeV energies the dominant contribution is from SCRs in SNRs that expand into a uniform interstellar medium. In fact, the sum of hadronic and Inverse Compton \grs would exceed the limits given by the existing experimental data, unless the confinement time $T_{SN}$, i.e. the time until which SNRs confine the main fraction of accelerated SCRs, is as small as $T_{SN}\sim 10^4$ yr and the typical magnetic field strength in SNRs as large as 30 $\mu$G. Both situations are however possible as a result of field amplification through CR backreaction in the acceleration process. It is pointed out that accurate measurements of the low-latitude diffuse Galactic \gr spectrum at TeV-energies can serve as a unique consistency test for CR origin from the Supernova Remnant population as a whole.

It is shown that the nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) describes the shell-type nonthermal X-ray morphology of Cas A, obtained in Chandra observations, in a satisfactory way. The set of empirical parameters, like distance, source size and total energy release, is the same which reproduces the dynamical properties of the SNR and the spectral characteristics of the emission produced by CRs. The extremely small spatial scales in the observed morphological structures at hard X-ray energies are due to a large effective magnetic field $B_d \sim 500$ $\mu$G in the interior which is at the same time not only required but also sufficient to achieve the excellent agreement between the spatially integrated radio and X-ray synchrotron spectra and their calculated form. The only reasonably thinkable condition for the production of such a large effective field strength is a very efficiently accelerated nuclear CR component. Therefore the Chandra data confirm first of all the inference that Cas A indeed accelerates nuclear CRs with the high efficiency required for Cas A to be considered as a member of the main class of Galactic CR sources and, secondly, that the nonlinear kinetic theory of CR acceleration in SNRs is a reliable method to determine the magnetic field value in SNRs.

It is shown that the nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) fits the shell-type nonthermal X-ray morphology, obtained in Chandra observations, in a satisfactory way. The set of empirical parameters is the same which reproduces the dynamical properties of the SNR and the spectral characteristics of the emission produced by CRs. The extremely small spatial scales of the observed X-ray distribution are due to the large effective magnetic field $B_{\mathrm d}\sim 100$ $\mu$G in the interior, which is also required to give a good fit for the spatially integrated radio and X-ray synchrotron spectra. The only reasonably thinkable condition for the production of such a large effective field strength is an efficiently accelerated nuclear CR component. Therefore the Chandra data confirm the inference that SN 1006 indeed accelerates nuclear CRs with the high efficiency required for SNRs to be considered as the main Galactic CR sources.

The injection rate of suprathermal protons into the diffusive shock acceleration process should vary strongly over the surface of supernova remnant shocks. These variations and the absolute value of the injection rate are investigated. In the simplest case, like for SN 1006, the shock can be approximated as being spherical in a uniform large-scale magnetic field. The injection rate depends strongly on the shock obliquity and diminishes as the angle between the ambient field and the shock normal increases. Therefore efficient particle injection, which leads to conversion of a significant fraction of the kinetic energy at a shock surface element, arises only in relatively small regions near the "poles", reducing the overall CR production. The sizes of these regions depend strongly on the random background field and the Alfven wave turbulence generated due to CR streaming instability. For the cases of SN 1006 and Tycho's SNR they correspond to about 20, and for Cas A to between 10 and 20 percent of the entire shock surface. In first approximation, the CR production rate, calculated under the assumption of spherical symmetry, has therefore to be renormalized by this factor, while the shock as such remains roughly spherical.

We study the production of cosmic rays (CRs) in supernova remnants (SNRs), including the reacceleration of background galactic cosmic rays (GCRs) - thus refining the early considerations by Blandford & Ostriker (1980) and Wandel et al. (1987) - and the effects of the nuclear spallation inside the sources (the SNRs). This combines for the first time nuclear spallation inside CR sources and in the diffuse interstellar medium, as well as reacceleration, with the injection and subsequent acceleration of suprathermal particles from the postshock thermal pool. Selfconsistent CR spectra are calculated on the basis of the nonlinear kinetic model. It is shown that GCR reacceleration and CR spallation produce a measurable effect at high energies, especially in the secondary to primary (s/p) ratio, making its energy-dependence substantially flatter than predicted by the standard model. Quantitatively, the effect depends strongly upon the density of the surrounding circumstellar matter. GCR reacceleration dominates secondary CR production for a low circumstellar density. It increases the expected s/p ratio substantially and flattens its spectrum to an almost energy-independent form for energies larger than 100 GeV/n if the supernovae explode on average into a hot dilute medium with hydrogen number density $N_H=0.003$cm$^{-3}$. The contribution of CR spallation inside SNRs to the s/p ratio increases with increasing circumstellar density and becomes dominant for $N_H\gsim 1$ cm$^{-3}$, leading at high energies to a flat s/p ratio which is only by a factor of three lower than in the case of the hot medium. Measurements of the boron to carbon ratio at energies above 100 GeV/n could be used in comparison with the values predicted here as a consistency test for the supernova origin of the GCRs.

The nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to describe the relevant properties of Cassiopeia A (Cas A). In order to reproduce the SNR's observed size, expansion rate and thermal X-ray emission we employ a piecewise homogeneous model for the progenitor's circumstellar medium developed by Borkowski et al. (1996). It consists of a tenuous inner wind bubble, a dense shell of swept-up red supergiant wind material, and a subsequent red supergiant wind region. A quite large SNR interior magnetic field B_d approx 1 mG is required to give a good fit for the radio and X-ray synchrotron emission. The steep radio spectrum is consistent with efficient proton acceleration which produces a significant shock modification and leads to a steep electron spectrum at energies E_e < 1 GeV. The calculated integral gamma-ray flux from Cas A, F propto E_gamma^-1, is dominated by pi^0-decay gamma-rays due to relativistic protons. It extends up to roughly 30 TeV if CR diffusion is as strong as the Bohm limit. At TeV energies it satisfactorily agrees with the value 5.8 x 10^-13 cm^-2 s^-1 detected by the HEGRA collaboration.

A nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to describe the properties of Tycho's SNR. Observations of the expansion characteristics and of the nonthermal radio and X-ray emission spectra, assumed to be of synchrotron origin, are used to constrain the overall dynamical evolution and the particle acceleration parameters of the system, in addition to what is known from independent estimates of the distance and thermal X-ray observations. It is shown that a very efficient production of nuclear cosmic rays, leading to strong shock modification, and a large downstream magnetic field strength B_d approx 240muG are required to reproduce the observed synchrotron emission from radio to X-ray frequencies. This field strength is still well within the upper bound for the effective magnetic field, consistent with the acceleration process. The pi^0-decay gamma-ray flux turns out to be somewhat greater than the inverse Compton (IC) flux off the Cosmic Microwave Background at energies below 1 TeV, dominating it strongly at 10 TeV. The predicted TeV gamma-ray flux is consistent with but close to the very low upper limit recently obtained by HEGRA. A future detection at E_gamma ~ 10 TeV would clearly indicate hadronic emission.

The nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to describe the properties of the remnant of SN 1006. It is shown, that the theory fits the existing data in a satisfactory way within a set of parameters which is consistent with the idea that SN 1006 is a typical source of Galactic CR nucleons, although not necessarily of CR electrons. The adjusted parameters are those that are not very well determined by present theory or not directly amenable to astronomical observations. The calculated expansion law and the radio-, X-ray and gamma-ray emissions produced by the accelerated CRs in SN 1006 agree quite well with the observations. A rather large interior magnetic field B_d\approx 100 \muG is required to give a good fit for the radio and X-ray synchrotron emission. In the predicted TeV gamma-ray flux from SN 1006, the pi^0-decay gamma-rays, generated by the nuclear CR component, dominate over the inverse Compton (IC) gamma-rays, generated by the CR electrons in the cosmic microwave background. The predicted source morphology in high energy gamma-rays roughly corresponds to that of the synchrotron emission. The predicted integral gamma-ray flux F \propto e^-1 extends up to energies 100 TeV if CR diffusion is as strong as the Bohm limit. Only if the interior magnetic field is much lower in the SNR, B_d\approx 10 \muG, then the observed gamma-ray emission can be due to the accelerated electron component alone. In this case, not plausible physically in our view, the lowest permissible value of the electron to proton ratio is rather high, and the maximum individual energy and total energy content of accelerated nucleons so small, that SN 1006 can not be considered as a typical source of the nuclear Galactic CRs.

This report discusses the possibilities for observational tests of current models and ideas regarding the origin of the Galactic Cosmic Rays using facilities which are already available.

The contribution of the Source Cosmic Rays (SCRs), confined in Supernova Remnants, to the diffuse high energy \gr emission above 1 GeV from the Galactic disk is studied. \grs produced by the SCRs have a much harder spectrum compared with those generated by the Galactic Cosmic Rays which occupy a much larger residence volume uniformly. SCRs contribute less than 10% at GeV energies and become dominant at \gr energies above 100 GeV. The contributions from $\pi^0$-decay and Inverse Compton \grs have comparable magnitude and spectral shape, whereas the Bremsstrahlung component is negligible. At TeV energies the contribution from SCRs increases the expected diffuse \gr flux almost by an order of magnitude. It is shown that for the inner Galaxy the discrepancy between the observed diffuse intensity and previous model predictions at energies above a few GeV can be attributed to the SCR contribution.

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.

We have extended a simple model of nonlinear diffusive shock acceleration (Berezhko and Ellison 1999) to include the injection and acceleration of electrons and the production of photons from bremsstrahlung, synchrotron, inverse-Compton, and pion-decay processes. We argue that the results of this model, which is simpler to use than more elaborate ones, offer a significant improvement over test-particle, power-law spectra which are often used in astrophysical applications of diffusive shock acceleration. With an evolutionary supernova remnant (SNR) model to obtain shock parameters as functions of ambient interstellar medium parameters and time, we predict broad-band continuum photon emission from supernova remnants in general, and SN1006 in particular, showing that our results compare well with the more complete time-dependent and spherically symmetric nonlinear model of Berezhko, Ksenofontov, and Petukhov (1999). We discuss the implications nonlinear shock acceleration has for X-ray line emission, and use our model to describe how ambient conditions determine the TeV/radio flux ratio, an important parameter for gamma-ray observations of radio SNRs.

Abbreviated Abstract: A kinetic model of particle acceleration in supernova remnants (SNRs) is extended to study the cosmic ray (CR) and associated high energy gamma-ray production during SN shock propagation through the inhomogeneous circumstellar medium of a progenitor star that emits a strong wind. The wind forms a low-density bubble surrounded by a swept-up shell of interstellar matter. Gamma rays are produced as a result of decay of pions which in turn are the result of collisions of CRs with nuclei of the thermal plasma. The time evolution of the SNRs is followed numerically, taking into account the nonlinear backreaction of the accelerated CRs. Examples typical for SN type Ib and SN type II explosions are considered. Apart from the confirmation of the known fact that acceleration is extremely rapid and that the upper momentum cutoff is reached almost immediately after the explosion due to the high wind magnetic field, it is also shown that the CRs are accelerated with high efficiency. The CR energy spectrum, ultimately produced in the SNRs, has a power law form with index = 2.0 to 2.1 in a wide energy range up to a maximum energy which is about 10^14 eV if the CR diffusion coefficient is as small as the Bohm limiting value. However, the expected pi^0-decay gamma-ray flux is considerably lower than in the case of a uniform interstellar medium; a relatively high gamma-ray luminosity in the TeV band, detectable at distances of several kpc, is only expected in the case of a relatively dense ISM with a number density above 10/cc.

Spectrum and chemical composition of cosmic rays accelerated in supernova remnants are studied on a basis of kinetic approach. The cosmic ray transport equation with the Bohm diffusion coefficient has been numerically solved self-consistently with gas dynamic equations for the underlying flow. Comparison with observational results gives some indications that galactic cosmic rays are produced by supernovae shocks which expanded in the relatively low temperature partially ionized interstellar medium.