Search

The Laser Interferometer Space Antenna (LISA) will observe gravitational-wave signals from a wide range of sources, including massive black hole binaries. Although numerous techniques have been developed to perform Bayesian inference for LISA, they are often computationally expensive; analyses often take at least $\sim 1$ month on a single CPU, even when using accelerated techniques. Not only does this make it difficult to concurrently analyse more than one gravitational-wave signal, it also makes it challenging to rapidly produce parameter estimates for possible electromagnetic follow-up campaigns. simple-pe was recently developed to produce rapid parameter estimates for gravitational-wave signals observed with ground-based gravitational-wave detectors. In this work, we extend simple-pe to produce rapid parameter estimates for LISA sources, including the effects of higher order multipole moments. We show that simple-pe infers the source properties of massive black hole binaries in zero-noise at least $\sim 100\times$ faster than existing techniques; $\sim 12$ hours on a single CPU. We further demonstrate that simple-pe can be applied before existing Bayesian techniques to mitigate biases in multi-modal parameter estimation analyses of MBHBs., Comment: 28 pages, 8 figures

Gravitational wave searches are crucial for studying compact sources like neutron stars and black holes. Many sensitive modeled searches use matched filtering to compare gravitational strain data to a set of waveform models known as template banks. We introduce a new stochastic placement method for constructing template banks, offering efficiency and flexibility to handle arbitrary parameter spaces, including orbital eccentricity, tidal deformability, and other extrinsic parameters. This method can be computationally limited by the ability to compare proposal templates with the accepted templates in the bank. To alleviate this computational load, we introduce the use of inner product inequalities to reduce the number of required comparisons. We also introduce a novel application of Gaussian Kernel Density Estimation to enhance waveform coverage in sparser regions. Our approach has been employed to search for eccentric binary neutron stars, low-mass neutron stars, primordial black holes, supermassive black hole binaries. We demonstrate that our method produces self-consistent banks that recover the required minimum fraction of signals. For common parameter spaces, our method shows comparable computational performance and similar template bank sizes to geometric placement methods and stochastic methods, while easily extending to higher-dimensional problems. The time to run a search exceeds the time to generate the bank by a factor of $\mathcal{O}(10^5)$ for dedicated template banks, such as geometric, mass-only stochastic, and aligned spin cases, $\mathcal{O}(10^4)$ for eccentric and $\mathcal{O}(10^3)$ for the tidal deformable bank. With the advent of efficient template bank generation, the primary area for improvement is developing more efficient search methodologies., Comment: 14 Pages, 4 Figures, 1 Table

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., Comment: 18 pages, 20 figures, accepted in Phys. Rev. D

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents

Multi-messenger searches for BNS and NSBH mergers are currently one of the most exciting areas of astronomy. The search for joint electromagnetic and neutrino counterparts to GWs has resumed with O4. To support this effort, public semi-automated 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 MDC in the form of a real-time replay to optimize and profile the software infrastructure and scientific deliverables. End-to-end 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.5s) and show consistency and accuracy of released data products using the MDC. For the first time, we report the expected median latency for triggers from early warning searches (-3.1s), which are new in O4 and target neutron star mergers during inspiral phase. This paper provides a performance overview for LVK low-latency alert infrastructure and data products using the MDC and serves as a useful reference for the interpretation of O4 detections.

The Laser Interferometer Space Antenna (LISA), due for launch in the mid 2030s, is expected to observe gravitational waves (GW)s from merging massive black hole binaries (MBHB)s. These signals can last from days to months, depending on the masses of the black holes, and are expected to be observed with high signal to noise ratios (SNR)s out to high redshifts. We have adapted the PyCBC software package to enable a template bank search and inference of GWs from MBHBs. The pipeline is tested on the LISA data challenge (LDC)'s Challenge 2a (\enquote{Sangria}), which contains MBHBs and thousands of galactic binaries (GBs) in simulated instrumental LISA noise. Our search identifies all 6 MBHB signals with more than $92\%$ of the optimal SNR. The subsequent parameter inference step recovers the masses and spins within their $90\%$ confidence interval. Sky position parameters have 8 high likelihood modes which are recovered but often our posteriors favour the incorrect sky mode. We observe that the addition of GBs biases the parameter recovery of masses and spins away from the injected values, reinforcing the need for a global fit pipeline which will simultaneously fit the parameters of the GB signals before estimating the parameters of MBHBs., Comment: 27 pages, 10 figures and 6 tables

Along their path from source to observer, gravitational waves may be gravitationally lensed by massive objects. This results in distortions of the observed signal which can be used to extract new information about fundamental physics, astrophysics, and cosmology. Searches for these distortions amongst the observed signals from the current detector network have already been carried out, though there have as yet been no confident detections. However, predictions of the observation rate of lensing suggest detection in the future is a realistic possibility. Therefore, preparations need to be made to thoroughly investigate the candidate lensed signals. In this work, we present some of the follow-up analyses and strategies that could be applied to assess the significance of such events and ascertain what information may be extracted about the lens-source system from such candidate signals by applying them to a number of O3 candidate events, even if these signals did not yield a high significance for any of the lensing hypotheses. For strongly-lensed candidates, we verify their significance using a background of simulated unlensed events and statistics computed from lensing catalogs. We also look for potential electromagnetic counterparts. In addition, we analyse in detail a candidate for a strongly-lensed sub-threshold counterpart that is identified by a new method. For microlensing candidates, we perform model selection using a number of lens models to investigate our ability to determine the mass density profile of the lens and constrain the lens parameters. We also look for millilensing signatures in one of the lensed candidates. Applying these additional analyses does not lead to any additional evidence for lensing in the candidates that have been examined. However, it does provide important insight into potential avenues to deal with high-significance candidates in future observations., Comment: 29 pages, 27 figures

