White Dwarf Mergers Research Articles

Lower-mass-gap Black Holes in Dense Star Clusters

The existence of compact stellar remnants in the mass range 2–5 M ⊙ has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5 M ⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte Carlo N-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc−3 yr−1.

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.

Blue large-amplitude pulsators formed from the merger of low-mass white dwarfs

Blue large-amplitude pulsators (BLAPs) are a recently discovered group of hot stars pulsating in radial modes. Their origin needs to be explained, and several scenarios for their formation have already been proposed. We investigate whether BLAPs can originate as the product of a merger of two low-mass white dwarfs (WDs) and estimate how many BLAPs can be formed in this evolutionary channel. We used the Modules for Experiments in Stellar Astrophysics (MESA) code to model the merger of three different double extremely low-mass (DELM) WDs and the subsequent evolution of the merger product. We also performed a population synthesis of Galactic DELM WDs using the COSMIC code. We find that BLAPs can be formed from DELM WDs provided that the total mass of the system ranges between 0.32 and 0.7\,M$_ BLAPs born in this scenario either do not have any thermonuclear fusion at all or show off-centre He burning. The final product evolves to hot subdwarfs and eventually finishes its evolution either as a cooling He WD or a hybrid He/CO WD. The merger products become BLAPs only a few thousand years after coalescence, and it takes them 20\,--\,70 thousand years to pass the BLAP region. We found the instability of the fundamental radial mode to be in fair agreement with observations, but we also observed instability of the radial first overtone. The calculated evolutionary rates of period change can be both positive and negative. From the population synthesis, we found that up to a few hundred BLAPs born in this scenario can exist at present in the Galaxy. Given the estimated number of BLAPs formed in the studied DELM WD merger scenario, there is a good chance to observe BLAPs that originated through this scenario. Since strong magnetic fields can be generated during mergers, this scenario could lead to the formation of magnetic BLAPs. This fits well with the discovery of two likely magnetic BLAPs whose pulsations can be explained in terms of the oblique rotator model.

Large-scale ordered magnetic fields generated in mergers of helium white dwarfs

Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1<!PCT!>$. We present two high-resolution magnetohydrodynamical (MHD) simulations of the merger of two helium white dwarfs in a binary system with the same total mass of $0.6\,M_ We analysed an equal-mass merger with two $0.3\,M_ white dwarfs, and an unequal-mass merger with white dwarfs of $0.25\,M_ and $0.35\,M_ We simulated the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations. We found efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at a similar strength for both simulations. We then identified a second phase of magnetic field amplification in both merger remnants that happens on a timescale of several tens of rotational periods of the merger remnant. This phase generates a large-scale ordered azimuthal field via a large-scale dynamo driven by the magneto-rotational instability. Finally, we speculate that in the unequal-mass merger remnant, helium burning will initially start in a shell around a cold core, rather than in the centre. This forms a convection zone that coincides with the region that contains most of the magnetic energy, and likely destroys the strong, ordered field. Ohmic resistivity might then quickly erase the remaining small-scale field. Therefore, the mass ratio of the initial merger could be the selecting factor that decides if a merger remnant will stay highly magnetised long after the merger.

Modelling hydrogen-deficient carbon stars in mesa – the effects of total mass and mass ratio

ABSTRACT Hydrogen-deficient carbon (HdC) stars are rare, low-mass, chemically peculiar, supergiant variables believed to be formed by a double white dwarf (DWD) merger, specifically of a carbon/oxygen- (CO-) and a helium-white dwarf (He-WD). They consist of two subclasses – the dust-producing R Coronae Borealis (RCB) variables and their dustless counterparts the dustless HdCs (dLHdCs). Additionally, there is another, slightly cooler set of potentially related carbon stars, the DY Persei type variables which have some, but not conclusive, evidence of hydrogen-deficiency. Recent works have begun to explore the relationship between these three classes of stars, theorizing that they share an evolutionary pathway (a DWD merger) but come from different binary populations, specifically different total masses ($M_\text{tot}$) and mass ratios (q). In this work, we use the mesa modelling framework that has previously been used to model RCB stars and vary the merger parameters, $M_\text{tot}$ and q, to explore how those parameters affect the abundances, temperatures, and luminosities of the resultant post-merger stars. We find that lower $M_\text{tot}$ and larger q’s both decrease the luminosity and temperatures of post-merger models to the region of the Hertzsprung–Russell diagram populated by the dLHdCs. These lower $M_\text{tot}$ and larger q models also have smaller oxygen isotopic ratios ($^{16}$O/$^{18}$O) which is consistent with recent observations of dLHdCs compared to RCBs. None of the models generated in this work can explain the existence of the DY Persei type variables, however this may arise from the assumed metallicity of the models.

