Chad Hanna

The first two years of electromagnetic follow-up with advanced LIGO and Virgo

We anticipate the first direct detections of gravitational waves (GWs) with Advanced LIGO and Virgo later this decade. Though this groundbreaking technical achievement will be its own reward, a still greater prize could be observations of compact binary mergers in both gravitational and electromagnetic channels simultaneously. During Advanced LIGO and Virgos first two years of operation, 2015 through 2016, we expect the global GW detector array to improve in sensitivity and livetime and expand from two to three detectors. We model the detection rate and the sky localization accuracy for binary neutron star (BNS) mergers across this transition. We have analyzed a large, astrophysically motivated source population using real-time detection and sky localization codes and higher-latency parameter estimation codes that have been expressly built for operation in the Advanced LIGO/Virgo era. We show that for most BNS events, the rapid sky localization, available about a minute after a detection, is as accurate as the full parameter estimation. We demonstrate that Advanced Virgo will play an important role in sky localization, even though it is anticipated to come online with only one-third as much sensitivity as the Advanced LIGO detectors. We find that the median 90% confidence region shrinks from ~500 deg^2 in 2015 to ~200 deg^2 in 2016. A few distinct scenarios for the first LIGO/Virgo detections emerge from our simulations.

TOWARD EARLY-WARNING DETECTION OF GRAVITATIONAL WAVES FROM COMPACT BINARY COALESCENCE

Rapid detection of compact binary coalescence (CBC) with a network of advanced gravitational-wave detectors will offer a unique opportunity for multi-messenger astronomy. Prompt detection alerts for the astronomical community might make it possible to observe the onset of electromagnetic emission from CBC. We demonstrate a computationally practical filtering strategy that could produce early-warning triggers before gravitational radiation from the final merger has arrived at the detectors.

PARAMETER ESTIMATION FOR BINARY NEUTRON-STAR COALESCENCES WITH REALISTIC NOISE DURING THE ADVANCED LIGO ERA

Advanced ground-based gravitational-wave (GW) detectors begin operation imminently. Their intended goal is not only to make the first direct detection of GWs, but also to make inferences about the source systems. Binary neutron-star mergers are among the most promising sources. We investigate the performance of the parameter-estimation (PE) pipeline that will be used during the first observing run of the Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) in 2015: we concentrate on the ability to reconstruct the source location on the sky, but also consider the ability to measure masses and the distance. Accurate, rapid sky localization is necessary to alert electromagnetic (EM) observatories so that they can perform follow-up searches for counterpart transient events. We consider PE accuracy in the presence of non-stationary, non-Gaussian noise. We find that the character of the noise makes negligible difference to the PE performance at a given signal-to-noise ratio. The source luminosity distance can only be poorly constrained, since the median 90% (50%) credible interval scaled with respect to the true distance is 0.85 (0.38). However, the chirp mass is well measured. Our chirp-mass estimates are subject to systematic error because we used gravitational-waveform templates without component spin to carry out inference on signals with moderate spins, but the total error is typically less than 10^(-3) M_☉. The median 90% (50%) credible region for sky localization is ~ 600 deg^2 (~150 deg^2), with 3% (30%) of detected events localized within 100 deg^2. Early aLIGO, with only two detectors, will have a sky-localization accuracy for binary neutron stars of hundreds of square degrees; this makes EM follow-up challenging, but not impossible.

Searching for gravitational waves from binary coalescence

We describe the implementation of a search for gravitational waves from compact binary coalescences in LIGO and Virgo data. This all-sky, all-time, multidetector search for binary coalescence has been used to search data taken in recent LIGO and Virgo runs. The search is built around a matched filter analysis of the data, augmented by numerous signal consistency tests designed to distinguish artifacts of non-Gaussian detector noise from potential detections. We demonstrate the search performance using Gaussian noise and data from the fifth LIGO science run and demonstrate that the signal consistency tests are capable of mitigating the effect of non-Gaussian noise and providing a sensitivity comparable to that achieved in Gaussian noise.

Analysis Framework for the Prompt Discovery of Compact Binary Mergers in Gravitational-wave Data

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 min 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.

Improving the sensitivity of a search for coalescing binary black holes with nonprecessing spins in gravitational wave data

