skip to main content


Title: Hydrodynamic simulations of white dwarf–white dwarf mergers and the origin of R Coronae Borealis stars
ABSTRACT

We study the properties of double white dwarf (DWD) mergers by performing hydrodynamic simulations using the new and improved adaptive mesh refinement code octo-tiger. We follow the orbital evolution of DWD systems of mass ratio $q=0.7$ for tens of orbits until and after the merger to investigate them as a possible origin for R Coronae Borealis (RCB) type stars. We reproduce previous results, finding that during the merger, the helium WD donor star is tidally disrupted within 20–80 min since the beginning of the simulation onto the accretor carbon–oxygen WD, creating a high temperature shell around the accretor. We investigate the possible helium burning in this shell and the merged object’s general structure. Specifically, we are interested in the amount of oxygen-16 dredged-up from the accretor to the hot shell and the amount of oxygen-18 produced. This is critical as the discovery of very low oxygen-16 to oxygen-18 ratios in RCB stars pointed out the merger scenario as a favourable explanation for their origin. A small amount of hydrogen in the donor may help keep the oxygen-16 to oxygen-18 ratios within observational bounds, even if moderate dredge-up from the accretor occurs. In addition, we perform a resolution study to reconcile the difference found in the amount of oxygen-16 dredge-up between smoothed-particle hydrodynamics and grid-based simulations.

 
more » « less
PAR ID:
10555327
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
535
Issue:
2
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 1914-1943
Size(s):
p. 1914-1943
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. RCB stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cmJ-band telescope with a 25 deg2camera that surveys 18,000 deg2of the northern sky (δ> −28°) at a cadence of 2 days. Using PGIRJ-band lightcurves for ∼60 million stars together with mid-IR colors from WISE, we selected a sample of 530 candidate RCB stars. We obtained near-IR spectra for these candidates and identified 53 RCB stars in our sample. Accounting for our selection criteria, we find that there are a total of350100+150RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8–5 × 10−3yr−1, consistent with observational and theoretical estimates of the He-CO WD merger rate. We searched for quasi-periodic pulsations in the PGIR lightcurves of RCB stars and present pulsation periods for 16 RCB stars. We also examined high-cadenced TESS lightcurves for RCB and the chemically similar, but dustless hydrogen-deficient carbon (dLHdC) stars. We find that dLHdC stars show variations on timescales shorter than RCB stars, suggesting that they may have lower masses than RCB stars. Finally, we identified 3 new spectroscopically confirmed and 12 candidate Galactic DY Per type stars—believed to be colder cousins of RCB star—doubling the sample of Galactic DY Per type stars.

     
    more » « less
  2. Abstract

    The origin of the bright and hard X-ray emission flux among theγCas subgroup of B-emission line (Be) stars may be caused by gas accretion onto an orbiting white dwarf (WD) companion. Such Be+WD binaries are the predicted outcome of a second stage of mass transfer from a helium star mass donor to a rapidly rotating mass gainer star. The stripped donor stars become small and hot white dwarfs that are extremely faint compared to their Be star companions. Here we discuss model predictions about the physical and orbital properties of Be+WD binaries, and we show that current observational results onγCas systems are consistent with the expected large binary frequency, companion faintness and small mass, and relatively high mass range of the Be star hosts. We determine that the companions are probably not stripped helium stars (hot subdwarf sdO stars), because these are bright enough to detect in ultraviolet spectroscopy, yet their spectroscopic signatures are not observed in studies ofγCas binaries. Interferometry of relatively nearby systems provides the means to detect very faint companions including hot subdwarf and cooler main-sequence stars. Preliminary observations of fiveγCas binaries with the CHARA Array interferometer show no evidence of the companion flux, leaving white dwarfs as the only viable candidates for the companions.

     
    more » « less
  3. Abstract

    The detonation of an overlying helium layer on a 0.8–1.1Mcarbon–oxygen (CO) white dwarf (WD) can detonate the CO WD and create a thermonuclear supernova (SN). Many authors have recently shown that when the mass of the He layer is low (≲0.03M), the ashes from its detonation minimally impact the spectra and light curve from the CO detonation, allowing the explosion to appear remarkably similar to Type Ia SNe. These new insights motivate our investigation of dynamical He shell burning and our search for a binary scenario that stably accumulates thermally unstable He shells in the 0.01–0.08Mrange, thick enough to detonate, but also often thin enough for minimal impact on the observables. We first show that our improved nonadiabatic evolution of convective He shell burning in this range of shell mass leads to conditions ripe for a He detonation. We also find that a stable mass transfer scenario with a high-entropy He WD donor of mass 0.15–0.25Myields the He shell masses needed to achieve the double detonations. This scenario also predicts that the surviving He donor leaves with a spatial velocity consistent with the unusual runaway object, D6-2. We find that hot He WD donors originate in common-envelope events when a 1.3–2.0Mstar fills its Roche lobe at the base of the red giant branch at orbital periods of 1–10 days with the CO WD.

     
    more » « less
  4. Abstract We seek signatures of the current experimental 12 C α , γ 16 O reaction rate probability distribution function in the pulsation periods of carbon–oxygen white dwarf (WD) models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium-burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental 12 C α , γ 16 O reaction rate probability distribution function. We find an average spread in relative period shifts of Δ P / P ≃ ±2% for the identified trapped g-modes over the ±3 σ uncertainty in the 12 C α , γ 16 O reaction rate probability distribution function—across the effective temperature range of observed DAV and DBV WDs and for different WD masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed WDs are typically given to six to seven significant figures of precision. This suggests that an astrophysical constraint on the 12 C α , γ 16 O reaction rate could, in principle, be extractable from the period spectrum of observed variable WDs. 
    more » « less
  5. The formation channels and predicted populations of double white dwarfs (DWDs) are important because a subset will evolve to be gravitational-wave sources and/or progenitors of Type Ia supernovae. Given the observed population of short-period DWDs, we calculate the outcomes of common envelope (CE) evolution when convective effects are included. For each observed white dwarf (WD) in a DWD system, we identify all progenitor stars with an equivalent proto-WD core mass from a comprehensive suite of stellar evolution models. With the second observed WD as the companion, we calculate the conditions under which convection can accommodate the energy released as the orbit decays, including (if necessary) how much the envelope must spin-up during the CE phase. The predicted post-CE final separations closely track the observed DWD orbital parameter space, further strengthening the view that convection is a key ingredient in CE evolution. 
    more » « less