by Cornish, Neil J.
4 pages, no figures

Theoretical studies in gravitational wave astronomy often require the calculation of Fisher Information Matrices and Likelihood functions, which in a direct approach entail the costly step of computing gravitational waveforms. Here I describe an alternative technique that sidesteps the need to compute full waveforms, resulting in significant computational savings. I describe how related techniques can be used to speed up Bayesian inference applied to real gravitational wave data.

by Hannam, Mark and Husa, Sascha and Ohme, Frank and Mueller, Doreen and Bruegmann, Bernd
20 pages, 9 figures, 6 tables

We present gravitational waveforms for the last orbits and merger of black-hole-binary (BBH) systems along two branches of the BBH parameter space: equal-mass binaries with equal non-precessing spins, and nonspinning unequal-mass binaries. The waveforms are calculated from numerical solutions of Einstein’s equations for black-hole binaries that complete between six and ten orbits before merger. Along the equal-mass spinning branch, the spin parameter of each BH is $latex \chi_i = S_i/M_i^2 \in [-0.85,0.85]$, and along the unequal-mass branch the mass ratio is $latex q =M_2/M_1 \in [1,4]$. We discuss the construction of low-eccentricity puncture initial data for these cases, the properties of the final merged BH, and compare the last 8-10 GW cycles up to $latex M\omega = 0.1$ with the phase and amplitude predicted by standard post-Newtonian (PN) approximants. As in previous studies, we find that the phase from the 3.5PN TaylorT4 approximant is most accurate for nonspinning binaries. For equal-mass spinning binaries the 3.5PN TaylorT1 approximant (including spin terms up to only 2.5PN order) gives the most robust performance, but it is possible to treat TaylorT4 in such a way that it gives the best accuracy for spins $latex \chi_i > -0.75$. When high-order amplitude corrections are included, the PN amplitude of the $latex (\ell=2,m=\pm2)$ modes is larger than the NR amplitude by between 2-4%.

by Nichols, David A. and Chen, Yanbei
13 pages, 6 figures, 1 table, revtex4 format

Black-hole-binary coalescence is often divided into three stages: inspiral, merger and ringdown. The post-Newtonian (PN) approximation treats the inspiral phase, black-hole perturbation (BHP) theory describes the ringdown, and the nonlinear dynamics of spacetime characterize the merger. In this paper, we introduce a hybrid method that incorporates elements of PN and BHP theories, and we apply it to the head-on collision of black holes with transverse, anti-parallel spins. We compare our approximation technique with a full numerical-relativity simulation, and we find good agreement between the gravitational waveforms and the radiated energy and momentum. Our results suggest that PN and BHP theories may suffice to explain the main features of outgoing gravitational radiation for head-on mergers. This would further imply that linear perturbations to exact black-hole solutions can capture the nonlinear aspects of head-on binary-black-hole mergers accessible to observers far from the collision.

by Tessmer, Manuel and Schaefer, Gerhard
28 pages

The article provides ready-to-use 1pN-accurate frequency-domain gravitational wave forms for eccentric nonspinning compact binaries of arbitrary mass ratio including the first post-Newtonian (1pN) point particle corrections to the far-zone gravitational wave amplitude, given in terms of tensor spherical harmonics. The averaged equations for the decay of the eccentricity and growth of radial frequency due to radiation reaction are used to provide stationary phase approximations to the frequency-domain wave forms.

by Klein, Antoine and Jetzer, Philippe
7 pages, 1 figure; Accepted for publication in Phys. Rev. D

We compute here the spin-orbit and spin-spin couplings needed for an accurate computation of the phasing of gravitational waves emitted by comparable-mass binaries on eccentric orbits at the second post-Newtonian (PN) order. We use a quasi-Keplerian parametrization of the orbit free of divergencies in the zero eccentricity limit. We find that spin-spin couplings induce a residual eccentricity for coalescing binaries at 2PN, of the order of $latex 10^{-4}$-$latex 10^{-3}$ for supermassive black hole binaries in the LISA band. Spin-orbit precession also induces a non-trivial pattern in the evolution of the eccentricity, which could help to reduce the errors on the determination of the eccentricity and spins in a gravitational wave measurement.

by Finn, Lee Samuel and Lommen, Andrea N.
43 pages, 13 figures, submitted to ApJ.

Efforts to detect gravitational waves by timing an array of pulsars have focused traditionally on stationary gravitational waves: e.g., stochastic or periodic signals. Gravitational wave bursts — signals whose duration is much shorter than the observation period — will also arise in the pulsar timing array waveband. Sources that give rise to detectable bursts include the formation or coalescence of supermassive black holes (SMBHs), the periapsis passage of compact objects in highly elliptic or unbound orbits about a SMBH, or cusps on cosmic strings. Here we describe how pulsar timing array data may be analyzed to detect and characterize these bursts. Our analysis addresses, in a mutually consistent manner, a hierarchy of three questions: \emph{i}) What are the odds that a dataset includes the signal from a gravitational wave burst? \emph{ii}) Assuming the presence of a burst, what is the direction to its source? and \emph{iii}) Assuming the burst propagation direction, what is the burst waveform’s time dependence in each of its polarization states? Applying our analysis to synthetic data sets we find that we can \emph{detect} gravitational waves even when the radiation is too weak to either localize the source of infer the waveform, and \emph{detect} and \emph{localize} sources even when the radiation amplitude is too weak to permit the waveform to be determined. While the context of our discussion is gravitational wave detection via pulsar timing arrays, the analysis itself is directly applicable to gravitational wave detection using either ground or space-based detector data.

