No abstract

We report the results of a series of three-dimensional (3D) simulations of the deflagration phase of the gravitationally confined detonation mechanism for Type Ia supernovae. In this mechanism, ignition occurs at one or several off-center points, resulting in a burning bubble of hot ash that rises rapidly, breaks through the surface of the star, and collides at a point opposite the breakout on the stellar surface. We find that detonation conditions are robustly reached in our 3D simulations for a range of initial conditions and resolutions. Detonation conditions are achieved as the result of an inwardly directed jet that is produced by the compression of unburnt surface material when the surface flow collides with itself. A high-velocity outwardly directed jet is also produced. The initial conditions explored in this paper lead to conditions at detonation that can be expected to produce large amounts of 56 Ni and small amounts of intermediate-mass elements. These particular simulations are therefore relevant only to high-luminosity Type Ia supernovae. Recent observations of Type Ia supernovae imply a compositional structure that is qualitatively consistent with that expected from these simulations.

We study the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia) through the detonation phase and into homologous expansion. In the GCD model, a detonation is triggered by the surface flow due to single-point, off-center flame ignition in carbon-oxygen white dwarfs (WDs). The simulations are unique in terms of the degree to which nonidealized physics is used to treat the reactive flow, including weak reaction rates and a time-dependent treatment of material in nuclear statistical equilibrium (NSE). Careful attention is paid to accurately calculating the final composition of material which is burned to NSE and frozen out in the rapid expansion following the passage of a detonation wave over the high-density core of the WD; and an efficient method for nucleosynthesis postprocessing is developed which obviates the need for costly network calculations along tracer particle thermodynamic trajectories. Observational diagnostics are presented for the explosion models, including abundance stratifications and integrated yields. We find that for all of the ignition conditions studied here a self-regulating process comprised of neutronization and stellar expansion results in final 56 Ni masses of ∼1.1 M . But, more energetic models result in larger total NSE and stable Fe-peak yields. The total yield of intermediate mass elements is ∼ 0.1 M and the explosion energies are all around 1.5 × 10 51 erg. The explosion models are briefly compared to the inferred properties of recent SN Ia observations. The potential for surface detonation models to produce lower-luminosity (lower 56 Ni mass) SNe is discussed.

We explore Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. DDMC accelerates IMC in diffusive domains. Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally grey DDMC method. We rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. This formulation includes an analysis that yields an additional factor in the standard IMC-to-DDMC spatial interface condition. To our knowledge the new boundary condition is distinct from others presented in prior DDMC literature. The method is suitable for a variety of opacity distributions and may be applied to semi-relativistic radiation transport in simple fluids and geometries. Additionally, we test the code, called SuperNu, using an analytic solution having static material, as well as with a manufactured solution for moving material with structured opacities. Finally, we demonstrate with a simple source and 10 group logarithmic wavelength grid that IMC-DDMC performs better than pure IMC in terms of accuracy and speed when there are large disparities between the magnitudes of opacities in adjacent groups. We also present and test our implementation of the new boundary condition.

Type Ia supernovae (SNe Ia) play a crucial role as standardizable cosmological candles, though the nature of their progenitors is a subject of active investigation. Recent observational and theoretical work has pointed to merging white dwarf binaries, referred to as the double-degenerate channel, as the possible progenitor systems for some SNe Ia. Additionally, recent theoretical work suggests that mergers which fail to detonate may produce magnetized, rapidly-rotating white dwarfs. In this paper, we present the first multidimensional simulations of the post-merger evolution of white dwarf binaries to include the effect of the magnetic field. In these systems, the two white dwarfs complete a final merger on a dynamical timescale, and are tidally disrupted, producing a rapidly-rotating white dwarf merger surrounded by a hot corona and a thick, differentially-rotating disk. The disk is strongly susceptible to the magnetorotational instability (MRI), and we demonstrate that this leads to the rapid growth of an initially dynamically weak magnetic field in the disk, the spin-down of the white dwarf merger, and to the subsequent central ignition of the white dwarf merger. Additionally, these magnetized models exhibit new features not present in prior hydrodynamic studies of white dwarf mergers, including the development of MRI turbulence in the hot disk, magnetized outflows carrying a significant fraction of the disk mass, and the magnetization of the white dwarf merger to field strengths ∼ 2 × 10 8 G. We discuss the impact of our findings on the origins, circumstellar media, and observed properties of SNe Ia and magnetized white dwarfs.

No abstract

We present a systematic study of the pure deflagration model of Type Ia supernovae using three-dimensional, high-resolution, full-star hydrodynamical simulations, nucleosynthetic yields calculated using Lagrangian tracer particles, and light curves calculated using radiation transport. We evaluate the simulations by comparing their predicted light curves with many observed SNe Ia using the SALT2 data-driven model and find that the simulations may correspond to under-luminous SNe Iax. We explore the effects of the initial conditions on our results by varying the number of randomly selected ignition points from 63 to 3500, and the radius of the centered sphere they are confined in from 128 to 384 km. We find that the rate of nuclear burning depends on the number of ignition points at early times, the density of ignition points at intermediate times, and the radius of the confining sphere at late times. The results depend primarily on the number of ignition points, but we do not expect this to be the case in general. The simulations with few ignition points release more nuclear energy E nuc , have larger kinetic energies E K , and produce more 56 Ni than those with many ignition points, and differ in the distribution of 56 Ni, Si, and C/O in the ejecta. For these reasons, the simulations with few ignition points exhibit higher peak B-band absolute magnitudes M B and light curves that rise and decline more quickly; their M B and light curves resemble those of under-luminous SNe Iax, while those for simulations with many ignition points are not.

We demonstrate that nucleosynthesis in rapid, high-entropy expansions of proton-rich matter from high temperature and density can result in a wider variety of abundance patterns than heretofore appreciated. In particular, such expansions can produce iron-group nuclides, p-process nuclei, or even heavy, neutron-rich isotopes. Such diversity arises because the nucleosynthesis enters a little explored regime in which the free nucleons are not in equilibrium with the abundant 4 He. This allows nuclei significantly heavier than iron to form in the presence of abundant free nucleons early in the expansion. As the temperature drops, nucleons increasingly assemble into 4 He and heavier nuclei. If the assembly is efficient, the resulting depletion of free neutrons allows disintegration flows to drive nuclei back down to iron and nickel. If this assembly is inefficient, then the large abundance of free nucleons prevents the disintegration flows and leaves a distribution of heavy nuclei after reaction freezeout. For cases in between, an intermediate abundance distribution, enriched in p-process isotopes, is frozen out. These last expansions may contribute to the solar system's supply of the p-process nuclides if mildly proton-rich, high-entropy matter is ejected from proto-neutron stars winds or other astrophysical sites. Also significant is the fact that, because the nucleosynthesis is primary, the signature of this nucleosynthesis may be evident in metal-poor stars.