Template Banks for Binary black hole searches with Numerical Relativity waveforms

Gravitational waves (GW) from coalescing stellar-mass black hole binaries (BBH) are expected to be detected by the Advanced Laser Interferometer Gravitational-wave Observatory and Advanced Virgo. Detection searches operate by matched-filtering the detector data using a bank of waveform templates. Traditionally, template banks for BBH are constructed from intermediary analytical waveform models which are calibrated against numerical relativity simulations and which can be aluated for any choice of BBH parameters. This paper explores an alternative to the traditional approach, namely the construction of template banks directly from numerical BBH simulations. Using non-spinning BBH systems as an example, we demonstrate which regions of the mass-parameter plane can be covered with existing numerical BBH waveforms. We estimate the required number and required length of BBH simulations to cover the entire non-spinning BBH parameter plane up to mass-ratio 10, thus illustrating that our approach can be used to guide parameter placement of future numerical simulations. We derive error bounds which are independent of analytical waveform models; therefore, our formalism can be used to independently test the accuracy of such waveform models. The resulting template banks are suitable for advanced LIGO searches.

Investigating the effect of precession on searches for neutron-star-black-hole binaries with Advanced LIGO

The first direct detection of neutron-star-black-hole binaries will likely be made with gravitational-wave observatories. Advanced LIGO and Advanced Virgo will be able to observe neutron-star-black-hole mergers at a maximum distance of 900Mpc. To acheive this sensitivity, gravitational-wave searches will rely on using a bank of filter waveforms that accurately model the expected gravitational-wave signal. The angular momentum of the black hole is expected to be comparable to the orbital angular momentum. This angular momentum will affect the dynamics of the inspiralling system and alter the phase evolution of the emitted gravitational-wave signal. In addition, if the black holes angular momentum is not aligned with the orbital angular momentum it will cause the orbital plane of the system to precess. In this work we demonstrate that if the effect of the black holes angular momentum is neglected in the waveform models used in gravitational-wave searches, the detection rate of $(10+1.4)M_{odot}$ neutron-star--black-hole systems would be reduced by $33 - 37%$. The error in this measurement is due to uncertainty in the Post-Newtonian approximations that are used to model the gravitational-wave signal of neutron-star-black-hole inspiralling binaries. We describe a new method for creating a bank of filter waveforms where the black hole has non-zero angular momentum, but is aligned with the orbital angular momentum. With this bank we find that the detection rate of $(10+1.4)M_{odot}$ neutron-star-black-hole systems would be reduced by $26-33%$. Systems that will not be detected are ones where the precession of the orbital plane causes the gravitational-wave signal to match poorly with non-precessing filter waveforms. We identify the regions of parameter space where such systems occur and suggest methods for searching for highly precessing neutron-star-black-hole binaries.

Accuracy of gravitational waveform models for observing neutron-star--black-hole binaries in Advanced LIGO

Gravitational waves radiated by the coalescence of compact-object binaries containing a neutron star and a black hole are one of the most interesting sources for the ground-based gravitational-wave observatories Advanced LIGO and Advanced Virgo. Advanced LIGO will be sensitive to the inspiral of a $1.4, M_odot$ neutron star into a $10,M_odot$ black hole to a maximum distance of $sim 900$ Mpc. Achieving this sensitivity and extracting the physics imprinted in observed signals requires accurate modeling of the binary to construct template waveforms. In a NSBH binary, the black hole may have significant angular momentum (spin), which affects the phase evolution of the emitted gravitational waves. We investigate the ability of post-Newtonian (PN) templates to model the gravitational waves emitted during the inspiral phase of NSBH binaries. We restrict the black holes spin to be aligned with the orbital angular momentum and compare several approximants. We examine restricted amplitude waveforms that are accurate to 3.5PN order in the orbital dynamics and complete to 2.5PN order in the spin dynamics. We also consider PN waveforms with the recently derived 3.5PN spin-orbit and 3PN spin-orbit tail corrections. We compare these approximants to the effective-one-body model. For all these models, large disagreements start at low to moderate black hole spins, particularly for binaries where the spin is anti-aligned with the orbital angular momentum. We show that this divergence begins in the early inspiral at $v sim 0.2$ for $chi_{BH} sim 0.4$. PN spin corrections beyond those currently known will be required for optimal detection searches and to measure the parameters of neutron star--black hole binaries. While this complicates searches, the strong dependence of the gravitational-wave signal on the spin dynamics will make it possible to extract significant astrophysical information.

Detecting binary neutron star systems with spin in advanced gravitational-wave detectors

