Double White Dwarf Research Articles

CIRCUMSTELLAR ABSORPTION IN DOUBLE DETONATION TYPE Ia SUPERNOVAE

Upon formation, degenerate He core white dwarfs are surrounded by a radiative H-rich layer primarily supported by ideal gas pressure. In this Letter, we examine the effect of this H-rich layer on mass transfer in He+C/O double white dwarf binaries that will eventually merge and possibly yield a Type Ia supernova (SN Ia) in the double detonation scenario. Because its thermal profile and equation of state differ from the underlying He core, the H-rich layer is transferred stably onto the C/O white dwarf prior to the He core's tidal disruption. We find that this material is ejected from the binary system and sweeps up the surrounding interstellar medium hundreds to thousands of years before the SN Ia. The close match between the resulting circumstellar medium profiles and values inferred from recent observations of circumstellar absorption in SNe Ia gives further credence to the resurgent double detonation scenario.

Constraints on the pairing properties of main-sequence stars from observations of white dwarfs in binary systems

Observations of main sequence stars conducted over the last several decades have clearly shown that something like 50 per cent of stars of spectral types G and F occur in multiple systems. For earlier spectral types, the incidence of multiplicity is even higher. Thus, a volume limited sample of white dwarfs should reflect the percentage of binarity observed in stars of F to late B spectral types, which are their Main Sequence progenitors. However, a study of the local volume limited sample of white dwarfs (20 pc from the Sun) conducted by Holberg has shown that a white dwarf has a probability of only $\sim 32$ per cent of occurring in a binary system, in stark contrast to the observations of multiplicity of Main Sequence stars. Others studies have also led to the same conclusion. In this paper, we argue that the "hidden" white dwarfs are either in double white dwarf systems or in Sirius-like systems. We also show that the white dwarf progenitors of the SDSS white dwarf - M dwarf wide binaries are distributed according to Salpeter's IMF. However, they cannot be paired with secondary stars which are also drawn from this IMF, since such a pairing would produce a percentage of white dwarf-M dwarfs systems that is several times larger than observed.

Supernova Type Ia progenitors from merging double white dwarfs

The study of Type Ia supernovae (SNIa) has lead to greatly improved insights into many fields in astrophysics, however a theoretical explanation of the origin of these events is still lacking. We investigate the potential contribution to the SNIa rate from the population of merging double carbon-oxygen white dwarfs. We aim to develope a model that fits the observed SNIa progenitors as well as the observed close double white dwarf population. We differentiate between two scenarios for the common envelope (CE) evolution; the alpha-formalism based on the energy equation and the gamma-formalism that is based on the angular momentum equation. In one model we apply the alpha-formalism always. In the second model the gamma-formalism is applied, unless the binary contains a compact object or the CE is triggered by a tidal instability for which the alpha-formalism is used. The binary population synthesis code SeBa was used to evolve binary systems from the zero-age main sequence to the formation of double white dwarfs and subsequent mergers. SeBa has been thoroughly updated since the last publication of the content of the code. The limited sample of observed double white dwarfs is better represented by the simulated population using the gamma-formalism than the alpha-formalism. For both CE formalisms, we find that although the morphology of the simulated delay time distribution matches that of the observations within the errors, the normalisation and time-integrated rate per stellar mass are a factor 7-12 lower than observed. Furthermore, the characteristics of the simulated populations of merging double carbon-oxygen white dwarfs are discussed and put in the context of alternative SNIa models for merging double white dwarfs.

ORBITAL EVOLUTION OF COMPACT WHITE DWARF BINARIES

The newfound prevalence of extremely low mass (ELM, MHe < 0.2 M☉) helium white dwarfs (WDs) in tight binaries with more massive WDs has raised our interest in understanding the nature of their mass transfer. Possessing small (Menv ∼ 10−3 M☉) but thick hydrogen envelopes, these objects have larger radii than cold WDs and so initiate mass transfer of H-rich material at orbital periods of 6–10 minutes. Building on the original work of D'Antona et al., we confirm the 106 yr period of continued inspiral with mass transfer of H-rich matter and highlight the fact that the inspiraling direct-impact double WD binary HM Cancri likely has an ELM WD donor. The ELM WDs have less of a radius expansion under mass loss, thus enabling a larger range of donor masses that can stably transfer matter and become a He mass transferring AM CVn binary. Even once in the long-lived AM CVn mass transferring stage, these He WDs have larger radii due to their higher entropy from the prolonged H-burning stage.