The measurement of eccentricity would provide strong constraints on the formation channel of stellar-mass binary black holes. However, current ground-based gravitational wave detectors will, in most cases, not be able to measure eccentricity due to orbital circularization. Space-based observatories, in contrast, can determine binary eccentricity at 0.01Hz to $e_{0.01}\gtrsim\mathcal{O}(10^{-4}) $. Directly observing stellar-mass binary black holes with space-based observatories remains a challenging problem. However, observing such systems with ground-based detectors allows the possibility to identify the same signal in archival data from space-based observatories in the years previous. Since ground-based detectors provide little constraints on eccentricity, including eccentricity in the archival search will increase the required number of filter waveforms for the archival search by 5 orders of magnitudes [from $\sim \mathcal{O}(10^3)$ to $\sim \mathcal{O}(10^8)$], and will correspondingly need $ \sim8\times10^5 $ core hours (and $ \sim 10^5$ GB of memory), even for a mild upper limit on eccentricity of $0.1$. In this work, we have constructed the first template bank for an archival search of space-based gravitational wave detectors, including eccentricity. We have demonstrated that even though the inclusion of eccentricity brings extra computational burden, an archival search including eccentricity will be feasible in the time frame of planned space-based observatories, and will provide strong constraints on the eccentricities of stellar-mass binary black holes., Comment: 8 pages, 4 figures, comments welcome

Gravitational-wave signals from compact binary coalescences are most efficiently identified through matched filter searches, which match the data against a pre-generated bank of gravitational-wave templates. Although different techniques for performing the matched filter, as well as generating the template bank, exist, currently all modelled gravitational-wave searches use templates that restrict the component spins to be aligned (or anti-aligned) with the orbital angular momentum. This means that current searches are less sensitive to gravitational-wave signals generated from binaries with generic spins (precessing), suggesting that, potentially, a significant fraction of signals may remain undetected. In this work we introduce a matched filter search that is sensitive to signals generated from precessing binaries and can realistically be used during a gravitational-wave observing run. We take advantage of the fact that a gravitational-wave signal from a precessing binary can be decomposed into a power series of five harmonics, to show that a generic-spin template bank, which is only $\sim 3\times$ larger than existing aligned-spin banks, is needed to increase our sensitive volume by $\sim 100\%$ for neutron star black hole binaries with total mass larger than $17.5\, M_{\odot}$ and in-plane spins $>0.67$. In fact, our generic spin search performs as well as existing aligned-spin searches for neutron star black hole signals with insignificant in-plane spins, but improves sensitivity by $\sim60\%$ on average across the full generic spin parameter space. We anticipate that this improved technique will identify significantly more gravitational-wave signals, and, ultimately, help shed light on the unknown spin distribution of binaries in the universe., Comment: 23 pages, 12 figures. For data release, see https://icg-gravwaves.github.io/precessing_search_paper/

Data recorded by gravitational wave detectors includes many non-astrophysical transient noise bursts, the most common of which is caused by scattered-light within the detectors. These so-called ``glitches'' in the data impact the ability to both observe and characterize incoming gravitational wave signals. In this work we use a scattered-light glitch waveform model to identify and characterize scattered-light glitches in a representative stretch of gravitational wave data. We identify $2749$ scattered-light glitches in $5.96$ days of LIGO-Hanford data and $1306$ glitches in $5.93$ days of LIGO-Livingston data taken from the third LIGO-Virgo observing run. By subtracting identified scattered-light glitches we demonstrate an increase in the sensitive volume of the gravitational wave search for binary black hole signals by $\sim1\%$., Comment: 31 pages + acknowledgements, references and appendices, 13 figures

Standard detection and analysis techniques for transient gravitational waves make the assumption that detector data contains, at most, one signal at any time. As detectors improve in sensitivity, this assumption will no longer be valid. In this paper we examine how current search techniques for transient gravitational waves will behave under the presence of more than one signal. We perform searches on data sets containing time-overlapping compact binary coalescences. This includes a modelled, matched filter search (PyCBC), and an unmodelled coherent search, coherent WaveBurst (cWB). Both of these searches are used by the LIGO-Virgo-KAGRA collaboration. We find that both searches are capable of identifying both signals correctly when the signals are dissimilar in merger time, $|\Delta t_c| \geq 1$ second, with PyCBC losing only $1\%$ of signals for overlapping binary black hole mergers. Both pipelines can find signal pairings within the region $|\Delta t_c| < 1$ second. However, clustering routines in the pipelines will cause only one of the two signals to be recovered, as such the efficiency is reduced. Within this region, we find that cWB can identify both signals. We also find that matched filter searches can be modified to provide estimates of the correct parameters for each signal., Comment: 17 pages, 12 figures, 4 tables

Theories of dark energy that affect the speed of gravitational waves $c_{\rm GW}$ on cosmological scales naturally lead to a frequency-dependent transition of that speed close to the LIGO/Virgo/KAGRA (LVK) band. While observations such as GW170817 assure us that $c_{\rm GW}$ is extremely close to the speed of light in the LVK band, a frequency-dependent transition below the LVK band is a smoking-gun signal for large classes of dynamical dark energy theories. Here we discuss 1) how the remnants of such a transition can be constrained with observations in the LVK band, 2) what signatures are associated with such a transition in the LISA band, and 3) how joint observations in the LVK and LISA bands allow us to place tight constraints on this transition and the underlying theories. We find that deviations of $c_{\rm GW}$ can be constrained down to a level of $\sim 10^{-17}$ in the LVK ${\textit and}$ LISA bands even for mild frequency-dependence, much stronger than existing bounds for frequency-independent $c_{\rm GW} \neq c$. We use the strain data from GW170817 to bound the deviation of $c_{\rm GW}$ to be less than $10^{-17}$ at 100 Hz and less than $10^{-18}$ at 500 Hz. We also identify a particularly interesting type of transition in between the LVK and LISA bands and show how multi-band observations can constrain this further. Finally, we discuss what these current and forecasted constraints imply for the underlying dark energy theories., Comment: 12 pages + appendices and references, 5 figures; v3: minor clarifications added, matches published version

