Nonprecessing Spins Research Articles

Effective-one-body waveforms for extreme-mass-ratio binaries: Consistency with second-order gravitational self-force quasicircular results and extension to nonprecessing spins and eccentricity

Effective-one-body waveforms for extreme-mass-ratio binaries: Consistency with second-order gravitational self-force quasicircular results and extension to nonprecessing spins and eccentricity

Numerical-relativity-informed effective-one-body model for black-hole–neutron-star mergers with higher modes and spin precession

We present the first effective-one-body (EOB) model for generic-spins quasi-circular black-hole--neutron-star (BHNS) inspiral-merger-ringdown gravitational waveforms (GWs). Our model is based on a new numerical-relativity (NR) informed expression of the BH remnant and its ringdown. It reproduces the NR $(\ell,m)=(2,2)$ waveform with typical phase agreement of ${\lesssim0.5}\,$rad (${\lesssim1}\,$rad) to merger (ringdown). The maximum (minimum) mismatch between the $(2,2)$ and the NR data is 4% (0.6%). Higher modes (HMs) $(2,1)$, $(3,2)$, $(3,3)$, $(4,4)$ and $(5,5)$ are included and their mismatch with the available NR waveforms are up to (down to) a 60% (1%) depending on the inclination. Phase comparison with a 16 orbit precessing simulation shows differences within the NR uncertainties. We demonstrate the applicability of the model in GW parameter estimation by perfoming the first BHNS Bayesian analysis with HMs (and non-precessing spins) of the event GW190814, together with new $(2,2)$-mode analysis of GW200105 and GW200115. For the GW190814 study, the inclusion of HMs gives tighter parameter posteriors. The Bayes factors of our analyses on this event show decisive evidence for the presence of HMs, but no clear preference for a BHNS or a binary black hole (BBH) source. Similarly, we confirm GW200105 and GW200115 show no evidence for tidal effects.

Cost minimization in acquisition for gravitational wave surrogate modeling

Gravitational wave science is dependent upon expensive numerical simulations, which provide the foundational understanding of binary merger radiation needed to interpret observations of massive binary black holes. The high cost of these simulations limits large-scale campaigns to explore the binary black hole parameter space. Surrogate models have been developed to efficiently interpolate between simulation results, but these models require a sufficiently comprehensive sample to train on. Acquisition functions can be used to identify points in the domain for simulation. We develop a new acquisition function which accounts for the cost of simulating new points. We show that when applied to a 3D domain of binary mass ratio and dimensionless nonprecessing spins, the accumulated cost of simulation is reduced by a factor of about 10.

Effective-one-body multipolar waveforms for eccentric binary black holes with nonprecessing spins

We construct an inspiral-merger-ringdown eccentric gravitational-wave (GW) model for binary black holes with non-precessing spins within the effective-one-body formalism. This waveform model, SEOBNRv4EHM, extends the accurate quasi-circular SEOBNRv4HM model to eccentric binaries by including recently computed eccentric corrections up to 2PN order in the gravitational waveform modes, notably the $(l,|m|)=(2,2),(2,1),(3,3),(4,4),(5,5)$ multipoles. The waveform model reproduces the zero eccentricity limit with an accuracy comparable to the underlying quasi-circular model, with the unfaithfulness of $\lesssim1\%$ against quasi-circular numerical-relativity (NR) simulations. When compared against 28 public eccentric NR simulations from the Simulating eXtreme Spacetimes catalog with initial orbital eccentricities up to $e\simeq0.3$ and dimensionless spin magnitudes up to $+0.7$, the model provides unfaithfulness $<1\%$, showing that both the $(2,|2|)$-modes and the higher-order modes are reliably described without calibration to NR datasets in the eccentric sector. The waveform model SEOBNRv4EHM is able to qualitatively reproduce the phenomenology of dynamical captures, and can be extended to include spin-precession effects. It can be employed for upcoming observing runs with the LIGO-Virgo-KAGRA detectors and used to re-analyze existing GW catalogs to infer the eccentricity parameters for binaries with $e\lesssim0.3$ (at 20 Hz or lower) and spins up to $\lesssim 0.9-0.95$. The latter is a promising region of the parameter space where some astrophysical formation scenarios of binaries predict mild eccentricity in the ground-based detectors' bandwidth. Assessing the accuracy and robustness of the eccentric waveform model SEOBNRv4EHM for larger eccentricities and spins will require comparisons with, and, likely, calibration to eccentric NR waveforms in a larger region of the parameter space.

