The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to ∼40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020–2021. Bridging the large sky areas with small field-of-view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground-based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to ∼200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7, respectively.

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 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 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 populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7 Moreover, <57 assuming flat (fading) evolution, at 90 the online terrestrial probability for each GW trigger, we find that <68 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)

Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the Universe. The duration and hardness distribution of GRBs has two clusters 1 , now understood to reflect (at least) two different progenitors 2 . Short-hard GRBs (SGRBs; T 90 < 2 s) arise from compact binary mergers, and long-soft GRBs (LGRBs; T 90 > 2 s) have been attributed to the collapse of peculiar massive stars (collapsars) 3 . The discovery of SN 1998bw/GRB 980425 (ref. 4 ) marked the first association of an LGRB with a collapsar, and AT 2017gfo (ref. 5 )/GRB 170817A/GW170817 (ref. 6 ) marked the first association of an SGRB with a binary neutron star merger, which also produced a gravitational wave. Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi satellite and the Interplanetary Network localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova, but which is consistent with being the supernova. Although the GRB duration is short (rest-frame T 90 of 0.65 s), our panchromatic follow-up data confirm a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets.

Cosmic neutrinos provide a unique window into the otherwise hidden mechanism of particle acceleration in astrophysical objects. The IceCube Collaboration recently reported the likely association of 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 likely association of a radio-emitting tidal disruption event, AT2019dsg, with a second high-energy neutrino. AT2019dsg was identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility. The probability of finding any coincident 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 multizone 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 petaelectronvolt neutrino production. Assuming that the association is genuine, our observations suggest that tidal disruption events with mildly relativistic outflows contribute to the cosmic neutrino flux. The tidal disruption event AT2019dsg is probably associated with a high-energy neutrino, suggesting that such events can contribute to the cosmic neutrino flux. The electromagnetic emission is explained in terms of a central engine, a photosphere and an extended synchrotron-emitting outflow.

Using selection criteria based on amplitude, time and color, we have identified 329 objects as known or candidate cataclysmic variable (CVs) during the first year of testing and operation of the Zwicky Transient Facility (ZTF). Of these, 90 are previously confirmed CVs, 218 are strong candidates based on the shape and color of their light curves obtained during 3-562 days of observations, and the remaining 21 are possible CVs but with too few data points to be listed as good candidates. Almost half the strong candidates are within 10 deg of the galactic plane, in contrast to most other large surveys which have avoided crowded fields. The available Gaia parallaxes are consistent with sampling the low mass transfer CVs, as predicted by population models. Our followup spectra have confirmed Balmer/helium emission lines in 27 objects, with four showing high excitation HeII emission, including candidates for an AM CVn, a polar and an intermediate polar. Our results demonstrate that a complete survey of the galactic plane is needed to accomplish an accurate determination of the number of CVs existing in the Milky Way.

We present multiwavelength modeling of the afterglow from the long γ -ray burst (GRB) 160625B using Markov Chain Monte Carlo techniques of the afterglowpy Python package. GRB 160625B is an extremely bright burst with a rich set of observations spanning from radio to γ -ray frequencies. These observations range from ∼0.1 days to >1000 days, thus making this event extremely well suited to such modeling. In this work we compare top-hat and Gaussian jet structure types in order to find best-fit values for the GRB jet collimation angle, viewing angle, and other physical parameters. We find that a Gaussian-shaped jet is preferred (2.7 σ –5.3 σ ) over the traditional top-hat model. Our estimate for the opening angle of the burst ranges from 1.°26 to 3.°90, depending on jet-shape model. We also discuss the implications that assumptions on jet shape, viewing angle, and particularly the participation a fraction of electrons have on the final estimation of GRB intrinsic energy release and the resulting energy budget of the relativistic outflow. Most notably, allowing the participation fraction to vary results in an estimated total relativistic energy of ∼10 53 erg. This is two orders of magnitude higher than when the total fraction is assumed to be unity; thus, this parameter has strong relevance for placing constraints on long GRB central engines, details of the circumburst media, and host environment.

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(2), 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 <57% (<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 <68% of putative kilonovae could be brighter than -16.6 mag. Comparing to model grids, we find that some kilonovae must have M-ej M, X-lan > 10(-4), or > 30 degrees 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%.

High cadence transient surveys are able to capture supernovae closer to their first light than before. Applying analytical models to such early emission, we can constrain the progenitor stars properties. In this paper, we present observations of SN2018fif (ZTF18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN2018fif was surrounded by relatively small amounts of circumstellar material (CSM) compared to all previous cases. This particularity, coupled with the high cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir Waxman and its extension to early times by Morag, Sapir Waxman. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN2018fif was a large red supergiant, with a radius of R=744.0_-128.0^+183.0 solar radii and an ejected mass of Mej=9.3_-5.8^+0.4 solar masses. Our model also gives information on the explosion epoch, the progenitor inner structure, the shock velocity and the extinction. The distribution of radii is double-peaked, with lower radii corresponding to lower values of the extinction, earlier recombination times and better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g. with ULTRASAT), will help break the radius-extinction degeneracy and independently determine both.

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of similar to 3760 square degrees/hour to median depths of g similar to 20.8 and r similar to 20.6 mag (AB, 5 sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The real-time pipeline employs a novel image-differencing algorithm, optimized for the detection of point-source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. We present some preliminary results of the calibration performance delivered by the real-time pipeline. The reconstructed astrometric accuracy per science image with respect to Gaia DR1 is typically 45 to 85 milliarcsec. This is the RMS per-axis on the sky for sources extracted with photometric S/N >= 10 and hence corresponds to the typical astrometric uncertainty down to this limit. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. The high end of these ranges corresponds to an airmass approaching similar to 2-the limit of the public survey. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae; rare flux transients; variable stars; eclipsing binaries; variability from active galactic nuclei;counterparts to gravitational wave sources; a more complete census of Type Ia supernovae; and solar-system objects.

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

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time-domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg(2) field of view and an 8 second readout time. It is well positioned in the development of time-domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities that provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r similar to 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF, including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei, and tidal disruption events, stellar variability, and solar system objects.

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

We report on a search for high-energy counterparts to fast radio bursts (FRBs) with the Fermi Gamma-ray Burst Monitor, Fermi Large Area Telescope, and the Neil Gehrels Swift Observatory Burst Alert Telescope. We find no significant associations for any of the 23 FRBs in our sample, but report upper limits to the high-energy fluence for each on timescales of 0.1, 1, 10, and 100 s. We report lower limits on the ratio of the radio to high-energy fluence, f r / f γ , for timescales of 0.1 and 100 s. We discuss the implications of our non-detections on various proposed progenitor models for FRBs, including analogs of giant pulses from the Crab pulsar and hyperflares from magnetars. This work demonstrates the utility of analyses of high-energy data for FRBs in tracking down the nature of these elusive sources.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48 inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg(2) field of view and 8 s readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory. We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.