Focus on the Search for Electromagnetic Counterparts to Gravitational-wave Sources

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline, and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

This Supplement provides supporting material for Abbott et al. (2016a). We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

The Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) discovered gravitational waves (GWs) from a binary black hole merger in 2015 September and may soon observe signals from neutron star mergers. There is considerable interest in searching for their faint and rapidly fading electromagnetic (EM) counterparts, though GW position uncertainties are as coarse as hundreds of square degrees. Because LIGO's sensitivity to binary neutron stars is limited to the local universe, the area on the sky that must be searched could be reduced by weighting positions by mass, luminosity, or star formation in nearby galaxies. Since GW observations provide information about luminosity distance, combining the reconstructed volume with positions and redshifts of galaxies could reduce the area even more dramatically. A key missing ingredient has been a rapid GW parameter estimation algorithm that reconstructs the full distribution of sky location and distance. We demonstrate the first such algorithm, which takes under a minute, fast enough to enable immediate EM follow-up. By combining the three-dimensional posterior with a galaxy catalog, we can reduce the number of galaxies that could conceivably host the event by a factor of 1.4, the total exposure time for the Swift X-ray Telescope by a factor of 2, the total exposure time for a synoptic optical survey by a factor of 2, and the total exposure time for a narrow-field optical telescope by a factor of 3. This encourages us to suggest a new role for small field of view optical instruments in performing targeted searches of the most massive galaxies within the reconstructed volumes.

This is a supplement to the Letter of Singer et al., in which we demonstrated a rapid algorithm for obtaining joint 3D estimates of sky location and luminosity distance from observations of binary neutron star mergers with Advanced LIGO and Virgo. We argued that combining the reconstructed volumes with positions and redshifts of possible host galaxies can provide large-aperture but small field of view instruments with a manageable list of targets to search for optical or infrared emission. In this Supplement, we document the new HEALPix-based file format for 3D localizations of gravitational-wave transients. We include Python sample code to show the reader how to perform simple manipulations of the 3D sky maps and extract ranked lists of likely host galaxies. Finally, we include mathematical details of the rapid volume reconstruction algorithm.

Using observations of the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL), we place upper limits on the gamma-ray and hard X-ray prompt emission associated with the gravitational wave event GW150914, which was discovered by the LIGO/Virgo Collaboration. The omnidirectional view of the INTEGRAL/SPI-ACS has allowed us to constrain the fraction of energy emitted in the hard X-ray electromagnetic component for the full high-probability sky region of LIGO triggers. Our upper limits on the hard X-ray fluence at the time of the event range from ${F}_{\gamma }=2\times {10}^{-8}$ erg cm⁻² to ${F}_{\gamma }={10}^{-6}$ erg cm⁻² in the 75 keV–2 MeV energy range for typical spectral models. Our results constrain the ratio of the energy promptly released in gamma-rays in the direction of the observer to the gravitational wave energy E${}_{\gamma }/$E${}_{\mathrm{GW}}\lt {10}^{-6}$. We discuss the implication of gamma-ray limits for the characteristics of the gravitational wave source, based on the available predictions for prompt electromagnetic emission.

We report the results of a deep search for an optical counterpart to the gravitational wave (GW) event GW150914, the first trigger from the Advanced LIGO GW detectors. We used the Dark Energy Camera (DECam) to image a 102 deg² area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in the i and z bands at 4–5, 7, and 24 days after the trigger. The median 5σ point-source limiting magnitudes of our search images are i = 22.5 and z = 21.8 mag. We processed the images through a difference-imaging pipeline using templates from pre-existing Dark Energy Survey data and publicly available DECam data. Due to missing template observations and other losses, our effective search area subtends 40 deg², corresponding to a 12% total probability in the initial map and 3% in the final map. In this area, we search for objects that decline significantly between days 4–5 and day 7, and are undetectable by day 24, finding none to typical magnitude limits of i = 21.5, 21.1, 20.1 for object colors (i − z) = 1, 0, −1, respectively. Our search demonstrates the feasibility of a dedicated search program with DECam and bodes well for future research in this emerging field.

The collapse of a stellar core is expected to produce gravitational waves (GWs), neutrinos, and in most cases a luminous supernova. Sometimes, however, the optical event could be significantly less luminous than a supernova and a direct collapse to a black hole, where the star just disappears, is possible. The GW event GW150914 was detected by the LIGO Virgo Collaboration via a burst analysis that gave localization contours enclosing the Large Magellanic Cloud (LMC). Shortly thereafter, we used DECam to observe 102 deg² of the localization area, including 38 deg² on the LMC for a missing supergiant search. We construct a complete catalog of LMC luminous red supergiants, the best candidates to undergo invisible core collapse, and collected catalogs of other candidates: less luminous red supergiants, yellow supergiants, blue supergiants, luminous blue variable stars, and Wolf–Rayet stars. Of the objects in the imaging region, all are recovered in the images. The timescale for stellar disappearance is set by the free-fall time, which is a function of the stellar radius. Our observations at 4 and 13 days after the event result in a search sensitive to objects of up to about 200 solar radii. We conclude that it is unlikely that GW150914 was caused by the core collapse of a relatively compact supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational waveform as best interpreted as a high mass binary black hole merger. We discuss how to generalize this search for future very nearby core-collapse candidates.

