The construction of a model of the gravitational-wave (GW) signal from generic configurations of spinning-black-hole binaries, through inspiral, merger, and ringdown, is one of the most pressing theoretical problems in the buildup to the era of GW astronomy. We present the first such model in the frequency domain, PhenomP, which captures the basic phenomenology of the seven-dimensional parameter space of binary configurations with only three key physical parameters. Two of these (the binary’s mass ratio and an effective total spin parallel to the orbital angular momentum, which determines the inspiral rate) define an underlying nonprecessing-binary model. The nonprecessing-binary waveforms are then twisted up with approximate expressions for the precessional motion, which require only one additional physical parameter, an effective precession spin, ${\chi }_p$. All other parameters (total mass, sky location, orientation and polarization, and initial phase) can be specified trivially. The model is constructed in the frequency domain, which will be essential for efficient GW searches and source measurements. We have tested the model’s fidelity for GW applications by comparison against hybrid post-Newtonian-numerical-relativity waveforms at a variety of configurations—although we did not use these numerical simulations in the construction of the model. Our model can be used to develop GW searches, to study the implications for astrophysical measurements, and as a simple conceptual framework to form the basis of generic-binary waveform modeling in the advanced-detector era.

We present imrphenomxphm, a phenomenological frequency-domain model for the gravitational-wave signal emitted by quasicircular precessing binary black holes, which incorporates multipoles beyond the dominant quadrupole in the precessing frame. The model is a precessing extension of imrphenomxhm, [C. García-Quirós et al., Phys. Rev. D 102, 064002 (2020)] based on approximate maps between aligned-spin waveform modes in the coprecessing frame and precessing waveform modes in the inertial frame, which is commonly referred to as “twisting up” the nonprecessing waveforms. imrphenomxphm includes imrphenomxp as a special case, the restriction to the dominant quadrupole contribution in the coprecessing frame. We implement two alternative mappings, one based on a single-spin post-Newtonian approximation, as used in imrphenompv2 [M. Hannam et al., Phys. Rev. Lett. 113, 151101 (2014).], and one based on the double-spin multiple scale analysis approach of [K. Chatziioannou et al., Phys. Rev. D 95, 104004 (2017).]. We include a detailed discussion of conventions used in the description of precessing binaries and of all choices made in constructing the model. The computational cost of imrphenomxphm is further reduced by extending the interpolation technique of [C. García-Quirós et al., Classical Quant. Grav. 38, 015006 (2021).] to the Euler angles. The accuracy, speed, robustness, and modularity of the imrphenomx family will make these models productive tools for gravitational wave astronomy in the current era of greatly increased number and diversity of detected events.

The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in 2019 April and lasting six months, O3b starting in 2019 November and lasting five months, and O3GK starting in 2020 April and lasting two weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org. The main data set, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages.

In this paper we present imrphenomxas, a thorough overhaul of the imrphenomd [S. Husa et al., Phys. Rev. D 93, 044006 (2016); S. Khan et al., Phys. Rev. D 93, 044007 (2016)] waveform model, which describes the dominant $l=2$, $|m|=2$ spherical harmonic mode of nonprecessing coalescing black holes in terms of piecewise closed form expressions in the frequency domain. Improvements include in particular the accurate treatment of unequal spin effects, and the inclusion of extreme mass ratio waveforms. imrphenomd has previously been extended to approximately include spin precession [M. Hannam et al., Phys. Rev. Lett. 113, 151101 (2014)] and subdominant spherical harmonics [L. London et al., Phys. Rev. Lett. 120, 161102 (2018)], and with its extensions it has become a standard tool in gravitational wave parameter estimation. Improved extensions of imrphenomxas are discussed in companion papers [C. García-Quirós et al., Phys. Rev. D 102, 064002 (2020); G. Pratten et al., arXiv:2004.06503].

We present TEOBResumS, a new effective-one-body (EOB) waveform model for nonprecessing (spin-aligned) and tidally interacting compact binaries. Spin-orbit and spin-spin effects are blended together by making use of the concept of centrifugal EOB radius. The point-mass sector through merger and ringdown is informed by numerical relativity (NR) simulations of binary black holes (BBHs) computed with the SpEC and bam codes. An improved, NR-based phenomenological description of the postmerger waveform is developed. The tidal sector of TEOBResumS describes the dynamics of neutron star binaries up to merger and incorporates a resummed attractive potential motivated by recent advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. Equation-of-state-dependent self-spin interactions (monopole-quadrupole effects) are incorporated in the model using leading order post-Newtonian results in a new expression of the centrifugal radius. TEOBResumS is compared to 135 SpEC and 19 bam BBH waveforms. The maximum unfaithfulness to SpEC data $\overline{F}$—at design Advanced LIGO sensitivity and evaluated with total mass $M$ with a variance of $10M_{\odot }\le M\le 200M_{\odot }$—is always below $2.5\times {10}^{-3}$ except for a single outlier that grazes the $7.1\times {10}^{-3}$ level. When compared to bam data, $\overline{F}$ is smaller than 0.01 except for a single outlier in one of the corners of the NR-covered parameter space that reaches the 0.052 level. TEOBResumS is also compatible, up to merger, to high-end NR waveforms from binary neutron stars with spin effects and reduced initial eccentricity computed with the bam and thc codes. The data quality of binary neutron star waveforms is assessed via rigorous convergence tests from multiple resolution runs and takes into account systematic effects estimated by using the two independent high-order NR codes. The model is designed to generate accurate templates for the analysis of LIGO-Virgo data through merger and ringdown. We demonstrate its use by analyzing the publicly available data for GW150914.