Current matched-filter searches for gravitational waves from binary black-hole mergers compare the calibrated detector data to waveform templates that omit the higher-order mode content of the signals predicted by General Relativity. However, higher-order emission modes become important for highly inclined asymmetric sources with masses above $\simeq 100 M_\odot$, causing current searches to be ill-suited at detecting them. We present a new gravitational-wave search that implements templates including higher-order modes, adapted signal-glitch discriminators, and trigger-ranking statistics to specifically target signals displaying strong higher modes, corresponding to nearly edge-on sources with total redshifted masses in the intermediate-mass black-hole range $M_T(1+z) \in (100,500) M_\odot$ and mass-ratios $q\in(1,10)$. Our search shows a volumetric sensitivity gain of up to $450\%$ to these signals compared to existing searches omitting higher-order modes. We deploy our search on public data from the third observing run of Advanced LIGO. While we find no statistically significant candidates beyond those already reported elsewhere, our search sets the stage to search for higher-mode rich signals in future observing runs. The efficient detection of such signals is crucial to performing detailed tests of General Relativity, observing strong-field phenomena, and maximizing the chances of observing the yet uncharted realm of intermediate-mass black hole binaries., Comment: 15 pages, 8 figures, matches published version

Gravitational-wave observations of compact binary coalescences are allowing us to see black holes and neutron stars further into the universe and recent results represent the most sensitive searches for compact objects ever undertaken. Most searches for gravitational waves from compact binary coalescence currently rely on detecting coincident triggers from multiple detectors. In this paper, we describe a new method for extrapolating significance of single-detector signals beyond the live-time of the analysis. Using this method, we can recover loud signals which only triggered in a single detector. We demonstrate this method in a search of O3 data, and recover seven single-detector events with a false alarm rate less than two per year. These were the same events as discovered in the GWTC-2.1 and GWTC-3 searches in a single detector, and all but one event from 3-OGC and 4-OGC. Through a campaign of injected signals, we estimate that the total time--volume sensitivity increases by a factor of up to $1.20 \pm 0.02$ at a false alarm rate of one per two years compared to completely ignoring single-detector events., Comment: 25 pages, 9 figures, 6 tables

The sensitivity of gravitational wave searches is reduced by the presence of non-Gaussian noise in the detector data. These non-Gaussianities often match well with the template waveforms used in matched filter searches, and require signal-consistency tests to distinguish them from astrophysical signals. However, empirically tuning these tests for maximum efficacy is time consuming and limits the complexity of these tests. In this work we demonstrate a framework to use machine-learning techniques to automatically tune signal-consistency tests. We implement a new $\chi^2$ signal-consistency test targeting the large population of noise found in searches for intermediate mass black hole binaries, training the new test using the framework set out in this paper. We find that this method effectively trains a complex model to down-weight the noise, while leaving the signal population relatively unaffected. This improves the sensitivity of the search by $\sim 11\%$ for signals with masses $> 300 M_\odot$. In the future this framework could be used to implement new tests in any of the commonly used matched-filter search algorithms, further improving the sensitivity of our searches., Comment: 10 pages, 5 figures. Supplementary data: https://icg-gravwaves.github.io/chisqnet/ . Version accepted for publication in PRD. Various updates made during review process. Typos corrected

The next generation of ground-based gravitational-wave detectors will observe coalescences of black holes and neutron stars throughout the cosmos, thousands of them with exceptional fidelity. The Science Book is the result of a 3-year effort to study the science capabilities of networks of next generation detectors. Such networks would make it possible to address unsolved problems in numerous areas of physics and astronomy, from Cosmology to Beyond the Standard Model of particle physics, and how they could provide insights into workings of strongly gravitating systems, astrophysics of compact objects and the nature of dense matter. It is inevitable that observatories of such depth and finesse will make new discoveries inaccessible to other windows of observation. In addition to laying out the rich science potential of the next generation of detectors, this report provides specific science targets in five different areas in physics and astronomy and the sensitivity requirements to accomplish those science goals. This report is the second in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book (this report), iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network., Comment: 69 pages, 18 figures

SkyPy is an open-source Python package for simulating the astrophysical sky. It comprises a library of physical and empirical models across a range of observables and a command-line script to run end-to-end simulations. The library provides functions that sample realisations of sources and their associated properties from probability distributions. Simulation pipelines are constructed from these models using a YAML-based configuration syntax, while task scheduling and data dependencies are handled internally and the modular design allows users to interface with external software. SkyPy is developed and maintained by a diverse community of domain experts with a focus on software sustainability and interoperability. By fostering development, it provides a framework for correlated simulations of a range of cosmological probes including galaxy populations, large scale structure, the cosmic microwave background, supernovae and gravitational waves. Version 0.4 implements functions that model various properties of galaxies including luminosity functions, redshift distributions and optical photometry from spectral energy distribution templates. Future releases will provide additional modules, for example, to simulate populations of dark matter halos and model the galaxy-halo connection, making use of existing software packages from the astrophysics community where appropriate., Comment: Published by JOSS. The package is available at https://github.com/skypyproject/skypy. Comments, issues and pull requests are welcome