by Cornish, Neil J. and Key, Joey Shapiro
9 pages

Several scenarios have been proposed in which the orbits of binary black holes enter the band of a gravitational wave detector with significant eccentricity. To avoid missing these signals or biasing the parameter estimation it is important that we consider waveform models that account for eccentricity. The ingredients needed to compute post-Newtonian (PN) waveforms produced by spinning black holes inspiralling on quasi-eccentric orbits have been available for almost two decades at 2 PN order, and this work has recently been extended to 2.5 PN order. However, the computational cost of directly implementing these waveforms is high, requiring many steps per orbit to evolve the system of coupled differential equations. Here we employ a separation of timescales and a generalized Keplarian parameterization of the orbits to produce efficient waveforms describing spinning black hole binaries with arbitrary spin orientations on quasi-eccentric orbits to 1.5 PN order. Our solution includes the spin contributions to the decay of the semi-major axis and eccentricity. We outline a scheme for extending our approach to higher post-Newtonian order.

by Mueller, Doreen and Grigsby, Jason and Bruegmann, Bernd
10 pages, 14 figures

Certain numerical frameworks used for the evolution of binary black holes make use of a gamma driver, which includes a damping factor. Such simulations typically use a constant value for damping. However, it has been found that very specific values of the damping factor are needed for the calculation of unequal mass binaries. We examine carefully the role this damping plays, and provide two explicit, non-constant forms for the damping to be used with mass-ratios further from one. Our analysis of the resultant waveforms compares well against the constant damping case.

by Favata, Marc
11 pages, 2 figures; proceedings of the 8th Amaldi Conference on Gravitational Waves (New York, June 2009); accepted for publication in special issue of Classical and Quantum Gravity

The nonlinear memory effect is a slowly-growing, non-oscillatory contribution to the gravitational-wave amplitude. It originates from gravitational waves that are sourced by the previously emitted waves. In an ideal gravitational-wave interferometer a gravitational-wave with memory causes a permanent displacement of the test masses that persists after the wave has passed. Surprisingly, the nonlinear memory affects the signal amplitude starting at leading (Newtonian-quadrupole) order. Despite this fact, the nonlinear memory is not easily extracted from current numerical relativity simulations. After reviewing the linear and nonlinear memory I summarize some recent work, including: (1) computations of the memory contribution to the inspiral waveform amplitude (thus completing the waveform to third post-Newtonian order); (2) the first calculations of the nonlinear memory that include all phases of binary black hole coalescence (inspiral, merger, ringdown); and (3) realistic estimates of the detectability of the memory with LISA.

by Tessmer, Manuel and Hartung, Johannes and Schaefer, Gerhard
26 pages, 1 figure

A quasi-Keplerian parameterisation for the solutions of second post-Newtonian (PN) accurate equations of motion for spinning compact binaries is obtained including leading order spin-spin and next-to-leading order spin-orbit interactions. Rotational deformation of the compact objects is incorporated. For arbitrary mass ratios the spin orientations are taken to be parallel or anti-parallel to the orbital angular momentum vector. The emitted gravitational wave forms are given in analytic form up to 2PN point particle, 1.5PN spin orbit and 1PN spin-spin contributions, where the spins are counted of 0PN order.

by Canizares, Priscilla and Sopuerta, Carlos F.
3 pages. To appear in Proceedings of the Twelfth Marcel Grossmann Meeting on General Relativity, edited by Thibault Damour, Robert T Jantzen and Remo Ruffini, World Scientific, Singapore, 2010

We introduce a new time-domain method for computing the self-force acting on a scalar particle in a Schwarzschild geometry. The principal feature of our method consists in the division of the spatial domain into several subdomains and locating the particle at the interface betweem two them. In this way, we avoid the need of resolving a small length scale associated with the presence of a particle in the computational domain and, at the same time, we avoid numerical problems due to the low differentiability of solutions of equations with point-like singular behaviour.

by Barausse, Enrico and Buonanno, Alessandra
22 pages, 9 figures

Building on a recent paper in which we computed the canonical Hamiltonian of a spinning test particle in curved spacetime, at linear order in the particle’s spin, we work out an improved effective-one-body (EOB) Hamiltonian for spinning black-hole binaries. As in previous descriptions, we endow the effective particle not only with a mass m, but also with a spin S*. Thus, the effective particle interacts with the effective Kerr background (having spin S_Kerr) through a geodesic-type interaction and an additional spin-dependent interaction proportional to S*. When expanded in post-Newtonian (PN) orders, the EOB Hamiltonian reproduces the leading order spin-spin coupling and the spin-orbit coupling through 2.5PN order, for any mass-ratio. Also, it reproduces all spin-orbit couplings in the test-particle limit. Similarly to the test-particle limit case, when we restrict the EOB dynamics to spins aligned or antialigned with the orbital angular momentum, for which circular orbits exist, the EOB dynamics has several interesting features, such as the existence of an innermost stable circular orbit, a photon circular orbit, and a maximum in the orbital frequency during the plunge subsequent to the inspiral. These properties are crucial for reproducing the dynamics and gravitational-wave emission of spinning black-hole binaries, as calculated in numerical relativity simulations.

