Common envelope (CE) phases in binary systems where the primary star reaches the tip of the red giant branch are discussed as a formation scenario for hot subluminous B-type (sdB) stars. For some of these objects, observations point to very low-mass companions. In hydrodynamical CE simulations with the moving-mesh code AREPO, we test whether low-mass objects can successfully unbind the envelope. The success of envelope removal in our simulations critically depends on whether or not the ionization energy released by recombination processes in the expanding material is taken into account. If this energy is thermalized locally, envelope ejection eventually leading to the formation of an sdB star is possible with companion masses down to the brown dwarf range. For even lower companion masses approaching the regime of giant planets, however, envelope removal becomes increasingly difficult or impossible to achieve. Our results are consistent with current observational constraints on companion masses of sdB stars. Based on a semianalytic model, we suggest a new criterion for the lowest companion mass that is capable of triggering a dynamical response of the primary star thus potentially facilitating the ejection of a common envelope. This gives an estimate consistent with the findings of our hydrodynamical simulations.

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.

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.

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 $\tau$ 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.

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.

Sub-Chandrasekhar mass white dwarfs accreting a helium shell on a carbon-oxygen core are potential progenitors of normal Type Ia supernovae. This work focuses on the details of the onset of the carbon detonation in the double detonation sub-Chandrasekhar model. In order to simulate the influence of core-shell mixing on the carbon ignition mechanism, the helium shell and its detonation are followed with an increased resolution compared to the rest of the star treating the propagation of the detonation wave more accurately. This significantly improves the predictions of the nucleosynthetic yields from the helium burning. The simulations were carried out with the AREPO code. A carbon-oxygen core with a helium shell was set up in one dimension and mapped to three dimensions. We ensured the stability of the white dwarf with a relaxation step before the hydrodynamic detonation simulation started. Synthetic observables were calculated with the radiative transfer code ARTIS. An ignition mechanism of the carbon detonation was observed, which received little attention before. In this "scissors mechanism", the impact the helium detonation wave has on unburnt material when converging opposite to its ignition spot is strong enough to ignite a carbon detonation. This is possible in a carbon enriched transition region between the core and shell. The detonation mechanism is found to be sensitive to details of the core-shell transition and our models illustrate the need to consider core-shell mixing taking place during the accretion process. Even though the detonation ignition mechanism differs form the converging shock mechanism, the differences in the synthetic observables are not significant. Though they do not fit observations better than previous simulations, they illustrate the need for multi-dimensional simulations.

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.

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.

We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about 15% are partial thermonuclear explosions, the model is able to reproduce the solar abundances of several important and problematic isotopes including $^{48}$Ca, $^{50}$Ti and $^{54}$Cr together with $^{58}$Fe, $^{64}$Ni, $^{82}$Se and $^{86}$Kr and several of the Zn--Zr isotopes. A model in which no cECSNe occur, only tECSNe with low-mass FeCCSNe or rotating massive stars, proves also very successful at reproducing the solar abundances for these isotopes. Despite the small mass range for the progenitors of ECSNe and low-mass FeCCSNe, the large production factors suffice for the solar inventory of the above isotopes. Our model is compelling because it introduces no new tensions with the solar abundance distribution for a Milky Way model -- only tending to improve the model predictions for several isotopes. The proposed astrophysical production model thus provides a natural and elegant way to explain one of the last uncharted territories on the periodic table of astrophysical element production.

Recent progress in the three-dimensional modeling of supernovae (SN) has shown the importance of asymmetries for the explosion. This calls for a reconsideration of the modeling of the subsequent phase, the supernova remnant (SNR), which has commonly relied on simplified ejecta models. In this paper we bridge SN and SNR studies by using the output of a SN simulation as the input of a SNR simulation carried on until 500~yr. We consider the case of a thermonuclear explosion of a carbon-oxygen white dwarf star as a model for a Type Ia SN; specifically we use the N100 delayed detonation model of Seitenzahl et al 2013. In order to analyze the morphology of the SNR, we locate the three discontinuities that delineate the shell of shocked matter: the forward shock, the contact discontinuity, and the reverse shock, and we decompose their radial variations as a function of angular scale and time. Assuming a uniform ambient medium, we find that the impact of the SN on the SNR may still be visible after hundreds of years. Previous 3D simulations aiming at reproducing Tycho's SNR, that started out from spherically symmetric initial conditions, failed to reproduce structures at the largest angular scales observed in X-rays. Our new simulations strongly suggest that the missing ingredient was the initial asymmetries from the SN itself. With this work we establish a way of assessing the viability of SN models based on the resulting morphology of the SNR.

(Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ratio measurements and the white dwarf mass-radius relation. Tracer particles from 3d hydrodynamic simulations were post-processed with a large nuclear reaction network in order to determine the complete compositional state of the bound ONeFe remnant and the ejecta, and population synthesis simulations were performed in order to estimate the ECSN rate with respect to the CCSN rate. The 3d deflagration simulations drastically overproduce the neutron-rich isotopes $^{48}$Ca, $^{50}$Ti, $^{54}$Cr, $^{60}$Fe and several of the Zn isotopes relative to their solar abundances. Using the solar abundance distribution as our constraint, we place an upper limit on the frequency of thermonuclear ECSNe as 1$-$3~\% the frequency at which core-collapse supernovae (FeCCSNe) occur. This is on par with or 1~dex lower than the estimates for ECSNe from single stars. The upper limit from the yields is also in relatively good agreement with the predictions from our population synthesis simulations. The $^{54}$Cr/$^{52}$Cr and $^{50}$Ti/$^{48}$Ti isotopic ratios in the ejecta are a near-perfect match with recent measurements of extreme pre-solar meteoritc oxide grains, and $^{53}$Cr/$^{52}$Cr can also be matched if the ejecta condenses before mixing with the interstellar medium. Theoretical mass-radius relations for the bound ONeFe WD remnants of these explosions are apparently consistent with several observational WD candidates.

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.

The observed sub-class of "superluminous" Type Ia supernovae lacks a convincing theoretical explanation. If the emission of such objects were powered exclusively by radioactive decay of 56Ni formed in the explosion, a progenitor mass close to or even above the Chandrasekhar limit for a non-rotating white dwarf star would be required. Masses significantly exceeding this limit can be supported by differential rotation. We, therefore, explore explosions and predict observables for various scenarios resulting from differentially rotating carbon-oxygen white dwarfs close to their respective limit of stability. Specifically, we have investigated a prompt detonation model, detonations following an initial deflagration phase ("delayed detonation" models), and a pure deflagration model. In postprocessing steps, we performed nucleosynthesis and three-dimensional radiative transfer calculations, that allow us, for the first time, to consistently derive synthetic observables from our models. We find that all explosion scenarios involving detonations produce very bright events. The observables predicted for them, however, are inconsistent with any known subclass of Type Ia supernovae. Pure deflagrations resemble 2002cx-like supernovae and may contribute to this class. We discuss implications of our findings for the explosion mechanism and for the existence of differentially rotating white dwarfs as supernova progenitors.

We give an overview of recent efforts to model Type Ia supernovae and related astrophysical transients resulting from thermonuclear explosions in white dwarfs. In particular we point out the challenges resulting from the multi-physics multi-scale nature of the problem and discuss possible numerical approaches to meet them in hydrodynamical explosion simulations and radiative transfer modeling. We give examples of how these methods are applied to several explosion scenarios that have been proposed to explain distinct subsets or, in some cases, the majority of the observed events. In case we comment on some of the successes and shortcoming of these scenarios and highlight important outstanding issues.

In studies on intermediate- and high-resolution spectra of Type Ia supernovae (SNe Ia), some objects exhibit narrow Na-I D absorptions often blueshifted with respect to the rest wavelength within the host galaxy. The absence of these in other SNe Ia may reflect that the explosions have different progenitors: blueshifted Na-I D features might be explained by the outflows of 'single-degenerate' systems (binaries of a white dwarf with a non-degenerate companion). In this work, we search for systematic differences among SNe Ia for which the Na-I D characteristics have been clearly established in previous studies. We perform an analysis of the chemical abundances in the outer ejecta of 13 'spectroscopically normal' SNe Ia (five of which show blueshifted Na lines), modelling time series of photospheric spectra with a radiative-transfer code. We find only moderate differences between 'blueshifted-Na', 'redshifted-Na' and 'no-Na' SNe Ia, so that we can neither conclusively confirm a 'one-scenario' nor a 'two-scenario' theory for normal SNe Ia. Yet, some of the trends we see should be further studied using larger observed samples: Models for blueshifted-Na SNe tend to show higher photospheric velocities than no-Na SNe, corresponding to a higher opacity of the envelope. Consistently, blueshifted-Na SNe show hints of a somewhat larger iron-group content in the outer layers with respect to the no-Na subsample (and also to the redshifted-Na subsample). This agrees with earlier work where it was found that the light curves of no-Na SNe - often appearing in elliptical galaxies - are narrower, i.e. decline more rapidly.

Despite the importance of Type Ia supernovae (SNe Ia) for modern astrophysics, their detailed mechanism is still not fully understood. In this contribution, we present recent findings from numerical explosion models in the context of the observed diversity of SNe Ia and we discuss how these models can help to shed light on the explosion mechanism and the progenitor stars of SNe Ia. In addition, we introduce the Heidelberg Supernova Model Archive (HESMA), a new online data base where we provide integrated isotopic abundances and radially averaged ejecta profiles and synthetic observables for a wide range of state-of-the-art explosion models.

Despite their use as cosmological distance indicators and their importance in the chemical evolution of galaxies, the unequivocal identification of the progenitor systems and explosion mechanism of normal type Ia supernova (SN Ia) remains elusive. The leading hypothesis is that such a supernova is a thermonuclear explosion of a carbon-oxygen white dwarf, but the exact explosion mechanism is still a matter of debate. Observation of a galactic SN Ia would be of immense value in answering the many open questions related to these events. One potentially useful source of information about the explosion mechanism and progenitor is the neutrino signal. In this paper we compute the expected neutrino signal from a gravitationally confined detonation (GCD) explosion scenario for a SN~Ia and show how the flux at Earth contains features in time and energy unique to this scenario. We then calculate the expected event rates in the Super-K, Hyper-K, JUNO, DUNE, and IceCube detectors and find both Hyper-K and IceCube would see a few events for a GCD supernova at 1 kpc or closer, while Super-K, JUNO, and DUNE would see a events if the supernova were closer than ${\sim}0.3$ kpc. The distance and detector criteria needed to resolve the time and spectral features arising from the explosion mechanism, neutrino production, and neutrino oscillation processes are also discussed. The neutrino signal from the GCD is then compared with the signal from a deflagration-to-detonation transition (DDT) explosion model computed previously. We find the overall event rate is the most discriminating feature between the two scenarios followed by the event rate time structure. Using the event rate in the Hyper-K detector alone, the DDT can be distinguished from the GCD at 2$\sigma$ if the distance to the supernova is less than $2.3\;{\rm kpc}$ for a normal mass ordering and $3.6\;{\rm kpc}$ for an inverted ordering.

The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of the diffusion approximation to dynamical shear in hydrostatic stellar evolution models by comparing 1D models to a first-principle hydrodynamics simulation starting from the same initial conditions. We chose an initial model calculated with the stellar evolution code GENEC that is just at the onset of a dynamical shear instability but does not show any other instabilities (e.g., convection). This was mapped to the hydrodynamics code SLH to perform a 2D simulation in the equatorial plane. We compare the resulting profiles in the two codes and compute an effective diffusion coefficient for the hydro simulation. Shear instabilities develop in the 2D simulation in the regions predicted by linear theory to become unstable in the 1D model. Angular velocity and chemical composition is redistributed in the unstable region, thereby creating new unstable regions. Eventually the 2D simulation settles in a symmetric, steady state, which is Richardson stable everywhere, whereas the instability remains for longer in the 1D model due to current limitations in the 1D code. A spatially resolved diffusion coefficient is extracted by comparing the initial and final profiles of mean atomic mass. The presented simulation gives a first insight on hydrodynamics of shear instabilities in a real stellar environment and even allows us to directly extract an effective diffusion coefficient. We see evidence for a critical Richardson number of 0.25 as regions above this value remain stable for the course of the simulation.

Type Ia supernovae are associated with thermonuclear explosions of white dwarf stars. Combustion processes convert material in nuclear reactions and release the energy required to explode the stars. At the same time, they produce the radioactive species that power radiation and give rise to the formation of the observables. Therefore, the physical mechanism of the combustion processes, as reviewed here, is the key to understand these astrophysical events. Theory establishes two distinct modes of propagation for combustion fronts: subsonic deflagrations and supersonic detonations. Both are assumed to play an important role in thermonuclear supernovae. The physical nature and theoretical models of deflagrations and detonations are discussed together with numerical implementations. A particular challenge arises due to the wide range of spatial scales involved in these phenomena. Neither the combustion waves nor their interaction with fluid flow and instabilities can be directly resolved in simulations. Substantial modeling effort is required to consistently capture such effects and the corresponding techniques are discussed in detail. They form the basis of modern multidimensional hydrodynamical simulations of thermonuclear supernova explosions. The problem of deflagration-to-detonation transitions in thermonuclear supernova explosions is briefly mentioned.

Over the past years type Ia supernovae (SNe Ia) have become a major tool to determine the expansion history of the Universe, and considerable attention has been given to, both, observations and models of these events. However, until now, their progenitors are not known. The observed diversity of light curves and spectra seems to point at different progenitor channels and explosion mechanisms. Here, we present a new way to compare model predictions with observations in a systematic way. Our method is based on the construction of a metric space for SN Ia spectra by means of linear Principal Component Analysis (PCA), taking care of missing and/or noisy data, and making use of Partial Least Square regression (PLS) to find correlations between spectral properties and photometric data. We investigate realizations of the three major classes of explosion models that are presently discussed: delayed-detonation Chandrasekhar-mass explosions, sub-Chandrasekhar-mass detonations, and double-degenerate mergers, and compare them with data. We show that in the PC space all scenarios have observed counterparts, supporting the idea that different progenitors are likely. However, all classes of models face problems in reproducing the observed correlations between spectral properties and light curves and colors. Possible reasons are briefly discussed.

Based on the Roe solver a new technique that allows to correctly represent low Mach number flows with a discretization of the compressible Euler equations was proposed in Miczek et al.: New numerical solver for flows at various mach numbers. A&A 576, A50 (2015). We analyze properties of this scheme and demonstrate that its limit yields a discretization of the continuous limit system. Furthermore we perform a linear stability analysis for the case of explicit time integration and study the performance of the scheme under implicit time integration via the evolution of its condition number. A numerical implementation demonstrates the capabilities of the scheme on the example of the Gresho vortex which can be accurately followed down to Mach numbers of ~1e-10 .

Hydrodynamical simulations of stellar interactions require stable models of stars as initial conditions. Such initial models, however, are difficult to construct for giant stars because of the wide range in spatial scales of the hydrostatic equilibrium and in dynamical timescales between the core and the envelope of the giant. They are needed for, e.g., modeling the common envelope phase where a giant envelope encompasses both the giant core and a companion star. Here, we present a new method of approximating and reconstructing giant profiles from a stellar evolution code to produce stable models for multi-dimensional hydrodynamical simulations. We determine typical stellar stratification profiles with the 1D stellar evolution code MESA. After an appropriate mapping, hydrodynamical simulations are conducted using the moving-mesh code AREPO. The giant profiles are approximated by replacing the core of the giant with a point mass and by constructing a suitable continuation of the profile to the center. Different reconstruction methods are tested that can specifically control the convective behaviour of the model. After mapping to a grid, a relaxation procedure that includes damping of spurious velocities yields stable models in three-dimensional hydrodynamical simulations. Initially convectively stable configurations lead to stable hydrodynamical models while for stratifications that are convectively unstable in the stellar evolution code, simulations recover the convective behaviour of the initial model and show large convective plumes with Mach numbers up to 0.8. Examples are shown for a $2M_\odot$ red giant and a $0.67M_\odot$ asymptotic giant branch star. A detailed analysis shows that the improved method reliably provides stable models of giant envelopes that can be used as initial conditions for subsequent hydrodynamical simulations of stellar interactions involving giant stars.

During the common envelope (CE) phase, a giant star in a binary system overflows its Roche lobe and unstable mass transfer leads to a spiral-in of the companion, resulting in a close binary system or in a merger of the stellar cores. Dynamo processes during the CE phase have been proposed as a mechanism to generate magnetic fields that are important for forming magnetic white dwarfs (MWDs) and for shaping planetary nebulae. Here, we present the first magnetohydrodynamics simulations of the dynamical spiral-in during a CE phase. We find that magnetic fields are strongly amplified in the accretion stream around the $1M_\odot$ companion as it spirals into the envelope of a $2M_\odot$ RG. This leads to field strengths of 10 to 100 kG throughout the envelope after 120 d. The magnetic field amplification is consistent with being driven by the magnetorotational instability. The field strengths reached in our simulation make the magnetic field interesting for diagnostic purposes, but they are dynamically irrelevant. They are also too small to explain the formation of the highest fields found in MWDs, but may be relevant for luminous red novae, and detecting magnetic fields in these events would support the scenario as proposed here.

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 {\lambda}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 {\lambda}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.

In the classical picture, electron-capture supernovae and the accretion-induced collapse of oxygen-neon white dwarfs (ONeWDs) undergo an oxygen deflagration phase before gravitational collapse produces a neutron star (NS). These types of core collapse events are postulated to explain several astronomical phenomena. In this work, the deflagration phase is simulated for the first time using multidimensional hydrodynamics, with the aim of gaining new insight into the explosive deaths of $8-10~M_\odot$ stars and ONeWDs that accrete material from a binary companion star. The main aim is to determine whether these events are thermonuclear or core-collapse supernova explosions, and hence whether NSs are formed by such phenomena. The deflagration is simulated in ONe cores with three different central ignition densities. The intermediate density case is perhaps the most realistic, being based on recent nuclear physics calculations and 1D stellar models. The 3D hydrodynamic simulations presented in this work begin from a centrally confined flame structure using a level-set-based flame approach and are performed in $256^3$ and $512^3$ numerical resolutions. In the simulations with intermediate and low ignition density, the cores do not appear to collapse into NSs. Instead, almost a solar mass of material becomes unbound from the cores, leaving bound remnants. These simulations represent the case in which semiconvective mixing during the electron-capture phase preceding the deflagration is inefficient. The masses of the bound remnants double when Coulomb corrections are included in the EoS, however they still do not exceed the effective Chandrasekhar mass and, hence, would not collapse into NSs. The simulations with the highest ignition density ($\log_{10}\rho_{\rm c}=10.3$), representing the case where semiconvective mixing is very efficient, show clear signs that the core will collapse into a NS.

The gravitationally confined detonation (GCD) model has been proposed as a possible explosion mechanism for Type Ia supernovae in the single-degenerate evolution channel. Driven by buoyancy, a deflagration flame rises in a narrow cone towards the surface. For the most part, the flow of the expanding ashes remains radial, but upon reaching the outer, low-pressure layers of the white dwarf, an additional lateral component develops. This makes the deflagration ashes converge again at the opposite side, where the compression heats fuel and a detonation may be launched. To test the GCD explosion model, we perform a 3D simulation for a model with an ignition spot offset near the upper limit of what is still justifiable, 200 km. This simulation meets our deliberately optimistic detonation criteria and we initiate a detonation. The detonation burns through the white dwarf and leads to its complete disruption. We determine nucleosynthetic yields by post-processing 10^6 tracer particles with a 384 nuclide reaction network and we present multi-band light curves and time-dependent optical spectra. We find that our synthetic observables show a prominent viewing-angle sensitivity in UV and blue bands, which is in tension with observed SNe Ia. The strong dependence on viewing-angle is caused by the asymmetric distribution of the deflagration ashes in the outer ejecta layers. Finally, we perform a comparison of our model to SN 1991T. The overall flux-level of the model is slightly too low and the model predicts pre-maximum light spectral features due to Ca, S, and Si that are too strong. Furthermore, the model chemical abundance stratification qualitatively disagrees with recent abundance tomography results in two key areas: our model lacks low velocity stable Fe and instead has copious amounts of high-velocity 56Ni and stable Fe. We therefore do not find good agreement of the model with SN 1991T.

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.

The hard X-ray detector (HXD) onboard {\it Suzaku} measured soft $\gamma$-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 \sigma$), the $3 \sigma$ 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 $\gamma$-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 \sigma$ error on the 170-250 keV X-ray flux of $(1.2\pm0.7) \times10^{-4}$ ph s$^{-1}$ cm$^{-2}$.

