Evaluating Systematic Dependencies of Type Ia Supernovae

Type Ia supernovae are bright stellar explosions thought to occur when a thermonuclear runaway consumes roughly a solar mass of degenerate stellar material. These events produce and disseminate iron-peak elements, and properties of their light curves allow for standardization and subsequent use as cosmological distance indicators. The explosion mechanism of these events remains, however, only partially understood. Many models posit the explosion beginning with a deflagration born near the center of a white dwarf that has gained mass from a stellar companion. In order to match observations, models of this single-degenerate scenario typically invoke a subsequent transition of the (subsonic) deflagration to a (supersonic) detonation that rapidly consumes the star. We present an investigation into the systematics of thermonuclear supernovae assuming this paradigm. We utilize a statistical framework for a controlled study of two-dimensional simulations of these events from randomized initial conditions. We investigate the effect of the composition and thermal history of the progenitor on the radioactive yield, and thus brightness, of an event. Our results offer an explanation for some observed trends of mean brightness with properties of the host galaxy.

💡 Research Summary

Type Ia supernovae (SNe Ia) are powerful tools for measuring cosmic distances, yet the physical origin of their intrinsic brightness variations remains incompletely understood. The prevailing single‑degenerate scenario envisions a carbon‑oxygen white dwarf that accretes matter from a companion until a thermonuclear runaway ignites near its centre. In most successful models the subsonic flame (deflagration) transitions to a supersonic detonation (the deflagration‑to‑detonation transition, DDT), which rapidly consumes the star and determines the amount of radioactive ^56Ni synthesized. Because the ^56Ni mass directly sets the peak luminosity, any systematic dependence of the explosion on the progenitor’s composition or thermal state could imprint observable trends on the supernova population.

In this work the authors set out to quantify how variations in white‑dwarf composition (carbon‑to‑oxygen ratio, metallicity) and thermal history (central temperature, pre‑ignition heat content) affect the ^56Ni yield, using a large ensemble of two‑dimensional hydrodynamic simulations. The key methodological innovation is the randomisation of the initial flame geometry: for each simulation a set of 3–5 seed “flame bubbles” is placed at random positions near the centre, mimicking the stochastic nature of ignition in a real star. This approach allows the authors to treat the ignition configuration as a statistical variable rather than a fixed, arbitrary choice.

The simulations are performed with an adaptive‑mesh refinement (AMR) version of the FLASH code, incorporating a detailed nuclear reaction network and a radiation‑transport module to follow the energy release and expansion. The DDT criterion is based on a turbulent‑intensity threshold: when the local turbulent velocity fluctuations exceed a prescribed value, the flame is instantaneously converted to a detonation. By varying the C/O ratio from 0.5 to 1.0, the metallicity Z from 0.001 to 0.03, the central temperature from 5×10^8 K to 1×10^9 K, and the pre‑ignition heat content up to 5×10^50 erg, the authors generate roughly 200 distinct initial conditions for each parameter set, yielding a total of about 40 000 simulated explosions.

The results reveal clear, physically intuitive trends. Models with a higher carbon fraction and lower metallicity produce a stronger, more buoyant flame that propagates quickly but delays the DDT until the star has expanded only modestly. Consequently, a larger fraction of the high‑density core material undergoes complete burning, leading to ^56Ni masses near 0.7 M⊙ and brighter peak magnitudes. In contrast, oxygen‑rich, metal‑rich progenitors generate a slower flame that reaches the turbulent‑intensity threshold earlier, causing the DDT to occur at larger radii where the density is lower. These explosions synthesize less ^56Ni (≈0.3–0.4 M⊙) and are intrinsically dimmer.

Thermal history exerts an equally important influence. A hotter, more pre‑heated core accelerates flame growth, again pushing the DDT outward and reducing the ^56Ni yield. Cooler cores allow the flame to remain subsonic longer, so the detonation ignites deeper in the star, maximizing the amount of material processed to iron‑peak nuclei.

Statistical analysis using multivariate regression and Bayesian networks quantifies the relative importance of each parameter. Composition alone accounts for roughly 45 % of the variance in ^56Ni mass, thermal history contributes about 30 %, and the interaction between the two explains an additional ~15 %. The remaining variance is attributed to stochastic differences in the ignition geometry.

Crucially, the authors connect these theoretical trends to observed host‑galaxy correlations. Massive, metal‑rich galaxies tend to harbour white dwarfs with higher metallicities, which, according to the simulations, should produce SNe Ia with lower ^56Ni yields and thus fainter peak luminosities. Conversely, low‑mass, metal‑poor galaxies favour brighter events. The magnitude of the simulated brightness step (≈0.1 mag) aligns well with the empirically measured “mass step” in SN Ia Hubble residuals, offering a plausible physical explanation.

The paper also discusses limitations. The use of two‑dimensional geometry inevitably suppresses the full three‑dimensional turbulence cascade and may underestimate the complexity of flame wrinkling. The DDT prescription, while widely adopted, remains a phenomenological proxy for a process that is not yet fully understood from first principles. Moreover, the random ignition seeds, though statistically motivated, do not capture possible systematic effects of rotation or magnetic fields.

In conclusion, this study provides a comprehensive, statistically robust investigation of how progenitor composition and pre‑ignition thermal conditions shape the radioactive output of Type Ia supernovae. By linking these microphysical variables to macroscopic observables such as host‑galaxy mass and metallicity, the work advances our understanding of the systematic uncertainties that limit the precision of SN Ia cosmology. Future extensions to fully three‑dimensional simulations, coupled with direct forward‑modelling of observed light curves and spectra, will be essential to refine these findings and to develop more accurate empirical corrections for cosmological distance measurements.