No Arabic abstract
We construct a new factorized waveform including $(l,|m|)=(2,2),(2,1),(3,3),(4,4)$ modes based on effective-one-body (EOB) formalism, which is valid for spinning binary black holes (BBH) in general equatorial orbit. When combined with the dynamics of $texttt{SEOBNRv4}$, the $(l,|m|)=(2,2)$ mode waveform generated by this new waveform can fit the original $texttt{SEOBNRv4}$ waveform very well in the case of a quasi-circular orbit. We have calibrated our new waveform model to the Simulating eXtreme Spacetimes (SXS) catalog. The comparison is done for BBH with total mass in $(20,200)M_odot$ using Advanced LIGO designed sensitivity. For the quasi-circular cases we have compared our $(2,2)$ mode waveforms to the 281 numerical relativity (NR) simulations of BBH along quasi-circular orbits. All of the matching factors are bigger than 98%. For the elliptical cases, 24 numerical relativity simulations of BBH along an elliptic orbit are used. For each elliptical BBH system, we compare our modeled gravitational polarizations against the NR results for different combinations of the inclination angle, the initial orbit phase and the source localization in the sky. We use the the minimal matching factor respect to the inclination angle, the initial orbit phase and the source localization to quantify the performance of the higher modes waveform. We found that after introducing the high modes, the minimum of the minimal matching factor among the 24 tested elliptical BBHs increases from 90% to 98%. Following our previous $texttt{SEOBNRE}$ waveform model, we call our new waveform model $texttt{SEOBNREHM}$. Our $texttt{SEOBNREHM}$ waveform model can match all tested 305 SXS waveforms better than 98% including highly spinning ($chi=0.99$) BBH, highly eccentric ($eapprox0.15$) BBH and large mass ratio ($q=10$) BBH.
Gravitational wave astrophysics relies heavily on the use of matched filtering both to detect signals in noisy data from detectors, and to perform parameter estimation on those signals. Matched filtering relies upon prior knowledge of the signals expected to be produced by a range of astrophysical systems, such as binary black holes. These waveform signals can be computed using numerical relativity techniques, where the Einstein field equations are solved numerically, and the signal is extracted from the simulation. Numerical relativity simulations are, however, computationally expensive, leading to the need for a surrogate model which can predict waveform signals in regions of the physical parameter space which have not been probed directly by simulation. We present a method for producing such a surrogate using Gaussian process regression which is trained directly on waveforms generated by numerical relativity. This model returns not just a single interpolated value for the waveform at a new point, but a full posterior probability distribution on the predicted value. This model is therefore an ideal component in a Bayesian analysis framework, through which the uncertainty in the interpolation can be taken into account when performing parameter estimation of signals.
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 Einsteins 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.
We develop the foundations of an effective-one-body (EOB) model for eccentric binary coalescences that includes the conservative dynamics, radiation reaction, and gravitational waveform modes from the inspiral and the merger-ringdown signals. We use the same approach as is commonly employed in black-hole perturbation theory by introducing a relativistic parameterization of the dynamics that is defined by the orbital geometry and consists of a set of phase variables and quantities that evolve only due to gravitational radiation reaction. Specializing to nonspinning binaries, we derive the EOB evolution equations and compute the binarys radiative multipole moments that determine the gravitational waves through a decomposition into the fundamental frequencies of the motion. The major differences between our treatment and the quasi-Keplerian approach often used in post-Newtonian (PN) calculations are that the orbital parameters describe strong-field dynamics, and that expressing the multipole moments in terms of the frequencies simplifies the calculations and also results in an unambiguous orbit-averaging operation. While our description of the conservative dynamics is fully relativistic, we limit explicit derivations in the radiative sector to 1.5PN order for simplicity. This already enables us to establish methods for computing both instantaneous and hereditary contributions to the gravitational radiation in EOB coordinates that have straightforward extensions to higher PN order. The weak-field, small eccentricity limit of our results for the orbit-averaged fluxes of energy and angular momentum agrees with known PN results when expressed in terms of gauge-invariant quantities. We further address considerations for the numerical implementation of the model and the completion of the waveforms to include the merger and ringdown signals, and provide illustrative results.
As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binarys orbital plane to precess. To accurately model such systems, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the inspiral, merger and ringdown (IMR) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to new 118 precessing NR waveforms, which span mass ratios 1-4 and spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for $94%$ ($57%$) of the cases, the maximum value, in the total mass range $20-200 M_odot$, is below $3%$ ($1%$). Those numbers become $83%$ ($20%$) when using the IMR, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models and identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM.
We compute the periastron advance using the effective-one-body formalism for binary black holes moving on quasi-circular 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 (PN) corrections linear in the spin of the particle. In fact, even when only the leading spin PN 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 1% 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.