The progenitor system(s) and the explosion mechanism(s) of Type Ia supernovae (SNe Ia) are still under debate. Non-electromagnetic observables, in particular gravitational waves and neutrino emission, of thermonuclear supernovae are a complementary window to light curves and spectra for studying these enigmatic objects. A leading model for SNe Ia is the thermonuclear incineration of a near-Chandrasekhar mass carbon-oxygen white dwarf star in a "delayed-detonation". We calculate a three-dimensional hydrodynamic explosion for the N100 delayed-detonation model extensively discussed in the literature, taking the dynamical effects of neutrino emission from all important contributing source terms into account. Although neutrinos carry away $2 \times 10^{49}$ erg of energy, we confirm the common view that neutrino energy losses are dynamically not very important, resulting in only a modest reduction of the final kinetic energy by two per cent. We then calculate the gravitational wave signal from the time evolution of the quadrupole moment. Our model radiates $7 \times 10^{39}$ erg in gravitational waves and the spectrum has a pronounced peak around 0.4 Hz. Depending on viewing angle and polarization, we find that the future space-based gravitational wave missions DECIGO and BBO would be able to detect our source to a distance of 1.3 Mpc. We predict a clear signature of the deflagration-to-detonation transition in the neutrino and the gravitational wave signals. If observed, such a feature would be a strong indicator of the realization of delayed-detonations in near-Chandrasekhar mass white dwarfs.

The common envelope (CE) phase is an important stage in binary stellar evolution. It is needed to explain many close binary stellar systems, such as cataclysmic variables, Type Ia supernova progenitors, or X-ray binaries. To form the resulting close binary, the initial orbit has to shrink, thereby transferring energy to the primary giant's envelope that is hence ejected. The details of this interaction, however, are still not understood. Here, we present new hydrodynamic simulations of the dynamical spiral-in forming a CE system. We apply the moving-mesh code AREPO to follow the interaction of a $1M_\odot$ compact star with a $2M_\odot$ red giant possessing a $0.4M_\odot$ core. The nearly Lagrangian scheme combines advantages of smoothed particle hydrodynamics and traditional grid-based hydrodynamic codes and allows us to capture also small flow features at high spatial resolution. Our simulations reproduce the initial transfer of energy and angular momentum from the binary core to the envelope by spiral shocks seen in previous studies, but after about 20 orbits a new phenomenon is observed. Large-scale flow instabilities are triggered by shear flows between adjacent shock layers. These indicate the onset of turbulent convection in the common envelope, thus altering the transport of energy on longer time scales. At the end of our simulation, only 8% of the envelope mass is ejected. The failure to unbind the envelope completely may be caused by processes on thermal time scales or unresolved microphysics.

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.

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.

The progenitors of many CCSNe are expected to be in binary systems. After the SN explosion, the companion may suffer from mass stripping and be shock heated as a result of the impact of the SN ejecta. If the binary system is disrupted, the companion is ejected as a runaway and hypervelocity star. By performing a series of 3D hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their companions. We use the BEC code to construct the detailed companion structure at the time of SN explosion. The impact of the SN blast wave on the companion is followed by means of 3D SPH simulations using the Stellar GADGET code. For main-sequence (MS) companions, we find that the amount of removed mass, impact velocity, and chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent with the results obtained from impact simulations of SNe~Ia. However, we find that the impact velocity is sensitive to the momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modeling of the ejecta. Because most companions to Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the strongly decaying impact effects with increasing binary separation, we argue that the majority of these SNe lead to inefficient mass stripping and shock heating of the companion star following the impact of the ejecta. Our simulations show that the impact effects of Ib/c SN ejecta on the structure of MS companions, and thus their long-term post-explosion evolution, is in general not dramatic. We find that at most 10% of their mass is lost, and their resulting impact velocities are less than 100 km/s.

