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.

We present very-high-resolution 1D MHD simulations of the late-stage supernova remnants (SNR). In the post-adiabatic stage the magnetic field has an important and significant dynamical effect on the shock dynamics, the flow structure, and hence the acceleration and emission of cosmic rays. We find that the tangential component of the magnetic field provides pressure support that to a fair degree prevents the collapse of the radiative shell and thus limits the total compression ratio of the partially or fully radiative forward shock. A consequence is that the spectra of cosmic rays would not be as hard as in hydrodynamic simulations. We also investigated the effect on the flow profiles of the magnetic-field inclination and a large-scale gradient in the gas density and/or the magnetic field. A positive density gradient shortens the evolutionary stages whereas a shock obliquity lowers the shock compression. The compression of the tangential component of the magnetic field leads to its dominance in the downstream region of post-adiabatic shocks for a wide range of orientation of the upstream field, which may explain why one preferentially observes tangential radio polarization in old SNRs. As most cosmic rays are produced at late stages of SNR evolution, the post-adiabatic phase and the influence of the magnetic field during it are most important for modeling the cosmic-ray acceleration at old SNRs and the gamma-ray emission from late-stage SNRs interacting with clouds.

Gamma-rays from hadronic collisions are expected from supernova remnants (SNRs) located near molecular clouds. The temperature on the shock interacting with the dense environment quickly reaches $10^5$ K. The radiative losses of plasma become essential in the evolution of SNRs. They decrease the thermal pressure and essentially increase the density behind the shock. The presence of ambient magnetic field may considerably alter the behavior of the post-adiabatic SNRs comparing to hydrodynamic scenario. In the present paper, the magneto-hydrodynamic simulations of radiative shocks in magnetic field are performed. High plasma compression due to the radiative losses results also in the prominent increase of the strength of the tangential component of magnetic field behind the shock and the decrease of the parallel one. If the strength of the tangential field before the shock is higher than about $3\mathrm{\mu G}$ it prevents formation of the very dense thin shell. The higher the strength of the tangential magnetic field the larger the thickness and the lower the maximum density in the radiative shell. Parallel magnetic field does not affect the distribution of the hydrodynamic parameters behind the shock. There are almost independent channels of energy transformations: radiative losses are due to the thermal energy, the magnetic energy increases by reducing the kinetic energy. The large density and high strength of the perpendicular magnetic field in the radiative shells of SNRs should result in considerable increase of the hadronic gamma-ray flux comparing to the leptonic one.

Images of SN1006 have a number of important properties. For instance, the bright limbs coincide spatially in various bands, they have different brightness, and the contrast of brightness varies from radio to gamma-rays. The reasons of such properties and the role of the magnetic field strength are discussed. Simple, almost model-independent methods and analytical approximations for the derivation of the strength of magnetic field from the comparison of radio, X-rays and TeV images of SNR are presented. The methods require the TeV image to be well resolved and accurate, at least to the level of the radio and X-ray maps, in order to put reasonable constraints on magnetic field. If we apply it to the present HESS data, they limit the strength of magnetic field in limbs of SN1006 to values lower than few hundreds micro Gauss. If applied to the Fermi-LAT band, the model predicts same position and same ratio of the surface brightness for GeV photons as for the radio band. We conclude that TeV and GeV future high-resolution data may be very informative about the MF of SN1006.