Recent gravitational-wave observations have begun to constrain the internal physics of neutron stars. However, current detection searches for neutron star systems assume that potential neutron stars are low-mass black holes, ignoring any affect on the gravitational-wave signal due to the internal neutron-star physics. We wish to create a template bank of binary neutron star waveforms including the effect of tidal deformability. However, we find that the Fisher matrix, which is commonly used to approximate match calculations when placing template banks, is unsuitable to predict the match between two binary neutron star waveforms. We find that the Fisher matrix can predict errors on the mismatch that are larger than $100\%$ when attempting to identify waveforms with a match of $0.97$. We explore the regime in which the Fisher matrix cannot be trusted and examine why it breaks down. We demonstrate that including higher-order terms in the Taylor series expansion of the match can reliably compute matches for these examples, but that it is prohibitively computationally expensive to do so. Finally, we demonstrate that stochastic placement can still be used to construct a template bank of tidally deformed neutron-star waveforms., Comment: 11 pages, 5 figures

The third observing run of Advanced LIGO and Advanced Virgo took place between April 2019 and March 2020 and resulted in dozens of gravitational-wave candidates, many of which are now published as confident detections. A crucial requirement of the third observing run has been the rapid identification and public reporting of compact binary mergers, which enabled massive followup observation campaigns with electromagnetic and neutrino observatories. PyCBC Live is a low-latency search for compact binary mergers based on frequency-domain matched filtering, which has been used during the second and third observing runs, together with other low-latency analyses, to generate these rapid alerts from the data acquired by LIGO and Virgo. This paper describes and evaluates the improvements made to PyCBC Live after the second observing run, which defined its operation and performance during the third observing run., Comment: 14 pages, 6 figures. Matches version accepted by ApJ

We present a search for merging compact binary gravitational-wave sources that produce a signal appearing solely or primarily in a single detector. Past analyses have heavily relied on coincidence between multiple detectors to reduce non-astrophysical background. However, for $\sim40\%$ of the total time of the 2015-2017 LIGO-Virgo observing runs only a single detector was operating. We discuss the difficulties in assigning significance and calculating the probability of astrophysical origin for candidates observed primarily by a single detector, and suggest a straightforward resolution using a noise model designed to provide a conservative assessment given the observed data. We also describe a procedure to assess candidates observed in a single detector when multiple detectors are observing. We apply these methods to search for binary black hole (BBH) and binary neutron star (BNS) mergers in the open LIGO data spanning 2015-2017. The most promising candidate from our search is 170817+03:02:46UTC (probability of astrophysical origin $p_{\rm astro} \sim 0.4$): if astrophysical, this is consistent with a BBH merger with primary mass $67_{-15}^{+21}\,M_{\odot}$, suggestive of a hierarchical merger origin. We also apply our method to the analysis of GW190425 and find $p_{\rm astro} \sim 0.5$, though this value is highly dependent on assumptions about the noise and signal models., Comment: 11 pages, 5 figures, 2 tables. Updated to match ApJ version. Supplementary materials at https://github.com/gwastro/single-search

The worldwide advanced gravitational-wave (GW) detector network has so far primarily consisted of the two Advanced LIGO observatories at Hanford and Livingston, with Advanced Virgo joining the 2016-7 O2 observation run at a relatively late stage. However Virgo has been observing alongside the LIGO detectors since the start of the O3 run; in the near future, the KAGRA detector will join the global network and a further LIGO detector in India is under construction. Gravitational-wave search methods would therefore benefit from the ability to analyse data from an arbitrary network of detectors. In this paper, we extend the PyCBC offline compact binary coalescence (CBC) search analysis to three or more detectors, and describe resulting updates to the coincident search and event ranking statistic. For a three-detector network, our improved multi-detector search finds 23% more simulated signals at fixed false alarm rate in idealized colored Gaussian noise, and up to 40% more in real data, compared to the two-detector analysis previously used during O2., Comment: 12 pages, 6 figures

Advanced LIGO and Advanced Virgo could observe the first lensed gravitational wave sources in the coming years, while the future Einstein Telescope could observe hundreds of lensed events. It is, therefore, crucial to develop methodologies to distinguish between lensed from unlensed gravitational-wave observations. A lensed signal not identified as such will lead to biases during the interpretation of the source. In particular, sources will appear to have intrinsically higher masses. No robust method currently exists to distinguish between the magnification bias caused by lensing and intrinsically high-mass sources. In this work, we show how to recognize lensed and unlensed binary neutron star systems through the measurement of their tidal effects for highly magnified sources as a proof-of-principle. The proposed method could be used to identify lensed binary neutron stars, which are the chief candidate for lensing cosmography studies. We apply our method on GW190425, finding no evidence in favor of lensing, mainly due to the poor measurement of the event's tidal effects. However, we expect that future detections with better tidal measurements can yield better constraints., Comment: 12 pages, 7 figures

Strong gravitational lensing can produce multiple images of the same gravitational-wave signal, each arriving at different times and with different magnification. Previous work has explored if lensed pairs exist among the known high-significance events and found no evidence of this. However, the possibility remains that weaker counterparts of these events are present in the data, unrecovered by previous searches. We conduct a targeted search specifically looking for sub-threshold lensed images of known binary black hole (BBH) observations. We recover candidates matching three of the additional events first reported by Venumadhav et al. (2019), but find no evidence for additional BBH events. We also find no evidence that any of the Venumadhav et al. observations are lensed counterparts. We demonstrate how this type of counterpart search can constrain hypotheses about the overall source and lens populations and we rule out at very high confidence the extreme hypothesis that all heavy BBH detections are in fact lensed systems at high redshift with intrinsic masses < 15Msun., Comment: 16 pages. Comments welcome. Data release: https://github.com/icg-gravwaves/lensed-o1-o2-data-release. Version accepted for publication in PRD. Supplement merged into main paper following PRD style. Various updates/improvements during review process. Typos corrected