The intermediate Palomar Transient Factory (iPTF) autonomously responded to and promptly tiled the error region of the first gravitational-wave event GW150914 to search for an optical counterpart. Only a small fraction of the total localized region was immediately visible in the northern night sky, due both to Sun-angle and elevation constraints. Here, we report on the transient candidates identified and rapid follow-up undertaken to determine the nature of each candidate. Even in the small area imaged of 126 deg², after extensive filtering, eight candidates were deemed worthy of additional follow-up. Within two hours, all eight were spectroscopically classified by the Keck II telescope. Curiously, even though such events are rare, one of our candidates was a superluminous supernova. We obtained radio data with the Jansky Very Large Array and X-ray follow-up with the Swift satellite for this transient. None of our candidates appear to be associated with the gravitational-wave trigger, which is unsurprising given that GW150914 came from the merger of two stellar-mass black holes. This end-to-end discovery and follow-up campaign bodes well for future searches in this post-detection era of gravitational waves.

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σ). 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}_{-1.0}^{+1.5}\times {10}^{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 detection of the first gravitational wave (GW) transient GW150914 prompted an extensive campaign of follow-up observations at all wavelengths. Although no dedicated XMM-Newton observations have been performed, the satellite passed through the GW150914 error region during normal operations. Here we report the analysis of the data taken during these satellite slews performed two hours and two weeks after the GW event. Our data cover 1.1 and 4.8 deg² of the final GW localization region. No X-ray counterpart to GW150914 is found down to a sensitivity of 6 × 10⁻¹³ erg cm⁻² s⁻¹ in the 0.2–2 keV band. Nevertheless, these observations show the great potential of XMM-Newton slew observations for searching for the electromagnetic counterparts of GW events. A series of adjacent slews performed in response to a GW trigger would take ≲1.5 days to cover most of the typical GW credible region. We discuss this scenario and its prospects for detecting the X-ray counterpart of future GW detections.

The Fermi Large Area Telescope (LAT) has an instantaneous field of view (FoV) covering $\sim 1/5$ of the sky and it completes a survey of the entire sky in high-energy gamma-rays every 3 hr. It enables searches for transient phenomena over timescales from milliseconds to years. Among these phenomena could be electromagnetic counterparts to gravitational wave (GW) sources. In this paper, we present a detailed study of the LAT observations relevant to Laser Interferometer Gravitational-wave Observatory (LIGO) event GW150914, which is the first direct detection of gravitational waves and has been interpreted as being due to the coalescence of two stellar-mass black holes. The localization region for GW150914 was outside the LAT FoV at the time of the GW signal. However, as part of routine survey observations, the LAT observed the entire LIGO localization region within ∼70 minutes of the trigger and thus enabled a comprehensive search for a γ-ray counterpart to GW150914. The study of the LAT data presented here did not find any potential counterparts to GW150914, but it did provide limits on the presence of a transient counterpart above 100 MeV on timescales of hours to days over the entire GW150914 localization region.

We report the results of an extensive search through the AGILE data for a gamma-ray counterpart to the LIGO gravitational-wave (GW) event GW150914. Currently in spinning mode, AGILE has the potential of cover 80% of the sky with its gamma-ray instrument, more than 100 times a day. It turns out that AGILE came within a minute of the event time of observing the accessible GW150914 localization region. Interestingly, the gamma-ray detector exposed ∼65% of this region during the 100 s time intervals centered at −100 and +300 s from the event time. We determine a 2σ flux upper limit in the band 50 MeV–10 GeV, UL = 1.9 × 10⁻⁸ erg cm⁻² s⁻¹, obtained ∼300 s after the event. The timing of this measurement is the fastest ever obtained for GW150914, and significantly constrains the electromagnetic emission of a possible high-energy counterpart. We also carried out a search for a gamma-ray precursor and delayed emission over five timescales ranging from minutes to days: in particular, we obtained an optimal exposure during the interval −150/−30 s. In all these observations, we do not detect a significant signal associated with GW150914. We do not reveal the weak transient source reported by Fermi-GBM 0.4 s after the event time. However, even though a gamma-ray counterpart of the GW150914 event was not detected, the prospects for future AGILE observations of GW sources are decidedly promising.

We present the results of the optical follow-up conducted by the TOROS collaboration of the first gravitational-wave event GW150914. We conducted unfiltered CCD observations (0.35–1 μm) with the 1.5 m telescope at Bosque Alegre starting ∼2.5 days after the alarm. Given our limited field of view (∼100 arcmin²), we targeted 14 nearby galaxies that were observable from the site and were located within the area of higher localization probability. We analyzed the observations using two independent implementations of difference-imaging algorithms, followed by a Random-Forest-based algorithm to discriminate between real and bogus transients. We did not find any bona fide transient event in the surveyed area down to a 5σ limiting magnitude of r = 21.7 mag (AB). Our result is consistent with the LIGO detection of a binary black hole merger, for which no electromagnetic counterparts are expected, and with the expected rates of other astrophysical transients.

In view of the recent report by Connaughton et al., we analyze continuous time-tagged event (TTE) data of Fermi-gamma-ray burst monitor (GBM) around the time of the gravitational-wave event GW 150914. We find that after proper accounting for low-count statistics, the GBM transient event at 0.4 s after GW 150914 is likely not due to an astrophysical source, but consistent with a background fluctuation, removing the tension between the INTEGRAL/ACS non-detection and GBM. Additionally, reanalysis of other short GRBs shows that without proper statistical modeling the fluence of faint events is over-predicted, as verified for some joint GBM–ACS detections of short GRBs. We detail the statistical procedure to correct these biases. As a result, faint short GRBs, verified by ACS detections, with significances in the broadband light curve even smaller than that of the GBM–GW150914 event are recovered as proper non-zero source, while the GBM–GW150914 event is consistent with zero fluence.

We present radio follow-up observations carried out with the Karl G. Jansky Very Large Array during the first observing run (O1) of the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). A total of three gravitational-wave triggers were followed-up during the $\approx 4$ months of O1, from 2015 September to 2016 January. Two of these triggers, GW150914 and GW151226, are binary black hole (BH) merger events of high significance. A third trigger, G194575, was subsequently declared as an event of no interest (i.e., a false alarm). Our observations targeted selected optical transients identified by the intermediate Palomar Transient Factory in the Advanced LIGO error regions of the three triggers, and a limited region of the gravitational-wave localization area of G194575 not accessible to optical telescopes due to Sun constraints, where a possible high-energy transient was identified. No plausible radio counterparts to GW150914 and GW151226 were found, in agreement with expectations for binary BH mergers. We show that combining optical and radio observations is key to identifying contaminating radio sources that may be found in the follow-up of gravitational-wave triggers, such as emission associated with star formation and active galactic nuclei. We discuss our results in the context of the theoretical predictions for radio counterparts to gravitational-wave transients, and describe our future plans for the radio follow-up of Advanced LIGO (and Virgo) triggers.

