skip to content

Radiative feedback in molecular clouds from cooling supernova shock waves

1 / 2
  • Enlarged view: Figure 1
    Figure 1: Example of the cooling radiation emission plot (projection) in optical (E < 13.6 eV), UV (13.6 eV < E < 100 eV) and X-ray (E > 100 eV) energy bands. We used high-resolution ( ~ 0.1 pc) magnetohydrodynamical simulations with the adaptive mesh refinement code FLASH 4.3 from Seifried et al. 2018. The MC is formed self-consistently within its large-scale galactic environment following the non-equilibrium formation of H2 and CO including (self-) shielding and important heating and cooling processes. After the cloud is fully refined (11.9 Myr) we explode single SNe (1.53 Myr later) at distance d = 25 pc from the center-of-mass of cloud along the x-direction.
  • Enlarged view: Figure 2
    Figure 2: Example of the line ratios emission map (projection) at time. Left: S[II] λ6717/S [II] λ6731, right: S[II]/Hα. After around 0.3 Myr after the SN event the temperature behind the shock front drops to 10^5 K and the energy lost more due to recombination, creating a cooling region behind the shock front and producing collisionally ionized species like as [SII].

Ekaterina Makarenko

Massive stars (with mass > 8 Solar masses) end their lives dramatically with an explosion of incredible power: a supernova. Supernovae (SNe) inject ~ 1051 erg per SN into the interstellar medium (ISM) shaping the chemical, thermal and also dynamic evolution of galaxies. The gas which is shocked by SN is heated to 107 K and cools by emitting radiation at different energies. As such, SNe do not interact with the ISM purely through the physical act of shocking gas but also through the production of high energy photons (and cosmic rays) which are re-absorbed by the surrounding dense molecular gas. Typically, about 70% of the initially injected SN energy is lost by radiative cooling. The question is how to make simulations taking into account all the physics of this process?

Commonly, the former process (the shocking and shock heating) is always treated in numerical simulations of the SN-driven ISM. But the fate of the emitted cooling photons is usually neglected as the surrounding medium is treated to be optically thin with respect to the cooling radiation. That is why we are developing a post-processing module for FLASH code to calculate the cooling radiation from shock-heated gas in different energy bands ranging from X-rays to UV and Optical emission as shown on Figure 1. Using modern collision excitation data from MAPPING V (Sutherland et al. 2018), we can produce realistic maps of SNe emission taking into account angular resolution, radiation transfer effects, etc. The results of this work can be compared with real observations of supernova remnants. For instance, the time-dependent emission in different energy bands is analyzed for supernovae, which have exploded in different ISM environments. Considering that supernova remnants (SNRs) are key objects to study the chemical enrichment of the ISM, we can extend our model to use them as an excitation diagnostic tool for distinguishing among HII regions  objects photoionized by a harder radiation field (SNe and SNRs), see Figure 2.