Context. Supernova remnants are known as sources of Galactic cosmic rays for their nonthermal emission of radio waves, X-rays, and gamma rays. However, the observed soft broken power-law spectra are ...hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvénic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a comoving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces nonlinearity in addition to that imposed by cosmic-ray feedback.
Context. Supernova remnants (SNRs) are thought to be the primary candidates for the sources of Galactic cosmic rays. According to the diffusive shock acceleration theory, SNR shocks produce a ...power-law spectrum with an index of s = 2, perhaps nonlinearly modified to harder spectra at high energy. Observations of SNRs often indicate particle spectra that are softer than that and show features not expected from classical theory. Known drawbacks of the standard approach are the assumption that SNRs evolve in a uniform environment, and that the reverse shock does not accelerate particles. Relaxing these assumptions increases the complexity of the problem, because one needs reliable hydrodynamical data for the plasma flow as well as good estimates for the magnetic field (MF) at the reverse shock. Aims. We show that these two factors are especially important when modeling young core-collapse SNRs that evolve in a complicated circumstellar medium shaped by the winds of progenitor stars. Methods. We used high-resolution numerical simulations for the hydrodynamical evolution of the SNR. Instead of parametrizations of the MF profiles inside the SNR, we followed the advection of the frozen-in MF inside the SNR, and thus obtained the B-field value at all locations, in particular at the reverse shock. To model cosmic-ray acceleration we solved the cosmic-ray transport equation in test-particle approximation. Results. We find that the complex plasma-flow profiles of core-collapse SNRs significantly modify the particle spectra. Additionally, the reverse shock strongly affects the emission spectra and the surface brightness.
Context. The radio spectra of many shell-type supernova remnants show deviations from those expected on theoretical grounds. Aims. In this paper we determine the effect of stochastic reacceleration ...on the spectra of electrons in the GeV band and at lower energies, and we investigate whether reacceleration can explain the observed variation in radio spectral indices. Methods. We explicitely calculated the momentum diffusion coefficient for 3 types of turbulence expected downstream of the forward shock: fast-mode waves, small-scale non-resonant modes, and large-scale modes arising from turbulent dynamo activity. After noting that low-energy particles are efficiently coupled to the quasi-thermal plasma, a simplified cosmic-ray transport equation can be formulated and is numerically solved. Results. Only fast-mode waves can provide momentum diffusion fast enough to significantly modify the spectra of particles. Using a synchrotron emissivity that accurately reflects a highly turbulent magnetic field, we calculated the radio spectral index and find that soft spectra with index α ≲ − 0.6 can be maintained over more than 2 decades in radio frequency, even if the electrons experience reacceleration for only one acceleration time. A spectral hardening is possible but considerably more frequency-dependent. The spectral modification imposed by stochastic reacceleration downstream of the forward shock depends only weakly on the initial spectrum provided by, e.g., diffusive shock acceleration at the shock itself.
► We study cosmic ray acceleration during free expansion phase of the SNR. ► We consider the acceleration at both the reverse and the forward shock of the SNR. ► We show the expected observable ...features of SNRs with acceleration at the RS. ► We present a wider study of the particle spectra and their emission in type-Ia SNe.
We study cosmic-ray acceleration in young Type Ia Supernova Remnants (SNRs) by means of test-particle diffusive shock acceleration theory and 1-D hydrodynamical simulations of their evolution. In addition to acceleration at the forward shock, we explore the
particle acceleration at the reverse shock in the presence of a possible substantial magnetic field, and consequently the impact of this acceleration on the particle spectra in the remnant. We investigate the time evolution of the spectra for various time-dependent profiles of the magnetic field in the shocked region of the remnant. We test a possible influence on particle spectra of the Alfvénic drift of scattering centers in the precursor regions of the shocks. In addition, we study the radiation spectra and morphology in a broad band from radio to gamma-rays. It is demonstrated that the reverse shock contribution to the cosmic-ray particle population of young Type Ia SNRs
may be significant, modifying the spatial distribution of particles and noticeably affecting the volume-integrated particle spectra
in young SNRs. In particular spectral structures may arise in test-particle calculations that are often discussed as signatures of non-linear cosmic-ray modification of shocks. Therefore, the spectrum and morphology of emission, and their time evolution,
differ from pure forward-shock solutions.
Context. Supernova remnants (SNRs) are thought to be the main source of Galactic cosmic rays (CRs) up to the “knee” in CR spectrum. During the evolution of a SNR, the bulk of the CRs are confined ...inside the SNR shell. The highest-energy particles leave the system continuously, while the remaining adiabatically cooled particles are released when the SNR has expanded sufficiently and decelerated so that the magnetic field at the shock is no longer able to confine them. Particles escaping from the parent system may interact with nearby molecular clouds, producing γ-rays in the process via pion decay. The soft gamma-ray spectra observed for a number of SNRs interacting with molecular clouds, however, challenge current theories of non-linear particle acceleration that predict harder spectra. Aims. We study how the spectrum of escaped particles depends on the time-dependent acceleration history in both Type Ia and corecollapse SNRs, as well as on different assumptions about the diffusion coefficient in the vicinity of the SNR. Methods. We solve the CR transport equation in a test-particle approach combined with numerical simulations of SNR evolution. Results. We extend our method for calculating the CR acceleration in SNRs to trace the escaped particles in a large volume around SNRs. We calculate the evolution of the spectra of CRs that have escaped from a SNR into a molecular cloud or dense shell for two diffusion models. We find a strong confinement of CRs in a close region around the SNR, and a strong dilution effect for CRs that were able to propagate out as far as a few SNR radii.
Context. RX J1713.7–3946 is the brightest shell-type supernova remnant (SNR) of the TeV γ-ray sky. Earlier Fermi-LAT results on low energy γ-ray emission suggested that, despite large uncertainties ...in the background determination, the spectrum is inconsistent with a hadronic origin. Aims. We update the GeV-band spectra using improved estimates for the diffuse Galactic γ-ray emission and more than double the volume of data. We further investigate the viability of hadronic emission models for RX J1713.7–3946. Methods. We produced a high-resolution map of the diffuse Galactic γ-ray background corrected for the HI self-absorption and used it in the analysis of more than five years worth of Fermi-LAT data. We used hydrodynamic scaling relations and a kinetic transport equation to calculate the acceleration and propagation of cosmic rays in SNR. We then determined spectra of hadronic γ-ray emission from RX J1713.7–3946, separately for the SNR interior and the cosmic-ray precursor region of the forward shock, and computed flux variations that would allow us to test the model with observations. Results. We find that RX J1713.7–3946 is now detected by Fermi-LAT with very high statistical significance, and the source morphology is best described by that seen in the TeV band. The measured spectrum of RX J1713.7–3946 is hard with index γ = 1.53 ± 0.07, and the integral flux above 500 MeV is F = (5.5 ± 1.1) × 10-9 photons cm-2 s-1. We demonstrate that scenarios based on hadronic emission from the cosmic-ray precursor region are acceptable for RX J1713.7–3946, and we predict a secular flux increase at a few hundred GeV at the level of around 15% over ten years, which may be detectable with the upcoming Cherenkov Telescope Array (CTA) observatory.
High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it ...is a young, type Ia SNR that has been well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho's SNR with 147 hr of VERITAS and 84 months of Fermi-LAT observations, which represent about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is 2.92 0.42stat 0.20sys. It is also softer than the spectral index in the GeV energy range, 2.14 0.09stat 0.02sys, measured in this study using Fermi-LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi-LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously.
We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGNs), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when ...multi-TeV gamma-rays from an AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background photons via inverse-Compton scattering to produce gamma-rays. Due to the deflection of the electron-positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218+304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around 10−14 G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.
We consider Vela Jr. as being the old Supernova Remnant (SNR) at the beginning of the transition from adiabatic to radiative stage of evolution. According to our model, Vela Jr. is situated outside ...Vela SNR at the distance of ∼600
pc and its age is 17500
yr. We model the high energy fluxes from Vela Jr. and its broadband spectrum. We find our results compatible with experimental data in radio waves, X- and γ-rays. Our hydrodynamical model of Vela Jr. explains the observed TeV γ-ray flux by hadronic mechanism. The proposed model does not contradict to the low density environment of the SNR and does not need extreme fraction of the explosion energy to be transferred to Cosmic Rays.
Context. Efficient particle acceleration can modify the structure of supernova remnants. We present the results of a combined analysis of the XMM-Newton EPIC archive observations of SN 1006. Aims. We ...attempt to describe the spatial distribution of the physical and chemical properties of the X-ray emitting plasma at the shock front. We investigate the contribution of thermal and non-thermal emission to the X-ray spectrum at the rim of the remnant to study how the acceleration processes affect the X-ray emitting plasma. Methods. We perform a spatially resolved spectral analysis of a set of regions covering the entire rim of the shell and we apply our results in producing a count-rate image of the “pure” thermal emission of SN 1006 in the 0.5-0.8 keV energy band (subtracting the non-thermal contribution). This image differs significantly from the total image in the same band, especially close to the bright limbs. Results. We find that thermal X-ray emission can be associated with the ejecta and study the azimuthal variation in the physical and chemical properties of the ejecta by identifying anisotropies in the temperature and chemical composition. By employing our thermal image, we trace the position of the contact discontinuity over the entire shell and compare it with that expected from 3D MHD models of SNRs with an unmodified shock. Conclusions. We conclude that the shock is modified everywhere in the rim and that the aspect angle between the interstellar magnetic field and the line of sight is significantly lower than $90^\circ$.