Search

The gravitational-wave signal GW250114 was observed by the two LIGO detectors with a network matched-filter signal-to-noise ratio of 80. The signal was emitted by the coalescence of two black holes with near-equal masses $m_1={33.6}_{-0.8}^{+1.2}{M}_{\odot }$ and $m_2={32.2}_{-1.3}^{+0.8}{M}_{\odot }$, and small spins ${\chi }_{1,2}\le 0.26$ (90% credibility) and negligible eccentricity $e\le 0.03$. Postmerger data excluding the peak region are consistent with the dominant quadrupolar $(\ell =|m|=2)$ mode of a Kerr black hole and its first overtone. We constrain the modes’ frequencies to $\pm 30\%$ of the Kerr spectrum, providing a test of the remnant’s Kerr nature. We also examine Hawking’s area law, also known as the second law of black hole mechanics, which states that the total area of the black hole event horizons cannot decrease with time. A range of analyses that exclude up to five of the strongest merger cycles confirm that the remnant area is larger than the sum of the initial areas to high credibility.

The binary black hole signal GW250114, the loudest gravitational wave detected to date, offers a unique opportunity to test Einstein’s general relativity (GR) in the high-velocity, strong-gravity regime and probe whether the remnant conforms to the Kerr metric. Upon perturbation, black holes emit a spectrum of damped sinusoids with specific, complex frequencies. Our analysis of the postmerger signal shows that at least two quasinormal modes are required to explain the data, with the most damped remaining statistically significant for about one cycle. We probe the remnant’s Kerr nature by constraining the spectroscopic pattern of the dominant quadrupolar ($\ell =m=2$) mode and its first overtone to match the Kerr prediction to tens of percent at multiple postpeak times. The measured mode amplitudes and phases agree with a numerical-relativity simulation having parameters close to GW250114. By fitting a parametrized waveform that incorporates the full inspiral-merger-ringdown sequence, we constrain the fundamental $(\ell =m=4)$ mode to tens of percent and bound the quadrupolar frequency to within a few percent of the GR prediction. We perform a suite of tests—spanning inspiral, merger, and ringdown—finding constraints that are comparable to, and in some cases 2–3 times more stringent than those obtained by combining dozens of events in the fourth Gravitational-Wave Transient Catalog. These results constitute the most stringent single-event verification of GR and the Kerr nature of black holes to date, and outline the power of black-hole spectroscopy for future gravitational-wave observations.

The Bondi-van der Burg-Metzner-Sachs (BMS) group, which uniquely describes the symmetries of asymptotic infinity and therefore of the gravitational waves that propagate there, has become increasingly important for accurate modeling of waveforms. In particular, waveform models, such as post-Newtonian (PN) expressions, numerical relativity (NR), and black hole perturbation theory, produce results that are in different BMS frames. Consequently, to build a model for the waveforms produced during the merging of compact objects, which ideally would be a hybridization of PN, NR, and black hole perturbation theory, one needs a fast and robust method for fixing the BMS freedoms. In this work, we present the first means of fixing the entire BMS freedom of NR waveforms to match the frame of either PN waveforms or black hole perturbation theory. We achieve this by finding the BMS transformations that change certain charges in a prescribed way—e.g., finding the center-of-mass transformation that maps the center-of-mass charge to a mean of zero. We find that this new method is 20 times faster, and more correct when mapping to the superrest frame, than previous methods that relied on optimization algorithms. Furthermore, in the course of developing this charge-based frame fixing method, we compute the PN expression for the Moreschi supermomentum to 3PN order without spins and 2PN order with spins. This Moreschi supermomentum is effectively equivalent to the energy flux or the null memory contribution at future null infinity ${\mathcal{I}}^{+}$. From this PN calculation, we also compute oscillatory ($m\ne 0$ modes) and spin-dependent memory terms that have not been identified previously or have been missing from strain expressions in the post-Newtonian literature.