by Pan, Yi and Buonanno, Alessandra and Buchman, Luisa T. and Chu, Tony and Kidder, Lawrence E. and Pfeiffer, Harald P. and Scheel, Mark A.
15 pages, 8 figures

We present the first attempt at calibrating the effective-one-body (EOB) model to accurate numerical-relativity simulations of spinning, non-precessing black-hole binaries. Aligning the EOB and numerical waveforms at low frequency over a time interval of 1000M, we first estimate the phase and amplitude errors in the numerical waveforms and then minimize the difference between numerical and EOB waveforms by calibrating a handful of EOB-adjustable parameters. In the equal-mass, spin aligned case, we find that phase and fractional amplitude differences between the numerical and EOB (2,2) mode can be reduced to 0.01 radians and 1%, respectively, over the entire inspiral waveforms. In the equal-mass, spin anti-aligned case, these differences can be reduced to 0.13 radians and 1% during inspiral and plunge, and to 0.4 radians and 10% during merger and ringdown. The waveform agreement is within numerical errors in the spin aligned case while slightly over numerical errors in the spin anti-aligned case. Using Enhanced LIGO and Advanced LIGO noise curves, we find that the overlap between the EOB and the numerical (2,2) mode, maximized over the initial phase and time of arrival, is larger than 0.999 for binaries with total mass 30-200Ms. In addition to the leading (2,2) mode, we compare four subleading modes. We find good amplitude and frequency agreements between the EOB and numerical modes for both spin configurations considered, except for the (3,2) mode in the spin anti-aligned case. We believe that the larger difference in the (3,2) mode is due to the lack of knowledge of post-Newtonian spin effects in the higher modes.

by Hadar, Shahar and Kol, Barak
17 pages, 2 figures. v3: minor changes

An extreme mass ratio inspiral consists of two parts: adiabatic inspiral and plunge. The plunge trajectory from the innermost stable circular orbit (ISCO) is special (somewhat independent of initial conditions). We write an expression for its solution in closed-form and for the emitted waveform. In particular we extract an expression for the associated black-hole ringdown amplitudes, which reduces to solving certain ordinary (radial) differential equations in the Schwarzschild background.

by Damour, Thibault and Nagar, Alessandro
21 pages, 5 figures. Submitted to Phys. Rev. D

The late part of the gravitational wave signal of binary neutron star inspirals can in principle yield crucial information on the nuclear equation of state via its dependence on relativistic tidal parameters. In the hope of analytically describing the gravitational wave phasing during the late inspiral (essentially up to contact) we propose an extension of the effective one body (EOB) formalism which includes tidal effects. We compare the prediction of this tidal-EOB formalism to recently computed nonconformally flat quasi-equilibrium circular sequences of binary neutron star systems. Our analysis suggests the importance of higher-order (post-Newtonian) corrections to tidal effects, even beyond the first post-Newtonian order, and their tendency to {\it significantly} increase the “effective tidal polarizability” of neutron stars. We compare the EOB predictions to some recently advocated, nonresummed, post-Newtonian based (“Taylor-T4”) description of the phasing of inspiralling systems. This comparison shows the strong sensitivity of the late-inspiral phasing to the choice of the analytical model, but raises the hope that a sufficiently accurate numerical–relativity–“calibrated” EOB model might give us a reliable handle on the nuclear equation of state

by Tessmer, Manuel
13 pages, 2 figures, submitted to PRD on 11 Sep. 2009

The paper generalizes the structure of gravitational waves from orbiting spinning binaries under leading order spin-orbit coupling, as given in the work by K\”onigsd\”orffer and Gopakumar [PRD 71, 024039 (2005)] for single-spin and equal-mass binaries, to unequal-mass binaries and arbitrary spin configurations. The orbital motion is taken to be quasi-circular and the fractional mass difference is assumed to be small against one. The emitted gravitational waveforms are given in analytic form.

by Tiec, Alexandre Le and Blanchet, Luc
34 pages, 6 figures

The ringdown phase of a black hole formed from the merger of two orbiting black holes is described by means of the close-limit (CL) approximation starting from second-post-Newtonian (2PN) initial conditions. The 2PN metric of point-particle binaries is formally expanded in CL form and identified with that of a perturbed Schwarzschild black hole. The multipolar coefficients describing the even-parity (or polar) and odd-parity (axial) components of the linear perturbation consistently satisfy the 2PN-accurate perturbative field equations. We use these coefficients to build initial conditions for the Regge-Wheeler and Zerilli wave equations, which we then evolve numerically. The ringdown waveform is obtained in two cases: head-on collision with zero-angular momentum, composed only of even modes, and circular orbits, for which both even and odd modes contribute. In a separate work, this formalism is applied to the study of the gravitational recoil produced during the ringdown phase of coalescing binary black holes.

by Amaro-Seoane, Pau and Santamaria, Lucia
Submitted to ApJ; abstract abridged, figure 1 has a lower resolution