The progenitor problem of Type Ia supernovae (SNe Ia) is still unsolved. Most of these events are thought to be explosions of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion scenarios, particularly those involving the externally triggered detonation of a sub-Chandrasekhar mass WD (sub-M Ch WD), there is also a possibility of having an oxygen-neon (ONe) WD as progenitor. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with detonations in CO WDs of similar mass and observational data of SNe Ia. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses below the Chandrasekhar mass (M Ch), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction network. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe Ia. We also perform binary evolution calculations to determine the number of SNe Ia involving ONe WDs relative to the number of other promising progenitor channels. The ejecta structures of our simulated detonations in sub-M Ch ONe WDs are similar to those from CO WDs. There are, however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are systematically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO models. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (1.2 Msun) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.

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.

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.

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

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.

Binary stellar evolution calculations predict that Chandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radially varying profile for the composition with a carbon depleted core. Many recent multi-dimensional simulations of Type Ia supernovae (SNe Ia), however, assume the progenitor WD has a homogeneous chemical composition. In this work, we explore the impact of different initial carbon profiles of the progenitor WD on the explosion phase and on synthetic observables in the Chandrasekhar-mass delayed detonation model. Spectra and light curves are compared to observations to judge the validity of the model. The explosion phase is simulated using the finite volume supernova code LEAFS, which is extended to treat different compositions of the progenitor WD. The synthetic observables are computed with the Monte Carlo radiative transfer code ARTIS. Differences in binding energies of carbon and oxygen lead to a lower nuclear energy release for carbon depleted material; thus, the burning fronts that develop are weaker and the total nuclear energy release is smaller. For otherwise identical conditions, carbon depleted models produce less Ni-56. Comparing different models with similar Ni-56 yields shows lower kinetic energies in the ejecta for carbon depleted models, but only small differences in velocity distributions and line velocities in spectra. The light curve width-luminosity relation (WLR) obtained for models with differing carbon depletion is roughly perpendicular to the observed WLR, hence the carbon mass fraction is probably only a secondary parameter in the family of SNe Ia.

The nucleosynthesis of proton-rich isotopes is calculated for multi-dimensional Chandrasekhar-mass models of Type Ia supernovae with different metallicities. The predicted abundances of the short-lived radioactive isotopes 92Nb, 97Tc, 98Tc and 146Sm are given in this framework. The abundance seeds are obtained by calculating s-process nucleosynthesis in the material accreted onto a carbon-oxygen white dwarf from a binary companion. A fine grid of s-seeds at different metallicities and 13C-pocket efficiencies is considered. A galactic chemical evolution model is used to predict the contribution of SNIa to the solar system p-nuclei composition measured in meteorites. Nuclear physics uncertainties are critical to determine the role of SNeIa in the production of 92Nb and 146Sm. We find that, if standard Chandrasekhar-mass SNeIa are at least 50% of all SNIa, they are strong candidates for reproducing the radiogenic p-process signature observed in meteorites.

Type-Ia supernovae result from binary systems that include a carbon-oxygen white dwarf, and these thermonuclear explosions typically produce 0.5 M_solar of radioactive 56Ni. The 56Ni is commonly believed to be buried deeply in the expanding supernova cloud. Surprisingly, in SN2014J we detected the lines at 158 and 812 keV from 56Ni decay ({\tau}~8.8 days) earlier than the expected several-week time scale, only ~20 days after the explosion, and with flux levels corresponding to roughly 10% of the total expected amount of 56Ni. Some mechanism must break the spherical symmetry of the supernova, and at the same time create a major amount of 56Ni at the outskirts. A plausible explanation is that a belt of helium from the companion star is accreted by the white dwarf, where this material explodes and then triggers the supernova event.

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.

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.

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.

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.

Context: Manganese is predominantly synthesised in Type Ia supernova (SN Ia) explosions. Owing to the entropy dependence of the Mn yield in explosive thermonuclear burning, SNe Ia involving near Chandrasekhar-mass white dwarfs (WDs) are predicted to produce Mn to Fe ratios significantly exceeding those of SN Ia explosions involving sub-Chandrasekhar mass primary WDs. Of all current supernova explosion models, only SN Ia models involving near-Chandrasekhar mass WDs produce [Mn/Fe] > 0.0. Aims: Using the specific yields for competing SN Ia scenarios, we aim to constrain the relative fractions of exploding near-Chandrasekhar mass to sub-Chandrasekhar mass primary WDs in the Galaxy. Methods: We extract the Mn yields from three-dimensional thermonuclear supernova simulations referring to different initial setups and progenitor channels. We then compute the chemical evolution of Mn in the Solar neighborhood, assuming SNe Ia are made up of different relative fractions of the considered explosion models. Results: We find that due to the entropy dependence of freeze-out yields from nuclear statistical equilibrium, [Mn/Fe] strongly depends on the mass of the exploding WD, with near-Chandraskher mass WDs producing substantially higher [Mn/Fe] than sub-Chandrasekhar mass WDs. Of all nucleosynthetic sources potentially influencing the chemical evolution of Mn, only explosion models involving the thermonuclear incineration of near-Chandrasekhar mass WDs predict solar or super-solar [Mn/Fe]. Consequently, we find in our chemical evolution calculations that the observed [Mn/Fe] in the Solar neighborhood at [Fe/H] > 0.0 cannot be reproduced without near-Chandrasekhar mass SN Ia primaries. Assuming that 50 per cent of all SNe Ia stem from explosive thermonuclear burning in near-Chandrasekhar mass WDs results in a good match to data.

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.

In the version of the SD scenario of SNe Ia studied here, a CO WD explodes close to the Chandrasekhar limit after accreting material from a non-degenerate He companion. In the present study, we employ the Stellar GADGET code to perform 3D hydrodynamical simulations of the interaction of the SN Ia ejecta with the He companion taking into account its orbital motion and spin. It is found that only 2%--5% of the initial companion mass are stripped off from the outer layers of He companions due to the SN impact. The dependence of the unbound mass (or the kick velocity) on the orbital separation can be fitted in good approximation by a power law for a given companion model. After the SN impact, the outer layers of a He donor star are significantly enriched with heavy elements from the low-expansion-velocity tail of SN Ia ejecta. The total mass of accumulated SN-ejecta material on the companion surface reaches about > 10e-3 M_sun for different companion models. This enrichment with heavy elements provides a potential way to observationally identify the surviving companion star in SN remnants. Finally, by artificially adjusting the explosion energy of the W7 explosion model, we find that the total accumulation of SN ejecta on the companion surface is also dependent on the explosion energy with a power law relation in good approximation.

The impact of nuclear physics uncertainties on nucleosynthesis in thermonuclear supernovae has not been fully explored using comprehensive and systematic studies with multiple models. To better constrain predictions of yields from these phenomena, we have performed a sensitivity study by post-processing thermodynamic histories from two different hydrodynamic, Chandrasekhar-mass explosion models. We have individually varied all input reaction and, for the first time, weak interaction rates by a factor of ten and compared the yields in each case to yields using standard rates. Of the 2305 nuclear reactions in our network, we find that the rates of only 53 reactions affect the yield of any species with an abundance of at least 10^-8 M_sun by at least a factor of two, in either model. The rates of the 12C(a,g), 12C+12C, 20Ne(a,p), 20Ne(a,g) and 30Si(p,g) reactions are among those that modify the most yields when varied by a factor of ten. From the individual variation of 658 weak interaction rates in our network by a factor of ten, only the stellar 28Si(b+)28Al, 32S(b+)32P and 36Ar(b+)36Cl rates significantly affect the yields of species in a model. Additional tests reveal that reaction rate changes over temperatures T > 1.5 GK have the greatest impact, and that ratios of radionuclides that may be used as explosion diagnostics change by a factor of less than two from the variation of individual rates by a factor of 10. Nucleosynthesis in the two adopted models is relatively robust to variations in individual nuclear reaction and weak interaction rates. Laboratory measurements of a limited number of reactions would help to further constrain predictions. As well, we confirm the need for a consistent treatment for relevant stellar weak interaction rates since simultaneous variation of these rates (as opposed to individual variation) has a significant effect on yields in our models.

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.

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.