We present the second Open Gravitational-wave Catalog (2-OGC) of compact-binary coalescences, obtained from the complete set of public data from Advanced LIGO's first and second observing runs. For the first time we also search public data from the Virgo observatory. The sensitivity of our search benefits from updated methods of ranking candidate events including the effects of non-stationary detector noise and varying network sensitivity; in a separate targeted binary black hole merger search we also impose a prior distribution of binary component masses. We identify a population of 14 binary black hole merger events with probability of astrophysical origin $> 0.5$ as well as the binary neutron star merger GW170817. We confirm the previously reported events GW170121, GW170304, and GW170727 and also report GW151205, a new marginal binary black hole merger with a primary mass of $67^{+28}_{-17}\,\mathrm{M}_{\odot}$ that may have formed through hierarchical merger. We find no additional significant binary neutron star merger or neutron star--black hole merger events. To enable deeper follow-up as our understanding of the underlying populations evolves, we make available our comprehensive catalog of events, including the sub-threshold population of candidates and posterior samples from parameter inference of the 30 most significant binary black hole candidates., Comment: 14 pages, 4 Figures, 4 tables, supplemental materials at https://github.com/gwastro/2-ogc, updated to match version accepted by ApJ

In the last few years, machine learning techniques, in particular convolutional neural networks, have been investigated as a method to replace or complement traditional matched filtering techniques that are used to detect the gravitational-wave signature of merging black holes. However, to date, these methods have not yet been successfully applied to the analysis of long stretches of data recorded by the Advanced LIGO and Virgo gravitational-wave observatories. In this work, we critically examine the use of convolutional neural networks as a tool to search for merging black holes. We identify the strengths and limitations of this approach, highlight some common pitfalls in translating between machine learning and gravitational-wave astronomy, and discuss the interdisciplinary challenges. In particular, we explain in detail why convolutional neural networks alone cannot be used to claim a statistically significant gravitational-wave detection. However, we demonstrate how they can still be used to rapidly flag the times of potential signals in the data for a more detailed follow-up. Our convolutional neural network architecture as well as the proposed performance metrics are better suited for this task than a standard binary classifications scheme. A detailed evaluation of our approach on Advanced LIGO data demonstrates the potential of such systems as trigger generators. Finally, we sound a note of caution by constructing adversarial examples, which showcase interesting "failure modes" of our model, where inputs with no visible resemblance to real gravitational-wave signals are identified as such by the network with high confidence., Comment: First two authors contributed equally; appeared at Phys. Rev. D

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects., Comment: 14 pages, 2 figures, White Paper submitted to the Astro-2020 (2020 Astronomy and Astrophysics Decadal Survey) by GWIC-3G Science Case Team (GWIC: Gravitational-Wave International Committee); replaced Figure 1 with a revised version that incorporates Gravity IR Flare data point

We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario., Comment: 14 pages, 8 figures, final published version, samples available at https://git.ligo.org/katerina.chatziioannou/gw170729hm_datarelease

LIGO and Virgo recently observed the first binary neutron star merger, demonstrating that gravitational-waves offer the ability to probe how matter behaves in one of the most extreme environments in the Universe. However, the gravitational-wave signal emitted by an inspiraling binary neutron star system is only weakly dependent on the equation of state and extracting this information is challenging. Previous studies have focused mainly on binary systems where the neutron stars are spinning slowly and the main imprint of neutron star matter in the inspiral signal is due to tidal effects. For binaries with non-negligible neutron-star spin the deformation of the neutron star due to its own rotation introduces additional variations in the emitted gravitational-wave signal. Here we explore whether highly spinning binary neutron-star systems offer a better chance to measure the equation-of-state than weakly spinning binary-neutron star systems. We focus on the dominant adiabatic quadrupolar effects and consider three main questions. First, we show that equation-of-state effects can be significant in the inspiral waveforms, and that the spin-quadrupole effect dominates for rapidly rotating neutron stars. Second, we show that variations in the spin-quadrupole phasing are strongly degenerate with changes in the component masses and spins, and neglecting these terms has a negligible impact on the number of observations with second generation observatories. Finally, we explore the bias in the masses and spins that would be introduced by using incorrect equation-of-state terms. Using a novel method to rapidly evaluate an approximation of the likelihood we show that assuming the incorrect equation-of-state when measuring source parameters can lead to a significant bias. We also find that the ability to measure the equation-of-state is improved when considering spinning systems., Comment: 29 pages (CQG formatting), 9 figures

Current searches for the gravitational-wave signature of compact binary mergers rely on matched-filtering data from interferometric observatories with sets of modelled gravitational waveforms. These searches currently use model waveforms that do not include the higher-order mode content of the gravitational-wave signal. Higher-order modes are important for many compact binary mergers and their omission reduces the sensitivity to such sources. In this work we explore the sensitivity loss incurred from omitting higher-order modes. We present a new method for searching for compact binary mergers using waveforms that include higher-order mode effects, and evaluate the sensitivity increase that using our new method would allow. We find that, when evaluating sensitivity at a constant rate-of-false alarm, and when including the fact that signal-consistency tests can reject some signals that include higher-order mode content, we observe a sensitivity increase of up to a factor of 2 in volume for high mass ratio, high total-mass systems. For systems with equal mass, or with total mass $\sim 50 M_{\odot}$, we see more modest sensitivity increases, $< 10\%$, which indicates that the existing search is already performing well. Our new search method is also directly applicable in searches for generic compact binaries., Comment: 17 pages, 7 figures. Updating to match the published version (published in 2018)

