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

Synoptic searches for the optical counterpart to a binary neutron star (BNS) or neutron star-black hole (NSBH) merger can pose significant challenges towards the discovery of kilonovae and performing multi-messenger science. In this work, we describe the advantage of a global multi-telescope network towards this end, with a particular focus on the key and complementary role the Dark Energy Camera (DECam) plays in multi-facility follow-up. We describe the Global Relay of Observatories Watching Transients Happen (GROWTH) Target-of-Opportunity (ToO) Marshal, a common web application we built to ingest events, plan observations, search for transient candidates, and retrieve performance summary statistics for all of the telescopes in our network. Our infrastructure enabled us to conduct observations of two events during O3a, S190426c and S190510g. Furthermore, our analysis of deep DECam observations of S190814bv conducted by the DESGW team, and access to a variety of global follow-up facilities allowed us to place meaningful constraints on the parameters of the kilonova and the merging binary. We emphasize the importance of a global telescope network in conjunction with a power telescope like DECam in performing searches for the counterparts to gravitational-wave sources.

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 deg2, 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.

Merging neutron stars offer an excellent laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart (EM170817) with gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic data set, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultrarelativistic jets. Instead, we suggest that breakout of a wide-angle, mildly relativistic cocoon engulfing the jet explains the low-luminosity gamma rays, the high-luminosity ultraviolet-optical-infrared, and the delayed radio and x-ray emission. We posit that all neutron star mergers may lead to a wide-angle cocoon breakout, sometimes accompanied by a successful jet and sometimes by a choked jet.

The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package, as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of interoperable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy Project.

We present a search for subsolar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational-wave candidates consistent with subsolar mass ultracompact binaries with at least one component between 0.2  M_{⊙}-1.0  M_{⊙}. We use the null result to constrain the binary merger rate of (0.2  M_{⊙}, 0.2  M_{⊙}) binaries to be less than 3.7×10^{5}  Gpc^{-3} yr^{-1} and the binary merger rate of (1.0  M_{⊙}, 1.0  M_{⊙}) binaries to be less than 5.2×10^{3}  Gpc^{-3} yr^{-1}. Subsolar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with subsolar masses. Assuming a particular primordial black hole (PBH) formation model, we constrain a population of merging 0.2  M_{⊙} black holes to account for less than 16% of the dark matter density and a population of merging 1.0  M_{⊙} black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict subsolar mass binaries.

On 17 August 2017, the Advanced LIGO and Virgo detectors observed the gravitational-wave event GW170817-a strong signal from the merger of a binary neutron-star system. Less than two seconds after the merger, a γ-ray burst (GRB 170817A) was detected within a region of the sky consistent with the LIGO-Virgo-derived location of the gravitational-wave source. This sky region was subsequently observed by optical astronomy facilities, resulting in the identification of an optical transient signal within about ten arcseconds of the galaxy NGC 4993. This detection of GW170817 in both gravitational waves and electromagnetic waves represents the first 'multi-messenger' astronomical observation. Such observations enable GW170817 to be used as a 'standard siren' (meaning that the absolute distance to the source can be determined directly from the gravitational-wave measurements) to measure the Hubble constant. This quantity represents the local expansion rate of the Universe, sets the overall scale of the Universe and is of fundamental importance to cosmology. Here we report a measurement of the Hubble constant that combines the distance to the source inferred purely from the gravitational-wave signal with the recession velocity inferred from measurements of the redshift using the electromagnetic data. In contrast to previous measurements, ours does not require the use of a cosmic 'distance ladder': the gravitational-wave analysis can be used to estimate the luminosity distance out to cosmological scales directly, without the use of intermediate astronomical distance measurements. We determine the Hubble constant to be about 70 kilometres per second per megaparsec. This value is consistent with existing measurements, while being completely independent of them. Additional standard siren measurements from future gravitational-wave sources will enable the Hubble constant to be constrained to high precision.