Search

The LIGO-Virgo-KAGRA collaboration has announced the detection of almost 100 binary black holes so far, which have been used in several studies to infer the features of the underlying binary black hole population. From these, it is possible to predict the overall gravitational-wave (GW) fractional energy density contributed by black holes throughout the Universe, and thus estimate the gravitational-wave background (GWB) spectrum emitted in the current GW detector band. These predictions are fundamental in our forecasts for background detection and characterization, with both present and future instruments. The uncertainties in the inferred population strongly impact the predicted energy spectrum, and in this paper we present a new, flexible method to quickly calculate the energy spectrum for varying black hole population features such as the mass spectrum and redshift distribution. We implement this method in an open-access package, popstock, and extensively test its capabilities. Using popstock, we investigate how uncertainties in these distributions impact our detection capabilities and present several caveats for background estimation. In particular, we find that the standard assumption that the background signal follows a 2/3 power-law at low frequencies is both waveform and mass-model dependent, and that the signal power-law is likely shallower than previously modelled, given the current waveform and population knowledge., Comment: 14 pages, 6 figures

Using an optimal quadratic estimator, we measure the large-scale cross-correlation between maps of the stochastic gravitational-wave intensity, constructed from the first three LIGO-Virgo observing runs, and a suite of tomographic samples of galaxies covering the redshift range $z\lesssim 2$. We do not detect any statistically significant cross-correlation, but the tomographic nature of the data allows us to place constraints on the (bias-weighted) production rate density of gravitational waves by astrophysical sources as a function of cosmic time. Our constraints range from $\langle b\dot{\Omega}_{\rm GW}\rangle<3.0\times10^{-9}\,{\rm Gyr}^{-1}$ at $z\sim0.06$ to $\langle b\dot{\Omega}_{\rm GW}\rangle<2.7\times10^{-7}\,{\rm Gyr}^{-1}$ at $z\sim1.5$ (95\% C.L.), assuming a frequency spectrum of the form $f^{2/3}$ (corresponding to an astrophysical background of binary mergers), and a reference frequency $f_{\rm ref}=25\,{\rm Hz}$. Although these constraints are $\sim2$ orders of magnitude higher than the expected signal, we show that a detection may be possible with future experiments., Comment: 14 pages, 11 figures

The vast majority of gravitational-wave signals from stellar-mass compact binary mergers are too weak to be individually detected with present-day instruments and instead contribute to a faint, persistent background. This astrophysical background is targeted by searches that model the gravitational-wave ensemble collectively with a small set of parameters. The traditional search models the background as a stochastic field and estimates its amplitude by cross-correlating data from multiple interferometers. A different search uses gravitational-wave templates to marginalize over all individual event parameters and measure the duty cycle and population properties of binary mergers. Both searches ultimately estimate the total merger rate of compact binaries and are expected to yield a detection in the coming years. Given the conceptual and methodological differences between them, though, it is not well understood how their results should be mutually interpreted. In this paper, we use the Fisher information to study the implications of a background detection in terms of which region of the Universe each approach probes. Specifically, we quantify how information about the compact binary merger rate is accumulated by each search as a function of the event redshift. For the LIGO Design sensitivity and a uniform-in-comoving-volume distribution of equal-mass 30M_sol binaries, the traditional cross-correlation search obtains 99% of its information from binaries up to redshift 2.5 (average signal-to-noise-ratio <8), and the template-based search from binaries up to redshift 1.0 (average signal-to-noise-ratio ~8). While we do not calculate the total information accumulated by each search, our analysis emphasizes the need to pair any claimed detection of the stochastic background with an assessment of which binaries contribute to said detection., Comment: 16 pages, 6 figures

One of the most exciting targets of current and future gravitational-wave observations is the angular power spectrum of the astrophysical GW background. This cumulative signal encodes information about the large-scale structure of the Universe, as well as the formation and evolution of compact binaries throughout cosmic time. However, the finite rate of compact binary mergers gives rise to temporal shot noise, which introduces a significant bias in measurements of the angular power spectrum if not explicitly accounted for. Previous work showed that this bias can be removed by cross-correlating GW sky maps constructed from different observing times. However, this work considered an idealised measurement scenario, ignoring detector specifics and in particular noise contributions. Here we extend this temporal cross-correlation method to account for these difficulties, allowing us to implement the first unbiased anisotropic search pipeline for LIGO-Virgo-KAGRA data. In doing so, we show that the existing pipeline is biased even in the absence of shot noise, due to previously neglected sub-leading contributions to the noise covariance. We apply our pipeline to mock LIGO data, and find that our improved analysis will be crucial for stochastic searches from the current observing run (O4) onwards., Comment: 12 pages, 9 figures

We present a new, simulation-based inference method to compute the angular power spectrum of the distribution of foreground gravitational-wave transient events. As a first application of this method, we use the binary black hole mergers observed during the LIGO, Virgo, and KAGRA third observation run to test the spatial distribution of these sources. We find no evidence for anisotropy in their angular distribution. We discuss further applications of this method to investigate other gravitational-wave source populations and their correlations to the cosmological large-scale structure., Comment: 6 pages, 5 figures

