An official website of the United States government
Detections of gravitational waves emitted from binary black hole coalescences allow us to probe the strong-field dynamics of general relativity (GR). One can compare the observed gravitational-wave signals with theoretical waveform models to constrain possible deviations from GR. Any physics that is not included in these waveform models might show up as apparent GR deviations. The waveform models used in current tests of GR describe binaries on quasicircular orbits, since most of the binaries detected by ground-based gravitational-wave detectors are expected to have negligible eccentricities. Thus, a signal from an eccentric binary in GR is likely to show up as a deviation from GR in the current implementation of these tests. We study the response of four standard tests of GR to eccentric binary black hole signals with the forecast O4 sensitivity of the LIGO-Virgo network. Specifically, we consider two parametrized tests (TIGER and FTI), the modified dispersion relation test, and the inspiral-merger-ringdown consistency test. To model eccentric signals, we use nonspinning numerical relativity simulations from the SXS catalog with three mass ratios (1, 2, 3), which we scale to a redshifted total mass of 80M⊙ and luminosity distance of 400 Mpc. For each of these mass ratios, we consider signals with eccentricities of ∼0.05 and ∼0.1 at 17 Hz. We find that signals with larger eccentricity lead to very significant false GR deviations in most tests while signals having smaller eccentricity lead to significant deviations in some tests. For the larger eccentricity cases, one would even get a deviation from GR with TIGER at ∼90% credibility at a distance of ≳1.5 Gpc. Thus, it will be necessary to exclude the possibility of an eccentric binary in order to make any claim about detecting a deviation from GR.
more »
« less
The SXS Collaboration’s third catalog of binary black hole simulations
Abstract We present a major update to the Simulating eXtreme Spacetimes (SXS) Collaboration’s catalog of binary black hole simulations. Using highly efficient spectral methods implemented in the Spectral Einstein Code (SpEC), we have nearly doubled the total number of binary configurations from 2,018 to 3,756. The catalog now more densely covers the parameter space with precessing simulations up to mass ratio q = 8 and dimensionless spins up to |χ⃗| ≤ 0.8 with near-zero eccentricity. The catalog also includes some simulations at higher mass ratios with moderate spin and more than 250 eccentric simulations. We have also deprecated and rerun some simulations from our previous catalog (e.g., simulations run with a much older version of SpEC or that had anomalously high errors in the waveform). The median waveform difference (which is similar to the mismatch) between resolutions over the simulations in the catalog is 4 × 10−4. The simulations have a median of 22 orbits, while the longest simulation has 148 orbits. We have corrected each waveform in the catalog to be in the binary’s center-of-mass frame and exhibit gravitational-wave memory. We estimate the total CPU cost of all simulations in the catalog to be 480,000,000 core-hours. We find that using spectral methods for binary black hole simulations is over 1,000 times more efficient than previously published finite-difference simulations. The full catalog is publicly available through the sxs Python package and at https://data.black-holes.org .
more »
« less
- PAR ID:
- 10629213
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Classical and Quantum Gravity
- ISSN:
- 0264-9381
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Simulating eXtreme Spacetimes Collaboration's code \texttt{SpEC} can now routinely simulate binary black hole mergers undergoing $$\sim25$$ orbits, with the longest simulations undergoing nearly $$\sim180$$ orbits. While this sounds impressive, the mismatch between the highest resolutions for this long simulation is $$\mathcal{O}(10^{-1})$$. Meanwhile, the mismatch between resolutions for the more typical simulations tends to be $$\mathcal{O}(10^{-4})$$, despite the resolutions being similar to the long simulations'. In this note, we explain why mismatch alone gives an incomplete picture of code---and waveform---quality, especially in the context of providing waveform templates for LISA and 3G detectors, which require templates with $$\mathcal{O}(10^{3}) - \mathcal{O}(10^{5})$$ orbits. We argue that to ready the GW community for the sensitivity of future detectors, numerical relativity groups must be aware of this caveat, and also run future simulations with at least three resolutions to properly assess waveform accuracy.more » « less
-
Higher-order gravitational wave modes from quasi-circular, spinning, non-precessing binary black hole mergers encode key information about these systems' nonlinear dynamics. We model these waveforms using transformer architectures, targeting the evolution from late inspiral through ringdown. Our data is derived from the \texttt{NRHybSur3dq8} surrogate model, which includes spherical harmonic modes up to ℓ≤4 (excluding (4,0), (4,±1) and including (5,5) modes). These waveforms span mass ratios q≤8, spin components sz1,2∈[−0.8,0.8], and inclination angles θ∈[0,π]. The model processes input data over the time interval t∈[−5000M,−100M) and generates predictions for the plus and cross polarizations, (h+,h×), over the interval t∈[−100M,130M]. Utilizing 16 NVIDIA A100 GPUs on the Delta supercomputer, we trained the transformer model in 15 hours on over 14 million samples. The model's performance was evaluated on a test dataset of 840,000 samples, achieving mean and median overlap scores of 0.996 and 0.997, respectively, relative to the surrogate-based ground truth signals. We further benchmark the model on numerical relativity waveforms from the SXS catalog, finding that it generalizes well to out-of-distribution systems, capable of reproducing the dynamics of systems with mass ratios up to q=15 and spin magnitudes up to 0.998, with a median overlap of 0.969 across 521 NR waveforms and up to 0.998 in face-on/off configurations. These results demonstrate that transformer-based models can capture the nonlinear dynamics of binary black hole mergers with high accuracy, even outside the surrogate training domain, enabling fast sequence modeling of higher-order wave modes.more » « less
-
Abstract Binary black holes are the most abundant source of gravitational-wave observations. Gravitational-wave observatories in the next decade will require tremendous increases in the accuracy of numerical waveforms modeling binary black holes, compared to today’s state of the art. One approach to achieving the required accuracy is using spectral-type methods that scale to many processors. Using theSpECTREnumerical-relativity (NR) code, we present the first simulations of a binary black hole inspiral, merger, and ringdown using discontinuous Galerkin (DG) methods. The efficiency of DG methods allows us to evolve the binary through ∼ 18 orbits at reasonable computational cost. We then useSpECTRE’s Cauchy Characteristic Evolution (CCE) code to extract the gravitational waves at future null infinity. The open-source nature ofSpECTREmeans this is the first time a spectral-type method for simulating binary black hole evolutions is available to the entire NR community.more » « less
-
null (Ed.)Accurate models of gravitational waves from merging binary black holes are crucial for detectors to measure events and extract new science. One important feature that is currently missing from the Simulating eXtreme Spacetimes (SXS) Collaboration’s catalog of waveforms for merging black holes, and other waveform catalogs, is the gravitational memory effect: a persistent, physical change to spacetime that is induced by the passage of transient radiation. We find, however, that by exploiting the Bondi-van der Burg-Metzner-Sachs (BMS) balance laws, which come from the extended BMS transformations, we can correct the strain waveforms in the SXS catalog to include the missing displacement memory. Our results show that these corrected waveforms satisfy the BMS balance laws to a much higher degree of accuracy. Furthermore, we find that these corrected strain waveforms coincide especially well with the waveforms obtained from Cauchy-characteristic extraction (CCE) that already exhibit memory effects. These corrected strain waveforms also evade the transient junk effects that are currently present in CCE waveforms. Last, we make our code for computing these contributions to the BMS balance laws and memory publicly available as a part of the python package sxs, thus enabling anyone to evaluate the expected memory effects and violation of the BMS balance laws.more » « less
