Inspiral-merger-ringdown waveforms for black-hole binaries with non-precessing spins

We present the first analytical inspiral-merger-ringdown gravitational waveforms from binary black holes (BBHs) with non-precessing 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.Comment: To appear in Phys. Rev. Lett. Significant new results. One figure removed due to page limitatio

Implementation of standard testbeds for numerical relativity

We discuss results that have been obtained from the implementation of the initial round of testbeds for numerical relativity which was proposed in the first paper of the Apples with Apples Alliance. We present benchmark results for various codes which provide templates for analyzing the testbeds and to draw conclusions about various features of the codes. This allows us to sharpen the initial test specifications, design a new test and add theoretical insight.Comment: Corrected versio

Radiation from low-momentum zoom-whirl orbits

We study zoom-whirl behaviour of equal mass, non-spinning black hole binaries in full general relativity. The magnitude of the linear momentum of the initial data is fixed to that of a quasi-circular orbit, and its direction is varied. We find a global maximum in radiated energy for a configuration which completes roughly one orbit. The radiated energy in this case exceeds the value of a quasi-circular binary with the same momentum by 15%. The direction parameter only requires minor tuning for the localization of the maximum. There is non-trivial dependence of the energy radiated on eccentricity (several local maxima and minima). Correlations with orbital dynamics shortly before merger are discussed. While being strongly gauge dependent, these findings are intuitive from a physical point of view and support basic ideas about the efficiency of gravitational radiation from a binary system.Comment: 9 pages, 6 figures, Amaldi8 conference proceedings as publishe

Method to estimate ISCO and ring-down frequencies in binary systems and consequences for gravitational wave data analysis

Recent advances in the description of compact binary systems have produced gravitational waveforms that include inspiral, merger and ring-down phases. Comparing results from numerical simulations with those of post-Newtonian (PN), and related, expansions has provided motivation for employing PN waveforms in near merger epochs when searching for gravitational waves and has encouraged the development of analytic fits to full numerical waveforms. The models and simulations do not yet cover the full binary coalescence parameter space. For these yet un-simulated regions, data analysts can still conduct separate inspiral, merger and ring-down searches. Improved knowledge about the end of the inspiral phase, the beginning of the merger, and the ring-down frequencies could increase the efficiency of both coherent inspiral-merger-ring-down (IMR) searches and searches over each phase separately. Insight can be gained for all three cases through a recently presented theoretical calculation, which, corroborated by the numerical results, provides an implicit formula for the final spin of the merged black holes, accurate to within 10% over a large parameter space. Knowledge of the final spin allows one to predict the end of the inspiral phase and the quasinormal mode ring-down frequencies, and in turn provides information about the bandwidth and duration of the merger. In this work we will discuss a few of the implications of this calculation for data analysis.Comment: Added references to section 3 14 pages 5 figures. Submitted to Classical and Quantum Gravit

Ninja data analysis with a detection pipeline based on the Hilbert-Huang Transform

The Ninja data analysis challenge allowed the study of the sensitivity of data analysis pipelines to binary black hole numerical relativity waveforms in simulated Gaussian noise at the design level of the LIGO observatory and the VIRGO observatory. We analyzed NINJA data with a pipeline based on the Hilbert Huang Transform, utilizing a detection stage and a characterization stage: detection is performed by triggering on excess instantaneous power, characterization is performed by displaying the kernel density enhanced (KD) time-frequency trace of the signal. Using the simulated data based on the two LIGO detectors, we were able to detect 77 signals out of 126 above SNR 5 in coincidence, with 43 missed events characterized by signal to noise ratio SNR less than 10. Characterization of the detected signals revealed the merger part of the waveform in high time and frequency resolution, free from time-frequency uncertainty. We estimated the timelag of the signals between the detectors based on the optimal overlap of the individual KD time-frequency maps, yielding estimates accurate within a fraction of a millisecond for half of the events. A coherent addition of the data sets according to the estimated timelag eventually was used in a characterization of the event.Comment: Accepted for publication in CQG, special issue NRDA proceedings 200

Matching post-Newtonian and numerical relativity waveforms: systematic errors and a new phenomenological model for non-precessing black hole binaries

