skip to main content


Title: Predicting the LISA white dwarf binary population in the Milky Way with cosmological simulations
ABSTRACT

White dwarf binaries with orbital periods below 1 h will be the most numerous sources for the space-based gravitational wave detector Laser Interferometer Space Antenna (LISA). Based on thousands of individually resolved systems, we will be able to constrain binary evolution and provide a new map of the Milky Way and its close surroundings. In this paper we predict the main properties of populations of different types of detached white dwarf binaries detected by LISA over time. For the first time, we combine a high-resolution cosmological simulation of a Milky Way-mass galaxy (taken from the FIRE project) with a binary population synthesis model for low- and intermediate-mass stars. Our Galaxy model therefore provides a cosmologically realistic star formation and metallicity history for the Galaxy and naturally produces its different components such as the thin and thick disc, the bulge, the stellar halo, and satellite galaxies and streams. Thanks to the simulation, we show how different Galactic components contribute differently to the gravitational wave signal, mostly due to their typical age and distance distributions. We find that the dominant LISA sources will be He–He double white dwarfs (DWDs) and He–CO DWDs with important contributions from the thick disc and bulge. The resulting sky map of the sources is different from previous models, with important consequences for the searches for electromagnetic counterparts and data analysis. We also emphasize that much of the science-enabling information regarding white dwarf binaries, such as the chirp mass and the sky localization, becomes increasingly rich with long observations, including an extended mission up to 8 yr.

 
more » « less
NSF-PAR ID:
10124586
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
490
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
p. 5888-5903
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Compact-object binaries including a white dwarf component are unique among gravitational-wave sources because their evolution is governed not just by general relativity and tides, but also by mass transfer. While the black hole and neutron star binaries observed with ground-based gravitational-wave detectors are driven to inspiral due to the emission of gravitational radiation—manifesting as a “chirp-like” gravitational-wave signal—the astrophysical processes at work in double white dwarf (DWD) systems can cause the inspiral to stall and even reverse into an outspiral. The dynamics of the DWD outspiral thus encode information about tides, which tell us about the behavior of electron-degenerate matter. We carry out a population study to determine the effect of the strength of tides on the distributions of the DWD binary parameters that the Laser Interferometer Space Antenna (LISA) will be able to constrain. We find that the strength of tidal coupling parameterized via the tidal synchronization timescale at the onset of mass transfer affects the distribution of gravitational-wave frequencies and frequency derivatives for detectably mass-transferring DWD systems. Using a hierarchical Bayesian framework informed by binary population synthesis simulations, we demonstrate how this parameter can be inferred using LISA observations. By measuring the population properties of DWDs, LISA will be able to probe the behavior of electron-degenerate matter.

     
    more » « less
  2. Abstract

    Future searches for gravitational waves from space will be sensitive to double compact objects in our Milky Way. We present new simulations of the populations of double black holes (BHBHs), BH neutron stars (BHNSs), and double neutron stars (NSNSs) that will be detectable by the planned space-based gravitational-wave detector called Laser Interferometer Space Antenna (LISA). For our estimates, we use an empirically informed model of the metallicity-dependent star formation history of the Milky Way. We populate it using an extensive suite of binary population-synthesis predictions for varying assumptions relating to mass transfer, common-envelope, supernova kicks, remnant masses, and wind mass-loss physics. For a 4(10) yr LISA mission, we predict between 30–370(50–550) detections over these variations, out of which 6–154 (9–238) are BHBHs, 2–198 (3–289) are BHNSs, and 3–35 (4–57) are NSNSs. We expect that about 50% (60%) can be distinguished from double white dwarf sources based on their mass or eccentricity and localization. Specifically, for about 10% (15%), we expect to be able to determine chirp masses better than 10%. For 13% (13%), we expect sky-localizations better than 1°. We discuss how the variations in the physics assumptions alter the distribution of properties of the detectable systems, even when the detection rates are unchanged. We further discuss the possibility of multimessenger observations of pulsar populations with the Square Kilometre Array and assess the benefits of extending the LISA mission.

     
    more » « less
  3. null (Ed.)
    ABSTRACT Compact white dwarf (WD) binaries are important sources for space-based gravitational-wave (GW) observatories, and an increasing number of them are being identified by surveys like Extremely Low Mass (ELM) and Zwicky Transient Facility (ZTF). We study the effects of non-linear dynamical tides in such binaries. We focus on the global three-mode parametric instability and show that it has a much lower threshold energy than the local wave-breaking condition studied previously. By integrating networks of coupled modes, we calculate the tidal dissipation rate as a function of orbital period. We construct phenomenological models that match these numerical results and use them to evaluate the spin and luminosity evolution of a WD binary. While in linear theory the WD’s spin frequency can lock to the orbital frequency, we find that such a lock cannot be maintained when non-linear effects are taken into account. Instead, as the orbit decays, the spin and orbit go in and out of synchronization. Each time they go out of synchronization, there is a brief but significant dip in the tidal heating rate. While most WDs in compact binaries should have luminosities that are similar to previous traveling-wave estimates, a few per cent should be about 10 times dimmer because they reside in heating rate dips. This offers a potential explanation for the low luminosity of the CO WD in J0651. Lastly, we consider the impact of tides on the GW signal and show that the Laser Interferometer Space Antenna (LISA) and TianGO can constrain the WD’s moment of inertia to better than $1{{\ \rm per\ cent}}$ for centi-Hz systems. 
    more » « less
  4. ABSTRACT

    Carbon-enhanced metal-poor (CEMP) stars comprise almost a third of stars with [Fe/H] < −2, although their origins are still poorly understood. It is highly likely that one sub-class (CEMP-s stars) is tied to mass-transfer events in binary stars, while another sub-class (CEMP-no stars) are enriched by the nucleosynthetic yields of the first generations of stars. Previous studies of CEMP stars have primarily concentrated on the Galactic halo, but more recently they have also been detected in the thick disc and bulge components of the Milky Way. Gaia DR3 has provided an unprecedented sample of over 200 million low-resolution (R ≈ 50) spectra from the BP and RP photometers. Training on the CEMP catalogue from the SDSS/SEGUE database, we use XGBoost to identify the largest all-sky sample of CEMP candidate stars to date. In total, we find 58 872 CEMP star candidates, with an estimated contamination rate of 12 per cent. When comparing to literature high-resolution catalogues, we positively identify 60–68 per cent of the CEMP stars in the data, validating our results and indicating a high completeness rate. Our final catalogue of CEMP candidates spans from the inner to outer Milky Way, with distances as close as r ∼ 0.8 kpc from the Galactic centre, and as far as r > 30 kpc. Future higher resolution spectroscopic follow-up of these candidates will provide validations of their classification and enable investigations of the frequency of CEMP-s and CEMP-no stars throughout the Galaxy, to further constrain the nature of their progenitors.

     
    more » « less
  5. ABSTRACT

    The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localization, multiwavelength, time-dependent EM models are needed to identify the right host galaxy. We studied merging MBHBs embedded in a circumbinary disc (CBD) using high-resolution two-dimensional simulations, with a Γ-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the CBD. We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multiband EM signature: we find that the MBHB produces strong thermal X-ray emission until 1–2 d prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days, regardless of any post-merger effects such as GW recoil or mass-loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring.

     
    more » « less