COMMON PROPER-MOTION WIDE WHITE DWARF BINARIES SELECTED FROM THE SLOAN DIGITAL SKY SURVEY

Wide binaries made up of two white dwarfs (WDs) receive far less attention than their tight counterparts. However, our tests using the binary population synthesis code {\tt StarTrack} indicate that, for any set of reasonable initial conditions, there exists a significant observable population of double white dwarfs (WDWDs) with orbital separations of 10$^2$ to 10$^5$ AU. We adapt the technique of Dhital et al.\ to search for candidate common proper motion WD companions separated by $<10\amin$ around the $>$12,000 spectroscopically confirmed hydrogen-atmosphere WDs recently identified in the Sloan Digital Sky Survey. Using two techniques to separate random alignments from high-confidence pairs, we find nine new high-probability wide WDWDs and confirm three previously identified candidate wide WDWDs. This brings the number of known wide WDWDs to 45; our new pairs are a significant addition to the sample, especially at small proper motions ($<$200 mas/yr) and large angular separations ($>$10\asec). Spectroscopic follow-up and an extension of this method to a larger, photometrically selected set of SDSS WDs may eventually produce a large enough dataset for WDWDs to realize their full potential as testbeds for theories of stellar evolution.

DO R CORONAE BOREALIS STARS FORM FROM DOUBLE WHITE DWARF MERGERS?

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

SOME CONSTRAINTS ON THE LOWER MASS LIMIT FOR DOUBLE-DEGENERATE PROGENITORS OF TYPE Ia SUPERNOVAE

Recent theoretical and observational studies both argue that the merging of double carbon-oxygen white dwarfs (WDs) is responsible for at least some Type Ia supernovae (SNe Ia). Previous (standard) studies of the anticipated SN birthrate from this channel have assumed that the merger process is conservative and that the primary criterion for explosion is that the merged mass exceeds the Chandrasekhar mass. Han & Webbink (1999) demonstrated that mass transfer and merger in close double WDs will in many cases be non-conservative. Pakmor et al. (2011) further suggested that the merger process should be violent in order to initiate an explosion. We have therefore investigated how the SN Ia birthrate from the double-degenerate (DD) channel is affected by these constraints. Using the binary-star population-synthesis method, we have calculated the DD SN Ia birthrate under conservative and non-conservative approximations, and including lower mass and mass-ratio limits indicated by recent smoothed-particle-hydrodynamic calculations. The predicted DD SN Ia rate is significantly reduced by all of these constraints. With dynamical mass loss alone (violent merger) the birthrate is reduced to 56% of the conservative rate. Requiring that the mass ratio $q>2/3$ further reduces the birthrate to 18% that of the standard assumption. An upper limit of 0.0061 SNuM, or a Galactic rate of $4.6 \times 10^{-4}{\rm yr}^{-1}$, might be realistic.

A NUMERICAL METHOD FOR STUDYING SUPER-EDDINGTON MASS TRANSFER IN DOUBLE WHITE DWARF BINARIES

We present a numerical method for the study of double white dwarf (DWD) binary systems at the onset of super-Eddington mass transfer. We incorporate the physics of ideal inviscid hydrodynamical flow, Newtonian self-gravity, and radiation transport on a three-dimensional uniformly rotating cylindrical Eulerian grid. Care has been taken to conserve the key physical quantities such as angular momentum and energy. Our new method conserves total energy to a higher degree of accuracy than other codes that are presently being used to model mass transfer in DWD systems. We present the results of verification tests and simulate the first 20 + orbits of a binary system of mass ratio q = 0.7 at the onset of dynamically unstable direct impact mass transfer. The mass transfer rate quickly exceeds the critical Eddington limit by many orders of magnitude, and thus we are unable to model a trans-Eddington phase. It appears that radiation pressure does not significantly affect the accretion flow in the highly super-Eddington regime. An optically thick common envelope forms around the binary within a few orbits. Although this envelope quickly exceeds the spatial domain of the computational grid, the fraction of the common envelope that exceeds zero gravitational binding energy is extremely small, suggesting that radiation-driven mass loss is insignificant in this regime. It remains to be seen whether simulations that capture the trans-Eddington phase of such flows will lead to the same conclusion or show that substantial material gets expelled.