In the SD scenario of SNe Ia the companion survives the SN explosion and thus should be visible near the center of the SN remnant and may show some unusual features. A promising approach to test progenitor models of SNe Ia is to search for the companion in SNRs. Here we present the results of 3D hydrodynamics simulations of the interaction between the SN Ia blast wave and a MS companion taking into consideration its orbital motion and spin. The primary goal of this work is to investigate the rotation of surviving companions after SN Ia explosions in the WD+MS scenario. We use Eggleton's code including the optically thick accretion wind model to obtain realistic models of companions. The impact of the SN blast wave on these companions is followed in 3D hydrodynamic simulations employing the SPH code GADGET3. We find that the rotation of the companion does not significantly affect the amount of stripped mass and the kick velocity caused by the SN impact. However, in our simulations, the rotational velocity of the companion is significantly reduced to about 14% to 32% of its pre-explosion value due to the expansion of the companion and the fact that 55%-89% of the initial angular momentum is carried away by the stripped matter. Compared with the observed rotational velocity of the presumed companion star of Tycho's SN, Tycho G, of 6 km/s the final rotational velocity we obtain is still higher by at least a factor of two. Whether this difference is significant, and may cast doubts on the suggestion that Tycho G is the companion of SN 1572, has to be investigated in future studies. Based on binary population synthesis results we present, for the first time, the expected distribution of rotational velocities of companions after the explosion which may provide useful information for the identification of the surviving companion in observational searches in other historical SNRs.

We explore SNIa as p-process sources in the framework of two-dimensional SNIa models using enhanced s-seed distributions as directly obtained from a sequence of thermal pulse instabilities. The SNIa WD precursor is assumed to have reached the Chandrasekhar mass limit in a binary system by mass accretion from a giant/main sequence companion. We apply the tracer-particle method to reconstruct the nucleosynthesis from the thermal histories of Lagrangian particles, passively advected in the hydrodynamic calculations. For each particle we follow the explosive nucleosynthesis with a detailed nuclear reaction network. We select tracers within the typical temperature range for p-process production, 1.5-3.7 109K, and analyse in detail their behaviour, exploring the influence of different s-process distributions on the p-process nucleosynthesis. We find that SNIa contribute to a large fraction of p-nuclei, both the light p-nuclei and the heavy-p nuclei at a quite flat average production factor. For the first time, the very abundant Ru and Mo p-isotopes are reproduced at the same level as the heavy p-nuclei. We investigate the metallicity effect on the p-process production. Starting with a range of s-seeds distributions obtained for different metallicities, running SNIa two-dimensional models and using a simple chemical evolution code, we give estimates of the SNIa contribution to the solar p-process composition. We find that SNIa contribute for at least 50% at the solar p-nuclei composition, in a primary way.

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.

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

We present results for a suite of fourteen three-dimensional, high resolution hydrodynamical simulations of delayed-detonation modelsof Type Ia supernova (SN Ia) explosions. This model suite comprises the first set of three-dimensional SN Ia simulations with detailed isotopic yield information. As such, it may serve as a database for Chandrasekhar-mass delayed-detonation 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 the 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 x 10^9 gcc, plus in addition one high central density (5.5 x 10^9 gcc), and one low central density (1.0 x 10^9 gcc) rendition of the 100 ignition kernel configuration. For each simulation 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 56Ni masses from 0.32 to 1.11 solar masses. 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) in a shell surrounding a 56Ni-rich core. The models further predict relatively low velocity oxygen and carbon, with typical minimum velocities around 4,000 and 10,000 km/s, respectively.

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.

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.

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.

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.

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.

We present the first maximum-light ultraviolet (UV) through near-infrared (NIR) Type Ia supernova (SN Ia) spectrum. This spectrum of SN 2011iv was obtained nearly simultaneously by the Hubble Space Telescope at UV/optical wavelengths and the Magellan Baade telescope at NIR wavelengths. These data provide the opportunity to examine the entire maximum-light SN Ia spectral-energy distribution. Since the UV region of a SN Ia spectrum is extremely sensitive to the composition of the outer layers of the explosion, which are transparent at longer wavelengths, this unprecedented spectrum can provide strong constraints on the composition of the SN ejecta, and similarly the SN explosion and progenitor system. SN 2011iv is spectroscopically normal, but has a relatively fast decline (Delta m_15 (B) = 1.69 +/- 0.05 mag). We compare SN 2011iv to other SNe Ia with UV spectra near maximum light and examine trends between UV spectral properties, light-curve shape, and ejecta velocity. We tentatively find that SNe with similar light-curve shapes but different ejecta velocities have similar UV spectra, while those with similar ejecta velocities but different light-curve shapes have very different UV spectra. Through a comparison with explosion models, we find that both a solar-metallicity W7 and a zero-metallicity delayed-detonation model provide a reasonable fit to the spectrum of SN 2011iv from the UV to the NIR.

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.

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.

Based on the largest homogeneously reduced set of Type Ia supernova luminosity data currently available -- the Union2 sample -- we reconstruct the expansion history of the Universe in a model-independent approach. Our method tests the geometry of the Universe directly without reverting to any assumptions made on its energy content. This allows us to constrain Dark Energy models and non-standard cosmologies in a straightforward way. The applicability of the presented method is not restricted to testing cosmological models. It can be a valuable tool for pointing out systematic errors hidden in the supernova data and planning future Type Ia supernova cosmology campaigns.

We present three-dimensional numerical simulations of turbulent combustion converting a neutron star into a quark star. Hadronic matter, described by a micro-physical finite-temperature equation of state, is converted into strange quark matter. We assume this phase, represented by a bag-model equation of state, to be absolutely stable. Following the example of thermonuclear burning in white dwarfs leading to Type Ia supernovae, we treat the conversion process as a potentially turbulent deflagration. Solving the non-relativistic Euler equations using established numerical methods we conduct large eddy simulations including an elaborate subgrid scale model, while the propagation of the conversion front is modeled with a level-set method. Our results show that for large parts of the parameter space the conversion becomes turbulent and therefore significantly faster than in the laminar case. Despite assuming absolutely stable strange quark matter, in our hydrodynamic approximation an outer layer remains in the hadronic phase, because the conversion front stops when it reaches conditions under which the combustion is no longer exothermic.

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.

We explore SNe Ia as p-process sites in the framework of two-dimensional SN Ia delayed detonation and pure deflagration models. The WD precursor is assumed to have reached the Chandrasekhar mass in a binary system by mass accretion from a giant/main sequence companion. We use enhanced s-seed distributions, obtained from a sequence of thermal pulse instabilities both in the AGB phase and in the accreted material. We apply the tracer-particle method to reconstruct the nucleosynthesis by the thermal histories of Lagrangian particles, passively advected in the hydrodynamic calculations. For each particle we follow the explosive nucleosynthesis with a detailed network for all isotopes up to 209Bi. We find that SNe Ia can produce a large amount of p-nuclei, both the light p-nuclei below A=120 and the heavy-p nuclei, at quite flat average production factors, tightly related to the s-process seed distribution. For the first time, we find a stellar source able to produce both, light and heavy p-nuclei almost at the same level as 56Fe, including the very debated neutron magic 92,94Mo and 96,98Ru. We also find that there is an important contribution from p-process nucleosynthesis to the s-only nuclei 80Kr, 86Sr, to the neutron magic 90Zr, and to the neutron-rich 96Zr. Finally, we investigate the metallicity effect on p-process. Starting with different s-process seed distributions, for two metallicities Z = 0.02 and Z = 0.001, running SNe Ia models with different initial composition, we estimate that SNe Ia can contribute to, at least, 50% of the solar p-process composition.

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.

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)

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.

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.

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.

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 {\gamma}-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.

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.

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.

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

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.

The theory of radiative transfer provides the link between the physical conditions in an astrophysical object and the observable radiation which it emits. Thus accurately modelling radiative transfer is often a necessary part of testing theoretical models by comparison with observations. We describe a new radiative transfer code which employs Monte Carlo methods for the numerical simulation of radiation transport in expanding media. We discuss the application of this code to the calculation of synthetic spectra and light curves for a Type Ia supernova explosion model and describe the sensitivity of the results to certain approximations made in the simulations.

Type Ia supernovae (SNe Ia) are thought to result from thermonuclear explosions of carbon-oxygen white dwarf stars. Existing models generally explain the observed properties, with the exception of the sub-luminous 1991-bg-like supernovae. It has long been suspected that the merger of two white dwarfs could give rise to a type Ia event, but hitherto simulations have failed to produce an explosion. Here we report a simulation of the merger of two equal-mass white dwarfs that leads to an underluminous explosion, though at the expense of requiring a single common-envelope phase, and component masses of ~0.9 M_sun. The light curve is too broad, but the synthesized spectra, red colour and low expansion velocities are all close to what is observed for sub-luminous 1991bg-like events. While mass ratios can be slightly less than one and still produce an underluminous event, the masses have to be in the range 0.83-0.9 M_sun.