An Infrared Census of R Coronae Borealis Stars II—Spectroscopic Classifications and Implications for the Rate of Low-mass White Dwarf Mergers

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 cm J-band telescope with a 25 deg2 camera that surveys 18,000 deg2 of the northern sky (δ > −28°) at a cadence of 2 days. Using PGIR J-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 of ≈350−100+150 RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8–5 × 10−3 yr−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.

Light Curves of the Explosion of ONe White Dwarf + CO White Dwarf Merger Remnant and Type Icn Supernovae

Type Icn supernovae (SNe Icn) are a newly detected, rare subtype of interacting stripped-envelope supernovae that show narrow P Cygni lines of highly ionized carbon, oxygen, and neon in their early spectra due to the interactions of the SNe ejecta with dense hydrogen- and helium-deficient circumstellar material (CSM). It has been suggested that SNe Icn may have multiple progenitor channels, such as the explosion of carbon-rich Wolf–Rayet stars or the explosion of stripped-envelope SNe, which undergo binary interactions. Among the SNe Icn, SN 2019jc shows unique properties, and previous work inferred that it may stem from the ultrastripped supernova, but other possibilities still exist. In this work, we aim to simulate the light curves from the explosions of oxygen-neon and carbon-oxygen double white dwarf (WD) merger remnants and to further investigate whether the corresponding explosions can appear as some particular SNe Icn. We generate the light curves from the explosive remnants and analyze the influence of different parameters on the light curves, such as the ejecta mass, explosion energy, mass of 56Ni, and CSM properties. Comparing our results with some SNe Icn, we found that the light curves from the explosions of double WD merger remnants can explain the observable properties of SN 2019jc, from which we infer that this special SN Icn may have a different progenitor. Our results indicate that double WD merger may be an alternative model in producing at least one of the SNe Icn.

Ground-based and JWST Observations of SN 2022pul. II. Evidence from Nebular Spectroscopy for a Violent Merger in a Peculiar Type Ia Supernova

We present an analysis of ground-based and JWST observations of SN 2022pul, a peculiar “03fg-like” (or “super-Chandrasekhar”) Type Ia supernova (SN Ia), in the nebular phase at 338 days postexplosion. Our combined spectrum continuously covers 0.4–14 μm and includes the first mid-infrared spectrum of a 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. A strong, broad, centrally peaked [Ne ii] line at 12.81 μm was previously predicted as a hallmark of “violent merger” SN Ia models, where dynamical interaction between two sub-M Ch white dwarfs (WDs) causes disruption of the lower-mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to reproduce the optical iron emission better, and add mass in the innermost region (<2000 km s−1) to account for the observed narrow [O i] λ λ6300, 6364 emission. A violent WD–WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SNe Ia.

A low-mass helium star progenitor model for the Type Ibn SN 2020nxt

ABSTRACT A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterized by both rapidly evolving light curves and persistent narrow He i lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor system and mass-loss mechanism. In this paper, we present multiwavelength data of the Type Ibn SN 2020nxt, including HST/STIS ultraviolet spectra. We fit the data with recently updated CMFGEN models designed to handle configurations for SNe Ibn. The UV coverage yields strong constraints on the energetics and, when combined with the CMFGEN models, offer new insight on potential progenitor systems. We find the most successful model is a ≲4 M⊙ helium star that lost its $\sim 1\, {\rm M}_\odot$ He-rich envelope in the years preceding core collapse. We also consider viable alternatives, such as a He white dwarf merger. Ultimately, we conclude at least some SNe Ibn do not arise from single, massive (&amp;gt;30 M⊙) Wolf–Rayet-like stars.

Radio emission from SN 1181 hosting a white dwarf merger product

