We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers.

Advanced LIGO and Advanced Virgo have detected gravitational waves from astronomical sources to open a new window on the Universe. To explore this new realm requires an exquisite level of detector sensitivity, meaning that the much stronger signal from instrumental and environmental noise must be rejected. Selected examples of unwanted noise in Advanced LIGO are presented. The initial focus is on how the existence of this noise (characterized by particular frequencies or time intervals) was discovered. Then, a variety of methods are used to track down the source of the noise, e.g., a fault within the instruments or coupling from an external source. The ultimate goal of this effort is to mitigate the noise by either fixing equipment or by augmenting methods to suppress the coupling to the environment.

KAGRA, the underground and cryogenic gravitational-wave detector, was operated for its solo observation from February 25 to March 10, 2020, and its first joint observation with the GEO 600 detector from April 7 to April 21, 2020 (O3GK). This study presents an overview of the input optics systems of the KAGRA detector, which consist of various optical systems, such as a laser source, its intensity and frequency stabilization systems, modulators, a Faraday isolator, mode-matching telescopes, and a high-power beam dump. These optics were successfully delivered to the KAGRA interferometer and operated stably during the observations. The laser frequency noise was observed to limit the detector sensitivity above a few kilohertz, whereas the laser intensity did not significantly limit the detector sensitivity.

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.

We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivative from $-10^{-8}$ to $10^{-9}$ Hz/s. No statistically-significant periodic gravitational-wave signal is observed by any of the four searches. As a result, upper limits on the gravitational-wave strain amplitude $h_0$ are calculated. The best upper limits are obtained in the frequency range of 100 to 200 Hz and they are ${\sim}1.1\times10^{-25}$ at 95\% confidence-level. The minimum upper limit of $1.10\times10^{-25}$ is achieved at a frequency 111.5 Hz. We also place constraints on the rates and abundances of nearby planetary- and asteroid-mass primordial black holes that could give rise to continuous gravitational-wave signals.

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO and Advanced Virgo. This is a semicoherent search that uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25 to 1600 Hz, as well as ranges in orbital speed, frequency, and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100 and 200 Hz, correspond to an amplitude h 0 of about 10 −25 when marginalized isotropically over the unknown inclination angle of the neutron star’s rotation axis, or less than 4 × 10 −26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically marginalized upper limits are close to the predicted amplitude from about 70 to 100 Hz; the limits assuming that the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40 to 200 Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500 Hz or more.

We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived from multidimensional numerical simulations and phenomenological extreme emission models. For neutrino-driven explosions the distance at which we reach 50% detection efficiency is approaching 5 kpc, and for magnetorotationally-driven explosions is up to 54 kpc. However, waveforms for extreme emission models are detectable up to 28 Mpc. For the first time, the gravitational-wave data enabled us to exclude part of the parameter spaces of two extreme emission models with confidence up to 83%, limited by coincident data coverage. Besides, using ad hoc harmonic signals windowed with Gaussian envelopes we constrained the gravitational-wave energy emitted during core-collapse at the levels of $4.27\times 10^{-4}\,M_\odot c^2$ and $1.28\times 10^{-1}\,M_\odot c^2$ for emissions at 235 Hz and 1304 Hz respectively. These constraints are two orders of magnitude more stringent than previously derived in the corresponding analysis using initial LIGO, initial Virgo and GEO 600 data.

The gravitational-wave signal GW190521 is consistent with a binary black hole merger source at redshift 0.8 with unusually high component masses, $85^{+21}_{-14}\,M_{\odot}$ and $66^{+17}_{-18}\,M_{\odot}$, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range $65 - 120\,M_{\odot}$. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger $(142^{+28}_{-16}\,M_{\odot})$ classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular binary black hole coalescence, we detail the physical properties of GW190521's source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be $0.13^{+0.30}_{-0.11}\,{\rm Gpc}^{-3}\,\rm{yr}^{-1}$. We discuss the astrophysical implications of GW190521 for stellar collapse, and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescence, or via hierarchical merger of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