Effective-one-body waveforms from dynamical captures in black hole binaries

Dynamical capture is a possible formation channel for BBH mergers leading to highly eccentric merger dynamics and to gravitational wave (GW) signals that are morphologically different from those of quasi-circular mergers. The future detection of these mergers by ground-based or space-based GW interferometers can provide invaluable insights on astrophysical black holes, but it requires precise predictions and dedicated waveform models for the analysis. We present a state-of-the-art effective-one-body (EOB) model for the multipolar merger-ringdown waveform from dynamical capture black-hole mergers with arbitrary mass-ratio and nonprecessing spins. The model relies on analytical descriptions of the radiation reaction and waveform along generic orbits that are obtained by incorporating generic Newtonian prefactors in the expressions used in the quasi-circular case. It provides a tool for generating waveforms for generic binary black hole coalescences and for GW data analysis. We demonstrate that the model reliably accounts for the rich phenomenology of dynamical captures, from direct plunge to successive close encounters up to merger. The parameter space is fully characterized in terms of the initial energy and angular momentum. Our model reproduces to few percent the scattering angle from ten equal-mass, nonspinning, hyperbolic encounter numerical-relativity (NR) simulations. The agreement can be further improved to the by incorporating 6PN-results in one of the EOB potentials and tuning currently unknown analytical parameters. Our results suggest that NR simulations of hyperbolic encounters (and dynamical captures) can be used to inform EOB waveform models for generic BBH mergers/encounters for present and future GW detectors.

Testing the multipole structure and conservative dynamics of compact binaries using gravitational wave observations: The spinning case

In an earlier work [S. Kastha et al., Phys. Rev. D 98, 124033 (2018)], we developed the parametrized multipolar gravitational wave phasing formula to test general relativity, for the nonspinning compact binaries in quasicircular orbit. In this paper, we extend the method and include the important effect of spins in the inspiral dynamics. Furthermore, we consider parametric scaling of post-Newtonian (PN) coefficients of the conserved energy for the compact binary, resulting in the parametrized phasing formula for nonprecessing spinning compact binaries in quasicircular orbit. We also compute the projected accuracies with which the second and third generation ground-based gravitational wave detector networks as well as the planned space-based detector LISA will be able to measure the multipole deformation parameters and the binding energy parameters. Based on different source configurations, we find that a network of third-generation detectors would have comparable ability to that of LISA in constraining the conservative and dissipative dynamics of the compact binary systems. This parametrized multipolar waveform would be extremely useful not only in deriving the first upper limits on any deviations of the multipole and the binding energy coefficients from general relativity using the gravitational wave detections, but also for science case studies of next generation gravitational wave detectors.

Time-domain effective-one-body gravitational waveforms for coalescing compact binaries with nonprecessing spins, tides, and self-spin effects

We present TEOBResumS, a new effective-one-body (EOB) waveform model for nonprecessing (spin-aligned) and tidally interacting compact binaries.Spin-orbit and spin-spin effects are blended together by making use of the concept of centrifugal EOB radius. The point-mass sector through merger and ringdown is informed by numerical relativity (NR) simulations of binary black holes (BBH) computed with the SpEC and BAM codes. An improved, NR-based phenomenological description of the postmerger waveform is developed.The tidal sector of TEOBResumS describes the dynamics of neutron star binaries up to merger and incorporates a resummed attractive potential motivated by recent advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. Equation-of-state dependent self-spin interactions (monopole-quadrupole effects) are incorporated in the model using leading-order post-Newtonian results in a new expression of the centrifugal radius. TEOBResumS is compared to 135 SpEC and 19 BAM BBH waveforms. The maximum unfaithfulness to SpEC data $\bar{F}$ -- at design Advanced-LIGO sensitivity and evaluated with total mass $M$ varying between $10M_\odot \leq M \leq 200 M_\odot$ --is always below $2.5 \times 10^{-3}$ except for a single outlier that grazes the $7.1 \times 10^{-3}$ level. When compared to BAM data, $\bar{F}$ is smaller than $0.01$ except for a single outlier in one of the corners of the NR-covered parameter space, that reaches the $0.052$ level.TEOBResumS is also compatible, up to merger, to high end NR waveforms from binary neutron stars with spin effects and reduced initial eccentricity computed with the BAM and THC codes. The model is designed to generate accurate templates for the analysis of LIGO-Virgo data through merger and ringdown. We demonstrate its use by analyzing the publicly available data for GW150914.

