Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant = m = 2 gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of ∼ 10 −4 . The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% (90 th percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms . black holes and of the surrounding spacetime [31,32]. Simulations have also been used for visualizations of curved spacetime [33][34][35][36][37][38][39][40], investigations of spin quantities [41], and the relaxation of spacetime to the Kerr solution following merger [42][43][44]. The motion of the black hole horizons and horizon curvature quantities have been used to explore eccentric dynamics [45][46][47][48], spin precession [49][50][51][52], and the first law of binary black hole mechanics [53][54][55][56][57]. These in turn have been compared to analytic post-Newtonian and self-force approximations (see also [58][59][60]), mapping out the bounds of validity of these approximations.A key application of BBH simulations is the accurate modeling of gravitational waves emitted by these systems during their late inspiral, merger, and final ringdown. Waveforms extracted from BBH simulations are essential for analyzing observed gravitational-wave signals from black hole binaries. Indeed, all BBH observations by LIGO and Virgo were analyzed using waveform families that rely on numerical relativity for their construction, most notably effective-one-body waveform models [61-65] and phenomenological waveform models [66][67][68]. Numerical simulations are also central in validating such waveform models [69][70][71][72][73][74][75][76], and were used to validate GW searches [77][78][79]. Waveforms from numerical relativity are also used directly in parameter estimation [80,81], to construct template banks [82], and to construct waveform families without intermediate analytical models, using methods such as reduced order modeling [83][84][85][86]. Today's simula...

As gravitational-wave detectors become more sensitive and broaden their frequency bandwidth, 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 different formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To accurately model such precessing signals, especially when masses and spins vary in the wide astrophysical range, 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 entire coalescence stage (i.e., inspiral, merger and ringdown) 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 black-hole 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 118 SXS precessing NR waveforms, produced as part of this project, which span mass ratios 1-4 and (dimensionless) black-hole's 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 ⊙ , is below 3% (1%). Those numbers change to 83% (20%) when using the inspiral-merger-ringdown, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two Bayesian, parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models as a function of the mass and spin parameters to 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 detailed study of the center-of-mass (c.m.) motion seen in simulations produced by the Simulating eXtreme Spacetimes (SXS) collaboration. We investigate potential physical sources for the large c.m. motion in binary black hole simulations and find that a significant fraction of the c.m. motion cannot be explained physically, thus concluding that it is largely a gauge effect. These large c.m. displacements cause mode mixing in the gravitational waveform, most easily recognized as amplitude oscillations caused by the dominant (2, ±2) modes mixing into subdominant modes. This mixing does not diminish with increasing distance from the source; it is present even in asymptotic waveforms, regardless of the method of data extraction. We describe the current c.m.-correction method used by the SXS collaboration, which is based on counteracting the motion of the c.m. as measured by the trajectories of the apparent horizons in the simulations, and investigate potential methods to improve that correction to the waveform. We also present a complementary method for computing an optimal c.m. correction or evaluating any other c.m. transformation based solely on the asymptotic waveform data.