[Paper Review] NuGrid: Toward High Precision Double-Degenerate Merger Simulations with SPH in 3D

Steven Diehl, Chris L. Fryer+12 more|ArXiv.org|Nov 28, 2008

Gamma-ray bursts and supernovae3 citations

TL;DR

This paper presents high-precision 3D SPH simulations of double-degenerate white dwarf mergers, emphasizing rigorous verification and validation through improved initial condition setup using weighted Voronoi tesselations and the inclusion of shocks via an ideal gas equation of state. Key results show that shocks dramatically alter merger dynamics, forming a hot, expanding halo and ejecting mass with angular momentum, which may impact the viability of such mergers as Type Ia supernova progenitors.

ABSTRACT

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

Motivation & Objective

• To conduct high-precision, verified, and validated 3D SPH simulations of double-degenerate white dwarf mergers to assess their viability as Type Ia supernova progenitors.
• To investigate the impact of initial condition setup—particularly particle distribution and co-rotation—on simulation accuracy and dynamics.
• To evaluate the role of shocks in merger dynamics by comparing simulations with and without shock physics via different equations of state.
• To enable post-processing of SPH particle trajectories with a full nuclear network to model dynamic nuclear reactions and compare with R Coronae Borealis star abundances.
• To establish a foundation for future simulations using realistic equations of state and stellar evolution models within the NuGrid collaboration framework.

Proposed method

• Utilizes a modified version of the SNSPH code with high-resolution smoothed particle hydrodynamics (SPH) in 3D to simulate binary white dwarf mergers.
• Employs the self-consistent field method to generate equilibrium initial conditions for co-rotating, unequal-mass white dwarf binaries with minimal numerical artifacts.
• Introduces a novel weighted Voronoi tesselation (WVT)-based particle setup to achieve high resolution in outer layers while minimizing particle noise and density fluctuations.
• Compares simulations using a polytropic equation of state (no shocks) versus an ideal gas equation of state (with shocks) to isolate the dynamical impact of shock heating.
• Applies post-processing with the tppnp tool to track SPH particles through a complete nuclear network, enabling synthesis of isotopes and comparison with observed R Coronae Borealis star abundances.
• Conducts verification via code comparison with grid-based simulations and numerical convergence studies, and validation against observational constraints from R Coronae Borealis stars.

Experimental results

Research questions

• RQ1 How do different initial condition setups—particularly particle distribution and rotational profiles—affect the outcome of double-degenerate merger simulations?
• RQ2 What is the dynamical impact of including shocks in the equation of state on merger evolution and mass ejection?
• RQ3 How do the resulting merger remnants, including their structure and angular momentum loss, compare with observational constraints from R Coronae Borealis stars?
• RQ4 To what extent do shock-heated, high-entropy halos formed during mergers alter the merger timescale and energy budget?
• RQ5 Can high-precision SPH simulations with realistic equations of state and nuclear network post-processing confirm the viability of double-degenerate mergers as progenitors of Type Ia supernovae?

Key findings

• The inclusion of shocks via an ideal gas equation of state leads to the formation of a hot, expanding halo of shocked material that engulfs both white dwarfs, fundamentally altering merger dynamics compared to shock-free polytropic simulations.
• Shocked gas is efficiently ejected from the backside of both stars, carrying away significant angular momentum and potentially accelerating orbital decay.
• The merger remnant settles into a differentially rotating, nearly spherically symmetric configuration with a fast-rotating core and a hot envelope extending to approximately 100 solar radii at the end of the simulation.
• The simulation shows that the accretion stream is resolved with over 1,000 SPH particles on average, indicating high resolution and reduced numerical noise due to the WVT-based particle setup.
• Non-corotating systems exhibit artificial orbital shrinkage due to angular momentum transfer to spin, highlighting the importance of proper initial rotational setup.
• The results demonstrate that the equation of state and initial condition quality are critical factors in determining merger outcomes, with shocks playing a central dynamical role.