Gravitational waves provide a unique tool for observational astronomy. While the first LIGO--Virgo catalogue of gravitational-wave transients (GWTC-1) contains eleven 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 eleven 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.

Machine learning has emerged as a popular and powerful approach for solving problems in astrophysics. We review applications of machine learning techniques for the analysis of ground-based gravitational-wave detector data. Examples include techniques for improving the sensitivity of Advanced LIGO and Advanced Virgo gravitational-wave searches, methods for fast measurements of the astrophysical parameters of gravitational-wave sources, and algorithms for reduction and characterization of non-astrophysical detector noise. These applications demonstrate how machine learning techniques may be harnessed to enhance the science that is possible with current and future gravitational-wave detectors.

We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

The combined observation of gravitational and electromagnetic waves from the coalescence of two neutron stars marks the beginning of multi-messenger astronomy with gravitational waves (GWs). The development of accurate gravitational waveform models is a crucial prerequisite to extract information about the properties of the binary system that generated a detected GW signal. In binary neutron star systems (BNS), tidal effects also need to be incorporated in the modeling for an accurate waveform representation. Building on previous work [Phys.Rev.D96 121501], we explore the performance of inspiral-merger waveform models that are obtained by adding a numerical relativity (NR) based approximant for the tidal part of the phasing (NRTidal) to existing models for nonprecessing and precessing binary black hole systems (SEOBNRv4, PhenomD and PhenomPv2), as implemented in the LSC Algorithm Library Suite. The resulting BNS waveforms are compared and contrasted to target waveforms hybridizing NR waveforms, covering the last approx. 10 orbits up to merger and extending through the postmerger phase, with inspiral waveforms calculated from 30Hz obtained with TEOBResumS. The latter is a state-of-the-art effective-one-body waveform model that blends together tidal and spin effects. We probe that the combination of the PN-based self-spin terms and of the NRTidal description is necessary to obtain minimal mismatches (< 0.01) and phase differences (< 1 rad) with respect to the target waveforms. However, we also discuss possible improvements and drawbacks of the NRTidal approximant in its current form, since we find that it tends to overestimate the tidal interaction with respect to the TEOBResumS model during the inspiral.