Even though the existence of intermediate-mass black holes has not yet been corroborated observationally, these objects are of high interest for astrophysics. Our understanding of formation and evolution of supermassive black holes (SMBHs), as well as galaxy evolution modeling and cosmography would dramatically change if an IMBH was observed. The prospect of detection and, possibly, observation and characterization of an IMBH has good chances in lower-frequency gravitational-wave (GW) astrophysics with ground-based detectors such as LIGO, Virgo and the future Einstein Telescope (ET). We present an analysis of the signal of a system of a binary of IMBHs based on a waveform model obtained with numerical relativity simulations coupled with post-Newtonian calculations at the highest available order so as to extend the waveform to lower frequencies. We find that initial LIGO and Virgo are in the position of detecting IMBHs with a signal-to-noise ratio (SNR) of $latex \sim 10$ for systems with total mass between 100 and $latex 500 M_{\odot}$ situated at a distance of 100 Mpc. Nevertheless, the event rate is too low and the possibility that these signals are mistaken with a glitch is, unfortunately, non-negligible. When going to second- and third-generation detectors, such as Advanced LIGO or the proposed ET, the event rate becomes much more promising (tens per year for the first and thousands per year for the latter) and the SNR at 100 Mpc is as high as 100 — 1000 and 1000 — $latex 10^{5}$ respectively. The prospects for IMBH detection and characterization with ground-based GW observatories would not only provide us with a robust test of general relativity, but would also corroborate the existence of these systems. Such detections would be a probe to the stellar environments of IMBHs and their formation.

by Yunes, Nicolas and Buonanno, Alessandra and Hughes, Scott A. and Miller, M. Coleman and Pan, Yi
4 pages, 3 figures, submitted to Phys. Rev. Letters

We present the first models of extreme-mass-ratio inspirals within the effective-one-body (EOB) formalism, focusing on quasi-circular orbits into non-rotating black holes. We show that the phase difference and (Newtonian normalized) amplitude difference between EOB and Teukolsky-based gravitational waveforms can be reduced to < 0.1 rads and < 0.002, respectively, after a 2-year evolution. The inclusion of post-Newtonian self-force terms in the EOB approach leads to a phase disagreement of ~6-27 rads after a 2-year evolution. Such inclusion could also allow for the EOB modeling of waveforms from intermediate-mass ratio, quasi-circular inspirals.

by Khanna, Gaurav and McKennon, Justin
8 pages, 4 figures, Accepted for publication in Parallel and Distributed Computing and Systems (PDCS 2009)

In this paper, we accelerate a gravitational physics numerical modelling application using hardware accelerators — Cell processor and Tesla CUDA GPU. We describe these new technologies and our approach in detail, and then present our final performance results. We obtain well over an order-of-magnitude performance gain in our application by making use of these many-core architectures.

by Ajith, P. and Hannam, M. and Husa, S. and Chen, Y. and Bruegmann, B. and Dorband, N. and Mueller, D. and Ohme, F. and Pollney, D. and Reisswig, C. and Santamaria, L. and Seiler, J.
4 pages, 5 figures

We present the first analytical inspiral-merger-ringdown gravitational waveforms from black-hole (BH) binaries with non-precessing spins. By matching a post-Newtonian description of the inspiral to a set of numerical calculations performed in full general relativity, we obtain a waveform family with a conveniently small number of physical parameters. The physical content of these waveforms includes the “orbital hang-up” effect, when BHs are spinning rapidly along the direction of the orbital angular momentum. These waveforms will allow us to detect a larger parameter space of BH binary coalescence, to explore various scientific questions related to GW astronomy, and could dramatically improve the expected detection rates of GW detectors.

by Capozziello, S. and De Laurentis, M. and Forte, L. and Garufi, F. and Milano, L.
14 pages, 7 figures

Gravitational waveforms and production could be considerably affected by gravitomagnetic corrections considered in relativistic theory of orbits. Beside the standard periastron effect of General Relativity, new nutation effects come out when c^{-3} corrections are taken into account. Such corrections emerge as soon as matter-current densities and vector gravitational potentials cannot be discarded into dynamics. We study the gravitational waves emitted through the capture, in the gravitational field of massive binary systems (e.g. a very massive black hole on which a stellar object is inspiralling) via the quadrupole approximation, considering precession and nutation effects. We present a numerical study to obtain the gravitational wave luminosity, the total energy output and the gravitational radiation amplitude. From a crude estimate of the expected number of events towards peculiar targets (e.g. globular clusters) and in particular, the rate of events per year for dense stellar clusters at the Galactic Center (SgrA*), we conclude that this type of capture could give signatures to be revealed by interferometric GW antennas, in particular by the forthcoming laser interferometer space antenna LISA.

by Whelan, John T. and Prix, Reinhard and Khurana, Deepak
26 pages, 11 figures

We describe the application of an F-statistic search for continuous gravitational waves to the search for galactic white-dwarf binaries in the Second Mock LISA Data Challenge. The search method employs a hierarchical template-grid based exploration of the parameter space, using a coincidence step to distinguish between primary (”true”) and secondary maxima, followed by a final (multi-TDI) “zoom” stage to provide an accurate parameter estimation of the final candidates. Suitably tuned, the pipeline is able to extract 1989 true signals with only 5 false alarms. The use of the rigid adiabatic approximation allows recovery of signal parameters comparable to statistical expectations, although there is still some systematic excess above expected statistical errors due to Gaussian noise. An experimental iterative pipeline with seven rounds of subtraction and re-analysis allows us to increase the number of signals recovered, up to a total of 3419 with 29 false alarms.

by Harry, Ian and Allen, Bruce and Sathyaprakash, B. S.
14 pages, 11 figures