Hangup effect in unequal mass binary black hole mergers and further studies of their gravitational radiation and remnant properties

We present the results of 74 new simulations of nonprecessing spinning black hole binaries with mass ratios $q=m_1/m_2$ in the range $1/7\leq q\leq1$ and individual spins covering the parameter space $-0.95\leq\alpha_{1,2}\leq0.95$ with one runs with spins of $\pm0.95$. We supplement those runs with 107 previous simulations to study the hangup effect in black hole mergers, i.e. the delay or prompt merger of spinning holes with respect to non spinning binaries. We perform the numerical evolution for typically the last ten orbits before the merger and down to the formation of the final remnant black hole. This allows us to study the hangup effect for unequal mass binaries leading us to identify the spin variable that controls the number of orbits before merger as $\vec{S}_{hu}\cdot{\hat{L}},$ where $\vec{S}_{hu}=(1+\frac12\frac{m_2}{m_1})\vec{S}_1+(1+\frac12\frac{m_1}{m_2})\vec{S}_2$. We also combine the total results of those 181 simulations to obtain improved fitting formulae for the remnant final black hole mass, spin and recoil velocity as well as for the peak luminosity and peak frequency of the gravitational strain, and find new correlations among them. This accurate new set of simulations enhances the number of available numerical relativity waveforms available for parameter estimation of gravitational wave observations.

Effects of nonquadrupole modes in the detection and parameter estimation of black hole binaries with nonprecessing spins

We study the effect of nonquadrupolar modes in the detection and parameter estimation of gravitational waves (GWs) from black hole binaries with nonprecessing spins, using Advanced LIGO. We evaluate the loss of the signal-to-noise ratio (SNR) and the systematic errors in the estimated parameters when a quadrupole-mode template family is used to detect GW signals with all the relevant modes. Target signals including nonquadrupole modes are constructed by matching numerical-relativity simulations of nonprecessing black hole binaries describing the late inspiral, merger and ringdown with post-Newtonian/effective-one-body waveforms describing the early inspiral. We find that neglecting nonquadrupole modes will, in general, cause unacceptable loss in the detection rate and unacceptably large systematic errors in the estimated parameters, for the case of massive binaries with large mass ratios. For a given mass ratio, neglecting subdominant modes will result in a larger loss in the detection rate for binaries with aligned spins. For binaries with antialigned spins, quadrupole-mode templates are more effectual in detection, at the cost of introducing a larger systematic bias in the parameter estimation. We provide a summary of the regions in the parameter space where neglecting nonquadrupole modes will cause an unacceptable loss of detection rates and unacceptably large systematic biases in the estimated parameters.

Remnant of binary black-hole mergers: New simulations and peak luminosity studies

We present the results of 61 new simulations of nonprecessing spinning black hole binaries with mass ratios $q={m}_{1}/{m}_{2}$ in the range $1/3\ensuremath{\le}q\ensuremath{\le}1$ and individual spins covering the parameter space $\ensuremath{-}0.85\ensuremath{\le}{\ensuremath{\alpha}}_{1,2}\ensuremath{\le}0.85$. We additionally perform ten new simulations of nonspinning black hole binaries with mass ratios covering the range $1/6\ensuremath{\le}q&lt;1$. We follow the evolution for typically the last ten orbits before merger down to the formation of the final remnant black hole. This allows for assessment of the accuracy of our previous empirical formulas for relating the binary parameters to the remnant final black hole mass, spin and recoil. We use the new simulation to improve the fit to the above remnant formulas and add a formula for the peak luminosity of gravitational waves, produced around the merger of the two horizons into one. We find excellent agreement (typical errors $\ensuremath{\sim}0.1%--0.2%$) for the mass and spin, and within $\ensuremath{\sim}5%$ for the recoil and peak luminosity. These formulas have direct application to parameter estimation techniques applied to LIGO observations of gravitational waves from binary black hole mergers.

Ready-to-use post-Newtonian gravitational waveforms for binary black holes with nonprecessing spins: An update

