Ivo R. Seitenzahl

Double-detonation sub-Chandrasekhar supernovae: can minimum helium shell masses detonate the core?

The explosion of sub-Chandrasekhar mass white dwarfs via the double detonation scenario is a potential explanation for type Ia supernovae. In this scenario, a surface detonation in a heli um layer initiates a detonation in the underlying carbon/oxygen core leading to an explosion. For a given core mass, a lower bound has been determined on the mass of the helium shell required for dynamical burning during a helium flash, which is a necessary prerequis ite for detonation. For a range of core and corresponding minimum helium shell masses, we investigate whether an assumed surface helium detonation is capable of triggering a subsequent detonation in the core even for this limiting case. We carried out hydrodynamic simulations on a co-expanding Eulerian grid in two dimensions assuming rotational symmetry. The detonations are propagated using the level-set approach and a simplified scheme for n uclear reactions that has been calibrated with a large nuclear network. The same network is used to determine detailed nucleosynthetic abundances in a post-processing step. Based on approximate detonation initiation criteria in the literature, we find th at secondary core

Three-dimensional delayed-detonation models with nucleosynthesis for Type Ia supernovae

We present results for a suite of fourteen 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 I a simulations with detailed isotopic yield information. As such, it may serve as a database for Chandrasekhar-mass delayeddetonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ a physically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ t he deflagration to detonation transition (DDT) 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 central density of 2.9× 10 9 g cm −3 , plus in addition one high central density (5.5× 10 9 g cm −3 ) and one low central density (1.0× 10 9 g cm −3 ) rendition of the 100 ignition kernel configuration. For each simulatio n we determined detailed nucleosynthetic yields by post-processing 10 6 tracer particles with a 384 nuclide reaction network. All delayed detonation models result in explosions unbinding the white dwarf, producing a range of 56 Ni masses from 0.32 to 1.11 M⊙. As a general trend, the models predict that the stable neutron-rich iron group isotopes are not found at the lowest velocities, but rather at intermediate velocities (∼3, 000− 10, 000 km s −1 ) in a shell surrounding a 56 Ni-rich core. The models further predict relatively low velocity oxygen and carbon, with typical minimum velocities around 4, 000 and 10, 000 km s −1 , respectively.

DETONATIONS IN SUB-CHANDRASEKHAR-MASS C+O WHITE DWARFS

Explosions of sub-Chandrasekhar-mass white dwarfs (WDs) are one alternative to the standard Chandrasekhar-mass model of Type Ia supernovae (SNe Ia). They are interesting since binary systems with sub-Chandrasekhar-mass primary WDs should be common and this scenario would suggest a simple physical parameter which determines the explosion brightness, namely the mass of the exploding WD. Here we perform one-dimensional hydrodynamical simulations, associated post-processing nucleosynthesis, and multi-wavelength radiation transport calculations for pure detonations of carbon-oxygen WDs. The light curves and spectra we obtain from these simulations are in good agreement with observed properties of SNe Ia. In particular, for WD masses from 0.97 to 1.15 M ☉ we obtain 56Ni masses between 0.3 and 0.8 M ☉, sufficient to capture almost the complete range of SN Ia brightnesses. Our optical light curve rise times, peak colors, and decline timescales display trends which are generally consistent with observed characteristics although the range of B-band decline timescales displayed by our current set of models is somewhat too narrow. In agreement with observations, the maximum light spectra of the models show clear features associated with intermediate-mass elements and reproduce the sense of the observed correlation between explosion luminosity and the ratio of the Si II lines at λ6355 and λ5972. We therefore suggest that sub-Chandrasekhar-mass explosions are a viable model for SNe Ia for any binary evolution scenario leading to explosions in which the optical display is dominated by the material produced in a detonation of the primary WD.

Constraining type Ia supernova models: SN 2011fe as A test case

The nearby supernova SN 2011fe can be observed in unprecedented detail. Therefore, it is an important test case for Type Ia supernova (SN Ia) models, which may bring us closer to understanding the physical nature of these objects. Here, we explore how available and expected future observations of SN 2011fe can be used to constrain SN Ia explosion scenarios. We base our discussion on three-dimensional simulations of a delayed detonation in a Chandrasekhar-mass white dwarf and of a violent merger of two white dwarfs (WDs)—realizations of explosion models appropriate for two of the most widely discussed progenitor channels that may give rise to SNe Ia. Although both models have their shortcomings in reproducing details of the early and near-maximum spectra of SN 2011fe obtained by the Nearby Supernova Factory (SNfactory), the overall match with the observations is reasonable. The level of agreement is slightly better for the merger, in particular around maximum, but a clear preference for one model over the other is still not justified. Observations at late epochs, however, hold promise for discriminating the explosion scenarios in a straightforward way, as a nucleosynthesis effect leads to differences in the 55Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect.

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

cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae. Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ejecta. Pure turbulent deagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative transfer calculations for a 3D full-star deagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a proto-typical 2002cx-like supernova), not only in the optical, but also in the near- infrared. For that purpose we present, for the rst time, ve near-infrared spectra of SN 2005hk from 0:2 to 26:6 days 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.03 M { consisting 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 inuence the late-time observables of 2002cx-like SNe.

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

We investigate whether pure deflagration models of Chandrasekhar-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-nuc less than or similar to 1.1 x 10(51) 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-nuc similar to 0.5 x 10(51) 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.

SPONTANEOUS INITIATION OF DETONATIONS IN WHITE DWARF ENVIRONMENTS: DETERMINATION OF CRITICAL SIZES

Some explosion models for Type Ia supernovae (SNe Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate / material of a white dwarf (WD). The length scales pertinent to the initiation of the detonation are notoriously unresolved in multidimensional stellar simulations, prompting the use of results of one-dimensional simulations at higher resolution, such as those performed for this work, as guidelines for deciding whether or not conditions reached in the higher dimensional full star simulations successfully would lead to the onset of a detonation. Spontaneous initiation relies on the existence of a suitable gradient in self-ignition (induction) times of the fuel, which we set up with a spatially localized nonuniformity of temperature?a hot spot. We determine the critical (smallest) sizes of such hot spots that still marginally result in a detonation in WD matter by integrating the reactive Euler equations with the hydrodynamics code FLASH. We quantify the dependences of the critical sizes of such hot spots on composition, background temperature, peak temperature, geometry, and functional form of the temperature disturbance, many of which were hitherto largely unexplored in the literature. We discuss the implications of our results in the context of modeling of SNe Ia.

Extensive HST ultraviolet spectra and multiwavelength observations of SN 2014J in M82 indicate reddening and circumstellar scattering by typical dust

SN 2014J in M82 is the closest detected Type Ia supernova (SN Ia) in at least 28 years and perhaps in 410 years. Despite its small distance of 3.3 Mpc, SN 2014J is surprisingly faint, peaking at V = 10.6 mag, and assuming a typical SN Ia luminosity, we infer an observed visual extinction of AV = 2.0 ± 0.1 mag. But this picture, with RV = 1.6 ± 0.2, is too simple to account for all observations. We combine 10 epochs (spanning a month) of HST/STIS ultraviolet through near-infrared spectroscopy with HST/WFC3, KAIT, and FanCam photometry from the optical to the infrared and 9 epochs of high-resolution TRES spectroscopy to investigate the sources of extinction and reddening for SN 2014J. We argue that the wide range of observed properties for SN 2014J is caused by a combination of dust reddening, likely originating in the interstellar medium of M82, and scattering off circumstellar material. For this model, roughly half of the extinction is caused by reddening from typical dust (E(B V ) = 0.45 mag and RV = 2.6) and roughly half by scattering off LMC-like dust in the circumstellar environment of SN 2014J.

Study of the Detonation Phase in the Gravitationally Confined Detonation Model of Type Ia Supernovae

We study the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia) through the detonation phase and into homologous expansion. In the GCD model, a detonation is triggered by the surface flow due to single-point, off-center flame ignition in carbon-oxygen white dwarfs (WDs). The simulations are unique in terms of the degree to which nonidealized physics is used to treat the reactive flow, including weak reaction rates and a time-dependent treatment of material in nuclear statistical equilibrium (NSE). Careful attention is paid to accurately calculating the final composition of material which is burned to NSE and frozen out in the rapid expansion following the passage of a detonation wave over the high-density core of the WD; and an efficient method for nucleosynthesis postprocessing is developed which obviates the need for costly network calculations along tracer particle thermodynamic trajectories. Observational diagnostics are presented for the explosion models, including abundance stratifications and integrated yields. We find that for all of the ignition conditions studied here a self-regulating process comprised of neutronization and stellar expansion results in final 56Ni masses of ~1.1?M ?. But, more energetic models result in larger total NSE and stable Fe-peak yields. The total yield of intermediate mass elements is ~0.1?M ? and the explosion energies are all around 1.5 ? 1051 erg. The explosion models are briefly compared to the inferred properties of recent SN Ia observations. The potential for surface detonation models to produce lower-luminosity (lower 56Ni mass) SNe is discussed.

On the brightness distribution of type Ia supernovae from violent white dwarf mergers

We investigate the brightness distribution expected for th ermonuclear explosions that might result from the ignition of a detonation during the violent merger of white dwarf (WD) binaries. Violent WD mergers are a subclass of the canonical double degenerate scenario where two carbon-oxygen (CO) WDs merge when the larger WD fills its Roche-lobe. Determining their brightness distribution is critical for evaluating w hether such an explosion model could be responsible for a significant fraction of the observed pop ulation of Type Ia supernovae (SNe Ia). We argue that the brightness of an explosion realised via the violent merger model is mainly determined by the mass of 56 Ni produced in the detonation of the primary CO WD. To quantify this link, we use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass (mWD) and expected peak bolometric brightness (Mbol). We use this mWD-Mbol relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We also investigate the sensitivity of our results to assumptions abo ut the conditions required to realise a detonation during violent mergers of WDs. We find a striking s imilarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a Mbol distribution that roughly covers the range and matches the shape of the one observed for SNe Ia. However, this agreement hinges on a particular phase of mass accretion during binary evolution: the primary WD gains �0.15 0.35 M⊙ from a slightly-evolved helium star companion. In our standard binary evolution model, such an accretion phase is predicted to occur for about 43% of all binary systems that ultimately give rise to binary CO WD mergers. We also find that with high probability, violent WD merge rs involving the most massive primaries ( �1.3M⊙, which should produce bright SNe) have delay times �500 Myr.