Search

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10–500 s in a frequency band of 40–1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between $3.4\times 10^{-5}$ and $9.4\times 10^{-4}{\mathrm{Mpc}}^{-3}{\mathrm{yr}}^{-1}$ at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves.

We present the results of a search for long-duration gravitational-wave transients in the data from the Advanced LIGO second observation run; we search for gravitational-wave transients of 2–500 s duration in the 24–2048 Hz frequency band with minimal assumptions about signal properties such as waveform morphologies, polarization, sky location or time of occurrence. Signal families covered by these search algorithms include fallback accretion onto neutron stars, broadband chirps from innermost stable circular orbit waves around rotating black holes, eccentric inspiral-merger-ringdown compact binary coalescence waveforms, and other models. The second observation run totals about 118.3 days of coincident data between November 2016 and August 2017. We find no significant events within the parameter space that we searched, apart from the already-reported binary neutron star merger GW170817. We thus report sensitivity limits on the root-sum-square strain amplitude $h_{\mathrm{rss}}$ at 50% efficiency. These sensitivity estimates are an improvement relative to the first observing run and also done with an enlarged set of gravitational-wave transient waveforms. Overall, the best search sensitivity is $h_{\mathrm{rss}}^{50\%}=2.7\times {10}^{-22}{\mathrm{Hz}}^{-1/2}$ for a millisecond magnetar model. For eccentric compact binary coalescence signals, the search sensitivity reaches $h_{\mathrm{rss}}^{50\%}=9.6\times {10}^{-22}{\mathrm{Hz}}^{-1/2}$.

The detection of gravitational waves by Advanced LIGO and Advanced Virgo provides an opportunity to test general relativity in a regime that is inaccessible to traditional astronomical observations and laboratory tests. We present four tests of the consistency of the data with binary black hole gravitational waveforms predicted by general relativity. One test subtracts the best-fit waveform from the data and checks the consistency of the residual with detector noise. The second test checks the consistency of the low- and high-frequency parts of the observed signals. The third test checks that phenomenological deviations introduced in the waveform model (including in the post-Newtonian coefficients) are consistent with 0. The fourth test constrains modifications to the propagation of gravitational waves due to a modified dispersion relation, including that from a massive graviton. We present results both for individual events and also results obtained by combining together particularly strong events from the first and second observing runs of Advanced LIGO and Advanced Virgo, as collected in the catalog GWTC-1. We do not find any inconsistency of the data with the predictions of general relativity and improve our previously presented combined constraints by factors of 1.1 to 2.5. In particular, we bound the mass of the graviton to be $m_g\le 4.7\times {10}^{-23}\mathrm{eV}/c^2$ (90% credible level), an improvement of a factor of 1.6 over our previously presented results. Additionally, we check that the four gravitational-wave events published for the first time in GWTC-1 do not lead to stronger constraints on alternative polarizations than those published previously.

Searches are under way in Advanced LIGO and Virgo data for persistent gravitational waves from continuous sources, e.g. rapidly rotating galactic neutron stars, and stochastic sources, e.g. relic gravitational waves from the Big Bang or superposition of distant astrophysical events such as mergers of black holes or neutron stars. These searches can be degraded by the presence of narrow spectral artifacts (lines) due to instrumental or environmental disturbances. We describe a variety of methods used for finding, identifying and mitigating these artifacts, illustrated with particular examples. Results are provided in the form of lists of line artifacts that can safely be treated as non-astrophysical. Such lists are used to improve the efficiencies and sensitivities of continuous and stochastic gravitational wave searches by allowing vetoes of false outliers and permitting data cleaning.

Cosmic strings are topological defects which can be formed in grand unified theory scale phase transitions in the early universe. They are also predicted to form in the context of string theory. The main mechanism for a network of Nambu-Goto cosmic strings to lose energy is through the production of loops and the subsequent emission of gravitational waves, thus offering an experimental signature for the existence of cosmic strings. Here we report on the analysis conducted to specifically search for gravitational-wave bursts from cosmic string loops in the data of Advanced LIGO 2015-2016 observing run (O1). No evidence of such signals was found in the data, and as a result we set upper limits on the cosmic string parameters for three recent loop distribution models. In this paper, we initially derive constraints on the string tension $G\mu$ and the intercommutation probability, using not only the burst analysis performed on the O1 data set but also results from the previously published LIGO stochastic O1 analysis, pulsar timing arrays, cosmic microwave background and big-bang nucleosynthesis experiments. We show that these data sets are complementary in that they probe gravitational waves produced by cosmic string loops during very different epochs. Finally, we show that the data sets exclude large parts of the parameter space of the three loop distribution models we consider.

