Search for Subsolar-Mass Binaries in the First Half of Advanced LIGO's and Advanced Virgo's Third Observing Run

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2  M_{⊙} and 1.0  M_{⊙} in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend our 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≥0.1. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14  yr^{-1}. This implies an upper limit on the merger rate of subsolar binaries in the range [220-24200]  Gpc^{-3} 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 in the mass range 0.2  M_{⊙}<m_{PBH}<1.0  M_{⊙} is f_{PBH}≡Ω_{PBH}/Ω_{DM}≲6%. This improves existing constraints on primordial black hole abundance by a factor of ∼3. 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_{min}=1  M_{⊙}, where f_{DBH}≡Ω_{DBH}/Ω_{DM}≲0.003%. These are the first constraints placed on dissipative dark models by subsolar-mass analyses.

GW190521: A Binary Black Hole Merger with a Total Mass of 150 M_{⊙}

On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of two black holes with masses of 85_{-14}^{+21}  M_{⊙} and 66_{-18}^{+17}  M_{⊙} (90% credible intervals). We infer that the primary black hole mass lies within the gap produced by (pulsational) pair-instability supernova processes, with only a 0.32% probability of being below 65  M_{⊙}. We calculate the mass of the remnant to be 142_{-16}^{+28}  M_{⊙}, which can be considered an intermediate mass black hole (IMBH). The luminosity distance of the source is 5.3_{-2.6}^{+2.4}  Gpc, corresponding to a redshift of 0.82_{-0.34}^{+0.28}. The inferred rate of mergers similar to GW190521 is 0.13_{-0.11}^{+0.30}  Gpc^{-3} yr^{-1}.

Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Advanced Virgo and KAGRA

We present our current best estimate of the plausible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next several years, with the intention of providing information to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals for the third (O3), fourth (O4) and fifth observing (O5) runs, including the planned upgrades of the Advanced LIGO and Advanced Virgo detectors. We study the capability of the network to determine the sky location of the source for gravitational-wave signals from the inspiral of binary systems of compact objects, that is binary neutron star, neutron star-black hole, and binary black hole systems. The ability to localize the sources is given as a sky-area probability, luminosity distance, and comoving volume. The median sky localization area (90% credible region) is expected to be a few hundreds of square degrees for all types of binary systems during O3 with the Advanced LIGO and Virgo (HLV) network. The median sky localization area will improve to a few tens of square degrees during O4 with the Advanced LIGO, Virgo, and KAGRA (HLVK) network. During O3, the median localization volume (90% credible region) is expected to be on the order of 10 5 , 10 6 , 10 7 Mpc 3 for binary neutron star, neutron star-black hole, and binary black hole systems, respectively. The localization volume in O4 is expected to be about a factor two smaller than in O3. We predict a detection count of 1 - 1 + 12 ( 10 - 10 + 52 ) for binary neutron star mergers, of 0 - 0 + 19 ( 1 - 1 + 91 ) for neutron star-black hole mergers, and 17 - 11 + 22 ( 79 - 44 + 89 ) for binary black hole mergers in a one-calendar-year observing run of the HLV network during O3 (HLVK network during O4). We evaluate sensitivity and localization expectations for unmodeled signal searches, including the search for intermediate mass black hole binary mergers.

Search for Subsolar Mass Ultracompact Binaries in Advanced LIGO's Second Observing Run

We present a search for subsolar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational-wave candidates consistent with subsolar mass ultracompact binaries with at least one component between 0.2  M_{⊙}-1.0  M_{⊙}. We use the null result to constrain the binary merger rate of (0.2  M_{⊙}, 0.2  M_{⊙}) binaries to be less than 3.7×10^{5}  Gpc^{-3} yr^{-1} and the binary merger rate of (1.0  M_{⊙}, 1.0  M_{⊙}) binaries to be less than 5.2×10^{3}  Gpc^{-3} yr^{-1}. Subsolar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with subsolar masses. Assuming a particular primordial black hole (PBH) formation model, we constrain a population of merging 0.2  M_{⊙} black holes to account for less than 16% of the dark matter density and a population of merging 1.0  M_{⊙} black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict subsolar mass binaries.

