We present a novel relaxation method to build three-dimensional rotating structures of barotropic bodies using the SPH technique. The method is able to relax gaseous structures in rigid as well as differential rotation. The relaxation procedure strongly relies on the excellent conservation of angular momentum that characterizes the SPH technique. The method has been successfully applied to a variety of zero-temperature white dwarfs and polytropic self-gravitating structures. Our SPH results have been validated by comparing the main features (energies, central densities and the polar to equatorial radius ratio) to those obtained with independent, albeit grid-based methods, as for example, the self-consistent field method, showing that both methods agree within few percents.

Type Ia supernovae (SNe Ia) originate from the thermonuclear explosion of carbon-oxygen white dwarfs (CO WDs), giving rise to luminous optical transients. A relatively common variety of subluminous SNe Ia events, referred to as SNe Iax, are believed to arise from the failed detonation of a CO WD. In this paper, we explore failed detonation SNe Ia in the context of the double-degenerate channel of merging white dwarfs. In particular, we have carried out the first fully three-dimensional simulation of the merger of a ONe WD with a CO WD. While the hot, tidally-disrupted carbon-rich disk material originating from the CO WD secondary is readily susceptible to detonation in such a merger, the ONe WD primary core is not. This merger yields a failed detonation, resulting in the ejection of a small amount of mass, and leaving behind a kicked, super-Chandrasekhar ONe WD remnant enriched by the fallback of the products of nuclear burning. The resulting outburst is a rapidly-fading optical transient with a small amount of radioactive $^{56}$Ni powering the light curve. Consequently, the ONe-CO WD merger naturally produces a very faint and rapidly-fading transient, fainter even than the faintest Type Iax events observed to date, such as SN $2008$ha and SN 2010ae. More massive ONe primaries than considered here may produce brighter and longer-duration transients.

Increasing observational and theoretical evidence points to binary white dwarf mergers as the origin of some if not most normal Type Ia supernovae (SNe Ia). In this paper, we discuss the post-merger evolution of binary white dwarf (WD) mergers, and their relevance to the double-degenerate channel of SNe Ia. We present 3D simulations of carbon-oxygen (C/O) WD binary systems undergoing unstable mass transfer, varying both the total mass and the mass ratio. We demonstrate that these systems generally give rise to a one-armed gravitational spiral instability. The spiral density modes transport mass and angular momentum in the disk even in the absence of a magnetic field, and are most pronounced for secondary-to-primary mass ratios larger than $0.6$. We further analyze carbon burning in these systems to assess the possibility of detonation. Unlike the case of a $1.1 + 1.0 M_{\odot}$ C/O WD binary, we find that WD binary systems with lower mass and smaller mass ratios do not detonate as SNe Ia up to $\sim8-22$ outer dynamical times. Two additional models do however undergo net heating, and their secular increase in temperature could possibly result in a detonation on timescales longer than those considered here.

Recent advances in time-domain astronomy have uncovered a new class of optical transients with timescales shorter than typical supernovae and a wide range of peak luminosities. Several subtypes have been identified within this broad class, including Ca-rich transients, .Ia supernovae, and fast/bright transients. We examine the predic- tions from a state-of-the-art grid of three-dimensional simulations of dynamical white dwarf interactions in the context of these fast optical transients. We find that for colli- sions involving carbon-oxygen or oxygen-neon white dwarfs the peak luminosities and durations of the light curves in our models are in good agreement with the properties of fast/bright transients. When one of the colliding white dwarfs is made of helium the properties of the light curves are similar to those of Ca-rich gap transients. The model lightcurves from our white dwarf collisions are too slow to reproduce those of .Ia SNe, and too fast to match any normal or peculiar Type Ia supernova.

Isolated magnetic white dwarfs have field strengths ranging from kilogauss to gigagauss. However, the origin of the magnetic field has not been hitherto elucidated. Whether these fields are fossil, hence the remnants of original weak magnetic fields amplified during the course of the evolution of their progenitor stars, or are the result of binary interactions or, finally, they are produced by other internal physical mechanisms during the cooling of the white dwarf itself, remains a mystery. At sufficiently low temperatures white dwarfs crystallize. Upon solidification, phase separation of its main constituents, 12C and 16O, and of the impurities left by previous evolution occurs. This process leads to the formation of a Rayleigh-Taylor unstable liquid mantle on top of a solid core. This convective region, as it occurs in Solar System planets like the Earth and Jupiter, can produce a dynamo able to yield magnetic fields of strengths of up to 0.1 MG, thus providing a mechanism that could explain magnetism in single white dwarfs.