For black-hole binaries whose spins are (anti-) aligned with respect to the orbital angular momentum of the binary, we compute the frequency domain phasing coefficients including the quadratic-in-spin terms up to the third post-Newtonian (3PN) order, the cubic-in-spin terms at the leading order, 3.5PN, and the spin-orbit effects up to the 4PN order. In addition, we obtain the 2PN spin contributions to the amplitude of the frequency-domain gravitational waveforms for non-precessing binaries, using recently derived expressions for the time-domain polarization amplitudes of binaries with generic spins, complete at that accuracy level. These two results are updates to Arun et al. (2009) [1] for amplitude and Wade et al. (2013) [2] for phasing. They should be useful to construct banks of templates that model accurately non-precessing inspiraling binaries, for parameter estimation studies, and or constructing analytical template families that accounts for the inspiral-merger-ringdown phases of the binary.

Energetics and phasing of nonprecessing spinning coalescing black hole binaries

We present an improved numerical relativity (NR) calibration of the new effective-one-body (EOB) model for coalescing non precessing spinning black hole binaries recently introduced by Damour and Nagar [Physical Review D 90, 044018 (2014)]. We do so by comparing the EOB predictions to both the phasing and the energetics provided by two independent sets of NR data covering mass ratios $1\leq q \leq 9.989$ and dimensionless spin range $-0.95\leq \chi\leq +0.994$. One set of data is a subset of the Simulating eXtreme Spacetimes (SXS) catalog of public waveforms; the other set consists of new simulations obtained with the Llama code plus Cauchy Characteristic Evolution. We present the first systematic computation of the gauge-invariant relation between the binding energy and the total angular momentum, $E_{b}(j)$, for a large sample of, spin-aligned, SXS and Llama data. The dynamics of the EOB model presented here involves only two free functional parameters, one ($a_6^c(\nu)$) entering the non spinning sector, as a 5PN effective correction to the interaction potential, and one ($c_3(\tilde{a}_1,\tilde{a}_2,\nu))$ in the spinning sector, as an effective next-to-next-to-next-to-leading order correction to the spin-orbit coupling. These parameters are determined (together with a third functional parameter $\Delta t_{\rm NQC}(\chi)$ entering the waveform) by comparing the EOB phasing with the SXS phasing, the consistency of the energetics being checked afterwards. The quality of the analytical model for gravitational wave data analysis purposes is assessed by computing the EOB/NR faithfulness. Over the NR data sample and when varying the total mass between 20 and 200~$M_\odot$ the EOB/NR unfaithfulness (integrated over the NR frequency range) is found to vary between $99.493\%$ and $99.984\%$ with a median value of $99.944\%$.

Testing the validity of the phenomenological gravitational waveform models for nonspinning binary black hole searches at low masses

The phenomenological gravitational waveform models, which we refer to as PhenomA, PhenomB, and PhenomC, generate full inspiral, merger, and ringdown (IMR) waveforms of coalescing binary back holes (BBHs). These models are defined in the Fourier domain, thus can be used for fast matched filtering in the gravitational wave search. PhenomA has been developed for nonspinning BBH waveforms, while PhenomB and PhenomC were designed to model the waveforms of BBH systems with nonprecessing (aligned) spins, but can also be used for nonspinning systems. In this work, we study the validity of the phenomenological models for nonspinning BBH searches at low masses, and with Advanced LIGO. As our complete signal waveform model, we adopt EOBNRv2, which is a time-domain IMR waveform model. To investigate the search efficiency of the phenomenological template models, we calculate fitting factors (FFs) by exploring overlap surfaces. We find that only PhenomC is valid to obtain FFs better than 0.97 in the mass range of Above PhenomA is most efficient in symmetric mass region, PhenomB is most efficient in highly asymmetric mass region, and PhenomC is most efficient in the intermediate region. Specifically, we propose an effective phenomenological template family that can be constructed by employing the phenomenological models in four subregions individually. We find that FFs of the effective templates are better than 0.97 in our entire mass region and mostly greater than 0.99.

New effective-one-body description of coalescing nonprecessing spinning black-hole binaries