Gravitational waves from merging neutron stars are expected to be observed in the next 5 years. We explore the potential impact of matter effects on gravitational waves from merging double neutron-star binaries. If neutron star binaries exist with chirp masses less than roughly 1 solar mass and typical neutron-star radii are larger than roughly 14 km, or if neutron-star radii are larger than 15-16 km for the chirp masses of galactic neutron-star binaries, then matter will have a significant impact on the effectiveness of a point-particle-based search at Advanced LIGO design sensitivity (roughly 5% additional loss of signals). In a configuration typical of LIGO's first observing run, extreme matter effects lead to up to 10% potential loss in the most extreme cases.

We describe the methodology and novel techniques used to construct a set of waveforms, or template bank, applicable to searches for compact binary coalescences in Advanced LIGO's second observing run. This template bank is suitable for observing systems composed of two neutron stars, two black holes, or a neutron star and a black hole. The Post-Newtonian formulation is used to model waveforms with total mass less than 4 $M_{\odot}$ and the most recent effective-one-body model, calibrated to numerical relativity to include the merger and ringdown, is used for total masses greater than 4 $M_{\odot}$. The effects of spin precession, matter, orbital eccentricity and radiation modes beyond the quadrupole are neglected. In contrast to the template bank used to search for compact binary mergers in Advanced LIGO's first observing run, here we are including binary-black-hole systems with total mass up to several hundreds of solar masses, thereby improving the ability to observe such systems. We introduce a technique to vary the starting frequency of waveform filters so that our bank can simultaneously contain binary-neutron-star and high-mass binary-black hole waveforms. We also introduce a lower-bound on the filter waveform length, to exclude very short-duration, high-mass templates whose sensitivity is strongly reduced by the characteristics and performance of the interferometers., Comment: 10 pages, 8 figures

This document describes the new Numerical Relativity (NR) injection infrastructure in the LIGO Algorithms Library (LAL), which henceforth allows for the usage of NR waveforms as a discrete waveform approximant in LAL. With this new interface, NR waveforms provided in the described format can directly be used as simulated GW signals ("injections") for data analyses, which include parameter estimation, searches, hardware injections etc. As opposed to the previous infrastructure, this new interface natively handles sub-dominant modes and waveforms from numerical simulations of precessing binary black holes, making them directly accessible to LIGO analyses. To correctly handle precessing simulations, the new NR injection infrastructure internally transforms the NR data into the coordinate frame convention used in LAL., Comment: 20 pages, 2 figures, technical report

We improve the accuracy of the effective-one-body (EOB) waveforms that were employed during the first observing run of Advanced LIGO for binaries of spinning, nonprecessing black holes by calibrating them to a set of 141 numerical-relativity (NR) waveforms. The NR simulations expand the domain of calibration towards larger mass ratios and spins, as compared to the previous EOBNR model. Merger-ringdown waveforms computed in black-hole perturbation theory for Kerr spins close to extremal provide additional inputs to the calibration. For the inspiral-plunge phase, we use a Markov-chain Monte Carlo algorithm to efficiently explore the calibration space. For the merger-ringdown phase, we fit the NR signals with phenomenological formulae. After extrapolation of the calibrated model to arbitrary mass ratios and spins, the (dominant-mode) EOBNR waveforms have faithfulness --- at design Advanced-LIGO sensitivity --- above $99\%$ against all the NR waveforms, including 16 additional waveforms used for validation, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling for these binary configurations. We find that future NR simulations at mass ratios $\gtrsim 4$ and double spin $\gtrsim 0.8$ will be crucial to resolve discrepancies between different ways of extrapolating waveform models. We also find that some of the NR simulations that already exist in such region of parameter space are too short to constrain the low-frequency portion of the models. Finally, we build a reduced-order version of the EOBNR model to speed up waveform generation by orders of magnitude, thus enabling intensive data-analysis applications during the upcoming observation runs of Advanced LIGO., Comment: 27 pages, 15 figures

Current searches for gravitational waves from compact-object binaries with the LIGO and Virgo observatories employ waveform models with spins aligned (or anti-aligned) with the orbital angular momentum. Here, we derive a new statistic to search for compact objects carrying generic (precessing) spins. Applying this statistic, we construct banks of both aligned- and generic-spin templates for binary black holes and neutron-star--black-hole binaries, and compare the effectualness of these banks towards simulated populations of generic-spin systems. We then use these banks in a pipeline analysis of Gaussian noise to measure the increase in background incurred by using generic- instead of aligned-spin banks. Although the generic-spin banks have a factor of ten to twenty more templates than the aligned-spin banks, we find an overall improvement in signal recovery at fixed false-alarm rate for systems with high-mass ratio and highly precessing spins ---up to 60\% for neutron-star--black-hole mergers. This gain in sensitivity comes at a small loss of sensitivity ($\lesssim$4\%) for systems that are already well-covered by aligned-spin templates. Since the observation of even a single binary merger with misalinged spins could provide unique astrophysical insights into the formation of these sources, we recommend that the method described here be developed further to mount a viable search for generic-spin binary mergers in LIGO/Virgo data., Comment: 21 pages, 9 figures. Updated to reflect published version in PRD