The collection of gravitational waves (GWs) that are either too weak or too numerous to be individually resolved is commonly referred to as the gravitational-wave background (GWB). A confident detection and model-driven characterization of such a signal will provide invaluable information about the evolution of the Universe and the population of GW sources within it. We present a new, user-friendly Python--based package for gravitational-wave data analysis to search for an isotropic GWB in ground--based interferometer data. We employ cross-correlation spectra of GW detector pairs to construct an optimal estimator of the Gaussian and isotropic GWB, and Bayesian parameter estimation to constrain GWB models. The modularity and clarity of the code allow for both a shallow learning curve and flexibility in adjusting the analysis to one's own needs. We describe the individual modules which make up {\tt pygwb}, following the traditional steps of stochastic analyses carried out within the LIGO, Virgo, and KAGRA Collaboration. We then describe the built-in pipeline which combines the different modules and validate it with both mock data and real GW data from the O3 Advanced LIGO and Virgo observing run. We successfully recover all mock data injections and reproduce published results., Comment: 32 pages, 14 figures

A likely source of a gravitational-wave background (GWB) in the frequency band of the Advanced LIGO, Virgo and KAGRA detectors is the superposition of signals from the population of unresolvable stellar-mass binary-black-hole (BBH) mergers throughout the Universe. Since the duration of a BBH merger in band ($\sim\!1~{\rm s}$) is much shorter than the expected separation between neighboring mergers ($\sim\!10^3~{\rm s}$), the observed signal will be "popcorn-like" or intermittent with duty cycles of order $10^{-3}$. However, the standard cross-correlation search for stochastic GWBs currently performed by the LIGO-Virgo-KAGRA collaboration is based on a continuous-Gaussian signal model, which does not take into account the intermittent nature of the background. The latter is better described by a Gaussian mixture-model, which includes a duty cycle parameter that quantifies the degree of intermittence. Building on an earlier paper by Drasco and Flanagan, we propose a stochastic-signal-based search for intermittent GWBs. For such signals, this search performs better than the standard continuous cross-correlation search. We present results of our stochastic-signal-based approach for intermittent GWBs applied to simulated data for some simple models, and compare its performance to the other search methods, both in terms of detection and signal characterization. Additional testing on more realistic simulated data sets, e.g., consisting of astrophysically-motivated BBH merger signals injected into colored detector noise containing noise transients, will be needed before this method can be applied with confidence on real gravitational-wave data., Comment: 23 pages, 8 figures, 1 table

A stochastic gravitational-wave (GW) background consists of a large number of weak, independent and uncorrelated events of astrophysical or cosmological origin. The GW power on the sky is assumed to contain anisotropies on top of an isotropic component, i.e., the angular monopole. Complementary to the LIGO--Virgo--KAGRA (LVK) searches, we develop an efficient analysis pipeline to compute the maximum-likelihood anisotropic sky maps in stochastic backgrounds directly in the sky pixel domain using data folded over one sidereal day. We invert the full pixel-pixel correlation matrix in map-making of the GW sky, up to an optimal eigenmode cutoff decided systematically using simulations. In addition to modeled mapping, we implement a model-independent method to probe spectral shapes of stochastic backgrounds. Using data from LIGO--Virgo's first three observing runs, we obtain upper limits on anisotropies as well as the isotropic monopole as a limiting case, consistent with the LVK results. We also set constraints on the spectral shape of the stochastic background using this novel model-independent method., Comment: 16 pages, 9 figures

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe.

The collection of individually resolvable gravitational wave (GW) events makes up a tiny fraction of all GW signals which reach our detectors, while most lie below the confusion limit and go undetected. Like voices in a crowded room, the collection of unresolved signals gives rise to a background which is well-described via stochastic variables, and hence referred to as the stochastic GW background (SGWB). In this review, we provide an overview of stochastic GW signals, and characterise them based on features of interest such as generation processes and observational properties. We then review the current detection strategies for stochastic backgrounds, offering a ready-to-use manual for stochastic GW searches in real data. In the process, we distinguish between interferometric measurements of GWs, either by ground-based or space-based laser interferometers, and timing-residuals analyses with pulsar timing arrays (PTAs). These detection methods have been applied to real data both by the large GW collaborations and smaller research groups, and the most recent and instructive results are reported here. We close this review with an outlook on future observations with third generation detectors, space-based interferometers, and potential non-interferometric detection methods proposed in the literature., Comment: 71 pages, 14 figures

