Search

We develop a robust and self-consistent framework to extract and classify gravitational wave candidates from noisy data, for the purpose of assisting in real-time multi-messenger follow-ups during LIGO-Virgo-KAGRA's fourth observing run~(O4). Our formalism implements several improvements to the low latency calculation of the probability of astrophysical origin~(\PASTRO{}), so as to correctly account for various factors such as the sensitivity change between observing runs, and the deviation of the recovered template waveform from the true gravitational wave signal that can strongly bias said calculation. We demonstrate the high accuracy with which our new formalism recovers and classifies gravitational wave triggers, by analyzing replay data from previous observing runs injected with simulated sources of different categories. We show that these improvements enable the correct identification of the majority of simulated sources, many of which would have otherwise been misclassified. We carry out the aforementioned analysis by implementing our formalism through the \GSTLAL{} search pipeline even though it can be used in conjunction with potentially any matched filtering pipeline. Armed with robust and self-consistent \PASTRO{} values, the \GSTLAL{} pipeline can be expected to provide accurate source classification information for assisting in multi-messenger follow-up observations to gravitational wave alerts sent out during O4.

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., 19 pages, 21 figures

Starting from May 2023, the LIGO Scientific, Virgo and KAGRA Collaboration is planning to conduct the fourth observing run with improved detector sensitivities and an expanded detector network including KAGRA. Accordingly, it is vital to optimize the detection algorithm of low-latency search pipelines, increasing their sensitivities to gravitational waves from compact binary coalescences. In this work, we discuss several new features developed for ranking statistics of GstLAL-based inspiral pipeline, which mainly consist of: the signal contamination removal, the bank-$\xi^2$ incorporation, the upgraded $\rho-\xi^2$ signal model and the integration of KAGRA. An injection study demonstrates that these new features improve the pipeline's sensitivity by approximately 15% to 20%, paving the way to further multi-messenger observations during the upcoming observing run., Comment: 13pages, 6figures

Template banks containing gravitational wave (GW) waveforms are essential for matched-filtering GW search pipelines. We describe the generation method, the design, and validation of the template bank used by the GstLAL-based inspiral pipeline to analyze data from the fourth observing run of LIGO scientific, Virgo, and KAGRA collaboration. This paper presents a template bank containing $1.8 \times 10^6$ templates that include merging neutron star - neutron star, neutron star - black hole, and black hole - black hole systems up to a total mass of $400$ $M_\odot$. Motivated by observations, component masses below $3$ $M_\odot$ have dimensionless spins ranging between $\pm 0.05$, while component masses between $3$ to $200$ $M_\odot$ have dimensionless spins ranging between $\pm 0.99$, where we assume spin-aligned systems. The low-frequency cutoff is $15$ Hz. The templates are placed in the parameter space according to the metric via a binary tree approach which took $\mathcal{O}\left(10\right)$ minutes when jobs were parallelized. The template bank generated with this method has a $98\%$ match or higher for $90\%$ of the injections, thus being as effective as the template placement method used for the previous observation runs. The volumes of the templates are computed prior to template placement and the nearby templates have similar volumes in the coordinate space, henceforth, enabling a more efficient and less biased implementation of population models. SVD sorting of the O4 template bank has been renewed to use post-Newtonian phase terms, which improved the computational efficiency of SVD by nearly $4 \sim 5$ times as compared to conventional SVD sorting schemes. Template banks and searches focusing on the sub-solar mass parameter space and intermediate-mass black hole parameter space are conducted separately.

We demonstrate a new geometric method for fast template placement for searches for gravitational waves from the inspiral, merger and ringdown of compact binaries. The method is based on a binary tree decomposition of the template bank parameter space into non-overlapping hypercubes. We use a numerical approximation of the signal overlap metric at the center of each hypercube to estimate the number of templates required to cover the hypercube and determine whether to further split the hypercube. As long as the expected number of templates in a given cube is greater than a given threshold, we split the cube along its longest edge according to the metric. When the expected number of templates in a given hypercube drops below this threshold, the splitting stops and a template is placed at the center of the hypercube. Using this method, we generate aligned-spin template banks covering the mass range suitable for a search of Advanced LIGO data. The aligned-spin bank required ~24 CPU-hours and produced 2 million templates. In general, we find that other methods, namely stochastic placement, produces a more strictly bounded loss in match between waveforms, with the same minimal match between waveforms requiring about twice as many templates with our proposed algorithm. Though we note that the average match is higher, which would lead to a higher detection efficiency. Our primary motivation is not to strictly minimize the number of templates with this algorithm, but rather to produce a bank with useful geometric properties in the physical parameter space coordinates. Such properties are useful for population modeling and parameter estimation.

Observations of the properties of multiple coalescing neutron stars will simultaneously provide insight into neutron star mass and spin distribution, the neutron star merger rate, and the nuclear equation of state. Not all merging binaries containing neutron stars are expected to be identical. Plausible sources of diversity in these coalescing binaries can arise from a broad or multi-peaked NS mass distribution; the effect of different and extreme NS natal spins; the possibility of NS-BH mergers; or even the possibility of phase transitions, allowing for NS with similar mass but strongly divergent radius. In this work, we provide a concrete algorithm to combine all information obtained from GW measurements into a joint constraint on the NS merger rate, the distribution of NS properties, and the nuclear equation of state. Using a concrete example, we show how biased mass distribution inferences can significantly impact the recovered equation of state, even in the small-$N$ limit. With the same concrete example, we show how small-$N$ observations could identify a bimodal mass and spin distribution for merging NS simultaneously with the EOS. Our concordance approach can be immediately generalized to incorporate other observational constraints.

After their successful first observing run (September 12, 2015 - January 12, 2016), the Advanced LIGO detectors were upgraded to increase their sensitivity for the second observing run (November 30, 2016 - August 26, 2017). The Advanced Virgo detector joined the second observing run on August 1, 2017. We discuss the updates that happened during this period in the GstLAL-based inspiral pipeline, which is used to detect gravitational waves from the coalescence of compact binaries both in low latency and an offline configuration. These updates include deployment of a zero-latency whitening filter to reduce the over-all latency of the pipeline by up to 32 seconds, incorporation of the Virgo data stream in the analysis, introduction of a single-detector search to analyze data from the periods when only one of the detectors is running, addition of new parameters to the likelihood ratio ranking statistic, increase in the parameter space of the search, and introduction of a template mass-dependent glitch-excision thresholding method., 12 pages, 7 figures, to be submitted to Phys. Rev. D, comments welcome