skip to main content

Title: Posterior samples of the parameters of binary black holes from Advanced LIGO, Virgo’s second observing run

This paper presents a parameter estimation analysis of the seven binary black hole mergers—GW170104, GW170608, GW170729, GW170809, GW170814, GW170818, and GW170823—detected during the second observing run of the Advanced LIGO and Virgo observatories using the gravitational-wave open data. We describe the methodology for parameter estimation of compact binaries using gravitational-wave data, and we present the posterior distributions of the inferred astrophysical parameters. We release our samples of the posterior probability density function with tutorials on using and replicating our results presented in this paper.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Data
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Since the very first detection of gravitational waves from the coalescence of two black holes in 2015, Bayesian statistical methods have been routinely applied by LIGO and Virgo to extract the signal out of noisy interferometric measurements, obtain point estimates of the physical parameters responsible for producing the signal, and rigorously quantify their uncertainties. Different computational techniques have been devised depending on the source of the gravitational radiation and the gravitational waveform model used. Prominent sources of gravitational waves are binary black hole or neutron star mergers, the only objects that have been observed by detectors to date. But also gravitational waves from core‐collapse supernovae, rapidly rotating neutron stars, and the stochastic gravitational‐wave background are in the sensitivity band of the ground‐based interferometers and expected to be observable in future observation runs. As nonlinearities of the complex waveforms and the high‐dimensional parameter spaces preclude analytic evaluation of the posterior distribution, posterior inference for all these sources relies on computer‐intensive simulation techniques such as Markov chain Monte Carlo methods. A review of state‐of‐the‐art Bayesian statistical parameter estimation methods will be given for researchers in this cross‐disciplinary area of gravitational wave data analysis.

    This article is categorized under:

    Applications of Computational Statistics > Signal and Image Processing and Coding

    Statistical and Graphical Methods of Data Analysis > Markov Chain Monte Carlo (MCMC)

    Statistical Models > Time Series Models

    more » « less
  2. Abstract We introduce deep learning models to estimate the masses of the binary components of black hole mergers, ( m 1 , m 2 ) , and three astrophysical properties of the post-merger compact remnant, namely, the final spin, a f , and the frequency and damping time of the ringdown oscillations of the fundamental ℓ = m = 2 bar mode, ( ω R , ω I ) . Our neural networks combine a modified WaveNet architecture with contrastive learning and normalizing flow. We validate these models against a Gaussian conjugate prior family whose posterior distribution is described by a closed analytical expression. Upon confirming that our models produce statistically consistent results, we used them to estimate the astrophysical parameters ( m 1 , m 2 , a f , ω R , ω I ) of five binary black holes: GW150914 , GW170104 , GW170814 , GW190521 and GW190630 . We use PyCBC Inference to directly compare traditional Bayesian methodologies for parameter estimation with our deep learning based posterior distributions. Our results show that our neural network models predict posterior distributions that encode physical correlations, and that our data-driven median results and 90% confidence intervals are similar to those produced with gravitational wave Bayesian analyses. This methodology requires a single V100 NVIDIA GPU to produce median values and posterior distributions within two milliseconds for each event. This neural network, and a tutorial for its use, are available at the Data and Learning Hub for Science . 
    more » « less

    On 2017 August 17, Advanced LIGO and Virgo observed GW170817, the first gravitational-wave (GW) signal from a binary neutron star merger. It was followed by a short-duration gamma-ray burst, GRB 170817A, and by a non-thermal afterglow emission. In this work, a combined simultaneous fit of the electromagnetic (EM, specifically, afterglow) and GW domains is implemented, both using the posterior distribution of a GW standalone analysis as prior distribution to separately process the EM data, and fitting the EM and GW domains simultaneously. These approaches coincide mathematically, as long as the actual posterior of the GW analysis, and not an approximation, is used as prior for the EM analysis. We treat the viewing angle, θv, as shared parameter across the two domains. In the afterglow modelling with a Gaussian structured jet this parameter and the jet core angle, θc, are correlated, leading to high uncertainties on their values. The joint EM + GW analysis relaxes this degeneracy, reducing the uncertainty compared to an EM-only fit. We also apply our methodology to hypothetical GW170817-like events occurring in the next GW observing run at ∼140 and 70 Mpc. At 70 Mpc the existing EM degeneracy is broken, thanks to the inclusion of the GW domain in the analysis. At 140 Mpc, the EM-only fit cannot constrain θv nor θc because of the lack of detections in the afterglow rising phase. Folding the GW data into the analysis leads to tighter constraints on θv, still leaving θc unconstrained, requiring instruments with higher sensitivities, such as Athena.

    more » « less
  4. Abstract

    The detection of GW170817, the first neutron star-neutron star merger observed by Advanced LIGO and Virgo, and its following analyses represent the first contributions of gravitational wave data to understanding dense matter. Parameterizing the high density section of the equation of state of both neutron stars through spectral decomposition, and imposing a lower limit on the maximum mass value, led to an estimate of the stars’ radii ofkm andkm (Abbottet al2018Phys. Rev. Lett.121161101). These values do not, however, take into account any uncertainty owed to the choice of the crust low-density equation of state, which was fixed to reproduce the SLy equation of state model (Douchin and Haensel 2001Astron. Astrophys.380151). We here re-analyze GW170817 data and establish that different crust models do not strongly impact the mass or tidal deformability of a neutron star—it is impossible to distinguish between low-density models with gravitational wave analysis. However, the crust does have an effect on inferred radius. We predict the systematic error due to this effect using neutron star structure equations, and compare the prediction to results from full parameter estimation runs. For GW170817, this systematic error affects the radius estimate by 0.3 km, approximatelyof the neutron stars’ radii.

    more » « less
  5. Abstract

    Integrated population models (IPMs) have become increasingly popular for the modelling of populations, as investigators seek to combine survey and demographic data to understand processes governing population dynamics. These models are particularly useful for identifying and exploring knowledge gaps within life histories, because they allow investigators to estimate biologically meaningful parameters, such as immigration or reproduction, that were previously unidentifiable without additional data. AsIPMs have been developed relatively recently, there is much to learn about model behaviour. Behaviour of parameters, such as estimates near boundaries, and the consequences of varying degrees of dependency among datasets, has been explored. However, the reliability of parameter estimates remains underexamined, particularly when models include parameters that are not identifiable from one data source, but are indirectly identifiable from multiple datasets and a presumed model structure, such as the estimation of immigration using capture‐recapture, fecundity and count data, combined with a life‐history model.

    To examine the behaviour of model parameter estimates, we simulated stable populations closed to immigration and emigration. We simulated two scenarios that might induce error into survival estimates: marker induced bias in the capture–mark–recapture data and heterogeneity in the mortality process. We subsequently fit capture–mark–recapture, state‐space and fecundity models, as well asIPMs that estimated additional parameters.

    Simulation results suggested that when model assumptions are violated, estimation of additional, previously unidentifiable, parameters usingIPMs may be extremely sensitive to these violations of model assumption. For example, when annual marker loss was simulated, estimates of survival rates were low and estimates of immigration rate from anIPMwere high. When heterogeneity in the mortality process was induced, there were substantial relative differences between the medians of posterior distributions and truth for juvenile survival and fecundity.

    Our results have important implications for biological inference when usingIPMs, as well as future model development and implementation. Specifically, using multiple datasets to identify additional parameters resulted in the posterior distributions of additional parameters directly reflecting the effects of the violations of model assumptions in integrated modelling frameworks. We suggest that investigators interpret posterior distributions of these parameters as a combination of biological process and systematic error.

    more » « less