Search

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]

Coalescing binaries of neutron stars and black holes are one of the most important sources of gravitational waves for the upcoming network of ground-based detectors. Detection and extraction of astrophysical information from gravitational-wave signals requires accurate waveform models. The effective-one-body and other phenomenological models interpolate between analytic results and numerical relativity simulations, that typically span O(10) orbits before coalescence. In this paper we study the faithfulness of these models for neutron star-black hole binaries. We investigate their accuracy using new numerical relativity (NR) simulations that span 36-88 orbits, with mass ratios q and black hole spins j BH of (q,Xbh) = (7, ±0.4), (7, ±0.6), and (5, -0.9). These simulations were performed treating the neutron star as a low-mass black hole, ignoring its matter effects. We find that (i) the recently published SEOBNRvl and SEOBNRv2 models of the effective-one-body family disagree with each other (mismatches of a few percent) for black hole spins XBH ≥ 0.5 or XBH ≤ -0.3, with waveform mismatch accumulating during early inspiral; (ii) comparison with numerical waveforms indicates that this disagreement is due to phasing errors of SEOBNRvl, with SEOBNRv2 in good agreement with all of our simulations; (iii) phenomenological waveforms agree with SEOBNRv2 only for comparable-mass low-spin binaries, with overlaps below 0.7 elsewhere in the neutron star-black hole binary parameter space; (iv) comparison with numerical waveforms shows that most of this model's dephasing accumulates near the frequency interval where it switches to a phenomenological phasing prescription; and finally (v) both SEOBNR and post-Newtonian models are effectual for neutron star-black hole systems, but post-Newtonian waveforms will give a significant bias in parameter recovery. Our results suggest that future gravitational-wave detection searches and parameter estimation efforts would benefit from using SEOBNRv2 waveform templates when focused on neutron star-black hole systems with q ≲ 7 and XBH ≈ [-0.9, +0.6]. For larger black hole spins and/or binary mass ratios, we recommend the models be further investigated as NR simulations in that region of the parameter space become available. [ABSTRACT FROM AUTHOR]

In this work we study the dynamics of spinning binary black hole systems in the strong field regime. For this purpose we extract from numerical relativity simulations the binding energy, specific orbital angular momentum, and gauge-invariant frequency. The goal of our work is threefold: First, we extract the individual spin contributions to the binding energy, in particular the spin-orbit, spin-spin, and cubic-in-spin terms. Second, we compare our results with predictions from waveform models and find that while post-Newtonian approximants are not capable of representing the dynamics during the last few orbits before merger, there is good agreement between our data and effective-one-body approximants as well as the numerical relativity surrogate models. Finally, we present phenomenological representations for the binding energy for nonspinning systems with mass ratios up to q=10 and for the spin-orbit interaction for mass ratios up to q=8 obtaining accuracies of ≲0.1% and ≲6%, respectively. [ABSTRACT FROM AUTHOR]

Spin-orbit resonances have important astrophysical implications as the evolution and subsequent coalescence of supermassive black hole binaries in one of these configurations may lead to low recoil velocity of merger remnants. It has also been shown that black hole spins in comparable mass stellar-mass black hole binaries could preferentially lie in a resonant plane when their gravitational waves (GWs) enter the advanced LIGO frequency band [1]. Therefore, it is highly desirable to investigate the possibility of detection and subsequent characterization of such GW sources in the advanced detector era, which can, in turn, improve our perception of their high mass counterparts. The current detection pipelines involve only nonprecessing templates for compact binary searches whereas parameter estimation pipelines can afford to use approximate precessing templates. In this paper, we test the performance of these templates in detection and characterization of spin-orbit resonant binaries. We use fully precessing time-domain SEOBNRv3 waveforms as well as four numerical relativity (NR) waveforms to model GWs from spin-orbit resonant binaries and filter them through IMRPhenomD, SEOBNRv4 and IMRPhenomPv2 approximants. We find that the nonprecessing approximants IMRPhenomD and SEOBNRv4 recover only ~70% of injections with fitting factor (FF) higher than 0.97 (or 90% of injections with FF>0.9). This loss in signal-to-noise ratio is mainly due to the missing physics in these approximants in terms of precession and nonquadrupole modes. However, if we use a new statistic, i.e., maximizing the matched filter output over the sky-location parameters as well, the precessing approximant IMRPhenomPv2 performs magnificently better than their nonprecessing counterparts with recovering 99% of the injections with FFs higher than 0.97. Interestingly, injections with Δϕ=180° have higher FFs (Δϕ is the angle between the components of the black hole spins in the plane orthogonal to the orbital angular momentum) as compared to their Δϕ=0° and generic counterparts. This is because Δϕ=180° binaries are not as strongly precessing as Δϕ=0° and generic binaries. This implies that we will have a slight observation bias towards Δϕ=180° and away from Δϕ=0° resonant binaries while using nonprecessing templates for searches. Moreover, all template approximants are able to recover most of the injected NR waveforms with FFs >0.95. For all the injections including NR, the systematic error in estimating chirp mass remains below <10% with minimum error for Δϕ=180° resonant binaries. The symmetric mass-ratio can be estimated with errors below 15%. The effective spin parameter χeff is measured with maximum absolute error of 0.13. The in-plane spin parameter χp is mostly underestimated indicating that a precessing signal will be recovered as a relatively less precessing signal. Based on our findings, we conclude that we not only need improvements in waveform models towards precession and nonquadrupole modes but also better search strategies for precessing GW signals. [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]