Gravitational wave science is a dynamical, fast-expanding research field founded on results, tools and methodologies drawn from different research areas and communities. Early career scientists entering this field must learn and combine knowledge and techniques from a range of disciplines. The Workshop on Gravitational-Wave Astrophysics for Early Career Scientists (GWAECS), held virtually in May 2021, planted the seeds of an interdisciplinary, well-connected and all-inclusive community of early career scientists working on gravitational waves, able to exchange relevant information and ideas, build a healthy professional and international environment, share and learn valuable skills, and ensure that ongoing research efforts are perpetuated and expanded in order to attain the main scientific goals envisioned by the whole community. GWAECS was the first event unifying early career scientists belonging to different communities, historically associated with different large-scale gravitational wave experiments. It provided a broad perspective on the future of gravitational waves, offered training on soft and transferable skills and allowed ample time for informal discussions between early career scientists and well-known research experts. The essence of those activities is summarised and collected in the present document, which presents a recap of each session of the workshop and aims to provide all early career scientists with a long-lasting, useful reference which constitutes the legacy of all the ideas that circulated at GWAECS., Comment: 60 pages, 12 figures; Legacy document of the Workshop on Gravitational Wave Astrophysics for Early Career Scientists (GW{\AE}CS), Lorentz Center Workshop, online, 3-7 May 2021

Detection of a stochastic background of gravitational waves is likely to occur in the next few years. Beyond searches for the isotropic component of SGWBs, there have been various mapping methods proposed to target anisotropic backgrounds. Some of these methods have been applied to data taken by the Laser Interferometer Gravitational-wave Observatories (LIGO) and Virgo. Specifically, these directional searches have focused on mapping the intensity of the signal on the sky via maximum likelihood solutions. We compare this intensity mapping approach to a previously proposed, but never employed, amplitude-phase mapping method to understand whether this latter approach may be employed in future searches. We build up our understanding of the differences between these two approaches by analysing simple toy models of time-stream data, and run mock-data mapping tests for the two methods. We find that the amplitude-phase method is only applicable to the case of a background which is phase-coherent on large scales or, at the very least, has an intrinsic coherence scale that is larger than that of the detector. Otherwise, the amplitude-phase mapping method leads to a loss of overall information, with respect to both phase and amplitude. Since we do not expect these phase-coherent properties to hold for any of the gravitational-wave background signals we hope to detect in the near future, we conclude that intensity mapping is the preferred method for such backgrounds., Comment: 15 pages, 4 figures

Given the recent advances in gravitational-wave detection technologies, the detection and characterisation of gravitational-wave backgrounds (GWBs) with the Laser Interferometer Space Antenna (LISA) is a real possibility. To assess the abilities of the LISA satellite network to reconstruct anisotropies of different angular scales and in different directions on the sky, we develop a map-maker based on an optimal quadratic estimator. The resulting maps are maximum likelihood representations of the GWB intensity on the sky integrated over a broad range of frequencies. We test the algorithm by reconstructing known input maps with different input distributions and over different frequency ranges. We find that, in an optimal scenario of well understood noise and high frequency, high SNR signals, the maximum scales LISA may probe are $\ell_{\rm max} \lesssim 15$. The map-maker also allows to test the directional dependence of LISA noise, providing insight on the directional sky sensitivity we may expect., Comment: Corresponding author: Arianna I. Renzini, as a part of their PhD thesis. 13 pages, 9 figures

We construct a model for the angular power spectrum of the instrumental noise in interferometer networks mapping gravitational wave backgrounds (GWBs) as a function of detector noise properties, network configuration and scan strategy. We use the model to calculate the noise power spectrum for current and future ground-based experiments, as well as for planned space missions. We present our results in a language similar to that used in cosmic microwave background and intensity mapping experiments, and connect the formalism with the sensitivity curves that are common lore in GWB analyses. Our formalism is implemented in a lightweight python module that we make publicly available at https://github.com/damonge/schNell., Comment: Comments welcome. Code available at https://github.com/damonge/schNell

We integrate the publicly available O1 LIGO time-domain data to obtain maximum-likelihood constraints on the Gravitational Wave Background (GWB) arising from stochastic, persistent signals. Our method produces sky-maps of the strain intensity $I$ as a function of direction on the sky at a reference frequency $f_0$. The data is integrated assuming a set of fixed power-law spectra for the signal. The maps provide upper limits on the amplitude of the GWB density $\Omega_{\rm GW}(f_0)$ and any anisotropy around the background. We find 95\% confidence upper limits of $\Omega_{\rm GW} < 4.8\times 10^{-7}$ at $f_0=50$ Hz with similar constraints on a dipole modulation for the inspiral-dominated stochastic background case., Comment: 5 pages, 3 figures, 3 tables

Einstein's theory of General Relativity implies that energy, i.e. matter, curves space-time and thus deforms lightlike geodesics, giving rise to gravitational lensing. This phenomenon is well understood in the case of the Schwarzschild metric, and has been accurately described in the past; however, lensing in the Kerr space-time has received less attention in the literature despite potential practical observational applications. In particular, lensing in such space is not expressible as the gradient of a scalar potential and as such is a source of curl-like signatures and an asymmetric shear pattern. In this paper, we develop a differentiable lensing map in the Kerr metric, reworking and extending previous approaches. By using standard tools of weak gravitational lensing, we isolate and quantify the distortion that is uniquely induced by the presence of angular momentum in the metric. We apply this framework to the distortion induced by a Kerr-like foreground object on a distribution of background of sources. We verify that the new unique lensing signature is orders of magnitude below current observational bounds for a range of lens configurations., Comment: 16 pages, 11 figures