Quasinormal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing general relativity. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity strain waveforms. As gravitational wave observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically important, and in turn require a more thorough treatment of the asymptotic frame at ${\mathcal{I}}^{+}$. The first main result of this work is a method for systematically fitting a QNM model containing many modes to a numerical waveform produced using Cauchy-characteristic extraction (CCE), a waveform extraction technique which is known to resolve memory effects. We choose the modes to model based on their power contribution to the residual between numerical and model waveforms. We show that the all-mode strain mismatch improves by a factor of $\sim {10}^5$ when using multimode fitting as opposed to only fitting the $(2,\pm 2,n)$ modes. Our most significant result addresses a critical point that has been overlooked in the QNM literature: the importance of matching the Bondi-van der Burg-Metzner-Sachs (BMS) frame of the numerical waveform to that of the QNM model. We show that by mapping the numerical waveforms—which exhibit the memory effect—to a BMS frame known as the super rest frame, there is an improvement of $\sim {10}^5$ in the all-mode strain mismatch compared to using a strain waveform whose BMS frame is not fixed. Furthermore, we find that by mapping CCE waveforms to the super rest frame, we can obtain all-mode mismatches that are, on average, a factor of $\sim 4$ better than using the publicly available extrapolated waveforms. We illustrate the effectiveness of these modeling enhancements by applying them to families of waveforms produced by numerical relativity and comparing our results to previous QNM studies.

Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools have not existed. In this paper, we show how to analyze and map to a suitable BMS frame for numerical waveforms calculated with the Spectral Einstein Code (SpEC). However, the methods and tools that we present are general and can be applied to any numerical waveforms. We present an extensive study of 13 binary black hole systems that broadly span parameter space. From these simulations, we extract the strain and also the Weyl scalars using both SpECTRE’s Cauchy-characteristic extraction module and also the standard extrapolation procedure with a displacement memory correction applied during postprocessing. First, we show that the current center-of-mass correction used to map these waveforms to the center-of-mass frame is not as effective as previously thought. Consequently, we also develop an improved correction that utilizes asymptotic Poincaré charges instead of a Newtonian center-of-mass trajectory. Next, we map our waveforms to the post-Newtonian (PN) BMS frame using a PN strain waveform. This helps us find the unique BMS transformation that minimizes the $L^2$ norm of the difference between the numerical and PN strain waveforms during the early inspiral phase. We find that once the waveforms are mapped to the PN BMS frame, they can be hybridized with a PN strain waveform much more effectively than if one used any of the previous alignment schemes, which only utilize the Poincaré transformations.

Numerical relativists can now produce gravitational waveforms with memory effects routinely and accurately. The gravitational-wave memory effect contains very low-frequency components, including a persistent offset. The presence of these components violates basic assumptions about time-shift behavior underpinning standard data-analysis techniques in gravitational-wave astronomy. This poses a challenge to the analysis of waveform spectra: how to preserve the low-frequency characteristics when transforming a time-domain waveform to the frequency domain. To tackle this challenge, we revisit the preprocessing procedures applied to the waveforms that contain memory effects. We find inconsistency between the zero-frequency limit of displacement memory and the low-frequency spectrum of the same memory preprocessed using the common scheme in literature. To resolve the inconsistency, we propose a new robust preprocessing scheme that produces the spectra of memory waveforms more faithfully. Using this new scheme, we inspect several characteristics of the spectrum of a memory waveform. In particular, we find a discernible beating pattern formed by the dominant oscillatory mode and the displacement memory. This pattern is absent in the spectrum of a waveform without memory. The difference between the memory and no-memory waveforms is too small to be observed by current-generation detectors in a single binary-black-hole event. Detecting the memory in a single event is likely to occur in the era of next-generation detectors.