The core-degenerate (CD) scenario for type Ia supernovae (SN Ia) involves the merger of the hot core of an asymptotic giant branch (AGB) star and a white dwarf, and might contribute a non-negligible fraction of all thermonuclear supernovae. Despite its potential interest, very few studies, and based on only crude simplifications, have been devoted to investigate this possible scenario, compared with the large efforts invested to study some other scenarios. Here we perform the first three-dimensional simulations of the merger phase, and find that this process can lead to the formation of a massive white dwarf, as required by this scenario. We consider two situations, according to the mass of the circumbinary disk formed around the system during the final stages of the common envelope phase. If the disk is massive enough, the stars merge on a highly eccentric orbit. Otherwise, the merger occurs after the circumbinary disk has been ejected and gravitational wave radiation has brought the stars close to the Roche lobe radius on a nearly circular orbit. Not surprisingly, the overall characteristics of the merger remnants are similar to those found for the double-degenerate (DD) scenario, independently of the very different core temperature and of the orbits of the merging stars. They consist of a central massive white dwarf, surrounded by a hot, rapidly rotating corona and a thick debris region.

Thermonuclear, or Type Ia supernovae (SNe Ia), originate from the explosion of carbon--oxygen white dwarfs, and serve as standardizable cosmological candles. However, despite their importance, the nature of the progenitor systems that give rise to SNe Ia has not been hitherto elucidated. Observational evidence favors the double-degenerate channel in which merging white dwarf binaries lead to SNe Ia. Furthermore, significant discrepancies exist between observations and theory, and to date, there has been no self-consistent merger model that yields a SNe Ia. Here we show that a spiral mode instability in the accretion disk formed during a binary white dwarf merger leads to a detonation on a dynamical timescale. This mechanism sheds light on how white dwarf mergers may frequently yield SNe Ia.

We compute the possible observational signatures of white dwarf dynamical interactions in dense stellar environments. Specifically, we compute the emission of gravitational waves, and we compare it with the sensitivity curves of planned space-borne gravitational wave detectors. We also compute the light curves for those interactions in which a detonation occurs, and one of the stars is destroyed, as well as the corresponding neutrino luminosities. We find that for the three possible outcomes of these interactions - which are the formation of an eccentric binary system, a lateral collision in which several mass transfer episodes occur, and a direct one in which just a single mass transfer episode takes place - only those in which an eccentric binary are formed are likely to be detected by the planned gravitational wave mission eLISA, while more sensitive detectors would be able to detect the signals emitted in lateral collisions. On the other hand, the light curves (and the thermal neutrino emission) of these interactions are considerably different, producing both very powerful outbursts and low luminosity events. Finally, we also calculate the X-ray signature produced in the aftermath of those interactions for which a merger occurs. We find that the temporal evolution follows a power law with the same exponent found in the case of the mergers of two neutron stars, although the total energy released is smaller.

A secular variation of $G$ modifies the structure and evolutionary time scales of white dwarfs. Using an state-of-the-art stellar evolutionary code, an up-to-date pulsational code, and a detailed population synthesis code we demonstrate that the effects of a running $G$ are obvious both in the properties of individual white dwarfs, and in those of the white dwarf populations in clusters. Specifically, we show that the white dwarf evolutionary sequences depend on both the value of $\dot G/G$, and on the value of $G$ when the white dwarf was born. We show as well that the pulsational properties of variable white dwarfs can be used to constrain $\dot G/G$. Finally, we also show that the ensemble properties of of white dwarfs in clusters can also be used to set upper bounds to $\dot G/G$. Precisely, the tightest bound --- $\dot G/G \sim -1.8 10^{-12}$ yr$^{-1}$ --- is obtained studying the population of the old, metal-rich, well populated, open cluster NGC 6791. Less stringent upper limits can be obtained comparing the theoretical results obtained taking into account the effects of a running $G$ with the measured rates of change of the periods of two well studied pulsating white dwarfs, G117--B15A and R548. Using these white dwarfs we obtain $\dot G/G\sim -1.8\times 10^{-10}$ yr$^{-1}$, and $\dot G/G\sim -1.3\times 10^{-10}$ yr$^{-1}$, respectively, which although less restrictive than the previous bound, can be improved measuring the rate of change of the period of massive white dwarfs.