Type Ia supernovae result when carbon-oxygen white dwarfs in binary systems accrete mass from companion stars, reach a critical mass, and explode. The near uniformity of their light curves makes these supernovae good standard candles for measuring cosmic expansion, but a correction must be applied to account for the fact that the brighter supernovae have broader light curves. One-dimensional modelling, with a certain choice of parameters, can reproduce this general trend in the width-luminosity relation, but the processes of ignition and detonation have recently been shown to be intrinsically asymmetric. Here we report on multi-dimensional modelling of the explosion physics and radiative transfer that reveals that the breaking of spherical symmetry is a critical factor in determining both the width luminosity relation and the observed scatter about it. The deviation from sphericity can also explain the finite polarization detected in the light from some supernovae. The slope and normalization of the width-luminosity relation has a weak dependence on certain properties of the white dwarf progenitor, in particular the trace abundances of elements other than carbon and oxygen. Failing to correct for this effect could lead to systematic overestimates of up to 2% in the distance to remote supernovae.

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.

The nature of carbon burning flames in Type Ia supernovae is explored as they interact with Kolmogorov turbulence. One-dimensional calculations using the Linear Eddy Model of Kerstein (1991) elucidate three regimes of turbulent burning. In the simplest case, large scale turbulence folds and deforms thin laminar flamelets to produce a flame brush with a total burning rate given approximately by the speed of turbulent fluctuations on the integral scale, U_L. This is the regime where the supernova explosion begins and where most of its pre-detonation burning occurs. As the density declines, turbulence starts to tear the individual flamelets, making broader structures that move faster. For a brief time, these turbulent flamelets are still narrow compared to their spacing and the concept of a flame brush moving with an overall speed of U_L remains valid. However, the typical width of the individual flamelets, which is given by the condition that their turnover time equals their burning time, continues to increase as the density declines. Eventually, mixed regions almost as large as the integral scale itself are transiently formed. At that point, a transition to detonation can occur. The conditions for such a transition are explored numerically and it is estimated that the transition will occur for densities near 1 x 10**7 g/cm**3, provided the turbulent speed on the integral scale exceeds about 15% sonic. An example calculation shows the details of a detonation actually developing.

The nature of Type Ia supernova progenitors is still unclear. The outstanding characteristic of the single-degenerate scenario is that it contains hydrogen in the binary companion of the exploding white dwarf star, which, if mixed into the ejecta of the supernova in large amounts may lead to conflicts with the observations thus ruling out the scenario. We investigate the effect of the impact of Type Ia supernova ejecta on a main sequence companion star of the progenitor system. With a series of simulations we investigate how different parameters of this system affect the amount of hydrogen stripped from the companion by the impact. The stellar evolution code GARSTEC is used to set up the structure of the companion stars mimicking the effect of a binary evolution phase. The impact itself is simulated with the smoothed particle hydrodynamics code GADGET2. We reproduce and confirm the results of earlier grid-based hydrodynamical simulation. Parameter studies of the progenitor system are extended to include the results of recent binary evolution studies. The more compact structure of the companion star found here significantly reduces the stripped hydrogen mass. The low hydrogen masses resulting from a more realistic companion structure are consistent with current observational constraints. Therefore, the single-degenerate scenario remains a valid possibility for Type Ia supernova progenitors. These new results are not a numerical effect, but the outcome of different initial conditions.

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.

The application of Type Ia supernovae (SNe Ia) as distance indicators in cosmology calls for a sound understanding of these objects. Recent years have seen a brisk development of astrophysical models which explain SNe Ia as thermonuclear explosions of white dwarf stars. While the evolution of the progenitor is still uncertain, the explosion mechanism certainly involves the propagation of a thermonuclear flame through the white dwarf star. Three-dimensional hydrodynamical simulations allowed to study a wide variety of possibilities involving subsonic flame propagation (deflagrations), flames accelerated by turbulence, and supersonic detonations. These possibilities lead to a variety of scenarios. I review the currently discussed approaches and present some recent results from simulations of the turbulent deflagration model and the delayed detonation model.

Vastly different time and length scales are a common problem in numerical simulations of astrophysical phenomena. Here, we present an approach to numerical modeling of such objects on the example of Type Ia supernova simulations. The evolution towards the explosion proceeds on much longer time scales than the explosion process itself. The physical length scales relevant in the explosion process cover 11 orders of magnitude and turbulent effects dominate the physical mechanism. Despite these challenges, three-dimensional simulations of Type Ia supernova explosions have recently become possible and pave the way to a better understanding of these important astrophysical objects.

In the "double-detonation sub-Chandrasekhar" model for type Ia supernovae, a carbon-oxygen (C + O) white dwarf accumulates sufficient amounts of helium such that a detonation ignites in that layer before the Chandrasekhar mass is reached. This detonation is thought to trigger a secondary detonation in the C + O core. By means of one- and two-dimensional hydrodynamic simulations, we investigate the robustness of this explosion mechanism for generic 1-M_sun models and analyze its observable predictions. Also a resolution dependence in numerical simulations is analyzed. The propagation of thermonuclear detonation fronts, both in helium and in the carbon-oxygen mixture, is computed by means of both a level-set function and a simplified description for nuclear reactions. The decision whether a secondary detonation is triggered in the white dwarf's core or not is made based on criteria given in the literature. In a parameter study involving different initial flame geometries for He-shell masses of 0.2 and 0.1 M_sun, we find that a secondary detonation ignition is a very robust process. Converging shock waves originating from the detonation in the He shell generate the conditions for a detonation near the center of the white dwarf in most of the cases considered. Finally, we follow the complete evolution of three selected models with 0.2 M_sun of He through the C/O-detonation phase and obtain nickel-masses of about 0.40 to 0.45 M_sun. Although we have not done a complete scan of the possible parameter space, our results show that sub-Chandrasekhar models are not good candidates for normal or sub-luminous type Ia supernovae. The chemical composition of the ejecta features significant amounts of nickel in the outer layers at high expansion velocities, which is inconsistent with near-maximum spectra. (abbreviated)

Although delayed detonation models of thermonuclear explosions of white dwarfs seem promising for reproducing Type Ia supernovae, the transition of the flame propagation mode from subsonic deflagration to supersonic detonation remains hypothetical. A potential instant for this transition to occur is the onset of the distributed burning regime, i.e. the moment when turbulence first affects the internal flame structure. Some studies of the burning microphysics indicate that a deflagration-to-detonation transition may be possible here, provided the turbulent intensities are strong enough. Consequently, the magnitude of turbulent velocity fluctuations generated by the deflagration flame is analyzed at the onset of the distributed burning regime in several three-dimensional simulations of deflagrations in thermonuclear supernovae. It is shown that the corresponding probability density functions fall off towards high turbulent velocity fluctuations much more slowly than a Gaussian distribution. Thus, values claimed to be necessary for triggering a detonation are likely to be found in sufficiently large patches of the flame. Although the microphysical evolution of the burning is not followed and a successful deflagration-to-detonation transition cannot be guaranteed from simulations presented here, the results still indicate that such events may be possible in Type Ia supernova explosions.