Searching for gravitational waves (GWs) from binary black holes (BBHs) with LIGO and Virgo involves matched-filtering data against a set of representative signal waveforms --- a template bank --- chosen to cover the full signal space of interest with as few template waveforms as possible. Although the component black holes may have significant angular momenta (spin), previous searches for BBHs have filtered LIGO and Virgo data using only waveforms where both component spins are zero. This leads to a loss of signal-to-noise ratio for signals where this is not the case. Combining the best available template placement techniques and waveform models, we construct a template bank of GW signals from BBHs with component spins $\chi_{1,2}\in [-0.99, 0.99]$ aligned with the orbital angular momentum, component masses $m_{1,2}\in [2, 48]\,\mathrm{M}_\odot$, and total mass $M_\mathrm{total} \leq 50\,\mathrm{M}_\odot$. Using effective-one-body waveforms with spin effects, we show that less than $3\%$ of the maximum signal-to-noise ratio (SNR) of these signals is lost due to the discreetness of the bank, using the early advanced LIGO noise curve. We use simulated advanced LIGO noise to compare the sensitivity of this bank to a non-spinning bank covering the same parameter space. In doing so, we consider the competing effects between improved SNR and signal-based vetoes, and the increase in the rate of false alarms of the aligned-spin bank due to covering a larger parameter space. We find that the aligned-spin bank can be a factor of $1.3$ -- $5$ more sensitive than a non-spinning bank to BBHs with dimensionless spins $> +0.6$ and component masses $\gtrsim 20\,\mathrm{M}_\odot$. Even larger gains are obtained for systems with equally high spins but smaller component masses., Comment: 17 pages, 8 figures. Version submitted to PRD

We describe the PyCBC search for gravitational waves from compact-object binary coalescences in advanced gravitational-wave detector data. The search was used in the first Advanced LIGO observing run and unambiguously identified two black hole binary mergers, GW150914 and GW151226. At its core, the PyCBC search performs a matched-filter search for binary merger signals using a bank of gravitational-wave template waveforms. We provide a complete description of the search pipeline including the steps used to mitigate the effects of noise transients in the data, identify candidate events and measure their statistical significance. The analysis is able to measure false-alarm rates as low as one per million years, required for confident detection of signals. Using data from initial LIGO's sixth science run, we show that the new analysis reduces the background noise in the search, giving a 30% increase in sensitive volume for binary neutron star systems over previous searches., Comment: 29 pages, 7 figures, accepted by Classical and Quantum Gravity

We apply a machine learning algorithm, the artificial neural network, to the search for gravitational-wave signals associated with short gamma-ray bursts. The multi-dimensional samples consisting of data corresponding to the statistical and physical quantities from the coherent search pipeline are fed into the artificial neural network to distinguish simulated gravitational-wave signals from background noise artifacts. Our result shows that the data classification efficiency at a fixed false alarm probability is improved by the artificial neural network in comparison to the conventional detection statistic. Therefore, this algorithm increases the distance at which a gravitational-wave signal could be observed in coincidence with a gamma-ray burst. In order to demonstrate the performance, we also evaluate a few seconds of gravitational-wave data segment using the trained networks and obtain the false alarm probability. We suggest that the artificial neural network can be a complementary method to the conventional detection statistic for identifying gravitational-wave signals related to the short gamma-ray bursts., Comment: 30 pages, 10 figures

We study the effect of spins on searches for gravitational waves from compact binary coalescences in realistic simulated early advanced LIGO data. We construct a detection pipeline including matched filtering, signal-based vetoes, a coincidence test between different detectors, and an estimate of the rate of background events. We restrict attention to neutron star--black hole (NS-BH) binary systems, and we compare a search using non-spinning templates to one using templates that include spins aligned with the orbital angular momentum. To run the searches we implement the binary inspiral matched-filter computation in PyCBC, a new software toolkit for gravitational-wave data analysis. We find that the inclusion of aligned-spin effects significantly increases the astrophysical reach of the search. Considering astrophysical NS-BH systems with non-precessing black hole spins, for dimensionless spin components along the orbital angular momentum uniformly distributed in $(-1, 1)$, the sensitive volume of the search with aligned-spin templates is increased by $\sim 50\%$ compared to the non-spinning search; for signals with aligned spins uniformly distributed in the range $(0.7,1)$, the increase in sensitive volume is a factor of $\sim 10$., Comment: 17 pages, 12 figures; version accepted by PRD

Gravitational waves radiated by the coalescence of compact-object binaries containing a neutron star and a black hole are one of the most interesting sources for the ground-based gravitational-wave observatories Advanced LIGO and Advanced Virgo. Advanced LIGO will be sensitive to the inspiral of a $1.4\, M_\odot$ neutron star into a $10\,M_\odot$ black hole to a maximum distance of $\sim 900$ Mpc. Achieving this sensitivity and extracting the physics imprinted in observed signals requires accurate modeling of the binary to construct template waveforms. In a NSBH binary, the black hole may have significant angular momentum (spin), which affects the phase evolution of the emitted gravitational waves. We investigate the ability of post-Newtonian (PN) templates to model the gravitational waves emitted during the inspiral phase of NSBH binaries. We restrict the black hole's spin to be aligned with the orbital angular momentum and compare several approximants. We examine restricted amplitude waveforms that are accurate to 3.5PN order in the orbital dynamics and complete to 2.5PN order in the spin dynamics. We also consider PN waveforms with the recently derived 3.5PN spin-orbit and 3PN spin-orbit tail corrections. We compare these approximants to the effective-one-body model. For all these models, large disagreements start at low to moderate black hole spins, particularly for binaries where the spin is anti-aligned with the orbital angular momentum. We show that this divergence begins in the early inspiral at $v \sim 0.2$ for $\chi_{BH} \sim 0.4$. PN spin corrections beyond those currently known will be required for optimal detection searches and to measure the parameters of neutron star--black hole binaries. While this complicates searches, the strong dependence of the gravitational-wave signal on the spin dynamics will make it possible to extract significant astrophysical information., Comment: 17 pages, 15 figures, abstract shortened to the 1920 character limit

