Source author record

F. K. Roepke

F. K. Roepke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.HE astro-ph astro-ph.CO astro-ph.GA astro-ph.IM nucl-th physics.flu-dyn hep-ph nucl-ex physics.comp-ph

Catalog footprint

What is connected

55works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

NLTE spectral modelling for a carbon-oxygen and helium white-dwarf merger as a Ca-rich transient candidate

We carry out NLTE (non local thermodynamic equilibrium) radiative transfer simulations to determine whether explosion during the merger of a carbon-oxygen (CO) white dwarf (WD) with a helium (He) WD can reproduce the characteristic Ca II/[Ca II] and He I lines observed in Ca-rich transients. Our study is based on a 1D representation of a hydrodynamic simulation of a 0.6 $M_{\odot}$ CO + 0.4 $M_{\odot}$ He WD merger. We calculate both photospheric and nebular-phase spectra including treatment for non-thermal electrons, as is required for accurate modelling of He I and [Ca II]. Consistent with Ca-rich transients, our simulation predicts a nebular spectrum dominated by emission from [Ca II] 7291, 7324 angstrom and the Ca II near-infrared (NIR) triplet. The photospheric-phase synthetic spectrum also exhibits a strong Ca II NIR triplet, prominent optical absorption due to He I 5876 angstrom and He I 10830 angstrom in the NIR, as is commonly observed for Ca-rich transients. Overall, our results therefore suggest that CO+He WD mergers are a promising channel for Ca-rich transients. However, the current simulation overpredicts some He I features, in particular both He I 6678 and 7065 angstrom and shows a significant contribution from Ti II, which results in a spectral energy distribution that is substantially redder than most Ca-rich transients at peak. Additionally the Ca II nebular emission features are too broad. Future work should investigate if these discrepancies can be resolved by considering full 3D models and exploring a range of CO+He WD binary configurations.

preprint2022arXiv

Bipolar planetary nebulae from common envelope evolution of binary stars

Asymmetric shapes and evidence for binary central stars suggest a common-envelope origin for many bipolar planetary nebulae. The bipolar components of the nebulae are observed to expand faster than the rest and the more slowly expanding material has been associated with the bulk of the envelope ejected during the common-envelope phase of a stellar binary system. Common-envelope evolution in general remains one of the biggest uncertainties in binary star evolution and the origin of the fast outflow has not been explained satisfactorily. We perform three-dimensional magnetohydrodynamic simulations of common-envelope interaction with the moving-mesh code AREPO. Starting from the plunge-in of the companion into the envelope of an asymptotic giant branch star and covering hundreds of orbits of the binary star system, we are able to follow the evolution to complete envelope ejection. We find that magnetic fields are strongly amplified in two consecutive episodes. First, when the companion spirals in the envelope and, second, when it forms a contact binary with the core of the former giant star. In the second episode, a magnetically-driven, high-velocity outflow of gas is launched self-consistently in our simulations. The outflow is bipolar and the gas is additionally collimated by the ejected common envelope. The resulting structure reproduces typical morphologies and velocities observed in young planetary nebulae. We propose that the magnetic driving mechanism is a universal consequence of common envelope interaction responsible for a substantial fraction of observed planetary nebulae. Such a mechanism likely also exists in the common-envelope phase of other binary stars that lead to the formation of Type Ia supernovae, X-ray binaries and gravitational-wave merger events.

preprint2022arXiv

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

preprint2022arXiv

Models of pulsationally assisted gravitationally confined detonations with different ignition conditions

Over the past decades, many explosion scenarios for Type Ia supernovae have been proposed and investigated including various combinations of deflagrations and detonations in white dwarfs of different masses up to the Chandrasekhar mass. One of these is the gravitationally confined detonation model. In this case a weak deflagration burns to the surface, wraps around the bound core, and collides at the antipode. A subsequent detonation is then initiated in the collision area. Since the parameter space for this scenario, that is, varying central densities and ignition geometries, has not been studied in detail, we used pure deflagration models of a previous parameter study dedicated to Type Iax supernovae as initial models to investigate the gravitationally confined detonation scenario. We aim to judge whether this channel can account for one of the many subgroups of Type Ia supernovae, or even normal events. To this end, we employed a comprehensive pipeline for three-dimensional Type Ia supernova modeling that consists of hydrodynamic explosion simulations, nuclear network calculations, and radiative transfer. The observables extracted from the radiative transfer are then compared to observed light curves and spectra. The study produces a wide range in masses of synthesized 56 Ni ranging from 0.257 to 1.057 $M_\odot$ , and, thus, can potentially account for subluminous as well as overluminous Type Ia supernovae in terms of brightness. However, a rough agreement with observed light curves and spectra can only be found for 91T-like objects. Although several discrepancies remain, we conclude that the gravitationally confined detonation model cannot be ruled out as a mechanism to produce 91T-like objects. However, the models do not provide a good explanation for either normal Type Ia supernovae or Type Iax supernovae.

preprint2021arXiv

Type Iax supernovae from deflagrations in Chandrasekhar mass white dwarfs

Due to the increasing number of observations Type Ia supernovae are nowadays regarded as a heterogeneous class of objects consisting of several subclasses. One of the largest of these is the class of Type Iax supernovae (SNe Iax) which have been suggested to originate from pure deflagrations in CO Chandrasekhar-mass white dwarfs (WDs). Although a few deflagration studies have been carried out, the full diversity of the class is not captured yet. We therefore present a parameter study of single-spot ignited deflagrations with varying ignition locations, central densities, metallicities and compositions. We also explore a rigidly rotating progenitor and carry out 3D hydrodynamic simulations, nuclear network calculations and radiative transfer. The new models extend the range in brightness covered by previous studies to the lower end. Our explosions produce $^{56}$Ni masses from $5.8 \times 10^{-3}$ to $9.2 \times 10^{-2}\,M_\odot$. In spite of the wide exploration of the parameter space the main characteristics of the models are primarily driven by the mass of $^{56}$Ni. Secondary parameters have too little impact to explain the observed trend among faint SNe~Iax. We report kick velocities of the bound explosion remnants from $6.9$ to $369.8\,$km$\,s^{-1}$. The wide exploration of the parameter space and viewing-angle effects in the radiative transfer lead to a significant spread in the synthetic observables. The trends towards the faint end of the class are, however, not reproduced. This motivates a quantification of the systematic uncertainties in the modeling procedure and the influence of the $^{56}$Ni-rich bound remnant. While the pure deflagration scenario remains a favorable explanation for bright and intermediate luminosity SNe~Iax, the possibility that SNe~Iax do not consist of a single explosion scenario needs to be considered.

preprint2020arXiv

Common-envelope evolution with an asymptotic giant branch star

