3D deflagration simulations leaving bound remnants:A model for 2002cx-like Type Ia supernovae

2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae (SNe). Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ejecta. Pure turbulent deflagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative-transfer calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a prototypical 2002cx-like supernova), not only in the optical, but also in the near-infrared. For that purpose we present, for the first time, five near-infrared spectra of SN 2005hk from -0.2 to 26.6 d with respect to B-band maximum. Since our model burns only small parts of the initial white dwarf, it fails to completely unbind the white dwarf and leaves behind a bound remnant of ~1.03Mconsisting mainly of unburned carbon and oxygen, but also enriched by some amount of intermediate-mass and iron-group elements from the explosion products that fall back on the remnant.We discuss possibilities for detecting this bound remnant and how it might influence the late-time observables of 2002cx-like SNe. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

Three-dimensional pure deflagration models with nucleosynthesis and synthetic observables for type ia supernovae

We investigate whether pure deflagration models ofChandrasekhar-mass carbon-oxygen white dwarf stars can account for one or more subclass of the observed population of Type Ia supernova (SN Ia) explosions. We compute a set of 3D full-star hydrodynamic explosion models, in which the deflagration strength is parametrized using the multispot ignition approach. For each model, we calculate detailed nucleosynthesis yields in a post-processing step with a 384 nuclide nuclear network. We also compute synthetic observables with our 3D Monte Carlo radiative transfer code for comparison with observations. For weak and intermediate deflagration strengths (energy release E {less-than or approximate} 1.1 × 10 erg), we find that the explosion leaves behind a bound remnant enriched with 3 to 10 per cent (by mass) of deflagration ashes. However, we do not obtain the large kick velocities recently reported in the literature. We find that weak deflagrations with E ~ 0.5 × 10 erg fit well both the light curves and spectra of 2002cx-like SNe Ia, and models with even lower explosion energies could explain some of the fainter members of this subclass. By comparing our synthetic observables with the properties of SNe Ia, we can exclude the brightest, most vigorously ignited models as candidates for any observed class of SN Ia: their B-V colours deviate significantly from both normal and 2002cx-like SNe Ia and they are too bright to be candidates for other subclasses.

The type Iax supernova, SN 2015H: A white dwarf deflagration candidate

We present results based on observations of SN 2015H which belongs to the small group of objects similar to SN 2002cx, otherwise known as type Iax supernovae. The availability of deep pre-explosion imaging allowed us to place tight constraints on the explosion epoch. Our observational campaign began approximately one day post-explosion, and extended over a period of about 150 days post maximum light, making it one of the best observed objects of this class to date. We find a peak magnitude of M_r = -17.27± 0.07, and a (Δm₁₅)_r = 0.69 ± 0.04. Comparing our observations to synthetic spectra generated from simulations of deflagrations of Chandrasekhar mass carbon-oxygen white dwarfs, we find reasonable agreement with models of weak deflagrations that result in the ejection of ∼0.2 M_⊙ of material containing ∼0.07 M_⊙ of ⁵⁶Ni. The model light curve however, evolves more rapidly than observations, suggesting that a higher ejecta mass is to be favoured. Nevertheless, empirical modelling of the pseudo-bolometric light curve suggests that ≲ 0.6 M_⊙ of material was ejected, implying that the white dwarf is not completely disrupted, and that a bound remnant is a likely outcome.

Signatures of delayed detonation, asymmetry, and electron capture in the mid-infrared spectra of supernovae 2003hv and 2005df

We present mid-infrared (5.2-15.2 mu m) spectra of the Type Ia supernovae (SNe Ia) 2003hv and 2005df observed with the Spitzer Space Telescope. These are the first observed mid-infrared spectra of thermonuclear supernovae, and show strong emission from fine-structure lines of Ni, Co, S, and Ar. The detection of Ni emission in SN 2005df 135 days after the explosion provides direct observational evidence of high-density nuclear burning forming a significant amount of stable Ni in a SN Ia. The SN 2005df Ar lines also exhibit a two-pronged emission profile, implying that the Ar emission deviates significantly from spherical symmetry. The spectrum of SN 2003hv also shows signs of asymmetry, exhibiting blueshifted [Co (III)], which matches the blueshift of [Fe (II)] lines in nearly coeval near-infrared spectra. Finally, local thermodynamic equilibrium abundance estimates for the yield of radioactive Ni-56 give M-56Ni approximate to 0.5 M-circle dot, for SN 2003hv, but only M-56Ni approximate to 0.13-0.22 M-circle dot for the apparently subluminous SN 2005df, supporting the notion that the luminosity of SNe Ia is primarily a function of the radioactive 56Ni yield. The observed emission-line profiles in the SN 2005df spectrum indicate a chemically stratified ejecta structure, which matches the predictions of delayed detonation (DD) models, but is entirely incompatible with current three-dimensional deflagration models. Furthermore, the degree that this layering persists to the innermost regions of the supernova is difficult to explain even in a DD scenario, where the innermost ejecta are still the product of deflagration burning. Thus, while these results are roughly consistent with a delayed detonation, it is clear that a key piece of physics is still missing from our understanding of the earliest phases of SN Ia explosions.

