skip to main content

Title: Gravitational Wave Sources in Our Galactic Backyard: Predictions for BHBH, BHNS, and NSNS Binaries Detectable with LISA

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 more » 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.

« less
; ; ; ; ;
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Article No. 118
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this

    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 ofmore »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.

    « less

    A variety of observational campaigns seek to test dark matter models by measuring dark matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within Lambda cold dark matter, we examine 11 Milky Way (MW)-mass haloes from the FIRE-2 baryonic simulations, quantifying the counts and orbital fluxes for subhaloes with properties relevant to stellar stream interactions: masses down to $10^{6}\, \text{M}_\odot$, distances ≲50 kpc of the galactic centre, across z = 0 − 1 (tlookback = 0–8 Gyr). We provide fits to our results and their dependence on subhalo mass, distance, and lookback time, for use in (semi)analytical models. A typical MW-mass halo contains ≈16 subhaloes $\gt 10^{7}\, \text{M}_\odot$ (≈1 subhalo $\gt 10^{8}\, \text{M}_\odot$) within 50 kpc at z ≈ 0. We compare our results with dark matter-only versions of the same simulations: because they lack a central galaxy potential, they overpredict subhalo counts by 2–10×, more so at smaller distances. Subhalo counts around a given MW-mass galaxy declined over time, being ≈10× higher at z = 1 than at z ≈ 0. Subhaloes have nearly isotropic orbital velocity distributions at z ≈more »0. Across our simulations, we also identified 4 analogues of Large Magellanic Cloud satellite passages; these analogues enhance subhalo counts by 1.4–2.1 times, significantly increasing the expected subhalo population around the MW today. Our results imply an interaction rate of ∼5 per Gyr for a stream like GD-1, sufficient to make subhalo–stream interactions a promising method of measuring dark subhaloes.

    « less
  3. Abstract Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly 50 M ⊙ and 100 M ⊙ , while above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational-wave emission with the current LIGO/Virgo/Kagra interferometers and future detectors such as LISA or the Einstein Telescope, can form both mass-gap BHs and IMBHs. Here we explore the possibility that mass-gap BHs and IMBHs are born as a result of successive BH mergers in dense star clusters. In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products after they receive significant recoil kicks due to anisotropic emission of gravitational radiation. Using for the first time simulations that include full stellar evolution, we show that a massive stellar BH seed can easily grow to ∼10 3 –10 4 M ⊙ as a result of repeated mergers with other smaller BHs. We find that lowering the cluster metallicity leads to larger final BH masses. We also show that the growing BH spin tends to decrease in magnitude with the number ofmore »mergers so that a negative correlation exists between the final mass and spin of the resulting IMBHs. Assumptions about the birth spins of stellar BHs affect our results significantly, with low birth spins leading to the production of a larger population of massive BHs.« less

    Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work, we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, Z, and redshift z, $\mathcal {S}(Z,z)$. Considering these uncertainties, we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal {S}(Z,z)$ variations can combined impact the predicted intrinsic and detectable merger rates by factors in the range 102–104. We find that BHBH rates are dominantly impacted by $\mathcal {S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses, and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS, we find that both uncertainties are important. Wemore »also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal {S}(Z,z)$ for BHBH, BHNS, and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95 per cent of BHBH detections to contain a BH $\gtrsim 8\, \rm {M}_{\odot }$ and have mass ratios ≲ 4. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties. Conversely, larger observed samples could allow us to decipher currently unconstrained stages of stellar and binary evolution.

    « less
  5. ABSTRACT The progenitor system of the compact binary merger GW190425 had a total mass of $3.4^{+0.3}_{-0.1}$ M⊙ (90th-percentile confidence region) as measured from its gravitational wave signal. This mass is significantly different from the Milky Way (MW) population of binary neutron stars (BNSs) that are expected to merge in a Hubble time and from that of the first BNS merger, GW170817. Here, we explore the expected electromagnetic (EM) signatures of such a system. We make several astrophysically motivated assumptions to further constrain the parameters of GW190425. By simply assuming that both components were NSs, we reduce the possible component masses significantly, finding $m_{1}=1.85^{+0.27}_{-0.19}$ M⊙ and $m_{2}=1.47^{+0.16}_{-0.18}$ M⊙. However, if the GW190425 progenitor system was an NS–black hole (BH) merger, we find best-fitting parameters $m_{1}=2.19^{+0.21}_{-0.17}$ M⊙ and $m_{2}=1.26^{+0.10}_{-0.08}$ M⊙. For a well-motivated BNS system where the lighter NS has a mass similar to the mass of non-recycled NSs in MW BNS systems, we find $m_{1}=2.03^{+0.15}_{-0.14}$ M⊙ and m2 = 1.35 ± 0.09 M⊙, corresponding to only 7 per cent mass uncertainties. For all scenarios, we expect a prompt collapse of the resulting remnant to a BH. Examining detailed models with component masses similar to our best-fitting results, we find the EM counterpart to GW190425 is expected to be significantly redder and fainter thanmore »that of GW170817. We find that almost all reported search observations were too shallow to detect the expected counterpart to GW190425. If the LIGO–Virgo Collaboration promptly provides the chirp mass, the astronomical community can adapt their observations to improve the likelihood of detecting a counterpart for similarly ‘high-mass’ BNS systems.« less