We present a new phenomenological gravitational waveform model for the inspiral and coalescence of non-precessing spinning black hole binaries. Our approach is based on a frequency domain matching of post-Newtonian inspiral waveforms with numerical relativity based binary black hole coalescence waveforms. We quantify the various possible sources of systematic errors that arise in matching post-Newtonian and numerical relativity waveforms, and we use a matching criteria based on minimizing these errors; we find that the dominant source of errors are those in the post-Newtonian waveforms near the merger. An analytical formula for the dominant mode of the gravitational radiation of non-precessing black hole binaries is presented that captures the phenomenology of the hybrid waveforms. Its implementation in the current searches for gravitational waves should allow cross-checks of other inspiral-merger-ringdown waveform families and improve the reach of gravitational wave searches.Comment: 22 pages, 11 figure

Mode coupling in the nonlinear response of black holes

We study the properties of the outgoing gravitational wave produced when a non-spinning black hole is excited by an ingoing gravitational wave. Simulations using a numerical code for solving Einstein's equations allow the study to be extended from the linearized approximation, where the system is treated as a perturbed Schwarzschild black hole, to the fully nonlinear regime. Several nonlinear features are found which bear importance to the data analysis of gravitational waves. When compared to the results obtained in the linearized approximation, we observe large phase shifts, a stronger than linear generation of gravitational wave output and considerable generation of radiation in polarization states which are not found in the linearized approximation. In terms of a spherical harmonic decomposition, the nonlinear properties of the harmonic amplitudes have simple scaling properties which offer an economical way to catalog the details of the waves produced in such black hole processes.Comment: 17 pages, 20 figures, abstract and introduction re-writte

A template bank for gravitational waveforms from coalescing binary black holes: non-spinning binaries

Gravitational waveforms from the inspiral and ring-down stages of the binary black hole coalescences can be modelled accurately by approximation/perturbation techniques in general relativity. Recent progress in numerical relativity has enabled us to model also the non-perturbative merger phase of the binary black-hole coalescence problem. This enables us to \emph{coherently} search for all three stages of the coalescence of non-spinning binary black holes using a single template bank. Taking our motivation from these results, we propose a family of template waveforms which can model the inspiral, merger, and ring-down stages of the coalescence of non-spinning binary black holes that follow quasi-circular inspiral. This two-dimensional template family is explicitly parametrized by the physical parameters of the binary. We show that the template family is not only \emph{effectual} in detecting the signals from black hole coalescences, but also \emph{faithful} in estimating the parameters of the binary. We compare the sensitivity of a search (in the context of different ground-based interferometers) using all three stages of the black hole coalescence with other template-based searches which look for individual stages separately. We find that the proposed search is significantly more sensitive than other template-based searches for a substantial mass-range, potentially bringing about remarkable improvement in the event-rate of ground-based interferometers. As part of this work, we also prescribe a general procedure to construct interpolated template banks using non-spinning black hole waveforms produced by numerical relativity.Comment: A typo fixed in Eq.(B11

Implementation of higher-order absorbing boundary conditions for the Einstein equations

We present an implementation of absorbing boundary conditions for the Einstein equations based on the recent work of Buchman and Sarbach. In this paper, we assume that spacetime may be linearized about Minkowski space close to the outer boundary, which is taken to be a coordinate sphere. We reformulate the boundary conditions as conditions on the gauge-invariant Regge-Wheeler-Zerilli scalars. Higher-order radial derivatives are eliminated by rewriting the boundary conditions as a system of ODEs for a set of auxiliary variables intrinsic to the boundary. From these we construct boundary data for a set of well-posed constraint-preserving boundary conditions for the Einstein equations in a first-order generalized harmonic formulation. This construction has direct applications to outer boundary conditions in simulations of isolated systems (e.g., binary black holes) as well as to the problem of Cauchy-perturbative matching. As a test problem for our numerical implementation, we consider linearized multipolar gravitational waves in TT gauge, with angular momentum numbers l=2 (Teukolsky waves), 3 and 4. We demonstrate that the perfectly absorbing boundary condition B_L of order L=l yields no spurious reflections to linear order in perturbation theory. This is in contrast to the lower-order absorbing boundary conditions B_L with L<l, which include the widely used freezing-Psi_0 boundary condition that imposes the vanishing of the Newman-Penrose scalar Psi_0.Comment: 25 pages, 9 figures. Minor clarifications. Final version to appear in Class. Quantum Grav

Conformally flat black hole initial data, with one cylindrical end

We give a complete analytical proof of existence and uniqueness of extreme-like black hole initial data for Einstein equations, which possess a cilindrical end, analogous to extreme Kerr, extreme Reissner Nordstrom, and extreme Bowen-York's initial data. This extends and refines a previous result \cite{dain-gabach09} to a general case of conformally flat, maximal initial data with angular momentum, linear momentum and matter.Comment: Minor changes and formula (21) revised according to the published version in Class. Quantum Grav. (2010). Results unchange