High-velocity features: A ubiquitous property of Type Ia supernovae

Evidence of high-velocity features (HVFs) such as those seen in the near-maximum spectra of some Type Ia supernovae (SNe Ia; e. g., SN 2000cx) has been searched for in the available SN Ia spectra observed earlier than 1 week before B maximum. Recent observational efforts have doubled the number of SNe Ia with very early spectra. Remarkably, all SNe Ia with early data ( seven in our Research Training Network sample and 10 from other programs) show signs of such features, to a greater or lesser degree, in Ca II IR and some also in the Si II lambda 6355 line. HVFs may be interpreted as abundance or density enhancements. Abundance enhancements would imply an outer region dominated by Si and Ca. Density enhancements may result from the sweeping up of circumstellar material (CSM) by the highest velocity SN ejecta. In this scenario, the high incidence of HVFs suggests that a thick disk and/or a high-density companion wind surrounds the exploding white dwarf, as may be the case in single degenerate systems. Large-scale angular fluctuations in the radial density and abundance distribution may also be responsible: this could originate in the explosion and would suggest a deflagration as the more likely explosion mechanism. CSM interaction and surface fluctuations may coexist, possibly leaving different signatures on the spectrum. In some SNe, the HVFs are narrowly confined in velocity, suggesting the ejection of blobs of burned material.

Three-dimensional delayed-detonation models with nucleosynthesis for type ia supernovae

We present results for a suite of 14 three-dimensional, high-resolution hydrodynamical simulations of delayed-detonation models of Type Ia supernova (SN Ia) explosions. This model suite comprises the first set of three-dimensional SN Ia simulations with detailed isotopic yield information. As such, it may serve as a data base for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ aphysically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration-to-detonation transition probabilities. To obtain different strengths of the deflagration phase and thereby different degrees of pre-expansion, we have chosen a sequence of initial models with 1, 3, 5, 10, 20, 40, 100, 150, 200, 300 and 1600 (two different realizations) ignition kernels in a hydrostatic white dwarf with a central density of 2.9 × 10 g cm, as well as one high central density (5.5 × 10 g cm) and one low central density (1.0 × 10 g cm) rendition of the 100 ignition kernel configuration. For each simulation, we determined detailed nucleosynthetic yields by postprocessing10 tracer particles with a 384 nuclide reaction network. All delayed-detonation models result in explosions unbinding thewhite dwarf, producing a range of 56Ni masses from 0.32 to 1.11M. As a general trend, the models predict that the stableneutron-rich iron-group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3000×10 000 km s) in a shell surrounding a Ni-rich core. The models further predict relatively low-velocity oxygen and carbon, with typical minimum velocities around 4000 and 10 000 km s, respectively. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Synthetic light curves and spectra for three-dimensional delayed-detonation models of type ia supernovae

In a companion paper, Seitenzahl et al. have presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar-mass white dwarfs (WDs). Here,we present multidimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (which is controlled by the ignition configuration in our models) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to the observed properties of normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time-scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are systematically too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width- luminosity relation (WLR) of SNe Ia: all our models show rather similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.

Deflagrations in hybrid CONe white dwarfs: a route to explain the faint Type Iax supernova 2008ha

