No Arabic abstract
We present the first numerical-relativity simulation of a compact-object binary whose gravitational waveform is long enough to cover the entire frequency band of advanced gravitational-wave detectors, such as LIGO, Virgo and KAGRA, for mass ratio 7 and total mass as low as $45.5,M_odot$. We find that effective-one-body models, either uncalibrated or calibrated against substantially shorter numerical-relativity waveforms at smaller mass ratios, reproduce our new waveform remarkably well, with a negligible loss in detection rate due to modeling error. In contrast, post-Newtonian inspiral waveforms and existing calibrated phenomenological inspiral-merger-ringdown waveforms display greater disagreement with our new simulation. The disagreement varies substantially depending on the specific post-Newtonian approximant used.
Gravitational waveforms which describe the inspiral, merger and ringdown of coalescing binaries are usually constructed by synthesising information from perturbative descriptions, in particular post-Newtonian theory and black-hole perturbation theory, with numerical solutions of the full Einstein equations. In this paper we discuss the glueing of numerical and post-Newtonian waveforms to produce hybrid waveforms which include subdominant spherical harmonics (higher order modes), and focus in particular on the process of consistently aligning the waveforms, which requires a comparison of both descriptions and a discussion of their imprecisions. We restrict to the non-precessing case, and illustrate the process using numerical waveforms of up to mass ratio $q=18$ produced with the BAM code, and publicly available waveforms from the SXS catalogue. The results also suggest new ways of analysing finite radius errors in numerical simulations.
Numerical simulations of 15 orbits of an equal-mass binary black hole system are presented. Gravitational waveforms from these simulations, covering more than 30 cycles and ending about 1.5 cycles before merger, are compared with those from quasi-circular zero-spin post-Newtonian (PN) formulae. The cumulative phase uncertainty of these comparisons is about 0.05 radians, dominated by effects arising from the small residual spins of the black holes and the small residual orbital eccentricity in the simulations. Matching numerical results to PN waveforms early in the run yields excellent agreement (within 0.05 radians) over the first $sim 15$ cycles, thus validating the numerical simulation and establishing a regime where PN theory is accurate. In the last 15 cycles to merger, however, {em generic} time-domain Taylor approximants build up phase differences of several radians. But, apparently by coincidence, one specific post-Newtonian approximant, TaylorT4 at 3.5PN order, agrees much better with the numerical simulations, with accumulated phase differences of less than 0.05 radians over the 30-cycle waveform. Gravitational-wave amplitude comparisons are also done between numerical simulations and post-Newtonian, and the agreement depends on the post-Newtonian order of the amplitude expansion: the amplitude difference is about 6--7% for zeroth order and becomes smaller for increasing order. A newly derived 3.0PN amplitude correction improves agreement significantly ($<1%$ amplitude difference throughout most of the run, increasing to 4% near merger) over the previously known 2.5PN amplitude terms.
The Christodoulou memory is a nonlinear contribution to the gravitational-wave field that is sourced by the gravitational-wave stress-energy tensor. For quasicircular, inspiralling binaries, the Christodoulou memory produces a growing, nonoscillatory change in the gravitational-wave plus polarization, resulting in the permanent displacement of a pair of freely-falling test masses after the wave has passed. In addition to its nonoscillatory behavior, the Christodoulou memory is interesting because even though it originates from 2.5 post-Newtonian (PN) order multipole interactions, it affects the waveform at leading (Newtonian/quadrupole) order. The memory is also potentially detectable in binary black-hole mergers. While the oscillatory pieces of the gravitational-wave polarizations for quasicircular, inspiralling compact binaries have been computed to 3PN order, the memory contribution to the polarizations has only been calculated to leading order (the next-to-leading order 0.5PN term has previously been shown to vanish). Here the calculation of the memory for quasicircular, inspiralling binaries is extended to 3PN order. While the angular dependence of the memory remains qualitatively unchanged, the PN correction terms tend to reduce the memorys magnitude. Explicit expressions are given for the memory contributions to the plus polarization and the spin-weighted spherical-harmonic modes of the metric and curvature perturbations. Combined with the recent results of Blanchet et al.(2008), this completes the waveform to 3PN order. This paper also discusses: (i) the difficulties in extracting the memory from numerical simulations, (ii) other nonoscillatory effects that enter the waveform at high PN orders, and (iii) issues concerning the observability of the memory.
In the adiabatic post-Newtonian (PN) approximation, the phase evolution of gravitational waves (GWs) from inspiralling compact binaries in quasicircular orbits is computed by equating the change in binding energy with the GW flux. This energy balance equation can be solved in different ways, which result in multiple approximants of the PN waveforms. Due to the poor convergence of the PN expansion, these approximants tend to differ from each other during the late inspiral. Which of these approximants should be chosen as templates for detection and parameter estimation of GWs from inspiraling compact binaries is not obvious. In this paper, we present estimates of the effective higher order (beyond the currently available 4PN and 3.5PN) non-spinning terms in the PN expansion of the binding energy and the GW flux that minimize the difference of multiple PN approximants (TaylorT1, TaylorT2, TaylorT4, TaylorF2) with effective one body waveforms calibrated to numerical relativity (EOBNR). We show that PN approximants constructed using the effective higher order terms show significantly better agreement (as compared to 3.5PN) with the inspiral part of the EOBNR. For non-spinning binaries with component masses $m_{1,2} in [1.4 M_odot, 15 M_odot]$, most of the approximants have a match (faithfulness) of better than 99% with both EOBNR and each other.
We estimate the probability of detecting a gravitational wave signal from coalescing compact binaries in simulated data from a ground-based interferometer detector of gravitational radiation using Bayesian model selection. The simulated waveform of the chirp signal is assumed to be a spin-less Post-Newtonian (PN) waveform of a given expansion order, while the searching template is assumed to be either of the same Post-Newtonian family as the simulated signal or one level below its Post-Newtonian expansion order. Within the Bayesian framework, and by applying a reversible jump Markov chain Monte Carlo simulation algorithm, we compare PN1.5 vs. PN2.0 and PN3.0 vs. PN3.5 wave forms by deriving the detection probabilities, the statistical uncertainties due to noise as a function of the SNR, and the posterior distributions of the parameters. Our analysis indicates that the detection probabilities are not compromised when simplified models are used for the comparison, while the accuracies in the determination of the parameters characterizing these signals can be significantly worsened, no matter what the considered Post-Newtonian order expansion comparison is.