Title: Automated determination of the end time of junk radiation in binary black hole simulations
Abstract When numerically solving Einstein’s equations for the evolution of binary black holes, physical imperfections in the initial data manifest as a transient, high-frequency pulse of ‘junk radiation.’ This unphysical signal must be removed before the waveform can be used. Improvements in the efficiency of numerical simulations now allow waveform catalogs containing thousands of waveforms to be produced. Thus, an automated procedure for identifying junk radiation is required. To this end, we present a new algorithm based on the empirical mode decomposition (EMD) from the Hilbert–Huang transform. This approach allows us to isolate and measure the high-frequency oscillations present in the measured irreducible masses of the black holes. The decay of these oscillations allows us to estimate the time from which the junk radiation can be ignored. To make this procedure more precise, we propose three distinct threshold criteria that specify how small the contribution of junk radiation has to be before it can be considered negligible. We apply this algorithm to 3403 BBH simulations from the Simulating eXtreme Spacetime catalog to find appropriate values for the thresholds in the three criteria. We find that this approach yields reliable decay time estimates, i.e. when to consider the simulation physical, for > 98.5% of the simulations studied. This demonstrates the efficacy of the EMD as a suitable tool to automatically isolate and characterize junk radiation in the simulation of binary black hole systems.  more » « less
Award ID(s):
2308615 2209655 2309211 2209656
PAR ID:
10589158
Author(s) / Creator(s):
; ;
Publisher / Repository:
IOP
Date Published:
Journal Name:
Classical and Quantum Gravity
Volume:
42
Issue:
7
ISSN:
0264-9381
Page Range / eLocation ID:
075004
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The Simulating eXtreme Spacetimes Collaboration’s codeSpECcan now routinely simulate binary black hole mergers undergoing 25 orbits, with the longest simulations undergoing nearly 180 orbits. While this sounds impressive, the mismatch between the highest resolutions for this long simulation is O ( 10 1 ) . Meanwhile, the mismatch between resolutions for the more typical simulations tends to be 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 O ( 10 3 ) 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
  2. Abstract We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from 47 numerical relativity simulations. We show that the ( = 2 , m = 2 ) and ( = 2 , m = 1 ) multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratioqand reduced tidal parameter Λ ~ . Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories. 
    more » « less
  3. Abstract Gravitational wave searches are crucial for studying compact sources such as neutron stars and black holes. Many sensitive modeled searches use matched filtering to compare gravitational strain data to a set of waveform models known as template banks. We introduce a new stochastic placement method for constructing template banks, offering efficiency and flexibility to handle arbitrary parameter spaces, including orbital eccentricity, tidal deformability, and other extrinsic parameters. This method can be computationally limited by the ability to compare proposal templates with the accepted templates in the bank. To alleviate this computational load, we introduce the use of inner product inequalities to reduce the number of required comparisons. We also introduce a novel application of Gaussian Kernel Density Estimation to enhance waveform coverage in sparser regions. Our approach has been employed to search for eccentric binary neutron stars, low-mass neutron stars, primordial black holes, and supermassive black hole binaries. We demonstrate that our method produces self-consistent banks that recover the required minimum fraction of signals. For common parameter spaces, our method shows comparable computational performance and similar template bank sizes to geometric placement methods and stochastic methods, while easily extending to higher-dimensional problems. The time to run a search exceeds the time to generate the bank by a factor of O ( 10 5 ) for dedicated template banks, such as geometric, mass-only stochastic, and aligned spin cases, O ( 10 4 ) for eccentric and O ( 10 3 ) for the tidal deformable bank. With the advent of efficient template bank generation, the primary area for improvement is developing more efficient search methodologies. 
    more » « less
  4. Abstract We present a major update to the Simulating eXtreme Spacetimes (SXSs) Collaboration’s catalog of binary black hole (BBH) simulations. Using highly efficient spectral methods implemented in the Spectral Einstein Code (SpEC), we have nearly doubled the total number of binary configurations from 2018 to 3756. The catalog now more densely covers the parameter space with precessing simulations up to mass ratioq = 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 ofSpECor 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 BBH simulations is over 1000 times more efficient than previously published finite-difference simulations. The full catalog is publicly available through thesxsPython package and athttps://data.black-holes.org . 
    more » « less
  5. Abstract On 2023 November 23, the two LIGO observatories both detected GW231123, a gravitational-wave signal consistent with the merger of two black holes with masses 13 7 18 + 23 M and 10 1 50 + 22 M (90% credible intervals), at a luminosity distance of 0.7–4.1 Gpc, a redshift of 0.4 0 0.25 + 0.27 , and with a network signal-to-noise ratio of ∼20.7. Both black holes exhibit high spins— 0.9 0 0.19 + 0.10 and 0.8 0 0.52 + 0.20 , respectively. A massive black hole remnant is supported by an independent ringdown analysis. Some properties of GW231123 are subject to large systematic uncertainties, as indicated by differences in the inferred parameters between signal models. The primary black hole lies within or above the theorized mass gap where black holes between 60–130Mshould be rare, due to pair-instability mechanisms, while the secondary spans the gap. The observation of GW231123 therefore suggests the formation of black holes from channels beyond standard stellar collapse and that intermediate-mass black holes of mass ∼200Mform through gravitational-wave-driven mergers. 
    more » « less