No Arabic abstract
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).
Within our Galaxy, supernova remnants are believed to be the major sources of cosmic rays up to the knee. However important questions remain regarding the share of the hadronic and leptonic components, and the fraction of the supernova energy channelled into these components. We address such question by the means of numerical simulations that combine a hydrodynamic treatment of the shock wave with a kinetic treatment of particle acceleration. Performing 3D simulations allows us to produce synthetic projected maps and spectra of the thermal and non-thermal emission, that can be compared with multi-wavelength observations (in radio, X-rays, and gamma-rays). Supernovae come in different types, and although their energy budget is of the same order, their remnants have different properties, and so may contribute in different ways to the pool of Galactic cosmic-rays. Our first simulations were focused on thermonuclear supernovae, like Tychos SNR, that usually occur in a mostly undisturbed medium. Here we present our 3D simulations of core-collapse supernovae, like the Cas A SNR, that occur in a more complex medium bearing the imprint of the wind of the progenitor star.
We present a model for the radio evolution of supernova remnants (SNRs) obtained by using three-dimensional (3D) hydrodynamic simulations, coupled with nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. We model the radio evolution of SNRs on a global level, by performing simulations for wide range of the relevant physical parameters, such as the ambient density, the supernova (SN) explosion energy, the acceleration efficiency and the magnetic field amplification (MFA) efficiency. We attribute the observed spread of radio surface brightnesses for corresponding SNR diameters to the spread of these parameters. In addition to our simulations of type Ia SNRs, we also considered SNR radio evolution in denser, nonuniform circumstellar environment, modified by the progenitor star wind. These simulations start with the mass of the ejecta substantially higher than in the case of a type Ia SN and presumably lower shock speed. The magnetic field is understandably seen as very important for the radio evolution of SNRs. In terms of MFA, we include both resonant and non-resonant modes in our large scale simulations, by implementing models obtained from first-principles, particle-in-cell (PIC) simulations and non-linear magnetohydrodynamical (MHD) simulations. We test the quality and reliability of our models on a sample consisting of Galactic and extragalactic SNRs. Our simulations give $Sigma-D$ slopes between -4 and -6 for the full Sedov regime. Recent empirical slopes obtained for the Galactic samples are around -5, while for the extragalactic samples are around -4.
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 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.