The ringdown portion of a binary black hole merger consists of a sum of modes, each containing an infinite number of tones that are exponentially damped sinusoids. In principle, these can be measured as gravitational-waves with observatories like LIGO/Virgo/KAGRA, however in practice it is unclear how many tones can be meaningfully resolved. We investigate the consistency and resolvability of the overtones of the quadrupolar $\ell =m=2$ mode by starting at late times when the gravitational waveform is expected to be well approximated by the $\ell mn=220$ tone alone. We present a Bayesian inference framework to measure the tones in numerical relativity data. We measure tones at different start times, checking for consistency: we classify a tone as stably recovered if and only if the 95% credible intervals for amplitude and phase at time $t$ overlap with the credible intervals at all subsequent times. We test a set of tones including the first four overtones of the fundamental mode and the 320 tone and find that the 220 and 221 tones can be measured consistently with the inclusion of additional overtones. The 222 tone measurements can be stabilized when we include the 223 tone, but only in a narrow time window, after which it is too weak to measure. The 223 tone recovery appears to be unstable, and does not become stable with the introduction of the 224 tone. We find that $N=3$ tones can be stably recovered simultaneously. However, when analyzing $N\ge 4$ tones, the amplitude of one tone is consistent with zero. Thus, within our framework, one can identify only $N=3$ tones with nonzero amplitude that are simultaneously stable.

Using gravitational waves to probe the geometry of the ringing remnant black hole formed in a binary black hole coalescence is a well-established way to test Einstein’s theory of general relativity. However, doing so requires knowledge of when the predictions of black hole perturbation theory, i.e., quasinormal modes (QNMs), are a valid description of the emitted gravitational wave as well as what the amplitudes of these excitations are. In this work, we develop an algorithm to systematically extract QNMs from the ringdown of black hole merger simulations. Our algorithm improves upon previous ones in three ways: it fits over the two-sphere, enabling a complete model of the strain; it performs a reverse search in time for QNMs using a more robust nonlinear least squares routine called varpro; and it checks the variance of QNM amplitudes, which we refer to as “stability,” over an interval matching the natural timescale of each QNM. Using this algorithm, we not only demonstrate the stability of a multitude of QNMs and their overtones across the parameter space of quasicircular, nonprecessing binary black holes, but we also identify new quadratic QNMs that may be detectable in the near future using ground-based interferometers. Furthermore, we provide evidence which suggests that the source of remnant black hole perturbations is roughly independent of the overtone index in a given angular harmonic across binary parameter space, at least for overtones with $n\lesssim 2$. This finding may hint at the spatiotemporal structure of ringdown perturbations in black hole coalescences, as well as the regime of validity of perturbation theory in the ringdown of these events.

We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical experiments of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell =|m|=2$ quadrupolar modes, we find good agreement (within $\sim 10\%$) between these approaches, with discrepancies attributed to truncation error and uncertainties from mode fitting. Our results align with earlier studies extracting the coupling coefficients from select binary black hole merger simulations, showing consistency for the same remnant spins. Notably, the coupling coefficient is insensitive to a diverse range of initial data, including configurations that led to a significant (up to 5%) increase in the remnant black hole mass. These findings present opportunities for testing the nonlinear dynamics of general relativity with ground-based gravitational wave observatories. Lastly, we provide evidence of a bifurcation in coupling coefficients between counterrotating and corotating quasinormal modes as black hole spin increases.

The Bondi–van der Burg–Metzner–Sachs (BMS) frame of gravitational waves produced by numerical relativity simulations is crucial for building accurate waveform models. A proper comparison of numerical relativity waveforms with other models requires fixing the arbitrary BMS frame. In this work we improve the center-of-mass (c.m.) frame fixing for quasicircular, nonprecessing binary systems. Past work approximated the c.m. motion with just a linear fit. We compute a post-Newtonian result of the boosted c.m. charge to also capture its physical outspiraling oscillations. We show that using the analytical results improves the robustness of the fit parameters—translation and boost vectors—to the choice of duration and time of the fitting window. Our analysis demonstrates a maximum improvement in robustness when the window is placed at the center of the inspiral. We quantified this improvement by computing the ratio of variances of fit parameters when the fit window size is varied. The largest improvement in robustness of parameters is by a factor of $\sim 25$ for the boost vector and $\sim 20$ for the translation vector. Finally, we incorporate this method into the BMS frame-fixing routine of the python package scri for waveforms produced with Cauchy-characteristic evolution.