THE LONG-TERM EVOLUTION OF DOUBLE WHITE DWARF MERGERS

In this paper, we present a model for the long-term evolution of the merger of two unequal mass C/O white dwarfs (WDs). After the dynamical phase of the merger, magnetic stresses rapidly redistribute angular momentum, leading to nearly solid-body rotation on a viscous timescale of 1e-4 to 1 yr, long before significant cooling can occur. Due to heating during the dynamical and viscous phases, the less massive WD is transformed into a hot, slowly rotating, and radially extended envelope supported by thermal pressure. Following the viscous phase of evolution, the maximum temperature near the envelope base may already be high enough to begin off-center convective carbon-burning. If not, Kelvin-Helmholtz contraction of the inner region of the envelope on a thermal timescale of 1e3-1e4 yr compresses the base of the envelope, again yielding off-center burning. As a result, the long-term evolution of the merger remnant is similar to that seen in previous calculations: the burning shell diffuses inwards over ~1e4 yr, eventually yielding a high-mass O/Ne WD or a collapse to a neutron star. During the cooling and shell-burning phases, the merger remnant radiates near the Eddington limit. Given the double WD merger rate of a few per 1000 yr, tens of these ~1e38 erg/s sources should exist in a Milky Way-type galaxy. While the end result is similar to that of previous studies, the physical picture and the dynamical state of the matter in our model differ from previous work. Furthermore, remaining uncertainties related to the convective structure near the photosphere and mass loss during the thermal evolution may significantly affect our conclusions. Thus, future work within the context of the physical model presented here is required to better address the eventual fate of double WD mergers, including those for which one or both of the components is a He WD.

TIDALLY INDUCED APSIDAL PRECESSION IN DOUBLE WHITE DWARFS: A NEW MASS MEASUREMENT TOOL WITHLISA

Galactic interacting double white dwarfs (DWD) are guaranteed gravitational wave (GW) sources for the GW detector LISA, with more than 10^4 binaries expected to be detected over the mission's lifetime. Part of this population is expected to be eccentric, and here we investigate the potential for constraining the white dwarf (WD) properties through apsidal precession in these binaries. We analyze the tidal, rotational, and general relativistic contributions to apsidal precession by using detailed He WD models, where the evolution of the star's interior is followed throughout the cooling phase. In agreement with previous studies of zero-temperature WDs, we find that apsidal precession in eccentric DWDs can lead to a detectable shift in the emitted GW signal when binaries with cool (old) components are considered. This shift increases significantly for hot (young) WDs. We find that apsidal motion in hot (cool) DWDs is dominated by tides at orbital frequencies above ~10^{-4}Hz (10^{- 3}$Hz). The analysis of apsidal precession in these sources while ignoring the tidal component would lead to an extreme bias in the mass determination, and could lead us to misidentify WDs as neutron stars or black holes. We use the detailed WD models to show that for older, cold WDs, there is a unique relationship that ties the radius and apsidal precession constant to the WD masses, therefore allowing tides to be used as a tool to constrain the source masses.

ROTATIONAL DOPPLER BEAMING IN ECLIPSING BINARIES

In eclipsing binaries the stellar rotation of the two components will cause a rotational Doppler beaming during eclipse ingress and egress when only part of the eclipsed component is covered. For eclipsing binaries with fast spinning components this photometric analogue of the well-known spectroscopic Rossiter-McLaughlin effect can exceed the strength of the orbital effect. Example light curves are shown for a detached double white dwarf binary, a massive O-star binary and a transiting exoplanet case, similar to WASP-33b. Inclusion of the rotational Doppler beaming in eclipsing systems is a prerequisite for deriving the correct stellar parameters from fitting high quality photometric light curves and can be used to determine stellar obliquities as well as e.g. an independent measure of the rotational velocity in those systems that may be expected to be fully synchronized.

