Modeled gravitational wave searches correlate the strain data with a bank of gravitational wave template waveforms to make detections of gravitational wave candidates, and these results are processed by downstream tools to calculate the likely sky location and distance of the source of the candidates. This is crucial for multi-messenger efforts, since it informs astronomers where to point their telescopes to facilitate electromagnetic follow-up of the gravitational wave candidates. We present a novel method to improve the low-latency results of the GstLAL gravitational wave search pipeline, and thus improving sky location estimates of low-latency candidates. This method involves ingesting the GstLAL low-latency results, and performing a small targeted hierarchical search to recover the candidates with more accurate parameters, in a medium-latency timescale (few seconds to five minutes). To test our method, we perform a GstLAL low-latency analysis on forty days of data from the third observing run of LIGO, Virgo, and KAGRA, and show that our method improves the GstLAL results by 5.38
We present the Stream Graph Navigator (SGN), a lightweight Python framework for building streaming data applications. In SGN, stream-processing pipelines are built by connecting computational components into directed acyclic graphs that run within an event loop. The time-series extension of the SGN library, SGN-TS, introduces signal processing methods to handle time series data. Together, SGN and SGN-TS provide the foundation for SGNL, a matched-filtering gravitational-wave search pipeline, and are being adopted by multiple projects across the low-latency gravitational-wave data analysis infrastructure as an extensible and maintainable framework for future gravitational-wave observations.
Low-latency gravitational-wave search pipelines provide early-warning alerts for multimessenger astrophysical transients. Current pipelines whiten the data stream using acausal, linear-phase filters, which require a look-ahead buffer that introduces several seconds of algorithmic latency. Eliminating this latency requires causal, minimum-phase whitening filters using only past data. However, operating causal filters under non-stationary noise is non-trivial: the drifting power spectral density must be tracked without degrading the matched-filter signal-to-noise ratio (SNR), filter updates must preserve the minimum-phase condition, and the altered phase response must be compensated to maintain sky-localization accuracy. In Paper I we introduced a gauge theoretic signal processing framework and showed that the minimum-phase connection on the manifold of power spectra provides a geometrically exact update rule for causal filters. Here we validate that framework numerically and operationally, demonstrating that parallel transport along this connection strictly preserves the minimum-phase property while exactly conserving the matched-filter SNR. We numerically certify the flatness of this connection, showing that the optimal filter is a path-independent state function of the instantaneous noise. Through an injection campaign on O3 data with 15,347 binary black hole signals across the LIGO-Virgo network, we confirm that this architecture preserves the detection sensitivity and inter-detector timing and phase accuracy of the linear-phase baseline. Implementing the framework in the production sgnl pipeline reduces whitening latency by 1.0 s (33
Gravitational-wave observations of merging binary neutron stars and black holes are now routinely made by detectors in the Advanced LIGO-Virgo-KAGRA network. Neutron star binary systems may also produce detectable electromagnetic and particle emission over times scales ranging from seconds to years. Real-time gravitational-wave searches play a central role in enabling time-critical electromagnetic and/or neutrino follow-up observations. During the fourth observing run (O4) of the Advanced LIGO-Virgo-KAGRA network, multiple real-time searches operated continuously to identify candidate gravitational-wave events and publicly disseminate information about these discoveries. Here, the performance and results of the GstLAL real-time analysis are reported. The analysis is designed to identify candidates with low latency, high detection efficiency, and sustained operational uptime over long observing periods. Across O4, it produced initial candidate uploads with a median latency of 15.8 s while maintaining an effective uptime of 98
We present the results of a search for gravitational-wave transients associated with core-collapse supernova SN 2023ixf, which was observed in the galaxy Messier 101 via optical emission on 2023 May 19, during the LIGO–Virgo–KAGRA 15th Engineering Run. We define a five-day on-source window during which an accompanying gravitational-wave signal may have occurred. No gravitational waves have been identified in data when at least two gravitational-wave observatories were operating, which covered ∼14% of this five-day window. We report the search detection efficiency for various possible gravitational-wave emission models. Considering the distance to M101 (6.7 Mpc), we derive constraints on the gravitational-wave emission mechanism of core-collapse supernovae across a broad frequency spectrum, ranging from 50 Hz to 2 kHz, where we assume the gravitational-wave emission occurred when coincident data are available in the on-source window. Considering an ellipsoid model for a rotating proto-neutron star, our search is sensitive to gravitational-wave energy 1 × 10 −4 M ⊙ c 2 and luminosity 2.6 × 10 −4 M ⊙ c 2 s −1 for a source emitting at 82 Hz. These constraints are around an order of magnitude more stringent than those obtained so far with gravitational-wave data. The constraint on the ellipticity of the proto-neutron star that is formed is as low as 1.08, at frequencies above 1200 Hz, surpassing past results.
Matched-filtering searches for gravitational-wave signals from compact binary mergers employ template banks which are a collection of modeled waveforms described by unique intrinsic parameters. We present two banks designed for low-latency and archive subsolar mass (SSM) searches in data from the fourth observing run of LIGO-Virgo-KAGRA, and demonstrate the efficacy of the banks via simulated signals. Further, we introduce several new methodological developments to the geometric placement algorithm [C. Hanna et al., Phys. Rev. D 108, 042003 (2023); C. Hanna, MANIFOLD (2024)], including a more stable metric-estimation technique and a boundary-padding procedure. Together, these enhancements enable the method to work for exceedingly low component masses necessary for SSM bank production, representing the first successful application of the MANIFOLD-method to lower-mass regions of the compact binary parameter space. The archive search bank contains a total of 3,452,006 templates, and covers a mass parameter space of 0.2 to 10Mo in the larger component and 0.2 to 1.0Mo in the smaller component, the spin parameter space of-0.9 to 0.9 for masses above 0.5Mo and-0.1 to 0.1 for masses below 0.5Mo, and the mass ratio parameter space of 1 to 10. The power spectral density used was from a week of the first half of the fourth observing run of Advanced LIGO, Advanced Virgo, and KAGRA, and the low frequency cutoff was set to 45 Hz with a maximum waveform duration of 128 secs. The bank simulations performed using SBank have shown that the banks presented in this paper have sufficient efficacy for use in their respective searches.
We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received with low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum-likelihood Non-imaging Transient Reconstruction and Temporal Search pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10(-3) Hz, we compute the GW-BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers.
Searches for gravitational waves from compact binary coalescences employ a process called matched filtering, in which gravitational wave strain data is cross-correlated against a bank of waveform templates. Data from every observing run of the LIGO, Virgo, and KAGRA collaboration is typically analyzed in this way twice, first in a low-latency mode in which gravitational wave candidates are identified in near-real time, and later in a high-latency mode. Such high-latency analyses have traditionally been considered more sensitive, since background data from the full observing run is available for assigning significance to all candidates, as well as more robust, since they do not need to worry about keeping up with live data. In this work, we present a novel technique to use the matched filtering data products from a low-latency analysis and re-process them by assigning significances in a high-latency way, effectively removing the need to perform matched filtering a second time. To demonstrate the efficacy of our method, we analyze 38 days of LIGO and Virgo data from the third observing run (O3) using the GstLAL pipeline, and show that our method is as sensitive and reliable as a traditional high-latency analysis. Since matched filtering represents the vast majority of computing time for a traditional analysis, our method greatly reduces the time and computational burden required to produce the same results as a traditional high-latency analysis. Consequently, it has already been adopted by GstLAL for the fourth observing run (O4) of the LIGO, Virgo, and KAGRA collaboration.
We present SGNL, a scalable, low-latency gravitational-wave search pipeline. It reimplements the core matched-filtering principles of the GstLAL pipeline within a modernized framework. The Stream Graph Navigator library, a lightweight Python streaming framework, replaces GstLAL's GStreamer infrastructure, simplifying pipeline construction and enabling flexible, modular graph design. The filtering core is reimplemented in PyTorch, allowing SGNL to leverage GPU acceleration for improved computational scalability. We describe the pipeline architecture and introduce a novel implementation of the Low-Latency Online Inspiral Detection algorithm in which components are pre-synchronized to reduce latency. Results from 40 days of data show that SGNL's event recovery and sensitivity are consistent with GstLAL's within statistical and systematic uncertainties. Notably, SGNL achieves a median latency of 4.7 seconds, compared to 9.0 seconds for GstLAL.
In this work, we present new methods implemented in the GstLAL offline gravitational wave search. These include a technique to reuse the matched filtering data products from a GstLAL online analysis, which hugely reduces the time and computational resources required to obtain offline results; a technique to combine these results with a separate search for heavier black hole mergers, enabling detections from a larger set of gravitational wave sources; changes to the likelihood ratio which increases the sensitivity of the analysis; and two separate changes to the background estimation, allowing more precise significance estimation of gravitational wave candidates. Some of these methods increase the sensitivity of the analysis, whereas others correct previous mis-estimations of sensitivity by eliminating false positives. These methods have been adopted for GstLAL's offline results during the fourth observing run of LIGO, Virgo, and KAGRA (O4). To test these new methods, we perform an offline analysis over one chunk of O3 data, lasting from May 12 19:36:42 UTC 2019 to May 21 14:45:08 UTC 2019, and compare it with previous GstLAL results over the same period of time. We show that cumulatively these methods afford around a 50
We present a subthreshold search for gravitational-wave inspirals from binary neutron stars using data from the first part of the fourth observing run of the LIGO-Virgo-KAGRA Collaboration. To enhance sensitivity to this targeted population, we incorporate a redshift-corrected population model informed by radio observations of Galactic double neutron star systems. The search identifies a significant trigger with a false-alarm rate of one per fifty years and a network signal-to-noise ratio of 9.7. This trigger was first reported in low-latency processing as S231109ci and subsequently listed in the GWTC-4.0 catalog as GW231109_235456, a subthreshold candidate. Accounting for a trials factor of five arising from four previous searches in GWTC-4.0 and this new search, the false-alarm rate of the candidate is approximately one per ten years. If the event is of astrophysical origin, the inferred source properties indicate component masses of 1.40 to 2.24 solar masses for the primary and 0.97 to 1.49 solar masses for the secondary, yielding a total mass of 2.95 solar masses with an uncertainty of plus 0.38 and minus 0.07 solar masses. The event is localized to a region of 450 square degrees enclosing ninety percent probability at a luminosity distance of 165 megaparsecs with an uncertainty of plus 70 and minus 69 megaparsecs. Assuming the signal arises from a binary neutron star merger, we estimate the local merger rate to lie between 53 and 342 per cubic gigaparsec per year.
In astronomy, we frequently face the decision problem: does this data contain a signal? Typically, a statistical approach is used, which requires a threshold. The choice of threshold presents a common challenge in settings where signals and noise must be delineated, but their distributions overlap. Gravitational-wave astronomy, which has gone from the first discovery to catalogues of hundreds of events in less than a decade, presents a fascinating case study. For signals from colliding compact objects, the field has evolved from a frequentist to a Bayesian methodology. However, the issue of choosing a threshold and validating noise contamination in a catalogue persists. Confusion and debate often arise due to the misapplication of statistical concepts, the complicated nature of the detection statistics, and the inclusion of astrophysical background models. We introduce Conformal Prediction (CP), a framework developed in Machine Learning to provide distribution-free uncertainty quantification to point predictors. We show that CP can be viewed as an extension of the traditional statistical frameworks whereby thresholds are calibrated such that the uncertainty intervals are statistically rigorous and the error rate can be validated. Moreover, we discuss how CP offers a framework to optimally build a meta-pipeline combining the outputs from multiple independent searches. We introduce CP with a toy cosmic-ray detector, which captures the salient features of most astrophysical search problems and allows us to demonstrate the features of CP in a simple context. We then apply the approach to a recent gravitational-wave Mock Data Challenge using multiple search algorithms for compact binary coalescence signals in interferometric gravitational-wave data. Finally, we conclude with a discussion on the future potential of the method for gravitational-wave astronomy.
We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers.
We introduce a new capability of the Neil Gehrels Swift Observatory, dubbed "continuous commanding," that achieves 10 s latency response time on orbit to unscheduled target-of-opportunity requests received on the ground. We show that this will allow Swift to respond to premerger (early-warning) gravitational-wave (GW) detections, rapidly slewing the Burst Alert Telescope (BAT) across the sky to place the GW origin in the BAT field of view at or before merger time. This will dramatically increase the GW/gamma-ray burst (GRB) codetection rate and enable prompt arcminute localization of a neutron star merger. We simulate the full Swift response to a GW early-warning alert, including input sky maps produced at different early-warning times, a complete model of the Swift attitude control system, and a full accounting of the latency between the GW detectors and the spacecraft. 60 s of early warning can double the rate of a prompt GRB detection with arcminute localization, and 140 s guarantees observation anywhere on the unocculted sky, even with localization areas >> 1000 deg(2). While 140 s is beyond current GW detector sensitivities, 30-70 s is achievable today. We show that the detection yield is now limited by the latency of LIGO/Virgo cyberinfrastructure and motivate a focus on its reduction. Continuous commanding has been integrated as a general capability of Swift, significantly increasing its versatility in response to the growing demands of time-domain astrophysics. We demonstrate this potential on an externally triggered fast radio burst (FRB), slewing 81 degrees across the sky, and collecting X-ray and UV photons from the source position <150 s after the trigger was received from the Canadian Hydrogen Intensity Mapping Experiment, thereby setting the earliest and deepest such constraints on high-energy activity from nonrepeating FRBs. The Swift Team invites the community to consider and propose novel scientific applications of ultra-low-latency UV, X-ray, and gamma-ray observations.
Leveraging the features of the GstLAL pipeline, we present the results of a matched filtering search for asymmetric binary black hole systems with heavily misaligned spins in LIGO and Virgo data taken during the third observing run. Our target systems show strong imprints of precession whereas current searches have nonoptimal sensitivity in detecting them. After measuring the sensitivity improvement brought by our search over standard spin-aligned searches, we report the detection of 30 gravitational wave events already discovered in the latest version of the Gravitational Wave Transient Catalog. However, we do not find any additional significant gravitational wave candidates. Our results allow us to place an upper limit of R_{90%}=0.28_{-0.04}^{+0.33} Gpc^{-3} yr^{-1} on the merger rate of a hypothetical subpopulation of asymmetric, heavily precessing signals, not identified by other searches. Since our upper limit is consistent with the latest rate estimates from the LIGO-Virgo-KAGRA Collaboration, our findings rule out the existence of a yet-to-be-discovered population of precessing binaries.
Precession in binary black holes (BBH) is caused by the failure of the black hole spins to be aligned and its study can open up new perspectives in gravitational wave astronomy, providing, among other advancements, a precise measure of distance and an accurate characterization of the BBH spins. However, detecting precessing signals is a highly nontrivial task, as standard matched filtering pipelines for gravitational wave searches are built on many assumptions that do not hold in the precessing case. This work details the upgrades made to the GstLAL pipeline to facilitate the search for precessing BBH signals. The implemented changes in the search statistics and in the signal consistency test are then described in detail. The performance of the upgraded pipeline is evaluated through two extensive searches of precessing signals, targeting two different regions in the mass space, and the consistency of the results is examined. Additionally, the benefits of the upgrades are assessed by comparing the sensitive volume of the precessing searches with two corresponding traditional aligned-spin searches. While no significant sensitivity improvement is observed for precessing binaries with mass ratio q less than or similar to 6, a volume increase of up to 100% is attainable for heavily asymmetric systems with largely misaligned spins. Furthermore, our findings suggest that the primary cause of degraded performance in an aligned-spin search targeting precessing signals is not a poor signal-to-noise-ratio recovery but rather the failure of the xi 2 2 signal- consistency test. Our work paves the way for a large-scale search for precessing signals, which could potentially result in exciting future detections.
GstLAL is a stream-based matched-filtering search pipeline aiming at the prompt discovery of gravitational waves from compact binary coalescences such as the mergers of black holes and neutron stars. Over the past three observation runs by the LIGO, Virgo, and KAGRA (LVK) collaboration, the GstLAL search pipeline has participated in several tens of gravitational wave discoveries. The fourth observing run (O4) is set to begin in May 2023 and is expected to see the discovery of many new and interesting gravitational wave signals which will inform our understanding of astrophysics and cosmology. We describe the current configuration of the GstLAL low-latency search and show its readiness for the upcoming observation run by presenting its performance on a mock data challenge. The mock data challenge includes 40 days of LIGO Hanford, LIGO Livingston, and Virgo strain data along with an injection campaign in order to fully characterize the performance of the search. We find an improved performance in terms of detection rate and significance estimation as compared to that observed in the O3 online analysis. The improvements are attributed to several incremental advances in the likelihood ratio ranking statistic computation and the method of background estimation.
The Soft X-ray Imager (SXI) is part of the scientific payload of the Solar wind Magnetosphere Ionosphere Link Explorer (SMILE) mission. SMILE is a joint science mission between the European Space Agency (ESA) and the Chinese Academy of Sciences (CAS) and is due for launch in 2025. SXI is a compact X-ray telescope with a wide field-of-view (FOV) capable of encompassing large portions of Earth's magnetosphere from the vantage point of the SMILE orbit. SXI is sensitive to the soft X-rays produced by the Solar Wind Charge eXchange (SWCX) process produced when heavy ions of solar wind origin interact with neutral particles in Earth's exosphere. SWCX provides a mechanism for boundary detection within the magnetosphere, such as the position of Earth's magnetopause, because the solar wind heavy ions have a very low density in regions of closed magnetic field lines. The sensitivity of the SXI is such that it can potentially track movements of the magnetopause on timescales of a few minutes and the orbit of SMILE will enable such movements to be tracked for segments lasting many hours. SXI is led by the University of Leicester in the United Kingdom (UK) with collaborating organisations on hardware, software and science support within the UK, Europe, China and the United States.
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.
Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level.