Common-envelope phases are decisive for the evolution of many binary systems. Of particular interest are cases with asymptotic giant branch (AGB) primary stars, because they are thought to be progenitors of various astrophysical transients. In three-dimensional hydrodynamic simulations with the moving-mesh code AREPO, we study the common-envelope evolution of a $1.0\,M_{\odot}$ early-AGB star with companions of different masses. Although the stellar envelope of the AGB star is less tightly bound than that of a red giant, we find that the release of orbital energy of the core binary is insufficient to eject more than about twenty percent of the envelope mass. Ionization energy released in the expanding envelope, however, can lead to complete envelope ejection. Because recombination proceeds largely at high optical depths in our simulations, it is likely that this effect indeed plays a significant role in the considered systems. The efficiency of mass loss and the final orbital separation of the core binary system depend on the mass ratio between the companion and the primary star. Our results suggest a linear relation between the ratio of final to initial orbital separation and this parameter.

preprint2020arXiv

Fully compressible simulations of waves and core convection in main-sequence stars

Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional hydrodynamics code SLH to simulate stellar oscillations. Methods. We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic $3\mathrm{M}_\odot$ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results. Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with observed spectra of massive stars. Conclusion. The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows the coupling of gravity and pressure waves to be studied too.

preprint2020arXiv

Long-term evolution of a magnetic massive merger product

About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $τ$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars.

preprint2020arXiv

Nucleosynthesis imprints from different Type Ia Supernova explosion scenarios and implications for galactic chemical evolution

We analyze the nucleosynthesis yields of various Type Ia supernova explosion simulations including pure detonations in sub- Chandrasekhar mass white dwarfs, double detonations and pure helium detonations of sub-Chandrasekhar mass white dwarfs with an accreted helium envelope, a violent merger model of two white dwarfs and deflagrations as well as delayed detonations in Chandrasekhar mass white dwarfs. We focus on the iron peak elements Mn, Zn and Cu. To this end, we also briefly review the different burning regimes and production sites of these elements as well as the results of abundance measurements and several galactic chemical evolution studies. We find that super-solar values of [Mn/Fe] are not restricted to Chandrasekhar mass explosion models. Scenarios including a helium detonation can significantly contribute to the production of Mn, in particular the models proposed for calcium-rich transients. Although Type Ia supernovae are often not accounted for as production sites of Zn and Cu, our models involving helium shell detonations can produce these elements in super-solar ratios relative to Fe. Our results suggest a re-consideration of Type Ia supernova yields in galactic chemical evolution models. A detailed comparison with observations can provide new insight into the progenitor and explosion channels of these events.

preprint2020arXiv

Photometric detection of internal gravity waves in upper main-sequence stars. II. Combined TESS photometry and high-resolution spectroscopy

Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detection of IGWs across a large number of O and early-B stars in the Hertzsprung-Russell diagram, and explain the ubiquitous detection of stochastic variability in the photospheres of massive stars. Methods. We combined high-precision time-series photometry from the NASA Transiting Exoplanet Survey Satellite with high-resolution ground-based spectroscopy of 70 stars with spectral types O and B to probe the relationship between the photometric signatures of IGWs and parameters such as spectroscopic mass, luminosity, and macroturbulence. Results. A relationship is found between the location of a star in the spectroscopic Hertzsprung-Russell diagram and the amplitudes and frequencies of stochastic photometric variability in the light curves of massive stars. Furthermore, the properties of the stochastic variability are statistically correlated with macroturbulent velocity broadening in the spectral lines of massive stars. Conclusions. The common ensemble morphology for the stochastic low-frequency variability detected in space photometry and its relationship to macroturbulence is strong evidence for IGWs in massive stars, since these types of waves are unique in providing the dominant tangential velocity field required to explain the observed spectroscopy.

preprint2020arXiv

White dwarf deflagrations for Type Iax supernovae: Polarisation signatures from the explosion and companion interaction

Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.

preprint2019arXiv

Discovery of an Exceptionally Strong $β$-Decay Transition of $^{20}$F and Implications for the Fate of Intermediate-Mass Stars

A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fully constrain the rate. By measuring the transition, we have established that its strength is exceptionally large and enhances the capture rate by several orders of magnitude. This has a decisive impact on the evolution of the core, increasing the likelihood that the star is (partially) disrupted by a thermonuclear explosion rather than collapsing to form a neutron star. Importantly, our measurement resolves the last remaining nuclear physics uncertainty in the final evolution of degenerate oxygen-neon stellar cores, allowing future studies to address the critical role of convection, which at present is poorly understood.

preprint2019arXiv

Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae

We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.

preprint2018arXiv

The interaction of core-collapse supernova ejecta with a stellar companion

The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. We find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, is in general not be dramatic.

preprint2016arXiv

Measurements of the Soft Gamma-ray Emission from SN2014J with Suzaku

The hard X-ray detector (HXD) onboard {\it Suzaku} measured soft $γ$-rays from the Type Ia supernova SN2014J at $77\pm2$ days after the explosion. Although the confidence level of the signal is about 90\% (i.e., $2 σ$), the $3 σ$ upper limit has been derived at $< 2.2 \times10^{-4}$ ph s$^{-1}$ cm$^{-2}$ in the 170 -- 250 keV band as the first independent measurement of soft $γ$-rays with an instrument other than {\it INTEGRAL}. For this analysis, we have examined the reproducibility of the NXB model of HXD/GSO using blank sky data. We find that the residual count rate in the 90 -- 500 keV band is distributed around an average of 0.19\% with a standard deviation of 0.42\% relative to the NXB rate. The averaged residual signals are consistent with that expected from the cosmic X-ray background. The flux of SN2014J derived from {\it Suzaku} measurements taken in one snapshot at $t=77\pm2$ days after the explosion is consistent with the {\it INTEGRAL} values averaged over the period between $t=$50 and 100 days and also with explosion models of single or double degenerate scenarios. Being sensitive to the total ejecta mass surrounding the radioactive material, the ratio between continuum and line flux in the soft gamma-ray regime might distinguish different progenitor models. The {\it Suzaku} data have been examined with this relation at $t=77\pm2$ days, but could not distinguish models between single and double degenerate-progenitors. We disfavor explosion models with larger $^{56}$Ni masses than 1 $M_\odot$, from our $1 σ$ error on the 170-250 keV X-ray flux of $(1.2\pm0.7) \times10^{-4}$ ph s$^{-1}$ cm$^{-2}$.

preprint2016arXiv

Predicting polarization signatures for double-detonation and delayed-detonation models of Type Ia supernovae

Calculations of synthetic spectropolarimetry are one means to test multi-dimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M$_{\odot}$ white dwarf binary system is too asymmetric to explain the low polarization levels commonly observed in normal Type Ia supernovae. Here, we present polarization simulations for two alternative scenarios: the sub-Chandrasekhar mass double-detonation and the Chandrasekhar mass delayed-detonation model. Specifically, we study a two-dimensional double-detonation model and a three-dimensional delayed-detonation model, and calculate polarization spectra for multiple observer orientations in both cases. We find modest polarization levels ($<$ 1 per cent) for both explosion models. Polarization in the continuum peaks at $\sim$ 0.1$-$0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si ii λ6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O i λ7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt.

preprint2016arXiv

The peculiar Type Ia supernova iPTF14atg: Chandrasekhar-mass explosion or violent merger?

