Geoffrey Lovelace

Effective-one-body model for black-hole binaries with generic mass ratios and spins

Gravitational waves emitted by black-hole binary systems have the highest signal-to-noise ratio in LIGO and Virgo detectors when black-hole spins are aligned with the orbital angular momentum and extremal. For such systems, we extend the effective-one-body inspiral-merger-ringdown waveforms to generic mass ratios and spins calibrating them to 38 numerical-relativity nonprecessing waveforms produced by the SXS Collaboration. The numerical-relativity simulations span mass ratios from 1 to 8, spin magnitudes up to 98% of extremality, and last for 40 to 60 gravitational-wave cycles. When the total mass of the binary is between 20 and 200M_⊙, the effective-one-body nonprecessing (dominant mode) waveforms have overlap above 99% (using the advanced-LIGO design noise spectral density) with all of the 38 nonprecessing numerical waveforms, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling. We also show that—without further calibration— the precessing effective-one-body (dominant mode) waveforms have overlap above 97% with two very long, strongly precessing numerical-relativity waveforms, when maximizing only on the initial phase and time.

Catalog of 174 Binary Black Hole Simulations for Gravitational Wave Astronomy

This Letter presents a publicly available catalog of 174 numerical binary black hole simulations following up to 35 orbits. The catalog includes 91 precessing binaries, mass ratios up to 8∶1, orbital eccentricities from a few percent to 10(-5), black hole spins up to 98% of the theoretical maximum, and radiated energies up to 11.1% of the initial mass. We establish remarkably good agreement with post-Newtonian precession of orbital and spin directions for two new precessing simulations, and we discuss other applications of this catalog. Formidable challenges remain: e.g., precession complicates the connection of numerical and approximate analytical waveforms, and vast regions of the parameter space remain unexplored.

Binary-black-hole initial data with nearly-extremal spins

There is a significant possibility that astrophysical black holes with nearly extremal spins exist. Numerical simulations of such systems require suitable initial data. In this paper, we examine three methods of constructing binary-black-hole initial data, focusing on their ability to generate black holes with nearly extremal spins: (i) Bowen-York initial data, including standard puncture data (based on conformal flatness and Bowen-York extrinsic curvature), (ii) standard quasiequilibrium initial data (based on the extended-conformal-thin-sandwich equations, conformal flatness, and maximal slicing), and (iii) quasiequilibrium data based on the superposition of Kerr-Schild metrics. We find that the two conformally flat methods (i) and (ii) perform similarly, with spins up to about 0.99 obtainable at the initial time. However, in an evolution, we expect the spin to quickly relax to a significantly smaller value around 0.93 as the initial geometry relaxes. For quasiequilibrium superposed Kerr-Schild data [method (iii)], we construct initial data with initial spins as large as 0.9997. We evolve superposed Kerr-Schild data sets with spins of 0.93 and 0.97 and find that the spin drops by only a few parts in 10^4 during the initial relaxation; therefore, we expect that superposed Kerr-Schild initial data will allow evolutions of binary black holes with relaxed spins above 0.99. Along the way to these conclusions, we also present several secondary results: the power-law coefficients with which the spin of puncture initial data approaches its maximal possible value; approximate analytic solutions for large spin puncture data; embedding diagrams for single spinning black holes in methods (i) and (ii); nonunique solutions for method (ii). All of the initial-data sets that we construct contain subextremal black holes, and when we are able to push the spin of the excision boundary surface into the superextremal regime, the excision surface is always enclosed by a second, subextremal apparent horizon. The quasilocal spin is measured by using approximate rotational Killing vectors, and the spin is also inferred from the extrema of the intrinsic scalar curvature of the apparent horizon. Both approaches are found to give consistent results, with the approximate-Killing-vector spin showing the least variation during the initial relaxation.

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

The Numerical–Relativity–Analytical–Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binarys total mass is ~100–200M⊙, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios ≤4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.

Prototype effective-one-body model for nonprecessing spinning inspiral-merger-ringdown waveforms

This paper presents a tunable effective-one-body (EOB) model for black-hole (BH) binaries of arbitrary mass ratio and aligned spins. This new EOB model incorporates recent results of small-mass-ratio simulations based on Teukolsky’s perturbative formalism. The free parameters of the model are calibrated to numerical-relativity simulations of nonspinning BH-BH systems of five different mass ratios and to equal-mass nonprecessing BH-BH systems with dimensionless BH spins χ_i≃±0.44. The present analysis focuses on the orbital dynamics of the resulting EOB model, and on the dominant (l,m)=(2,2) gravitational-wave mode. The calibrated EOB model can generate inspiral-merger-ringdown waveforms for nonprecessing, spinning BH binaries with any mass ratio and with individual BH spins -1≤χ_i≲0.7. Extremizing only over time and phase shifts, the calibrated EOB model has overlaps larger than 0.997 with each of the seven numerical-relativity waveforms for total masses between 20M_⊙ and 200M_⊙, using the Advanced LIGO noise curve. We compare the calibrated EOB model with two additional equal-mass highly spinning (χ_i≃-0.95,+0.97) numerical-relativity waveforms, which were not used during calibration. We find that the calibrated model has an overlap larger than 0.995 with the simulation with nearly extremal antialigned spins. Extension of this model to black holes with aligned spins χ_i≳0.7 requires improvements of our modeling of the plunge dynamics and inclusion of higher-order PN spin terms in the gravitational-wave modes and radiation-reaction force.