Gravitational-wave observations of compact binaries have the potential to uncover the distribution of masses and angular momenta of black holes and neutron stars in the universe. The binary components' physical parameters can be inferred from their effect on the phasing of the gravitational-wave signal, but a partial degeneracy between the components' mass ratio and their angular momenta limits our ability to measure the individual component masses. At the typical signal amplitudes expected by the Advanced Laser Interferometer Gravitational-wave Observatory (signal-to-noise ratios between 10 and 20), we show that it will in many cases be difficult to distinguish whether the components are neutron stars or black holes. We identify when the masses of the binary components could be unambiguously measured outside the range of current observations: a system with a chirp mass $\mathcal{M} \le 0.871 $ M$_\odot$ would unambiguously contain the smallest-mass neutron star observed, and a system with $\mathcal{M} \ge 2.786 \Msun$ must contain a black hole. However, additional information would be needed to distinguish between a binary containing two 1.35 M$_\odot$ neutron stars and an exotic neutron-star--black-hole binary. We also identify those configurations that could be unambiguously identified as black-hole binaries, and show how the observation of an electromagnetic counterpart to a neutron-star--black-hole binary could be used to constrain the black-hole spin., Comment: 5 pages, 4 figures. Final version to be published in Ap.J.Lett

The detection of gravitational waves from binary neutron stars is a major goal of the gravitational-wave observatories Advanced LIGO and Advanced Virgo. Previous searches for binary neutron stars with LIGO and Virgo neglected the component stars' angular momentum (spin). We demonstrate that neglecting spin in matched-filter searches causes advanced detectors to lose more than 3% of the possible signal-to-noise ratio for 59% (6%) of sources, assuming that neutron star dimensionless spins, $c\mathbf{J}/GM^2$, are uniformly distributed with magnitudes between 0 and 0.4 (0.05) and that the neutron stars have isotropically distributed spin orientations. We present a new method for constructing template banks for gravitational wave searches for systems with spin. We present a new metric in a parameter space in which the template placement metric is globally flat. This new method can create template banks of signals with non-zero spins that are (anti-)aligned with the orbital angular momentum. We show that this search loses more than 3% of the maximium signal-to-noise for only 9% (0.2%) of BNS sources with dimensionless spins between 0 and 0.4 (0.05) and isotropic spin orientations. Use of this template bank will prevent selection bias in gravitational-wave searches and allow a more accurate exploration of the distribution of spins in binary neutron stars., Comment: 11 pages, 10 figures

In this paper we present a method for conducting a coherent search for single spin compact binary coalescences in gravitational wave data and compare this search to the existing coincidence method for single spin searches. We propose a method to characterize the regions of the parameter space where the single spin search, both coincident and coherent, will increase detection efficiency over the existing non-precessing search. We also show example results of the coherent search on a stretch of data from LIGO's fourth science run but note that a set of signal based vetoes will be needed before this search can be run to try to make detections., Comment: 14 pages, 4 figures

We present the details of a method for conducting a targeted, coherent search for compact binary coalescences. The search is tailored to be used as a followup to electromagnetic transients such as Gamma Ray Bursts. We derive the coherent search statistic for Gaussian detector noise and discuss the benefits of a coherent, multi-detector search over coincidence methods. To mitigate the effects of non-stationary data, we introduce a number of signal consistency tests, including the null SNR, amplitude consistency and several $\chi^{2}$ tests. We demonstrate the search performance on Gaussian noise and on data from LIGO's fourth science run and verify that the signal consistency tests are capable of removing the majority of noise transients and the search gives an efficiency comparable to that achieved in Gaussian noise., Comment: 21 pages, 14 figures

This paper presents an algorithm for constructing matched-filter template banks in an arbitrary parameter space. The method places templates at random, then removes those which are "too close" together. The properties and optimality of stochastic template banks generated in this manner are investigated for some simple models. The effectiveness of these template banks for gravitational wave searches for binary inspiral waveforms is also examined. The properties of a stochastic template bank are then compared to the deterministically placed template banks that are currently used in gravitational wave data analysis., Comment: 14 pages, 11 figures

Gravitational waves from coalescing compact binaries are one of the most promising sources for detectors such as LIGO, Virgo and GEO600. If the components of the binary posess significant angular momentum (spin), as is likely to be the case if one component is a black hole, spin-induced precession of a binary's orbital plane causes modulation of the gravitational-wave amplitude and phase. If the templates used in a matched-filter search do not accurately model these effects then the sensitivity, and hence the detection rate, will be reduced. We investigate the ability of several search pipelines to detect gravitational waves from compact binaries with spin. We use the post-Newtonian approximation to model the inspiral phase of the signal and construct two new template banks using the phenomenological waveforms of Buonanno, Chen and Vallisneri. We compare the performance of these template banks to that of banks constructed using the stationary phase approximation to the non-spinning post-Newtonian inspiral waveform currently used by LIGO and Virgo in the search for compact binary coalescence. We find that, at the same false alarm rate, a search pipeline using phenomenological templates is no more effective than a pipeline which uses non-spinning templates. We recommend the continued use of the non-spinning stationary phase template bank until the false alarm rate associated with templates which include spin effects can be substantially reduced., Comment: 11 pages, 3 figures

The Mock LISA Data Challenges are a programme to demonstrate and encourage the development of LISA data-analysis capabilities, tools and techniques. At the time of this workshop, three rounds of challenges had been completed, and the next was about to start. In this article we provide a critical analysis of entries to the latest completed round, Challenge 1B. The entries confirm the consolidation of a range of data-analysis techniques for Galactic and massive--black-hole binaries, and they include the first convincing examples of detection and parameter estimation of extreme--mass-ratio inspiral sources. In this article we also introduce the next round, Challenge 3. Its data sets feature more realistic waveform models (e.g., Galactic binaries may now chirp, and massive--black-hole binaries may precess due to spin interactions), as well as new source classes (bursts from cosmic strings, isotropic stochastic backgrounds) and more complicated nonsymmetric instrument noise., Comment: 20 pages, 3 EPS figures. Proceedings of the 12th Gravitational Wave Data Analysis Workshop, Cambridge MA, 13--16 December 2007. Typos corrected