iPTF14atg, a subluminous peculiar Type Ia supernova (SN Ia) similar to SN 2002es, is the first SN Ia for which a strong UV flash was observed in the early-time light curves. This has been interpreted as evidence for a single-degenerate (SD) progenitor system where such a signal is expected from interactions between the SN ejecta and the non-degenerate companion star. Here, we compare synthetic observables of multi-dimensional state-of-the-art explosion models for different progenitor scenarios to the light curves and spectra of iPTF14atg. From our models, we have difficulties explaining the spectral evolution of iPTF14atg within the SD progenitor channel. In contrast, we find that a violent merger of two carbon-oxygen white dwarfs with 0.9 and 0.76 solar masses, respectively, provides an excellent match to the spectral evolution of iPTF14atg from 10d before to several weeks after maximum light. Our merger model does not naturally explain the initial UV flash of iPTF14atg. We discuss several possibilities like interactions of the SN ejecta with the circum-stellar medium and surface radioactivity from a He ignited merger that may be able to account for the early UV emission in violent merger models.

preprint2015arXiv

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 ~1.1 Msun, hybrid WDs in a binary system may easily approach the Chandrasekhar mass (MCh) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-MCh 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 Msun 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-MCh WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-MCh 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 MCh to be on the order of 1 percent of the Galactic SN Ia rate.

preprint2015arXiv

New numerical solver for flows at various Mach numbers

Many problems in stellar astrophysics feature flows at low Mach numbers. Conventional compressible hydrodynamics schemes frequently used in the field have been developed for the transonic regime and exhibit excessive numerical dissipation for these flows. While schemes were proposed that solve hydrodynamics strictly in the low Mach regime and thus restrict their applicability, we aim at developing a scheme that correctly operates in a wide range of Mach numbers. Based on an analysis of the asymptotic behavior of the Euler equations in the low Mach limit we propose a novel scheme that is able to maintain a low Mach number flow setup while retaining all effects of compressibility. This is achieved by a suitable modification of the well-known Roe solver. Numerical tests demonstrate the capability of this new scheme to reproduce slow flow structures even in moderate numerical resolution. Our scheme provides a promising approach to a consistent multidimensional hydrodynamical treatment of astrophysical low Mach number problems such as convection, instabilities, and mixing in stellar evolution.

preprint2015arXiv

Quantitative spectral analysis of the sdB star HD 188112: a helium-core white dwarf progenitor

HD 188112 is a bright (V = 10.2 mag) hot subdwarf B (sdB) star with a mass too low to ignite core helium burning and is therefore considered as a pre-extremely low mass (ELM) white dwarf (WD). ELM WDs (M $\le$ 0.3 Msun) are He-core objects produced by the evolution of compact binary systems. We present in this paper a detailed abundance analysis of HD 188112 based on high-resolution Hubble Space Telescope (HST) near and far-ultraviolet spectroscopy. We also constrain the mass of the star's companion. We use hybrid non-LTE model atmospheres to fit the observed spectral lines and derive the abundances of more than a dozen elements as well as the rotational broadening of metallic lines. We confirm the previous binary system parameters by combining radial velocities measured in our UV spectra with the already published ones. The system has a period of 0.60658584 days and a WD companion with M $\geq$ 0.70 Msun. By assuming a tidally locked rotation, combined with the projected rotational velocity (v sin i = 7.9 $\pm$ 0.3 km s$^{-1}$) we constrain the companion mass to be between 0.9 and 1.3 Msun. We further discuss the future evolution of the system as a potential progenitor of a (underluminous) type Ia supernova. We measure abundances for Mg, Al, Si, P, S, Ca, Ti, Cr, Mn, Fe, Ni, and Zn, as well as for the trans-iron elements Ga, Sn, and Pb. In addition, we derive upper limits for the C, N, O elements and find HD 188112 to be strongly depleted in carbon. We find evidence of non-LTE effects on the line strength of some ionic species such as Si II and Ni II. The metallic abundances indicate that the star is metal-poor, with an abundance pattern most likely produced by diffusion effects.

preprint2015arXiv

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarisation signatures

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarisation typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. (2012) for an explosion triggered during the merger of a 1.1 M$_{\odot}$ and 0.9 M$_{\odot}$ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarisation signal varies significantly with viewing angle. We find that polarisation levels for observers in the equatorial plane are modest ($\lesssim$ 1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarisation and departures from a dominant axis. While the particular model studied here gives a good match to highly-polarised events such as SN 2004dt, it has difficulties in reproducing the low polarisation levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarisation features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

preprint2014arXiv

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).

preprint2014arXiv

Extensive HST Ultraviolet Spectra and Multi-wavelength 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 A_V = 2.0 +/- 0.1 mag. But this picture, with R_V = 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 R_V = 2.6) and roughly half by scattering off LMC-like dust in the circumstellar environment of SN 2014J.

preprint2014arXiv

Testing the role of SNe Ia for Galactic chemical evolution of p-nuclei with 2D models and with s-process seeds at different metallicities

The bulk of p isotopes is created in the 'gamma processes' mainly by sequences of photodisintegrations and beta decays in explosive conditions in Type Ia supernovae (SNIa) or in core collapse supernovae (ccSN). The contribution of different stellar sources to the observed distribution of p-nuclei in the Solar System is still under debate. We explore single degenerate Type Ia supernovae in the framework of two-dimensional SNIa delayed-detonation explosion models. Travaglio et al. (2011, TRV11) discussed the sensitivity of p-nuclei production to different SNIa models, i.e. delayed detonations of different strength, deflagrations, and the dependence on selected s-process seed distributions. Here we present a detailed study of p-process nucleosynthesis occuring in SNIa with s-process seeds at different metallicities. Based on the delayed-detonation model DDT-a of TRV11, we analyze the dependence of p-nucleosynthesis on the s-seed distribution obtained from different strengths of the 13C-pocket. Using a Galactic chemical evolution code (see Travaglio et al. 2004) we study the contribution of SNIa to the solar stable p-nuclei. We find that explosions of Chandrasekhar-mass single degenerate systems produce a large amount of p-nuclei in our Galaxy, both in the range of light (A < 120) and heavy p-nuclei, at almost flat average production factors (within a factor of about 3). We discussed in details p-isotopes such as 94Mo with a behavior diverging from the average, which we attribute to uncertainties in the nuclear data or in SNIa modelling.

preprint2014arXiv

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 <~ 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 ~ 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.

preprint2013arXiv

A Subgrid-scale Model for Deflagration-to-Detonation Transitions in Type Ia Supernova Explosion Simulations - Numerical implementation

