skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on December 1, 2025

Title: Phenomenological gravitational waveform model of binary black holes incorporating horizon fluxes
Subjected to the tidal field of its companion, each component of a coalescing binary suffers a slow change in its mass (tidal heating) and spin (tidal torquing) during the inspiral and merger. For black holes, these changes are associated with their absorption of energy and angular momentum fluxes. This effect modifies the inspiral rate of the binary, and consequently, the phase and amplitude of its gravitational waveform. Numerical relativity (NR) waveforms contain these effects inherently, whereas analytical approximants for the early inspiral phase have to include them manually in the energy balance equation. In this work, we construct IMRPhenomD_Horizon, a frequency-domain gravitational waveform model that incorporates the effects of tidal heating of black holes. This is achieved by recalibrating the inspiral phase of the waveform model IMRPhenomD to incorporate the phase corrections for tidal heating. We also include corrections to the amplitude, but add them directly to the inspiral amplitude model of IMRPhenomD. First we demonstrate that the inclusion of the corrections, especially in the phase, confers an overall improvement in the phase agreement between the analytical inspiral model (uncalibrated SEOBNRv2) and NR data. The model presented here is faithful, with less than 1% mismatches against a set of hybrid waveforms (except for one outlier that barely breaches this limit). The recalibrated model shows mismatches of up to ∼14% with IMRPhenomD for high mass ratios and spins. Amplitude corrections become less significant for higher mass ratios, whereas the phase corrections leave more impact—suggesting that the former is practically irrelevant for gravitational wave data analysis in Advanced LIGO (aLIGO), Virgo and KAGRA. Comparing with a set of 219 numerical relativity waveforms, we find that the median of mismatches decreases by ∼4% in aLIGO zero-detuned high power noise curve, and by ∼1.5% with a flat noise curve. This implies a modest but notable improvement in waveform accuracy.  more » « less
Award ID(s):
2309352
PAR ID:
10578933
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical Review D
Volume:
110
Issue:
12
ISSN:
2470-0010
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Numerical relativity (NR) simulations of binary black hole (BBH) systems provide the most accurate gravitational wave predictions, but at a high computational cost—especially when the black holes have nearly extremal spins (i.e. spins near the theoretical upper limit) or very unequal masses. Recently, the technique of reduced order modeling has enabled the construction of ‘surrogate models’ trained on an existing set of NR waveforms. Surrogate models enable the rapid computation of the gravitational waves emitted by BBHs. Typically these models are used for interpolation to compute gravitational waveforms for BBHs with mass ratios and spins within the bounds of the training set. Because simulations with nearly extremal spins are so technically challenging, surrogate models almost always rely on training sets with only moderate spins. In this paper, we explore how well surrogate models can extrapolate to nearly extremal spins when the training set only includes moderate spins. For simplicity, we focus on one-dimensional surrogate models trained on NR simulations of BBHs with equal masses and equal, aligned spins. We assess the performance of the surrogate models at higher spin magnitudes by calculating the mismatches between extrapolated surrogate model waveforms and NR waveforms, by calculating the differences between extrapolated and NR measurements of the remnant black-hole mass, and by testing how the surrogate model improves as the training set extends to higher spins. We find that while extrapolation in this one-dimensional case is viable for current detector sensitivities, surrogate models for next-generation detectors should use training sets that extend to nearly extremal spins. 
    more » « less
  2. Abstract Gravitational waves emitted by black hole binary inspiral and mergers enable unprecedented strong-field tests of gravity, requiring accurate theoretical modeling of the expected signals in extensions of general relativity. In this paper we model the gravitational wave emission of inspiralling binaries in scalar Gauss–Bonnet gravity theories. Going beyond the weak-coupling approximation, we derive the gravitational waveform to relative first post-Newtonian order beyond the quadrupole approximation and calculate new contributions from nonlinear curvature terms. We also compute the scalar waveform to relative 0.5PN order beyond the leading −0.5PN order terms. We quantify the effect of these terms and provide ready-to-implement gravitational wave and scalar waveforms as well as the Fourier domain phase for quasi-circular binaries. We also perform a parameter space study, which indicates that the values of black hole scalar charges play a crucial role in the detectability of deviation from general relativity. We also compare the scalar waveforms to numerical relativity simulations to assess the impact of the relativistic corrections to the scalar radiation. Our results provide important foundations for future precision tests of gravity. 
    more » « less
  3. 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
  4. The inspiral-merger-ringdown (IMR) consistency test checks the consistency of the final mass and final spin of a binary black hole merger remnant, independently inferred via the inspiral and merger-ringdown parts of the waveform. As binaries are expected to be nearly circularized when entering the frequency band of ground-based detectors, tests of general relativity (GR) currently employ quasicircular waveforms. We quantify the effect of residual orbital eccentricity on the IMR consistency test. We find that eccentricity causes a significant systematic bias in the inferred final mass and spin of the remnant black hole at an orbital eccentricity (defined at 10 Hz) of e0≳0.1 in the LIGO band (for a total binary mass in the range 65-200M⊙). For binary black holes observed by Cosmic Explorer (CE), the systematic bias becomes significant for e0≳0.015 (for 200-600M⊙ systems). This eccentricity-induced bias on the final mass and spin leads to an apparent inconsistency in the IMR consistency test, manifesting as a false violation of GR. Hence, eccentric corrections to waveform models are important for constructing a robust test of GR, especially for third-generation detectors. We also estimate the eccentric corrections to the relationship between the inspiral parameters and the final mass and final spin; they are shown to be quite small. 
    more » « less
  5. Binary black holes with misaligned spins will generically induce both precession and nutation of the orbital angular momentum 𝐋 about the total angular momentum 𝐉. These phenomena modulate the phase and amplitude of the gravitational waves emitted as the binary inspirals to merger. We introduce a “taxonomy” of binary black-hole spin precession that encompasses all the known phenomenology, then present five new phenomenological parameters that describe generic precession and constitute potential building blocks for future gravitational waveform models. These are the precession amplitude ⟨𝜃𝐿⟩, the precession frequency ⟨Ω𝐿⟩, the nutation amplitude Δ⁢𝜃𝐿, the nutation frequency 𝜔, and the precession-frequency variation Δ⁢Ω𝐿. We investigate the evolution of these five parameters during the inspiral and explore their statistical properties for sources with isotropic spins. In particular, we find that nutation of 𝐋 is most prominent for binaries with high spins (𝜒≳0.5) and moderate mass ratios (𝑞∼0.6). 
    more » « less