The LIGO, Virgo, and KAGRA (LVK) gravitational-wave observatories have opened new scientific research in astrophysics, fundamental physics, and cosmology. The collaborations that build and operate these observatories release the interferometric strain data as well as a catalogue of observed signals with accompanying Bayesian posterior distributions. These posteriors, in the form of equally-weighted samples, form a dataset that is used by a multitude of further analyses seeking to constrain the population of merging black holes, identify lensed pairs of signals, and much more. However, many of these analyses rely, often implicitly, on the ability to reconstruct the likelihood and prior from the inputs to the analysis and apply resampling (a statistical technique to generate new samples varying the underlying analysis assumptions). In this work, we first provide a guide on how to reconstruct and modify the posterior density accurately from the inputs for analyses performed with the Bilby inference library. We then demonstrate and compare resampling techniques to produce new posterior sample sets and discuss Pareto-smoothing to improve the efficiency. Finally, we provide examples of how to use resampling to study observed gravitational-wave signals. We hope this guide provides a useful resource for those wishing to use open data products from the LVK for gravitational-wave astronomy.

Continuous gravitational waves (CWs) emission from neutron stars carries information about their internal structure and equation of state, and it can provide tests of General Relativity. We present a search for CWs from a set of 45 known pulsars in the first part of the fourth LIGO--Virgo--KAGRA observing run, known as O4a. We conducted a targeted search for each pulsar using three independent analysis methods considering the single-harmonic and the dual-harmonic emission models. We find no evidence of a CW signal in O4a data for both models and set upper limits on the signal amplitude and on the ellipticity, which quantifies the asymmetry in the neutron star mass distribution. For the single-harmonic emission model, 29 targets have the upper limit on the amplitude below the theoretical spin-down limit. The lowest upper limit on the amplitude is $6.4\!\times\!10^{-27}$ for the young energetic pulsar J0537-6910, while the lowest constraint on the ellipticity is $8.8\!\times\!10^{-9}$ for the bright nearby millisecond pulsar J0437-4715. Additionally, for a subset of 16 targets we performed a narrowband search that is more robust regarding the emission model, with no evidence of a signal. We also found no evidence of non-standard polarizations as predicted by the Brans-Dicke theory.

Progress in gravitational-wave astronomy depends upon having sensitive detectors with good data quality. Since the end of the LIGO-Virgo-KAGRA third Observing run in March 2020, detector-characterization efforts have lead to increased sensitivity of the detectors, swifter validation of gravitational-wave candidates and improved tools used for data-quality products. In this article, we discuss these efforts in detail and their impact on our ability to detect and study gravitational-waves. These include the multiple instrumental investigations that led to reduction in transient noise, along with the work to improve software tools used to examine the detectors data-quality. We end with a brief discussion on the role and requirements of detector characterization as the sensitivity of our detectors further improves in the future Observing runs.

Primordial gravitational waves are expected to create a stochastic background encoding information about the early Universe that may not be accessible by other means. However, the primordial background is obscured by an astrophysical foreground consisting of gravitational waves from compact binaries. We demonstrate a Bayesian method for estimating the primordial background in the presence of an astrophysical foreground. Since the background and foreground signal parameters are estimated simultaneously, there is no subtraction step, and therefore we avoid astrophysical contamination of the primordial measurement, sometimes referred to as "residuals". Additionally, since we include the non-Gaussianity of the astrophysical foreground in our model, this method represents the statistically optimal approach to the simultaneous detection of a multi-component stochastic background.

The second Gravitational-Wave Transient Catalog reported on 39 compact binary coalescences observed by the Advanced LIGO and Advanced Virgo detectors between 1 April 2019 15:00 UTC and 1 October 2019 15:00 UTC. We present GWTC-2.1, which reports on a deeper list of candidate events observed over the same period. We analyze the final version of the strain data over this period with improved calibration and better subtraction of excess noise, which has been publicly released. We employ three matched-filter search pipelines for candidate identification, and estimate the astrophysical probability for each candidate event. While GWTC-2 used a false alarm rate threshold of 2 per year, we include in GWTC-2.1, 1201 candidates that pass a false alarm rate threshold of 2 per day. We calculate the source properties of a subset of 44 high-significance candidates that have an astrophysical probability greater than 0.5. Of these candidates, 36 have been reported in GWTC-2. If the 8 additional high-significance candidates presented here are astrophysical, the mass range of events that are unambiguously identified as binary black holes (both objects $\geq 3M_\odot$) is increased compared to GWTC-2, with total masses from $\sim 14 M_\odot$for GW190924_021846 to$\sim 182 M_\odot$ for GW190426_190642. The primary components of two new candidate events (GW190403_051519 and GW190426_190642) fall in the mass gap predicted by pair instability supernova theory. We also expand the population of binaries with significantly asymmetric mass ratios reported in GWTC-2 by an additional two events (the mass ratio is less than $0.65$ and $0.44$ at $90\%$ probability for GW190403_051519 and GW190917_114630 respectively), and find that 2 of the 8 new events have effective inspiral spins \chi_\mathrm{eff} > 0 (at $90\%$ credibility), while no binary is consistent with \chi_\mathrm{eff} < 0 at the same significance.

