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.

Search for Subsolar-Mass Ultracompact Binaries in Advanced LIGO's First Observing Run.

We present the first Advanced LIGO and Advanced Virgo search for ultracompact binary systems with component masses between 0.2 M_{⊙}-1.0 M_{⊙} using data taken between September 12, 2015 and January 19, 2016. We find no viable gravitational wave candidates. Our null result constrains the coalescence rate of monochromatic (delta function) distributions of nonspinning (0.2 M_{⊙}, 0.2 M_{⊙}) ultracompact binaries to be less than 1.0×10^{6} Gpc^{-3} yr^{-1} and the coalescence rate of a similar distribution of (1.0 M_{⊙}, 1.0 M_{⊙}) ultracompact binaries to be less than 1.9×10^{4} Gpc^{-3} yr^{-1} (at 90% confidence). Neither black holes nor neutron stars are expected to form below ∼1 M_{⊙} through conventional stellar evolution, though it has been proposed that similarly low mass black holes could be formed primordially through density fluctuations in the early Universe and contribute to the dark matter density. The interpretation of our constraints in the primordial black hole dark matter paradigm is highly model dependent; however, under a particular primordial black hole binary formation scenario we constrain monochromatic primordial black hole populations of 0.2 M_{⊙} to be less than 33% of the total dark matter density and monochromatic populations of 1.0 M_{⊙} to be less than 5% of the dark matter density. The latter strengthens the presently placed bounds from microlensing surveys of massive compact halo objects (MACHOs) provided by the MACHO and EROS Collaborations.

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.

A dust-enshrouded tidal disruption event with a resolved radio jet in a galaxy merger.

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5 × 1052 erg at infrared and radio wavelengths but was not luminous at optical or x-ray wavelengths. We interpret this as a TDE with much of its emission reradiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH.

Estimation of the mechanical properties of the eye through the study of its vibrational modes.

Measuring the eye's mechanical properties in vivo and with minimally invasive techniques can be the key for individualized solutions to a number of eye pathologies. The development of such techniques largely relies on a computational modelling of the eyeball and, it optimally requires the synergic interplay between experimentation and numerical simulation. In Astrophysics and Geophysics the remote measurement of structural properties of the systems of their realm is performed on the basis of (helio-)seismic techniques. As a biomechanical system, the eyeball possesses normal vibrational modes encompassing rich information about its structure and mechanical properties. However, the integral analysis of the eyeball vibrational modes has not been performed yet. Here we develop a new finite difference method to compute both the spheroidal and, specially, the toroidal eigenfrequencies of the human eye. Using this numerical model, we show that the vibrational eigenfrequencies of the human eye fall in the interval 100 Hz-10 MHz. We find that compressible vibrational modes may release a trace on high frequency changes of the intraocular pressure, while incompressible normal modes could be registered analyzing the scattering pattern that the motions of the vitreous humour leave on the retina. Existing contact lenses with embebed devices operating at high sampling frequency could be used to register the microfluctuations of the eyeball shape we obtain. We advance that an inverse problem to obtain the mechanical properties of a given eye (e.g., Young's modulus, Poisson ratio) measuring its normal frequencies is doable. These measurements can be done using non-invasive techniques, opening very interesting perspectives to estimate the mechanical properties of eyes in vivo. Future research might relate various ocular pathologies with anomalies in measured vibrational frequencies of the eye.

The unusual γ-ray burst GRB 101225A from a helium star/neutron star merger at redshift 0.33.

Long γ-ray bursts (GRBs) are the most dramatic examples of massive stellar deaths, often associated with supernovae. They release ultra-relativistic jets, which produce non-thermal emission through synchrotron radiation as they interact with the surrounding medium. Here we report observations of the unusual GRB 101225A. Its γ-ray emission was exceptionally long-lived and was followed by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the first 10 days, the optical emission evolved as an expanding, cooling black body, after which an additional component, consistent with a faint supernova, emerged. We estimate its redshift to be z = 0.33 by fitting the spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a merger of a helium star with a neutron star that underwent a common envelope phase, expelling its hydrogen envelope. The resulting explosion created a GRB-like jet which became thermalized by interacting with the dense, previously ejected material, thus creating the observed black body, until finally the emission from the supernova dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy.