George C. Jordan

Three-Dimensional Simulations of the Deflagration Phase of the Gravitationally Confined Detonation Model of Type Ia Supernovae

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 56Ni 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.

FAILED-DETONATION SUPERNOVAE: SUBLUMINOUS LOW-VELOCITY Ia SUPERNOVAE AND THEIR KICKED REMNANT WHITE DWARFS WITH IRON-RICH CORES

Type Ia supernovae (SNe Ia) originate from the thermonuclear explosions of carbon-oxygen (C-O) white dwarfs (WDs). The single-degenerate scenario is a well-explored model of SNe Ia where unstable thermonuclear burning initiates in an accreting, Chandrasekhar-mass WD and forms an advancing flame. By several proposed physical processes, the rising, burning material triggers a detonation, which subsequently consumes and unbinds the WD. However, if a detonation is not triggered and the deflagration is too weak to unbind the star, a completely different scenario unfolds. We explore the failure of the gravitationally confined detonation mechanism of SNe Ia, and demonstrate through two-dimensional and three-dimensional simulations the properties of failed-detonation SNe. We show that failed-detonation SNe expel a few 0.1 M ☉ of burned and partially burned material and that a fraction of the material falls back onto the WD, polluting the remnant WD with intermediate-mass and iron-group elements that likely segregate to the core forming a WD whose core is iron rich. The remaining material is asymmetrically ejected at velocities comparable to the escape velocity from the WD, and in response, the WD is kicked to velocities of a few hundred km s–1. These kicks may unbind the binary and eject a runaway/hypervelocity WD. Although the energy and ejected mass of the failed-detonation SN are a fraction of typical thermonuclear SNe, they are likely to appear as subluminous low-velocity SNe Ia. Such failed detonations might therefore explain or are related to the observed branch of peculiar SNe Ia, such as the family of low-velocity subluminous SNe (SN 2002cx/SN 2008ha-like SNe).

Study of the Detonation Phase in the Gravitationally Confined Detonation Model of Type Ia Supernovae

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 56Ni 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 ? 1051 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 56Ni mass) SNe is discussed.

THE POST-MERGER MAGNETIZED EVOLUTION OF WHITE DWARF BINARIES: THE DOUBLE-DEGENERATE CHANNEL OF SUB-CHANDRASEKHAR TYPE Ia SUPERNOVAE AND THE FORMATION OF MAGNETIZED WHITE DWARFS

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 × 108 G. We discuss the impact of our findings on the origins, circumstellar media, and observed properties of SNe Ia and magnetized white dwarfs.

RADIATION TRANSPORT FOR EXPLOSIVE OUTFLOWS: A MULTIGROUP HYBRID MONTE CARLO METHOD

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 MC 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 gray 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.

Pragmatic optimizations for better scientific utilization of large supercomputers

Advances in modeling and algorithms, combined with growth in computing resources, have enabled simulations of multiphysics–multiscale phenomena that can greatly enhance our scientific understanding. However, on currently available high-performance computing (HPC) resources, maximizing the scientific outcome of simulations requires many trade-offs. In this paper we describe our experiences in running simulations of the explosion phase of Type Ia supernovae on the largest available platforms. The simulations use FLASH, a modular, adaptive mesh, parallel simulation code with a wide user base. The simulations use multiple physics components: hydrodynamics, gravity, a sub-grid flame model, a three-stage burning model, and a degenerate equation of state. They also use Lagrangian tracer particles, which are then post-processed to determine the nucleosynthetic yields. We describe the simulation planning process, and the algorithmic optimizations and trade-offs that were found to be necessary. Several of the optimizations and trade-offs were made during the course of the simulations as our understanding of the challenges evolved, or when simulations went into previously unexplored physical regimes. We also briefly outline the anticipated challenges of, and our preparations for, the next-generation computing platforms.