This paper presents an algorithm for constructing matched-filter template banks in an arbitrary parameter space. The method places templates at random, then removes those which are “too close” together. The properties and optimality of stochastic template banks generated in this manner are investigated for some simple models. The effectiveness of these template banks for gravitational wave searches for binary inspiral waveforms is also examined. The properties of a stochastic template bank are then compared to the deterministically placed template banks that are currently used in gravitational wave data analysis.

by Asada, Hideki
21 pages, 3 figures, 1 table

Continuing work initiated in an earlier publication [Torigoe et al. Phys. Rev. Lett. {\bf 102}, 251101 (2009)], gravitational wave forms for a three-body system in Lagrange’s orbit are considered especially in an analytic method. First, we derive an expression of the three-body wave forms at the mass quadrupole, octupole and current quadrupole orders. By using the expressions, we solve a gravitational-wave {\it inverse} problem of determining the source parameters to this particular configuration (three masses, a distance of the source to an observer, and the orbital inclination angle to the line of sight) through observations of the gravitational wave forms alone. For this purpose, the chirp mass to a three-body system in the particular configuration is expressed in terms of only the mass ratios by deleting initial angle positions. We discuss also whether and how a binary source can be distinguished from a three-body system in Lagrange’s orbit or others.

by Klein, Antoine and Jetzer, Philippe and Sereno, Mauro
20 pages, 20 figures. Submitted to Phys. Rev. D

Gravitational waves emitted by binary systems in the inspiral phase carry a complicated structure, consisting in a superposition of different harmonics of the orbital frequency, the amplitude of each of them taking the form of a Post-Newtonian series. In addition to that, spinning binaries experience couplings which induce a precession of the orbital angular momentum and of the individual spins. So far, all studies of the measurement accuracy of gravitational wave experiments for comparable-mass binary systems have considered either spinless binaries, or spinning binaries without subdominant harmonics in the waveform as well as no amplitude modulations. In this paper, we consider supermassive black hole binaries as expected to be observed with the planned space-based interferometer LISA, and study the measurement accuracy for several astrophysically interesting parameters obtainable taking into account the full 2PN waveform for spinning bodies, as well as spin-precession effects. We find that for binaries with a total mass in the range 10^5 M_Sun < M 10 for M < 10^7 M_Sun, 1.5 – 5 times higher than with the restricted waveform. We computed that the full waveform allows to use supermassive black hole binaries as standard sirens up to a redshift of z = 1.6, about 0.4 larger than what previous studies allowed.

by Reisswig, C. and Bishop, N. T. and Pollney, D. and Szilagyi, B.
4 pages, 3 figures

The accurate calculation of gravitational radiation emitted during black hole merger events is important for both improving the chances of detecting such events, as well as for making an astrophysical interpretation of the data once an event has been detected. Gravitational radiation is properly defined only at future null infinity (denoted by Scri), but in practice it is estimated from data calculated at a finite radius. We have used characteristic extraction to calculate gravitational radiation at Scri for the real astrophysical problem of the inspiral and merger of two equal mass non-spinning black holes (this problem has become a standard benchmark case in numerical relativity). Thus we have determined the first unambiguous merger waveforms for this problem. The implementation that has been developed is general purpose, and can be applied to calculate the gravitational radiation, at Scri, given data at a finite radius calculated in another computation.

by Johnson-McDaniel, Nathan K. and Yunes, Nicolas and Tichy, Wolfgang and Owen, Benjamin J.
48 pages, 23 figures

(Abridged) By asymptotically matching a post-Newtonian (PN) metric to two tidally perturbed Schwarzschild metrics, we generate approximate initial data (in the form of a 4-metric) for a nonspinning black hole binary in a circular orbit. We carry out this matching through O(v^4) in the binary’s orbital velocity v, so the resulting data are conformally curved. Far from the holes, we use the appropriate PN metric that accounts for retardation, which we construct using the highest-order PN expressions available to compute the binary’s past history. The data set’s uncontrolled remainders are thus O(v^5) throughout the timeslice; we also generate an extension to the data set that has uncontrolled remainders of O(v^6) in the purely PN portion of the timeslice (i.e., not too close to the holes). The resulting data are smooth, since we join all the metrics together by smoothly interpolating between them. We perform this interpolation using transition functions constructed to avoid introducing excessive additional constraint violations. Due to their inclusion of tidal deformations and outgoing radiation, these data should substantially reduce the initial spurious (”junk”) radiation observed in current simulations that use conformally flat initial data. Such reductions in the nonphysical components of the initial data will be necessary for simulations to achieve the accuracy required to supply Advanced LIGO and LISA with the templates necessary for parameter estimation.

by Buonanno, Alessandra and Iyer, Bala and Ochsner, Evan and Pan, Yi and Sathyaprakash, B. S.
27 pages, 8 figures, 4 tables, submitted to PRD

