Effective-one-body multipolar waveforms for eccentric binary black holes with nonprecessing spins

We construct an inspiral-merger-ringdown eccentric gravitational-wave (GW) model for binary black holes with non-precessing spins within the effective-one-body formalism. This waveform model, SEOBNRv4EHM, extends the accurate quasi-circular SEOBNRv4HM model to eccentric binaries by including recently computed eccentric corrections up to 2PN order in the gravitational waveform modes, notably the $(l,|m|)=(2,2),(2,1),(3,3),(4,4),(5,5)$ multipoles. The waveform model reproduces the zero eccentricity limit with an accuracy comparable to the underlying quasi-circular model, with the unfaithfulness of $\lesssim1\%$ against quasi-circular numerical-relativity (NR) simulations. When compared against 28 public eccentric NR simulations from the Simulating eXtreme Spacetimes catalog with initial orbital eccentricities up to $e\simeq0.3$ and dimensionless spin magnitudes up to $+0.7$, the model provides unfaithfulness $<1\%$, showing that both the $(2,|2|)$-modes and the higher-order modes are reliably described without calibration to NR datasets in the eccentric sector. The waveform model SEOBNRv4EHM is able to qualitatively reproduce the phenomenology of dynamical captures, and can be extended to include spin-precession effects. It can be employed for upcoming observing runs with the LIGO-Virgo-KAGRA detectors and used to re-analyze existing GW catalogs to infer the eccentricity parameters for binaries with $e\lesssim0.3$ (at 20 Hz or lower) and spins up to $\lesssim 0.9-0.95$. The latter is a promising region of the parameter space where some astrophysical formation scenarios of binaries predict mild eccentricity in the ground-based detectors' bandwidth. Assessing the accuracy and robustness of the eccentric waveform model SEOBNRv4EHM for larger eccentricities and spins will require comparisons with, and, likely, calibration to eccentric NR waveforms in a larger region of the parameter space

Initial Data and Eccentricity Reduction Toolkit for Binary Black Hole Numerical Relativity Waveforms

The production of numerical relativity waveforms that describe quasicircular binary black hole mergers requires high-quality initial data, and an algorithm to iteratively reduce residual eccentricity. To date, these tools remain closed source, or in commercial software that prevents their use in high performance computing platforms. To address these limitations, and to ensure that the broader numerical relativity community has access to these tools, herein we provide all the required elements to produce high-quality numerical relativity simulations in supercomputer platforms, namely: open source parameter files to numerical simulate spinning black hole binaries with asymmetric mass-ratios; open source $\texttt{Python}$ tools to produce high-quality initial data for numerical relativity simulations of spinning black hole binaries on quasi-circular orbits; open source $\texttt{Python}$ tools for eccentricity reduction, both as stand-alone software and deployed in the $\texttt{Einstein Toolkit}$'s software infrastructure. This open source toolkit fills in a critical void in the literature at a time when numerical relativity has an ever increasing role in the study and interpretation of gravitational wave sources. As part of our community building efforts, and to streamline and accelerate the use of these resources, we provide tutorials that describe, step by step, how to obtain and use these open source numerical relativity tools

Time-domain phenomenological model of gravitational-wave subdominant harmonics for quasicircular nonprecessing binary black hole coalescences

In this work we present an extension of the time domain phenomenological model IMRPhenomT for gravitational wave signals from binary black hole coalescences to include subdominant harmonics, specifically the $(l=2, m=\pm 1)$, $(l=3, m=\pm 3)$, $(l=4, m=\pm 4)$ and $(l=5, m=\pm 5)$ spherical harmonics. We also improve our model for the dominant $(l=2, m=\pm 2)$ mode and discuss mode mixing for the $(l=3, m=\pm 2)$ mode. The model is calibrated to numerical relativity solutions of the full Einstein equations up to mass ratio 18, and to numerical solutions of the Teukolsky equations for higher mass ratios. This work complements the latest generation of traditional frequency domain phenomenological models (IMRPhenomX), and provides new avenues to develop computationally efficient models for gravitational wave signals from generic compact binaries

Eccentric binary black holes: Comparing numerical relativity and small mass-ratio perturbation theory

The modelling of unequal mass binary black hole systems is of high importanceto detect and estimate parameters from these systems. Numerical relativity (NR)is well suited to study systems with comparable component masses, m_1\simm_2, whereas small mass ratio (SMR) perturbation theory applies to binarieswhere $q=m_2/m_1as a function of mass ratio for eccentric non-spinning binary black holes. Weproduce$52$NR simulations with mass ratios between$1:10$and$1:1$andinitial eccentricities up to$0.7$. From these we extract quantities likegravitational wave energy and angular momentum fluxes and periastron advance,and assess their accuracy. To facilitate comparison, we develop tools to mapbetween NR and SMR inspiral evolutions of eccentric binary black holes. Wederive post-Newtonian accurate relations between different definitions ofeccentricity. Based on these analyses, we introduce a new definition ofeccentricity based on the (2,2)-mode of the gravitational radiation, whichreduces to the Newtonian definition of eccentricity in the Newtonian limit.From the comparison between NR simulations and SMR results, we quantify theunknown next-to-leading order SMR contributions to the gravitational energy andangular momentum fluxes, and periastron advance. We show that in the comparablemass regime these contributions are subdominant and higher order SMRcontributions are negligible.<br