We demonstrate for the first time a search pipeline with improved sensitivity to gravitational waves from coalescing binary black holes with spins aligned to the orbital angular momentum by the inclusion of spin effects in the search templates. We study the pipeline recovery of simulated gravitational wave signals from aligned-spin binary black holes added to real detector noise, comparing the pipeline performance with aligned-spin filter templates to the same pipeline with nonspinning filter templates. Our results exploit a three-parameter phenomenological waveform family that models the full inspiral-merger-ringdown coalescence and treats the effect of aligned spins with a single effective spin parameter χ. We construct template banks from these waveforms by a stochastic placement method and use these banks as filters in the recently developed gstlal search pipeline. We measure the observable volume of the analysis pipeline for binary black hole signals with M_(total) and χ∈[0,0.85]. We find an increase in observable volume of up to 45% for systems with 0.2≤χ≤0.85 with almost no loss of sensitivity to signals with 0≤χ≤0.2. We also show that the use of spinning templates in the search pipeline provides for more accurate recovery of the binary mass parameters as well as an estimate of the effective spin parameter. We demonstrate this analysis on 25.9 days of data obtained from the Hanford and Livingston detectors in LIGO’s fifth observation run.

UTILITY OF GALAXY CATALOGS FOR FOLLOWING UP GRAVITATIONAL WAVES FROM BINARY NEUTRON STAR MERGERS WITH WIDE-FIELD TELESCOPES

The first detections of gravitational waves from binary neutron star mergers with advanced LIGO and Virgo observatories are anticipated in the next five years. These detections could pave the way for multi-messenger gravitational-wave (GW) and electromagnetic (EM) astronomy if GW triggers are successfully followed up with targeted EM observations. However, GW sky localization is relatively poor, with expected localization areas of ~10-100 deg2; this presents a challenge for following up GW signals from compact binary mergers. Even for wide-field instruments, tens or hundreds of pointings may be required. Prioritizing pointings based on the relative probability of successful imaging is important since it may not be possible to tile the entire gravitational-wave localization region in a timely fashion. Galaxy catalogs were effective at narrowing down regions of the sky to search in initial attempts at joint GW/EM observations. The relatively limited range of initial GW instruments meant that few galaxies were present per pointing and galaxy catalogs were complete within the search volume. The next generation of GW detectors will have a 10-fold increase in range thereby increasing the expected number of galaxies per unit solid angle by a factor of ~1000. As an additional complication, catalogs will be highly incomplete. Nevertheless, galaxy catalogs can still play an important role in prioritizing pointings for the next era of GW searches. We show how to quantify the advantages of using galaxy catalogs to prioritize wide-field follow-ups as a function of only two parameters: the three-dimensional volume within the field of view of a telescope after accounting for the GW distance measurement uncertainty, and the fraction of the GW sky localization uncertainty region that can be covered with telescope pointings. We find that the use of galaxy catalogs can improve the probability of successful imaging by ~10% to ~300% relative to follow-up strategies that do not utilize such catalogs for the scenarios we considered. We determine that catalogs with a 75% completeness perform comparably to complete catalogs in most cases, while 33%-complete catalogs can lead to lower follow-up success rates than complete catalogs for small fields of view, though still providing an advantage over strategies that do not use a catalog at all.

Intense electromagnetic outbursts from collapsing hypermassive neutron stars

We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts.

Reconstructing the Sky Location of Gravitational-Wave Detected Compact Binary Systems: Methodology for Testing and Comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localization have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high latency and aimed at in-depth studies of all the parameters of a source, including sky position, and “triangulation-based” techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyze methods for assessing the self consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localization. We apply these methods to two existing sky localization techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with a total mass of ≤20M_⊙ and nonspinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor ≈20 smaller for fully coherent techniques than for the specific variant of the triangulation-based technique used during the last science runs, at the expense of a factor ≈1000 longer processing time.

Parameter estimation on gravitational waves from neutron-star binaries with spinning components

Inspiraling binary neutron stars are expected to be one of the most significant sources of gravitational-wave signals for the new generation of advanced ground-based detectors. We investigate how well we could hope to measure properties of these binaries using the Advanced LIGO detectors, which began operation in September 2015. We study an astrophysically motivated population of sources (binary components with masses