Nicholas W. Taylor

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.

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.

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.

Dynamical excision boundaries in spectral evolutions of binary black hole spacetimes

Simulations of binary black hole systems using the Spectral Einstein Code (SpEC) are done on a computational domain that excises the regions inside the black holes. It is imperative that the excision boundaries are outflow boundaries with respect to the hyperbolic evolution equations used in the simulation. We employ a time-dependent mapping between the fixed computational frame and the inertial frame through which the black holes move. The time-dependent parameters of the mapping are adjusted throughout the simulation by a feedback control system in order to follow the motion of the black holes, to adjust the shape and size of the excision surfaces so that they remain outflow boundaries, and to prevent large distortions of the grid. We describe in detail the mappings and control systems that we use. We show how these techniques have been essential in the evolution of binary black hole systems with extreme configurations, such as large spin magnitudes and high mass ratios, especially during the merger, when apparent horizons are highly distorted and the computational domain becomes compressed. The techniques introduced here may be useful in other applications of partial differential equations that involve time-dependent mappings.

Comparing Gravitational Waveform Extrapolation to Cauchy-Characteristic Extraction in Binary Black Hole Simulations

We extract gravitational waveforms from numerical simulations of black hole binaries computed using the Spectral Einstein Code. We compare two extraction methods: direct construction of the Newman-Penrose (NP) scalar Ψ_4 at a finite distance from the source and Cauchy-characteristic extraction (CCE). The direct NP approach is simpler than CCE, but NP waveforms can be contaminated by near-zone effects—unless the waves are extracted at several distances from the source and extrapolated to infinity. Even then, the resulting waveforms can in principle be contaminated by gauge effects. In contrast, CCE directly provides, by construction, gauge-invariant waveforms at future null infinity. We verify the gauge invariance of CCE by running the same physical simulation using two different gauge conditions. We find that these two gauge conditions produce the same CCE waveforms but show differences in extrapolated-Ψ_4 waveforms. We examine data from several different binary configurations and measure the dominant sources of error in the extrapolated-Ψ4 and CCE waveforms. In some cases, we find that NP waveforms extrapolated to infinity agree with the corresponding CCE waveforms to within the estimated error bars. However, we find that in other cases extrapolated and CCE waveforms disagree, most notably for m=0 “memory” modes.

Final spin and radiated energy in numerical simulations of binary black holes with equal masses and equal, aligned or antialigned spins

The behavior of merging black holes (including the emitted gravitational waves and the properties of the remnant) can currently be computed only by numerical simulations. This paper introduces ten numerical relativity simulations of binary black holes with equal masses and equal spins aligned or antialigned with the orbital angular momentum. The initial spin magnitudes have |χ_i|≲0.95 and are more concentrated in the aligned direction because of the greater astrophysical interest of this case. We combine these data with five previously reported simulations of the same configuration, but with different spin magnitudes, including the highest spin simulated to date, χ_i≈0.97. This data set is sufficiently accurate to enable us to offer improved analytic fitting formulas for the final spin and for the energy radiated by gravitational waves as a function of initial spin. The improved fitting formulas can help to improve our understanding of the properties of binary black hole merger remnants and can be used to enhance future approximate waveforms for gravitational wave searches, such as effective-one-body waveforms.

What does a binary black hole merger look like

We present a method of calculating the strong-field gravitational lensing caused by many analytic and numerical spacetimes. We use this procedure to calculate the distortion caused by isolated black holes (BHs) and by numerically evolved BH binaries. We produce both demonstrative images illustrating details of the spatial distortion and realistic images of collections of stars taking both lensing amplification and redshift into account. On large scales the lensing from inspiraling binaries resembles that of single BHs, but on small scales the resulting images show complex and in some cases self-similar structure across different angular scales.

Suitability of hybrid gravitational waveforms for unequal-mass binaries

This article studies sufficient accuracy criteria of hybrid post-Newtonian (PN) and numerical relativity (NR) waveforms for parameter estimation of strong binary black-hole sources in second-generation ground-based gravitational-wave detectors. We investigate equal-mass nonspinning binaries with a new 33-orbit NR waveform, as well as unequal-mass binaries with mass ratios 2, 3, 4 and 6. For equal masses, the 33-orbit NR waveform allows us to recover previous results and to extend the analysis toward matching at lower frequencies. For unequal masses, the errors between different PN approximants increase with mass ratio. Thus, at 3.5 PN, hybrids for higher-mass-ratio systems would require NR waveforms with many more gravitational-wave cycles to guarantee no adverse impact on parameter estimation. Furthermore, we investigate the potential improvement in hybrid waveforms that can be expected from fourth-order post-Newtonian waveforms and find that knowledge of this fourth post-Newtonian order would significantly improve the accuracy of hybrid waveforms.

Periastron advance in spinning black hole binaries: comparing effective-one-body and Numerical Relativity

We compute the periastron advance using the effective-one-body formalism for binary black holes moving on quasicircular orbits and having spins collinear with the orbital angular momentum. We compare the predictions with the periastron advance recently computed in accurate numerical-relativity simulations and find remarkable agreement for a wide range of spins and mass ratios. These results do not use any numerical-relativity calibration of the effective-one-body model, and stem from two key ingredients in the effective-one-body Hamiltonian: (i) the mapping of the two-body dynamics of spinning particles onto the dynamics of an effective spinning particle in a (deformed) Kerr spacetime, fully symmetrized with respect to the two-body masses and spins, and (ii) the resummation, in the test-particle limit, of all post-Newtonian corrections linear in the spin of the particle. In fact, even when only the leading spin post-Newtonian corrections are included in the effective-one-body spinning Hamiltonian but all the test-particle corrections linear in the spin of the particle are resummed we find very good agreement with the numerical results (within the numerical error for equal-mass binaries and discrepancies of at most 1% for larger mass ratios). Furthermore, we specialize to the extreme mass-ratio limit and derive, using the equations of motion in the gravitational skeleton approach, analytical expressions for the periastron advance, the meridional Lense-Thirring precession and spin precession frequency in the case of a spinning particle on a nearly circular equatorial orbit in Kerr spacetime, including also terms quadratic in the spin.