No Arabic abstract
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.
Gravitational-wave observations of binary neutron star systems can provide information about the masses, spins, and structure of neutron stars. However, this requires accurate and computationally efficient waveform models that take <1s to evaluate for use in Bayesian parameter estimation codes that perform 10^7 - 10^8 waveform evaluations. We present a surrogate model of a nonspinning effective-one-body waveform model with l = 2, 3, and 4 tidal multipole moments that reproduces waveforms of binary neutron star numerical simulations up to merger. The surrogate is built from compact sets of effective-one-body waveform amplitude and phase data that each form a reduced basis. We find that 12 amplitude and 7 phase basis elements are sufficient to reconstruct any binary neutron star waveform with a starting frequency of 10Hz. The surrogate has maximum errors of 3.8% in amplitude (0.04% excluding the last 100M before merger) and 0.043 radians in phase. The version implemented in the LIGO Algorithm Library takes ~0.07s to evaluate for a starting frequency of 30Hz and ~0.8s for a starting frequency of 10Hz, resulting in a speed-up factor of ~10^3 - 10^4 relative to the original Matlab code. This allows parameter estimation codes to run in days to weeks rather than years, and we demonstrate this with a Nested Sampling run that recovers the masses and tidal parameters of a simulated binary neutron star system.
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).
We discuss the properties of the effective-one-body (EOB) multipolar gravitational waveform emitted by nonspinning black-hole binaries of masses $mu$ and $M$ in the extreme-mass-ratio limit, $mu/M= ull 1$. We focus on the transition from quasicircular inspiral to plunge, merger and ringdown.We compare the EOB waveform to a Regge-Wheeler-Zerilli (RWZ) waveform computed using the hyperboloidal layer method and extracted at null infinity. Because the EOB waveform keeps track analytically of most phase differences in the early inspiral, we do not allow for any arbitrary time or phase shift between the waveforms. The dynamics of the particle, common to both wave-generation formalisms, is driven by leading-order ${cal O}( u)$ analytically--resummed radiation reaction. The EOB and the RWZ waveforms have an initial dephasing of about $5times 10^{-4}$ rad and maintain then a remarkably accurate phase coherence during the long inspiral ($sim 33$ orbits), accumulating only about $-2times 10^{-3}$ rad until the last stable orbit, i.e. $Deltaphi/phisim -5.95times 10^{-6}$. We obtain such accuracy without calibrating the analytically-resummed EOB waveform to numerical data, which indicates the aptitude of the EOB waveform for LISA-oriented studies. We then improve the behavior of the EOB waveform around merger by introducing and tuning next-to-quasi-circular corrections both in the gravitational wave amplitude and phase. For each multipole we tune only four next-to-quasi-circular parameters by requiring compatibility between EOB and RWZ waveforms at the light-ring. The resulting phase difference around merger time is as small as $pm 0.015$ rad, with a fractional amplitude agreement of 2.5%. This suggest that next-to-quasi-circular corrections to the phase can be a useful ingredient in comparisons between EOB and numerical relativity waveforms.
The use of effective-one-body (EOB) waveforms for black hole binaries analysis in gravitational-wave astronomy requires faithful models and fast generation times. A key aspect to achieve faithfulness is the inclusion of numerical-relativity (NR) informed next-to-quasicircular corrections(NQC), dependent on the radial momentum, to the waveform and radiation reaction. A robust method to speed up the waveform generation is the post-adiabatic iteration to approximate the solution of the EOB Hamiltonian equations. In this work, we assess the performances of a fast NQC prescription in combination to the post-adiabatic method for generating multipolar gravitational waves. The outlined approach allows a consistent treatment of NQC in both the waveform and the radiation-reaction, does not require iterative procedures to achieve high faithfulness, and can be efficiently employed for parameter estimation. Comparing to 611 NR simulations, for total mass $10M_odotleq M leq 200M_odot$ and using the Advanded LIGO noise, the model has EOB/NR unfaithfulness well below $0.01$, with 78.5% of the cases below $0.001$. We apply the model to the parameter estimation of GW150914 exploring the impact of the new NQC and of the higher modes up to $ell=m=8$.
After eleven gravitational-wave detections from compact-binary mergers, we are yet to observe the striking general-relativistic phenomenon of orbital precession. Measurements of precession would provide valuable insights into the distribution of black-hole spins, and therefore into astrophysical binary formation mechanisms. Using our recent two-harmonic approximation of precessing-binary signals~cite{Fairhurst:2019_2harm}, we introduce the ``precession signal-to-noise ratio, $rho_p$. We demonstrate that this can be used to clearly identify whether precession was measured in an observation (by comparison with both current detections and simulated signals), and can immediately quantify the measurability of precession in a given signal, which currently requires computationally expensive parameter-estimation studies. $rho_p$ has numerous potential applications to signal searches, source-property measurements, and population studies. We give one example: assuming one possible astrophysical spin distribution, we predict that precession has a one in $sim 25$ chance of being observed in any detection.