Tests of General Relativity with GW170817

The recent discovery by Advanced LIGO and Advanced Virgo of a gravitational wave signal from a binary neutron star inspiral has enabled tests of general relativity (GR) with this new type of source. This source, for the first time, permits tests of strong-field dynamics of compact binaries in the presence of matter. In this Letter, we place constraints on the dipole radiation and possible deviations from GR in the post-Newtonian coefficients that govern the inspiral regime. Bounds on modified dispersion of gravitational waves are obtained; in combination with information from the observed electromagnetic counterpart we can also constrain effects due to large extra dimensions. Finally, the polarization content of the gravitational wave signal is studied. The results of all tests performed here show good agreement with GR.

Constraining the p-Mode-g-Mode Tidal Instability with GW170817

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.

GW170817: Measurements of Neutron Star Radii and Equation of State

On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function p(ρ) of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as R_{1}=10.8_{-1.7}^{+2.0}  km for the heavier star and R_{2}=10.7_{-1.5}^{+2.1}  km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than 1.97  M_{⊙} as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain R_{1}=11.9_{-1.4}^{+1.4}  km and R_{2}=11.9_{-1.4}^{+1.4}  km at the 90% credible level. Finally, we obtain constraints on p(ρ) at supranuclear densities, with pressure at twice nuclear saturation density measured at 3.5_{-1.7}^{+2.7}×10^{34}  dyn cm^{-2} at the 90% level.

GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral

On August 17, 2017 at 12∶41:04 UTC the Advanced LIGO and Advanced Virgo gravitational-wave detectors made their first observation of a binary neutron star inspiral. The signal, GW170817, was detected with a combined signal-to-noise ratio of 32.4 and a false-alarm-rate estimate of less than one per 8.0×10^{4}  years. We infer the component masses of the binary to be between 0.86 and 2.26  M_{⊙}, in agreement with masses of known neutron stars. Restricting the component spins to the range inferred in binary neutron stars, we find the component masses to be in the range 1.17-1.60  M_{⊙}, with the total mass of the system 2.74_{-0.01}^{+0.04}M_{⊙}. The source was localized within a sky region of 28  deg^{2} (90% probability) and had a luminosity distance of 40_{-14}^{+8}  Mpc, the closest and most precisely localized gravitational-wave signal yet. The association with the γ-ray burst GRB 170817A, detected by Fermi-GBM 1.7 s after the coalescence, corroborates the hypothesis of a neutron star merger and provides the first direct evidence of a link between these mergers and short γ-ray bursts. Subsequent identification of transient counterparts across the electromagnetic spectrum in the same location further supports the interpretation of this event as a neutron star merger. This unprecedented joint gravitational and electromagnetic observation provides insight into astrophysics, dense matter, gravitation, and cosmology.

GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence

On August 14, 2017 at 10∶30:43 UTC, the Advanced Virgo detector and the two Advanced LIGO detectors coherently observed a transient gravitational-wave signal produced by the coalescence of two stellar mass black holes, with a false-alarm rate of ≲1 in 27 000 years. The signal was observed with a three-detector network matched-filter signal-to-noise ratio of 18. The inferred masses of the initial black holes are 30.5_{-3.0}^{+5.7}M_{⊙} and 25.3_{-4.2}^{+2.8}M_{⊙} (at the 90% credible level). The luminosity distance of the source is 540_{-210}^{+130}  Mpc, corresponding to a redshift of z=0.11_{-0.04}^{+0.03}. A network of three detectors improves the sky localization of the source, reducing the area of the 90% credible region from 1160   deg^{2} using only the two LIGO detectors to 60  deg^{2} using all three detectors. For the first time, we can test the nature of gravitational-wave polarizations from the antenna response of the LIGO-Virgo network, thus enabling a new class of phenomenological tests of gravity.