148 research outputs found
Type Ia Supernova: Burning and Detonation in the Distributed Regime
A simple, semi-analytic representation is developed for nuclear burning in
Type Ia supernovae in the special case where turbulent eddies completely
disrupt the flame. The speed and width of the ``distributed'' flame front are
derived. For the conditions considered, the burning front can be considered as
a turbulent flame brush composed of corrugated sheets of well-mixed flames.
These flames are assumed to have a quasi-steady-state structure similar to the
laminar flame structure, but controlled by turbulent diffusion. Detonations
cannot appear in the system as long as distributed flames are still
quasi-steady-state, but this condition is violated when the distributed flame
width becomes comparable to the size of largest turbulent eddies. When this
happens, a transition to detonation may occur. For current best estimates of
the turbulent energy, the most likely density for the transition to detonation
is in the range 0.5 - 1.5 x 10^7 g cm^{-3}.Comment: 12 pages, 4 figure
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 12 high-resolution three-dimensional delayed detonation SN Ia explosion simulations that employ a new criterion to trigger the deflagration to detonation transition (DDT). The simulations fall into three ignition categories: relatively bright SNe with five ignition kernels and a weak deflagration phase (three different central densities); relatively dim SNe with 1600 ignition kernels and a strong deflagration phase (three different central densities) and intermediate SNe with 200 ignition kernels (six different central densities). 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 all three categories. 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 offsets the percental reduction of 56Ni due to neutronization
A subgrid-scale model for deflagration-to-detonation transitions in Type la supernova explosion simulations Numerical implementation
Delayed detonations of Chandrasekhar-mass white dwarfs are a promising model for normal Type Ia supernova explosions. In these white dwarfs, 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 recently determined. The region of detonation initiation cannot be spatially resolved in multidimensional full-star simulations of the explosion
On the small-scale stability of thermonuclear flames in Type Ia supernovae
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.Comment: 9 pages, 7 figures, Uses AASTEX, emulateapj5.sty, onecolfloat.sty.
Replaced with accepted version (ApJ), Figures 1 and 3 are ne
Constraints on the high-density nuclear equation of state from the phenomenology of compact stars and heavy-ion collisions
A new scheme for testing nuclear matter equations of state (EsoS) at high
densities using constraints from neutron star phenomenology and a flow data
analysis of heavy-ion collisions is suggested. An acceptable EoS shall not
allow the direct Urca process to occur in neutron stars with masses below
, and also shall not contradict flow and kaon production data of
heavy-ion collisions. Compact star constraints include the mass measurements of
2.1 +/- 0.2 M_sun (1 sigma level) for PSR J0751+1807, of 2.0 +/- 0.1 M_sun from
the innermost stable circular orbit for 4U 1636-536, the baryon mass -
gravitational mass relationships from Pulsar B in J0737-3039 and the
mass-radius relationships from quasiperiodic brightness oscillations in 4U
0614+09 and from the thermal emission of RX J1856-3754. This scheme is applied
to a set of relativistic EsoS constrained otherwise from nuclear matter
saturation properties with the result that no EoS can satisfy all constraints
simultaneously, but those with density-dependent masses and coupling constants
appear most promising.Comment: 15 pages, 8 figures, 5 table
Three-dimensional delayed-detonation models with nucleosynthesis for type ia supernovae
We present results for a suite of 14 three-dimensional, high-resolution hydrodynamical simulations of delayed-detonation models of Type Ia supernova (SN Ia) explosions. This model suite comprises the first set of three-dimensional SN Ia simulations with detailed isotopic yield information. As such, it may serve as a data base for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ a physically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration-to-detonation transition probabilities. To obtain different strengths of the deflagration phase and thereby different degrees of pre-expansion, we have chosen a sequence of initial models with 1, 3, 5, 10, 20, 40, 100, 150, 200, 300 and 1600 (two different realizations) ignition kernels in a hydrostatic white dwarf with a central density of 2.9 × 109 g cm−3, as well as one high central density (5.5 × 109 g cm−3) and one low central density (1.0 × 109 g cm−3) rendition of the 100 ignition kernel configuration. For each simulation, we determined detailed nucleosynthetic yields by post-processing 106 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 M⊙. As a general trend, the models predict that the stable neutron-rich iron-group isotopes are not found at the lowest velocities, but rather at intermediate velocities (∼3000–10 000 km s−1) in a shell surrounding a 56Ni-rich core. The models further predict relatively low-velocity oxygen and carbon, with typical minimum velocities around 4000 and 10 000 km s−1, respectively
Equation of State of Nuclear Matter at high baryon density
A central issue in the theory of astrophysical compact objects and heavy ion
reactions at intermediate and relativistic energies is the Nuclear Equation of
State (EoS). On one hand, the large and expanding set of experimental and
observational data is expected to constrain the behaviour of the nuclear EoS,
especially at density above saturation, where it is directly linked to
fundamental processes which can occur in dense matter. On the other hand,
theoretical predictions for the EoS at high density can be challenged by the
phenomenological findings. In this topical review paper we present the
many-body theory of nuclear matter as developed along different years and with
different methods. Only nucleonic degrees of freedom are considered. We compare
the different methods at formal level, as well as the final EoS calculated
within each one of the considered many-body schemes. The outcome of this
analysis should help in restricting the uncertainty of the theoretical
predictions for the nuclear EoS.Comment: 51 pages, to appear in J. Phys. G as Topical Revie
Predicting polarization signatures for double-detonation and delayed-detonation models of Type Ia supernovae
Calculations of synthetic spectropolarimetry are one means to test multidimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M⊙ 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 2D double-detonation model and a 3D 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 ˜0.1-0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si II λ6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O I λ7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt
500 Days of SN 2013dy: spectra and photometry from the ultraviolet to the infrared
SN 2013dy is a Type Ia supernova for which we have compiled an extraordinary
dataset spanning from 0.1 to ~ 500 days after explosion. We present 10 epochs
of ultraviolet (UV) through near-infrared (NIR) spectra with HST/STIS, 47
epochs of optical spectra (15 of them having high resolution), and more than
500 photometric observations in the BVrRiIZYJH bands. SN 2013dy has a broad and
slowly declining light curve (delta m(B) = 0.92 mag), shallow Si II 6355
absorption, and a low velocity gradient. We detect strong C II in our earliest
spectra, probing unburned progenitor material in the outermost layers of the SN
ejecta, but this feature fades within a few days. The UV continuum of SN
2013dy, which is strongly affected by the metal abundance of the progenitor
star, suggests that SN 2013dy had a relatively high-metallicity progenitor.
Examining one of the largest single set of high-resolution spectra for a SN Ia,
we find no evidence of variable absorption from circumstellar material.
Combining our UV spectra, NIR photometry, and high-cadence optical photometry,
we construct a bolometric light curve, showing that SN 2013dy had a maximum
luminosity of 10.0^{+4.8}_{-3.8} * 10^{42} erg/s. We compare the synthetic
light curves and spectra of several models to SN 2013dy, finding that SN 2013dy
is in good agreement with a solar-metallicity W7 model.Comment: 22 pages, 18 figures, replaced with version accecpted for publication
in MNRA
- …