Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical Stability of Detonations in White Dwarf Simulations
64
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Some simulations of Type Ia supernovae feature self-consistent thermonuclear detonations. However, these detonations are not meaningful if the simulations are not resolved, so it is important to establish the requirements for achieving a numerically converged detonation. In this study we examine a test detonation problem inspired by collisions of white dwarfs. This test problem demonstrates that achieving a converged thermonuclear ignition requires spatial resolution much finer than 1 km in the burning region. Current computational resource constraints place this stringent resolution requirement out of reach for multi-dimensional supernova simulations. Consequently, contemporary simulations that self-consistently demonstrate detonations are possibly not converged and should be treated with caution.
rate research
Read More
Type Ia supernovae (SNe Ia) span a range of luminosities and timescales, from rapidly evolving subluminous to slowly evolving overluminous subtypes. Previous theoretical work has, for the most part, been unable to match the entire breadth of observed SNe Ia with one progenitor scenario. Here, for the first time, we apply non-local thermodynamic equilibrium radiative transfer calculations to a range of accurate explosion models of sub-Chandrasekhar-mass white dwarf detonations. The resulting photometry and spectra are in excellent agreement with the range of observed non-peculiar SNe Ia through 15 d after the time of B-band maximum, yielding one of the first examples of a quantitative match to the entire Phillips (1993) relation. The intermediate-mass element velocities inferred from theoretical spectra at maximum light for the more massive white dwarf explosions are higher than those of bright observed SNe Ia, but these and other discrepancies likely stem from the one-dimensional nature of our explosion models and will be improved upon by future non-local thermodynamic equilibrium radiation transport calculations of multi-dimensional sub-Chandrasekhar-mass white dwarf detonations.
We study the violent phase of the merger of massive binary white dwarf systems. Our aim is to characterize the conditions for explosive burning to occur, and identify a possible explosion mechanism of Type Ia supernovae. The primary components of our model systems are carbon-oxygen (C/O) white dwarfs, while the secondaries are made either of C/O or of pure helium. We account for tidal effects in the initial conditions in a self-consistent way, and consider initially well-separated systems with slow inspiral rates. We study the merger evolution using an adaptive mesh refinement, reactive, Eulerian code in three dimensions, assuming symmetry across the orbital plane. We use a co-rotating reference frame to minimize the effects of numerical diffusion, and solve for self-gravity using a multi-grid approach. We find a novel detonation mechanism in C/O mergers with massive primaries. Here the detonation occurs in the primarys core and relies on the combined action of tidal heating, accretion heating, and self-heating due to nuclear burning. The exploding structure is compositionally stratified, with a reverse shock formed at the surface of the dense ejecta. The existence of such a shock has not been reported elsewhere. The explosion energy ($1.6times 10^{51}$ erg) and $^{56}$Ni mass (0.86 M$_odot$) are consistent with a SN Ia at the bright end of the luminosity distribution, with an approximated decline rate of $Delta m_{15}(B)approx 0.99$. Our study does not support double-detonation scenarios in the case of a system with a 0.6 M$_odot$ helium secondary and a 0.9 M$_odot$ primary. Although the accreted helium detonates, it fails to ignite carbon at the base of the boundary layer or in the primarys core.
Despite their unique astrophysical relevance, the outcome of white dwarf binary mergers has so far only been studied for a very restricted number of systems. Here we present the results of a survey with more than two hundred simulations systematically scanning the white dwarf binary parameter space. We consider white dwarf masses ranging from 0.2 to 1.2 $M_odot$ and account for their different chemical compositions. We find excellent agreement with the orbital evolution predicted by mass transfer stability analysis. Much of our effort in this paper is dedicated to determining which binary systems are prone to a thermonuclear explosion just prior to merger or at surface contact. We find that a large fraction of He-accreting binary systems explode: all dynamically unstable systems with accretor masses below 1.1 $M_odot$ and donor masses above $sim$ 0.4 $M_odot$ are found to trigger a helium detonation at surface contact. A substantial fraction of these systems could explode at earlier times via detonations induced by instabilities in the accretion stream, as we have demonstrated in our previous work. We do not find definitive evidence for an explosion prior to merger or at surface contact in any of the studied double carbon-oxygen systems. Although we cannot exclude their occurrence if some helium is present, the available parameter space for a successful detonation in a white dwarf binary of pure carbon-oxygen composition is small. We demonstrate that a wide variety of dynamically unstable systems are viable type Ia candidates. The next decade thus holds enormous promise for the study of these events, in particular with the advent of wide-field synoptic surveys allowing a detailed characterization of their explosive properties.
Merging white dwarfs are a possible progenitor of Type Ia supernovae (SNe Ia). While it is not entirely clear if and when an explosion is triggered in such systems, numerical models suggest that a detonation might be initiated before the stars have coalesced to form a single compact object. Here we study such peri-merger detonations by means of numerical simulations, modeling the disruption and nucleosynthesis of the stars until the ejecta reach the coasting phase. Synthetic light curves and spectra are generated for comparison with observations. Three models are considered with primary masses 0.96 Msun, 1.06 Msun, and 1.20 Msun. Of these, the 0.96 Msun dwarf merging with an 0.81 Msun companion, with a Ni56 yield of 0.58 Msun, is the most promising candidate for reproducing common SNe Ia. The more massive mergers produce unusually luminous SNe Ia with peak luminosities approaching those attributed to super-Chandrasekhar mass SNe Ia. While the synthetic light curves and spectra of some of the models resemble observed SNe Ia, the significant asymmetry of the ejecta leads to large orientation effects. The peak bolometric luminosity varies by more than a factor of 2 with the viewing angle, and the velocities of the spectral absorption features are lower when observed from angles where the light curve is brightest. The largest orientation effects are seen in the ultraviolet, where the flux varies by more than an order of magnitude. Despite the large variation with viewing angle, the set of three models roughly obeys a width-luminosity relation, with the brighter light curves declining more slowly in the B-band. Spectral features due to unburned carbon from the secondary star are also seen in some cases.
Merging carbon-oxygen (CO) white dwarfs are a promising progenitor system for Type Ia supernovae (SN Ia), but the underlying physics and timing of the detonation are still debated. If an explosion occurs after the secondary star is fully disrupted, the exploding primary will expand into a dense CO medium that may still have a disk-like structure. This interaction will decelerate and distort the ejecta. Here we carry out multi-dimensional simulations of ``tamped SN Ia models, using both particle and grid-based codes to study the merger and explosion dynamics, and a radiative transfer code to calculate synthetic spectra and light curves. We find that post-merger explosions exhibit an hourglass-shaped asymmetry, leading to strong variations in the light curves with viewing angle. The two most important factors affecting the outcome are the scale-height of the disk, which depends sensitively on the binary mass ratio, and the total ${}^{56}$Ni yield, which is governed by the central density of the remnant core. The synthetic broadband light curves rise and decline very slowly, and the spectra generally look peculiar, with weak features from intermediate mass elements but relatively strong carbon absorption. We also consider the effects of the viscous evolution of the remnant, and show that a longer time delay between merger and explosion probably leads to larger ${}^{56}$Ni yields and more symmetrical remnants. We discuss the relevance of this class of aspherical ``tamped SN Ia for explaining the class of ``super-Chandrasekhar SN Ia.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