The detection of gravitational waves from binary neutron stars is a major goal of the gravitational-wave observatories Advanced LIGO and Advanced Virgo. Previous searches for binary neutron stars with LIGO and Virgo neglected the component stars angular momentum (spin). We demonstrate that neglecting spin in matched-filter searches causes advanced detectors to lose more than 3% of the possible signal-to-noise ratio for 59% (6%) of sources, assuming that neutron star dimensionless spins, $cmathbf{J}/GM^2$, are uniformly distributed with magnitudes between 0 and 0.4 (0.05) and that the neutron stars have isotropically distributed spin orientations. We present a new method for constructing template banks for gravitational wave searches for systems with spin. We present a new metric in a parameter space in which the template placement metric is globally flat. This new method can create template banks of signals with non-zero spins that are (anti-)aligned with the orbital angular momentum. We show that this search loses more than 3% of the maximium signal-to-noise for only 9% (0.2%) of BNS sources with dimensionless spins between 0 and 0.4 (0.05) and isotropic spin orientations. Use of this template bank will prevent selection bias in gravitational-wave searches and allow a more accurate exploration of the distribution of spins in binary neutron stars.

Nonspinning searches for spinning binaries in ground-based detector data: Amplitude and mismatch predictions in the constant precession cone approximation

Current searches for compact binary mergers by ground-based gravitational-wave detectors assume for simplicity the two bodies are not spinning. If the binary contains compact objects with significant spin, then this can reduce the sensitivity of these searches, particularly for black hole--neutron star binaries. In this paper we investigate the effect of neglecting precession on the sensitivity of searches for spinning binaries using non-spinning waveform models. We demonstrate that in the sensitive band of Advanced LIGO, the angle between the binarys orbital angular momentum and its total angular momentum is approximately constant. Under this emph{constant precession cone} approximation, we show that the gravitational-wave phasing is modulated in two ways: a secular increase of the gravitational-wave phase due to precession and an oscillation around this secular increase. We show that this secular evolution occurs in precisely three ways, corresponding to physically different apparent evolutions of the binarys precession about the line of sight. We estimate the best possible fitting factor between emph{any} non-precessing template model and a single precessing signal, in the limit of a constant precession cone. Our closed form estimate of the fitting-factor depends only the geometry of the in-band precession cone; it does not depend explicitly on binary parameters, detector response, or details of either signal model. The precessing black hole--neutron star waveforms least accurately matched by nonspinning waveforms correspond to viewing geometries where the precession cone sweeps the orbital plane repeatedly across the line of sight, in an unfavorable polarization alignment.

The Effect of Eccentricity on Searches for Gravitational-Waves from Coalescing Compact Binaries in Ground-based Detectors

Inspiralling compact binaries are expected to circularize before their gravitational-wave signals reach the sensitive frequency band of ground-based detectors. Current searches for gravitational waves from compact binaries using the LIGO and Virgo detectors therefore use circular templates to construct matched filters. Binary formation models have been proposed which suggest that some systems detectable by the LIGO--Virgo network may have non-negligible eccentricity. We investigate the ability of the restricted 3.5 post-Newtonian order TaylorF2 template bank, used by LIGO and Virgo to search for gravitational waves from compact binaries with masses $M le 35 M_odot$, to detect binaries with non-zero eccentricity. We model the gravitational waves from eccentric binaries using the $x$-model post-Newtonian formalism proposed by Hinder emph{et. al.} [I. Hinder, F. Hermann, P. Laguna, and D. Shoemaker, arXiv:0806.1037v1]. We find that small residual eccentricities ($e_0 lesssim 0.05$ at 40 Hz) do not significantly affect the ability of current LIGO searches to detect gravitational waves from coalescing compact binaries with total mass $2 M_odot < M < 15 M_odot$. For eccentricities $e_0 gtrsim 0.1$, the loss in matched filter signal-to-noise ratio due to eccentricity can be significant and so templates which include eccentric effects will be required to perform optimal searches for such systems.

Comparison of high-accuracy numerical simulations of black-hole binaries with stationary phase post-Newtonian template waveforms for Initial and Advanced LIGO

