Following multi-dimensional Type Ia supernova explosion models to homologous expansion

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.Comment: 15 pages, 16 figures, resolution of some figures reduced to meet astro-ph file size restriction, submitted to A&

Double-detonation supernovae of sub-Chandrasekhar mass white dwarfs

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)Comment: 11 pages, 10 figures, PDFLaTeX, accepted for publication in A&

Full-star Type Ia supernova explosion 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.Comment: 11 pages, 10 figures, resolution of some figures reduced to meet astro-ph file size restriction, submitted to A&

Nucleosynthesis in thermonuclear supernovae with tracers: convergence and variable mass particles

Nucleosynthetic yield predictions for multi-dimensional simulations of thermonuclear supernovae generally rely on the tracer particle method to obtain isotopic information of the ejected material for a given supernova simulation. We investigate how many tracer particles are required to determine converged integrated total nucleosynthetic yields. For this purpose, we conduct a resolution study in the number of tracer particles for different hydrodynamical explosion models at fixed spatial resolution. We perform hydrodynamic simulations on a co-expanding Eulerian grid in two dimensions assuming rotational symmetry for both pure deflagration and delayed detonation Type Ia supernova explosions. Within a given explosion model, we vary the number of tracer particles to determine the minimum needed for the method to give a robust prediction of the integrated yields of the most abundant nuclides. For the first time, we relax the usual assumption of constant tracer particle mass and introduce a radially vary- ing distribution of tracer particle masses. We find that the nucleosynthetic 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 - more or less independent of the explosion model. We conclude that the number of tracer particles that were used in extant published works appear to have been sufficient as far as integrated yields are concerned for the most copiously produced nuclides. Additionally we find that a suitably chosen tracer mass distribution can improve convergence for nuclei produced in the outer layer of the supernova where the constant tracer mass prescription suffers from poor spatial resolution.Comment: 9 pages, 5 figures, accepted for publication in MNRA

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

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

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

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

Towards an understanding of Type Ia supernovae from a synthesis of theory and observations

Motivated by the fact that calibrated light curves of Type Ia supernovae (SNe Ia) have become a major tool to determine the expansion history of the Universe, considerable attention has been given to, both, observations and models of these events over the past 15 years. Here, we summarize new observational constraints, address recent progress in modeling Type Ia supernovae by means of three-dimensional hydrodynamic simulations, and discuss several of the still open questions. It will be be shown that the new models have considerable predictive power which allows us to study observable properties such as light curves and spectra without adjustable non-physical parameters. This is a necessary requisite to improve our understanding of the explosion mechanism and to settle the question of the applicability of SNe Ia as distance indicators for cosmology. We explore the capabilities of the models by comparing them with observations and we show how such models can be applied to study the origin of the diversity of SNe Ia.Comment: 26 pages, 13 figures, Frontiers of Physics, in prin

Towards a Unified Quark-Hadron Matter Equation of State for Applications in Astrophysics and Heavy-Ion Collisions

We outline an approach to a unified equation of state for quark-hadron matter on the basis of a $\Phi-$derivable approach to the generalized Beth-Uhlenbeck equation of state for a cluster decomposition of thermodynamic quantities like the density. To this end we summarize the cluster virial expansion for nuclear matter and demonstrate the equivalence of the Green's function approach and the $\Phi-$derivable formulation. For an example, the formation and dissociation of deuterons in nuclear matter is discussed. We formulate the cluster $\Phi-$derivable approach to quark-hadron matter which allows to take into account the specifics of chiral symmetry restoration and deconfinement in triggering the Mott-dissociation of hadrons. This approach unifies the description of a strongly coupled quark-gluon plasma with that of a medium-modified hadron resonance gas description which are contained as limiting cases. The developed formalism shall replace the common two-phase approach to the description of the deconfinement and chiral phase transition that requires a phase transition construction between separately developed equations of state for hadronic and quark matter phases. Applications to the phenomenology of heavy-ion collisions and astrophysics are outlined.Comment: 35 pages, 3 figures, Special Issue "Compact Stars in the QCD Phase Diagram

Delayed detonations in full-star models of Type Ia supernova explosions

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.Comment: 4 pages, 1 figur

Three-dimensional modeling of Type Ia supernovae - The power of late time spectra

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.Comment: 20 pages, 9 figures, resolution of Fig 1 is reduced to meet astro-ph file size restriction, submitted to A&