An official website of the United States government
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
Comparing Numerical Relativity and Perturbation Theory Waveforms for a Non-Spinning Equal-Mass Binary
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
- PAR ID:
- 10495477
- Publisher / Repository:
- MDPI
- Date Published:
- Journal Name:
- Universe
- Volume:
- 10
- Issue:
- 1
- ISSN:
- 2218-1997
- Page Range / eLocation ID:
- 25
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Simulations to calculate a single gravitational waveform (GW) can take several weeks. Yet, thousands of such simulations are needed for the detection and interpretation of gravitational waves. Future detectors will require even more accurate waveforms than those currently used. We present here the first large scale, adaptive mesh, multi-GPU numerical relativity (NR) code together with performance analysis and benchmarking. While comparisons are difficult to make, our GPU extension of the Dendro-GR NR code achieves a 6x speedup over existing state-of-the-art codes. We achieve 800 GFlops/s on a single NVIDIA A100 GPU with an overall 2.5x speedup over a two-socket, 128-core AMD EPYC 7763 CPU node with an equivalent CPU implementation. We present detailed performance analyses, parallel scalability results, and accuracy assessments for GWs computed for mass ratios q=1,2,4. We also present strong scalability up to 8 A100s and weak scaling up to 229,376 ×86 cores on the Texas Advanced Computing Center's Frontera system.more » « less
-
For a finite point set P⊂R^d, denote by diam(P) the ratio of the largest to the smallest distances between pairs of points in P. Let c_{d,α}(n) be the largest integer c such that any n-point set P⊂R^d in general position, satisfying diam(P)<αn^{1/d}, contains an c-point convex independent subset. We determine the asymptotics of c_{d,α}(n) as n→∞ by showing the existence of positive constants β=β(d,α) and γ=γ(d) such that βn^{(d−1)/(d+1)}≤c_{d,α}(n)≤γn^{(d−1)/(d+1)} for α≥2.more » « less
-
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
-
Abstract P-type point contact (PPC) HPGe detectors are a leading technology for rare event searches due to their excellent energy resolution, low thresholds, and multi-site event rejection capabilities. We have characterized a PPC detector’s response to $$\alpha $$ α particles incident on the sensitive passivated and p $$^+$$ + surfaces, a previously poorly-understood source of background. The detector studied is identical to those in the Majorana Demonstrator experiment, a search for neutrinoless double-beta decay ( $$0\nu \beta \beta $$ 0 ν β β ) in $$^{76}$$ 76 Ge. $$\alpha $$ α decays on most of the passivated surface exhibit significant energy loss due to charge trapping, with waveforms exhibiting a delayed charge recovery (DCR) signature caused by the slow collection of a fraction of the trapped charge. The DCR is found to be complementary to existing methods of $$\alpha $$ α identification, reliably identifying $$\alpha $$ α background events on the passivated surface of the detector. We demonstrate effective rejection of all surface $$\alpha $$ α events (to within statistical uncertainty) with a loss of only 0.2% of bulk events by combining the DCR discriminator with previously-used methods. The DCR discriminator has been used to reduce the background rate in the $$0\nu \beta \beta $$ 0 ν β β region of interest window by an order of magnitude in the Majorana Demonstrator and will be used in the upcoming LEGEND-200 experiment.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/universe10010025
