Search

Binary black-hole systems are expected to be important sources of gravitational waves for upcoming gravitational-wave detectors. If the spins are not colinear with each other or with the orbital angular momentum, these systems exhibit complicated precession dynamics that are imprinted on the gravitational waveform. We develop a new procedure to match the precession dynamics computed by post-Newtonian (PN) theory to those of numerical binary black-hole simulations in full general relativity. For numerical relativity (NR) simulations lasting approximately two precession cycles, we find that the PN and NR predictions for the directions of the orbital angular momentum and the spins agree to better than ~1° with NR during the inspiral, increasing to 5° near merger. Nutation of the orbital plane on the orbital time scale agrees well between NR and PN, whereas nutation of the spin direction shows qualitatively different behavior in PN and NR. We also examine how the PN equations for precession and orbital-phase evolution converge with PN order, and we quantify the impact of various choices for handling partially known PN terms. [ABSTRACT FROM AUTHOR]

The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries (BHNS), very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin nonprecessing binaries but limited accuracy, and a much smaller number of longer, more recent simulations limited to nonspinning black holes. In this paper, we present long, accurate numerical waveforms for three new systems that include rapidly spinning black holes, and one precessing configuration. We study in detail the accuracy of the simulations, and in particular perform for the first time in the context of BHNS binaries a detailed comparison of waveform extrapolation methods to the results of Cauchy characteristic extraction. The new waveforms have < 0.1     rad phase errors during inspiral, rising to ∼ (0.2-0.4)     rad errors at merger, and ≲1 % error in their amplitude. We compute the faithfulness of recent analytical models to these numerical results for the late inspiral and merger phases covered by the numerical simulations, and find that models specifically designed for BHNS binaries perform well (faithfulness F>0.99) for binaries seen face on. For edge-on observations, particularly for precessing systems, disagreements between models and simulations increase, and models that include precession and/or higher-order modes start to perform better than BHNS models that currently lack these features. [ABSTRACT FROM AUTHOR]

Characterizing the properties of black holes is one of the most important science objectives for gravitational-wave observations. Astrophysical evidence suggests that black holes that are nearly extremal (i.e., spins near the theoretical upper limit) might exist and, thus, might be among the merging black holes observed with gravitational waves. In this paper, we explore how well current gravitational wave parameter estimation methods can measure the spins of rapidly spinning black holes in binaries. We simulate gravitational-wave signals using numerical-relativity waveforms for nearly-extremal, merging black holes. For simplicity, we confine our attention to binaries with spins parallel or antiparallel with the orbital angular momentum. We find that recovering the holes' nearly extremal spins is challenging. When the spins are nearly extremal and parallel to each other, the resulting parameter estimates do recover spins that are large, though the recovered spin magnitudes are still significantly smaller than the true spin magnitudes. When the spins are nearly extremal and antiparallel to each other, the resulting parameter estimates recover the small effective spin but incorrectly estimate the individual spins as nearly zero. We study the effect of spin priors and argue that a commonly used prior (uniform in spin magnitude and direction) hinders unbiased recovery of large black-hole spins. [ABSTRACT FROM AUTHOR]

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 200AiQ, 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. [ABSTRACT FROM AUTHOR]