We present a search, using KamLAND, a kiloton-scale anti-neutrino detector, for low-energy anti-neutrino events that were coincident with the gravitational-wave (GW) events GW150914 and GW151226, and the candidate event LVT151012. We find no inverse beta-decay neutrino events within ±500 s of either GW signal. This non-detection is used to constrain the electron anti-neutrino fluence and the total integrated luminosity of the astrophysical sources.

We have performed a blind search for a gamma-ray transient of arbitrary duration and energy spectrum around the time of the LIGO gravitational-wave event GW150914 with the six-spacecraft interplanetary network (IPN). Four gamma-ray bursts were detected between 30 hr prior to the event and 6.1 hr after it, but none could convincingly be associated with GW150914. No other transients were detected down to limiting 15–150 keV fluences of roughly 5 ×10⁻⁸–5 × 10⁻⁷ erg cm⁻². We discuss the search strategies and temporal coverage of the IPN on the day of the event and compare the spatial coverage to the region where GW150914 originated. We also report the negative result of a targeted search for the Fermi-GBM event reported in conjunction with GW150914.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational-wave events, GW150914 and GW151226, as well as LVT151012, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of ±500 s around the gravitational-wave detection time. Four neutrino candidates are found for GW150914, and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90% confidence level upper limits on the combined neutrino fluence for both gravitational-wave events, which depends on event energy and topologies. Considering the upward-going muon data set (1.6 GeV–100 PeV), the neutrino fluence limit for each gravitational-wave event is 14–37 (19–50) cm⁻² for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational-wave events range from 2.4 × 10⁴ to 7.0 × 10⁹ cm⁻².

The weak transient detected by the Fermi Gamma-ray Burst Monitor (GBM) 0.4 s after GW150914 has generated much speculation regarding its possible association with the black hole binary merger. Investigation of the GBM data by Connaughton et al. revealed a source location consistent with GW150914 and a spectrum consistent with a weak, short gamma-ray burst. Greiner et al. present an alternative technique for fitting background-limited data in the low-count regime, and call into question the spectral analysis and the significance of the detection of GW150914-GBM presented in Connaughton et al. The spectral analysis of Connaughton et al. is not subject to the limitations of the low-count regime noted by Greiner et al. We find Greiner et al. used an inconsistent source position and did not follow the steps taken in Connaughton et al. to mitigate the statistical shortcomings of their software when analyzing this weak event. We use the approach of Greiner et al. to verify that our original spectral analysis is not biased. The detection significance of GW150914-GBM is established empirically, with a false-alarm rate (FAR) of $\sim {10}^{-4}$ Hz. A post-trials false-alarm probability (FAP) of $2.2\times {10}^{-3}$ ($2.9\sigma$) of this transient being associated with GW150914 is based on the proximity in time to the gravitational-wave event of a transient with that FAR. The FAR and the FAP are unaffected by the spectral analysis that is the focus of Greiner et al.