Parameter estimation with the current generation of phenomenological waveform models applied to the black hole mergers of GWTC-1

We consider the ten confidently detected gravitational-wave signals in theGWTC-1 catalog which are consistent with mergers of binary black hole systems,and perform a thorough parameter estimation re-analysis. This is made possibleby using computationally efficient waveform models of the current (fourth)generation of the IMRPhenom family of phenomenological waveform models, whichconsists of the IMRPhenomX frequency-domain modelsand the IMRPhenomTtime-domain models. The analysis is performed with both precessing andnon-precessing waveform models with and without subdominant spherical harmonicmodes. Results for all events are validated with convergence tests, discussingin particular the events GW170729 and GW151226. For the latter and the othertwo lowest-mass events, we also compare results between two independentsampling codes, Bilby and LALInference. We find overall consistent results withthe original GWTC-1 results, with all Jensen-Shannon divergences between theprevious results using IMRPhenomPv2 and our default IMRPhenomXPHM posteriorsbelow 0.045 bits, but we also discuss cases where including subdominantharmonics and/or precession influences the posteriors.<br

New twists in compact binary waveform modeling: A fast time-domain model for precession

We present IMRPhenomTPHM, a phenomenological model for the gravitational wave signals emitted by the coalescence of quasi-circular precessing binary black holes systems. The model is based on the "twisting up" approximation, which maps non-precessing signals to precessing ones in terms of a time dependent rotation described by three Euler angles, and which has been utilized in several frequency domain waveform models that have become standard tools in gravitational wave data analysis. Our model is however constructed in the time domain, which allows several improvements over the frequency domain models: we do not use the stationary phase approximation, we employ a simple approximation for the precessing Euler angles for the ringdown signal, and we implement a new method for computing the Euler angles through the evolution of the spin dynamics of the system, which is more accurate and also computationally efficient

Eccentric binary black holes: Comparing numerical relativity and small mass-ratio perturbation theory

The modelling of unequal mass binary black hole systems is of high importance to detect and estimate parameters from these systems. Numerical relativity (NR) is well suited to study systems with comparable component masses, $m_1\sim m_2$, whereas small mass ratio (SMR) perturbation theory applies to binaries where $q=m_2/m_1<< 1$. This work investigates the applicability for NR and SMR as a function of mass ratio for eccentric non-spinning binary black holes. We produce $52$ NR simulations with mass ratios between $1:10$ and $1:1$ and initial eccentricities up to $0.7$. From these we extract quantities like gravitational wave energy and angular momentum fluxes and periastron advance, and assess their accuracy. To facilitate comparison, we develop tools to map between NR and SMR inspiral evolutions of eccentric binary black holes. We derive post-Newtonian accurate relations between different definitions of eccentricity. Based on these analyses, we introduce a new definition of eccentricity based on the (2,2)-mode of the gravitational radiation, which reduces to the Newtonian definition of eccentricity in the Newtonian limit. From the comparison between NR simulations and SMR results, we quantify the unknown next-to-leading order SMR contributions to the gravitational energy and angular momentum fluxes, and periastron advance. We show that in the comparable mass regime these contributions are subdominant and higher order SMR contributions are negligible

A detailed analysis of GW190521 with phenomenological waveform models

In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolve approximate waveform degeneracies, which may result in multi-modal posterior distributions. The family of waveform models we use includes a new fast time-domain model IMRPhenomTPHM, which allows us extensive tests of different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper [Nitz&Capano], associated with possible support for a high-mass ratio source, but confirm their finding of a multi-modal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the PISN mass-gap, and the support for an intermediate mass ratio binary are drastically reduced with respect to what Nitz&Capano found. We also provide updated probabilities for associating GW190521 to the potential electromagnetic counterpart from ZTF

Extending black-hole remnant surrogate models to extreme mass ratios

Numerical-relativity surrogate models for both black-hole merger waveforms and remnants have emerged as important tools in gravitational-wave astronomy. While producing very accurate predictions, their applicability is limited to the region of the parameter space where numerical-relativity simulations are available and computationally feasible. Notably, this excludes extreme mass ratios. We present a machine-learning approach to extend the validity of existing and future numerical-relativity surrogate models toward the test-particle limit, targeting in particular the mass and spin of post-merger black-hole remnants. Our model is trained on both numerical-relativity simulations at comparable masses and analytical predictions at extreme mass ratios. We extend the gaussian-process-regression model NRSur7dq4Remnant, validate its performance via cross validation, and test its accuracy against additional numerical-relativity runs. Our fit, which we dub NRSur7dq4EmriRemnant, reaches an accuracy that is comparable to or higher than that of existing remnant models while providing robust predictions for arbitrary mass ratios.Comment: 10 pages, 3 figures. Model publicly available at https://pypi.org/project/surfinB