Improved effective-one-body model of spinning, nonprecessing binary black holes for the era of gravitational-wave astrophysics with advanced detectors

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 NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

The numerical injection analysis (NINJA) project is a collaborative effort between members of the numerical-relativity and gravitational wave data-analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms and to foster closer collaboration between the numerical-relativity and data-analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian—numerical-relativity hybrid waveforms, large numbers of injections and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical-relativity groups have contributed 56 hybrid waveforms consisting of a numerical portion modeling the late inspiral, merger and ringdown stitched to a post-Newtonian portion modeling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the (l, m) = (2, 2) mode. Based on these studies, we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.

Reducing orbital eccentricity in binary black hole simulations

Binary black hole simulations starting from quasi-circular (i.e., zero radial velocity) initial data have orbits with small but nonzero orbital eccentricities. In this paper, the quasi-equilibrium initial-data method is extended to allow nonzero radial velocities to be specified in binary black hole initial data. New low-eccentricity initial data are obtained by adjusting the orbital frequency and radial velocities to minimize the orbital eccentricity, and the resulting (~5 orbit) evolutions are compared with those of quasi-circular initial data. Evolutions of the quasi-circular data clearly show eccentric orbits, with eccentricity that decays over time. The precise decay rate depends on the definition of eccentricity; if defined in terms of variations in the orbital frequency, the decay rate agrees well with the prediction of Peters (1964 Phys. Rev. 136 1224–32). The gravitational waveforms, which contain ~8 cycles in the dominant l = m = 2 mode, are largely unaffected by the eccentricity of the quasi-circular initial data. The overlap between the dominant mode in the quasi-circular evolution and the same mode in the low-eccentricity evolution is about 0.99.

Simulating merging binary black holes with nearly extremal spins

Astrophysically realistic black holes may have spins that are nearly extremal (i.e., close to 1 in dimensionless units). Numerical simulations of binary black holes are important tools both for calibrating analytical templates for gravitational-wave detection and for exploring the nonlinear dynamics of curved spacetime. However, all previous simulations of binary-black-hole inspiral, merger, and ringdown have been limited by an apparently insurmountable barrier: the merging holes’ spins could not exceed 0.93, which is still a long way from the maximum possible value in terms of the physical effects of the spin. In this paper, we surpass this limit for the first time, opening the way to explore numerically the behavior of merging, nearly extremal black holes. Specifically, using an improved initial-data method suitable for binary black holes with nearly extremal spins, we simulate the inspiral (through 12.5 orbits), merger and ringdown of two equal-mass black holes with equal spins of magnitude 0.95 antialigned with the orbital angular momentum.

High-accuracy gravitational waveforms for binary black hole mergers with nearly extremal spins

Motivated by the possibility of observing gravitational waves from merging black holes whose spins are nearly extremal (i.e. 1 in dimensionless units), we present numerical waveforms from simulations of merging black holes with the highest spins simulated to date: (1) a 25.5-orbit inspiral, merger and ringdown of two holes with equal masses and spins of magnitude 0.97 aligned with the orbital angular momentum; and (2) a previously reported 12.5-orbit inspiral, merger and ringdown of two holes with equal masses and spins of magnitude 0.95 anti-aligned with the orbital angular momentum. First, we consider the horizon mass and spin evolution of the new aligned-spin simulation. During the inspiral, the horizon area and spin evolve in remarkably close agreement with Alvis analytic predictions, and the remnant holes final spin agrees reasonably well with several analytic predictions. We also find that the total energy emitted by a real astrophysical system with these parameters—almost all of which is radiated during the time included in this simulation—would be 10.952% of the initial mass at infinite separation. Second, we consider the gravitational waveforms for both simulations. After estimating their uncertainties, we compare the waveforms to several post-Newtonian approximants, finding significant disagreement well before merger, although the phase of the TaylorT4 approximant happens to agree remarkably well with the numerical prediction in the aligned-spin case. We find that the post-Newtonian waveforms have sufficient uncertainty that hybridized waveforms will require far longer numerical simulations (in the absence of improved post-Newtonian waveforms) for accurate parameter estimation of low-mass binary systems.