3D MHD modeling of the expanding remnant of SN 1987A : role of magnetic field and non-thermal radio emission

Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface to the expansion of its remnant through the surrounding nebula by 3D MHD simulations. The model considers the radiative cooling, the deviations from equilibrium of ionization, the deviation from temperature-equilibration between electrons and ions, and a plausible configuration of the pre-SN ambie…

The outer shock of SN1006

ROLE OF EJECTA CLUMPING AND BACK-REACTION OF ACCELERATED COSMIC RAYS IN THE EVOLUTION OF TYPE Ia SUPERNOVA REMNANTS

We investigate the role played by initial clumping of ejecta and by efficient acceleration of cosmic rays (CRs) in determining the density structure of the post-shock region of a Type Ia supernova remnant (SNR) through detailed 3D MHD modeling. Our model describes the expansion of a SNR through a magnetized interstellar medium (ISM), including the initial clumping of ejecta and the effects on shock dynamics due to back-reaction of accelerated CRs. The model predictions are compared to the observations of SN 1006. We found that the back-reaction of accelerated CRs alone cannot reproduce the observed separation between the forward shock (FS) and the contact discontinuity (CD) unless the energ…

Collisionless shock heating of heavy ions in SN 1987A

Astrophysical shocks at all scales, from those in the heliosphere up to the cosmological shock waves, are typically "collisionless", because the thickness of their jump region is much shorter than the collisional mean free path. Across these jumps, electrons, protons, and ions are expected to be heated at different temperatures. Supernova remnants (SNRs) are ideal targets to study collisionless processes because of their bright post-shock emission and fast shocks. Although optical observations of Balmer-dominated shocks in young SNRs showed that the post-shock proton temperature is higher than the electron temperature, the actual dependence of the post-shock temperature on the particle mass…

Linking gamma-ray spectra of supernova remnants to the cosmic ray injection properties in the aftermath of supernovae

The acceleration times of the highest-energy particles which emit gamma-rays in young and middle-age SNRs are comparable with SNR age. If the number of particles starting acceleration was varying during early times after the supernova explosion then this variation should be reflected in the shape of the gamma-ray spectrum. We use the solution of the non-stationary equation for particle acceleration in order to analyze this effect. As a test case, we apply our method to describe gamma-rays from IC443. As a proxy of the IC443 parent supernova we consider SN1987A. First, we infer the time dependence of injection efficiency from evolution of the radio spectral index in SN1987A. Then, we use the…

Radio polarization maps of shell-type SNRs II. Sedov models with evolution of turbulent magnetic field

Polarized radio emission has been mapped with great detail in several Galactic supernova remnants (SNRs), but has not yet been exploited to the extent it deserves. We have developed a method to model maps of the Stokes parameters for shell-like SNRs during their Sedov evolution phase. At first, 3-dimensional structure of a SNR has been computed, by modeling the distribution of the magnetohydrodynamic parameters and of the accelerated particles. The generation and dissipation of the turbulent component of magnetic field everywhere in SNR are also considered taking into account its interaction with accelerated particles. Then, in order to model the emission, we have used a generalization of t…

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Full list of authors: Abdalla, H.; Abe, H.; Acero, F.; Acharyya, A.; Adam, R.; Agudo, I; Aguirre-Santaella, A.; Alfaro, R.; Alfaro, J.; Alispach, C.; Aloisio, R.; Batista, R. Alves; Amati, L.; Amato, E.; Ambrosi, G.; Anguner, E. O.; Araudo, A.; Armstrong, T.; Arqueros, F.; Arrabito, L.; Asano, K.; Ascasibar, Y.; Ashley, M.; Backes, M.; Balazs, C.; Balbo, M.; Balmaverde, B.; Baquero Larriva, A.; Martins, V. Barbosa; Barkov, M.; Baroncelli, L.; de Almeida, U. Barres; Barrio, J. A.; Batista, P-, I; Becerra Gonzalez, J.; Becherini, Y.; Beck, G.; Tjus, J. Becker; Belmont, R.; Benbow, W.; Bernardini, E.; Berti, A.; Berton, M.; Bertucci, B.; Beshley, V; Bi, B.; Biasuzzi, B.; Biland, A.; Bissaldi, …