In Advanced LIGO, detection and astrophysical source parameter estimation of the binary black hole merger GW150914 requires a calibrated estimate of the gravitational-wave strain sensed by the detectors. Producing an estimate from each detector’s differential arm length control loop readout signals requires applying time domain filters, which are designed from a frequency domain model of the detector’s gravitational-wave response. The gravitational-wave response model is determined by the detector’s opto-mechanical response and the properties of its feedback control system. The measurements used to validate the model and characterize its uncertainty are derived primarily from a dedicated photon radiation pressure actuator, with cross-checks provided by optical and radio frequency references. We describe how the gravitational-wave readout signal is calibrated into equivalent gravitational-wave-induced strain and how the statistical uncertainties and systematic errors are assessed. Detector data collected over 38 calendar days, from September 12 to October 20, 2015, contain the event GW150914 and approximately 16 days of coincident data used to estimate the event false alarm probability. The calibration uncertainty is less than 10% in magnitude and 10° in phase across the relevant frequency band, 20 Hz to 1 kHz.

We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, including several performed specifically to reproduce this event. Our calculations go beyond existing semianalytic models, because for all simulations—including sources with two independent, precessing spins—we perform comparisons which account for all the spin-weighted quadrupolar modes, and separately which account for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported by Abbott et al. [Phys. Rev. Lett. 116, 241102 (2016)] (at the 90% credible level), we find the data are compatible with a wide range of nonprecessing and precessing simulations. Follow-up simulations performed using previously estimated binary parameters most resemble the data, even when all quadrupolar and octopolar modes are included. Comparisons including only the quadrupolar modes constrain the total redshifted mass $M_z\in [64M_{\odot }-82M_{\odot }]$, mass ratio $1/q=m_2/m_1\in [0.6,1]$, and effective aligned spin ${\chi }_{\mathrm{eff}}\in [-0.3,0.2]$, where ${\chi }_{\mathrm{eff}}=({\mathbf{S}}_1/m_1+{\mathbf{S}}_2/m_2)·\widehat{\mathbf{L}}/M$. Including both quadrupolar and octopolar modes, we find the mass ratio is even more tightly constrained. Even accounting for precession, simulations with extreme mass ratios and effective spins are highly inconsistent with the data, at any mass. Several nonprecessing and precessing simulations with similar mass ratio and ${\chi }_{\mathrm{eff}}$ are consistent with the data. Though correlated, the components’ spins (both in magnitude and directions) are not significantly constrained by the data: the data is consistent with simulations with component spin magnitudes $a_{1,2}$ up to at least 0.8, with random orientations. Further detailed follow-up calculations are needed to determine if the data contain a weak imprint from transverse (precessing) spins. For nonprecessing binaries, interpolating between simulations, we reconstruct a posterior distribution consistent with previous results. The final black hole’s redshifted mass is consistent with $M_{f,z}$ in the range $64.0M_{\odot }-73.5M_{\odot }$ and the final black hole’s dimensionless spin parameter is consistent with $a_f=0.62–0.73$. As our approach invokes no intermediate approximations to general relativity and can strongly reject binaries whose radiation is inconsistent with the data, our analysis provides a valuable complement to Abbott et al. [Phys. Rev. Lett. 116, 241102 (2016)].

On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of $1.0\times {10}^{-21}$. It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203 000 years, equivalent to a significance greater than $5.1\sigma$. The source lies at a luminosity distance of $410_{-180}^{+160}\mathrm{Mpc}$ corresponding to a redshift $z=0.09_{-0.04}^{+0.03}$. In the source frame, the initial black hole masses are $36_{-4}^{+5}{M}_{\odot }$ and $29_{-4}^{+4}{M}_{\odot }$, and the final black hole mass is $62_{-4}^{+4}{M}_{\odot }$, with $3.0_{-0.5}^{+0.5}{M}_{\odot }{c}^2$ radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