In old, dense stellar systems collisions of white dwarfs are a rather frequent phenomenon. Here we present the results of a comprehensive set of Smoothed Particle Hydrodynamics simulations of close encounters of white dwarfs aimed to explore the outcome of the interaction and the nature of the final remnants for different initial conditions. Depending on the initial conditions and the white dwarf masses, three different outcomes are possible. Specifically, the outcome of the interaction can be either a direct or a lateral collision or the interaction can result in the formation of an eccentric binary system. In those cases in which a collision occurs, the infalling material is compressed and heated such that the physical conditions for a detonation may be reached during the most violent phases of the merger. While we find that detonations occur in a significant number of our simulations, in some of them the temperature increase in the shocked region rapidly lifts degeneracy, leading to the quenching of the burning. We thus characterize under which circumstances a detonation is likely to occur as a result of the impact of the disrupted star on the surface of the more massive white dwarf. Finally, we also study which interactions result in bound systems, and in which ones the more massive white dwarf is also disrupted as a consequence of the dynamical interaction. The sizable number of simulations performed in this work allows to find how the outcome of the interaction depends on the distance at closest approach, and on the masses of the colliding white dwarfs, and which is the chemical pattern of the nuclearly processed material. Finally, we also discuss the influence of the masses and core chemical compositions of the interacting white dwarfs and the different kinds of impact in the properties of the remnants.

It has long been accepted that a possible mechanism for explaining the existence of magnetic white dwarfs is the merger of a binary white dwarf system, as there are viable mechanisms for producing sustainable magnetic fields within the merger product. However, the lack of rapid rotators in the magnetic white dwarf population has been always considered a problematic issue of this scenario. Smoothed Particle Hydrodynamics simulations show that in mergers in which the two white dwarfs have different masses a disc around the central compact object is formed. If the central object is magnetized it can interact with the disc through its magnetosphere. The torque applied by the disc changes the spin of the star, whereas the transferred angular momentum from the star to the disc determines the properties of the disc. In this work we build a model for the disc evolution under the effect of magnetic accretion, and for the angular momentum evolution of the star, which can be compared with the observations. Our model predicts that the magnetospheric interaction of magnetic white dwarfs with their discs results in a significant spin down, and we show that for magnetic white dwarfs with relatively strong fields (larger than 10 MG) the observed rotation periods of the magnetic white dwarf population can be reproduced. We also investigate whether turbulence can be sustained during the late phases of the evolution of the system. When a critical temperature below which turbulence is not sustained is introduced into the model, the periods of the three fast rotating, strongly magnetic, massive white dwarfs in the solar neighborhood are recovered.

It has long been accepted that a possible mechanism for explaining the existence of magnetic white dwarfs is the merger of a binary white dwarf system, as there are viable mechanisms for producing sustainable magnetism within the merger product. However, the lack of rapid rotators in the magnetic white dwarf population has been always considered a problematic issue of this scenario. In order to explain this discrepancy we build a model in which the interaction between the magnetosphere of the star and the disk induces angular momentum transfer. Our model predicts that the magnetospheric interaction of magnetic white dwarfs with their disks results in a significant spin down, and we show that the observed rotation period of REJ 0317-853, which is suggested to be a product of a double degenerate merger, can be reproduced.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the solar neighborhood.

The origin of R Coronae Borealis stars has been elusive for over 200 years. Currently, two theories for their formation have been presented. These are the Final Flash scenario, in which a dying asymptotic giant branch (AGB) star throws off its atmosphere to reveal the hydrogen poor, heavily processed material underneath, and the double degenerate scenario, in which two white dwarfs merge to produce a new star with renewed vigour. Some theories predict that the temperatures reached during the latter scenario would destroy any lithium originally present in the white dwarfs. The observed lithium content of some R Coronae Borealis stars, therefore, is often interpreted as an indication that the Final Flash scenario best describes their formation. In this paper, it is shown that lithium production can, indeed, occur in the merging of a helium white dwarf with a carbon-oxygen white dwarf if their chemical composition, particularly that of 3He, is fully considered. The production mechanism is described in detail, and the sensitivity of lithium production to the merger environment is investigated. Nucleosynthesis post-processing calculations of smoothed-particle hydrodynamics (SPH) tracer particles are performed to show that any lithium produced in these environments will be concentrated towards the cloud of material surrounding the R CrB star. Measurements of the lithium content of these clouds would, therefore, provide a valuable insight into the formation mechanism of R CrB stars.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the Solar neighborhood.

Type Ia supernovae are thought to be the outcome of the thermonuclear explosion of a white dwarf in a close binary system. Two possible scenarios, not necessarily incompatible, have been advanced. One assumes a white dwarf that accretes matter from a nondegenerate companion (the single degenerate scenario), the other assumes two white dwarfs that merge as a consequence of the emission of gravitational waves (the double degenerate scenario). The delay time distribution of star formation bursts strongly suggests that the DD scenario should be responsible of the late time explosions, but this contradicts the common wisdom that the outcome of the merging of two white dwarfs is an accretion induced collapse to a neutron star. In this contribution we review some of the most controversial issues of this problem.