We study the effectiveness of stationary-phase approximated post-Newtonian waveforms currently used by ground-based gravitational-wave detectors to search for the coalescence of binary black holes by comparing them to an accurate waveform obtained from numerical simulation of an equal-mass non-spinning binary black hole inspiral, merger and ringdown. We perform this study for the Initial- and Advanced-LIGO detectors. We find that overlaps between the templates and signal can be improved by integrating the match filter to higher frequencies than used currently. We propose simple analytic frequency cutoffs for both Initial and Advanced LIGO, which achieve nearly optimal matches, and can easily be extended to unequal-mass, spinning systems. We also find that templates that include terms in the phase evolution up to 3.5 pN order are nearly always better, and rarely significantly worse, than 2.0 pN templates currently in use. For Initial LIGO we recommend a strategy using templates that include a recently introduced pseudo-4.0 pN term in the low-mass ($M leq 35 MSun$) region, and 3.5 pN templates allowing unphysical values of the symmetric reduced mass $eta$ above this. This strategy always achieves overlaps within 0.3% of the optimum, for the data used here. For Advanced LIGO we recommend a strategy using 3.5 pN templates up to $M=12 MSun$, 2.0 pN templates up to $M=21 MSun$, pseudo-4.0 pN templates up to $65 MSun$, and 3.5 pN templates with unphysical $eta$ for higher masses. This strategy always achieves overlaps within 0.7% of the optimum for Advanced LIGO.

Rates and Characteristics of Intermediate Mass Ratio Inspirals Detectable by Advanced LIGO

Gravitational waves (GWs) from the inspiral of a neutron star (NS) or stellar-mass black hole (BH) into an intermediate-mass black hole (IMBH) with mass between ~50 and ~350 solar masses may be detectable by the planned advanced generation of ground-based GW interferometers. Such intermediate mass ratio inspirals (IMRIs) are most likely to be found in globular clusters. We analyze four possible IMRI formation mechanisms: (1) hardening of an NS-IMBH or BH-IMBH binary via three-body interactions, (2) hardening via Kozai resonance in a hierarchical triple system, (3) direct capture, and (4) inspiral of a compact object from a tidally captured main-sequence star; we also discuss tidal effects when the inspiraling object is an NS. For each mechanism we predict the typical eccentricities of the resulting IMRIs. We find that IMRIs will have largely circularized by the time they enter the sensitivity band of ground-based detectors. Hardening of a binary via three-body interactions, which is likely to be the dominant mechanism for IMRI formation, yields eccentricities under 10^-4 when the GW frequency reaches 10 Hz. Even among IMRIs formed via direct captures, which can have the highest eccentricities, around 90% will circularize to eccentricities under 0.1 before the GW frequency reaches 10 Hz. We estimate the rate of IMRI coalescences in globular clusters and the sensitivity of a network of three Advanced LIGO detectors to the resulting GWs. We show that this detector network may see up to tens of IMRIs per year, although rates of one to a few per year may be more plausible. We also estimate the loss in signal-to-noise ratio that will result from using circular IMRI templates for data analysis and find that, for the eccentricities we expect, this loss is negligible.

High-accuracy comparison of numerical relativity simulations with post-Newtonian expansions

Numerical simulations of 15 orbits of an equal-mass binary black hole system are presented. Gravitational waveforms from these simulations, covering more than 30 cycles and ending about 1.5 cycles before merger, are compared with those from quasi-circular zero-spin post-Newtonian (PN) formulae. The cumulative phase uncertainty of these comparisons is about 0.05 radians, dominated by effects arising from the small residual spins of the black holes and the small residual orbital eccentricity in the simulations. Matching numerical results to PN waveforms early in the run yields excellent agreement (within 0.05 radians) over the first $sim 15$ cycles, thus validating the numerical simulation and establishing a regime where PN theory is accurate. In the last 15 cycles to merger, however, {em generic} time-domain Taylor approximants build up phase differences of several radians. But, apparently by coincidence, one specific post-Newtonian approximant, TaylorT4 at 3.5PN order, agrees much better with the numerical simulations, with accumulated phase differences of less than 0.05 radians over the 30-cycle waveform. Gravitational-wave amplitude comparisons are also done between numerical simulations and post-Newtonian, and the agreement depends on the post-Newtonian order of the amplitude expansion: the amplitude difference is about 6--7% for zeroth order and becomes smaller for increasing order. A newly derived 3.0PN amplitude correction improves agreement significantly ($<1%$ amplitude difference throughout most of the run, increasing to 4% near merger) over the previously known 2.5PN amplitude terms.

Searching for Gravitational Waves from Binary Inspirals with LIGO

We describe the current status of the search for gravitational waves from inspiralling compact binary systems in LIGO data. We review the result from the first scientific run of LIGO (S1). We present the goals of the search of data taken in the second scientific run (S2) and describe the differences between the methods used in S1 and S2.