GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-infrared Counterpart to a Neutron Star-Black Hole Merger

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 BH merger. Due to the low false-alarm rate and the precise localization (23 deg2 at 90%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to an NSBH merger. To search for potential counterparts, the GROWTH Collaboration performed real-time image subtraction on six nights of public Dark Energy Camera images acquired in the 3 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 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 ej < 0.04 M o˙ at polar viewing angles, or M ej < 0.03 M o˙ if the opacity is κ < 2 cm2g-1. Assuming a tidal deformability for the NS 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 χ < 0.7 for mass ratios Q < 6, with weaker constraints for more compact NSs. © 2020. The American Astronomical Society

The First Post-Kepler Brightness Dips of KIC 8462852

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process

Kilonova Luminosity Function Constraints Based on Zwicky Transient Facility Searches for 13 Neutron Star Merger Triggers during O3

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 (O3). 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 area of 4480 deg², a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10⁻²⁵ yr⁻¹. The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (−16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than −16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day⁻¹ (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than −16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than 10⁻⁴, or φ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of −16 mag would constrain the maximum fraction of bright kilonovae to <25%

GROWTH on GW190425: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

The beginning of the third observing run by the network of gravitational-wave detectors has brought the discovery of many compact binary coalescences. Prompted by the detection of the first binary neutron star merger in this run (GW190425 / LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. As it was a single gravitational-wave detector discovery, the initial skymap spanned most of the sky observable from Palomar Observatory, the site of both instruments. Covering 8000 deg$^2$ of the inner 99\% of the initial skymap over the next two nights, corresponding to an integrated probability of 46\%, the ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered a total of 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% of the integrated probability based on the initial sky map. However, the revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. Out of the 338,646 ZTF transient "alerts" over the first two nights of observations, we narrowed this list to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod were particularly compelling given that their location, distance, and age were consistent with the gravitational-wave 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 photometrically or spectroscopically ruled-out as supernovae. Palomar Gattini-IR identified one fast evolving infrared transient after the merger, PGIR19bn, which was later spectroscopically classified as an M-dwarf flare. [abridged

GROWTH on S190425z: Searching Thousands of Square Degrees to Identify an Optical or Infrared Counterpart to a Binary Neutron Star Merger with the Zwicky Transient Facility and Palomar Gattini-IR

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 deg2 of the initial skymap over the next two nights, corresponding to 46% integrated probability, ZTF system achieved a depth of ≈21 m 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 skymap. The revised skymap issued the following day reduced these numbers to 21% for the ZTF 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 light curves 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

Target-of-opportunity observations of gravitational-wave events with Vera C. Rubin Observatory

The discovery of the electromagnetic counterpart to the binary neutron star (NS) merger GW170817 has opened the era of gravitational-wave multimessenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multiwavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of NS mergers and other gravitational-wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving NSs (∼tens per year) out to distances of several hundred megaparsecs. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of NS and other compact-object mergers, and yet unknown classes of gravitational-wave events

Kilonova Luminosity Function Constraints Based on Zwicky Transient Facility Searches for 13 Neutron Star Merger Triggers during O3

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 (O3). 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 area of 4480 deg2, a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10-25 yr-1. The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than-16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day-1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than-16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than &lt;57% (&lt;89%) of putative kilonovae could be brighter than-16.6 mag assuming flat evolution (fading by 1 mag day-1) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than &lt;68% of putative kilonovae could be brighter than-16.6 mag. Comparing to model grids, we find that some kilonovae must have M ej &lt; 0.03 M o˙, X lan &gt; 10-4, or φ &gt; 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of-16 mag would constrain the maximum fraction of bright kilonovae to &lt;25%