skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Rapid likelihood free inference of compact binary coalescences using accelerated hardware
Abstract We report a gravitational-wave parameter estimation algorithm,AMPLFI, based on likelihood-free inference using normalizing flows. The focus ofAMPLFIis to perform real-time parameter estimation for candidates detected by machine-learning based compact binary coalescence search,Aframe. We present details of our algorithm and optimizations done related to data-loading and pre-processing on accelerated hardware. We train our model using binary black-hole (BBH) simulations on real LIGO-Virgo detector noise. Our model has 6 million trainable parameters with training times 24 h. Based on online deployment on a mock data stream of LIGO-Virgo data,Aframe+AMPLFIis able to pick up BBH candidates and infer parameters for real-time alerts from data acquisition with a net latency of 6 s.  more » « less
Award ID(s):
2308862 2117997
PAR ID:
10598595
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Purpose-Led Publishing
Date Published:
Journal Name:
Machine Learning: Science and Technology
Volume:
5
Issue:
4
ISSN:
2632-2153
Page Range / eLocation ID:
045030
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract The cluster mass–richness relation (MRR) is an observationally efficient and potentially powerful cosmological tool for constraining the matter density Ωmand the amplitude of fluctuationsσ8using the cluster abundance technique. We derive the MRR relation usingGalWCat19, a publicly available galaxy cluster catalog we created from the Sloan Digital Sky Survey-DR13 spectroscopic data set. In the MRR, cluster mass scales with richness as log M 200 = α + β log N 200 . We find that the MRR we derive is consistent with both the IllustrisTNG and mini-Uchuu cosmological numerical simulations, with a slope ofβ≈ 1. We use the MRR we derived to estimate cluster masses from theGalWCat19catalog, which we then use to set constraints on Ωmandσ8. Utilizing the all-member MRR, we obtain constraints of Ωm= 0.31 0.03 + 0.04 andσ8= 0.82 0.04 + 0.05 , and utilizing the red member MRR only, we obtain Ωm= 0.31 0.03 + 0.04 andσ8= 0.81 0.04 + 0.05 . Our constraints on Ωmandσ8are consistent and very competitive with the Planck 2018 results. 
    more » « less
  4. 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
  5. Abstract Introduction: We present an extensive theoretical investigation of the electron impact excitation of doubly-ionized titanium (Ti III) to meet the needs of spectral analysis and plasma modeling. OBJECTIVES: The main objective of this work is to extend the currently scarce database of both structure and collision data for Ti III. METHODS: The calculation was performed in the close-coupling approximation using theB-splineR-matrix method. The multi-configuration Hartree–Fock method in combination withB-spline configuration interaction expansions and the non-orthogonal orbitals technique is employed for accurate descriptions of the target wave functions and adequate accounts of the various interactions between the target states. Relativistic effects are treated at the semi-relativistic Breit-Pauli approximation level. RESULTS: The present close-coupling expansion includes 138 fine-structure levels of Ti III belonging to the 3 d 2 , 4 s 2 , 4 s 4 p , 3 d 4 l ( l = 0 3 ), 3 d 5 l ( l = 0 3 ), 3 d 6 s , and 3 d 6 p configurations. Comprehensive sets of radiative and electron collisional data are reported for all of the possible transitions between the 138 fine-structure levels. Thermally averaged collision strengths are determined using a Maxwellian distribution for a wide range of temperatures from 10 2 K to 10 5 K. The accuracy of the calculated radiative parameters is validated by comparing with available values from the NIST database and previous literature. CONCLUSION: Given the lack of sufficient currently available experimental and theoretical data, the electron impact excitation cross sections of the Ti III fine-structure levels presented here are systematic, extensive, and internally consistent, thus making them suitable for many modeling applications. 
    more » « less