A promising model for normal Type Ia supernova (SN Ia) explosions are delayed detonations of Chandrasekhar-mass white dwarfs, in which the burning starts out as a subsonic deflagration and turns at a later phase of the explosion into a supersonic detonation. The mechanism of the underlying deflagration-to-detonation transition (DDT) is unknown in detail, but necessary conditions have been determined recently. The region of detonation initiation cannot be spatially resolved in multi-dimensional full-star simulations of the explosion. We develop a subgrid-scale (SGS) model for DDTs in thermonuclear supernova simulations that is consistent with the currently known constraints. The probability for a DDT to occur is calculated from the distribution of turbulent velocities measured on the grid scale in the vicinity of the flame and the fractal flame surface area that satisfies further physical constraints, such as fuel fraction and fuel density. The implementation of our DDT criterion provides a solid basis for simulations of thermonuclear supernova explosions in the delayed detonation scenario. It accounts for the currently known necessary conditions for the transition and avoids the inclusion of resolution-dependent quantities in the model. The functionality of our DDT criterion is demonstrated on the example of one three-dimensional thermonuclear supernova explosion simulation.

preprint2013arXiv

Combustion of a neutron star into a strange quark star: The neutrino signal

There are strong indications that the process of conversion of a neutron star into a strange quark star proceeds as a strong deflagration implying that in a few milliseconds almost the whole star is converted. Starting from the three-dimensional hydrodynamic simulations of the combustion process which provide the temperature profiles inside the newly born strange star, we calculate for the first time the neutrino signal that is to be expected if such a conversion process takes place. The neutrino emission is characterized by a luminosity and a duration that is typical for the signal expected from protoneutron stars and represents therefore a powerful source of neutrinos which could be possibly directly detected in case of events occurring close to our Galaxy. We discuss moreover possible connections between the birth of strange stars and explosive phenomena such as supernovae and gamma-ray-bursts.

preprint2013arXiv

Gamma-ray diagnostics of Type Ia supernovae: Predictions of observables from three-dimensional modeling

Besides the fact that the gamma-ray emission due to radioactive decays is responsible for powering the light curves of Type Ia supernovae (SNe Ia), gamma rays themselves are of particular interest as a diagnostic tool because they provide a direct way to obtain deeper insights into the nucleosynthesis and the kinematics of these explosion events. Focusing on two of the most broadly discussed SN Ia progenitor scenarios - a delayed detonation in a Chandrasekhar-mass white dwarf (WD) and a violent merger of two WDs - we use three-dimensional explosion models and perform radiative transfer simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of 56Ni and have similar optical properties that are in reasonable agreement with the recently observed supernova SN 2011fe. In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be rather different. The almost direct connection of the emission of gamma rays to fundamental physical processes occuring in SNe Ia permits additional constraints concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV missions such as GRIPS will resolve all spectral details only for nearby SNe Ia, but hardness ratio and light curve measurements still allow for a distinction of the two different models at 10 and 16 Mpc for an exposure time of 10^6 s, respectively. The possibility to detect the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN Ia detections in the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays.

preprint2013arXiv

Predicting the amount of hydrogen stripped by the supernova explosion for SN 2002cx-like SNe Ia

The most favored progenitor scenarios for Type Ia supernovae (SNe Ia) involve the single-degenerate (SD) scenario and the double-degenerate scenario. The absence of stripped hydrogen (H) in the nebular spectra of SNe Ia challenges the SD progenitor models. Recently, it was shown that pure deflagration explosion models of Chandrasekhar-mass white dwarfs ignited off-center reproduce the characteristic observational features of 2002cx-like SNe Ia very well. In this work we predict, for the first time, the amount of stripped H for the off-center pure deflagration explosions. We find that their low kinetic energies lead to inefficient H mass stripping (less than 0.01 M_sun), indicating that the stripped H may be hidden in (observed) late-time spectra of SN 2002cx-like SNe Ia.

preprint2013arXiv

SN 2010lp---a Type Ia supernova from a violent merger of two carbon--oxygen White Dwarfs

SN 2010lp is a subluminous Type Ia supernova (SN Ia) with slowly-evolving lightcurves. Moreover, it is the only subluminous SN Ia observed so far that shows narrow emission lines of OI in late-time spectra, indicating unburned oxygen close to the centre of the ejecta. Most explosion models for SNe Ia cannot explain the narrow OI emission. Here, we present hydrodynamic explosion and radiative transfer calculations showing that the violent merger of two carbon-oxygen white dwarfs of 0.9 and 0.76 solar masses, respectively, adequately reproduces the early-time observables of SN 2010lp. Moreover, our model predicts oxygen close to the centre of the explosion ejecta, a pre-requisite for narrow OI emission in nebular spectra as observed in SN 2010lp.

preprint2012arXiv

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. 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 proto-typical 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 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 solar masses -- 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 influence the late-time observables of 2002cx-like SNe.

preprint2012arXiv

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-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.

preprint2012arXiv

Monte Carlo radiation hydrodynamics: methods, tests and application to supernova Type Ia ejecta

In astrophysical systems, radiation-matter interactions are important in transferring energy and momentum between the radiation field and the surrounding material. This coupling often makes it necessary to consider the role of radiation when modelling the dynamics of astrophysical fluids. During the last few years, there have been rapid developments in the use of Monte Carlo methods for numerical radiative transfer simulations. Here, we present an approach to radiation hydrodynamics that is based on coupling Monte Carlo radiative transfer techniques with finite-volume hydrodynamical methods in an operator-split manner. In particular, we adopt an indivisible packet formalism to discretize the radiation field into an ensemble of Monte Carlo packets and employ volume-based estimators to reconstruct the radiation field characteristics. In this paper the numerical tools of this method are presented and their accuracy is verified in a series of test calculations. Finally, as a practical example, we use our approach to study the influence of the radiation-matter coupling on the homologous expansion phase and the bolometric light curve of Type Ia supernova explosions.

preprint2012arXiv

Normal Type Ia supernovae from violent mergers of white dwarf binaries

One of the most important questions regarding the progenitor systems of Type Ia supernovae (SNe Ia) is whether mergers of two white dwarfs can lead to explosions that reproduce observations of normal events. Here we present a fully three-dimensional simulation of a violent merger of two carbon-oxygen white dwarfs with masses of $0.9 \mathrm{M_\odot}$ and $1.1 \mathrm{M_\odot}$ combining very high resolution and exact initial conditions. A well-tested combination of codes is used to study the system. We start with the dynamical inspiral phase and follow the subsequent thermonuclear explosion under the plausible assumption that a detonation forms in the process of merging. We then perform detailed nucleosynthesis calculations and radiative transfer simulations to predict synthetic observables from the homologously expanding supernova ejecta. We find that synthetic color lightcurves of our merger, which produces about $0.62 \mathrm{M_\odot}$ of $^{56}\mathrm{Ni}$, show good agreement with those observed for normal SNe Ia in all wave bands from U to K. Line velocities in synthetic spectra around maximum light also agree well with observations. We conclude, that violent mergers of massive white dwarfs can closely resemble normal SNe Ia. Therefore, depending on the number of such massive systems available these mergers may contribute at least a small fraction to the observed population of normal SNe Ia.

preprint2012arXiv

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

We investigate the brightness distribution expected for thermonuclear explosions that might result from the ignition of a detonation during the violent merger of white dwarf (WD) binaries. Determining their brightness distribution is critical for evaluating whether such an explosion model could be responsible for a significant fraction of the observed population of SNe Ia. We argue that the brightness of an explosion realized via the violent merger model is determined by the mass of nickel 56 produced in the detonation of the primary WD. We use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass and expected peak bolometric brightness. We use this relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a brightness 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 Msun from a slightly-evolved helium star companion. We also find that with high probability, violent WD mergers involving the most massive primaries (which should produce bright SNe) have delay times <500 Myr (Abridged).