ON THE FORMATION OF DOUBLE WHITE DWARFS THROUGH STABLE MASS TRANSFER AND A COMMON ENVELOPE

Although several dozen double white dwarfs (DWDs) have been observed, for many the exact nature of the evolutionary channel(s) by which they form remains uncertain. The canonical explanation calls for the progenitor binary system to undergo two subsequent mass-transfer events, both of which are unstable and lead to a common envelope (CE). However, it has been shown that if both CE events obey the standard alpha-prescription (parametrizing energy loss), it is not possible to reproduce all of the observed systems. The gamma-prescription was proposed as an alternative to this description, instead parametrizing the fraction of angular momentum carried away in dynamical-timescale mass loss. In this paper, we analyze simultaneous energy and angular-momentum conservation, and show that the gamma-prescription cannot adequately describe a CE event for an arbitrary binary, nor can the first phase of mass loss always be understood in general as a dynamical-timescale event. We consider in detail the first episode of mass transfer in binary systems with initially low companion masses, with a primary mass in the range 1.0--1.3 solar masses and an initial mass ratio between the secondary and primary stars of 0.83--0.92. In these systems, the first episode of dramatic mass loss may be stable, non-conservative mass transfer. This strips the donor's envelope and dramatically raises the mass ratio; the considered progenitor binary systems can then evolve into DWDs after passing through a single CE during the second episode of mass loss. We find that such a mechanism reproduces the properties of the observed DWD systems which have an older component with less than approximately 0.46 solar masses and mass ratios between the younger and older WDs greater than 1.

SDSS J163030.58+423305.8: a 40-min orbital period detached white dwarf binary

ABSTRACT We report the discovery of a new detached, double white dwarf (WD) system with an orbital period of 39.8 min. We targeted SDSS J163030.58+423305.8 (hereafter J1630) as part of our radial velocity programme to search for companions around low-mass WDs using the 6.5-m MMT. We detect peak-to-peak radial velocity variations of 576 km s−1. The mass function and optical photometry rule out main-sequence companions. In addition, no millisecond pulsar companions are detected in radio observations. Thus the invisible companion is most likely another WD. Unlike the other 39-min binary SDSS J010657.39−100003.3, follow-up high-speed photometric observations of J1630 obtained at the McDonald 2.1-m telescope do not show significant ellipsoidal variations, indicating a higher primary mass and smaller radius. The absence of eclipses constrain the inclination angle to i≤ 82°. J1630 contains a pair of WDs, 0.3 M⊙ primary +≥0.3 M⊙ invisible secondary, at a separation of ≥0.32 R⊙. The two WDs will merge in less than 31 Myr. Depending on the core composition of the companion, the merger will form either a single core He-burning subdwarf star or a rapidly rotating massive WD. The gravitational wave strain from J1630 is detectable by instruments like the Laser Interferometer Space Antenna (LISA) within the first year of operation.

Evolutionary models for double helium white dwarf mergers and the formation of helium-rich hot subdwarfs

Recentsurveyshavedemonstratedtheexistenceofseveralshort-periodbinarysystemscontaining two white dwarfs. Following orbital decay by gravitational-wave radiation, such binaries are expected to merge at a rate of two or three per thousand years per galaxy. The consequences of such a merger depend on the individual white dwarf masses, but are believed to include helium-rich subdwarfs, RCrB stars, extreme helium stars and also AMCVn systems and possibly Type Ia supernovae. Whilst the hydrodynamics of the merger process remains difficult to compute, it is possible to compute the evolution of a double white dwarf merger following the destruction of one component. In this paper, we describe the evolution following the merger of two helium white dwarfs. We examine three sets of assumptions concerning the distribution of debris material between a disc and a corona. Our results demonstrate that a model comprising both fast accretion to form a (hot) corona andslow accretion from a (cold) debris disc can reproduce the observed distribution of helium-rich subdwarfs in terms of their surface temperatures, gravities, nitrogen and carbon