Constraints on the local interstellar magnetic field from non-thermal emission of SN1006

The synchrotron radio morphology of bilateral supernova remnants depends on the mechanisms of particle acceleration and on the viewing geometry. However, unlike X-ray and $\gamma$-ray morphologies, the radio emission does not depend on the cut-off region of the parent electron population, making it a simpler and more straightforward tool to investigate the physics of cosmic ray production in SNRs. Our aim is to derive from the radio morphology tight constraints on the direction of the local magnetic field and its gradient, and on the obliquity dependence of the electron injection efficiency. We perform a set of 3D MHD simulations describing the expansion of a spherical SNR through a magneti…

Post-adiabatic supernova remnants in an interstellar magnetic field: oblique shocks and non-uniform environment

Monthly notices of the Royal Astronomical Society 479(3), 4253 - 4270 (2018). doi:10.1093/mnras/sty1750

Three-dimensional Simulations from Supernovae to Their Supernova Remnants: The Dynamical and Chemical Evolution of Supernova 1987A

Modeling particle acceleration and non-thermal emission in supernova remnants

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 acceler…

On the Origin of Asymmetries in Bilateral Supernova Remnants

AIMS: We investigate whether the morphology of bilateral supernova remnants (BSNRs) observed in the radio band is determined mainly either by a non-uniform interstellar medium (ISM) or by a non-uniform ambient magnetic field. METHODS: We perform 3-D MHD simulations of a spherical SNR shock propagating through a magnetized ISM. Two cases of shock propagation are considered: 1) through a gradient of ambient density with a uniform ambient magnetic field; 2) through a homogeneous medium with a gradient of ambient magnetic field strength. From the simulations, we synthesize the synchrotron radio emission, making different assumptions about the details of acceleration and injection of relativisti…

Hydrodynamic simulations unravel the progenitor-supernova-remnant connection in SN 1987A

(Abridged) We aim at linking the dynamical and radiative properties of the remnant of SN 1987A to the geometrical and physical characteristics of the parent aspherical SN explosion and to the internal structure of its progenitor star. We performed 3D hydrodynamic simulations which describe the long-term evolution of SN 1987A from the onset of the SN to the full-fledged remnant at the age of 50 years, accounting for the pre-SN structure of the progenitor star. The simulations include all physical processes relevant for the complex phases of SN evolution and for the interaction of the SNR with the highly inhomogeneous ambient environment around SN 1987A. From the simulations, we synthesize ob…

Effects of non-uniform interstellar magnetic field on synchrotron X-ray and inverse-Compton γ-ray morphology of supernova remnants

Context. Observations of SuperNova Remnants (SNRs) in X-ray and γ-ray bands promise to contribute important information to our understanding of the kinematics of charged particles and magnetic fields in the vicinity of strong non-relativistic shocks and, therefore, the nature of Galactic cosmic rays. The accurate analysis of SNR images collected in different energy bands requires theoretical modeling of synchrotron and inverse Compton emission from SNRs. Aims. We develop a numerical code (remlight) to synthesize, from MHD simulations, the synchrotron radio, X-ray, and inverse Compton γ-ray emission originating in SNRs expanding in a non-uniform interstellar medium (ISM) and/or non-uniform i…

Predicted gamma-ray image of SN 1006 due to inverse Compton emission

We propose a method to synthesize the inverse Compton (IC) gamma-ray image of a supernova remnant starting from the radio (or hard X-ray) map and using results of the spatially resolved X-ray spectral analysis. The method is successfully applied to SN 1006. We found that synthesized IC gamma-ray images of SN 1006 show morphology in nice agreement with that reported by the H.E.S.S. collaboration. The good correlation found between the observed very-high energy gamma-ray and X-ray/radio appearance can be considered as an evidence that the gamma-ray emission of SN 1006 observed by H.E.S.S. is leptonic in origin, though the hadronic origin may not be excluded.

Thermal emission, shock modification, and X-ray emitting ejecta in SN 1006

Efficient particle acceleration can modify the structure of supernova remnants. In this context we present the results of the combined analysis of the XMM-Newton EPIC archive observations of SN 1006. We aim at describing 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. We perform a spatially resolved spectral analysis on a set of regions covering the whole rim of the shell and we exploit the results of the spectral analysis to produce a cou…