Abstract The remnant of the historical supernova 1181 is claimed to be associated with a white dwarf merger remnant J005311. The supernova remnant (SNR) shock, and a termination shock expected to be formed by the intense wind of J005311, are potential sites for radio emission via synchrotron emission from shock-accelerated electrons. In this paper, we estimate the radio emission from these two shocks, and find the peak radio flux to be 0.1–10 mJy (at 0.01–1 GHz) in the outer SNR shock and 0.01–0.1 mJy (at 1–10 GHz) in the inner termination shock. We also search for radio emission from this source in the archival data of the Karl G. Jansky Very Large Array (VLA) Sky Survey at 3 GHz, the NRAO VLA Sky Survey at 1.4 GHz and the Canadian Galactic Plane Survey at 408 MHz, finding no significant detection. While targeted observations with higher sensitivity are desired, we particularly encourage those at higher frequency and angular resolution to probe the inner termination shock and its evolution.

White Dwarf Merger Remnants: The DAQ Subclass

Four years after the discovery of a unique DAQ white dwarf with a hydrogen-dominated and carbon-rich atmosphere, we report the discovery of four new DAQ white dwarfs, including two that were not recognized properly in the literature. We find all five DAQs in a relatively narrow mass and temperature range of M = 1.14–1.19 M ⊙ and T eff = 13,000–17,000 K. In addition, at least two show photometric variations due to rapid rotation with ≈10 minute periods. All five are also kinematically old, but appear photometrically young, with estimated cooling ages of about 1 Gyr based on standard cooling tracks, and their masses are roughly twice the mass of the most common white dwarfs in the solar neighborhood. These characteristics are smoking gun signatures of white dwarf merger remnants. Comparing the DAQ sample with warm DQ white dwarfs, we demonstrate that there is a range of hydrogen abundances among the warm DQ population and that the distinction between DAQ and warm DQ white dwarfs is superficial. We discuss the potential evolutionary channels for the emergence of the DAQ subclass, suggesting that DAQ white dwarfs are trapped on the crystallization sequence and may remain there for a significant fraction of the Hubble time.

What Powered the Kilonova-like Emission after GRB 230307A in the Framework of a Neutron Star–White Dwarf Merger?

The second brightest gamma-ray burst, GRB 230307A (with a duration T 90 ∼ 40 s), exhibited characteristics indicative of a magnetar engine during the prompt emission phase. Notably, a suspected kilonova was identified in its follow-up optical and infrared observations. Here we propose that the origin of GRB 230307A is a neutron star–white dwarf (NS–WD) merger as this could naturally explain the long duration and the large physical offset from the center of its host galaxy. In the framework of such an NS–WD merger event, the late-time kilonova-like emission is very likely to be powered by the spin-down of the magnetar and the radioactive decay of 56Ni, rather than by the decay of r-process elements as these heavy elements may not easily be synthesized in an NS-WD merger. It is demonstrated that the above scenario can be supported by our fit to the late-time observational data, where a mass of ∼10−3 M ⊙ 56Ni is involved in the ejecta of a mass of ∼0.1 M ⊙. Particularly, the magnetar parameters required by the fit are consistent with those derived from the early X-ray observation.

The Optically Thick Rotating Magnetic Wind from a Massive White Dwarf Merger Product. II. Axisymmetric Magnetohydrodynamic Simulations

We numerically construct a series of axisymmetric rotating magnetic wind solutions, aiming at exploring the observation properties of massive white dwarf (WD) merger remnants with a strong magnetic field, a fast spin, and an intense mass loss, as inferred for WD J005311. We investigate the magnetospheric structure and the resultant spin-down torque exerted to the merger remnant with respect to the surface magnetic flux Φ*, spin angular frequency Ω* and the mass-loss rate Ṁ . We confirm that the wind properties for σ≡Φ*2Ω*2/Ṁvesc3≳1 significantly deviate from those of the spherical Parker wind, where v esc is the escape velocity at stellar surface. For such a rotating magnetic wind sequence, we find (i) a quasiperiodic mass eruption triggered by magnetic reconnection along with the equatorial plane and (ii) a scaling relation for the spin-down torque T≈(1/2)×ṀΩ*R*2σ1/4 . We apply our results to discuss the spin-down evolution and wind anisotropy of massive WD merger remnants, the latter of which could be probed by a successive observation of WD J005311 using Chandra.

Neutron star mergers as the dominant contributor to the production of heavy r-process elements