preprint2012arXiv

Prospect of Studying Hard X- and Gamma-Rays from Type Ia Supernovae

We perform multi-dimensional, time-dependent radiation transfer simulations for hard X-ray and gamma-ray emissions, following radioactive decays of 56Ni and 56Co, for two-dimensional delayed detonation models of Type Ia supernovae (SNe Ia). The synthetic spectra and light curves are compared with the sensitivities of current and future observatories for an exposure time of 10^6 seconds. The non-detection of the gamma-ray signal from SN 2011fe at 6.4 Mpc by SPI on board INTEGRAL places an upper limit for the mass of 56Ni of \lesssim 1.0 Msun, independently from observations in any other wavelengths. Signals from the newly formed radioactive species have not been convincingly measured yet from any SN Ia, but the future X-ray and gamma-ray missions are expected to deepen the observable horizon to provide the high energy emission data for a significant SN Ia sample. We predict that the hard X-ray detectors on board NuStar (launched in 2012) or ASTRO-H (scheduled for launch in 2014) will reach to SNe Ia at \sim15 Mpc, i.e., one SN every few years. Furthermore, according to the present results, the soft gamma-ray detector on board ASTRO-H will be able to detect the 158 keV line emission up to \sim25 Mpc, i.e., a few SNe Ia per year. Proposed next generation gamma-ray missions, e.g., GRIPS, could reach to SNe Ia at \sim20 - 35 Mpc by MeV observations. Those would provide new diagnostics and strong constraints on explosion models, detecting rather directly the main energy source of supernova light.

preprint2012arXiv

Stellar GADGET: A smooth particle hydrodynamics code for stellar astrophysics and its application to Type Ia supernovae from white dwarf mergers

Mergers of two carbon-oxygen white dwarfs have long been suspected to be progenitors of Type Ia Supernovae. Here we present our modifications to the cosmological smoothed particle hydrodynamics code Gadget to apply it to stellar physics including but not limited to mergers of white dwarfs. We demonstrate a new method to map a one-dimensional profile of an object in hydrostatic equilibrium to a stable particle distribution. We use the code to study the effect of initial conditions and resolution on the properties of the merger of two white dwarfs. We compare mergers with approximate and exact binary initial conditions and find that exact binary initial conditions lead to a much more stable binary system but there is no difference in the properties of the actual merger. In contrast, we find that resolution is a critical issue for simulations of white dwarf mergers. Carbon burning hotspots which may lead to a detonation in the so-called violent merger scenario emerge only in simulations with sufficient resolution but independent of the type of binary initial conditions. We conclude that simulations of white dwarf mergers which attempt to investigate their potential for Type Ia supernovae should be carried out with at least 10^6 particles.

preprint2012arXiv

Three-dimensional simulations of the interaction between Type Ia supernova ejecta and their main sequence companions

The identity of the progenitor systems of SNe Ia is still uncertain. In the single-degenerate (SD) scenario, the interaction between the SN blast wave and the outer layers of a main sequence (MS) companion star strips off H-rich material which is then mixed into the ejecta. Strong contamination of the SN ejecta with stripped material could lead to a conflict with observations of SNe Ia. This constrains the SD progenitor model. In this work, our previous simulations based on simplified progenitor donor stars have been updated by adopting more realistic progenitor-system models that result from fully detailed, state-of-the-art binary evolution calculations. We use Eggleton's stellar evolution code including the optically thick accretion wind model and the possibility of the effects of accretion disk instabilities to obtain realistic models of companions for different progenitor systems. The impact of the SN blast wave on these companion stars is followed in three-dimensional hydrodynamic simulations employing the SPH code GADGET3. We find that the stripped masses range from 0.11 to 0.18 M_sun. The kick velocity is between 51 and 105 km/s. We find that the stripped mass and kick velocity depend on the ratio of the orbital separation to the radius of a companion. They can be fitted by a power law for a given companion model. However, the structure of the companion star is also important for the amount of stripped material. With more realistic companion star models than in previous studies, our simulations show that the H masses stripped from companions are inconsistent with the best observational limits (< 0.01 M_sun) derived from nebular spectra. However, a rigorous forward modeling based on impact simulations with radiation transfer is required to reliably predict observable signatures of the stripped H and to conclusively assess the viability of the considered SN Ia progenitor scenario.

preprint2011arXiv

2D simulations of the double-detonation model for thermonuclear transients from low-mass carbon-oxygen white dwarfs

Thermonuclear explosions may arise in binaries in which a CO white dwarf (WD) accretes He from a companion. If the accretion rate allows a sufficiently large mass of He to accumulate prior to ignition of nuclear burning, the He surface layer may detonate, giving rise to an astrophysical transient. Detonation of the He layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation at the edge of the CO WD or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either ignition mechanism works, the two detonations can release sufficient energy to completely unbind the WD. Here we extend our 2D studies of this double-detonation model to low-mass CO WDs. We investigate the feasibility of triggering a secondary core detonation by shock convergence in low-mass CO WDs and the observable consequences of such a detonation. Our results suggest that core detonation is probable, even for the lowest CO core masses realized in nature. We compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur and compare these to models in which no CO detonation was allowed to occur. If significant shock compression of the CO WD occurs prior to detonation, explosion of the CO WD can produce a sufficiently large mass of radioactive iron-group nuclei to affect the light curves. In particular, this can lead to relatively slow post-maximum decline. If the secondary detonation is edge-lit, however, the CO WD explosion primarily yields intermediate-mass elements that affect the observables more subtly. In this case, NIR observations and detailed spectroscopic analysis would be needed to determine whether core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E.

preprint2011arXiv

Confronting 2D delayed-detonation models with light curves and spectra of Type Ia supernovae

We compare models for Type Ia supernova (SN Ia) light curves and spectra with an extensive set of observations. The models come from a recent survey of 44 two-dimensional delayed-detonation models computed by Kasen, Roepke & Woosley (2009), each viewed from multiple directions. The data include optical light curves of 251 SNe Ia and 2231 low-dispersion spectra from the Center for Astrophysics, plus data from the literature. The analysis uses standard techniques employed by observers, including MLCS2k2, SALT2, and SNooPy for light-curve analysis, and the Supernova Identification (SNID) code of Blondin & Tonry for spectroscopic comparisons to assess how well the models match the data. We show that the models that match observed spectra best lie systematically on the observed width-luminosity relation. Conversely, we reject six models with highly asymmetric ignition conditions and a large amount (>1 M_sun) of synthesized 56Ni that yield poor matches to observed SN Ia spectra. More subtle features of the comparison include the general difficulty of the models to match the U-band flux at early times, caused by a hot ionized ejecta that affect the subsequent redistribution of flux at longer wavelengths. We examine ways in which the asymptotic kinetic energy of the explosion affects both the predicted velocity and velocity gradient in the Si II and Ca II lines. Models with an asymmetric distribution of 56Ni are found to result in a larger variation of photometric and spectroscopic properties with viewing angle, regardless of the initial ignition setup. We discuss more generally whether highly anisotropic ignition conditions are ruled out by observations, and how detailed comparisons between models and observations involving both light curves and spectra can lead to a better understanding of SN Ia explosion mechanisms.

