Type Ia Supernova Explosion Models Research Articles

Type Ia supernova explosion models are inherently multidimensional

Theoretical and observational approaches to settling the important questions surrounding the progenitor systems and the explosion mechanism of normal Type Ia supernovae have thus far failed. With its unique capability to obtain continuous spectra through the near- and mid-infrared, JWST now offers completely new insights into Type Ia supernovae. In particular, observing them in the nebular phase allows us to directly see the central ejecta and thereby constrain the explosion mechanism. We aim to understand and quantify differences in the structure and composition of the central ejecta of various Type Ia supernova explosion models. We examined the currently most popular explosion scenarios using self-consistent multidimensional explosion simulations of delayed-detonation and pulsationally assisted, gravitationally confined delayed detonation Chandrasekhar-mass models and double-detonation sub-Chandrasekhar-mass and violent merger models. We find that the distribution of radioactive and stable nickel in the final ejecta, both observable in nebular spectra, are significantly different between different explosion scenarios. Therefore, comparing synthetic nebular spectra with JWST observations should allow us to distinguish between explosion models. We show that the explosion ejecta are inherently multidimensional for all models, and the Chandrasekhar-mass explosions simulated in spherical symmetry in particular lead to a fundamentally unphysical ejecta structure. Moreover, we show that radioactive and stable nickel cover a significant range of densities at a fixed velocity of the homologously expanding ejecta. Any radiation transfer postprocessing has to take these variations into account to obtain faithful synthetic observables; this will likely require multidimensional radiation transport simulations.

Nebular spectra from Type Ia supernova explosion models compared to JWST observations of SN 2021aefx

Context. Recent JWST observations of the Type Ia supernova (SN Ia) 2021aefx in the nebular phase have paved the way for late-time studies covering the full optical to mid-infrared (MIR) wavelength range, and with it the hope to better constrain SN Ia explosion mechanisms. Aims. We investigate whether public SN Ia models covering a broad range of progenitor scenarios and explosion mechanisms (Chandrasekhar-mass, or MCh, delayed detonations, pulsationally assisted gravitationally confined detonations, sub-MCh double detonations, and violent mergers) can reproduce the full optical-MIR spectrum of SN 2021aefx at ∼270 days post explosion. Methods. We consider spherically averaged 3D models available from the Heidelberg Supernova Model Archive with a 56Ni yield in the range 0.5–0.8 M⊙. We performed 1D steady-state non-local thermodynamic equilibrium simulations with the radiative-transfer code CMFGEN, and compared the predicted spectra to SN 2021aefx. Results. The models can explain the main features of SN 2021aefx over the full wavelength range. However, no single model, or mechanism, emerges as a preferred match, and the predicted spectra are similar to each other despite the very different explosion mechanisms. We discuss possible causes for the mismatch of the models, including ejecta asymmetries and ionisation effects. Our new calculations of the collisional strengths for Ni III have a major impact on the two prominent lines at 7.35 μm and 11.00 μm, and highlight the need for more accurate collisional data for forbidden transitions. Using updated atomic data, we identify a strong feature due to [Ca IV] 3.21 μm, attributed to [Ni I] in previous studies. We also provide a tentative identification of a forbidden line due to [Ne II] 12.81 μm, whose peaked profile indicates the presence of neon all the way to the innermost region of the ejecta, as predicted for instance in violent merger models. Contrary to previous claims, we show that the [Ar III] 8.99 μm line can be broader in sub-MCh models compared to near-MCh models. Last, the total luminosity in lines of Ni is found to correlate strongly with the stable nickel yield, although ionisation effects can bias the inferred abundance. Conclusions. Our models suggest that key physical ingredients are missing from either the explosion models, or the radiative-transfer post-processing, or both. Nonetheless, they also show the potential of the near- and MIR to uncover new spectroscopic diagnostics of SN Ia explosion mechanisms.

Early light curves for Type Ia supernova explosion models