ABSTRACT The discovery of the radioactively powered kilonova AT2017gfo, associated with the short-duration gamma-ray burst GRB 170817A and the gravitational wave source GW170817, has provided the first direct evidence supporting binary neutron star mergers as crucial astrophysical sites for the synthesis of heavy elements beyond iron through r-process nucleosysthesis in the universe. However, recent identification of kilonovae following long-duration gamma-ray bursts, such as GRB 211211A and GRB 230307A, has sparked discussions about the potential of neutron star–white dwarf mergers to also produce neutron-rich ejecta and contribute to the production of heavy r-process elements. In this work, we estimate the contribution of binary neutron star mergers to the total mass of r-process elements in the Milky Way and investigate the possibility of neutron star–white dwarf mergers as alternative astrophysical sites for r-process nucleosynthesis through an analysis of the total mass of the r-process elements in the Milky Way. Our results reveal that binary neutron star mergers can sufficiently account for the Galactic heavy r-process elements, suggesting that these events are the dominant contributor to the production of heavy r-process elements in the Milky Way. Considering the total mass of r-process elements in the Milky Way and the higher occurrence rate of neutron star–white dwarf mergers, it is unlikely that such mergers can produce a significant amount of neutron-rich ejecta, with the generated mass of r-process elements being lower than $0.005\, {\rm M}_{\odot }$.

Prospects for detecting neutron star–white dwarf mergers with decihertz gravitational-wave observatories

ABSTRACT Based on different neutron star–white dwarf (NS–WD) population models, we investigate the prospects of gravitational-wave (GW) detections for NS–WD mergers, with the help of early warnings from two space-borne decihertz GW observatories, DO-Optimal and DECIGO. We not only give quick assessments of the GW detection rates for NS–WD mergers with the two decihertz GW detectors, but also report systematic analyses on the characteristics of GW-detectable merger events using the method of Fisher matrix. With a sufficient 1-d early-warning time, the yearly GW detection number for DO-Optimal is in the range of (1.5–1.9) × 103, while it is (3.3–4.6) × 104 for DECIGO. More importantly, our results show that most NS–WD mergers can be localized with an uncertainty of $\mathcal {O}(10^{-2})\, \mathrm{deg}^2$. Given the NS–WD merger as a possible origin for a peculiar long-duration gamma-ray burst, GRB 211211A, followed with kilonova-like emissions, we further suggest that the GW early-warning detection would allow future electromagnetic telescopes to get prepared to follow up transients after some special NS–WD mergers. Based on our analyses, we emphasize that such a feasible ‘wait-for’ pattern can help to firmly identify the origin of GRB 211211A-like events in the future and bring excellent opportunities for the multimessenger astronomy.

Rotation plays a role in the generation of magnetic fields in single white dwarfs

ABSTRACT Recent surveys of close white dwarf binaries as well as single white dwarfs have provided evidence for the late appearance of magnetic fields in white dwarfs, and a possible generation mechanism, a crystallization and rotation-driven dynamo has been suggested. A key prediction of this dynamo is that magnetic white dwarfs rotate, at least on average, faster than their non-magnetic counterparts and/or that the magnetic field strength increases with rotation. Here we present rotation periods of ten white dwarfs within 40 pc measured using photometric variations. Eight of the light curves come from TESS observations and are thus not biased towards short periods, in contrast to most period estimates that have been reported previously in the literature. These TESS spin periods are indeed systematically shorter than those of non-magnetic white dwarfs. This means that the crystallization and rotation-driven dynamo could be responsible for a fraction of the magnetic fields in white dwarfs. However, the full sample of magnetic white dwarfs also contains slowly rotating strongly magnetic white dwarfs which indicates that another mechanism that leads to the late appearance of magnetic white dwarfs might be at work, either in addition to or instead of the dynamo. The fast-spinning and massive magnetic white dwarfs that appear in the literature form a small fraction of magnetic white dwarfs, and probably result from a channel related to white dwarf mergers.

Self-consistent magnetohydrodynamic simulation of jet launching in a neutron star – white dwarf merger