Delay times and rates for Type Ia supernovae and thermonuclear explosions from double-detonation sub-Chandrasekhar mass models

We present theoretical delay times and rates of thermonuclear explosions that are thought to produce Type Ia supernovae, including the double-detonation sub-Chandrasekhar mass model, using the population synthesis binary evolution code StarTrack. If detonations of sub-Chandrasekhar mass carbon-oxygen white dwarfs following a detonation in an accumulated layer of helium on the white dwarf's surface ("double-detonation" models) are able to produce thermonuclear explosions which are characteristically similar to those of SNe Ia, then these sub-Chandrasekhar mass explosions may account for at least some substantial fraction of the observed SN Ia rate. Regardless of whether all double-detonations look like 'normal' SNe Ia, in any case the explosions are expected to be bright and thus potentially detectable. Additionally, we find that the delay time distribution of double-detonation sub-Chandrasekhar mass SNe Ia can be divided into two distinct formation channels: the 'prompt' helium-star channel with delay times <500 Myr (~10% of all sub-Chandras), and the 'delayed' double white dwarf channel, with delay times >800 Myr spanning up to a Hubble time (~90%). These findings coincide with recent observationally-derived delay time distributions which have revealed that a large number of SNe Ia are prompt with delay times <500 Myr, while a significant fraction also have delay times spanning ~1 Gyr to a Hubble time.

DETECTION OF A WHITE DWARF COMPANION TO THE WHITE DWARF SDSSJ125733.63+542850.5

SDSSJ125733.63+542850.5 (hereafter SDSSJ1257+5428) is a compact white dwarf binary from the Sloan Digital Sky Survey that exhibits high-amplitude radial velocity variations on a period of 4.56 hours. While an initial analysis suggested the presence of a neutron star or black-hole binary companion, a follow-up study concluded that the spectrum was better understood as a combination of two white dwarfs. Here we present optical spectroscopy and ultraviolet fluxes which directly reveal the presence of the second white dwarf in the system. SDSSJ1257+5428's spectrum is a composite, dominated by the narrow-lined spectrum from a cool, low gravity white dwarf (Teff ~6300K, log g = 5 to 6.6) with broad wings from a hotter, high-mass white dwarf companion (11,000 to 14,000K; ~ 1Msun). The high-mass white dwarf has unusual line profiles which lack the narrow central core to Halpha that is usually seen in white dwarfs. This is consistent with rapid rotation with v sin i = 500 to 1750km/s, although other broadening mechanisms such as magnetic fields, pulsations or a helium-rich atmosphere could also be contributory factors. The cool component is a puzzle since no evolutionary model matches its combination of low gravity and temperature. Within the constraints set by our data, SDSSJ1257+5428 could have a total mass greater than the Chandrasekhar limit and thus be a potential Type Ia supernova progenitor. However, SDSSJ1257+5428's unusually low mass ratio q ~ 0.2 suggests that it is more likely that it will evolve into an accreting double white dwarf (AM CVn star).

Constraining White Dwarf Masses Via Apsidal Precession in Eccentric Double White Dwarf Binaries

AbstractGalactic short period double white dwarfs (DWD) are guaranteed gravitational wave (GW) sources for the next generation of space-based interferometers sensitive to low-frequency GWs (10−4− 1 Hz). Here we investigate the possibility of constraining the white dwarf (WD) properties through measurements of apsidal precession in eccentric binaries. We analyze the general relativistic (GR), tidal, and rotational contributions to apsidal precession by using detailed He WD models. We find that apsidal precession can lead to a detectable shift in the emitted GW signal, the effect being stronger (weaker) for binaries hosting hot (cool) WDs. We find that in hot (cool) DWDs tides dominate the precession at orbital frequencies above ~0.01 mHz (~1 mHz). Analyzing the apsidal precession of these sources only accounting for GR would potentially lead to an extreme overestimate of the component masses. Finally, we derive a relation that ties the radius and apsidal precession constant of cool WD components to their masses, therefore allowing tides to be used as an additional mass measurement tool.