Highly precise and robust waveform models are required as improvements in detector sensitivity enable us to test general relativity with more precision than ever before. In this work, we introduce a spin-aligned surrogate ringdown model. This ringdown surrogate, nrsur3dq8_rd, is built with numerical waveforms produced using Cauchy-characteristic evolution. In addition, these waveforms are in the super-rest frame of the remnant black hole allowing us to do a correct analysis of the ringdown spectrum. The novel prediction of our surrogate model is complex-valued quasinormal mode (QNM) amplitudes, with median relative errors of ${10}^{-2}-{10}^{-3}$ over the parameter space. Like previous remnant surrogates, we also predict the remnant black hole’s mass and spin. The QNM mode amplitude errors translate into median errors on ringdown waveforms of $\sim {10}^{-4}$. The high accuracy and QNM mode content provided by our surrogate will enable high-precision ringdown analyses such as tests of general relativity. Our ringdown model is publicly available through the python package surfinbh.

Using high-accuracy numerical relativity waveforms, we confirm the presence of numerous overtones of the $\ell =2$, $m=2$ quasinormal mode early in the ringdown of binary black hole mergers. We do this by demonstrating the stability of the mode amplitudes at different fit times, ruling out the possibility that a linear superposition of modes unphysically fits a highly nonlinear part of the waveform. We also find a number of previously unidentified subdominant second-order quasinormal modes in the (2, 2) mode. Even though these modes are mathematically nonlinear, they nevertheless confirm the validity of perturbation theory as a good approximation for describing much of the ringdown.

We present the first numerical relativity waveforms for binary black hole mergers produced using spectral methods that show both the displacement and the spin memory effects. Explicitly, we use the SXS (Simulating eXtreme Spacetimes) Collaboration’s spec code to run a Cauchy evolution of a binary black hole merger and then extract the gravitational wave strain using spectre’s version of a Cauchy-characteristic extraction. We find that we can accurately resolve the strain’s traditional $m=0$ memory modes and some of the $m\ne 0$ oscillatory memory modes that have previously only been theorized. We also perform a separate calculation of the memory using equations for the Bondi-Metzner-Sachs charges as well as the energy and angular momentum fluxes at asymptotic infinity. Our new calculation uses only the gravitational wave strain and two of the Weyl scalars at infinity. Also, this computation shows that the memory modes can be understood as a combination of a memory signal throughout the binary’s inspiral and merger phases, and a quasinormal mode signal near the ringdown phase. Additionally, we find that the magnetic memory, up to numerical error, is indeed zero as previously conjectured. Last, we find that signal-to-noise ratios of memory for LIGO, the Einstein Telescope, and the Laser Interferometer Space Antenna with these new waveforms and new memory calculation are larger than previous expectations based on post-Newtonian or minimal waveform models.

Numerical relativity (NR) simulations of binary black holes provide precise waveforms, but are typically too computationally expensive to produce waveforms with enough orbits to cover the whole frequency band of gravitational-wave observatories. Accordingly, it is important to be able to hybridize NR waveforms with analytic, post-Newtonian (PN) waveforms, which are accurate during the early inspiral phase. We show that to build such hybrids, it is crucial to both fix the Bondi-Metzner-Sachs (BMS) frame of the NR waveforms to match that of PN theory, and optimize over the PN parameters to mitigate the error caused by the discrepancy between NR and PN parameters. We test such a hybridization procedure including all spin-weighted spherical harmonic modes with $|m|\le \ell$ for $\ell \le 8$, using 29 NR waveforms with mass ratios $q\le 10$ and spin magnitudes $|{\chi }_1|,|{\chi }_2|\le 0.8$. We find that for spin-aligned systems, the PN and NR waveforms agree very well. The difference is limited by the small nonzero orbital eccentricity of the NR waveforms, or equivalently by the lack of eccentric terms in the PN waveforms. To maintain full accuracy of the simulations, the matching window for spin-aligned systems should be at least 5 orbits long and end at least 15 orbits before merger. For precessing systems, the errors are larger than for spin-aligned cases. The errors are likely limited by the absence of precession-related spin-spin PN terms. Using ${10}^{5}M$ long NR waveforms, we find that there is no optimal choice of the matching window within this time span, because the hybridization result for precessing cases is always better if using earlier or longer matching windows. We provide the mean orbital frequency of the smallest acceptable matching window as a function of the target error between the PN and NR waveforms and the black hole spins.

