Radiation reaction in the 2.5PN waveform from inspiralling binaries in circular orbits

In this Comment we compute the contributions of the radiation reaction force in the 2.5 post-Newtonian (PN) gravitational wave polarizations for compact binaries in circular orbits. (i) We point out and correct an inconsistency in the derivation of Arun, Blanchet, Iyer, and Qusailah. (ii) We prove that all contributions from radiation reaction in the 2.5PN waveform are actually negligible since they can be absorbed into a modification of the orbital phase at the 5PN order.Comment: 7 pages, no figures, submitted to CQ

Hyperboloidal evolution of test fields in three spatial dimensions

We present the numerical implementation of a clean solution to the outer boundary and radiation extraction problems within the 3+1 formalism for hyperbolic partial differential equations on a given background. Our approach is based on compactification at null infinity in hyperboloidal scri fixing coordinates. We report numerical tests for the particular example of a scalar wave equation on Minkowski and Schwarzschild backgrounds. We address issues related to the implementation of the hyperboloidal approach for the Einstein equations, such as nonlinear source functions, matching, and evaluation of formally singular terms at null infinity.Comment: 10 pages, 8 figure

Solving Einstein's Equations With Dual Coordinate Frames

A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad basis. These field components are then determined as functions of a second independent coordinate system. The transformation to the second coordinate system can be thought of as a mapping from the original ``inertial'' coordinate system to the computational domain. This dual-coordinate method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame at infinity; such evolutions are found to be generically unstable using a single rotating coordinate frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with a feedback control system that keeps the excision boundaries of the holes within their respective apparent horizons.Comment: Updated to agree with published versio

Evolution systems for non-linear perturbations of background geometries

The formulation of the initial value problem for the Einstein equations is at the heart of obtaining interesting new solutions using numerical relativity and still very much under theoretical and applied scrutiny. We develop a specialised background geometry approach, for systems where there is non-trivial a priori knowledge about the spacetime under study. The background three-geometry and associated connection are used to express the ADM evolution equations in terms of physical non-linear deviations from that background. Expressing the equations in first order form leads naturally to a system closely linked to the Einstein-Christoffel system, introduced by Anderson and York, and sharing its hyperbolicity properties. We illustrate the drastic alteration of the source structure of the equations, and discuss why this is likely to be numerically advantageous.Comment: 12 pages, 3 figures, Revtex v3.0. Revised version to appear in Physical Review

On the Performance Prediction of BLAS-based Tensor Contractions

Tensor operations are surging as the computational building blocks for a variety of scientific simulations and the development of high-performance kernels for such operations is known to be a challenging task. While for operations on one- and two-dimensional tensors there exist standardized interfaces and highly-optimized libraries (BLAS), for higher dimensional tensors neither standards nor highly-tuned implementations exist yet. In this paper, we consider contractions between two tensors of arbitrary dimensionality and take on the challenge of generating high-performance implementations by resorting to sequences of BLAS kernels. The approach consists in breaking the contraction down into operations that only involve matrices or vectors. Since in general there are many alternative ways of decomposing a contraction, we are able to methodically derive a large family of algorithms. The main contribution of this paper is a systematic methodology to accurately identify the fastest algorithms in the bunch, without executing them. The goal is instead accomplished with the help of a set of cache-aware micro-benchmarks for the underlying BLAS kernels. The predictions we construct from such benchmarks allow us to reliably single out the best-performing algorithms in a tiny fraction of the time taken by the direct execution of the algorithms.Comment: Submitted to PMBS1

Non-precessional spin-orbit effects on gravitational waves from inspiraling compact binaries to second post-Newtonian order

We derive all second post-Newtonian (2PN), non-precessional effects of spin- orbit coupling on the gravitational wave forms emitted by an inspiraling binary composed of spinning, compact bodies in a quasicircular orbit. Previous post- Newtonian calculations of spin-orbit effects (at 1.5PN order) relied on a fluid description of the spinning bodies. We simplify the calculations by introducing into post-Newtonian theory a delta-function description of the influence of the spins on the bodies' energy-momentum tensor. This description was recently used by Mino, Shibata, and Tanaka (MST) in Teukolsky-formalism analyses of particles orbiting massive black holes, and is based on prior work by Dixon. We compute the 2PN contributions to the wave forms by combining the MST energy-momentum tensor with the formalism of Blanchet, Damour, and Iyer for evaluating the binary's radiative multipoles, and with the well-known 1.5PN order equations of motion for the binary. Our results contribute at 2PN order only to the amplitudes of the wave forms. The secular evolution of the wave forms' phase, the quantity most accurately measurable by LIGO, is not affected by our results until 2.5PN order, at which point other spin-orbit effects also come into play. We plan to evaluate the entire 2.5PN spin-orbit contribution to the secular phase evolution in a future paper, using the techniques of this paper.Comment: 11 pages, submitted to Phys. Rev.