The merger of a white dwarf (WD) and a neutron star (NS) is a relatively common event that produces an observable electromagnetic signal. Furthermore, the compactness of these stellar objects makes them an interesting candidate for gravitational wave (GW) astronomy, potentially being in the frequency range of LISA and other missions. To date, three-dimensional simulations of these mergers have not fully modeled the WD disruption or have used lower resolutions and have not included magnetic fields even though they potentially shape the evolution of the merger remnant. In this work, we simulated the merger of a 1.4 M⊙ NS with a 1 M⊙ carbon oxygen WD in the magnetohydrodynamic moving mesh code AREPO. We find that the disruption of the WD forms an accretion disk around the NS, and the subsequent accretion by the NS powers the launch of strongly magnetized, mildly relativistic jets perpendicular to the orbital plane. Although the exact properties of the jets could be altered by unresolved physics around the NS, the event could result in a transient with a larger luminosity than kilonovae. We discuss possible connections to fast blue optical transients (FBOTs) and long-duration gamma-ray bursts. We find that the frequency of GWs released during the merger is too high to be detectable by the LISA mission, but suitable for deci-hertz observatories such as LGWA, BBO, or DECIGO.

Ground-based and JWST Observations of SN 2022pul. I. Unusual Signatures of Carbon, Oxygen, and Circumstellar Interaction in a Peculiar Type Ia Supernova

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground- and space-based follow-up campaign to characterize SN 2022pul, a super-Chandrasekhar mass SN Ia (alternatively “03fg-like” SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon–oxygen (C/O)-rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity (M B = −18.9 mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak B-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O i] λ λ6300, 6364 (FWHM ≈ 2000 km s−1), strong, broad emission from [Ca ii] λ λ7291, 7323 (FWHM ≈ 7300 km s−1), and a rapid Fe iii to Fe ii ionization change. Finally, we present the first ever optical-to-MIR nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with T ≈ 500 K), combined with prominent [O i] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within C/O-rich CSM.

On the Optical Transients from Double White-dwarf Mergers

Double white dwarf (DWD) mergers are relevant astrophysical sources expected to produce massive, highly magnetized white dwarfs (WDs), supernovae (SNe) Ia, and neutron stars (NSs). Although they are expected to be numerous sources in the sky, their detection has evaded the most advanced transient surveys. This article characterizes the optical transient expected from DWD mergers in which the central remnant is a stable (sub-Chandrasekhar) WD. We show that the expansion and cooling of the merger’s dynamical ejecta lead to an optical emission peaking at 1–10 days postmerger, with luminosities of 1040–1041 erg s−1. We present simulations of the light curves, spectra, and the color evolution of the transient. We show that these properties, together with the estimated rate of mergers, are consistent with the absence of detection, e.g., by the Zwicky Transient Facility. More importantly, we show that the Legacy Survey of Space and Time of the Vera C. Rubin Observatory will likely detect a few/several hundred per year, opening a new window to the physics of WDs, NSs, and SNe Ia.

The distribution, kinematics, and luminosities of extreme helium stars as probes of their origin and evolution

ABSTRACT Hydrogen-deficient stars include the cool R CrB variable (RCBs) and hydrogen-deficient carbon (HdCs) giants through extreme helium stars (EHes) to the very hot helium-rich subdwarfs (He-sdO and O(He) stars) and white dwarfs. With surfaces rich in helium, nitrogen, and carbon, their origins have been identified with the merger of two white dwarfs. Using Gaia to focus on the EHes, we aim to identify progenitor populations and test the evolution models. Gaia DR3 measurements and ground-based radial velocities have been used to compute Galactic orbits using galpy. Each orbit has been classified by population; EHe stars are found in all of the thin disc, thick disc, halo, and bulge, as are RCB, HdC, and He-sdO stars. Spectral energy distributions were constructed for all EHes, to provide angular diameters, and hence radii and luminosities. The EHes fall into two luminosity groups divided at $L \approx 2500 \, {\rm L_{\odot }}$. This supports theory for the origin of EHes, and is the strongest confirmation so far in terms of luminosity. The lower luminosity EHes correspond well with the post-merger evolution of a double helium white dwarf binary. Likewise, the higher luminosity EHes match the post-merger evolution of a carbon/oxygen plus helium white dwarf binary. In terms of parent populations, current models predict that double white dwarf mergers should occur in all Galactic populations, but favour mergers arising from recent star formation (i.e. thin disc), whereas the statistics favour an older epoch (i.e. thick disc).