The LIGO detection of GW150914 provides an unprecedented opportunity to study the two-body motion of a compact-object binary in the large-velocity, highly nonlinear regime, and to witness the final merger of the binary and the excitation of uniquely relativistic modes of the gravitational field. We carry out several investigations to determine whether GW150914 is consistent with a binary black-hole merger in general relativity. We find that the final remnant’s mass and spin, as determined from the low-frequency (inspiral) and high-frequency (postinspiral) phases of the signal, are mutually consistent with the binary black-hole solution in general relativity. Furthermore, the data following the peak of GW150914 are consistent with the least-damped quasinormal mode inferred from the mass and spin of the remnant black hole. By using waveform models that allow for parametrized general-relativity violations during the inspiral and merger phases, we perform quantitative tests on the gravitational-wave phase in the dynamical regime and we determine the first empirical bounds on several high-order post-Newtonian coefficients. We constrain the graviton Compton wavelength, assuming that gravitons are dispersed in vacuum in the same way as particles with mass, obtaining a 90%-confidence lower bound of ${10}^{13}\mathrm{km}$. In conclusion, within our statistical uncertainties, we find no evidence for violations of general relativity in the genuinely strong-field regime of gravity.

We report on a comprehensive all-sky search for periodic gravitational waves in the frequency band 100–1500 Hz and with a frequency time derivative in the range of $[-1.18,+1.00]\times 10^{-8}\mathrm{Hz}/\mathrm{s}$. Such a signal could be produced by a nearby spinning and slightly nonaxisymmetric isolated neutron star in our galaxy. This search uses the data from the initial LIGO sixth science run and covers a larger parameter space with respect to any past search. A Loosely Coherent detection pipeline was applied to follow up weak outliers in both Gaussian (95% recovery rate) and non-Gaussian (75% recovery rate) bands. No gravitational wave signals were observed, and upper limits were placed on their strength. Our smallest upper limit on worst-case (linearly polarized) strain amplitude $h_0$ is $9.7\times 10^{-25}$ near 169 Hz, while at the high end of our frequency range we achieve a worst-case upper limit of $5.5\times 10^{-24}$. Both cases refer to all sky locations and entire range of frequency derivative values.

The Advanced LIGO observatories detected gravitational waves from two binary black hole mergers during their first observation run (O1). We present a high-energy neutrino follow-up search for the second gravitational wave event, GW151226, as well as for gravitational wave candidate LVT151012. We find two and four neutrino candidates detected by IceCube, and one and zero detected by Antares, within $\pm 500\mathrm{s}$ around the respective gravitational wave signals, consistent with the expected background rate. None of these neutrino candidates are found to be directionally coincident with GW151226 or LVT151012. We use nondetection to constrain isotropic-equivalent high-energy neutrino emission from GW151226, adopting the GW event’s 3D localization, to less than $2\times 10^{51}–2\times 10^{54}\mathrm{erg}$.

We report on an all-sky search for periodic gravitational waves in the frequency band 20–475 Hz and with a frequency time derivative in the range of $[-1.0,+0.1]\times {10}^{-8}\mathrm{Hz}/\mathrm{s}$. Such a signal could be produced by a nearby spinning and slightly nonaxisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO’s first observational run, O1. No periodic gravitational wave signals were observed, and upper limits were placed on their strengths. The lowest upper limits on worst-case (linearly polarized) strain amplitude $h_0$ are $\sim 4\times {10}^{-25}$ near 170 Hz. For a circularly polarized source (most favorable orientation), the smallest upper limits obtained are $\sim 1.5\times {10}^{-25}$. These upper limits refer to all sky locations and the entire range of frequency derivative values. For a population-averaged ensemble of sky locations and stellar orientations, the lowest upper limits obtained for the strain amplitude are $\sim 2.5\times {10}^{-25}$.

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and $1700{\mathrm{Gpc}}^{-3}{\mathrm{yr}}^{-1}$ and the neutron star–black hole merger rate to be between 7.8 and $140{\mathrm{Gpc}}^{-3}{\mathrm{yr}}^{-1}$, assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and $44{\mathrm{Gpc}}^{-3}{\mathrm{yr}}^{-1}$ at a fiducial redshift ($z=0.2$). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to $(1+z{)}^{\kappa }$ with $\kappa =2.9_{-1.8}^{+1.7}$ for $z\lesssim 1$. Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from ${1.2}_{-0.2}^{+0.1}$ to ${2.0}_{-0.3}^{+0.3}{M}_{\odot }$. We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of ${8.3}_{-0.5}^{+0.3}$ and ${27.9}_{-1.8}^{+1.9}{M}_{\odot }$. While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately $60M_{\odot }$, which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below ${\chi }_i\approx 0.25$. While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum.

