skip to main content

Title: Bayesian inference for compact binary coalescences with BILBY: validation and application to the first LIGO-Virgo gravitational-wave transient catalogue
Gravitational waves provide a unique tool for observational astronomy. While the first LIGO–Virgo catalogue of gravitational wave transients (GWTC-1) contains 11 signals from black hole and neutron star binaries, the number of observations is increasing rapidly as detector sensitivity improves. To extract information from the observed signals, it is imperative to have fast, flexible, and scalable inference techniques. In a previous paper, we introduced BILBY: a modular and user-friendly Bayesian inference library adapted to address the needs of gravitational-wave inference. In this work, we demonstrate that BILBY produces reliable results for simulated gravitational-wave signals from compact binary mergers, and verify that it accurately reproduces results reported for the 11 GWTC-1 signals. Additionally, we provide configuration and output files for all analyses to allow for easy reproduction, modification, and future use. This work establishes that BILBY is primed and ready to analyse the rapidly growing population of compact binary coalescence gravitational-wave signals.  more » « less
Award ID(s):
1912594 1807046 1806630
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
Monthly notices of the Royal Astronomical Society
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The observation of gravitational waves from multiple compact binary coalescences by the LIGO–Virgo–KAGRA detector networks has enabled us to infer the underlying distribution of compact binaries across a wide range of masses, spins, and redshifts. In light of the new features found in the mass spectrum of binary black holes and the uncertainty regarding binary formation models, nonparametric population inference has become increasingly popular. In this work, we develop a data-driven clustering framework that can identify features in the component mass distribution of compact binaries simultaneously with those in the corresponding redshift distribution, from gravitational-wave data in the presence of significant measurement uncertainties, while making very few assumptions about the functional form of these distributions. Our generalized model is capable of inferring correlations among various population properties, such as the redshift evolution of the shape of the mass distribution itself, in contrast to most existing nonparametric inference schemes. We test our model on simulated data and demonstrate the accuracy with which it can reconstruct the underlying distributions of component masses and redshifts. We also reanalyze public LIGO–Virgo–KAGRA data from events in GWTC-3 using our model and compare our results with those from some alternative parametric and nonparametric population inference approaches. Finally, we investigate the potential presence of correlations between mass and redshift in the population of binary black holes in GWTC-3 (those observed by the LIGO–Virgo–KAGRA detector network in their first three observing runs), without making any assumptions about the specific nature of these correlations.

    more » « less
  2. Standard sirens have been the central paradigm in gravitational-wave cosmology so far. From the gravitational wave signature of compact star binaries, it is possible to measure the luminosity distance of the source directly, and if additional information on the source redshift is provided, a measurement of the cosmological expansion can be performed. This review article discusses several methodologies that have been proposed to use gravitational waves for cosmological studies. Methods that use only gravitational-wave signals and methods that use gravitational waves in conjunction with additional observations such as electromagnetic counterparts and galaxy catalogs will be discussed. The review also discusses the most recent results on gravitational-wave cosmology, starting from the binary neutron star merger GW170817 and its electromagnetic counterpart and finishing with the population of binary black holes, observed with the third Gravitational-wave Transient Catalog GWTC–3. 
    more » « less
  3. Abstract We present our current best estimate of the plausible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next several years, with the intention of providing information to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals for the third (O3), fourth (O4) and fifth observing (O5) runs, including the planned upgrades of the Advanced LIGO and Advanced Virgo detectors. We study the capability of the network to determine the sky location of the source for gravitational-wave signals from the inspiral of binary systems of compact objects, that is binary neutron star, neutron star–black hole, and binary black hole systems. The ability to localize the sources is given as a sky-area probability, luminosity distance, and comoving volume. The median sky localization area (90% credible region) is expected to be a few hundreds of square degrees for all types of binary systems during O3 with the Advanced LIGO and Virgo (HLV) network. The median sky localization area will improve to a few tens of square degrees during O4 with the Advanced LIGO, Virgo, and KAGRA (HLVK) network. During O3, the median localization volume (90% credible region) is expected to be on the order of $$10^{5}, 10^{6}, 10^{7}\mathrm {\ Mpc}^3$$ 10 5 , 10 6 , 10 7 Mpc 3 for binary neutron star, neutron star–black hole, and binary black hole systems, respectively. The localization volume in O4 is expected to be about a factor two smaller than in O3. We predict a detection count of $$1^{+12}_{-1}$$ 1 - 1 + 12 ( $$10^{+52}_{-10}$$ 10 - 10 + 52 ) for binary neutron star mergers, of $$0^{+19}_{-0}$$ 0 - 0 + 19 ( $$1^{+91}_{-1}$$ 1 - 1 + 91 ) for neutron star–black hole mergers, and $$17^{+22}_{-11}$$ 17 - 11 + 22 ( $$79^{+89}_{-44}$$ 79 - 44 + 89 ) for binary black hole mergers in a one-calendar-year observing run of the HLV network during O3 (HLVK network during O4). We evaluate sensitivity and localization expectations for unmodeled signal searches, including the search for intermediate mass black hole binary mergers. 
    more » « less

    The detection of an intermediate-mass black hole population (102–106 M⊙) will provide clues to their formation environments (e.g. discs of active galactic nuclei, globular clusters) and illuminate a potential pathway to produce supermassive black holes. Ground-based gravitational-wave detectors are sensitive to mergers that can form intermediate-mass black holes weighing up to ∼450 M⊙. However, ground-based detector data contain numerous incoherent short duration noise transients that can mimic the gravitational-wave signals from merging intermediate-mass black holes, limiting the sensitivity of searches. Here, we follow-up on binary black hole merger candidates using a ranking statistic that measures the coherence or incoherence of triggers in multiple-detector data. We use this statistic to rank candidate events, initially identified by all-sky search pipelines, with lab-frame total masses ≳ 55 M⊙ using data from LIGO’s second observing run. Our analysis does not yield evidence for new intermediate-mass black holes. However, we find support for eight stellar-mass binary black holes not reported in the first LIGO–Virgo gravitational wave transient catalogue GWTC-1, seven of which have been previously reported by other catalogues.

    more » « less

    We consider the 10 confidently detected gravitational-wave signals in the GWTC-1 catalog, which are consistent with mergers of binary black hole systems, and perform a thorough parameter estimation re-analysis. This is made possible by using computationally efficient waveform models of the current (fourth) generation of the IMRPhenom family of phenomenological waveform models, which consists of the IMRPhenomX frequency-domain models and the IMRPhenomT time-domain models.The analysis is performed with both precessing and non-precessing waveform models with and without subdominant spherical harmonic modes. Results for all events are validated with convergence tests, discussing in particular the events GW170729 and GW151226. For the latter and the other two lowest-mass events, we also compare results between two independent sampling codes, bilbyand lalinference. We find overall consistent results with the original GWTC-1 results, with all Jensen–Shannon divergences between the previous results using IMRPhenomPv2 and our default IMRPhenomXPHM posteriors below 0.045 bits. However, we also discuss cases where subdominant harmonics and/or precession influence the posteriors.

    more » « less