An official website of the United States government
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
Numerical-relativity surrogate modeling with nearly extremal black-hole spins
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
- PAR ID:
- 10435491
- Date Published:
- Journal Name:
- Classical and Quantum Gravity
- Volume:
- 40
- Issue:
- 5
- ISSN:
- 0264-9381
- Page Range / eLocation ID:
- 055003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Past studies have empirically demonstrated a surprising agreement between gravitational waveforms computed using adiabatic–driven–inspiral point–particle black hole perturbation theory (ppBHPT) and numerical relativity (NR) following a straightforward calibration step, sometimes referred to as α-β scaling. Specifically focusing on the quadrupole mode, this calibration technique necessitates only two time-independent parameters to scale the overall amplitude and time coordinate. In this article, part of a Special Issue, we investigate this scaling for non-spinning binaries at the equal-mass limit. Even without calibration, NR and ppBHPT waveforms exhibit an unexpected degree of similarity after accounting for different mass scale definitions. Post-calibration, good agreement between ppBHPT and NR waveforms extends nearly up to the point of the merger. We also assess the breakdown of the time-independent assumption of the scaling parameters, shedding light on current limitations and suggesting potential generalizations for the α-β scaling technique.more » « less
-
Abstract Accurate modelling of black hole binaries is critical to achieve the science goals of gravitational-wave detectors. Modelling such configurations relies strongly on calibration to numerical-relativity (NR) simulations. Binaries on quasi-circular orbits have been widely explored in NR, however, coverage of the broader 9-dimensional parameter space, including orbital eccentricity, remains sparse. This article develops a new procedure to control orbital eccentricity of binary black hole simulations that enables choosing initial data parameters with precise control over eccentricity and mean anomaly of the subsequent evolution, as well as the coalescence time. We then calculate several sequences of NR simulations that nearly uniformly cover the 2-dimensional eccentricity--mean anomaly space for equal mass, non-spinning binary black holes. We demonstrate that, for fixed eccentricity, many quantities related to the merger dynamics of binary black holes show an oscillatory dependence on mean anomaly. The amplitude of these oscillations scales nearly linearly with the eccentricity of the system. We find that for the eccentricities explored in this work, the magnitude of deviations in various quantities such as the merger amplitude and peak luminosity can approach $$\sim5\%$$ of their quasi-circular value. We use our findings to explain eccentric phenomena reported in other studies. We also show that methods for estimating the remnant mass employed in the effective-one-body approach exhibit similar deviations, roughly matching the amplitude of the oscillations we find in NR simulations. This work is an important step towards a complete description of eccentric binary black hole mergers, and demonstrates the importance of considering the entire 2-dimensional parameter subspace related to eccentricity.more » « less
-
Abstract The component black holes (BHs) observed in gravitational-wave (GW) binary black hole (BBH) events tend to be more massive and slower spinning than those observed in black hole X-ray binaries (BH-XRBs). Without modeling their evolutionary histories, we investigate whether these apparent tensions in the BH populations can be explained by GW observational selection effects alone. We find that this is indeed the case for the discrepancy between BH masses in BBHs and the observed high-mass X-ray binaries (HMXBs), when we account for statistical uncertainty from the small sample size of just three HMXBs. On the other hand, the BHs in observed low-mass X-ray binaries (LMXBs) are significantly lighter than the astrophysical BBH population, but this may just be due to a correlation between component masses in a binary system. Given their light stellar companions, we expect light BHs in LMXBs. The observed spins in HMXBs and LMXBs, however, are in tension with the inferred BBH spin distribution at the >99.9% level. We discuss possible scenarios behind the significantly larger spins in observed BH-XRBs. One possibility is that a small subpopulation (conservatively <30%) of BBHs have rapidly spinning primary components, indicating that they may have followed a similar evolutionary pathway to the observed HMXBs. In LMXBs, it has been suggested that BHs can spin up by accretion. If LMXB natal spins follow the BBH spin distribution, we find LMXBs must gain an average dimensionless spin of 0.47 − 0.11 + 0.10 , but if their natal spins follow the observed HMXB spins, the average spin-up must be <0.03.more » « less
-
Abstract The population-level distributions of the masses, spins, and redshifts of binary black holes (BBHs) observed using gravitational waves can shed light on how these systems form and evolve. Because of the complex astrophysical processes shaping the inferred BBH population, models allowing for correlations among these parameters will be necessary to fully characterize these sources. We hierarchically analyze the BBH population detected by LIGO and Virgo with a model allowing for correlations between the effective aligned spin and the primary mass and redshift. We find that the width of the effective spin distribution grows with redshift at 98.6% credibility. We determine this trend to be robust under the application of several alternative models and additionally verify that such a correlation is unlikely to be spuriously introduced using a simulated population. We discuss the possibility that this correlation could be due to a change in the natal black hole spin distribution with redshift.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1361-6382/acb3a7