The two-body dynamics in general relativity has been solved perturbatively using the post-Newtonian (PN) approximation. The evolution of the orbital phase and the emitted gravitational radiation are now known to a rather high order up to O(v^8), v being the characteristic velocity of the binary. The orbital evolution, however, cannot be specified uniquely due to the inherent freedom in the choice of parameter used in the PN expansion as well as the method pursued in solving the relevant differential equations. The goal of this paper is to determine the (dis)agreement between different PN waveform families in the context of initial and advanced gravitational-wave detectors. The waveforms employed in our analysis are those that are currently used by Initial LIGO/Virgo, that is the time-domain PN models TaylorT1, TaylorT2, TaylorT3, TaylorT4 and TaylorEt, the effective one-body (EOB) model, and the Fourier-domain representation TaylorF2. We examine the overlaps of these models with one another and with the prototype effective one-body model (calibrated to numerical relativity simulations, as currently used by initial LIGO) for a number of different binaries at 2PN, 3PN and 3.5PN orders to quantify their differences and to help us decide whether there exist preferred families that are the most appropriate as search templates. We conclude that as long as the total mass remains less than a certain upper limit M_crit, all template families at 3.5PN order (except TaylorT3 and TaylorEt) are equally good for the purpose of detection. The value of M_crit is found to be ~ 12M_Sun for Initial, Enhanced and Advanced LIGO. From a purely computational point of view we recommend that 3.5PN TaylorF2 be used below Mcrit and EOB calibrated to numerical relativity simulations be used for total binary mass M > Mcrit.

by Reisswig, Christian and Husa, Sascha and Rezzolla, Luciano and Dorband, Ernst Nils and Pollney, Denis and Seiler, Jennifer
18 pages, 11 figures

Binary black-hole systems with spins aligned or anti-aligned to the orbital angular momentum provide the natural ground to start detailed studies of the influence of strong-field spin effects on gravitational wave observations of coalescing binaries. Furthermore, such systems may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins parallel to the orbital angular momentum but otherwise arbitrary. Our attention is particularly focused on the gravitational-wave emission so as to quantify how much spin effects contribute to the signal-to-noise ratio, to the horizon distances, and to the relative event rates for the representative ranges in masses and detectors. As expected, the signal-to-noise ratio increases with the projection of the total black hole spin in the direction of the orbital momentum. We find that equal-spin binaries with maximum spin aligned with the orbital angular momentum are more than “three times as loud” as the corresponding binaries with anti-aligned spins, thus corresponding to event rates up to 30 times larger. We also consider the waveform mismatch between the different spinning configurations and find that, within our numerical accuracy, binaries with opposite spins S_1=-S_2 cannot be distinguished whereas binaries with spin S_1=S_2 have clearly distinct gravitational-wave emissions. Finally, we derive a simple expression for the energy radiated in gravitational waves and find that the binaries always have efficiencies E_rad/M > 3.6%, which can become as large as E_rad/M = 10% for maximally spinning binaries with spins aligned with the orbital angular momentum.

by Lindblom, Lee
9 pages, 7 figures

Accuracy standards have been developed to ensure that the waveforms used for gravitational-wave data analysis are good enough to serve their intended purposes. These standards place constraints on certain norms of the frequency-domain representations of the waveform errors. Examples are given here of possible misinterpretations and misapplications of these standards, whose effect could be to vitiate the quality control they were intended to enforce. Suggestions are given for ways to avoid these problems.

by Capozziello, Salvatore and De Laurentis, Mariafelicia and Forte, Luca and Garufi, Fabio and Milano, Leopoldo
6 pages, 10 figures; Multifrequency Behaviour of High-Energy Cosmic Sources, Vulcano Workshop 2009

Corrections to the relativistic theory of orbits are discussed considering higher order approximations induced by gravitomagnetic effects. Beside the standard periastron effect of General Relativity (GR), a new nutation effect was found due to the $latex {\displaystyle c^{-3}}$ orbital correction. According to the presence of that new nutation effect we studied the gravitational waveforms emitted through the capture in a gravitational field of a massive black hole (MBH) of a compact object (neutron star (NS) or BH) via the quadrupole approximation. We made a numerical study to obtain the emitted gravitational wave (GW) amplitudes. We conclude that the effects we studied could be of interest for the future space laser interferometric GW antenna LISA.

by Chung, Shin Kee and Wen, Linqing and Blair, David and Cannon, Kipp and Datta, Amitava
12 pages, 4 figures, submitted to CQG

We report a novel application of graphics processing units (GPUs) for the purpose of accelerating the search pipelines for gravitational waves from coalescing binaries of compact objects. A speed-up of 16 fold has been achieved compared with a single central processing unit (CPU). We show that substantial improvements are possible and discuss the reduction in CPU count required for the detection of inspiral sources afforded by the use of GPUs.

by Yagi, Kent and Tanaka, Takahiro
37 pages, 16 figures

We calculate how strong one can put constraints on the alternative theories of gravities such as Brans-Dicke and massive graviton theories with LISA. We consider the inspiral gravitational waves from NS/IMBH binaries in Brans-Dicke theory and SMBH/BH binaries in massive graviton theories. We use the 2PN waveforms including spins. We also take both precession and small eccentricity of the orbit into account. We neglect the spin of one of the binary object so that we can apply the so-called \textit{simple precession}. We perform the Monte Carlo simulations of $latex 10^4$ binaries, whose parameters include the Brans-Dicke parameter $latex \omega_{\mathrm{BD}}$ and the graviton Compton length $latex \lambda_g$. We find that including both the spin-spin coupling $latex \sigma$ and the small eccentricity into the binary parameters reduces the determination accuracy by an order of magnitude for the Brans-Dicke case, whilst it has less influence on massive graviton theories. On the other hand, including precession enhances the constraint on $latex \omega_{\mathrm{BD}}$ only 20% but it increases the constraint on $latex \lambda_g$ by several factors. For $latex (1.4+1000)M_{\odot}$ NS/BH binaries of SNR=10, one can put $latex \omega_{\mathrm{BD}}>7040$, whilst for $latex (10^7+10^6)M_{\odot}$ BH/BH binaries at 3Gpc, one can put $latex \lambda_g>4.24\times10^{21}$cm, on average. This is four orders of magnitude stronger than the one obtained from the solar system experiment. From these results, it is understood that the effects of precession and eccentricity cannot be neglected in the parameter estimation analysis.

