We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently associated one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the association of a radio-emitting tidal disruption event (AT2019dsg) with another high-energy neutrino, identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility (ZTF). The probability of finding any radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. The association suggests that tidal disruption events contribute to the cosmic neutrino flux. Unlike previous work which considered the rare subset of tidal disruption events with relativistic jets, our observations of AT2019dsg suggest an empirical model with a mildly-relativistic outflow.

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.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019-05-10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a BNS merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 +- 92 Mpc. We triggered Target of Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4m Blanco Telescope at CTIO in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hours of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was re-classified 1.7 days after the merger, after our observations were completed, with a "binary neutron star merger" probability reduced from 98% to 42% in favor of a "terrestrial" classification.

The discovery of a transient kilonova following the gravitational-wave event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up across the electromagnetic spectrum. In the Southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4-m telescope is well-suited to this task, as it is able to cover wide-fields quickly while still achieving the depths required to find kilonovae like the one accompanying GW170817 to $\sim$500 Mpc, the binary neutron star horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility followup by the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration, we describe the observations and automated data movement, data reduction, candidate discovery, and vetting pipeline of our target-of-opportunity DECam observations of S190426c, the first possible neutron star--black hole merger detected via gravitational waves. Starting 7.5hr after S190426c, over 11.28\,hr of observations, we imaged an area of 525\,deg$^2$ ($r$-band) and 437\,deg$^2$ ($z$-band); this was 16.3\% of the total original localization probability and nearly all of the probability density visible from the Southern hemisphere. The machine-learning based pipeline was optimized for fast turnaround, delivering transient candidates for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hours before the end of our observations. Our observations yielded no detection of a bona fide counterpart to $m_z = 22.5$ and $m_r = 22.9$ at the 5$\sigma$ level of significance, consistent with the refined LVC positioning. We view these observations and rapid inferencing as an important real-world test for this novel end-to-end wide-field pipeline.

There is significant interest in the models for production of short gamma-ray bursts. Until now, the number of known short gamma-ray bursts with multi-wavelength afterglows has been small. While the {\it Fermi} Gamma-Ray Burst Monitor detects many gamma-ray bursts relative to the Neil Gehrels {\it Swift} Observatory, the large localization regions makes the search for counterparts difficult. With the Zwicky Transient Facility recently achieving first light, it is now fruitful to use its combination of depth ($m_\textrm{AB} \sim 20.6$), field of view ($\approx$ 47 square degrees), and survey cadence (every $\sim 3$ days) to perform Target of Opportunity observations. We demonstrate this capability on GRB 180523B, which was recently announced by the {\it Fermi} Gamma-Ray Burst Monitor as a short gamma-ray burst. ZTF imaged $\approx$ 2900\,square degrees of the localization region, resulting in the coverage of 61.6\,\% of the enclosed probability over 2 nights to a depth of $m_\textrm{AB} \sim 20.5$. We characterized 14 previously unidentified transients, and none were found to be consistent with a short gamma-ray burst counterpart. This search with the Zwicky Transient Facility shows it is an efficient camera for searching for coarsely-localized short gamma-ray burst and gravitational-wave counterparts, allowing for a sensitive search with minimal interruption to its nominal cadence.

The Kitt Peak Electron Multiplying CCD (EMCCD) demonstrator is a new instrument that has been developed for use at the Kitt Peak National Observatory's 84-inch telescope. The EMCCD enables single-band optical imaging in the Sloan \textit{g} and \textit{r} bands and Johnson \textit{UVRI} filters. The EMCCD is chosen for its sub-electron effective read noise using large multiplicative gains. As a result, frame rates of greater than 1\,Hz are possible. The field-of-view is 4.4$^\prime$ $\times$ 4.4$^\prime$ and the pixel size is 0.259$^{\prime\prime}$. This camera, coupled with a fully roboticized telescope, is ideal for follow-up of short period, white dwarf binary candidates, as well as short duration transient and periodic sources identified by large field-of-view all-sky surveys such as the Zwicky Transient Facility.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.