Upcoming high-cadence transient survey programmes will produce a wealth of observational data for Type Ia supernovae. These data sets will contain numerous events detected very early in their evolution, shortly after explosion. Here, we present synthetic light curves, calculated with the radiation hydrodynamical approach Stella for a number of different explosion models, specifically focusing on these first few days after explosion. We show that overall the early light curve evolution is similar for most of the investigated models. Characteristic imprints are induced by radioactive material located close to the surface. However, these are very similar to the signatures expected from ejecta-CSM or ejecta-companion interaction. Apart from the pure deflagration explosion models, none of our synthetic light curves exhibit the commonly assumed power-law rise. We demonstrate that this can lead to substantial errors in the determination of the time of explosion. In summary, we illustrate with our calculations that even with very early data an identification of specific explosion scenarios is challenging, if only photometric observations are available.

Polarization spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multi-dimensional supernova explosion models. The approach utilises "virtual-packets" that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is vital for calculating synthetic spectropolarimetry (since the degree of polarisation is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealised test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarisation from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarisation spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

5.9-keV Mn K-shell X-ray luminosity from the decay of 55Fe in Type Ia supernova models

We show that the X-ray line flux of the Mn Kalpha line at 5.9 keV from the decay of 55Fe is a promising diagnostic to distinguish between Type Ia supernova (SN Ia) explosion models. Using radiation transport calculations, we compute the line flux for two 3D explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two white dwarfs. Both models are based on solar metallicity zero-age main sequence progenitors. Due to explosive nuclear burning at higher density, the delayed-detonation model synthesises 3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kalpha line flux of the delayed-detonation model exceeds that of the merger model by a factor of 4.5. Since in both models the 5.9 keV X-ray flux peaks five to six years after the explosion, a single measurement of the X-ray line emission at this time can place a constraint on the explosion physics that is complementary to those derived from earlier phase optical spectra or light curves. We perform detector simulations of current and future X-ray telescopes to investigate the possibilities of detecting the X-ray line at 5.9 keV. For the delayed-detonation scenario, a line detection is feasible with Chandra up to 3 Mpc for an exposure time of 10^6 s. We find that it should be possible with currently existing X-ray instruments (with exposure times 5x10^5 s) to detect both of our models at sufficiently high S/N to distinguish between them for hypothetical events within the Local Group. The prospects for detection will be better with future missions. For example, the proposed Athena/X-IFU instrument could detect our delayed-detonation model out to a distance of 5 Mpc. This would make it possible to study future events occurring during its operational life at distances comparable to those of the recent supernovae SN 2011fe (6.4 Mpc) and SN 2014J (3.5 Mpc).

Early X-ray emission from Type Ia supernovae originating from symbiotic progenitors or recurrent novae

One of the key observables for determining the progenitor nature of Type Ia supernovae is provided by their immediate circumstellar medium, which according to several models should be shaped by the progenitor binary system. So far, X-ray and radio observations indicate that the surroundings are very tenuous, producing severe upper-limits on the mass loss from winds of the progenitors. In this study, we perform numerical hydro-dynamical simulations of the interaction of the SN ejecta with circumstellar structures formed by possible mass outflows from the progenitor systems and we estimate numerically the expected numerical X-ray luminosity. We consider two kinds of circumstellar structures: a) A circumstellar medium formed by the donor star's stellar wind, in case of a symbiotic binary progenitor system; b) A circumstellar medium shaped by the interaction of the slow wind of the donor star with consecutive nova outbursts for the case of a symbiotic recurrent nova progenitor system. For the hydro-simulations we used well-known Type Ia supernova explosion models, as well as an approximation based on a power law model for the density structure of the outer ejecta. We confirm the strict upper limits on stellar wind mass loss, provided by simplified interpretations of X-ray upper limits of Type Ia supernovae. However, we show that supernova explosions going off in the cavities created by repeated nova explosions, provide a possible explanation for the lack of X-ray emission from supernovae originating from symbiotic binaries. Moreover, the velocity structure of circumstellar medium, shaped by a series of nova explosion matches well with the Na absorption features seen in absorption toward several Type Ia supernovae.

Probing thermonuclear supernova explosions with neutrinos