Two analysis errors have been identified that affect the results for a handful of the high-value pulsars given in Table 1 of Abbott et al. (2019). One affects the Bayesian analysis for the five pulsars that glitched during the analysis period, and the other affects the 5n-vector analysis for J0711-6830. Updated results after correcting the errors are shown in Table 1, which now supersedes the results given for those pulsars in Table 1 of Abbott et al. (2019). Updated versions of figures can be seen in Figures 1-4. Bayesian analysis.-For the glitching pulsars, the signal phase evolution caused by the glitch was wrongly applied twice and was therefore not consistent with our expected model of the pulsar phase. This error did not affect the F/G-statistic or 5n-vector analysis. Analyses of the five pulsars PSR J0205+6449, PSR J0534+2200, PSR J0835-4510, PSR J1028-5819, and PSR J1718-3825 have been repeated after correcting for the error. There are small quantitative differences in the results, but the changes do not affect the main conclusions of the paper. The largest differences are for PSR J0835-4510 (the Vela pulsar), for which the updated upper limits from the Bayesian method are found to be between 1.1 and 2 times larger than those obtained when the error was present. This appears primarily to be due to the error leading to the decohering of a strong spectral line in the LIGO Livingston detector and thus lowering the amplitude limit. 5n-vector analysis.-An error was also identified in the settings of the 5n-vector analysis, which affected the upper limit computation at the rotation frequency for C21 95% of J0711-6830. Specifically, we found an incorrect choice for the range of amplitudes used to inject simulated signals in the O2 data. The updated upper limit is about 2.5 times worse than that obtained when the error was present. This error did not affect the Bayesian or F/G-statistic results. (Table Presented) (Figure Presented).

Isolated spinning neutron stars, asymmetric with respect to their rotation axis, are expected to be sources of continuous gravitational waves. The most sensitive searches for these sources are based on accurate matched filtering techniques, that assume the continuous wave to be phase-locked with the pulsar beamed emission. While matched filtering maximizes the search sensitivity, a significant signal-to-noise ratio loss will happen in case of a mismatch between the assumed and the true signal phase evolution. Narrow-band algorithms allow for a small mismatch in the frequency and spin-down values of the pulsar while integrating coherently the entire data set. In this paper we describe a narrow-band search using LIGO O2 data for the continuous wave emission of 33 pulsars. No evidence for a continuous wave signal has been found and upper-limits on the gravitational wave amplitude, over the analyzed frequency and spin-down volume, have been computed for each of the targets. In this search we have surpassed the spin-down limit for some of the pulsars already present in the O1 LIGO narrow-band search, such as J1400\textminus6325 J1813\textminus1246, J1833\textminus1034, J1952+3252, and for new targets such as J0940\textminus5428 and J1747\textminus2809. For J1400\textminus6325, J1833\textminus1034 and J1747\textminus2809 this is the first time the spin-down limit is surpassed.

We present a search for gravitational waves from 222 pulsars with rotation frequencies greater than or similar to 10 Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the l = m = 2 mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and the l = 2, m = 1 mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the l = m = 2 mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of 10(38) kg m(2) for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of 1.2 x 10(-8). We also place new limits on the emission amplitude at the rotation frequency of the pulsars.

We analyze the impact of a proposed tidal instability coupling p modes and g modes within neutron stars on GW170817. This nonresonant instability transfers energy from the orbit of the binary to internal modes of the stars, accelerating the gravitational-wave driven inspiral. We model the impact of this instability on the phasing of the gravitational wave signal using three parameters per star: an overall amplitude, a saturation frequency, and a spectral index. Incorporating these additional parameters, we compute the Bayes factor (lnB_{!pg}^{pg}) comparing our p-g model to a standard one. We find that the observed signal is consistent with waveform models that neglect p-g effects, with lnB_{!pg}^{pg}=0.03_{-0.58}^{+0.70} (maximum a posteriori and 90% credible region). By injecting simulated signals that do not include p-g effects and recovering them with the p-g model, we show that there is a ≃50% probability of obtaining similar lnB_{!pg}^{pg} even when p-g effects are absent. We find that the p-g amplitude for 1.4  M_{⊙} neutron stars is constrained to less than a few tenths of the theoretical maximum, with maxima a posteriori near one-tenth this maximum and p-g saturation frequency ∼70  Hz. This suggests that there are less than a few hundred excited modes, assuming they all saturate by wave breaking. For comparison, theoretical upper bounds suggest ≲10^{3} modes saturate by wave breaking. Thus, the measured constraints only rule out extreme values of the p-g parameters. They also imply that the instability dissipates ≲10^{51}  erg over the entire inspiral, i.e., less than a few percent of the energy radiated as gravitational waves.

Advanced LIGO's second observing run (O2), conducted from November 30, 2016 to August 25, 2017, combined with Advanced Virgo's first observations in August 2017 witnessed the birth of gravitational-wave multi-messenger astronomy. The first ever gravitational-wave detection from the coalescence of two neutron stars, GW170817, and its gamma-ray counterpart, GRB 170817A, led to an electromagnetic follow-up of the event at an unprecedented scale. Several teams from across the world searched for EM/neutrino counterparts to GW170817, paving the way for the discovery of optical, X-ray, and radio counterparts. In this article, we describe the online identification of gravitational-wave transients and the distribution of gravitational-wave alerts by the LIGO and Virgo collaborations during O2. We also describe the gravitational-wave observables which were sent in the alerts to enable searches for their counterparts. Finally, we give an overview of the online candidate alerts shared with observing partners during O2. Alerts were issued for 14 candidates, six of which have been confirmed as gravitational-wave events associated with the merger of black holes or neutron stars. Eight of the 14 alerts were issued less than an hour after data acquisition.

