The merger of two white dwarfs (WDs) has for many years not been considered as the favoured model for the progenitor system of type Ia supernovae (SNe Ia). But recent years have seen a change of opinion as a number of studies, both observational and theoretical, have concluded that they should contribute significantly to the observed type Ia supernova rate. In this paper, we study the ignition and propagation of detonation through post-merger remnants and we follow the resulting nucleosynthesis up to the point where a homologous expansion is reached. In our study we cover the entire range of WD masses and compositions. For the emergence of a detonation we study three different setups. The first two are guided by the merger remnants from our earlier simulations (Dan et al. 2014}, while for the third one the ignitions were set by placing hotspots with properties determined by spatially resolved calculations taken from the literature. There are some caveats to our approach which we investigate. We carefully compare the nucleosynthetic yields of successful explosions with SN Ia observations. Only three of our models are consistent with all the imposed constraints and potentially lead to a standard type Ia event. The first one, a $0.45\,M_\odot\ {\rm helium\ (He)} + 0.9 \,M_\odot$ carbon-oxygen (CO) WD system produces a sub-luminous, SN 1991bg-like event while the other two, a $0.45\,M_\odot\ {\rm He} + 1.1\,M_\odot$ oxygen-neon (ONe) WD system and a $1.05+1.05\,M_\odot$ system with two CO WDs, are good candidates for common SNe Ia.

We present a large parameter study where we investigate the structure of white dwarf (WD) merger remnants after the dynamical phase. A wide range of WD masses and compositions are explored and we also probe the effect of different initial conditions. We investigated the degree of mixing between the WDs, the conditions for detonations as well as the amount of gas ejected. We find that systems with lower mass ratios have more total angular momentum and as a result more mass is flung out in a tidal tail. Nuclear burning can affect the amount of mass ejected. Many WD binaries that contain a helium-rich WD achieve the conditions to trigger a detonation. In contrast, for carbon-oxygen transferring systems only the most massive mergers with a total mass above ~2.1 solar masses detonate. Even systems with lower mass may detonate long after the merger if the remnant remains above the Chandrasekhar mass and carbon is ignited at the centre. Finally, our findings are discussed in the context of several possible observed astrophysical events and stellar systems, such as hot subdwarfs, R Coronae Borealis stars, single massive white dwarfs, supernovae of type Ia and other transient events. A large database containing 225 white dwarf merger remnants is made available via a dedicated web page.

The merger of two white dwarfs (WDs) creates a differentially rotating remnant which is unstable to magnetohydrodynamic instabilities. These instabilities can lead to viscous evolution on a time-scale short compared to the thermal evolution of the remnant. We present multi-dimensional hydrodynamic simulations of the evolution of WD merger remnants under the action of an $\alpha$-viscosity. We initialize our calculations using the output of eight WD merger simulations from Dan et al. (2011), which span a range of mass ratios and total masses. We generically find that the merger remnants evolve towards spherical states on time-scales of hours, even though a significant fraction of the mass is initially rotationally supported. The viscous evolution unbinds only a very small amount of mass $(< 10^{-5} M_\odot)$. Viscous heating causes some of the systems we study with He WD secondaries to reach conditions of nearly dynamical burning. It is thus possible that the post-merger viscous phase triggers detonation of the He envelope in some WD mergers, potentially producing a Type Ia supernova via a double detonation scenario. Our calculations provide the proper initial conditions for studying the long-term thermal evolution of WD merger remnants. This is important for understanding WD mergers as progenitors of Type Ia supernovae, neutron stars, R Coronae Borealis stars and other phenomena.

Double white dwarf binaries in the Galaxy dominate the gravitational wave sky and would be detectable for an instrument such as LISA. Most studies have calculated the expected gravitational wave signal under the assumption that the binary white dwarf system can be represented by two point masses in orbit. We discuss the accuracy of this approximation for real astrophysical systems. For non-relativistic binaries in circular orbit the gravitational wave signal can easily be calculated. We show that for these systems the point mass approximation is completely justified when the individual stars are axisymmetric irrespective of their size. We find that the signal obtained from Smoothed-Particle Hydrodynamics simulations of tidally deformed, Roche-lobe filling white dwarfs, including one case when an accretion disc is present, is consistent with the point mass approximation. The difference is typically at the level of one per cent or less in realistic cases, yielding small errors in the inferred parameters of the binaries.

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.