We present a new study of remnant black hole properties from 13 binary black hole systems, numerically evolved using the Spectral Einstein Code. The mass, spin, and recoil velocity of each remnant were determined quasilocally from apparent horizon data and asymptotically from Bondi data $(h,{\psi }_4,{\psi }_3,{\psi }_2,{\psi }_1)$ computed at future null infinity using SpECTRE’s Cauchy characteristic evolution. We compare these independent measurements of the remnant properties in the bulk and on the boundary of the spacetime, giving insight into how well asymptotic data are able to reproduce local properties of the remnant black hole in numerical relativity. We also discuss the theoretical framework for connecting horizon quantities to asymptotic quantities and how it relates to our results. This study recommends a simple improvement to the recoil velocities reported in the Simulating eXtreme Spacetimes waveform catalog, provides an improvement to future surrogate remnant models, and offers new analysis techniques for evaluating the physical accuracy of numerical simulations.

The spectroscopic study of black hole quasinormal modes in gravitational-wave ringdown observations is hindered by our ignorance of which modes should dominate astrophysical signals for different binary configurations, limiting tests of general relativity and astrophysics. In this work, we present a description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasicircular precessing binary black hole coalescences—a key region of parameter space that has yet to be fully explored within the framework of black hole spectroscopy. We suggest that the remnant perturbation for precessing and nonprecessing systems is approximately the same up to a rotation, which implies that the relative amplitudes of the quasinormal modes in both systems are also related by a rotation. We present evidence for this by analyzing an extensive catalog of numerical relativity simulations. Additional structure in the amplitudes is connected to the system’s kick velocity and other asymmetries in the orbital dynamics. We find that the ringdowns of precessing systems need not be dominated by the $(\ell ,m,n)=(2,\pm 2,0)$ quasinormal modes, and that instead the $(2,\pm 1,0)$ or (2, 0,0) quasinormal modes can dominate. Our results are consistent with a ringdown analysis of the LIGO-Virgo gravitational wave signal GW190521, and may also help in understanding phenomenological inspiral-merger-ringdown waveform model systematics.

Understanding the characteristics of the remnant black hole formed in a binary black hole merger is crucial for conducting gravitational wave astronomy. Typically, models of remnant black holes provide information about their mass, spin, and kick velocity. However, other information related to the supertranslation symmetries of the BMS group, such as the memory effect, is also important for characterizing the final state of the system. In this work, we build a model of the BMS transformation that maps a binary black hole’s inspiral frame to the remnant black hole’s canonical rest frame. Training data for this model are created using high-precision numerical relativity simulations of quasicircular systems with mass ratios $q\le 8$ and spins parallel to the orbital angular momentum with magnitudes ${\chi }_1$, ${\chi }_2\le 0.8$. We use Gaussian process regression to model the BMS transformations over the three-dimensional parameter space $(q,{\chi }_1^z,{\chi }_2^z)$. The physics captured by this model is strictly nonperturbative and cannot be obtained from post-Newtonian approximations alone, as it requires knowledge of the strong nonlinear effects that are sourced during the merger. Apart from providing the first model of the supertranslation induced by a binary black hole merger, we also find that the kick velocities predicted using Cauchy-characteristic evolution waveforms are, on average, $\sim 5\%$ larger than the ones obtained from extrapolated waveforms. Our work has broad implications for improving gravitational wave models and studying the large-scale impact of memory, such as on the cosmological background. The fits produced in this work are available through the python package surfinbh under the name nrsur3dq8bmsremnant.

