skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on February 26, 2026

Title: Recovering Injected Astrophysics from the LISA Double White Dwarf Binaries
Abstract We present the successful recovery of common-envelope ejection efficiency assumed in a simulated population of double white dwarf (DWD) binaries like those which may be observed by the future Laser Interferometer Space Antenna (LISA) mission. We simulate the formation of DWD binaries by using the COSMIC population synthesis code to sample binary formation conditions such as initial mass function, metallicity of star formation, initial orbital period, and initial eccentricity. These binaries are placed in the m12i synthetic Milky Way–like galaxy, and their signal-to-noise ratio (SNR) for the LISA instrument is estimated, considering a Galactic gravitational-wave foreground informed by the population. Through the use of Fisher estimates, we construct a likelihood function for the measurement error of the LISA-bright DWD binaries (≥20 SNR,fGW≥ 5 mHz), in their gravitational-wave frequency (fGW) and chirp mass. By repeating this process for different assumptions of the common-envelope ejection efficiency, we apply Bayesian hierarchical inference to find the best match to an injected astrophysical assumption for a fiducial population model. We conclude that the impact of common-envelope ejection efficiency on the mass-transfer processes involved in DWD formation may be statistically relevant in the future observed LISA population, and that constraints on binary formation may be found by comparing simulated populations to a future observed population.  more » « less
Award ID(s):
2012057 2206321 2309172
PAR ID:
10615793
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Astronomical Society
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
981
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
66
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The coalescence of two neutron stars was recently observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12 M⊙ red supergiant and a 1.4 M⊙ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and a final orbital separation of ≈1.1 - 2.8R⊙ , leading to a binary neutron star that will merge within 0.01-1 Gyr. We find an αCE -equivalent efficiency of ≈0.1 - 0.4 for the models we study, but this may be specific for these extended progenitors. We fully resolve the core of the star to ≲0.005R⊙ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources. 
    more » « less
  2. Abstract The dynamical formation of binary black holes (BBHs) in globular clusters (GCs) may contribute significantly to the observed gravitational-wave (GW) merger rate. Furthermore, the Laser Interferometer Space Antenna (LISA) may detect many BBH sources from GCs at mHz frequencies, enabling the characterization of such systems within the Milky Way and nearby Universe. In this work, we use Monte CarloN-body simulations to construct a realistic sample of Galactic clusters, thus estimating the population, detectability, and parameter measurement accuracy of BBHs hosted within them. In particular, we show that the GW signal from 0.7 ± 0.7, 2.0 ± 1.7, 3.6 ± 2.3, and 13.4 ± 4.7 BBHs in Milky Way GCs can exceed the signal-to-noise ratio (SNR) threshold of SNR = 30, 5, 3, and 1 for a 10 yr LISA observation, with ∼50% of detectable sources exhibiting high eccentricities (e ≳ 0.9). Moreover, the Fisher matrix and Bayesian analyses of the GW signals indicate that these systems typically feature highly resolved orbital frequencies (δforb/forb ∼ 10−7to 10−5) and eccentricities (δe/e ∼ 10−3to 0.1), as well as a measurable total mass when SNR exceeds ∼20. Notably, we show that high-SNR BBHs can be confidently localized to specific Milky Way GCs with a sky localization accuracy ofδΩ ∼ 1 deg2, and we address the large uncertainties in their distance measurement (δR ∼ 0.3–20 kpc). The detection and localization of even a single BBH in a Galactic GC would allow accurate tracking of its long-term orbital evolution, enable a direct test of the role of GCs in BBH formation, and provide a unique probe into the evolutionary history of Galactic clusters. 
    more » « less
  3. Abstract The dynamical formation channels of gravitational wave (GW) sources typically involve a stage when the compact object binary source interacts with the environment, which may excite its eccentricity, yielding efficient GW emission. For the wide eccentric compact object binaries, the GW emission happens mostly near the pericenter passage, creating a unique, burst-like signature in the waveform. This work examines the possibility of stellar-mass bursting sources in the mHz band for future LISA detections. Because of their long lifetime (∼107yr) and promising detectability, the number of mHz bursting sources can be large in the local Universe. For example, based on our estimates, there will be ∼3–45 bursting binary black holes in the Milky Way, with ∼102–104bursts detected during the LISA mission. Moreover, we find that the number of bursting sources strongly depends on their formation history. If certain regions undergo active formation of compact object binaries in the recent few million years, there will be a significantly higher bursting source fraction. Thus, the detection of mHz GW bursts not only serves as a clue for distinguishing different formation channels, but also helps us understand the star formation history in different regions of the Milky Way. 
    more » « less
  4. Abstract Gravitational-wave (GW) detectors are observing compact object mergers from increasingly far distances, revealing the redshift evolution of the binary black hole (BBH)—and soon the black hole–neutron star (BHNS) and binary neutron star (BNS)—merger rate. To help interpret these observations, we investigate the expected redshift evolution of the compact object merger rate from the isolated binary evolution channel. We present a publicly available catalog of compact object mergers and their accompanying cosmological merger rates from population synthesis simulations conducted with the COMPAS software. To explore the impact of uncertainties in stellar and binary evolution, our simulations use two-parameter grids of binary evolution models that vary the common-envelope efficiency with mass transfer accretion efficiency and supernova (SN) remnant mass prescription with SN natal kick velocity, respectively. We quantify the redshift evolution of our simulated merger rates using the local (z∼ 0) rate, the redshift at which the merger rate peaks, and the normalized differential rates (as a proxy for slope). We find that although the local rates span a range of ∼103across our model variations, their redshift evolutions are remarkably similar for BBHs, BHNSs, and BNSs, with differentials typically within a factor 3 and peaks ofz≈ 1.2–2.4 across models. Furthermore, several trends in our simulated rates are correlated with the model parameters we explore. We conclude that future observations of the redshift evolution of the compact object merger rate can help constrain binary models for stellar evolution and GW formation channels. 
    more » « less
  5. ABSTRACT The Laser Interferometer Space Antenna (LISA) is expected to detect a wide variety of gravitational wave sources in the mHz band. Some of these signals will elude individual detection, instead contributing as confusion noise to one of several stochastic gravitational-wave backgrounds (SGWBs) – notably including the ‘Galactic foreground’, a loud signal resulting from the superposition of millions of unresolved double white dwarf binaries (DWDs) in the Milky Way. It is possible that similar, weaker SGWBs will be detectable from other DWD populations in the local Universe, including the Large Magellanic Cloud (LMC). We use the Bayesian LISA Inference Package (blip) to investigate the possibility of an anisotropic SGWB generated by unresolved DWDs in the LMC. To do so, we compute the LMC SGWB from a realistic DWD population generated via binary population synthesis, simulate 4 years of time-domain data with blip comprised of stochastic contributions from the LMC SGWB and the LISA detector noise, and analyse this data with blip’s spherical harmonic anisotropic SGWB search. We also consider the case of spectral separation from the Galactic foreground. We present the results of these analyses and show, for the first time, that the unresolved DWDs in the LMC will comprise a significant SGWB for LISA. 
    more » « less