preprint2011arXiv

Type Ia supernova diversity: white dwarf central density as a secondary parameter in three-dimensional delayed detonation models

Delayed detonations of Chandrasekhar-mass white dwarfs (WDs) have been very successful in explaining the spectra, light curves, and the width-luminosity relation of spectroscopically normal Type Ia supernovae (SNe Ia). The ignition of the thermonuclear deflagration flame at the end of the convective carbon "simmering" phase in the core of the WD is still not well understood and much about the ignition kernel distribution remains unknown. Furthermore, the central density at the time of ignition depends on the still uncertain screened carbon fusion reaction rates, the accretion history and cooling time of the progenitor, and the composition. We present the results of twelve high-resolution three-dimensional delayed detonation SN Ia explosion simulations that employ a new criterion to trigger the deflagration to detonation transition (DDT). All simulations trigger our DDT criterion and the resulting delayed detonations unbind the star. We find a trend of increasing iron group element (IGE) production with increasing central density for bright, faint, and intermediate SNe. The total 56Ni yield, however, remains more or less constant, even though increased electron captures at high density result in a decreasing 56Ni mass fraction of the IGE material. We attribute this to an approximate balance of 56Ni producing and destroying effects. The deflagrations that were ignited at higher density initially have a faster growth rate of subgrid-scale turbulence. Hence, the effective flame speed increases faster, which triggers the DDT criterion earlier, at a time when the central density of the expanded star is higher. This leads to an overall increase of IGE production, which off-sets the percental reduction of 56Ni due to neutronization.

preprint2011arXiv

Violent mergers of nearly equal-mass white dwarf as progenitors of subluminous Type Ia supernovae

The origin of subluminous Type Ia supernovae (SNe Ia) has long eluded any explanation, as all Chandrasekhar-mass models have severe problems reproducing them. Recently, it has been proposed that violent mergers of two white dwarfs of 0.9 M_sun could lead to subluminous SNe Ia events that resemble 1991bg-like SNe~Ia. Here we investigate whether this scenario still works for mergers of two white dwarfs with a mass ratio smaller than one. We aim to determine the range of mass ratios for which a detonation still forms during the merger, as only those events will lead to a SN Ia. This range is an important ingredient for population synthesis and one decisive point to judge the viability of the scenario. In addition, we perform a resolution study of one of the models. Finally we discuss the connection between violent white dwarf mergers with a primary mass of 0.9 M_sun and 1991bg-like SNe Ia. The latest version of the smoothed particle hydrodynamics code Gadget3 is used to evolve binary systems with different mass ratios until they merge. We analyze the result and look for hot spots in which detonations can form. We show that mergers of two white dwarfs with a primary white dwarf mass of ~0.9 M_sun and a mass ratio larger than about $0.8$ robustly reach the conditions we require to ignite a detonation and thus produce thermonuclear explosions during the merger itself. We also find that while our simulations do not yet completely resolve the hot spots, increasing the resolution leads to conditions that are even more likely to ignite detonations. (abridged)

preprint2010arXiv

An asymmetric explosion as the origin of spectral evolution diversity in type Ia supernovae

