Numerical and Physical Challenges to Nebular Spectroscopy in Thermonuclear Supernovae

TL;DR

This paper discusses the complex physical and numerical challenges in simulating nebular spectra of Type Ia supernovae, emphasizing the importance of accurate physics modeling for interpreting JWST observations.

Contribution

It introduces a comprehensive non-LTE radiation code (HYDRA) for simulating nebular spectra and demonstrates its application to models matching JWST supernova observations.

Findings

Nebular spectra depend critically on the treatment of energy conversion processes.

Stimulated recombination limits over-ionization of high ions in supernova ejecta.

Hydrodynamical models can reproduce observed spectral line ratios and profiles.

Abstract

Thermodynamical explosions of White Dwarfs (WD)are one of the keys to high precision cosmology. Nebular spectra, namely mid-infrared (MIR) with JWST are an effective tool to probe for the multi-dimensional imprints of the explosion physics of WDs and their progenitor systems but also pose a challenge for simulations. What we observe as SNe Ia are low-energy photons, namely light curves, and spectra detected some days to years after the explosion. The light is emitted from a rapidly expanding envelope consisting of a low-density and low-temperature plasma with atomic population numbers far from thermodynamical equilibrium. SNe Ia are powered radioactive decays which produce hard X- and gamma-rays and MeV leptons which are converted within the ejecta to low-energy photons. We find that the optical and IR nebular spectra depend sensitively on the proper treatment of the physical conversion of high to low energies. The low-energy photons produced by forbidden line transitions originate from a mostly optically thin envelope. However, the UV is optically thick because of a quasi-continuum formed by allowed lines and bound-free transitions even several years after the explosion. We find that stimulated recombination limits the over-ionization of high ions with populations governed by the far UV. The requirements to simulate nebular spectra are well beyond both classical stellar atmospheres and nebulae. Using our full non-LTE HYDrodynamical RAdiation code (HYDRA) as a test-bed, the sensitivity on the physics on synthetic spectra are demonstrated using observations as a benchmark. At some examples, we establish the power of high-precision nebular spectroscopy as quantitative tool. Centrally ignited, off-center delayed-detonation near Chandrasekhar-mass models can reproduce line-ratios and line profiles of Branch-normal and underluminous SNe Ia observed with JWST.