A Numerical Relativity Waveform Surrogate Model for Generically Precessing Binary Black Hole Mergers

A generic, non-eccentric binary black hole (BBH) system emits gravitational waves (GWs) that are completely described by 7 intrinsic parameters: the black hole spin vectors and the ratio of their masses. Simulating a BBH coalescence by solving Einstein's equations numerically is computationally expensive, requiring days to months of computing resources for a single set of parameter values. Since theoretical predictions of the GWs are often needed for many different source parameters, a fast and accurate model is essential. We present the first surrogate model for GWs from the coalescence of BBHs including all $7$ dimensions of the intrinsic non-eccentric parameter space. The surrogate model, which we call NRSur7dq2, is built from the results of $744$ numerical relativity simulations. NRSur7dq2 covers spin magnitudes up to $0.8$ and mass ratios up to $2$, includes all $\ell \leq 4$ modes, begins about $20$ orbits before merger, and can be evaluated in $\sim~50\,\mathrm{ms}$. We find the largest NRSur7dq2 errors to be comparable to the largest errors in the numerical relativity simulations, and more than an order of magnitude smaller than the errors of other waveform models. Our model, and more broadly the methods developed here, will enable studies that would otherwise require millions of numerical relativity waveforms, such as parameter inference and tests of general relativity with GW observations.Comment: 10 pages, 5 figures; Added report numbe

On Estimation of the Post-Newtonian Parameters in the Gravitational-Wave Emission of a Coalescing Binary

The effect of the recently obtained 2nd post-Newtonian corrections on the accuracy of estimation of parameters of the gravitational-wave signal from a coalescing binary is investigated. It is shown that addition of this correction degrades considerably the accuracy of determination of individual masses of the members of the binary. However the chirp mass and the time parameter in the signal is still determined to a very good accuracy. The possibility of estimation of effects of other theories of gravity is investigated. The performance of the Newtonian filter is investigated and it is compared with performance of post-Newtonian search templates introduced recently. It is shown that both search templates can extract accurately useful information about the binary.Comment: 34 pages, 118Kb, LATEX format, submitted to Phys. Rev.

Einstein boundary conditions for the 3+1 Einstein equations

In the 3+1 framework of the Einstein equations for the case of vanishing shift vector and arbitrary lapse, we calculate explicitly the four boundary equations arising from the vanishing of the projection of the Einstein tensor along the normal to the boundary surface of the initial-boundary value problem. Such conditions take the form of evolution equations along (as opposed to across) the boundary for certain components of the extrinsic curvature and for certain space-derivatives of the intrinsic metric. We argue that, in general, such boundary conditions do not follow necessarily from the evolution equations and the initial data, but need to be imposed on the boundary values of the fundamental variables. Using the Einstein-Christoffel formulation, which is strongly hyperbolic, we show how three of the boundary equations should be used to prescribe the values of some incoming characteristic fields. Additionally, we show that the fourth one imposes conditions on some outgoing fields.Comment: Revtex 4, 6 pages, text and references added, typos corrected, to appear in Phys. Rev.

Geometrical optics analysis of the short-time stability properties of the Einstein evolution equations

Many alternative formulations of Einstein's evolution have lately been examined, in an effort to discover one which yields slow growth of constraint-violating errors. In this paper, rather than directly search for well-behaved formulations, we instead develop analytic tools to discover which formulations are particularly ill-behaved. Specifically, we examine the growth of approximate (geometric-optics) solutions, studied only in the future domain of dependence of the initial data slice (e.g. we study transients). By evaluating the amplification of transients a given formulation will produce, we may therefore eliminate from consideration the most pathological formulations (e.g. those with numerically-unacceptable amplification). This technique has the potential to provide surprisingly tight constraints on the set of formulations one can safely apply. To illustrate the application of these techniques to practical examples, we apply our technique to the 2-parameter family of evolution equations proposed by Kidder, Scheel, and Teukolsky, focusing in particular on flat space (in Rindler coordinates) and Schwarzchild (in Painleve-Gullstrand coordinates).Comment: Submitted to Phys. Rev.