We summarize the four tests of general relativity carried out using the ten binary black hole signals detected by Advanced LIGO and Advanced Virgo and included in their first catalog, GWTC-1. These events provide unprecedented opportunities for testing the predictions of general relativity for the gravitational waveforms from these highly dynamical, strong-field events. The first two tests check the consistency of the residuals with noise and the consistency of the low- and high-frequency parts of the signal. The other two tests check that parameterized deviations in the waveform model---including in the post-Newtonian coefficients---are consistent with zero, and that the propagation of the waves is nondispersive. These tests reveal no evidence for deviations from general relativity, and the combined constraints improve previous results by factors of up to 2.4. We also check that the binary black hole signals observed by all three detectors do not give stronger constraints on alternative polarizations than those obtained from GW170817.

We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above 1 M-circle dot during the first and second observing runs of the advanced gravitational-wave detector network. During the first observing run (O1), from September 12, 2015 to January 19, 2016, gravitational waves from three binary black hole mergers were detected. The second observing run (O2), which ran from November 30, 2016 to August 25, 2017, saw the first detection of gravitational waves from a binary neutron star inspiral, in addition to the observation of gravitational waves from a total of seven binary black hole mergers, four of which we report here for the first time: GW170729, GW170809, GW170818, and GW170823. For all significant gravitational-wave events, we provide estimates of the source properties. The detected binary black holes have total masses between 18.6(-0.7)(+3.2) M-circle dot and 84.4(-11.1)(+15.8) M-circle dot and range in distance between 320(-110)(+120) and 2840(-1360)(+1400) Mpc. No neutron star-black hole mergers were detected. In addition to highly significant gravitational-wave events, we also provide a list of marginal event candidates with an estimated false-alarm rate less than 1 per 30 days. From these results over the first two observing runs, which include approximately one gravitational-wave detection per 15 days of data searched, we infer merger rates at the 90% confidence intervals of 110 - 3840 Gpc(-3) y(-1) for binary neutron stars and 9.7 - 101 Gpc(-3) y(-1) for binary black holes assuming fixed population distributions and determine a neutron star-black hole merger rate 90% upper limit of 610 Gpc(-3) y(-1).

On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just ∼40 Mpc, consistent with the gravitational-wave measurement, and the merger was localized to be at a projected distance of ∼2 kpc away from the galaxy’s center. We use this minimal set of facts and the mass posteriors of the two neutron stars to derive the first constraints on the progenitor of GW170817 at the time of the second supernova (SN). We generate simulated progenitor populations and follow the three-dimensional kinematic evolution from binary neutron star (BNS) birth to the merger time, accounting for pre-SN galactic motion, for considerably different input distributions of the progenitor mass, pre-SN semimajor axis, and SN-kick velocity. Though not considerably tight, we find these constraints to be comparable to those for Galactic BNS progenitors. The derived constraints are very strongly influenced by the requirement of keeping the binary bound after the second SN and having the merger occur relatively close to the center of the galaxy. These constraints are insensitive to the galaxy’s star formation history, provided the stellar populations are older than 1 Gyr.

We present a multi-messenger measurement of the Hubble constant H0 using the binary–black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO)/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DES Year 3 data. Black hole mergers such as GW170814 are expected to lack bright electromagnetic emission to uniquely identify their host galaxies and build an object-by-object Hubble diagram. However, they are suitable for a statistical measurement, provided that a galaxy catalog of adequate depth and redshift completion is available. Here we present the first Hubble parameter measurement using a black hole merger. Our analysis results in H 75 km s Mpc 0 32 40 1 1 = + , which is consistent with both SN Ia and cosmic microwave background measurements of the Hubble constant. The quoted 68% credible region comprises 60% of the uniform prior range [20, 140] km s Mpc, and it depends on the assumed prior range. If we take a broader prior of [10, 220] kmsMpc, we find H 78 km s Mpc 0 24 96 1 1 = + (57% of the prior range). Although a weak constraint on the Hubble constant from a single event is expected using the dark siren method, a multifold increase in the LVC event rate is anticipated in the coming years and combinations of many sirens will lead to improved constraints on H0.