An implementation of Deflate in Coq

The widely-used compression format "Deflate" is defined in RFC 1951 and is based on prefix-free codings and backreferences. There are unclear points about the way these codings are specified, and several sources for confusion in the standard. We tried to fix this problem by giving a rigorous mathematical specification, which we formalized in Coq. We produced a verified implementation in Coq which achieves competitive performance on inputs of several megabytes. In this paper we present the several parts of our implementation: a fully verified implementation of canonical prefix-free codings, which can be used in other compression formats as well, and an elegant formalism for specifying sophisticated formats, which we used to implement both a compression and decompression algorithm in Coq which we formally prove inverse to each other -- the first time this has been achieved to our knowledge. The compatibility to other Deflate implementations can be shown empirically. We furthermore discuss some of the difficulties, specifically regarding memory and runtime requirements, and our approaches to overcome them

GW170817: Implications for the Stochastic Gravitational-Wave Background from Compact Binary Coalescences

The LIGO Scientific and Virgo Collaborations have announced the event GW170817, the first detection of gravitational waves from the coalescence of two neutron stars. The merger rate of binary neutron stars estimated from this event suggests that distant, unresolvable binary neutron stars create a significant astrophysical stochastic gravitational-wave background. The binary neutron star component will add to the contribution from binary black holes, increasing the amplitude of the total astrophysical background relative to previous expectations. In the Advanced LIGO-Virgo frequency band most sensitive to stochastic backgrounds (near 25 Hz), we predict a total astrophysical background with amplitude ΩGW(f=25 Hz)=1.8-1.3+2.7×10-9 with 90% confidence, compared with ΩGW(f=25 Hz)=1.1-0.7+1.2×10-9 from binary black holes alone. Assuming the most probable rate for compact binary mergers, we find that the total background may be detectable with a signal-to-noise-ratio of 3 after 40 months of total observation time, based on the expected timeline for Advanced LIGO and Virgo to reach their design sensitivity

GW170608: Observation of a 19 Solar-mass Binary Black Hole Coalescence

On 2017 June 8 at 02:01:16.49 UTC, a gravitational-wave (GW) signal from the merger of two stellar-mass blackholes was observed by the two Advanced Laser Interferometer Gravitational-Wave Observatory detectors with anetwork signal-to-noise ratio of 13. This system is the lightest black hole binary so far observed, with componentmasses of 12+7-2M⊙7+2-2 (90% credible intervals). These lie in the range of measured black hole masses inlow-mass X-ray binaries, thus allowing us to compare black holes detected through GWs with electromagneticobservations. The source\u27s luminosity distance is 340+140-140corresponding to redshift -0.07+0.03003. We verify thatthe signal waveform is consistent with the predictions of general relativity

First narrow-band search for continuous gravitational waves from known pulsars in advanced detector data

Spinning neutron stars asymmetric with respect to their rotation axis are potential sources of continuous gravitational waves for ground-based interferometric detectors. In the case of known pulsars a fully coherent search, based on matched filtering, which uses the position and rotational parameters obtained from electromagnetic observations, can be carried out. Matched filtering maximizes the signal-to-noise (SNR) ratio, but a large sensitivity loss is expected in case of even a very small mismatch between the assumed and the true signal parameters. For this reason, narrow-band analysis methods have been developed, allowing a fully coherent search for gravitational waves from known pulsars over a fraction of a hertz and several spin-down values. In this paper we describe a narrow-band search of 11 pulsars using data from Advanced LIGO\u27s first observing run. Although we have found several initial outliers, further studies show no significant evidence for the presence of a gravitational wave signal. Finally, we have placed upper limits on the signal strain amplitude lower than the spin-down limit for 5 of the 11 targets over the bands searched; in the case of J1813-1749 the spin-down limit has been beaten for the first time. For an additional 3 targets, the median upper limit across the search bands is below the spin-down limit. This is the most sensitive narrow-band search for continuous gravitational waves carried out so far

Estimating the Contribution of Dynamical Ejecta in the Kilonova Associated with GW170817