The ever-increasing number of detections of gravitational waves from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of those detectors. We restrict our analysis to the 15 confident signals that have false alarm rates $\le {10}^{-3}{\mathrm{yr}}^{-1}$. In addition to signals consistent with binary black hole mergers, the new events include GW200115_042309, a signal consistent with a neutron star–black hole merger. We find the residual power, after subtracting the best fit waveform from the data for each event, to be consistent with the detector noise. Additionally, we find all the post-Newtonian deformation coefficients to be consistent with the predictions from GR, with an improvement by a factor of $\sim 2$ in the $-1\mathrm{PN}$ parameter. We also find that the spin-induced quadrupole moments of the binary black hole constituents are consistent with those of Kerr black holes in GR. We find no evidence for dispersion of gravitational waves, non-GR modes of polarization, or post-merger echoes in the events that were analyzed. We update the bound on the mass of the graviton, at 90% credibility, to $m_g\le 2.42\times {10}^{-23}\mathrm{eV}/c^2$. The final mass and final spin as inferred from the premerger and postmerger parts of the waveform are consistent with each other. The studies of the properties of the remnant black holes, including deviations of the quasinormal mode frequencies and damping times, show consistency with the predictions of GR. In addition to considering signals individually, we also combine results from the catalog of gravitational waves signals to calculate more precise population constraints. We find no evidence in support of physics beyond general relativity.

Quantum fluctuations in the phase and amplitude quadratures of light set limitations on the sensitivity of modern optical instruments. The sensitivity of the interferometric gravitational-wave detectors, such as the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO), is limited by quantum shot noise, quantum radiation pressure noise, and a set of classical noises. We show how the quantum properties of light can be used to distinguish these noises using correlation techniques. Particularly, in the first part of the paper we show estimations of the coating thermal noise and gas phase noise, hidden below the quantum shot noise in the Advanced LIGO sensitivity curve. We also make projections on the observatory sensitivity during the next science runs. In the second part of the paper we discuss the correlation technique that reveals the quantum radiation pressure noise from the background of classical noises and shot noise. We apply this technique to the Advanced LIGO data, collected during the first science run, and experimentally estimate the quantum correlations and quantum radiation pressure noise in the interferometer.

In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz–1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO 600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyze GRB events with large sky localization uncertainties, for example the localizations provided by the Fermi Gamma-ray Burst Monitor (GBM). Coherent searches for gravitational waves (GWs) can be computationally intensive when the GRB sky position is not well localized, due to the corrections required for the difference in arrival time between detectors. Using a linear search grid we are able to reduce the computational cost of the analysis by a factor of $\mathcal{O}(10)$ for GBM events. Furthermore, we demonstrate that our analysis pipeline can improve upon the sky localization of GRBs detected by the GBM, if a high-frequency GW signal is observed in coincidence. We use the method of the linear grid in a search for GWs associated with 129 GRBs observed satellite-based gamma-ray experiments between 2006 and 2011. The GRBs in our sample had not been previously analyzed for GW counterparts. A fraction of our GRB events are analyzed using data from GEO 600 while the detector was using squeezed-light states to improve its sensitivity; this is the first search for GWs using data from a squeezed-light interferometric observatory. We find no evidence for GW signals, either with any individual GRB in this sample or with the population as a whole. For each GRB we place lower bounds on the distance to the progenitor, under an assumption of a fixed GW emission energy of ${10}^{-2}\mathrm{M}\odot c^2$, with a median exclusion distance of 0.8 Mpc for emission at 500 Hz and 0.3 Mpc at 1 kHz. The reduced computational cost associated with a linear search grid will enable rapid searches for GWs associated with Fermi GBM events once the advanced LIGO and Virgo detectors begin operation.

