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.

For the first time, we construct an inspiral-merger-ringdown waveform model within the effective-one-body formalism for spinning, nonprecessing binary black holes that includes gravitational modes beyond the dominant (l vertical bar m vertical bar) = (2, 2) mode, specifically (l vertical bar m vertical bar) = (2, 1), (3, 3), (4, 4), (5, 5). Our multipolar waveform model incorporates recent (resummed) post-Newtonian results for the inspiral and information from 157 numerical-relativity simulations, and 13 waveforms from black-hole perturbation theory for the (plunge-) merger and ringdown. We quantify the improvement in accuracy when including higher-order modes by computing the faithfulness of the waveform model against the numerical-relativity waveforms used to construct the model. We define the faithfulness as the match maximized over time, phase of arrival, gravitational-wave polarization and sky position of the waveform model, and averaged over binary orientation, gravitational-wave polarization and sky position of the numerical-relativity waveform. When the waveform model contains only the (2,2) mode, we find that the averaged faithfulness to numerical-relativity waveforms containing all modes with l <= 5 ranges from 90% to 99.9% for binaries with total mass 20-200 M-circle dot (using the Advanced LIGO's design noise curve). By contrast, when the (2,1), (3,3), (4,4), (5,5) modes are also included in the model, the faithfulness improves to 99% for all but four configurations in the numerical-relativity catalog, for which the faithfulness is greater than 98.5%. Starting from the complete inspiral-merger-ringdown model, we develop also a (stand-alone) waveform model for the merger-ringdown signal, calibrated to numerical-relativity waveforms, which can be used to measure multiple quasi-normal modes. The multipolar waveform model can be extended to include spinprecessional effects, and will be employed in upcoming observing runs of Advanced LIGO and Virgo.

We improve the accuracy of the effective-one-body (EOB) waveforms that were employed during the first observing run of Advanced LIGO for binaries of spinning, nonprecessing black holes by calibrating them to a set of 141 numerical-relativity (NR) waveforms. The NR simulations expand the domain of calibration toward larger mass ratios and spins, as compared to the previous EOBNR model. Merger-ringdown waveforms computed in black-hole perturbation theory for Kerr spins close to extremal provide additional inputs to the calibration. For the inspiral-plunge phase, we use a Markov-chain Monte Carlo algorithm to efficiently explore the calibration space. For the merger-ringdown phase, we fit the NR signals with phenomenological formulae. After extrapolation of the calibrated model to arbitrary mass ratios and spins, the (dominant-mode) EOBNR waveforms have faithfulness—at design Advanced-LIGO sensitivity—above 99% against all the NR waveforms, including 16 additional waveforms used for validation, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling for these binary configurations. We find that future NR simulations at mass ratios ≳4 and double spin ≳0.8 will be crucial to resolving discrepancies between different ways of extrapolating waveform models. We also find that some of the NR simulations that already exist in such region of parameter space are too short to constrain the low-frequency portion of the models. Finally, we build a reduced-order version of the EOBNR model to speed up waveform generation by orders of magnitude, thus enabling intensive data-analysis applications during the upcoming observation runs of Advanced LIGO.

The Fokker action of point-particle binaries at the fourth post-Newtonian (4PN) approximation of general relativity has been determined previously. However two ambiguity parameters associated with infra-red (IR) divergencies of spatial integrals had to be introduced. These two parameters were fixed by comparison with gravitational self-force (GSF) calculations of the conserved energy and periastron advance for circular orbits in the test-mass limit. In the present paper together with a companion paper, we determine both these ambiguities from first principle, by means of dimensional regularization. Our computation is thus entirely defined within the dimensional regularization scheme, for treating at once the IR and ultra-violet (UV) divergencies. In particular, we obtain crucial contributions coming from the Einstein-Hilbert part of the action and from the non-local tail term in arbitrary dimensions, which resolve the ambiguities.

In this paper, we revisit and complete our preceding work on the Fokker Lagrangian describing the dynamics of compact binary systems at the fourth post-Newtonian (4PN) order in harmonic coordinates. We clarify the impact of the nonlocal character of the Fokker Lagrangian or the associated Hamiltonian on both the conserved energy and the relativistic periastron precession for circular orbits. We show that the nonlocality of the action, due to the presence of the tail effect at the 4PN order, gives rise to an extra contribution to the conserved integral of energy with respect to the Hamiltonian computed on shell, which was not taken into account in our previous work. We also provide a direct derivation of the periastron advance by taking carefully into account this nonlocality. We then argue that the infrared (IR) divergences in the calculation of the gravitational part of the action are problematic, which motivates us to introduce a second ambiguity parameter, in addition to the one already assumed previously. After fixing these two ambiguity parameters by requiring that the conserved energy and the relativistic periastron precession for circular orbits be in agreement with numerical and analytical gravitational self-force calculations, valid in the limiting case of small mass ratio, we find that our resulting Lagrangian is physically equivalent to the one obtained in the ADM Hamiltonian approach.

We review our recent derivation of a Fokker action describing the conservative dynamics of a compact binary system at the fourth post-Newtonian (4PN) approximation of general relativity. The two bodies are modeled by point particles, which induces ultraviolet (UV) divergences that are cured by means of dimensional regularization combined with a renormalization of the particle's wordlines. Associated with the propagation of wave tails at infinity is the appearance of a non-local-in-time conservative tail effect at the 4PN order in the Lagrangian. In turn this implies the appearance of infrared (IR) divergent integrals which are also regularized by means of dimensional regularization. We compute the Noetherian conserved energy and periastron advance for circular orbits at 4PN order, paying special attention to the treatment of the non-local terms. One ambiguity parameter remaining in the current formalism is determined by comparing those quantities, expressed as functions of the orbital frequency, with self-force results valid in the small mass ratio limit.

We present a new frequency-domain phenomenological model of the gravitational-wave signal from the inspiral, merger and ringdown of nonprecessing (aligned-spin) black-hole binaries. The model is calibrated to 19 hybrid effective-one-body-numerical-relativity waveforms up to mass ratios of 1:18 and black-hole spins of vertical bar a/m vertical bar similar to 0.85 (0.98 for equal-mass systems). The inspiral part of the model consists of an extension of frequency-domain post-Newtonian expressions, using higher-order terms fit to the hybrids. The merger ringdown is based on a phenomenological ansatz that has been significantly improved over previous models. The model exhibits mismatches of typically less than 1% against all 19 calibration hybrids and an additional 29 verification hybrids, which provide strong evidence that, over the calibration region, the model is sufficiently accurate for all relevant gravitational-wave astronomy applications with the Advanced LIGO and Virgo detectors. Beyond the calibration region the model produces physically reasonable results, although we recommend caution in assuming that any merger-ringdown waveform model is accurate outside its calibration region. As an example, we note that an alternative nonprecessing model, SEOBNRv2 (calibrated up to spins of only 0.5 for unequal-mass systems), exhibits mismatch errors of up to 10% for high spins outside its calibration region. We conclude that waveform models would benefit most from a larger number of numerical-relativity simulations of high-aligned-spin unequal-mass binaries.

Current searches for gravitational waves from compact-object binaries with the LIGO and Virgo observatories employ waveform models with spins aligned (or antialigned) with the orbital angular momentum. Here, we derive a new statistic to search for compact objects carrying generic (precessing) spins. Applying this statistic, we construct banks of both aligned-and generic-spin templates for binary black holes and neutron star-black hole binaries, and compare the effectualness of these banks towards simulated populations of generic-spin systems. We then use these banks in a pipeline analysis of Gaussian noise to measure the increase in background incurred by using generic-instead of aligned-spin banks. Although the generic-spin banks have roughly a factor of ten more templates than the aligned-spin banks, we find an overall improvement in signal recovery at a fixed false-alarm rate for systems with high-mass ratio and highly precessing spins. This gain in sensitivity comes at a small loss of sensitivity (less than or similar to 4%) for systems that are already well covered by aligned-spin templates. Since the observation of even a single binary merger with misaligned spins could provide unique astrophysical insights into the formation of these sources, we recommend that the method described here be developed further to mount a viable search for generic-spin binary mergers in LIGO/Virgo data.

The Fokker action governing the motion of compact binary systems without spins is derived in harmonic coordinates at the fourth post-Newtonian approximation (4PN) of general relativity. Dimensional regularization is used for treating the local ultraviolet (UV) divergences associated with point particles, followed by a renormalization of the poles into a redefinition of the trajectories of the point masses. Effects at the 4PN order associated with wave tails propagating at infinity are included consistently at the level of the action. A finite part procedure based on analytic continuation deals with the infrared (IR) divergencies at spatial infinity, which are shown to be fully consistent with the presence of near-zone tails. Our end result at 4PN order is Lorentz invariant and has the correct self-force limit for the energy of circular orbits. However, we find that it differs from the recently published result derived within the ADM Hamiltonian formulation of general relativity [T. Damour, P. Jaranowski, and G. Schafer, Phys. Rev. D 89, 064058 (2014)]. More work is needed to understand this discrepancy.

In this paper we discuss the anatomy of frequency-domain gravitational-wave signals from nonprecessing black-hole coalescences with the goal of constructing accurate phenomenological waveform models. We first present new numerical-relativity simulations for mass ratios up to 18, including spins. From a comparison of different post-Newtonian approximants with numerical-relativity data we select the uncalibrated SEOBNRv2 model as the most appropriate for the purpose of constructing hybrid postNewtonian/ numerical-relativity waveforms, and we discuss how we prepare time-domain and frequencydomain hybrid data sets. We then use our data together with results in the literature to calibrate simple explicit expressions for the final spin and radiated energy. Equipped with our prediction for the final state we then develop a simple and accurate merger-ringdown model based on modified Lorentzians in the gravitational-wave amplitude and phase, and we discuss a simple method to represent the low frequency signal augmenting the TaylorF2 post-Newtonian approximant with terms corresponding to higher orders in the post-Newtonian expansion. We finally discuss different options for modelling the small intermediate frequency regime between inspiral and merger ringdown. A complete phenomenological model based on the present work is presented in a companion paper [S. Khan et al., following paper, Phys. Rev. D 93 044007 (2016)].

Gravitational waveforms which describe the inspiral, merger and ringdown of coalescing binaries are usually constructed by synthesising information from perturbative descriptions, in particular post-Newtonian theory and black-hole perturbation theory, with numerical solutions of the full Einstein equations. In this paper we discuss the "glueing" of numerical and post-Newtonian waveforms to produce hybrid waveforms which include subdominant spherical harmonics ("higher order modes"), and focus in particular on the process of consistently aligning the waveforms, which requires a comparison of both descriptions and a discussion of their imprecisions. We restrict to the non-precessing case, and illustrate the process using numerical waveforms of up to mass ratio $q=18$ produced with the BAM code, and publicly available waveforms from the SXS catalogue. The results also suggest new ways of analysing finite radius errors in numerical simulations.

We investigate the dynamics of spinning binaries of compact objects at the next-to-leading order in the quadratic-in-spin effects, which corresponds to the third post-Newtonian order (3PN). Using a Dixon-type multipolar formalism for spinning point particles endowed with spin-induced quadrupoles and computing iteratively in harmonic coordinates the relevant pieces of the PN metric within the near zone, we derive the post-Newtonian equations of motion as well as the equations of spin precession. We find full equivalence with available results. We then focus on the far-zone field produced by those systems and obtain the previously unknown 3PN spin contributions to the gravitational-wave energy flux by means of the multipolar post-Minkowskian wave generation formalism. Our results are presented in the center-of-mass frame for generic orbits, before being further specialized to the case of spin-aligned, circular orbits. Based on the energy balance equation, we derive the orbital phase of the binary at the order 3PN including all conservative spin effects and briefly discuss the relevance of the new terms.

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, χ(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 compute the next-to-next-to-leading order spin–orbit contributions in the total energy flux emitted in gravitational waves by compact binary systems. Such contributions correspond to the post-Newtonian order 3.5PN for maximally spinning compact objects. Continuing our recent work on the next-to-next-to-leading spin–orbit terms at the 3.5PN order in the equations of motion, we obtain the spin–orbit terms in the multipole moments of the compact binary system up to the same order within the multipolar post-Newtonian wave generation formalism. Our calculation of the multipole moments is valid for general orbits and in an arbitrary frame, the moments are then reduced to the center-of-mass frame and the resulting energy flux is specialized to quasi-circular orbits. The test-mass limit of our final result for the flux agrees with the already known Kerr black hole perturbation limit. Furthermore, the various multipole moments of the compact binary reduce in the one-body case to those of a single-boosted Kerr black hole. We briefly discuss the implications of our result for the gravitational wave flux in terms of the binary's phase evolution, and address its importance for the future detection and parameter estimation of signals in gravitational wave detectors.

We compute next-to-next-to-leading order spin contributions to the post-Newtonian equations of motion for binaries of compact objects, such as black holes or neutron stars. For maximally spinning black holes, those contributions are of third-and-a-half post-Newtonian (3.5 PN) order, improving our knowledge of the equations of motion, already known for non-spinning objects up to this order. Building on previous work, we represent the rotation of the two bodies using a pole–dipole matter stress–energy tensor and iterate Einstein's field equations for a set of potentials parametrizing the metric in harmonic coordinates. Checks of the result include the existence of a conserved energy, the approximate global Lorentz invariance of the equations of motion in harmonic coordinates and the recovery of the motion of a spinning object on a Kerr background in the test-mass limit. We verified the existence of a contact transformation, together with a redefinition of the spin variables that makes our result equivalent to a previously published reduced Hamiltonian, obtained from the Arnowitt–Deser–Misner formalism.

We report on recent results obtained in the post-Newtonian framework for the modelling of the gravitational waves emitted by binary systems of spinning compact objects (black holes and/or neutron stars). These new results are obtained at the spin-orbit (linear-in-spin) level and solving Einstein's field equations iteratively in harmonic coordinates as well as the multipolar post-Newtonian formalism. The dynamics of the binary was tackled at the next-to-next-to-leading order, corresponding to the 3.5 post-Newtonian (PN) order for maximally spinning objects, and the result is found to be consistent with a previously obtained reduced Hamiltonian in the ADM approach. The corresponding contribution to the energy flux emitted by the binary was obtained at the 3.5PN order, as well as the next-to-leading 4PN tail contribution to this flux, an imprint of the non-linearity in the propagation of the wave. These new terms can be used to build more accurate PN templates for the next generation of gravitational wave detectors. We give an illustrative estimate of the quantitative relevance of the new terms in the orbital phasing of the binary.

We extend our previous work devoted to the computation of the next-to-next-to-leading order spin–orbit correction (corresponding to 3.5PN order) in the equations of motion of spinning compact binaries by (i) deriving the corresponding spin–orbit terms in the evolution equations for the spins, the conserved integrals of the motion and the metric regularized at the location of the particles (obtaining also the metric all over the near zone but with some lower precision); (ii) performing the orbital reduction of the precession equations, near-zone metric and conserved integrals to the center-of-mass frame and then further assuming quasi-circular orbits (neglecting gravitational radiation reaction). The results are systematically expressed in terms of the spin variables with a conserved Euclidean norm instead of the original antisymmetric spin tensors of the pole–dipole formalism. This work paves the way to the future computation of the next-to-next-to-leading order spin–orbit terms in the gravitational-wave phasing of spinning compact binaries.

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 build-up to the era of GW astronomy. We present such a model, “PhenomP”, which captures the basic phenomenology of the seven-dimensional parameter space of binary configurations with only three physical parameters. Essentially, we simply “twist up” a twoparameter non-precessing-binary model with approximate expressions for the precessional motion, which require only one additional physical parameter, an effective precession spin. 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, perhaps most importantly, as a simple conceptual framework to form the basis of generic-binary waveform modelling in the advanced-detector era.

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.