We present the results of a systematic numerical study of the onset of mass transfer in double degenerate binary systems and its impact on the subsequent evolution. All investigated systems belong to the regime of direct impact, unstable mass transfer. In all of the investigated cases, even those considered unstable by conventional stability analysis, we find a long-lived mass transfer phase continuing for as many as several dozen orbital periods. This settles a recent debate sparked by a discrepancy between earlier SPH calculations that showed disruptions after a few orbital periods and newer grid-based studies in which mass transfer continued for tens of orbits. As we show that these binaries can survive at small separation for hundreds of orbital periods, their associated gravitational wave signal should be included when calculating the gravitational wave foreground (although expected to below LISA's sensitivity at these high frequencies). We also show that the inclusion of the entropy increase associated with shock-heating of the accreted material reduces the number of orbits a binary survives given the same initial conditions, although the effect is not as pronounced when using the appropriate initial conditions. The use of accurate initial conditions and a correct treatment of shock heating allows for a reliable time evolution of the temperature, density, and angular momentum, which are important when considering thermonuclear events that may occur during the mass transfer phase and/or after merger. Our treatment allows us to accurately identify when surface detonations may occur in the lead-up to the merger, as well as the properties of final merger products.

We present three-dimensional simulations on a new mechanism for the detonation of a sub-Chandrasekhar CO white dwarf in a dynamically unstable system where the secondary is either a pure He white dwarf or a He/CO hybrid. For dynamically unstable systems where the accretion stream directly impacts the surface of the primary, the final tens of orbits can have mass accretion rates that range from $10^{-5}$ to $10^{-3} M_{\odot}$ s$^{-1}$, leading to the rapid accumulation of helium on the surface of the primary. After $\sim 10^{-2}$ $M_{\odot}$ of helium has been accreted, the ram pressure of the hot helium torus can deflect the accretion stream such that the stream no longer directly impacts the surface. The velocity difference between the stream and the torus produces shearing which seeds large-scale Kelvin-Helmholtz instabilities along the interface between the two regions. These instabilities eventually grow into dense knots of material that periodically strike the surface of the primary, adiabatically compressing the underlying helium torus. If the temperature of the compressed material is raised above a critical temperature, the timescale for triple-$\alpha$ reactions becomes comparable to the dynamical timescale, leading to the detonation of the primary's helium envelope. This detonation drives shockwaves into the primary which tend to concentrate at one or more focal points within the primary's CO core. If a relatively small amount of mass is raised above a critical temperature and density at these focal points, the CO core may itself be detonated.

We present a numerical study of the mass transfer dynamics prior to the gravitational wave-driven merger of a double white dwarf system. Recently, there has been some discussion about the dynamics of these last stages, different methods seemed to provide qualitatively different results. While earlier SPH simulations indicated a very quick disruption of the binary on roughly the orbital time scale, more recent grid-based calculations find long-lived mass transfer for many orbital periods. Here we demonstrate how sensitive the dynamics of this last stage is to the exact initial conditions. We show that, after a careful preparation of the initial conditions, the reportedly short-lived systems undergo mass transfer for many dozens of orbits. The reported numbers of orbits are resolution-biased and therefore represent only lower limits to what is realized in nature. Nevertheless, the study shows convincingly the convergence of different methods to very similar results.

The existence of supermassive black holes lurking in the centers of galaxies and of stellar binary systems containing a black hole with a few solar masses has been established beyond reasonable doubt. The idea that black holes of intermediate masses ($\sim 1000$ \msun) may exist in globular star clusters has gained credence over recent years but no conclusive evidence has been established yet. An attractive feature of this hypothesis is the potential to not only disrupt solar-type stars but also compact white dwarf stars. In close encounters the white dwarfs can be sufficiently compressed to thermonuclearly explode. The detection of an underluminous thermonuclear explosion accompanied by a soft, transient X-ray signal would be compelling evidence for the presence of intermediate mass black holes in stellar clusters. In this paper we focus on the numerical techniques used to simulate the entire disruption process from the initial parabolic orbit, over the nuclear energy release during tidal compression, the subsequent ejection of freshly synthesized material and the formation process of an accretion disk around the black hole.