We carry out astrophysical inference for compact binary merger events in LIGO-Virgo's first gravitational-wave transient catalog (GWTC-1) using a physically motivated calibration model. We demonstrate that importance sampling can be used to reduce the cost of what would otherwise be a computationally challenging analysis. We show that including the physical estimate for the calibration error distribution has negligible impact on the inference of parameters for the events in GWTC-1. Studying a simulated signal with matched filter signal-to-noise ratio $\text{SNR}=200$, we project that a calibration error estimate typical of GWTC-1 is likely to be negligible for the current generation of gravitational-wave detectors. We argue that other sources of systematic error---from waveforms, prior distributions, and noise modelling---are likely to be more important. Finally, using the events in GWTC-1 as standard sirens, we infer an astrophysically-informed improvement on the estimate of the calibration error in the LIGO interferometers.

The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasting 2 weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at this https URL. The main dataset, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages.

Higher-order spatial laser modes have recently been investigated as candidates for reducing test-mass thermal noise in ground-based gravitational-wave detectors such as advanced LIGO. In particular, higher-order Hermite-Gauss (HG) modes have gained attention within the community for their more robust behaviors against random test-mass surface deformations and stronger sensing and control capacities. In this letter we offer experimental investigations on various aspects of HG mode interferometry. We have generated purified HG modes up to the 12-th order $\mathrm{HG}_{6,6}$ mode, with a power conversion efficiency of 38.8% and 27.7% for the $\mathrm{HG}_{3,3}$ and $\mathrm{HG}_{6,6}$ modes respectively. We demonstrated for the first time the misalignment and mode mismatch-induced power coupling loss measurements for HG modes up to the $\mathrm{HG}_{6,6}$. We report an excellent agreement with the extended numerical power loss factors that in the ``small power loss'' region converge to $2n+1$ or $n^2+n+1$ for a misaligned or mode mismatched $\mathrm{HG}_{n,n}$ mode. We also demonstrated the wavefront sensing (WFS) signal measurement for HG modes up to the $\mathrm{HG}_{6,6}$. The measurement result is accurately in accordance with theoretical WFS gain $\beta_{n,n-1}\sqrt{n} + \beta_{n,n+1}\sqrt{n+1}$ for an $\mathrm{HG}_{n,n}$ mode, with $\beta_{n,n-1}$ being the beat coefficient of the adjacent $\mathrm{HG}_{n,n}$ and $\mathrm{HG}_{n-1,n}$ modes on a split photodetector.

Hardware injections are simulated gravitational-wave signals added to the Laser Interferometer Gravitational-wave Observatory (LIGO). The detectors' test masses are physically displaced by an actuator in order to simulate the effects of a gravitational wave. The simulated signal initiates a control-system response which mimics that of a true gravitational wave. This provides an end-to-end test of LIGO's ability to observe gravitational waves. The gravitational-wave analyses used to detect and characterize signals are exercised with hardware injections. By looking for discrepancies between the injected and recovered signals, we are able to characterize the performance of analyses and the coupling of instrumental subsystems to the detectors' output channels. This paper describes the hardware injection system and the recovery of injected signals representing binary black hole mergers, a stochastic gravitational wave background, spinning neutron stars, and sine-Gaussians.

The LIGO, Virgo, and KAGRA (LVK) gravitational-wave observatories have opened new scientific research in astrophysics, fundamental physics, and cosmology. The collaborations that build and operate these observatories release the interferometric strain data as well as a catalogue of observed signals with accompanying Bayesian posterior distributions. These posteriors, in the form of equally-weighted samples, form a dataset that is used by a multitude of further analyses seeking to constrain the population of merging black holes, identify lensed pairs of signals, and much more. However, many of these analyses rely, often implicitly, on the ability to reconstruct the likelihood and prior from the inputs to the analysis and apply resampling (a statistical technique to generate new samples varying the underlying analysis assumptions). In this work, we first provide a guide on how to reconstruct and modify the posterior density accurately from the inputs for analyses performed with the Bilby inference library. We then demonstrate and compare resampling techniques to produce new posterior sample sets and discuss Pareto-smoothing to improve the efficiency. Finally, we provide examples of how to use resampling to study observed gravitational-wave signals. We hope this guide provides a useful resource for those wishing to use open data products from the LVK for gravitational-wave astronomy.

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv \Omega_\mathrm{PBH} / \Omega_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv \Omega_\mathrm{PBH} / \Omega_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date.