Constraining Type Ia Supernova Progenitors

AbstractWe present observational and theoretical studies constraining Type Ia supernova progenitors.First, we use a new observational technique to show that “prompt” SNe Ia that trace star-formation on cosmic timescales exhibit a significant delay time of 200-500 million years. This implies that either the majority of SNe Ia companion stars have main-sequence masses less than three solar masses, or that most SNe Ia arise from double-white dwarf binaries.Second we present a comprehensive study of white dwarf collisions as an avenue for creating SNe Ia. Using a smooth particle hydrodynamics code with a 13-isotope nuclear network, we show that several combinations of white dwarf masses and impact parameters produce enough 56Ni to result in luminosities ranging from those of sub-luminous to super-luminous SNe Ia, depending on the parameters of the collision.Finally, we conduct a simulation survey of double-degenerate white dwarf mergers with varying mass combinations. Unlike previous works, we do not add detonations by hand to our simulations, and we do not find any thermonuclear explosions during the mergers. Instead, all but one of our simulations forms a cold, degenerate core surrounded by a hot disk, while our least massive pair of stars forms only a hot disk. We characterize the remnants by core mass, rotational velocity, and half-mass radius, and discuss how we will evolve them further with simulations that incorporate dissipative processes. Such simulations may indeed lead to double-degenerate Type Ia explosions that occur many orbits after the mergers themselves.

A DEEPLY ECLIPSING DETACHED DOUBLE HELIUM WHITE DWARF BINARY

Using Liverpool Telescope+RISE photometry we identify the 2.78 hour period binary star CSS 41177 as a detached eclipsing double white dwarf binary with a 21,100K primary star and a 10,500K secondary star. This makes CSS 41177 only the second known eclipsing double white dwarf binary after NLTT 11748. The 2 minute long primary eclipse is 40% deep and the secondary eclipse 10% deep. From Gemini+GMOS spectroscopy we measure the radial velocities of both components of the binary from the H{\alpha} absorption line cores. These measurements, combined with the light curve information, yield white dwarf masses of M1 = 0.283\pm0.064Msun and M2 = 0.274\pm0.034Msun, making them both helium core white dwarfs. As an eclipsing, double-lined spectroscopic binary CSS 41177 is ideally suited to measuring precise, model-independent masses and radii. The two white dwarfs will merge in roughly 1.1 Gyr to form a single sdB star.

Double white dwarf mergers and elemental surface abundances in extreme helium and R Coronae Borealis stars

The surface abundances of extreme helium (EHe) and R Coronae Borealis (RCB) stars are discussed in terms of the merger of a carbon-oxygen white dwarf with a helium white dwarf. The model is expressed as a linear mixture of the individual layers of both constituent white dwarfs, taking account of the specific evolution of each star. In developing this recipe from previous versions, particular attention has been given to the inter-shell abundances of the asymptotic giant branch star which evolved to become the carbon-oxygen white dwarf. Thus the surface composition of the merged star is estimated as a function of the initial mass and metallicity of its progenitor. The question of whether additional nucleosynthesis occurs during the white dwarf merger has been examined. The high observed abundances of carbon and oxygen must either originate by dredge-up from the core of the carbon-oxygen white dwarf during a cold merger or be generated directly by alpha-burning during a hot merger. The presence of large quantities of O18 may be consistent with both scenarios, since a significant O18 pocket develops at the carbon/helium boundary in a number of our post-AGB models. The production of fluorine, neon and phosphorus in the AGB intershell produces n overabundance at the surface of the merged stars, but generally not in sufficient quantity. However, the evidence for an AGB origin for these elements points to progenitor stars with initial masses in the range 1.9 - 3 solar masses. There is not yet sufficient information to discriminate the origin (fossil or prompt) of all the abundance anomalies observed in EHe and RCB stars. Further work is required on argon and s-process elements in the AGB intershell, and on the predicted yields of all elements from a hot merger.