We present a new, tunable effective-one-body (EOB) model of the motion and radiation of coalescing black hole binaries with arbitrary mass ratio and aligned spins. The most novel feature of our formalism is the introduction, and systematic use, of the (gauge-invariant) concept of {\it centrifugal radius} $r_{c}$. In the spinning small mass-ratio limit, the main radial potential expressed in terms of $r_{c}$ differs very little (and only multiplicatively so) from the usual Schwarzschild potential $1-2M/r_{c}$. This motivates a new, multiplicative way of blending finite-mass-ratio deformations with spin-deformations. In the present exploratory work we consider a minimal version of our spinning EOB model (containing essentially only two adjustable parameters: one in the Hamiltonian and one in the waveform) and calibrate its (dominant mode) waveform against a sample of fifteen equal-mass, equal-spin waveforms produced by the SXS collaboration, and covering the dimensionless spin range $-0.95\leq \chi \leq +0.98$. The numerical relativity / EOB phasing disagreement remains remarkably small ($\lesssim \pm 0.15$ rad) over the entire spin range.

Optimizing gravitational-wave searches for a population of coalescing binaries: Intrinsic parameters

We revisit the problem of searching for gravitational waves from inspiralling compact binaries in Gaussian coloured noise. For binaries with quasicircular orbits and non-precessing component spins, considering dominant mode emission only, if the intrinsic parameters of the binary are known then the optimal statistic for a single detector is the well-known two-phase matched filter. However, the matched filter signal-to-noise ratio is /not/ in general an optimal statistic for an astrophysical population of signals, since their distribution over the intrinsic parameters will almost certainly not mirror that of noise events, which is determined by the (Fisher) information metric. Instead, the optimal statistic for a given astrophysical distribution will be the Bayes factor, which we approximate using the output of a standard template matched filter search. We then quantify the possible improvement in number of signals detected for various populations of non-spinning binaries: for a distribution of signals uniformly distributed in volume and with component masses distributed uniformly over the range $1\leq m_{1,2}/M_\odot\leq 24$, $(m_1+m_2) /M_\odot\leq 25$ at fixed expected SNR, we find $\gtrsim 20\%$ more signals at a false alarm threshold of $10^{-6}\,$Hz in a single detector. The method may easily be generalized to binaries with non-precessing spins.

Improving the sensitivity of a search for coalescing binary black holes with nonprecessing spins in gravitational wave data

We demonstrate for the first time a search pipeline with improved sensitivity to gravitational waves from coalescing binary black holes with spins aligned to the orbital angular momentum by the inclusion of spin effects in the search templates. We study the pipeline recovery of simulated gravitational wave signals from aligned-spin binary black holes added to real detector noise, comparing the pipeline performance with aligned-spin filter templates to the same pipeline with nonspinning filter templates. Our results exploit a three-parameter phenomenological waveform family that models the full inspiral-merger-ringdown coalescence and treats the effect of aligned spins with a single effective spin parameter $\ensuremath{\chi}$. We construct template banks from these waveforms by a stochastic placement method and use these banks as filters in the recently developed gstlal search pipeline. We measure the observable volume of the analysis pipeline for binary black hole signals with ${M}_{\text{total}}$ and $\ensuremath{\chi}\ensuremath{\in}[0,0.85]$. We find an increase in observable volume of up to 45% for systems with $0.2\ensuremath{\le}\ensuremath{\chi}\ensuremath{\le}0.85$ with almost no loss of sensitivity to signals with $0\ensuremath{\le}\ensuremath{\chi}\ensuremath{\le}0.2$. We also show that the use of spinning templates in the search pipeline provides for more accurate recovery of the binary mass parameters as well as an estimate of the effective spin parameter. We demonstrate this analysis on 25.9 days of data obtained from the Hanford and Livingston detectors in LIGO's fifth observation run.

Testing the validity of the single-spin approximation in inspiral-merger-ringdown waveforms

