Search

We extract gravitational waveforms from numerical simulations of black hole binaries computed using the Spectral Einstein Code. We compare two extraction methods: direct construction of the Newman-Penrose (NP) scalar ψ4 at a finite distance from the source and Cauchy-characteristic extraction (CCE). The direct NP approach is simpler than CCE, but NP waveforms can be contaminated by near-zone effects--unless the waves are extracted at several distances from the source and extrapolated to infinity. Even then, the resulting waveforms can in principle be contaminated by gauge effects. In contrast, CCE directly provides, by construction, gauge-invariant waveforms at future null infinity. We verify the gauge invariance of CCE by running the same physical simulation using two different gauge conditions. We find that these two gauge conditions produce the same CCE waveforms but show differences in extrapolated-ψ4 waveforms. We examine data from several different binary configurations and measure the dominant sources of error in the extrapolated-ψ4 and CCE waveforms. In some cases, we find that NP waveforms extrapolated to infinity agree with the corresponding CCE waveforms to within the estimated error bars. However, we find that in other cases extrapolated and CCE waveforms disagree, most notably for m = 0 "memory" modes. [ABSTRACT FROM AUTHOR]

This paper develops a method for solving Einstein's equation numerically on multicube representations of manifolds with arbitrary spatial topologies. This method is designed to provide a set of flexible, easy to use computational procedures that make it possible to explore the never before studied properties of solutions to Einstein's equation on manifolds with arbitrary toplogical structures. A new covariant, first-order symmetric-hyperbolic representation of Einstein's equation is developed for this purpose, along with the needed boundary conditions at the interfaces between adjoining cubic regions. Numerical tests are presented that demonstrate the long-term numerical stability of this method for evolutions of a complicated, time-dependent solution of Einstein's equation coupled to a complex scalar field on a manifold with spatial topology S³. The accuracy of these numerical test solutions is evaluated by performing convergence studies and by comparing the full nonlinear numerical results to the analytical perturbation solutions, which are also derived here. [ABSTRACT FROM AUTHOR]

Gravitational waves emitted by black-hole binary systems have the highest signal-to-noise ratio in LIGO and Virgo detectors when black-hole spins are aligned with the orbital angular momentum and extremal. For such systems, we extend the effective-one-body inspiral-merger-ringdown waveforms to generic mass ratios and spins calibrating them to 38 numerical-relativity nonprecessing waveforms produced by the SXS Collaboration. The numerical-relativity simulations span mass ratios from 1 to 8, spin magnitudes up to 98% of extremality, and last for 40 to 60 gravitational-wave cycles. When the total mass of the binary is between 20 and 200AiQ, the effective-one-body nonprecessing (dominant mode) waveforms have overlap above 99% (using the advanced-LIGO design noise spectral density) with all of the 38 nonprecessing numerical waveforms, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling. We also show that--without further calibration-- the precessing effective-one-body (dominant mode) waveforms have overlap above 97% with two very long, strongly precessing numerical-relativity waveforms, when maximizing only on the initial phase and time. [ABSTRACT FROM AUTHOR]

We study the general relativistic periastron advance in spinning black hole binaries on quasicircular orbits, with spins aligned or antialigned with the orbital angular momentum, using numerical-relativity simulations, the post-Newtonian approximation, and black hole perturbation theory. By imposing a symmetry by exchange of the bodies' labels, we devise an improved version of the perturbative result and use it as the leading term of a new type of expansion in powers of the symmetric mass ratio. This allows us to measure, for the first time, the gravitational self-force effect on the periastron advance of a nonspinning particle orbiting a Kerr black hole of mass M and spin S=-0.5M², down to separations of order 9M. Comparing the predictions of our improved perturbative expansion with the exact results from numerical simulations of equal-mass and equal-spin binaries, we find a remarkable agreement over a wide range of spins and orbital separations. [ABSTRACT FROM AUTHOR]

We compute the periastron advance using the effective-one-body formalism for binary black holes moving on quasicircular orbits and having spins collinear with the orbital angular momentum. We compare the predictions with the periastron advance recently computed in accurate numerical-relativity simulations and find remarkable agreement for a wide range of spins and mass ratios. These results do not use any numerical-relativity calibration of the effective-one-body model, and stem from two key ingredients in the effective-one-body Hamiltonian: (i) the mapping of the two-body dynamics of spinning particles onto the dynamics of an effective spinning particle in a (deformed) Kerr spacetime, fully symmetrized with respect to the two-body masses and spins, and (ii) the resummation, in the test-particle limit, of all post-Newtonian corrections linear in the spin of the particle. In fact, even when only the leading spin post-Newtonian corrections are included in the effective-one-body spinning Hamiltonian but all the test-particle corrections linear in the spin of the particle are resummed we find very good agreement with the numerical results (within the numerical error for equal-mass binaries and discrepancies of at most I % for larger mass ratios). Furthermore, we specialize to the extreme mass-ratio limit and derive, using the equations of motion in the gravitational skeleton approach, analytical expressions for the periastron advance, the meridional Lense-Thirring precession and spin precession frequency in the case of a spinning particle on a nearly circular equatorial orbit in Kerr spacetime, including also terms quadratic in the spin. [ABSTRACT FROM AUTHOR]

The behavior of merging black holes (including the emitted gravitational waves and the properties of the remnant) can currently be computed only by numerical simulations. This paper introduces ten numerical relativity simulations of binary black holes with equal masses and equal spins aligned or antialigned with the orbital angular momentum. The initial spin magnitudes have |Xi| ≤ 0.95 and are more concentrated in the aligned direction because of the greater astrophysical interest of this case. We combine these data with five previously reported simulations of the same configuration, but with different spin magnitudes, including the highest spin simulated to date, Xi ≈ 0.97. This data set is sufficiently accurate to enable us to offer improved analytic fitting formulas for the final spin and for the energy radiated by gravitational waves as a function of initial spin. The improved fitting formulas can help to improve our understanding of the properties of binary black hole merger remnants and can be used to enhance future approximate waveforms for gravitational wave searches, such as effective-one-body waveforms. [ABSTRACT FROM AUTHOR]

This article studies sufficient accuracy criteria of hybrid post-Newtonian (PN) and numerical relativity (NR) waveforms for parameter estimation of strong binary black-hole sources in second-generation ground-based gravitational-wave detectors. We investigate equal-mass nonspinning binaries with a new 33-orbit NR waveform, as well as unequal-mass binaries with mass ratios 2, 3, 4 and 6. For equal masses, the 33-orbit NR waveform allows us to recover previous results and to extend the analysis toward matching at lower frequencies. For unequal masses, the errors between different PN approximants increase with mass ratio. Thus, at 3.5 PN, hybrids for higher-mass-ratio systems would require NR waveforms with many more gravitational-wave cycles to guarantee no adverse impact on parameter estimation. Furthermore, we investigate the potential improvement in hybrid waveforms that can be expected from fourth-order post-Newtonian waveforms and find that knowledge of this fourth post- Newtonian order would significantly improve the accuracy of hybrid waveforms. [ABSTRACT FROM AUTHOR]