Type Ia Supernovae (SNe Ia) form an observationally uniform class of stellar explosions, in that more luminous objects have smaller decline-rates. This one-parameter behavior allows SNe Ia to be calibrated as cosmological `standard candles', and led to the discovery of an accelerating Universe. Recent investigations, however, have revealed that the true nature of SNe Ia is more complicated. Theoretically, it has been suggested that the initial thermonuclear sparks are ignited at an offset from the centre of the white-dwarf (WD) progenitor, possibly as a result of convection before the explosion. Observationally, the diversity seen in the spectral evolution of SNe Ia beyond the luminosity decline-rate relation is an unresolved issue. Here we report that the spectral diversity is a consequence of random directions from which an asymmetric explosion is viewed. Our findings suggest that the spectral evolution diversity is no longer a concern in using SNe Ia as cosmological standard candles. Furthermore, this indicates that ignition at an offset from the centre of is a generic feature of SNe Ia.

preprint2010arXiv

Detonations in Sub-Chandrasekhar Mass C+O White Dwarfs

Explosions of sub-Chandrasekhar-mass white dwarfs are one alternative to the standard Chandrasekhar-mass model of Type Ia supernovae. They are interesting since binary systems with sub-Chandrasekhar-mass primary white dwarfs should be common and this scenario would suggest a simple physical parameter which determines the explosion brightness, namely the mass of the exploding white dwarf. Here we perform one-dimensional hydrodynamical simulations, associated post-processing nucleosynthesis and multi-wavelength radiation transport calculations for pure detonations of carbon-oxygen white dwarfs. The light curves and spectra we obtain from these simulations are in good agreement with observed properties of Type Ia supernovae. In particular, for white dwarf masses from 0.97 - 1.15 Msun we obtain 56Ni masses between 0.3 and 0.8 Msun, sufficient to capture almost the complete range of Type Ia supernova brightnesses. Our optical light curve rise times, peak colours 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 Angstroms. We therefore suggest that sub-Chandrasekhar mass explosions are a viable model for Type Ia supernovae 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 white dwarf.

preprint2010arXiv

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 helium 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 prerequisite 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 nuclear 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 that secondary core detonations are triggered for all of the simulated models, ranging in core mass from 0.810 up to 1.385 M_solar with corresponding shell masses from 0.126 down to 0.0035 M_solar. This implies that, as soon as a detonation triggers in a helium shell covering a carbon/oxygen white dwarf, a subsequent core detonation is virtually inevitable.

preprint2010arXiv

Double-detonation sub-Chandrasekhar supernovae: synthetic observables for minimum helium shell mass models

Abridged. In the double detonation scenario for Type Ia supernovae (SNe Ia) a detonation initiates in a shell of He-rich material accreted from a companion star by a sub-Chandrasekhar-mass White Dwarf (WD). This shell detonation drives a shock front into the carbon-oxygen (C/O) WD that triggers a secondary detonation in the core. The core detonation results in a complete disruption of the WD. Earlier studies concluded that this scenario has difficulties in accounting for the observed properties of SNe Ia since the explosion ejecta are surrounded by the products of explosive He burning in the shell. Recently, it was proposed that detonations might be possible for much less massive He shells than previously assumed. Moreover, it was shown that even detonations of these minimum He shell masses robustly trigger detonations of the C/O core. Here we present time-dependent multi-wavelength radiative transfer calculations for models with minimum He shell mass and derive synthetic observables for both the optical and γ-ray spectral regions. These differ strongly from those found in earlier simulations of sub-Chandrasekhar-mass explosions in which more massive He shells were considered. Our models predict light curves which cover both the range of brightnesses and the rise and decline times of observed SNe Ia. However, their colours and spectra do not match the observations. In particular, their B-V colours are generally too red. We show that this discrepancy is mainly due to the composition of the burning products of the He shell of our models which contain significant amounts of Ti and Cr. Using a toy model, we also show that the burning products of the He shell depend crucially on its initial composition. This leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed SNe Ia may still be feasible but further study of the shell properties is required.

preprint2010arXiv

Modeling the Diversity of Type Ia Supernova Explosions

Type Ia supernovae (SNe Ia) are a prime tool in observational cosmology. A relation between their peak luminosities and the shapes of their light curves allows to infer their intrinsic luminosities and to use them as distance indicators. This relation has been established empirically. However, a theoretical understanding is necessary in order to get a handle on the systematics in SN Ia cosmology. Here, a model reproducing the observed diversity of normal SNe Ia is presented. The challenge in the numerical implementation arises from the vast range of scales involved in the physical mechanism. Simulating the supernova on scales of the exploding white dwarf requires specific models of the microphysics involved in the thermonuclear combustion process. Such techniques are discussed and results of simulations are presented.

preprint2010arXiv

Nucleosynthesis in Two-Dimensional Delayed Detonation Models of Type Ia Supernova Explosions

The nucleosynthetic characteristics of various explosion mechanisms of Type Ia supernovae (SNe Ia) is explored based on three two-dimensional explosion simulations representing extreme cases: a pure turbulent deflagration, a delayed detonation following an approximately spherical ignition of the initial deflagration, and a delayed detonation arising from a highly asymmetric deflagration ignition. Apart from this initial condition, the deflagration stage is treated in a parameter-free approach. The detonation is initiated when the turbulent burning enters the distributed burning regime. This occurs at densities around $10^{7}$ g cm$^{-3}$ -- relatively low as compared to existing nucleosynthesis studies for one-dimensional spherically symmetric models. The burning in these multidimensional models is different from that in one-dimensional simulations as the detonation wave propagates both into unburned material in the high density region near the center of a white dwarf and into the low density region near the surface. Thus, the resulting yield is a mixture of different explosive burning products, from carbon-burning products at low densities to complete silicon-burning products at the highest densities, as well as electron-capture products synthesized at the deflagration stage. In contrast to the deflagration model, the delayed detonations produce a characteristic layered structure and the yields largely satisfy constraints from Galactic chemical evolution. In the asymmetric delayed detonation model, the region filled with electron capture species (e.g., $^{58}$Ni, $^{54}$Fe) is within a shell, showing a large off-set, above the bulk of $^{56}$Ni distribution, while species produced by the detonation are distributed more spherically (abridged).

preprint2010arXiv

Nucleosynthetic post-processing of Type Ia supernovae with variable tracer masses

The post-processing of passively advected Lagrangian tracer particles is still the most common way for obtaining detailed nucleosynthetic yield predictions of Type Ia supernova (SN Ia) hydrodynamical simulations. Historically, tracer particles of constant mass are employed. However, intermediate mass elements, such as e.g. Ne, Mg, Al, or Si, are typically synthesized in the outer layers of SNe Ia, where due to the lower initial density a constant mass tracer distribution results in poor resolution of the spatial morphology of the abundance distribution. We show how to alleviate this problem with a suitably chosen distribution of variable tracer particle masses. We also present results of the convergence of integrated nucleosynthetic yields with increasing tracer particle number. We find that the yields of the most abundant species (mass fraction > 10E-5) are reasonably well predicted for a tracer number as small as 32 per axis and direction. Convergence for isotopes produced in regions where a constant tracer mass implementation results in poor spatial resolution can be improved by suitably choosing tracers of variable mass.

preprint2009arXiv

Turbulence in a 3D deflagration model for type Ia SNe: II. Intermittency and the deflagration-to-detonation transition probability

The delayed detonation model describes the observational properties of the majority of type Ia supernovae very well. Using numerical data from a three-dimensional deflagration model for type Ia supernovae, the intermittency of the turbulent velocity field and its implications on the probability of a deflagration-to-detonation (DDT) transition are investigated. From structure functions of the turbulent velocity fluctuations, we determine intermittency parameters based on the log-normal and the log-Poisson models. On the other hand, the analysis of the turbulent velocity fluctuations in the vicinity of the flame front by Roepke suggests a much higher probability of large velocity fluctuations on the grid scale in comparison to the log-normal intermittency model. Following Pan et al., we computed probability density functions for a DDT for the different distributions. Assuming that a DDT can occur in the stirred flame regime, as proposed by Woosley et al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 seconds after the beginning of the deflagration phase for the multi-spot ignition scenario used in the simulation. However, the probability drops to virtually zero if a DDT is further constrained by the requirement that the turbulent velocity fluctuations reach about 500 km/s. Under this condition, delayed detonations are only possible if the distribution of the velocity fluctuations is not log-normal. From our calculations follows that the distribution obtained by Roepke allow for multiple DDTs around 0.8 seconds after ignition at a transition density close to 1x10^7 g/cm^3.

preprint2009arXiv

Turbulence in a three-dimensional deflagration model for Type Ia supernovae: I. Scaling properties

We analyze the statistical properties of the turbulent velocity field in the deflagration model for Type Ia supernovae. In particular, we consider the question of whether turbulence is isotropic and consistent with the Kolmogorov theory at small length scales. Using numerical data from a high-resolution simulation of a thermonuclear supernova explosion, spectra of the turbulence energy and velocity structure functions are computed. We show that the turbulent velocity field is isotropic at small length scales and follows a scaling law that is consistent with the Kolmogorov theory until most of the nuclear fuel is burned. At length scales greater than a certain characteristic scale, turbulence becomes anisotropic. Here, the radial velocity fluctuations follow the scaling law of the Rayleigh-Taylor instability, whereas the angular component still obeys Kolmogorov scaling. In the late phase of the explosion, this characteristic scale drops below the numerical resolution of the simulation. The analysis confirms that a subgrid-scale model for the unresolved turbulence energy is required for the consistent calculation of the flame speed in deflagration models of Type Ia supernovae, and that the assumption of isotropy on these scales is appropriate.

preprint2008arXiv

Multi-dimensional numerical simulations of type Ia supernova explosions

The major role type Ia supernovae play in many fields of astrophysics and in particular in cosmological distance determinations calls for self-consistent models of these events. Since their mechanism is believed to crucially depend on phenomena that are inherently three-dimensional, self-consistent numerical models of type Ia supernovae must be multi-dimensional. This field has recently seen a rapid development, which is reviewed in this article. The different modeling approaches are discussed and as an illustration a particular explosion model -- the deflagration model -- in a specific numerical implementation is presented in greater detail. On this exemplary case, the procedure of validating the model on the basis of comparison with observations is discussed as well as its application to study questions arising from type Ia supernova cosmology.

preprint2006arXiv

Off-center ignition in type Ia supernova: I. Initial evolution and implications for delayed detonation

The explosion of a carbon-oxygen white dwarf as a Type Ia supernova is known to be sensitive to the manner in which the burning is ignited. Studies of the pre-supernova evolution suggest asymmetric, off-center ignition, and here we explore its consequences in two- and three-dimensional simulations. Compared with centrally ignited models, one-sided ignitions initially burn less and release less energy. For the distributions of ignition points studied, ignition within two hemispheres typically leads to the unbinding of the white dwarf, while ignition within a small fraction of one hemisphere does not. We also examine the spreading of the blast over the surface of the white dwarf that occurs as the first plumes of burning erupt from the star. In particular, our studies test whether the collision of strong compressional waves can trigger a detonation on the far side of the star as has been suggested by Plewa et al. (2004). The maximum temperature reached in these collisions is sensitive to how much burning and expansion has already gone on, and to the dimensionality of the calculation. Though detonations are sometimes observed in 2D models, none ever happens in the corresponding 3D calculations. Collisions between the expansion fronts of multiple bubbles also seem, in the usual case, unable to ignite a detonation. "Gravitationally confined detonation" is therefore not a robust mechanism for the explosion. Detonation may still be possible in these models however, either following a pulsation or by spontaneous detonation if the turbulent energy is high enough.

preprint2006arXiv

Surface detonation in type Ia supernova explosions?

We explore the evolution of thermonuclear supernova explosions when the progenitor white dwarf star ignites asymmetrically off-center. Several numerical simulations are carried out in two and three dimensions to test the consequences of different initial flame configurations such as spherical bubbles displaced from the center, more complex deformed configurations, and teardrop-shaped ignitions. The burning bubbles float towards the surface while releasing energy due to the nuclear reactions. If the energy release is too small to gravitationally unbind the star, the ash sweeps around it, once the burning bubble approaches the surface. Collisions in the fuel on the opposite side increase its temperature and density and may -- in some cases -- initiate a detonation wave which will then propagate inward burning the core of the star and leading to a strong explosion. However, for initial setups in two dimensions that seem realistic from pre-ignition evolution, as well as for all three-dimensional simulations the collimation of the surface material is found to be too weak to trigger a detonation.

preprint2004arXiv

On the Stability of Thermonuclear Burning Fronts in Type Ia Supernovae

The propagation of cellularly stabilized thermonuclear flames is investigated by means of numerical simulations. In Type Ia supernova explosions the corresponding burning regime establishes at scales below the Gibson length. The cellular flame stabilization - which is a result of an interplay between the Landau-Darrieus instability and a nonlinear stabilization mechanism - is studied for the case of propagation into quiescent fuel as well as interaction with vortical fuel flows. Our simulations indicate that in thermonuclear supernova explosions stable cellular flames develop around the Gibson scale and that deflagration-to-detonation transition is unlikely to be triggered from flame evolution effects here.

F. K. Roepke

What is connected

Connect this record

See the researcher in context

Building this map preview

55 published item(s)

NLTE spectral modelling for a carbon-oxygen and helium white-dwarf merger as a Ca-rich transient candidate

Bipolar planetary nebulae from common envelope evolution of binary stars

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Models of pulsationally assisted gravitationally confined detonations with different ignition conditions

Type Iax supernovae from deflagrations in Chandrasekhar mass white dwarfs

Common-envelope evolution with an asymptotic giant branch star

Fully compressible simulations of waves and core convection in main-sequence stars

Long-term evolution of a magnetic massive merger product

Nucleosynthesis imprints from different Type Ia Supernova explosion scenarios and implications for galactic chemical evolution

Photometric detection of internal gravity waves in upper main-sequence stars. II. Combined TESS photometry and high-resolution spectroscopy

White dwarf deflagrations for Type Iax supernovae: Polarisation signatures from the explosion and companion interaction

Discovery of an Exceptionally Strong $β$-Decay Transition of $^{20}$F and Implications for the Fate of Intermediate-Mass Stars

Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae

The interaction of core-collapse supernova ejecta with a stellar companion

Measurements of the Soft Gamma-ray Emission from SN2014J with Suzaku

Predicting polarization signatures for double-detonation and delayed-detonation models of Type Ia supernovae

The peculiar Type Ia supernova iPTF14atg: Chandrasekhar-mass explosion or violent merger?

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

New numerical solver for flows at various Mach numbers

Quantitative spectral analysis of the sdB star HD 188112: a helium-core white dwarf progenitor

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarisation signatures

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

Extensive HST Ultraviolet Spectra and Multi-wavelength Observations of SN 2014J in M82 Indicate Reddening and Circumstellar Scattering by Typical Dust

Testing the role of SNe Ia for Galactic chemical evolution of p-nuclei with 2D models and with s-process seeds at different metallicities

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

A Subgrid-scale Model for Deflagration-to-Detonation Transitions in Type Ia Supernova Explosion Simulations - Numerical implementation

Combustion of a neutron star into a strange quark star: The neutrino signal

Gamma-ray diagnostics of Type Ia supernovae: Predictions of observables from three-dimensional modeling

Predicting the amount of hydrogen stripped by the supernova explosion for SN 2002cx-like SNe Ia

SN 2010lp---a Type Ia supernova from a violent merger of two carbon--oxygen White Dwarfs

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

Constraining Type Ia supernova models: SN 2011fe as a test case

Monte Carlo radiation hydrodynamics: methods, tests and application to supernova Type Ia ejecta

Normal Type Ia supernovae from violent mergers of white dwarf binaries

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

Prospect of Studying Hard X- and Gamma-Rays from Type Ia Supernovae

Stellar GADGET: A smooth particle hydrodynamics code for stellar astrophysics and its application to Type Ia supernovae from white dwarf mergers

Three-dimensional simulations of the interaction between Type Ia supernova ejecta and their main sequence companions

2D simulations of the double-detonation model for thermonuclear transients from low-mass carbon-oxygen white dwarfs

Confronting 2D delayed-detonation models with light curves and spectra of Type Ia supernovae

Type Ia supernova diversity: white dwarf central density as a secondary parameter in three-dimensional delayed detonation models

Violent mergers of nearly equal-mass white dwarf as progenitors of subluminous Type Ia supernovae

An asymmetric explosion as the origin of spectral evolution diversity in type Ia supernovae

Detonations in Sub-Chandrasekhar Mass C+O White Dwarfs

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

Double-detonation sub-Chandrasekhar supernovae: synthetic observables for minimum helium shell mass models

Modeling the Diversity of Type Ia Supernova Explosions

Nucleosynthesis in Two-Dimensional Delayed Detonation Models of Type Ia Supernova Explosions

Nucleosynthetic post-processing of Type Ia supernovae with variable tracer masses

Turbulence in a 3D deflagration model for type Ia SNe: II. Intermittency and the deflagration-to-detonation transition probability

Turbulence in a three-dimensional deflagration model for Type Ia supernovae: I. Scaling properties

Multi-dimensional numerical simulations of type Ia supernova explosions

Off-center ignition in type Ia supernova: I. Initial evolution and implications for delayed detonation

Surface detonation in type Ia supernova explosions?

On the Stability of Thermonuclear Burning Fronts in Type Ia Supernovae