Stellar evolution models predict the existence of hybrid white dwarfs (WDs) with a carbon-oxygen core surrounded by an oxygen-neon mantle. Being born with masses similar to 1.1 M-aS (TM), hybrid WDs in a binary system may easily approach the Chandrasekhar mass (M-Ch) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-M-Ch hybrid WD under the assumption that nuclear burning only occurs in carbon-rich material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthesis post-processing calculations, we find that only 0.014 M-aS (TM) of material is ejected while the remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, O, C, Si and S. We also calculate synthetic observables for our model and find reasonable agreement with the faint Type Iax SN 2008ha. This shows for the first time that deflagrations in near-M-Ch WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-M-Ch bound remnant opens the possibility for recurrent explosions or a subsequent accretion-induced collapse in faint Type Iax SNe, if further accretion episodes occur. From binary population synthesis calculations, we find the rate of hybrid WDs approaching M-Ch to be of the order of 1 per cent of the Galactic SN Ia rate.

Polarization spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multidimensional supernova explosion models. The approach utilizes '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 not only vital for calculating synthetic spectropolarimetry (since the degree of polarization is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealized 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 polarization 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 polarization spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

Three-dimensional simulations of gravitationally confined detonations compared to observations of SN 1991T

The gravitationally confined detonation (GCD) model has been proposed as a possible explosion mechanism for Type Ia supernovae in the single-degenerate evolution channel. It starts with ignition of a deflagration in a single off-centre bubble in a near-Chandrasekhar-mass white dwarf. Driven by buoyancy, the deflagration flame rises in a narrow cone towards the surface. For the most part, the main component of the flow of the expanding ashes remains radial, but upon reaching the outer, low-pressure layers of the white dwarf, an additional lateral component develops. This causes the deflagration ashes to converge again at the opposite side, where the compression heats fuel and a detonation may be launched. We first performed five three-dimensional hydrodynamic simulations of the deflagration phase in 1.4 M⊙ carbon/oxygen white dwarfs at intermediate-resolution (2563computational zones). We confirm that the closer the initial deflagration is ignited to the centre, the slower the buoyant rise and the longer the deflagration ashes takes to break out and close in on the opposite pole to collide. To test the GCD explosion model, we then performed a high-resolution (5123 computational zones) simulation for a model with an ignition spot offset near the upper limit of what is still justifiable, 200 km. This high-resolution simulation met our deliberately optimistic detonation criteria, and we initiated a detonation. The detonation burned through the white dwarf and led to its complete disruption. For this model, we determined detailed nucleosynthetic yields by post-processing 106 tracer particles with a 384 nuclide reaction network, and we present multi-band light curves and time-dependent optical spectra. We find that our synthetic observables show a prominent viewing-angle sensitivity in ultraviolet and blue wavelength bands, which contradicts observed SNe Ia. The strong dependence on the viewing angle is caused by the asymmetric distribution of the deflagration ashes in the outer ejecta layers. Finally, we compared our model to SN 1991T. The overall flux level of the model is slightly too low, and the model predicts pre-maximum light spectral features due to Ca, S, and Si that are too strong. Furthermore, the model chemical abundance stratification qualitatively disagrees with recent abundance tomography results in two key areas: our model lacks low-velocity stable Fe and instead has copious amounts of high-velocity 56Ni and stable Fe. We therefore do not find good agreement of the model with SN 1991T.

LES of flame acceleration and DDT in hydrogen-air mixture using artificially thickened flame approach and detailed chemical kinetics

A large eddy simulation is performed to study the deflagration to detonation transition phenomenon in an obstructed channel containing premixed stoichiometric hydrogen–air mixture. Two-dimensional filtered reactive Navier–Stokes equations are solved utilizing the artificially thickened flame approach (ATF) for modeling sub-grid scale combustion. To include the effect of induction time, a 27-step detailed mechanism is utilized along with an in situ adaptive tabulation (ISAT) method to reduce the computational cost due to the detailed chemistry. The results show that in the slow flame propagation regime, the flame–vortex interaction and the resulting flame folding and wrinkling are the main mechanisms for the increase of the flame surface and consequently acceleration of the flame. Furthermore, at high speed, the major mechanisms responsible for flame propagation are repeated reflected shock–flame interactions and the resulting baroclinic vorticity. These interactions intensify the rate of heat release and maintain the turbulence and flame speed at high level. During the flame acceleration, it is seen that the turbulent flame enters the ‘thickened reaction zones’ regime. Therefore, it is necessary to utilize the chemistry based combustion model with detailed chemical kinetics to properly capture the salient features of the fast deflagration propagation.