The origins of merging compact binaries observed by gravitational-wave detectors remains highly uncertain. Several astrophysical channels may contribute to the overall merger rate, with distinct formation processes imprinted on the structure and correlations in the underlying distributions of binary source parameters. In the absence of confident theoretical models, the current understanding of this population mostly relies on simple parametric models that make strong assumptions and are prone to misspecification. Recent work has made progress using more flexible nonparametric models, but detailed measurement of the multidimensional population remains challenging. In pursuit of this, we present PixelPop-a high resolution Bayesian nonparametric model to infer joint distributions and parameter correlations with minimal assumptions. PixelPop densely bins the joint parameter space and directly infers the merger rate in each bin, assuming only that bins are coupled to their nearest neighbors. We demonstrate this method on mock populations with and without bivariate source correlations, employing several statistical metrics for information gain and correlation significance to quantify our nonparametric results. We show that PixelPop correctly recovers the true populations within posterior uncertainties and offers a conservative assessment of population-level features and parameter correlations. Its flexibility and tractability make it a useful data-driven tool to probe gravitational-wave populations in multiple dimensions.

Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM.

Gravitational lensing offers unique opportunities to learn about the astrophysical origin of distant sources, the abundance of intervening objects acting as lenses, and gravity and cosmology in general. However, all this information can only be retrieved as long as one can disentangle each effect from the finite number of observables. In the geometric optics regime, typical of electromagnetic radiation, when the wavelength of the lensed signal is small compared to the size of the lens, there are invariance transformations that change the mass of the lens and the source-lens configuration but leave the observables unchanged. Neglecting this ``mass-sheet degeneracy'' can lead to biased lens parameters or unrealistic low uncertainties, which could then transfer to an incorrect cosmography study. This might be different for gravitational waves as their long wavelengths can be comparable to the lens size and lensing enters into the wave-optics limit. We explore the existence of invariance transformations in the wave-optics regime of gravitational-wave lensing, extending previous work and examining the implications for astrophysical and cosmological studies. We study these invariance transformations using three different methods of increasing level of complexity: template mismatch, Fisher Matrix, and Bayesian parameter estimation. We find that, for a sufficiently loud signal, the degeneracy is partially broken and the lens and cosmological parameters, e.g. $H_0$, can be retrieved independently and unbiased. In current ground-based detectors, though, considering also population studies, a strong constraint on these parameters seems quite remote and the prevailing degeneracy implies a larger uncertainty in the lens model reconstruction. However, with better sensitivity of the third-generation ground-based detectors, a meaningful constraint on $H_0$ is possible to obtain.

Gravitational-wave observations became commonplace in Advanced LIGO-Virgo's recently concluded third observing run. 56 non-retracted candidates were identified and publicly announced in near real time. Gravitational waves from binary neutron star mergers, however, remain of special interest since they can be precursors to high-energy astrophysical phenomena like $\gamma$-ray bursts and kilonovae. While late-time electromagnetic emissions provide important information about the astrophysical processes within, the prompt emission along with gravitational waves uniquely reveals the extreme matter and gravity during - and in the seconds following - merger. Rapid communication of source location and properties from the gravitational-wave data is crucial to facilitate multi-messenger follow-up of such sources. This is especially enabled if the partner facilities are forewarned via an early-warning (pre-merger) alert. Here we describe the commissioning and performance of such a low-latency infrastructure within LIGO-Virgo. We present results from an end-to-end mock data challenge that detects binary neutron star mergers and alerts partner facilities before merger. We set expectations for these alerts in future observing runs.

Teleseismic, or distant, earthquakes regularly disrupt the operation of ground--based gravitational wave detectors such as Advanced LIGO. Here, we present \emph{EQ mode}, a new global control scheme, consisting of an automated sequence of optimized control filters that reduces and coordinates the motion of the seismic isolation platforms during earthquakes. This, in turn, suppresses the differential motion of the interferometer arms with respect to one another, resulting in a reduction of DARM signal at frequencies below 100\,mHz. Our method greatly improved the interferometers' capability to remain operational during earthquakes, with ground velocities up to 3.9\,$\mu \mbox{m/s}$ rms in the beam direction, setting a new record for both detectors. This sets a milestone in seismic controls of the Advanced LIGO detectors' ability to manage high ground motion induced by earthquakes, opening a path for further robust operation in other extreme environmental conditions.

With the discovery of the black hole binary (BBH) coalescence GW150914 the era of gravitational-wave (GW) astronomy has started. It has recently been shown that BBH with masses comparable to or higher than GW150914 would be visible in the eLISA band a few years before they finally merge in the band of ground-based detectors. This would allow for pre-merger electromagnetic alerts, dramatically increasing the chances of a joint detection, if BBH are indeed luminous in the electromagnetic band. In this paper we explore a quite different aspect of multi-band GW astronomy, and verify if, and to what extent, measurement of masses and sky position with eLISA could improve parameter estimation and tests of general relativity with ground-based detectors. We generate a catalog of 200 BBH and find that having prior information from eLISA can reduce the uncertainty in the measurement of source distance and primary black hole spin by up to factor of 2 in ground-based GW detectors. The component masses estimate from eLISA will not be refined by the ground based detectors, whereas joint analysis will yield precise characterization of the newly formed black hole and improve consistency tests of general relativity.

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results.