Gravitational waves from a variety of sources are predicted to superpose to create a stochastic background. This background is expected to contain unique information from throughout the history of the Universe that is unavailable through standard electromagnetic observations, making its study of fundamental importance to understanding the evolution of the Universe. We carry out a search for the stochastic background with the latest data from the LIGO and Virgo detectors. Consistent with predictions from most stochastic gravitational-wave background models, the data display no evidence of a stochastic gravitational-wave signal. Assuming a gravitational-wave spectrum of ${\mathrm{\Omega }}_{\mathrm{GW}}(f)={\mathrm{\Omega }}_{\alpha }{(f/f_{\mathrm{ref}})}^{\alpha }$, we place 95% confidence level upper limits on the energy density of the background in each of four frequency bands spanning 41.5–1726 Hz. In the frequency band of 41.5–169.25 Hz for a spectral index of $\alpha =0$, we constrain the energy density of the stochastic background to be ${\mathrm{\Omega }}_{\mathrm{GW}}(f)<5.6\times 10^{-6}$. For the 600–1000 Hz band, ${\mathrm{\Omega }}_{\mathrm{GW}}(f)<0.14(f/900\mathrm{Hz}{)}^3$, a factor of 2.5 lower than the best previously reported upper limits. We find ${\mathrm{\Omega }}_{\mathrm{GW}}(f)<1.8\times 10^{-4}$ using a spectral index of zero for 170–600 Hz and ${\mathrm{\Omega }}_{\mathrm{GW}}(f)<1.0(f/1300\mathrm{Hz}{)}^3$ for 1000–1726 Hz, bands in which no previous direct limits have been placed. The limits in these four bands are the lowest direct measurements to date on the stochastic background. We discuss the implications of these results in light of the recent claim by the BICEP2 experiment of the possible evidence for inflationary gravitational waves.

We present an archival search for transient gravitational-wave bursts in coincidence with 27 single-pulse triggers from Green Bank Telescope pulsar surveys, using the LIGO, Virgo, and GEO interferometer network. We also discuss a check for gravitational-wave signals in coincidence with Parkes fast radio bursts using similar methods. Data analyzed in these searches were collected between 2007 and 2013. Possible sources of emission of both short-duration radio signals and transient gravitational-wave emission include starquakes on neutron stars, binary coalescence of neutron stars, and cosmic string cusps. While no evidence for gravitational-wave emission in coincidence with these radio transients was found, the current analysis serves as a prototype for similar future searches using more sensitive second-generation interferometers.

We present results of a search for continuously emitted gravitational radiation, directed at the brightest low-mass x-ray binary, Scorpius X-1. Our semicoherent analysis covers 10 days of LIGO S5 data ranging from 50–550 Hz, and performs an incoherent sum of coherent $\mathcal{F}$-statistic power distributed amongst frequency-modulated orbital sidebands. All candidates not removed at the veto stage were found to be consistent with noise at a 1% false alarm rate. We present Bayesian 95% confidence upper limits on gravitational-wave strain amplitude using two different prior distributions: a standard one, with no a priori assumptions about the orientation of Scorpius X-1; and an angle-restricted one, using a prior derived from electromagnetic observations. Median strain upper limits of $1.3\times {10}^{-24}$ and $8\times {10}^{-25}$ are reported at 150 Hz for the standard and angle-restricted searches respectively. This proof-of-principle analysis was limited to a short observation time by unknown effects of accretion on the intrinsic spin frequency of the neutron star, but improves upon previous upper limits by factors of $\sim 1.4$ for the standard, and 2.3 for the angle-restricted search at the sensitive region of the detector.

We present the first results of an all-sky search for continuous gravitational waves from unknown spinning neutron stars in binary systems using LIGO and Virgo data. Using a specially developed analysis program, the TwoSpect algorithm, the search was carried out on data from the sixth LIGO science run and the second and third Virgo science runs. The search covers a range of frequencies from 20 Hz to 520 Hz, a range of orbital periods from 2 to $\sim 2,254\mathrm{h}$ and a frequency- and period-dependent range of frequency modulation depths from 0.277 to 100 mHz. This corresponds to a range of projected semimajor axes of the orbit from $\sim 0.6\times 10^{-3}\mathrm{ls}$ to $\sim 6,500\mathrm{ls}$ assuming the orbit of the binary is circular. While no plausible candidate gravitational wave events survive the pipeline, upper limits are set on the analyzed data. The most sensitive 95% confidence upper limit obtained on gravitational wave strain is $2.3\times 10^{-24}$ at 217 Hz, assuming the source waves are circularly polarized. Although this search has been optimized for circular binary orbits, the upper limits obtained remain valid for orbital eccentricities as large as 0.9. In addition, upper limits are placed on continuous gravitational wave emission from the low-mass x-ray binary Scorpius X-1 between 20 Hz and 57.25 Hz.