The detection of gravitational waves and the extraction of physical information from them requires the prediction of accurate waveforms to be used in template banks. For that purpose, the accuracy of effective-one-body (EOB) waveforms has been improved over the last years by calibrating them to numerical-relativity (NR) waveforms. So far, the calibration has employed a handful of NR waveforms with a total length of ~30 cycles, the length being limited by the computational cost of NR simulations. Here we address the outstanding problem of the stability of the EOB calibration with respect to the length of NR waveforms. Performing calibration studies against NR waveforms of nonspinning black-hole binaries with mass ratios 1, 1.5, 5, and 8, and with a total length of ~60 cycles, we find that EOB waveforms calibrated against either 30 or 60 cycles will be indistinguishable by the advanced detectors LIGO and Virgo when the signal-to-noise ratio (SNR) is below 110. When extrapolating to a very large number of cycles, using very conservative assumptions, we can conclude that state-of-the-art nonspinning EOB waveforms of any length are sufficiently accurate for parameter estimation with advanced detectors when the SNR is below 20, the mass ratio is below 5 and total mass is above 20 M☉. The results are not conclusive for the entire parameter space because of current NR errors. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 (ℓ, 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 20Mo and 2OOM0, 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 ant ¡aligned spins. Extension of this model to black holes with aligned spins x¡ --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. [ABSTRACT FROM AUTHOR]

We present a new study of remnant black hole properties from 13 binary black hole systems, numerically evolved using the Spectral Einstein Code. The mass, spin, and recoil velocity of each remnant were determined quasilocally from apparent horizon data and asymptotically from Bondi data (h,ψ4,ψ3,ψ2,ψ1) computed at future null infinity using SpECTRE's Cauchy characteristic evolution. We compare these independent measurements of the remnant properties in the bulk and on the boundary of the spacetime, giving insight into how well asymptotic data are able to reproduce local properties of the remnant black hole in numerical relativity. We also discuss the theoretical framework for connecting horizon quantities to asymptotic quantities and how it relates to our results. This study recommends a simple improvement to the recoil velocities reported in the Simulating eXtreme Spacetimes waveform catalog, provides an improvement to future surrogate remnant models, and offers new analysis techniques for evaluating the physical accuracy of numerical simulations. [ABSTRACT FROM AUTHOR]

After the discovery of gravitational waves from binary black holes (BBHs) and binary neutron stars (BNSs) with the LIGO and Virgo detectors, neutron-star black holes (NSBHs) are the natural next class of binary systems to be observed. In this work, we develop a waveform model for aligned-spin NSBHs combining a BBH baseline waveform (available in the effective-one-body approach) with a phenomenological description of tidal effects (extracted from numerical-relativity simulations) and correcting the amplitude during the late inspiral, merger and ringdown to account for the NS tidal disruption. In particular, we calibrate the amplitude corrections using NSBH waveforms obtained with the numerical-relativity spectral Einstein code (spec) and the sacra code. The model was calibrated using simulations with NS masses in the range 1.2-1.4 M⊙, tidal deformabilities up to 4200 (for a 1.2 M⊙ NS), and dimensionless BH spin magnitude up to 0.9. Based on the simulations used and on checking that sensible waveforms are produced, we recommend our model to be employed with a NS mass in the range 1-3 M⊙, tidal deformability 0-5000, and (dimensionless) BH spin magnitude up to 0.9. We also validate our model against two new, highly accurate NSBH waveforms with BH spin 0.9 and mass ratios 3 and 4, characterized by tidal disruption, produced with SpEC, and find very good agreement. Furthermore, we compute the unfaithfulness between waveforms from NSBH, BBH, and BNS systems, finding that it will be challenging for the Advanced LIGO-Virgo detector network at design sensitivity to distinguish different source classes. We perform a Bayesian parameter-estimation analysis on a synthetic numerical-relativity signal in zero noise to study parameter biases. Finally, we reanalyze GW170817, with the hypothesis that it is a NSBH. We do not find evidence to distinguish the BNS and NSBH hypotheses; however, the posterior for the mass ratio is shifted to less equal masses under the NSBH hypothesis. [ABSTRACT FROM AUTHOR]

Gravitational-wave (GW) detectors have begun to observe coalescences of heavy black hole binaries (M≳50 M⊙) at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters, such as masses and spins of component black holes. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations provide ideal models for constraining physical parameters describing these heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogate models as reliable parameter estimation tools. We show that within a fully Bayesian framework, NR surrogates can help extract additional information from GW observations that is inaccessible to traditional models. We demonstrate this by analyzing a set of synthetic signals with NR surrogate templates and comparing against contemporary approximate models. We then consider the case of two of the earliest binary black holes detected by the LIGO observatories, GW150914 and GW170104. We reanalyze their data with the generically precessing NR-based surrogate model and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. Our analysis estimates the sources of both events to be 20%-25% further away than was previously estimated. Our analysis also constrains their orbital orientations more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin χeff more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It strongly motivates that NR surrogates be extended to cover more of the binary black hole parameter space. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]