An official website of the United States government
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
Length dependence of waveform mismatch: a caveat on waveform accuracy
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
- PAR ID:
- 10590488
- 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 Objective. Accurate prediction of motor unit (MU) discharge activity from surface electromyogram (sEMG) signals is critical for understanding neuromuscular control and for enabling practical neural interface applications. However, current MU decomposition approaches rely on person-specific data, limiting their generalizability. Approach. We developed a cross-person decomposition framework and validated the algorithm using synthesized high-density sEMG data by convoluting simulated MU firing spike trains with action potential templates derived from human experimental data. We first obtained separation matrix from multiple training subjects and applied them to decompose sEMG signals from unseen test subjects. This allowed us to obtain MU spike trains. The predicted outcomes were then compared with the ground truth across multiple metrics, including spike detection accuracy, MU firing rate (FR), waveform similarity of motor unit action potentials (MUAP), and MU recruitment thresholds. Main results. Our results demonstrated strong agreement between predicted and true MU activity. Specifically, we found high R² values (≥0.95) for the populational FR, and the coefficient of variation of FR remained stable across different MU retention thresholds. The MU similarity analyses revealed that the predicted MUAPs closely matched ground truth counterparts both in waveform shape and spatial distribution. Furthermore, recruitment thresholds exhibited strong linear relation (R² = 0.98 ± 0.006) with minimal error. Significance. These findings demonstrate the feasibility of efficient cross-person MU decomposition with minimal accuracy loss, laying the groundwork for generalized, plug-and-play myoelectric systems in neurophysiology, neuroprosthetic, and rehabilitation applications.more » « less
-
Abstract LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the early inspirals of stellar-mass black holes that will ultimately venture into the ground-based detectors’ view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA’s discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This White Paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.more » « less
-
Abstract We demonstrate how to quantify the frequency-domain amplitude and phase accuracy of waveform models, andδφ, in a form that could be marginalized over in gravitational-wave inference using techniques currently applied for quantifying calibration uncertainty. For concreteness, waveform uncertainties affecting neutron-star inspiral measurements are considered, and post-hoc error estimates from a variety of waveform models are made by comparing time-domain and frequency-domain analytic models with multiple-resolution numerical simulations. These waveform uncertainty estimates can be compared to GW170817 calibration envelopes or to Advanced LIGO and Virgo calibration goals. Signal-specific calibration and waveform uncertainties are compared to statistical fluctuations in gravitational-wave observatories, giving frequency-dependent modeling requirements for detectors such as Advanced LIGO Plus, Cosmic Explorer, or Einstein Telescope. Finally, the distribution of waveform error for the GW170817 posterior is computed from tidal models and compared to the constraints onδφor from GWTC-1 by Edelmanet al.In general,δφand can also be interpreted in terms of unmodeled astrophysical energy transfer within or from the source system.more » « less
-
Abstract Next-generation gravitational wave detectors such as Cosmic Explorer, the Einstein Telescope, and LISA, demand highly accurate and extensive gravitational wave (GW) catalogs to faithfully extract physical parameters from observed signals. However, numerical relativity (NR) faces significant challenges in generating these catalogs at the required scale and accuracy on modern computers, as NR codes do not fully exploit modern GPU capabilities. In response, we extend NRPy, a Python-based NR code-generation framework, to develop NRPyEllipticGPU—a CUDA-optimized elliptic solver tailored for the binary black hole (BBH) initial data problem. NRPyEllipticGPU is the first GPU-enabled elliptic solver in the NR community, supporting a variety of coordinate systems and demonstrating substantial performance improvements on both consumer-grade and HPC-grade GPUs. We show that, when compared to a high-end CPU, NRPyEllipticGPU achieves on a high- end GPU up to a sixteenfold speedup in single precision while increasing double- precision performance by a factor of 2–4. This performance boost leverages the GPU’s superior parallelism and memory bandwidth to achieve a compute-bound application and enhancing the overall simulation efficiency. As NRPyEllipticGPU shares the core infrastructure common to NR codes, this work serves as a practical guide for developing full, CUDA-optimized NR codes.more » « less