by Berti, Emanuele and Cardoso, Vitor and Starinets, Andrei O.
112 pages, 20 figures, 4 tables. Invited Topical Review for CQG. Comments and suggestions are welcome. QNM data and Mathematica notebooks are available at http://www.phy.olemiss.edu/~berti/qnms.html and http://gamow.ist.utl.pt/~vitor/ringdown

Quasinormal modes are eigenmodes of dissipative systems. Perturbations of classical gravitational backgrounds involving black holes or branes naturally lead to quasinormal modes. The analysis and classification of the quasinormal spectra requires solving non-Hermitian eigenvalue problems for the associated linear differential equations. Within the recently developed gauge-gravity duality, these modes serve as an important tool for determining the near-equilibrium properties of strongly coupled quantum field theories, in particular their transport coefficients, such as viscosity, conductivity and diffusion constants. In astrophysics, the detection of quasinormal modes in gravitational wave experiments would allow precise measurements of the mass and spin of black holes as well as new tests of general relativity. This review is meant as an introduction to the subject, with a focus on the recent developments in the field.

by Yunes, Nicolas and Arun, K. G. and Berti, Emanuele and Will, Clifford M.
24 pages, 9 figures, submitted to Phys. Rev. D

We lay the foundations for the construction of analytic expressions for Fourier-domain gravitational waveforms produced by eccentric, inspiraling compact binaries in a post-circular or small-eccentricity approximation. The time-dependent, “plus” and “cross” polarizations are expanded in Bessel functions, which are then self-consistently re-expanded in a power series about zero initial eccentricity to eighth order. The stationary phase approximation is then employed to obtain explicit analytic expressions for the Fourier transform of the post-circular expanded, time-domain signal. We exemplify this framework by considering Newtonian-accurate waveforms, which in the post-circular scheme give rise to higher harmonics of the orbital phase and amplitude corrections both to the amplitude and the phase of the Fourier domain waveform. Such higher harmonics lead to an effective increase in the inspiral mass reach of a detector as a function of the binary’s eccentricity $latex e_0$ at the time when the binary enters the detector sensitivity band. Using the largest initial eccentricity allowed by our approximations ($latex e_0 < 0.4$), the mass reach is found to be enhanced up to factors of approximately 5 relative to that of circular binaries for Advanced LIGO, LISA, and the proposed Einstein Telescope at a signal-to-noise ratio of ten. A post-Newtonian generalization of the post circular scheme is also discussed, which holds the promise to provide “ready-to-use” Fourier-domain waveforms for data analysis of eccentric inspirals.

by Arun, K G and Will, Clifford M
12 pages, 4 figures

Observations by laser interferometric detectors of gravitational waves from inspiraling compact binary systems can be used to search for a dependence of the waves’ propagation speed on wavelength, and thereby to bound the mass or Compton wavelength of a putative graviton. We study the effect of including higher harmonics, as well as their post-Newtonian amplitude corrections, in the template gravitational waveforms employed in the process of parameter estimation using matched filtering. We consider the bounds that could be achieved using advanced LIGO, a proposed third generation instrument called Einstein Telescope, and the proposed space interferometer LISA. We find that in all cases, the bounds on the graviton Compton wavelength are improved by almost an order of magnitude for higher masses when amplitude corrections are included.

by Field, Scott E. and Hesthaven, Jan S. and Lau, Stephen R.
27 pages, 7 figures, 3 tables. Uses revtex4

Gravitational wave emission from extreme-mass-ratio binaries (EMRBs) should be detectable by the joint NASA-ESU LISA project, spurring interest in analytical and numerical methods for investigating EMRBs. We describe a discontinuous Galerkin (dG) method for solving the distributionally forced 1+1 wave equations which arise when modeling EMRBs via the perturbation theory of Schwarzschild black holes. Despite the presence of jump discontinuities in the relevant polar and axial gravitational “master functions”, our dG method achieves global spectral accuracy, provided that we know the instantaneous position, velocity, and acceleration of the small particle. Here these variables are known, since we assume that the particle follows a timelike geodesic of the Schwarzschild geometry. We document the results of several numerical experiments testing our method, and discuss the possible incorporation of radiation reaction in the model.

by Buonanno, Alessandra and Pan, Yi and Pfeiffer, Harald P. and Scheel, Mark A. and Buchman, Luisa T. and Kidder, Lawrence E.
19 pages, 19 figures

