No Arabic abstract
While most binary inspirals are expected to have circularized before they enter the LIGO/Virgo frequency band, a small fraction of those binaries could have non-negligible orbital eccentricity depending on their formation channel. Hence, it is important to accurately model eccentricity effects in waveform models used to detect those binaries, infer their properties, and shed light on their astrophysical environment. We develop a multipolar effective-one-body (EOB) eccentric waveform model for compact binaries whose components have spins aligned or anti-aligned with the orbital angular momentum. The waveform model contains eccentricity effects in the radiation-reaction force and gravitational modes through second post-Newtonian (PN) order, including tail effects, and spin-orbit and spin-spin couplings. We recast the PN-expanded, eccentric radiation-reaction force and modes in factorized form so that the newly derived terms can be directly included in the state-of-the-art, quasi-circular--orbit EOB model currently used in LIGO/Virgo analyses (i.e., the SEOBNRv4HM model).
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 present a frequency domain reduced order model (ROM) for the aligned-spin effective-one-body (EOB) model for binary black holes (BBHs) SEOBNRv4HM that includes the spherical harmonics modes $(ell, |m|) = (2,1),(3,3),(4,4),(5,5)$ beyond the dominant $(ell, |m|) = (2,2)$ mode. These higher modes are crucial to accurately represent the waveform emitted from asymmetric BBHs. We discuss a decomposition of the waveform, extending other methods in the literature, that allows us to accurately and efficiently capture the morphology of higher mode waveforms. We show that the ROM is very accurate with median (maximum) values of the unfaithfulness against SEOBNRv4HM lower than $0.001% (0.03%)$ for total masses in $[2.8,100] M_odot$. For a total mass of $M = 300 M_odot$ the median (maximum) value of the unfaithfulness increases up to $0.004% (0.17%)$. This is still at least an order of magnitude lower than the estimated accuracy of SEOBNRv4HM compared to numerical relativity simulations. The ROM is two orders of magnitude faster in generating a waveform compared to SEOBNRv4HM. Data analysis applications typically require $mathcal{O}(10^6-10^8)$ waveform evaluations for which SEOBNRv4HM is in general too slow. The ROM is therefore crucial to allow the SEOBNRv4HM waveform to be used in searches and Bayesian parameter inference. We present a targeted parameter estimation study that shows the improvements in measuring binary parameters when using waveforms that includes higher modes and compare against three other waveform models.
Spinning neutron stars acquire a quadrupole moment due to their own rotation. This quadratic-in-spin, self-spin effect depends on the equation of state (EOS) and affects the orbital motion and rate of inspiral of neutron star binaries. We incorporate the EOS-dependent self-spin (or monopole-quadrupole) terms in the spin-aligned effective-one-body (EOB) waveform model TEOBResumS at next-to-next-to-leading (NNLO) order, together with other (bilinear, cubic and quartic) nonlinear-in-spin effects (at leading order, LO). The structure of the Hamiltonian of TEOBResumS is such that it already incorporates, in the binary black hole case, the recently computed quartic-in-spin LO term. Using the gauge-invariant characterization of the phasing provided by the function $Q_omega=omega^2/dot{omega}$ of $omega=2pi f$ , where $f$ is the gravitational wave frequency, we study the EOS dependence of the self-spin effects and show that: (i) the next-to-leading order (NLO) and NNLO monopole-quadrupole corrections yield increasingly phase-accelerating effects compared to the corresponding LO contribution; (ii) the standard TaylorF2 post-Newtonian (PN) treatment of NLO (3PN) EOS-dependent self-spin effects makes their action stronger than the corresponding EOB description; (iii) the addition to the standard 3PN TaylorF2 post-Newtonian phasing description of self-spin tail effects at LO allows one to reconcile the self-spin part of the TaylorF2 PN phasing with the corresponding TEOBResumS one up to dimensionless frequencies $Momegasimeq 0.04-0.06$. By generating the inspiral dynamics using the post-adiabatic approximation, incorporated in a new implementation of TEOBResumS, one finds that the computational time needed to obtain a typical waveform (including all multipoles up to $ell=8$) from 10 Hz is of the order of 0.4 sec.
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.
We present TEOBResumS, a new effective-one-body (EOB) waveform model for nonprecessing (spin-aligned) and tidally interacting compact binaries.Spin-orbit and spin-spin effects are blended together by making use of the concept of centrifugal EOB radius. The point-mass sector through merger and ringdown is informed by numerical relativity (NR) simulations of binary black holes (BBH) computed with the SpEC and BAM codes. An improved, NR-based phenomenological description of the postmerger waveform is developed.The tidal sector of TEOBResumS describes the dynamics of neutron star binaries up to merger and incorporates a resummed attractive potential motivated by recent advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. Equation-of-state dependent self-spin interactions (monopole-quadrupole effects) are incorporated in the model using leading-order post-Newtonian results in a new expression of the centrifugal radius. TEOBResumS is compared to 135 SpEC and 19 BAM BBH waveforms. The maximum unfaithfulness to SpEC data $bar{F}$ -- at design Advanced-LIGO sensitivity and evaluated with total mass $M$ varying between $10M_odot leq M leq 200 M_odot$ --is always below $2.5 times 10^{-3}$ except for a single outlier that grazes the $7.1 times 10^{-3}$ level. When compared to BAM data, $bar{F}$ is smaller than $0.01$ except for a single outlier in one of the corners of the NR-covered parameter space, that reaches the $0.052$ level.TEOBResumS is also compatible, up to merger, to high end NR waveforms from binary neutron stars with spin effects and reduced initial eccentricity computed with the BAM and THC codes. The model is designed to generate accurate templates for the analysis of LIGO-Virgo data through merger and ringdown. We demonstrate its use by analyzing the publicly available data for GW150914.