skip to main content


This content will become publicly available on July 1, 2024

Title: pygwb: A Python-based Library for Gravitational-wave Background Searches
Abstract The collection of gravitational waves (GWs) that are either too weak or too numerous to be individually resolved is commonly referred to as the gravitational-wave background (GWB). A confident detection and model-driven characterization of such a signal will provide invaluable information about the evolution of the universe and the population of GW sources within it. We present a new, user-friendly, Python-based package for GW data analysis to search for an isotropic GWB in ground-based interferometer data. We employ cross-correlation spectra of GW detector pairs to construct an optimal estimator of the Gaussian and isotropic GWB, and Bayesian parameter estimation to constrain GWB models. The modularity and clarity of the code allow for both a shallow learning curve and flexibility in adjusting the analysis to one’s own needs. We describe the individual modules that make up pygwb , following the traditional steps of stochastic analyses carried out within the LIGO, Virgo, and KAGRA Collaboration. We then describe the built-in pipeline that combines the different modules and validate it with both mock data and real GW data from the O3 Advanced LIGO and Virgo observing run. We successfully recover all mock data injections and reproduce published results.  more » « less
Award ID(s):
2207758 2207945 1912594
NSF-PAR ID:
10439738
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
952
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
25
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The accretion disks of active galactic nuclei (AGNs) are promising locations for the merger of compact objects detected by gravitational wave (GW) observatories. Embedded within a baryon-rich, high-density environment, mergers within AGNs are the only GW channel where an electromagnetic (EM) counterpart must occur (whether detectable or not). Considering AGNs with unusual flaring activity observed by the Zwicky Transient Facility (ZTF), we describe a search for candidate EM counterparts to binary black hole (BBH) mergers detected by LIGO/Virgo in O3. After removing probable false positives, we find nine candidate counterparts to BBH mergers during O3 (seven in O3a, two in O3b) with ap-value of 0.0019. Based on ZTF sky coverage, AGN geometry, and merger geometry, we expect ≈3(NBBH/83)(fAGN/0.5) potentially detectable EM counterparts from O3, whereNBBHis the total number of observed BBH mergers andfAGNis the fraction originating in AGNs. Further modeling of breakout and flaring phenomena in AGN disks is required to reduce our false-positive rate. Two of the events are also associated with mergers with total masses >100M, which is the expected rate for O3 if hierarchical (large-mass) mergers occur in the AGN channel. Candidate EM counterparts in future GW observing runs can be better constrained by coverage of the Southern sky as well as spectral monitoring of unusual AGN flaring events in LIGO/Virgo alert volumes. A future set of reliable AGN EM counterparts to BBH mergers will yield an independent means of measuring cosmic expansion (H0) as a function of redshift.

     
    more » « less
  2. Multimessenger searches for binary neutron star (BNS) and neutron star-black hole (NSBH) mergers are currently one of the most exciting areas of astronomy. The search for joint electromagnetic and neutrino counterparts to gravitational wave (GW)s has resumed with ALIGO’s, AdVirgo’s and KAGRA’s fourth observing run (O4). To support this effort, public semiautomated data products are sent in near real-time and include localization and source properties to guide complementary observations. In preparation for O4, we have conducted a study using a simulated population of compact binaries and a mock data challenge (MDC) in the form of a real-time replay to optimize and profile the software infrastructure and scientific deliverables. End-toend performance was tested, including data ingestion, running online search pipelines, performing annotations, and issuing alerts to the astrophysics community. We present an overview of the low-latency infrastructure and the performance of the data products that are now being released during O4 based on the MDC. We report the expected median latency for the preliminary alert of full bandwidth searches (29.5 s) and show consistency and accuracy of released data products using the MDC. We report the expected median latency for triggers from early warning searches (−3.1 s), which are new in O4 and target neutron star mergers during inspiral phase. This paper provides a performance overview for LIGO-Virgo-KAGRA (LVK) low-latency alert infrastructure and data products using theMDCand serves as a useful reference for the interpretation of O4 detections.

     
    more » « less
  3. null (Ed.)
    Advanced LIGO and Advanced Virgo are monitoring the sky and collecting gravitational-wave strain data with sufficient sensitivity to detect signals routinely. In this paper we describe the data recorded by these instruments during their first and second observing runs. The main data products are gravitational-wave strain time series sampled at 16384 Hz. The datasets that include this strain measurement can be freely accessed through the Gravitational Wave Open Science Center at http://gw-openscience.org, together with data-quality information essential for the analysis of LIGO and Virgo data, documentation, tutorials, and supporting software. 
    more » « less
  4. Abstract

    The LIGO/Virgo collaboration published the catalogs GWTC-1, GWTC-2.1, and GWTC-3 containing candidate gravitational-wave (GW) events detected during its runs O1, O2, and O3. These GW events can be possible sites of neutrino emission. In this paper, we present a search for neutrino counterparts of 90 GW candidates using IceCube DeepCore, the low-energy infill array of the IceCube Neutrino Observatory. The search is conducted using an unbinned maximum likelihood method, within a time window of 1000 s, and uses the spatial and timing information from the GW events. The neutrinos used for the search have energies ranging from a few GeV to several tens of TeV. We do not find any significant emission of neutrinos, and place upper limits on the flux and the isotropic-equivalent energy emitted in low-energy neutrinos. We also conduct a binomial test to search for source populations potentially contributing to neutrino emission. We report a nondetection of a significant neutrino-source population with this test.

     
    more » « less
  5. Abstract

    We present a new class of AI models for the detection of quasi-circular, spinning, non-precessing binary black hole mergers whose waveforms include the higher order gravitational wave modes(,|m|)={(2,2),(2,1),(3,3),(3,2),(4,4)}, and mode mixing effects in the=3,|m|=2harmonics. These AI models combine hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational waves; and graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a three detector network encompassing the Advanced LIGO and Virgo detectors. We first trained these spatiotemporal-graph AI models using synthetic noise, using 1.2 million modeled waveforms to densely sample this signal manifold, within 1.7 h using 256 NVIDIA A100 GPUs in the Polaris supercomputer at the Argonne Leadership Computing Facility. This distributed training approach exhibited optimal classification performance, and strong scaling up to 512 NVIDIA A100 GPUs. With these AI ensembles we processed data from a three detector network, and found that an ensemble of 4 AI models achieves state-of-the-art performance for signal detection, and reports two misclassifications for every decade of searched data. We distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, and completed the processing of a decade of gravitational wave data from a three detector network within 3.5 h. Finally, we fine-tuned these AI ensembles to process the entire month of February 2020, which is part of the O3b LIGO/Virgo observation run, and found 6 gravitational waves, concurrently identified in Advanced LIGO and Advanced Virgo data, and zero false positives. This analysis was completed in one hour using one NVIDIA A100 GPU.

     
    more » « less