skip to main content

Title: LEGWORK: A Python Package for Computing the Evolution and Detectability of Stellar-origin Gravitational-wave Sources with Space-based Detectors

We present LEGWORK (LISA Evolution and Gravitational Wave Orbit Kit), an open-source Python package for making predictions about stellar-origin gravitational-wave sources and their detectability in LISA or other space-based gravitational-wave detectors. LEGWORK can be used to evolve the orbits of sources due to gravitational-wave emission, calculate gravitational-wave strains (using post-Newtonian approximations), compute signal-to-noise ratios, and visualize the results. It can be applied to a variety of potential sources, including binaries consisting of white dwarfs, neutron stars, and black holes. Although we focus on double compact objects, in principle LEGWORK can be used for any system with a user-specified orbital evolution, such as those affected by a third object or gas drag. We optimized the package to make it efficient for use in population studies, which can contain tens of millions of sources. This paper describes the package and presents several potential use cases. We explain in detail the derivations of the expressions behind the package as well as identify and clarify some discrepancies currently present in the literature. We hope that LEGWORK will enable and accelerate future studies triggered by the rapidly growing interest in gravitational-wave sources.

more » « less
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Supplement Series
Page Range / eLocation ID:
Article No. 52
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    With the inception of gravitational wave astronomy, astrophysical studies using interferometric techniques have begun to probe previously unknown parts of the Universe. In this work, we investigate the potential of a new interferometric experiment to study a unique group of gravitationally interacting sources within our Solar system: binary asteroids. We present the first study into binary asteroid detection via gravitational signals. We identify the interferometer sensitivity necessary for detecting a population of binary asteroids in the asteroid belt. We find that the space-based gravitational wave detector LISA will have negligible ability to detect these sources as these signals will be well below the LISA noise curve. Consequently, we propose a 4.6 au and a 1 au arm-length interferometer specialized for binary asteroid detection, targeting frequencies between 10−6 and 10−4 Hz. Our results demonstrate that the detection of binary asteroids with space-based gravitational wave interferometers is possible though very difficult, requiring substantially improved interferometric technology over what is presently proposed for space-based missions. If that threshold can be met, an interferometer may be used to map the asteroid belt, allowing for new studies into the evolution of our Solar system.

    more » « less

    The upcoming Laser Interferometer Space Antenna (LISA) will detect a large gravitational-wave foreground of Galactic white dwarf binaries. These sources are exceptional for their probable detection at electromagnetic wavelengths, some long before LISA flies. Studies in both gravitational and electromagnetic waves will yield strong constraints on system parameters not achievable through measurements of one messenger alone. In this work, we present a Bayesian inference pipeline and simulation suite in which we study potential constraints on binaries in a variety of configurations. We show how using LISA detections and parameter estimation can significantly improve constraints on system parameters when used as a prior for the electromagnetic analyses. We also provide rules of thumb for how current measurements will benefit from LISA measurements in the future.

    more » « less
  3. 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
  4. 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

    The launch of space-based gravitational-wave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays will extend the GW window to low frequencies, opening new investigations into dynamical processes involving massive black hole binaries (MBHBs) and their mergers across cosmic time. MBHBs are expected to be among the primary sources for the upcoming low-frequency (10−4–10−1 Hz) window probed by LISA. It is important to investigate the expected supermassive BH merger rates and associated signals, to determine how potential LISA events are affected by physics included in current models. To study this, we post-process the large population of MBHBs in the Illustris simulation to account for dynamical friction time delays associated with BH infall/inspiral. We show that merger delays associated with binary evolution have the potential to decrease the expected merger rates, with $M_{\rm {BH}}\ \gt\ 10^6\ \mathrm{M}_\odot$ MBHBs (the lowest mass in Illustris) decreasing from ∼3 to ∼0.1 yr−1, and shifting the merger peak from z ∼2 to ∼1.25. During this time, we estimate that accretion grows the total merging mass by as much as 7x the original mass. Importantly, however, dynamical friction-associated delays (which shift the mergers toward lower redshift and higher masses) lead to a stronger signal/strain for the emitted GWs in the LISA band, increasing mean frequency from 10−3.1 to 10−3.4–10−4.0 Hz, and mean strain from 10−17.2 to 10−16.3–10−15.3. Finally, we show that after including a merger delay and associated MBH growth, mergers still tend to lie on the typical MBH–M* relation, but with an increased likelihood of an undermassive BH.

    more » « less