The source of the gravitational-wave (GW) signal GW170817, very likely a binary neutron star merger, was also observed electromagnetically, providing the first multi-messenger observations of this type. The two-week-long electromagnetic (EM) counterpart had a signature indicative of an r-process-induced optical transient known as a kilonova. This Letter examines how the mass of the dynamical ejecta can be estimated without a direct electromagnetic observation of the kilonova, using GW measurements and a phenomenological model calibrated to numerical simulations of mergers with dynamical ejecta. Specifically, we apply the model to the binary masses inferred from the GW measurements, and use the resulting mass of the dynamical ejecta to estimate its contribution (without the effects of wind ejecta) to the corresponding kilonova light curves from various models. The distributions of dynamical ejecta mass range between Mej = 10-3 - 10-2 M⊙ for various equations of state, assuming that the neutron stars are rotating slowly. In addition, we use our estimates of the dynamical ejecta mass and the neutron star merger rates inferred from GW170817 to constrain the contribution of events like this to the r-process element abundance in the Galaxy when ejecta mass from post-merger winds is neglected. We find that if ≳10% of the matter dynamically ejected from binary neutron star (BNS) mergers is converted to r-process elements, GW170817-like BNS mergers could fully account for the amount of r-process material observed in the Milky Way

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.7M and 25.3-4.2+2.8M (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 deg2 using only the two LIGO detectors to 60 deg2 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

Search for Post-merger Gravitational Waves from the Remnant of the Binary Neutron Star Merger GW170817

The first observation of a binary neutron star (NS) coalescence by the Advanced LIGO and Advanced Virgo gravitational-wave (GW) detectors offers an unprecedented opportunity to study matter under the most extreme conditions. After such a merger, a compact remnant is left over whose nature depends primarily on the masses of the inspiraling objects and on the equation of state of nuclear matter. This could be either a black hole (BH) or an NS, with the latter being either long-lived or too massive for stability implying delayed collapse to a BH. Here, we present a search for GWs from the remnant of the binary NS merger GW170817 using data from Advanced LIGO and Advanced Virgo. We search for short- (≲1 s) and intermediate-duration (≲500 s) signals, which include GW emission from a hypermassive NS or supramassive NS, respectively. We find no signal from the post-merger remnant. Our derived strain upper limits are more than an order of magnitude larger than those predicted by most models. For short signals, our best upper limit on the root sum square of the GW strain emitted from 1-4 kHz is at 50% detection efficiency. For intermediate-duration signals, our best upper limit at 50% detection efficiency is for a millisecond magnetar model, and for a bar-mode model. These results indicate that post-merger emission from a similar event may be detectable when advanced detectors reach design sensitivity or with next-generation detectors

On the Progenitor of Binary Neutron Star Merger GW170817

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\u27s 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\u27s star formation history, provided the stellar populations are older than 1 Gyr

Search for Tensor Vector and Scalar Polarizations in the Stochastic Gravitational-Wave Background

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy densities of tensor, vector, and scalar modes at 95% credibility to Ω0T\u3c5.58×10-8, Ω0V\u3c6.35×10-8, and Ω0S\u3c1.08×10-7 at a reference frequency f0=25 Hz

GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2

We describe the observation of GW170104, a gravitational-wave signal produced by the coalescence of a pair of stellar-mass black holes. The signal was measured on January 4, 2017 at 10: 11: 58.6 UTC by the twin advanced detectors of the Laser Interferometer Gravitational-Wave Observatory during their second observing run, with a network signal-to-noise ratio of 13 and a false alarm rate less than 1 in 70 000 years. The inferred component black hole masses are 31.2(-6.0)(+8.4) M-circle dot and 19.4(-5.9)(+5.3)M(circle dot) (at the 90% credible level). The black hole spins are best constrained through measurement of the effective inspiral spin parameter, a mass-weighted combination of the spin components perpendicular to the orbital plane, chi(eff) = -0.12(-0.30)(+0.21) . This result implies that spin configurations with both component spins positively aligned with the orbital angular momentum are disfavored. The source luminosity distance is 880(-390)(+450) Mpc corresponding to a redshift of z = 0.18(-0.07)(+0.08) . We constrain the magnitude of modifications to the gravitational-wave dispersion relation and perform null tests of general relativity. Assuming that gravitons are dispersed in vacuum like massive particles, we bound the graviton mass to m(g) \u3c= 7.7 x 10(-23) eV/c(2). In all cases, we find that GW170104 is consistent with general relativity