We present SEOBNRv5HM, a more accurate and faster inspiral-merger-ringdown gravitational waveform model for quasicircular, spinning, nonprecessing binary black holes within the effective-one-body (EOB) formalism. Compared to its predecessor, SEOBNRv4HM, the waveform model (i) incorporates recent high-order post-Newtonian results in the inspiral, with improved resummations, (ii) includes the gravitational modes $(\ell ,|m|)=(3,2),(4,3)$, in addition to the (2,2), (3,3), (2,1), (4,4), (5,5) modes already implemented in SEOBNRv4HM, (iii) is calibrated to larger mass ratios and spins using a catalog of 442 numerical-relativity (NR) simulations and 13 additional waveforms from black-hole perturbation theory, and (iv) incorporates information from second-order gravitational self-force in the nonspinning modes and radiation-reaction force. Computing the unfaithfulness against NR simulations, we find that for the dominant (2,2) mode the maximum unfaithfulness in the total mass range $10–300M_{\odot }$ is below ${10}^{-3}$ for 90% of the cases (38% for SEOBNRv4HM). When including all modes up to $\ell =5$ we find 98% (49%) of the cases with unfaithfulness below ${10}^{-2}$ (${10}^{-3}$), while these numbers reduce to 88% (5%) when using SEOBNRv4HM. Furthermore, the model shows improved agreement with NR in other dynamical quantities (e.g., the angular momentum flux and binding energy), providing a powerful check of its physical robustness. We implemented the waveform model in a high-performance python package (pyseobnr), which leads to evaluation times faster than SEOBNRv4HM by a factor of 10 to 50, depending on the configuration, and provides the flexibility to easily include spin-precession and eccentric effects, thus making it the starting point for a new generation of EOBNR waveform models (SEOBNRv5) to be employed for upcoming observing runs of the LIGO-Virgo-KAGRA detectors.

Numerical relativity simulations provide a full description of the dynamics of binary systems, including gravitational radiation. The waveforms produced by these simulations have a number of applications in gravitational-wave detection and inference. In this work, we revisit the accuracy of the waveforms produced by the spectral Einstein code. Motivated by the wide range of waveform applications, we propose and explore three accuracy metrics between simulation resolutions: (i) the generalized frequency-weighted mismatch, (ii) the relative amplitude difference, and (iii) the phase difference at different times. We confirm that numerical errors accumulate over the binary evolution, but the error is not intrinsically larger during the latest, more dynamical stages. Studying errors across the parameter space, we identify a positive correlation between both the mismatch and the phase difference with precessing spin, but little correlation with aligned spin or eccentricity. Lastly, amplitude and phases differences are symmetric upon exchanging resolutions across the catalog, suggesting that the dominant source of error is random, rather than something systematic that affects all waveforms similarly.

With the multitude of gravitational wave observations that have been made in the past ten years, probing the dynamical and nonlinear nature of strong gravity is becoming more and more feasible. One promising way to test the nonlinear nature of Einstein’s theory of general relativity (GR) is through the gravitational wave null memory effect: a nonlinear prediction of GR which corresponds to initially comoving observers being permanently displaced due to a burst of gravitational radiation. Previous studies have shown that, while it is unlikely that the memory effect will be observed in a single event by the LIGO-Virgo-KAGRA (LVK) detectors, evidence for memory in the population of LVK events should be attainable after $\sim 2000$ gravitational wave detections. Many of these works, however, largely relied on Bayes factors to perform their memory analyses: an approach that can depend sensitively on the analysis priors and, when multiplied across many events, can even favor incorrect conclusions. In this work, using the GWTC-5.0 catalog of binary black hole observations, we instead perform hierarchical Bayesian inference—which is not subject to the issues associated with Bayes factors—to measure the evidence for memory in current LVK observations. We find that we can constrain what we call the memory enhancement factor—the constant appearing in front of the contribution to the strain from the supermomentum flux—to ${0.26}_{-4.08}^{+4.09}$ (with $\pm$ values denoting the 68% credible interval), consistent with its GR value of 1. We also forecast that $\sim 2000$ detections will be needed to constrain the memory enhancement factor away from zero at the $1\sigma$ level.