We present a nucleosynthesis study of the merger of a 0.4 solar masses helium white dwarf with a 0.8 solar masses carbon-oxygen white dwarf, coupling the thermodynamic history of Smoothed Particle Hydrodynamics particles with a post-processing code. The resulting chemical abundance pattern, particularly for oxygen and fluorine, is in qualitative agreement with the observed abundances in R Coronae Borealis stars.

Many hydrogen deficient stars are characterised by surface abundance patterns that are hard to reconcile with conventional stellar evolution. Instead, it has been suggested that they may represent the result of a merger episode between a helium and a carbon-oxygen white dwarf. In this Letter, we present a nucleosynthesis study of the merger of a 0.4 M_sol helium white dwarf with a 0.8 M_sol carbon-oxygen white dwarf, by coupling the thermodynamic history of Smoothed Particle Hydrodynamics particles with a post-processing code. The resulting chemical abundance pattern, particularly for oxygen and fluorine, is in qualitative agreement with the observed abundances in R Coronae Borealis stars.

A variation of the gravitational constant over cosmological ages modifies the main sequence lifetimes and white dwarf cooling ages. Using an state-of-the-art stellar evolutionary code we compute the effects of a secularly varying G on the main sequence ages and, employing white dwarf cooling ages computed taking into account the effects of a running G, we place constraints on the rate of variation of Newton's constant. This is done using the white dwarf luminosity function and the distance of the well studied open Galactic cluster NGC 6791. We derive an upper bound G'/G ~ -1.8 10^{-12} 1/yr. This upper limit for the secular variation of the gravitational constant compares favorably with those obtained using other stellar evolutionary properties, and can be easily improved if deep images of the cluster allow to obtain an improved white dwarf luminosity function.

NGC 6791 is a well studied open cluster1 that it is so close to us that can be imaged down to very faint luminosities. The main sequence turn-off age (~8 Gyr) and the age derived from the termination of the white dwarf cooling sequence (~6 Gyr) are significantly different. One possible explanation is that as white dwarfs cool, one of the ashes of helium burning, 22Ne, sinks in the deep interior of these stars. At lower temperatures, white dwarfs are expected to crystallise and phase separation of the main constituents of the core of a typical white dwarf, 12C and 16O, is expected to occur. This sequence of events is expected to introduce significant delays in the cooling times, but has not hitherto been proven. Here we report that, as theoretically anticipated, physical separation processes occur in the cores of white dwarfs, solving the age discrepancy for NGC 6791.

The collision of two white dwarfs is a quite frequent event in dense stellar systems, like globular clusters and galactic nuclei. In this paper we present the results of a set of simulations of the close encounters and collisions of two white dwarfs. We use an up- to-date smoothed particle hydrodynamics code that incorporates very detailed input physics and an improved treatment of the artificial viscosity. Our simulations have been done using a large number of particles (~ 4 \times 10^5) and covering a wide range of velocities and initial distances of the colliding white dwarfs. We discuss in detail when the initial eccentric binary white dwarf survives the closest approach, when a lateral collision in which several mass transfer episodes occur is the outcome of the newly formed binary system, and which range of input parameters leads to a direct collision, in which only one mass transfer episode occurs. We also discuss the characteristics of the final configuration and we assess the possible observational signatures of the merger, such as the associated gravitational waveforms and the fallback luminosities. We find that the overall evolution of the system and the main characteristics of the final object agree with those found in previous studies. We also find that the fallback luminosities are close to 10^48 erg/s. Finally, we find as well that in the case of lateral and direct collisions the gravitational waveforms are characterized by large-amplitude peaks which are followed by a ring-down phase, while in the case in which the binary white dwarf survives the closest approach, the gravitational pattern shows a distinctive behavior, typical of eccentric systems.

We present the results of a set of high-resolution simulations of the merging process of two white dwarfs. In order to do so, we use an up-to-date Smoothed Particle Hydrodynamics code which incorporates very detailed input physics and an improved treatment of the artificial viscosity. Our simulations have been done using a large number of particles (4x10^5) and cover the full range of masses and chemical compositions of the coalescing white dwarfs. We also compare the time evolution of the system during the first phases of the coalescence with that obtained using a simplified treatment of mass transfer, we discuss in detail the characteristics of the final configuration, we assess the possible observational signatures of the merger, like the associated gravitational waveforms and the fallback X-ray flares, and we study the long-term evolution of the coalescence.