Gravitational-wave signals from black-hole binaries with non-precessing spins are described by four parameters -- each black hole's mass and spin. It has been shown that the dominant spin effects can be modeled by a \emph{single} spin parameter, leading to the development of several \emph{three-parameter} waveform models. Previous studies indicate that these models should be adequate for gravitational-wave detection. In this paper we focus on the systematic biases that would result from using them to estimate binary parameters, and consider a one-parameter family of configurations at one choice of mass ratio and effective single spin. We find that for low-mass binaries within that family of configurations, where the observable waveform is dominated by the inspiral, the systematic bias in all physical parameters is smaller than the parameter uncertainty due to degeneracies between the mass ratio and the spins, at least up to signal-to-noise ratios (SNRs) of 50. For higher-mass binaries, where the merger and ringdown make a greater contribution to the observed signal, the bias in the mass ratio is comparable to its uncertainty at SNRs of only $\sim$30, and the bias in the measurement of the total spin is \emph{larger} than the uncertainty defined by the 90% confidence region even at an SNR of only 10. Although this bias may be mitigated in future models by a better choice of single-effective-spin parameter, these results suggest that it may be possible to accurately measure \emph{both} black-hole spins in intermediate-mass binaries.

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

We first use five non-spinning and two mildly spinning (chi_i \simeq -0.44, +0.44) numerical-relativity waveforms of black-hole binaries and calibrate an effective-one-body (EOB) model for non-precessing spinning binaries, notably its dynamics and the dominant (2,2) gravitational-wave mode. Then, we combine the above results with recent outcomes of small-mass-ratio simulations produced by the Teukolsky equation and build a prototype EOB model for detection purposes, which is capable of generating inspiral-merger-ringdown waveforms for non-precessing spinning black-hole binaries with any mass ratio and individual black-hole spins -1 \leq chi_i \lesssim 0.7. We compare the prototype EOB model to two equal-mass highly spinning numerical-relativity waveforms of black holes with spins chi_i = -0.95, +0.97, which were not available at the time the EOB model was calibrated. In the case of Advanced LIGO we find that the mismatch between prototype-EOB and numerical-relativity waveforms is always smaller than 0.003 for total mass 20-200 M_\odot, the mismatch being computed by maximizing only over the initial phase and time. To successfully generate merger waveforms for individual black-hole spins chi_i \gtrsim 0.7, the prototype-EOB model needs to be improved by (i) better modeling the plunge dynamics and (ii) including higher-order PN spin terms in the gravitational-wave modes and radiation-reaction force.

Addressing the spin question in gravitational-wave searches: Waveform templates for inspiralling compact binaries with nonprecessing spins

This paper presents a post-Newtonian (PN) template family of gravitational waveforms from inspiralling compact binaries with nonprecessing spins, where the spin effects are described by a single ``reduced-spin'' parameter. This template family, which reparametrizes all the spin-dependent PN terms in terms of the leading-order (1.5PN) spin-orbit coupling term in an approximate way, has very high overlaps (fitting factor $&gt;0.99$) with nonprecessing binaries with arbitrary mass ratios and spins. We also show that this template family is ``effectual'' for the detection of a significant fraction of generic spinning binaries in the comparable-mass regime (${m}_{2}/{m}_{1}\ensuremath{\lesssim}10$), providing an attractive and feasible way of searching for gravitational waves from spinning low-mass binaries. We also show that the secular (nonoscillatory) spin-dependent effects in the phase evolution (which are taken into account by the nonprecessing templates) are more important than the oscillatory effects of precession in the comparable-mass (${m}_{1}\ensuremath{\simeq}{m}_{2}$) regime. Hence the effectualness of nonspinning templates is particularly poor in this case, as compared to non-precessing-spin templates. For the case of binary neutron stars observable by Advanced LIGO, even moderate spins (${\stackrel{^}{\mathbf{L}}}_{N}\ifmmode\cdot\else\textperiodcentered\fi{}\mathbf{S}/{m}^{2}\ensuremath{\simeq}0.015--0.1$) will cause considerable mismatches ($\ensuremath{\sim}3%--25%$) with nonspinning templates. This is contrary to the expectation that neutron-star spins may not be relevant for gravitational wave detection.

Inspiral-Merger-Ringdown Waveforms for Black-Hole Binaries with Nonprecessing Spins

We present the first analytical inspiral-merger-ringdown gravitational waveforms from binary black holes (BBHs) with nonprecessing spins, that is based on a description of the late-inspiral, merger and ringdown in full general relativity. By matching a post-Newtonian description of the inspiral to a set of numerical-relativity simulations, we obtain a waveform family with a conveniently small number of physical parameters. These waveforms will allow us to detect a larger parameter space of BBH coalescence, including a considerable fraction of precessing binaries in the comparable-mass regime, thus significantly improving the expected detection rates.