The inspiral and merger of black-hole binary systems are a promising source of gravitational waves. The most effective method to look for a signal with a well understood waveform, such as the binary black hole signal, is matched filtering against a l
ibrary of model waveforms. Current model waveforms are comprised solely of the dominant radiation mode, the quadrupole mode, although it is known that there can be significant power in the higher-order modes. The binary black hole waveforms produced by numerical relativity are accurate through late inspiral, merger, and ringdown and include the higher-order modes. The available numerical-relativity waveforms span an increasing portion of the physical parameter space of unequal mass, spin and precession. In this paper, we investigate the degree to which gravitational-wave searches could be improved by the inclusion of higher modes in the model waveforms, for signals with a variety of initial mass ratios and generic spins. Our investigation studies how well the quadrupole-only waveform model matches the signal as a function of the inclination and orientation of the source and how the modes contribute to the distance reach into the Universe of Advanced LIGO for a fixed set of internal source parameters. The mismatch between signals and quadrupole-only waveform can be large, dropping below 0.97 for up to 65% of the source-sky for the non-precessing cases we studied, and over a larger area in one precessing case. There is a corresponding 30% increase in detection volume that could be achieved by adding higher modes to the search; however, this is mitigated by the fact that the mismatch is largest for signals which radiate the least energy and to which the search is therefore least sensitive. Likewise, the mismatch is largest in directions from the source along which the least energy is radiated.
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 fr
om 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.