We calibrate the effective-one-body (EOB) model to an accurate numerical simulation of an equal-mass, non-spinning binary black-hole coalescence produced by the Caltech-Cornell collaboration. Aligning the EOB and numerical waveforms at low frequency over a time interval of ~1000M, and taking into account the uncertainties in the numerical simulation, we investigate the significance and degeneracy of the EOB adjustable parameters during inspiral, plunge and merger, and determine the minimum number of EOB adjustable parameters that achieves phase and amplitude agreements on the order of the numerical error. We find that phase and fractional amplitude differences between the numerical and EOB values of the dominant gravitational wave mode h_{22} can be reduced to 0.02 radians and 2%, respectively, until a time 26 M before merger, and to 0.1 radians and 10%, at a time 16M after merger (during ringdown), respectively. Using LIGO, Enhanced LIGO and Advanced LIGO noise curves, we find that the overlap between the EOB and the numerical h_{22}, maximized only over the initial phase and time of arrival, is larger than 0.999 for equal-mass binary black holes with total mass 30-150 Msun. In addition to the leading gravitational mode (2,2), we compare the dominant subleading modes (4,4) and (3,2) and find phase and amplitude differences on the order of the numerical error. We also determine the mass-ratio dependence of one of the EOB adjustable parameters by fitting to numerical {\it inspiral} waveforms for black-hole binaries with mass ratios 2:1 and 3:1. These results improve and extend recent successful attempts aimed at providing gravitational-wave data analysts the best analytical EOB model capable of interpolating accurate numerical simulations.

by Damour, Thibault and Nagar, Alessandro
5 pages, 5 figures, to apper as a Phys. Rev. D Rapid Communication

We present an analytical formalism, within the Effective-One-Body framework, which predicts gravitational-wave signals from inspiralling and coalescing black-hole binaries that agree, within numerical errors, with the results of the currently most accurate numerical relativity simulations for several different mass ratios. In the equal-mass case, the gravitational wave energy flux predicted by our formalism agrees, within numerical errors, with the most accurate numerical-relativity energy flux. We think that our formalism opens a realistic possibility of constructing a sufficiently accurate, large bank of gravitational wave templates, as needed both for detection and data analysis of (non spinning) coalescing binary black holes.

by Santamaria, Lucia and Krishnan, Badri and Whelan, John T.
13 pages, 4 figures, to appear in the proceedings of the NRDA08 meeting, Syracuse, Aug. 11-14, 2008

Recent progress in numerical relativity now allows computation of the binary black hole merger, whereas post-Newtonian and perturbative techniques can be used to model the inspiral and ringdown phases. So far, most gravitational-wave searches have made use of various post-Newtonian-inspired templates to search for signals arising from the coalescence of compact binary objects. Ajith et al have produced hybrid waveforms for non-spinning binary black-hole systems which include the three stages of the coalescence process, and constructed from them phenomenological templates which capture the features of these waveforms in a parametrized form. As a first step towards extending the present inspiral searches to higher-mass binary black-hole systems, we have used these phenomenological waveforms in a search for numerically-simulated signals injected into synthetic LIGO data as part of the NINJA project.

by Aylott, Benjamin and Baker, John G. and Boggs, William D. and Boyle, Michael and Brady, Patrick R. and Brown, Duncan A. and Brügmann, Bernd and Buchman, Luisa T. and Buonanno, Alessandra and Cadonati, Laura and Camp, Jordan and Campanelli, Manuela and Centrella, Joan and Chatterji, Shourov and Christensen, Nelson and Chu, Tony and Diener, Peter and Dorband, Nils and Etienne, Zachariah B. and Faber, Joshua and Fairhurst, Stephen and Farr, Benjamin and Fischetti, Sebastian and Guidi, Gianluca and Goggin, Lisa M. and Hannam, Mark and Herrmann, Frank and Hinder, Ian and Husa, Sascha and Kalogera, Vicky and Keppel, Drew and Kidder, Lawrence E. and Kelly, Bernard J. and Krishnan, Badri and Laguna, Pablo and Lousto, Carlos O. and Mandel, Ilya and Marronetti, Pedro and Matzner, Richard and McWilliams, Sean T. and Matthews, Keith D. and Mercer, R. Adam and Mohapatra, Satyanarayan R. P. and Mroué, Abdul H. and Nakano, Hiroyuki and Ochsner, Evan and Pan, Yi and Pekowsky, Larne and Pfeiffer, Harald P. and Pollney, Denis and Pretorius, Frans and Raymond, Vivien and Reisswig, Christian and Rezzolla, Luciano and Rinne, Oliver and Robinson, Craig and Röver, Christian and Santamaría, Lucía and Sathyaprakash, Bangalore and Scheel, Mark A. and Schnetter, Erik and Seiler, Jennifer and Shapiro, Stuart L. and Shoemaker, Deirdre and Sperhake, Ulrich and Stroeer, Alexander and Sturani, Riccardo and Tichy, Wolfgang and Liu, Yuk Tung and van der Sluys, Marc and van Meter, James R. and Vaulin, Ruslan and Vecchio, Alberto and Veitch, John and Viceré, Andrea and Whelan, John T. and Zlochower, Yosef
53 pages, 25 figures

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

by Hannam, Mark
22 pages, 6 figures, Version to be published in CQG, NRDA 2008 Special Issue

It is now possible to theoretically calculate the gravitational-wave signal from the inspiral, merger and ringdown of a black-hole-binary system. The late inspiral, merger and ringdown can be calculated in full general relativity using numerical methods. The numerical waveforms can then be either stitched to inspiral waveforms predicted by approximation techniques (in particular post-Newtonian calculations) that start at an arbitrarily low frequency, or used to calibrate free parameters in analytic models of the full waveforms. In this review I summarize the status of numerical-relativity (NR) waveforms that include at least ten cycles of the dominant mode of the GW signal before merger, which should be long enough to produce accurate, complete waveforms for GW observations.