An official website of the United States government
Abstract Brownian coating thermal noise in detector test masses is limiting the sensitivity of current gravitational-wave detectors on Earth. Therefore, accurate numerical models can inform the ongoing effort to minimize Brownian coating thermal noise in current and future gravitational-wave detectors. Such numerical models typically require significant computational resources and time, and often involve closed-source commercial codes. In contrast, open-source codes give complete visibility and control of the simulated physics, enable direct assessment of the numerical accuracy, and support the reproducibility of results. In this article, we use the open-sourceSpECTREnumerical relativity code and adopt a novel discontinuous Galerkin numerical method to model Brownian coating thermal noise. We demonstrate thatSpECTREachieves significantly higher accuracy than a previous approach at a fraction of the computational cost. Furthermore, we numerically model Brownian coating thermal noise in multiple sub-wavelength crystalline coating layers for the first time. Our new numerical method has the potential to enable fast exploration of realistic mirror configurations, and hence to guide the search for optimal mirror geometries, beam shapes and coating materials for gravitational-wave detectors.
more »
« less
Simulating binary black hole mergers using discontinuous Galerkin methods
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
- PAR ID:
- 10565353
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Classical and Quantum Gravity
- Volume:
- 42
- Issue:
- 3
- ISSN:
- 0264-9381
- Format(s):
- Medium: X Size: Article No. 035001
- Size(s):
- Article No. 035001
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Errors due to imperfect boundary conditions in numerical relativity simulations of binary black holes (BBHs) can produce unphysical reflections of gravitational waves which compromise the accuracy of waveform predictions, especially for subdominant modes. A system of higher order absorbing boundary conditions which greatly reduces this problem was introduced in earlier work (Buchman and Sarbach 2006Class. Quantum Grav.236709). In this paper, we devise two new implementations of this boundary condition system in the Spectral Einstein Code (SpEC), and test them in both linear multipolar gravitational wave and inspiralling mass ratio 7:1 BBH simulations. One of our implementations in particular is shown to be extremely robust and to produce accuracy superior to the standard freezing-Ψ0boundary condition usually used bySpEC.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
-
ABSTRACT The detection of an intermediate-mass black hole population (102–106 M⊙) will provide clues to their formation environments (e.g. discs of active galactic nuclei, globular clusters) and illuminate a potential pathway to produce supermassive black holes. Ground-based gravitational-wave detectors are sensitive to mergers that can form intermediate-mass black holes weighing up to ∼450 M⊙. However, ground-based detector data contain numerous incoherent short duration noise transients that can mimic the gravitational-wave signals from merging intermediate-mass black holes, limiting the sensitivity of searches. Here, we follow-up on binary black hole merger candidates using a ranking statistic that measures the coherence or incoherence of triggers in multiple-detector data. We use this statistic to rank candidate events, initially identified by all-sky search pipelines, with lab-frame total masses ≳ 55 M⊙ using data from LIGO’s second observing run. Our analysis does not yield evidence for new intermediate-mass black holes. However, we find support for eight stellar-mass binary black holes not reported in the first LIGO–Virgo gravitational wave transient catalogue GWTC-1, seven of which have been previously reported by other catalogues.more » « less
-
Abstract The detection of orbital eccentricity for a binary black hole system via gravitational waves is a key signature to distinguish between the possible binary origins. The identification of eccentricity has been difficult so far due to the limited availability of eccentric gravitational waveforms over the full range of black hole masses and eccentricities. Here we evaluate the eccentricity of five black hole mergers detected by the LIGO and Virgo observatories using theTEOBResumS-DALI,TEOBResumS-GIOTTO, andTEOBResumSPmodels. This analysis studies eccentricities up to 0.6 at the reference frequency of 5 Hz and incorporates higher-order gravitational-wave modes critical to model emission from highly eccentric orbits. The binaries have been selected due to previous hints of eccentricity or due to their unusual mass and spin. While other studies found marginal evidence for eccentricity for some of these events, our analyses do not favor the incorporation of eccentricity compared to the quasi-circular case. While lacking the eccentric evidence of other analyses, we find our analyses marginally shifts the posterior in multiple parameters for several events when allowing eccentricity to be nonzero.more » « less
