The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies.We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star (BNS), a neutron star -black hole binary (NSBH) and a binary black hole (BBH), where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with nonspinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations.We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior.We discuss the computational cost of these algorithms, and describe the general and problemspecific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence (CBC) parameter space.

No abstract

We describe a stream-based analysis pipeline to detect gravitational waves from the merger of binary neutron stars, binary black holes, and neutron-star-black-hole binaries within ∼ 1 minute of the arrival of the merger signal at Earth. Such low-latency detection is crucial for the prompt response by electromagnetic facilities in order to observe any fading electromagnetic counterparts that might be produced by mergers involving at least one neutron star. Even for systems expected not to produce counterparts, low-latency analysis of the data is useful for deciding when not to point telescopes, and as feedback to observatory operations. Analysts using this pipeline were the first to identify GW151226, the second gravitational-wave event ever detected. The pipeline also operates in an offline mode, in which it incorporates more refined information about data quality and employs acausal methods that are inapplicable to the online mode. The pipeline's offline mode was used in the detection of the first two gravitational-wave events, GW150914 and GW151226, as well as the identification of a third candidate, LVT151012.

Inferring astrophysical information from gravitational waves emitted by compact binaries is one of the key science goals of gravitational-wave astronomy. In order to reach the full scientific potential of gravitational-wave experiments, we require techniques to mitigate the cost of Bayesian inference, especially as gravitational-wave signal models and analyses become increasingly sophisticated and detailed. Reduced-order models (ROMs) of gravitational waveforms can significantly reduce the computational cost of inference by removing redundant computations. In this paper, we construct the first reduced-order models of gravitational-wave signals that include the effects of spin precession, inspiral, merger, and ringdown in compact object binaries and that are valid for component masses describing binary neutron star, binary black hole, and mixed binary systems. This work utilizes the waveform model known as "IMRPhenomPv2." Our ROM enables the use of a fast reduced-order quadrature (ROQ) integration rule which allows us to approximate Bayesian probability density functions at a greatly reduced computational cost. We find that the ROQ rule can be used to speed-up inference by factors as high as 300 without introducing systematic bias. This corresponds to a reduction in computational time from around half a year to half a day for the longest duration and lowest mass signals. The ROM and ROQ rules are available with the main inference library of the LIGO Scientific Collaboration, LALInference.

No abstract

We report on an all-sky search with the LIGO detectors for periodic gravitational waves in the frequency range 50 -1000 Hz and with the frequency's time derivative in the range ÿ1 10 ÿ8 Hz s ÿ1 to zero. Data from the fourth LIGO science run (S4) have been used in this search. Three different semicoherent methods of transforming and summing strain power from short Fourier transforms (SFTs) of the calibrated data have been used. The first, known as StackSlide, averages normalized power from each SFT. A ''weighted Hough'' scheme is also developed and used, which also allows for a multiinterferometer search. The third method, known as PowerFlux, is a variant of the StackSlide method in which the power is weighted before summing. In both the weighted Hough and PowerFlux methods, the weights are chosen according to the noise and detector antenna-pattern to maximize the signal-to-noise ratio. The respective advantages and disadvantages of these methods are discussed. Observing no evidence of periodic gravitational radiation, we report upper limits; we interpret these as limits on this radiation from isolated rotating neutron stars. The best population-based upper limit with 95% confidence on the gravitational-wave strain amplitude, found for simulated sources distributed isotropically across the sky and with isotropically distributed spin axes, is 4:28 10 ÿ24 (near 140 Hz). Strict upper limits are also obtained for small patches on the sky for best-case and worst-case inclinations of the spin axes.

After their successful first observing run (September 12, 2015 -January 12, 2016, the Advanced LIGO detectors were upgraded to increase their sensitivity for the second observing run (November 30, 2016 -August 26, 2017. The Advanced Virgo detector joined the second observing run on August 1, 2017. We discuss the updates that happened during this period in the GstLAL-based inspiral pipeline, which is used to detect gravitational waves from the coalescence of compact binaries both in low latency and an offline configuration. These updates include deployment of a zero-latency whitening filter to reduce the over-all latency of the pipeline by up to 32 seconds, incorporation of the Virgo data stream in the analysis, introduction of a single-detector search to analyze data from the periods when only one of the detectors is running, addition of new parameters to the likelihood ratio ranking statistic, increase in the parameter space of the search, and introduction of a template mass-dependent glitch-excision thresholding method.

We carry out two searches for periodic gravitational waves using the most sensitive few hours of data from the second LIGO science run. Both searches exploit fully coherent matched filtering and cover wide areas of parameter space, an innovation over previous analyses which requires considerable algorithm development and computational power. The first search is targeted at isolated, previously unknown neutron stars, covers the entire sky in the frequency band 160 -728.8 Hz, and assumes a frequency derivative of less than 4 10 ÿ10 Hz=s. The second search targets the accreting neutron star in the lowmass x-ray binary Scorpius X-1 and covers the frequency bands 464-484 Hz and 604-624 Hz as well as the two relevant binary orbit parameters. Because of the high computational cost of these searches we limit the analyses to the most sensitive 10 hours and 6 hours of data, respectively. Given the limited sensitivity and duration of the analyzed data set, we do not attempt deep follow-up studies. Rather we concentrate on demonstrating the data analysis method on a real data set and present our results as upper limits over large volumes of the parameter space. In order to achieve this, we look for coincidences in parameter space between the Livingston and Hanford 4-km interferometers. For isolated neutron stars our 95% confidence level upper limits on the gravitational wave strain amplitude range from 6:6 10 ÿ23 to 1 10 ÿ21 across the frequency band; for Scorpius X-1 they range from 1:7 10 ÿ22 to 1:3 10 ÿ21 across the two 20-Hz frequency bands. The upper limits presented in this paper are the first broadband wide parameter space upper limits on periodic gravitational waves from coherent search techniques. The methods developed here lay the foundations for upcoming hierarchical searches of more sensitive data which may detect astrophysical signals.