Aims: We present neutrino light curves and energy spectra for two representative type Ia supernova explosion models: a pure deflagration and a delayed detonation. Methods: We calculate the neutrino flux from $\beta$ processes using nuclear statistical equilibrium abundances convoluted with approximate neutrino spectra of the individual nuclei and the thermal neutrino spectrum (pair+plasma). Results: Although the two considered thermonuclear supernova explosion scenarios are expected to produce almost identical electromagnetic output, their neutrino signatures appear vastly different, which allow an unambiguous identification of the explosion mechanism: a pure deflagration produces a single peak in the neutrino light curve, while the addition of the second maximum characterizes a delayed-detonation. We identified the following main contributors to the neutrino signal: (1) weak electron neutrino emission from electron captures (in particular on the protons Co55 and Ni56) and numerous beta-active nuclei produced by the thermonuclear flame and/or detonation front, (2) electron antineutrinos from positron captures on neutrons, and (3) the thermal emission from pair annihilation. We estimate that a pure deflagration supernova explosion at a distance of 1 kpc would trigger about 14 events in the future 50 kt liquid scintillator detector and some 19 events in a 0.5 Mt water Cherenkov-type detector. Conclusions: While in contrast to core-collapse supernovae neutrinos carry only a very small fraction of the energy produced in the thermonuclear supernova explosion, the SN Ia neutrino signal provides information that allows us to unambiguously distinguish between different possible explosion scenarios. These studies will become feasible with the next generation of proposed neutrino observatories.

On the road to consistent Type Ia supernova explosion models

AbstractKeeping up with ever more detailed observations, Type Ia supernova (SN Ia) explosion models have seen a brisk development over the past years. The aim is to construct a self-consistent picture of the physical processes in order to gain the predictive power necessary to answer questions arising from the application of SNe Ia as cosmological distance indicators. We review recent developments in modeling these objects focusing on three-dimensional simulations.

Multi-spot ignition in type Ia supernova models

We present a systematic survey of the capabilities of type Ia supernova explosion models starting from a number of flame seeds distributed around the center of the white dwarf star. To this end we greatly improved the resolution of the numerical simulations in the initial stages. This novel numerical approach facilitates a detailed study of multi-spot ignition scenarios with up to hundreds of ignition sparks. Two-dimensional simulations are shown to be inappropriate to study the effects of initial flame configurations. Based on a set of three-dimensional models, we conclude that multi-spot ignition scenarios may improve type Ia supernova models towards better agreement with observations. The achievable effect reaches a maximum at a limited number of flame ignition kernels as shown by the numerical models and corroborated by a simple dimensional analysis.

Following multi-dimensional type Ia supernova explosion models to homologous expansion

The last years have witnessed a rapid development of three-dimensional models of Type Ia supernova explosions. Consequently, the next step is to evaluate these models under variation of the initial parameters and to compare them with observations. To calculate synthetic lightcurves and spectra from numerical models, it is mandatory to follow the evolution up to homologous expansion. We report on methods to achieve this in our current implementation of multi-dimensional Type Ia supernova explosion models. The novel scheme is thoroughly tested in two dimensions and a simple example of a three-dimensional simulation is presented. We discuss to what degree the assumption of homologous expansion is justified in these models.

Full-star type Ia supernova explosion models

We present full-star simulations of type Ia supernova explosions on the basis of the standard Chandrasekhar-mass deflagration model. Most simulations so far considered only one spatial octant and assumed mirror symmetry to the other octants. Two full-star models are evolved to homologous expansion and compared with previous single-octant simulations. Therefrom we analyze the effect of abolishing the artificial symmetry constraint on the evolution of the flame surface. It turns out that the development of asymmetries depends on the chosen initial flame configuration. Such asymmetries of the explosion process could possibly contribute to the observed polarization of some type Ia supernova spectra.

A model grid for the spectral analysis of X-ray emission in young Type Ia supernova remnants

We address a new set of models for the spectral analysis of the X-ray emission from young, ejecta-dominated Type Ia supernova remnants. These models are based on hydrodynamic simulations of the interaction between Type Ia supernova explosion models and the surrounding ambient medium, coupled to self-consistent ionization and electron heating calculations in the shocked supernova ejecta, and the generation of synthetic spectra with an appropriate spectral code. The details are provided elsewhere, but in this paper we concentrate on a specific class of Type Ia explosion models (delayed detonations), commenting on the differences that arise between their synthetic X-ray spectra under a variety of conditions.