L. P. Singer

Astropy: A community Python package for astronomy

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as flexible image transport system (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

Fermi GBM Observations of LIGO Gravitational Wave event GW150914

With an instantaneous view of 70% of the sky, the Fermi Gamma-ray Burst Monitor (GBM) is an excellent partner in the search for electromagnetic counterparts to gravitational-wave (GW) events. GBM observations at the time of the Laser Interferometer Gravitational-wave Observatory (LIGO) event GW150914 reveal the presence of a weak transient above 50 keV, 0.4 s after the GW event, with a false-alarm probability of 0.0022 (2.9(sigma)). This weak transient lasting 1 s was not detected by any other instrument and does not appear to be connected with other previously known astrophysical, solar, terrestrial, or magnetospheric activity. Its localization is ill-constrained but consistent with the direction of GW150914. The duration and spectrum of the transient event are consistent with a weak short gamma-ray burst (GRB) arriving at a large angle to the direction in which Fermi was pointing where the GBM detector response is not optimal. If the GBM transient is associated with GW150914, then this electromagnetic signal from a stellar mass black hole binary merger is unexpected. We calculate a luminosity in hard X-ray emission between 1 keV and 10 MeV of 1.8(sup +1.5, sub -1.0) x 10(exp 49) erg/s. Future joint observations of GW events by LIGO/Virgo and Fermi GBM could reveal whether the weak transient reported here is a plausible counterpart to GW150914 or a chance coincidence, and will further probe the connection between compact binary mergers and short GRBs.

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.

A radio counterpart to a neutron star merger

GROWTH observations of GW170817 The gravitational wave event GW170817 was caused by the merger of two neutron stars (see the Introduction by Smith). In three papers, teams associated with the GROWTH (Global Relay of Observatories Watching Transients Happen) project present their observations of the event at wavelengths from x-rays to radio waves. Evans et al. used space telescopes to detect GW170817 in the ultraviolet and place limits on its x-ray flux, showing that the merger generated a hot explosion known as a blue kilonova. Hallinan et al. describe radio emissions generated as the explosion slammed into the surrounding gas within the host galaxy. Kasliwal et al. present additional observations in the optical and infrared and formulate a model for the event involving a cocoon of material expanding at close to the speed of light, matching the data at all observed wavelengths. Science, this issue p. 1565, p. 1579, p. 1559; see also p. 1554 Radio observations constrain the energy and geometry of relativistic material ejected from a binary neutron star merger. Gravitational waves have been detected from a binary neutron star merger event, GW170817. The detection of electromagnetic radiation from the same source has shown that the merger occurred in the outskirts of the galaxy NGC 4993, at a distance of 40 megaparsecs from Earth. We report the detection of a counterpart radio source that appears 16 days after the event, allowing us to diagnose the energetics and environment of the merger. The observed radio emission can be explained by either a collimated ultrarelativistic jet, viewed off-axis, or a cocoon of mildly relativistic ejecta. Within 100 days of the merger, the radio light curves will enable observers to distinguish between these models, and the angular velocity and geometry of the debris will be directly measurable by very long baseline interferometry.

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.

A mildly relativistic wide-angle outflow in the neutron-star merger event GW170817

GW170817 was the first gravitational-wave detection of a binary neutron-star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 megaparsecs. It has been proposed that the observed γ-ray, X-ray and radio emission is due to an ultra-relativistic jet being launched during the merger (and successfully breaking out of the surrounding material), directed away from our line of sight (off-axis). The presence of such a jet is predicted from models that posit neutron-star mergers as the drivers of short hard-γ-ray bursts. Here we report that the radio light curve of GW170817 has no direct signature of the afterglow of an off-axis jet. Although we cannot completely rule out the existence of a jet directed away from the line of sight, the observed γ-ray emission could not have originated from such a jet. Instead, the radio data require the existence of a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high-velocity tail of the neutron-rich material that was ejected dynamically during the merger, or a cocoon of material that breaks out when a jet launched during the merger transfers its energy to the dynamical ejecta. Because the cocoon model explains the radio light curve of GW170817, as well as the γ-ray and X-ray emission (and possibly also the ultraviolet and optical emission), it is the model that is most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron-star mergers, giving rise to a hitherto unidentified population of radio, ultraviolet, X-ray and γ-ray transients in the local Universe.

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.

Rapid Bayesian position reconstruction for gravitational-wave transients

Within the next few years, Advanced LIGO and Virgo should detect gravitational waves from binary neutron star and neutron star-black hole mergers. These sources are also predicted to power a broad array of electromagnetic transients. Because the electromagnetic signatures can be faint and fade rapidly, observing them hinges on rapidly inferring the sky location from the gravitational-wave observations. Markov chain Monte Carlo methods for gravitational-wave parameter estimation can take hours or more. We introduce BAYESTAR, a rapid, Bayesian, non-Markov chain Monte Carlo sky localization algorithm that takes just seconds to produce probability sky maps that are comparable in accuracy to the full analysis. Prompt localizations from BAYESTAR will make it possible to search electromagnetic counterparts of compact binary mergers.

DISCOVERY AND REDSHIFT OF AN OPTICAL AFTERGLOW IN 71 deg2: iPTF13bxl AND GRB 130702A

We report the discovery of the optical afterglow of the γ-ray burst (GRB) 130702A, identified upon searching 71 deg^2 surrounding the Fermi Gamma-ray Burst Monitor (GBM) localization. Discovered and characterized by the intermediate Palomar Transient Factory, iPTF13bxl is the first afterglow discovered solely based on a GBM localization. Real-time image subtraction, machine learning, human vetting, and rapid response multi-wavelength follow-up enabled us to quickly narrow a list of 27,004 optical transient candidates to a single afterglow-like source. Detection of a new, fading X-ray source by Swift and a radio counterpart by CARMA and the Very Large Array confirmed the association between iPTF13bxl and GRB 130702A. Spectroscopy with the Magellan and Palomar 200 inch telescopes showed the afterglow to be at a redshift of z = 0.145, placing GRB 130702A among the lowest redshift GRBs detected to date. The prompt γ-ray energy release and afterglow luminosity are intermediate between typical cosmological GRBs and nearby sub-luminous events such as GRB 980425 and GRB 060218. The bright afterglow and emerging supernova offer an opportunity for extensive panchromatic follow-up. Our discovery of iPTF13bxl demonstrates the first observational proof-of-principle for ~10 Fermi-iPTF localizations annually. Furthermore, it represents an important step toward overcoming the challenges inherent in uncovering faint optical counterparts to comparably localized gravitational wave events in the Advanced LIGO and Virgo era.

Galaxy Strategy for Ligo-Virgo Gravitational Wave Counterpart Searches

In this work we continue a line of inquiry begun in Kanner et al. which detailed a strategy for utilizing telescopes with narrow fields of view, such as the Swift X-ray Telescope (XRT), to localize gravity wave (GW) triggers from LIGO/Virgo. If one considers the brightest galaxies that produce ~50% of the light, then the number of galaxies inside typical GW error boxes will be several tens. We have found that this result applies both in the early years of Advanced LIGO when the range is small and the error boxes large, and in the later years when the error boxes will be small and the range large. This strategy has the beneficial property of reducing the number of telescope pointings by a factor 10 to 100 compared with tiling the entire error box. Additional galaxy count reduction will come from a GW rapid distance estimate which will restrict the radial slice in search volume. Combining the bright galaxy strategy with a convolution based on anticipated GW localizations, we find that the searches can be restricted to about 18±5 galaxies for 2015, about 23±4 for 2017, and about 11±2 for 2020. This assumes a distance localization at or near the putative NS-NS merger range for each target year, and these totals are integrated out to the range. Integrating out to the horizon would roughly double the totals. For nearer localizations the totals would decrease. The galaxy strategy we present in this work will enable numerous sensitive optical and X-ray telescopes with small fields of view to participate meaningfully in searches wherein the prospects for rapidly fading afterglow place a premium on a fast response time.