No Arabic abstract
Supernova remnants are known to accelerate cosmic rays (CRs) on account of their non-thermal emission of radio waves, X-rays, and gamma rays. However, the ability to accelerate CRs up to PeV-energies has yet to be demonstrated. The presence of cut-offs in the gamma-ray spectra of several young SNRs led to the idea that PeV energies might only be achieved during the very initial stages of a remnants evolution. We use the time-dependent acceleration code RATPaC to study the acceleration of cosmic rays in supernovae expanding into dense environments around massive stars, where the plentiful target material might offer a path to the detection of gamma-rays by current and future experiments. We performed spherically symmetric 1-D simulations in which we simultaneously solve the transport equations for cosmic rays, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in the test-particle limit. We investigated typical parameters of the circumstellar medium (CSM) in the freely expanding winds around red supergiant (RSG) and luminous blue variable (LBV) stars. The maximum achievable energy might be limited to sub-PeV energies despite strong magnetic fields close to the progenitor star that enhance turbulence-damping by cascading: we find a maximum CR energy of 100-200 TeV, reached within one month after explosion. The peak luminosity for a LBV progenitor is 1e43 erg/s (1e42 erg/s) at GeV (TeV) energies and, for a RSG progenitor, 1e41 erg/s (1e40 erg/s). All calculated SNe reach their peak gamma-ray luminosity after <~1 month and then fade at a rate ~1/t as long as the SN shock remains in the freely expanding wind of the progenitor. Potentially detectable gamma-ray signals can be expected in the Fermi-LAT waveband weeks to months after an explosion into a freely expanding wind.
We study the non-thermal emission from old shell-type supernova remnants (SNRs) on the frame of a time-dependent model. In this model, the time-dependent non-thermal spectra of both primary electrons and protons as well as secondary electron/positron ($e^{pm}$) pairs can be calculated numerically by taking into account the evolution of the secondary $e^{pm}$ pairs produced from proton-proton (p-p) interactions due to the accelerated protons collide with the ambient matter in an SNR. The multi-wavelength photon spectrum for a given SNR can be produced through leptonic processes such as electron/positron synchrotron radiation, bremsstrahlung and inverse Compton scattering as well as hadronic interaction. Our results indicate that the non-thermal emission of the secondary $e^{pm}$ pairs is becoming more and more prominent when the SNR ages in the radiative phase because the source of the primary electrons has been cut off and the electron synchrotron energy loss is significant for a radiative SNR, whereas the secondary $e^{pm}$ pairs can be produced continuously for a long time in the phase due to the large energy loss time for the p-p interaction. We apply the model to two old SNRs, G8.7$-$0.1 and G23.3$-$0.3, and the predicted results can explain the observed multi-wavelength photon spectra for the two sources.
According to the most popular model for the origin of cosmic rays (CRs), supernova remnants (SNRs) are the site where CRs are accelerated. Observations across the electromagnetic spectrum support this picture through the detection of non-thermal emission that is compatible with being synchrotron or inverse Compton radiation from high energy electrons, or pion decay due to proton-proton interactions. These observations of growing quantity and quality promise to unveil many aspects of CRs acceleration and require more and more accurate tools for their interpretation. Here, we show how multi-dimensional MHD models of SNRs, including the effects on shock dynamics due to back-reaction of accelerated CRs and the synthesis of non-thermal emission, turned out to be very useful to investigate the signatures of CRs acceleration and to put constraints on the acceleration mechanism of high energy particles. These models have been used to interpret accurately observations of SNRs in various bands (radio, X-ray and $gamma$-ray) and to extract from them key information about CRs acceleration.
Supernova remnants are believed to be the major contributors to Galactic cosmic rays. In this paper, we explore how the non-thermal emission from young remnants can be used to probe the production of energetic particles at the shock (both protons and electrons). Our model couples hydrodynamic simulations of a supernova remnant with a kinetic treatment of particle acceleration. We include two important back-reaction loops upstream of the shock: energetic particles can (i) modify the flow structure and (ii) amplify the magnetic field. As the latter process is not fully understood, we use different limit cases that encompass a wide range of possibilities. We follow the history of the shock dynamics and of the particle transport downstream of the shock, which allows us to compute the non-thermal emission from the remnant at any given age. We do this in 3D, in order to generate projected maps that can be compared with observations. We observe that completely different recipes for the magnetic field can lead to similar modifications of the shock structure, although to very different configurations of the field and particles. We show how this affects the emission patterns in different energy bands, from radio to X-rays and $gamma$-rays. High magnetic fields ($>100 mu$G) directly impact the synchrotron emission from electrons, by restricting their emission to thin rims, and indirectly impact the inverse Compton emission from electrons and also the pion decay emission from protons, mostly by shifting their cut-off energies to respectively lower and higher energies.
Supernova remnants (SNRs) are believed to be the major contributors to Galactic cosmic rays. The detection of non-thermal emission from SNRs demonstrates the presence of energetic particles, but direct signatures of protons and other ions remain elusive. If these particles receive a sizeable fraction of the explosion energy, the morphological and spectral evolution of the SNR must be modified. To assess this, we run 3D hydrodynamic simulations of a remnant coupled with a non-linear acceleration model. We obtain the time-dependent evolution of the shocked structure, impacted by the Rayleigh-Taylor hydrodynamic instabilities at the contact discontinuity and by the back-reaction of particles at the forward shock. We then compute the progressive temperature equilibration and non-equilibrium ionization state of the plasma, and its thermal emission in each cell. This allows us to produce the first realistic synthetic maps of the projected X-ray emission from the SNR. Plasma conditions (temperature, ionization age) can vary widely over the projected surface of the SNR, especially between the ejecta and the ambient medium owing to their different composition. This demonstrates the need for spatially-resolved spectroscopy. We find that the integrated emission is reduced with particle back-reaction, with the effect being more significant for the highest photon energies. Therefore different energy bands, corresponding to different emitting elements, probe different levels of the impact of particle acceleration. Our work provides a framework for the interpretation of SNR observations with current X-ray missions (Chandra, XMM-Newton, Suzaku) and with upcoming X-ray missions (such as Astro-H).
Shock waves associated with shell type supernova remnants are considered to be possible sites of cosmic ray acceleration. Since shocks are capable of accelerating electrons in addition to protons one anticipates both species to contribute to the high energy radiation expected from these objects. Adopting a simple model for particle acceleration we calculate in a self-consistent manner the time-dependent synchrotron and inverse Compton radiation of high energy electrons assumed either to be accelerated directly by the shock wave or to be injected at high energies as secondaries from the hadronic collisions of relativistic protons with the circumstellar material. We deduce that for standard supernova parameters the TeV flux produced from neutral pion decay is about the same order as the flux expected from directly accelerated electrons.