LIGO-Virgo searches for gravitational waves from coalescing binaries: a status update
Coalescing compact binaries of neutron stars and/or black holes are considered as one of the most promising sources for Earth based gravitational wave detectors. The LIGO-Virgo joint collaboration’s Compact Binary Coalescence (CBC) group is searching for gravitational waves emitted by these astrophysical systems by matched filtering the data against theoretically modeled template waveforms. A variety of waveform template families are employed depending on the mass range probed by the search and the stage of the inspiral phase targeted: restricted post-Newtonian for systems having total mass less than , numerical relativity inspired complete inspiral-merger-ringdown waveforms for more massive systems up to and ringdown templates for modeling perturbed black holes up to . We give a status update on CBC group’s current efforts and upcoming plans in detecting signatures of astrophysical gravitational waves.
LIGO Laboratory, California Institute of Technology, Pasadena CA 91125, USA.
Gravitational waves emitted by compact binaries of neutron stars and/or black holes in the last few minutes before merger lie within the sensitive bandwidth of present ground-based interferometric detectors like LIGO  and Virgo  and are considered likely candidates for first detection in these instruments. This is partly motivated by the fact that the phase and amplitude evolution of waveforms from such sources are well modeled analytically in the inspiral phase and ringdown phase when the excited final black hole emits gravitational waves and loses energy [3, 4]. This allows the construction of templates for use in the application of matched-filtering techniques  to filter the noisy data. The intermediate merger phase is not well modeled but can be obtained by numerical relativity . These have, in turn, inspired a whole generation of template waveforms which bridge the inspiral, merger and ringdown parts by various approximations [7, 8].
The three LIGO detectors at Livingston and Hanford operated at their design sensitivity in a continuous data acquisition mode between November 2005 and September 2007 in their fifth science run (S5), collecting a full year of triple coincident data. The first run of the Virgo detectors (VSR1) coincided with the last 6 months of S5 at a sensitivity comparable to the LIGO instruments above 400 Hz. The aim of this article is to provide a snap-shot of the current activities of the CBC group and present our plans for the near future.
2 Data processing
At the core of the CBC group’s data analysis pipeline  is the matched filtering engine that correlates data from detectors against theoretical template waveforms expected from the sources. A discrete lattice of templates spanning the intrinsic parameter space of the binary is constructed such that the maximum loss in signal-to-noise ratio (SNR) is no more than a prescribed limit. Filtering the data against a template yields the SNR given by
where denotes the noise-weighted scalar product of the data and the template. For each template, triggers that have SNR greater than a prescribed threshold are retained. The templates are typically parametrized by the component masses of the coalescing binary (or, for ringdowns, the mass and spin of the final blackhole). These triggers are tagged by the template parameters and coalescence time, which can be later used for clustering  and determining coincidence . Although the trigger generation by matched filtering is at the heart of the CBC pipelines, other important steps including: determining detection efficiency and accidental background estimation; vetoing glitchy times; determining coincidence across different instruments and following up interesting coincident triggers; all form important parts of the full data analysis pipeline and are used both for GW-triggered and externally triggered (eg, GRB events) searches.
3 Status of current activities of the CBC group
The template families are chosen to optimize the sensitivity of the searches to probe the part of the signal parameter space accessible by the ground-based detectors. The current activities of the group include: the analysis of coincident LIGO S5 and Virgo VSR1 data for low mass coalescing binaries (systems with total mass between and ); analysis of the S5 data for high-mass binaries (systems with total mass ) using inspiral-merger-ringdown (IMR) waveforms; those for ringdown of perturbed black holes using exponentially damped sinusoidal waveforms; and analysis of S5 data around times triggered by short, hard gamma-ray bursts.
The ‘S5 low-mass search’ covers the region of the total-mass parameter space employing frequency domain restricted post-Newtonian (PN) template waveforms accurate to 2PN order in phase evolution up to the innermost stable circular orbit (ISCO). Data is cross-correlated against such template banks consisting of 5-7 thousand templates to provide adequate coverage via semi-automated data analysis pipeline. Historically, it was divided into three separate searches to analyze different subsets of the S5 data: (a) The S5 first calendar year search  - sensitive to CBC as far as for systems with total mass at an SNR of 8. No gravitational waves were detected by this search. Improved upper limits on the rate of coalescence for neutron star and black hole binary systems were established. (b) The S5 12-18 month analysis covering 186 days beyond the first calendar year . No gravitational wave candidate events were reported in this search which yielded a -confidence rate upper limit of CBC of neutron stars at per year per where the latter is times the blue solar luminosity. (c) The last six months of S5 data was analyzed in coincidence with VSR1 allowing the first joint LIGO-Virgo data analysis. This work is nearing completion and publication. The upper limits are cumulative (i.e, they include the results from all previous analyses) and are combined in a Bayesian manner.
-  B. P. Abbott et al., Reports on Progress in Physics 72, 076901 (25pp) (2009).
-  F. Acernese et al., Class. Quant. Grav. 23, S635 (2006).
-  B. P. Abbott et al., Physical Review D 80, 062001 (2009).
-  A. Buonanno, B. R. Iyer, E. Ochsner, Y. Pan, and B. S. Sathyaprakash, Physical Review D (Particles, Fields, Gravitation, and Cosmology) 80, 084043 (2009).
-  B. S. Sathyaprakash and S. V. Dhurandhar, Phys. Rev D 44, 3819 (1991).
-  F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005).
-  P. Ajith, M. Hannam, S. Husa, Y. Chen, B. Bruegmann, N. Dorband, D. Mueller, F. Ohme, D. Pollney, C. Reisswig, L. Santamaria, and J. Seiler, ArXiv e-prints (2009).
-  A. Buonanno, Y. Pan, J. G. Baker, J. Centrella, B. J. Kelly, S. T. McWilliams, and J. R. van Meter, Physical Review D 76, 104049 (2007).
-  B. Allen, W. G. Anderson, P. R. Brady, D. A. Brown, and J. D. E. Creighton, FINDCHIRP: an algorithm for detection of gravitational waves from inspiraling compact binaries, 2005.
-  A. S. Sengupta, J. A. Gupchup, and C. A. K. Robinson, (2008), in Preparation.
-  C. A. K. Robinson, B. S. Sathyaprakash, and A. S. Sengupta, Physical Review D 78, 062002 (2008).
-  B. P. Abbott et al., Phys. Rev. D79, 122001 (2009).
-  B. P. Abbott et al., Phys. Rev. D80, 047101 (2009).
-  LIGO Scientific Collaboration and K. Hurley, The Astrophysical Journal 681, 1419 (2008).
-  M. van der Sluys, I. Mandel, V. Raymond, V. Kalogera, C. Roever, and N. Christensen, CLASS.QUANT.GRAV 204010, (2009).
-  P. Ajith, S. Babak, Y. Chen, M. Hewitson, B. Krishnan, A. M. Sintes, J. T. Whelan, B. Brügmann, P. Diener, N. Dorband, J. Gonzalez, M. Hannam, S. Husa, D. Pollney, L. Rezzolla, L. Santamaría, U. Sperhake, and J. Thornburg, Physical Review D 77, 104017 (2008).
-  T. Miller, A. S. Sengupta, K. Hodge, and A. Weinstein, LIGO Technical Document (LIGO-T0900454) from 2009 SURF program (2009).
-  J. Veitch and A. Vecchio, Classical and Quantum Gravity 25, 184010 (2008).
Apart from these searches, the CBC group is also developing techniques to follow-up interesting candidate events with a view to developing detection confidence. These include using Markov-chain Monte-Carlo methods  to estimate the parameters of detection candidates. Other activities include: developing better IMR waveforms that are faithful to the numerical waveforms over a wider range of parameters; including spin as a parameter in our detection templates - especially for high-mass end of the parameter space; detector characterization to understand the source of coincident glitches; and preparing the ground for upcoming, more sensitive science runs.
4 New data analysis techniques on the horizon
Since the beginning of S5, the basic data analysis pipeline has seen several improvements to its core functionality and many new features have been added. In this section, we highlight some of the most important ones. A key area of improvement has been the addition of numerical relativity inspired waveforms to the LSC and Virgo collaboration’s arsenal of existing templates: time-domain EOB waveforms calculated up to the light ring, calibrated to numerical simulations and then matched to the quasi-normal modes of the ringdown at the end, have extended the effective bandwidth of the high-mass search thereby increasing the distance reach and sensitivity of the high-mass search; see Figure 1. Phenomenological IMR waveforms [7, 16] that include spins (albeit restricted to non-precessing aligned/anti-aligned configurations) now allow us to not only probe significantly larger parts of the CBC parameter space but also to inject such signals into the data stream for more reliable estimates of the efficiency of our pipeline in detecting gravitational waves from CBCs. Several new analysis technologies are on the horizon: hardware-accelerated signal processing where the computer’s graphical processing units (GPU) are used to parallelize the fast-Fourier transform calculations have been tested and found to make our search pipeline an order of magnitude faster. Multi-variate statistical classification tools are being developed to bring in machine-learning algorithms that are able to classify our putative events and better distinguish them from background . New ranking statistics based on likelihood ratios have been found to perform better than traditional ‘effective SNR’ statistic and have been deployed in some of the searches. The coherent inspiral pipeline has matured to the point where it can be run in a follow-up mode for interesting candidates that pass the stringent mass and time coincidence tests. Efforts are also on towards a Bayesian analysis of data from arbitrary network of detectors for the extraction of parameters of a likely signal in the data . The full posterior probability density function for parameter estimation (using a wide range of template waveforms) can be computed, allowing estimation of likely parameters and the range of their uncertainty. These new techniques are at various stages of development and stem from a common desire of the group to work towards building confidence in our first detection.
5 Future plans: S6, VSR2 and beyond
Following S5 and VSR1 data taking, the detectors are undergoing scheduled upgrades which will improve the sensitivity in these instruments and the volume of space that these ‘enhanced’ detectors can probe. Data taking with the enhanced detectors (LIGO’s S6 and Virgo’s VSR2) began in June 2009 and will continue for approximately 1.5 years. At that point, the detectors will undergo more major upgrades. Advanced LIGO is scheduled to be fully commissioned in 2015. In these instruments, the improvement in sensitivity is expected to be a factor of 10 over initial LIGO and a 1000 fold increase in the event rates, making the detection of gravitational waves routine. S6/VSR2 data is being used as a test bed for several data analysis techniques developed by the CBC group in the recent past. A major paradigm shift has been to focus on rapid turn around of the results by analyzing a week’s data at a time. Low latency pipelines being developed by the group have been deployed and are being extensively tested. Many of the post-processing tasks and follow-up of interesting triggers from the main pipeline are being automated: among them, rapid sky localization and parameter estimation; and improved background estimation are being explored. The latter is motivated by a chance for electromagnetic followup by quick identification and localization of a transient corresponding to an interesting gravitational-wave candidates. This is also establishing the protocol for external collaborations with other ground and space based telescopes and the exciting opportunity to ‘see’ a gravitational wave trigger in different wavelengths.
I would like to thank Szabi Marka, Zzusa Marka and other organisers of the 8th E. Amaldi conference in New York City. I also thank Ajith Parameswaran, Craig Robinson and Alan Weinstein for reading and editing the manuscript, comments and feedback.