A simulation of the thermonuclear explosion of a Chandrasekhar-mass C+O white dwarf, the most popular scenario of a type Ia supernova (SN Ia), is presented. The underlying modeling is pursued in a self-consistent way, treating the combustion wave as a turbulent deflagration using well tested methods developed for laboratory combustion and based on the concept of `large eddy simulations' (LES). Such consistency requires to capture the onset of the turbulent cascade on resolved scales. This is achieved by computing the dynamical evolution on a 1024$^3$ moving grid, which resulted in the best-resolved three-dimensional SN Ia simulation carried out thus far, reaching the limits of what can be done on present supercomputers. Consequently, the model has no free parameters other than the initial conditions at the onset of the explosion, and therefore it has considerable predictive power. Our main objective is to determine to which extent such a simulation can account for the observations of normal SNe Ia. Guided by previous simulations with less resolution and a less sophisticated flame model, initial conditions were chosen that yield a reasonably strong explosion and a sufficient amount of radioactive nickel for a bright display. We show that observables are indeed matched to a reasonable degree. In particular, good agreement is found with the light curves of normal SNe Ia. Moreover, the model reproduces the general features of the abundance stratification as inferred from the analysis of spectra. This indicates that it captures the main features of the explosion mechanism of SNe Ia. However, we also show that even a seemingly best-choice pure deflagration model has shortcomings that indicate the need for a different mode of nuclear burning at late times, perhaps the transition to a detonation at low density.

Aims: We present the first full-star three-dimensional explosion simulations of thermonuclear supernovae including parameterized deflagration-to-detonation transitions that occur once the flame enters the distributed burning regime. Methods: Treating the propagation of both the deflagration and the detonation waves in a common front-tracking approach, the detonation is prevented from crossing ash regions. Results: Our criterion triggers the detonation wave at the outer edge of the deflagration flame and consequently it has to sweep around the complex structure and to compete with expansion. Despite the impeded detonation propagation, the obtained explosions show reasonable agreement with global quantities of observed type Ia supernovae. By igniting the flame in different numbers of kernels around the center of the exploding white dwarf, we set up three different models shifting the emphasis from the deflagration phase to the detonation phase. The resulting explosion energies and iron group element productions cover a large part of the diversity of type Ia supernovae. Conclusions: Flame-driven deflagration-to-detonation transitions, if hypothetical, remain a possibility deserving further investigation.

Type Ia supernovae, the thermonuclear explosions of white dwarf stars composed of carbon and oxygen, were instrumental as distance indicators in establishing the acceleration of the universe's expansion. However, the physics of the explosion are debated. Here we report a systematic spectral analysis of a large sample of well observed type Ia supernovae. Mapping the velocity distribution of the main products of nuclear burning, we constrain theoretical scenarios. We find that all supernovae have low-velocity cores of stable iron-group elements. Outside this core, nickel-56 dominates the supernova ejecta. The outer extent of the iron-group material depends on the amount of nickel-56 and coincides with the inner extent of silicon, the principal product of incomplete burning. The outer extent of the bulk of silicon is similar in all SNe, having an expansion velocity of ~11000 km/s and corresponding to a mass of slightly over one solar mass. This indicates that all the supernovae considered here burned similar masses, and suggests that their progenitors had the same mass. Synthetic light curve parameters and three-dimensional explosion simulations support this interpretation. A single explosion scenario, possibly a delayed detonation, may thus explain most type Ia supernovae.

Modeling type Ia supernova (SN Ia) explosions in three dimensions allows to eliminate any undetermined parameters and provides predictive power to simulations. This is necessary to improve the understanding of the explosion mechanism and to settle the question of the applicability of SNe Ia in cosmological distance measurements. Since the models contain no tunable parameters, it is also possible to directly assess their validity on the basis of a comparison with observations. Here, we describe the modeling of SNe Ia as thermonuclear explosions in which the flame after ignition near the center of the progenitor white dwarf star propagates outward in the sub-sonic deflagration mode accelerated by the interaction with turbulence. We explore the capabilities of this model by comparison with observations and show in a preliminary approach, how such a model can be applied to study the origin of the diversity of SNe Ia.

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.

Type Ia supernovae (SNe Ia) are one of the major tools to determine the cosmological parameters. Utilizing them as distance indicators, it is possible to geometrically survey the universe. To this end, the intrinsic scatter in the luminosities of these events needs to be calibrated using empirical relations between observables. A theoretical explanation for these relations is still lacking and can only be provided by a sound understanding of the mechanism of SNe Ia. Recently there has been significant progress in modeling SNe Ia. We report on numerical simulations of the explosion process, compare the results with observations of nearby SNe Ia, and discuss current uncertainties of the models. The presented simulations shed some light on the origin of the diversity of SNe Ia. Such simulations will pave the way towards an understanding of SN Ia diversities and correlations of their properties and ultimately provide a tool to validate the cosmological implications of SN Ia distance measurements.

We present results from a systematic study of the effects of initial parameters on three-dimensional thermonuclear supernova models.

We present a systematic study of the diversity of three-dimensional deflagration simulations of Type Ia supernova explosions arising from variations of the initial parameters. By altering the carbon mass fraction, the central density, and the metallicity of the progenitor white dwarf star, we set up a grid of numerical explosion models. While changing the central density has the largest impact on the explosion energy, the largest variation in the 56Ni production is found by changing the metallicity of the models. Varying the carbon mass fraction hardly affects the 56Ni synthesized although it alters the energetics of the explosion. Possible consequences for the shape of light curves of Type Ia supernovae are discussed.

We review the astrophysical modeling of type Ia supernova explosions and describe numerical methods to implement numerical simulations of these events. Some results of such simulations are discussed.

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.

We present synthetic bolometric and broad-band UBVRI light curves of SNe Ia, for four selected 3-D deflagration models of thermonuclear supernovae. The light curves are computed with the 1-D hydro code STELLA, which models (multi-group time-dependent) non-equilibrium radiative transfer inside SN ejecta. Angle-averaged results from 3-D hydrodynamical explosion simulations with the composition determined in a nucleosynthetic postprocessing step served as the input to the radiative transfer model. The predicted model UBV light curves do agree reasonably well with the observed ones for SNe Ia in the range of low to normal luminosities, although the underlying hydrodynamical explosion models produced only a modest amount of radioactive Ni56 and relatively low kinetic energy in the explosion. The evolution of predicted B and V fluxes in the model with a Ni56 mass of 0.42 M_sun follows the observed decline rate after the maximum very well, although the behavior of fluxes in other filters somewhat deviates from observations, and the bolometric decline rate is a bit slow. Using our models, we check the validity of Arnett's rule and the accuracy of the procedure for extracting the Ni56 mass from the observed light curves. We find that the comparison between theoretical light curves and observations provides a useful tool to validate SN Ia models. The steps necessary to improve the agreement between theory and observations are set out.

The dynamics of the explosive burning process is highly sensitive to the flame speed model in numerical simulations of type Ia supernovae. Based upon the hypothesis that the effective flame speed is determined by the unresolved turbulent velocity fluctuations, we employ a new subgrid scale model which includes a localised treatment of the energy transfer through the turbulence cascade in combination with semi-statistical closures for the dissipation and non-local transport of turbulence energy. In addition, subgrid scale buoyancy effects are included. In the limit of negligible energy transfer and transport, the dynamical model reduces to the Sharp-Wheeler relation. According to our findings, the Sharp-Wheeler relation is insuffcient to account for the complicated turbulent dynamics of flames in thermonuclear supernovae. The application of a co-moving grid technique enables us to achieve very high spatial resolution in the burning region. Turbulence is produced mostly at the flame surface and in the interior ash regions. Consequently, there is a pronounced anisotropy in the vicinity of the flame fronts. The localised subgrid scale model predicts significantly enhanced energy generation and less unburnt carbon and oxygen at low velocities compared to earlier simulations.

We present a systematic survey of the capabilities of type Ia supernova explosion models starting from a number of flame seeds distributed around the center of the white dwarf star. To this end we greatly improved the resolution of the numerical simulations in the initial stages. This novel numerical approach facilitates a detailed study of multi-spot ignition scenarios with up to hundreds of ignition sparks. Two-dimensional simulations are shown to be inappropriate to study the effects of initial flame configurations. Based on a set of three-dimensional models, we conclude that multi-spot ignition scenarios may improve type Ia supernova models towards better agreement with observations. The achievable effect reaches a maximum at a limited number of flame ignition kernels as shown by the numerical models and corroborated by a simple dimensional analysis.

The use of type Ia supernovae as distance indicators for cosmology has initiated a search for theoretical arguments supporting the empirical calibration methods applied. To this end, as a first step, a sound understanding of the origin of the observed diversity in type Ia supernova properties is needed. Here we present a first systematic study of effects resulting from changing some physical parameters of three-dimensional deflagration models of thermonuclear supernovae. In our study we vary the progenitor's carbon-to-oxygen ratio and its central density prior to ignition because both properties are not well determined by stellar evolution theory and they may change from supernova to supernova. Next we compute for these explosion models the nucleosynthesis yields in a post-processing step. This, in addition, allows us to study variations in the progenitor's metallicity by means of different $^{22}$Ne mass fractions in the initial composition. We find that the progenitor's carbon-to-oxygen ratio and its central density affect the energy release of the models and thus the expansion velocity of the supernova. Moreover, we find that changing the metallicity and the central density changes the production of radioactive $^{56}$Ni and thus affects the luminosity. In contrast, the carbon-to-oxygen ratio has little effect on the $^{56}$Ni production. Implications of the found variations of the explosion energy and the produced $^{56}$Ni mass for the type Ia supernova diversity are discussed.

Late time synthetic spectra of Type Ia supernovae, based on three-dimensional deflagration models, are presented. We mainly focus on one model,"c3_3d_256_10s", for which the hydrodynamics (Roepke 2005) and nucleosynthesis (Travaglio et al. 2004) was calculated up to the homologous phase of the explosion. Other models with different ignition conditions and different resolution are also briefly discussed. The synthetic spectra are compared to observed late time spectra. We find that while the model spectra after 300 to 500 days show a good agreement with the observed Fe II-III features, they also show too strong O I and C I lines compared to the observed late time spectra. The oxygen and carbon emission originates from the low-velocity unburned material in the central regions of these models. To get agreement between the models and observations we find that only a small mass of unburned material may be left in the center after the explosion. This may be a problem for pure deflagration models, although improved initial conditions, as well as higher resolution decrease the discrepancy. The relative intensity from the different ionization stages of iron is sensitive to the density of the emitting iron-rich material. We find that clumping, with the presence of low density regions, is needed to reproduce the observed iron emission, especially in the range between 4000 and 6000 AA. Both temperature and ionization depend sensitively on density, abundances and radioactive content. This work therefore illustrates the importance of including the inhomogeneous nature of realistic three-dimensional explosion models. We briefly discuss the implications of the spectral modeling for the nature of the explosion.

The deflagration mode of flame propagation in Type Ia supernova (SN Ia) models requires a correct description of the interaction of the flame with turbulent motions. It is well-known that turbulent combustion proceeds in different regimes. For most parts of the deflagration in SNe Ia the flamelet regime applies. This has been modeled in previous multi-dimensional simulations. However, at late stages of the explosion, the flame will propagate in the so-called distributed regime. We investigate the effects of this regime on SN Ia models in a first and simplified approach and show that the trend of effects seems capable of curing some problems of current pure deflagration models.

We present full-star simulations of Type Ia supernova explosions on the basis of the standard Chandrasekhar-mass deflagration model. Most simulations so far considered only one spatial octant and assumed mirror symmetry to the other octants. Two full-star models are evolved to homologous expansion and compared with previous single-octant simulations. Therefrom we analyze the effect of abolishing the artificial symmetry constraint on the evolution of the flame surface. It turns out that the development of asymmetries depends on the chosen initial flame configuration. Such asymmetries of the explosion process could possibly contribute to the observed polarization of some Type Ia supernova spectra.

The last years have witnessed a rapid development of three-dimensional models of Type Ia supernova explosions. Consequently, the next step is to evaluate these models under variation of the initial parameters and to compare them with observations. To calculate synthetic lightcurves and spectra from numerical models, it is mandatory to follow the evolution up to homologous expansion. We report on methods to achieve this in our current implementation of multi-dimensional Type Ia supernova explosion models. The novel scheme is thoroughly tested in two dimensions and a simple example of a three-dimensional simulation is presented. We discuss to what degree the assumption of homologous expansion is justified in these models.

Type Ia supernovae have recently received considerable attention because it appears that they can be used as "standard candles" to measure cosmic distances out to billions of light years away from us. Observations of type Ia supernovae seem to indicate that we are living in a universe that started to accelerate its expansion when it was about half its present age. These conclusions rest primarily on phenomenological models which, however, lack proper theoretical understanding, mainly because the explosion process, initiated by thermonuclear fusion of carbon and oxygen into heavier elements, is difficult to simulate even on supercomputers. Here, we investigate a new way of modeling turbulent thermonuclear deflagration fronts in white dwarfs undergoing a type Ia supernova explosion. Our approach is based on a level set method which treats the front as a mathematical discontinuity and allows for full coupling between the front geometry and the flow field. New results of the method applied to the problem of type Ia supernovae are obtained. It is shown that in 2-D with high spatial resolution and a physically motivated subgrid scale model for the nuclear flames numerically "converged" results can be obtained, but for most initial conditions the stars do not explode. In contrast, simulations in 3-D, do give the desired explosions and many of their properties, such as the explosion energies, lightcurves and nucleosynthesis products, are in very good agreement with observed 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.

One of the major challenges for theoretical modeling of Type Ia supernova explosions is to explain the diversity of these events and the empirically established correlation between their peak luminosity and light curve shape. In the framework of the so-called Chandrasekhar mass models, the progenitor's carbon-to-oxygen ratio has been suggested to be a principal source of peak luminosity variations due to a variation in the production of radioactive $^{56}$Ni during the explosion. The underlying idea is that an enhanced carbon mass fraction should result in a more vigorous explosion since here the energy release from nuclear reactions is increased. It was suspected that this would produce a higher amount of $^{56}$Ni in the ejecta. In this letter we describe a mechanism resulting from an interplay between nucleosynthesis and turbulent flame evolution which counteracts such an effect. Based on three-dimensional simulations we argue that it is nearly balanced and only minor differences in the amount of synthesized $^{56}$Ni with varying carbon mass fraction in the progenitor can be expected. Therefore this progenitor parameter is unlikely to account for the observed variations in Type Ia supernova luminosity. We discuss possible effects on the calibration of cosmological measurements.

We investigate the interaction of thermonuclear flames in Type Ia supernova explosions with vortical flows by means of numerical simulations. In our study, we focus on small scales, where the flame propagation is no longer dominated by the turbulent cascade originating from large-scale effects. Here, the flame propagation proceeds in the cellular burning regime, resulting from a balance between the Landau-Darrieus instability and its nonlinear stabilization. The interaction of a cellularly stabilized flame front with a vortical fuel flow is explored applying a variety of fuel densities and strengths of the velocity fluctuations. We find that the vortical flow can break up the cellular flame structure if it is sufficiently strong. In this case the flame structure adapts to the imprinted flow field. The transition from the cellularly stabilized front to the flame structure dominated by vortices of the flow proceeds in a smooth way. The implications of the results of our simulations for Type Ia Supernova explosion models are discussed.

We present a numerical investigation of the cellular burning regime in Type Ia supernova explosions. This regime holds at small scales (i.e. below the Gibson scale), which are unresolved in large-scale Type Ia supernova simulations. The fundamental effects that dominate the flame evolution here are the Landau-Darrieus instability and its nonlinear stabilization, leading to a stabilization of the flame in a cellular shape. The flame propagation into quiescent fuel is investigated addressing the dependence of the simulation results on the specific parameters of the numerical setup. Furthermore, we investigate the flame stability at a range of fuel densities. This is directly connected to the questions of active turbulent combustion (a mechanism of flame destabilization and subsequent self-turbulization) and a deflagration-to-detonation transition of the flame. In our simulations we find no substantial destabilization of the flame when propagating into quiescent fuels of densities down to ~10^7 g/cm^3, corroborating fundamental assumptions of large-scale SN Ia explosion models. For these models, however, we suggest an increased lower cutoff for the flame propagation velocity to take the cellular burning regime into account.

We present a numerical model which allows us to investigate thermonuclear flames in Type Ia supernova explosions. The model is based on a finite-volume explicit hydrodynamics solver employing PPM. Using the level-set technique combined with in-cell reconstruction and flux-splitting schemes we are able to describe the flame in the discontinuity approximation. We apply our implementation to flame propagation in Chandrasekhar-mass Type Ia supernova models. In particular we concentrate on intermediate scales between the flame width and the Gibson-scale, where the burning front is subject to the Landau-Darrieus instability. We are able to reproduce the theoretical prediction on the growth rates of perturbations in the linear regime and observe the stabilization of the flame in a cellular shape. The increase of the mean burning velocity due to the enlarged flame surface is measured. Results of our simulation are in agreement with semianalytical studies.

We present a numerical model to study the behavior of thermonuclear flames in the discontinuity approximation. This model is applied to investigate the Landau-Darrieus instability under conditions found in Type Ia supernova explosions of Chandrasekhar mass white dwarfs. This is a first step to explore the flame microstructure in these events. The model reproduces Landau's linearized stability analysis in early stages of the flame evolution and the stabilization in a cellular flame structure in the nonlinear stage.