M. Pürrer

Frequency-domain gravitational waves from non-precessing black-hole binaries. II. A phenomenological model for the advanced detector era

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 |a/m|∼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.

Simple model of complete precessing black-hole-binary gravitational waveforms

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 binarys 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 models 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.

Frequency-domain gravitational waves from nonprecessing black-hole binaries. I. New numerical waveforms and anatomy of the signal

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 post-Newtonian/numerical-relativity waveforms, and we discuss how we prepare time-domain and frequency-domain 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)].

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.

Frequency domain reduced order model of aligned-spin effective-one-body waveforms with generic mass ratios and spins

I provide a frequency domain reduced order model (ROM) for the aligned-spin effective-one-body (EOB) model “SEOBNRv2” for data analysis with second and third generation ground based gravitational wave (GW) detectors. SEOBNRv2 models the dominant mode of the GWs emitted by the coalescence of black hole (BH) binaries. The large physical parameter space (dimensionless spins −1 ≤ χi ≤ 0.99 and symmetric mass-ratios 0.01 ≤ η ≤ 0.25) requires sophisticated reduced order modeling techniques, including patching in the parameter space and in frequency. I find that the time window over which the inspiral-plunge and the merger-ringdown waveform in SEOBNRv2 are connected is discontinuous when the spin of the deformed Kerr BH χ = 0.8 or the symmetric mass-ratio η ∼ 0.083. This discontinuity increases resolution requirements for the ROM. The ROM can be used for compact binary systems with total masses of 2M or higher for the advanced LIGO (aLIGO) design sensitivity and a 10 Hz lower cutoff frequency. The ROM has a worst mismatch against SEOBNRv2 of ∼ 1%, but in general mismatches are better than ∼ 0.1%. The ROM is crucial for key data analysis applications for compact binaries, such as GW searches and parameter estimation carried out within the LIGO Scientific Collaboration (LSC). LIGO-P1500175

Fast and accurate inference on gravitational waves from precessing compact binaries

Inferring astrophysical information from gravitational waves emitted by compact binaries is one of the key science goals of gravitational-wave astronomy. In order to reach the full scientific potential of gravitational-wave experiments, we require techniques to mitigate the cost of Bayesian inference, especially as gravitational-wave signal models and analyses become increasingly sophisticated and detailed. Reduced-order models (ROMs) of gravitational waveforms can significantly reduce the computational cost of inference by removing redundant computations. In this paper, we construct the first reduced-order models of gravitational-wave signals that include the effects of spin precession, inspiral, merger, and ringdown in compact object binaries and that are valid for component masses describing binary neutron star, binary black hole, and mixed binary systems. This work utilizes the waveform model known as “IMRPhenomPv2.” Our ROM enables the use of a fast reduced-order quadrature (ROQ) integration rule which allows us to approximate Bayesian probability density functions at a greatly reduced computational cost. We find that the ROQ rule can be used to speed-up inference by factors as high as 300 without introducing systematic bias. This corresponds to a reduction in computational time from around half a year to half a day for the longest duration and lowest mass signals. The ROM and ROQ rules are available with the main inference library of the LIGO Scientific Collaboration, LALInference.

Hierarchical data-driven approach to fitting numerical relativity data for nonprecessing binary black holes with an application to final spin and radiated energy

Numerical relativity is an essential tool in studying the coalescence of binary black holes (BBHs). It is still computationally prohibitive to cover the BBH parameter space exhaustively, making phenomenological fitting formulas for BBH waveforms and final-state properties important for practical applications. We describe a general hierarchical bottom-up fitting methodology to design and calibrate fits to numerical relativity simulations for the three-dimensional parameter space of quasicircular nonprecessing merging BBHs, spanned by mass ratio and by the individual spin components orthogonal to the orbital plane. Particular attention is paid to incorporating the extreme-mass-ratio limit and to the subdominant unequal-spin effects. As an illustration of the method, we provide two applications, to the final spin and final mass (or equivalently: radiated energy) of the remnant black hole. Fitting to 427 numerical relativity simulations, we obtain results broadly consistent with previously published fits, but improving in overall accuracy and particularly in the approach to extremal limits and for unequal-spin configurations. We also discuss the importance of data quality studies when combining simulations from diverse sources, how detailed error budgets will be necessary for further improvements of these already highly accurate fits, and how this first detailed study of unequal-spin effects helps in choosing the most informative parameters for future numerical relativity runs.

Impact of gravitational radiation higher order modes on single aligned-spin gravitational wave searches for binary black holes

Current template-based gravitational wave searches for compact binary coalescences (CBC) use waveform models that neglect the higher order modes content of the gravitational radiation emitted, considering only the quadrupolar

Can we measure individual black-hole spins from gravitational-wave observations?

(\ell,|m|)=(2,2)

Parametrized tests of the strong-field dynamics of general relativity using gravitational wave signals from coalescing binary black holes: Fast likelihood calculations and sensitivity of the method

modes. We study the effect of such a neglection for the case of aligned-spin CBC searches for equal-spin (and non-spinning) binary black holes in the context of two versions of Advanced LIGO: the upcoming 2015 version, known as early Advanced LIGO (eaLIGO) and its Zero-Detuned High Energy Power version, that we will refer to as Advanced LIGO (AdvLIGO). In addition, we study the case of a non-spinning search for initial LIGO (iLIGO). We do this via computing the effectualness of the aligned-spin SEOBNRv1 ROM waveform family, which only considers quadrupolar modes, towards hybrid post-Newtonian/Numerical Relativity waveforms which contain higher order modes. We find that for all LIGO versions, losses of more than