Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

We report the discovery and panchromatic followup observations of the young Type Ic supernova, SN2020oi, in M100, a grand design spiral galaxy at a mere distance of $14$ Mpc. We followed up with observations at radio, X-ray and optical wavelengths from only a few days to several months after explosion. The optical behaviour of the supernova is similar to those of other normal Type Ic supernovae. The event was not detected in the X-ray band but our radio observation revealed a bright mJy source ($L_{\nu} \approx 1.2 \times 10^{27} {\rm erg\,s}^{-1} {\rm Hz}^{-1}$). Given, the relatively small number of stripped envelope SNe for which radio emission is detectable, we used this opportunity to perform a detailed analysis of the comprehensive radio dataset we obtained. The radio emitting electrons initially experience a phase of inverse Compton cooling which leads to steepening of the spectral index of the radio emission. Our analysis of the cooling frequency points to a large deviation from equipartition at the level of $\epsilon_e/\epsilon_B \gtrsim 200$, similar to a few other cases of stripped envelope SNe. Our modeling of the radio data suggests that the shockwave driven by the SN ejecta into the circumstellar matter (CSM) is moving at $\sim 3\times 10^{4}\,{\rm km\,s}^{-1}$. Assuming a constant mass-loss from the stellar progenitor, we find that the mass-loss rate is $\dot{M} \approx 1.4\times 10^{-4}\,{M}_{\odot}\,{\rm yr}^{-1}$, for an assumed wind velocity of $1000\,{\rm km\,s}^{-1}$. The temporal evolution of the radio emission suggests a radial CSM density structure steeper than the standard $r^{-2}$.

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). At $z=0.0866$, ZTF19aawfbtg (SN2019hge) has more than 10 optical spectra at phases from $-41$ to $+103$ days. Most of these spectra match well with that of PTF10hgi, suggesting SN2019hge as the second case of a He-rich SLSN-I. Confirmation comes from its near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features, yielding a lower limit to the He-rich fraction of $\sim9$\%. It is suggested that the excitation of He I atoms requires non-thermal radiation. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. However, how the magnetar energy is transported to the outer layers is unclear. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities, around $-20.06\pm0.16$ mag ($g$-band, rest-frame), in comparison to the full range of $-20$ to $-23$ mag covered by SLSN-I.

We present ZTF20aajnksq (AT2020blt), a fast-fading ($\Delta r=2.4$ mag in $\Delta t=1.3$ days) red ($g-r\approx0.6$ mag) and luminous ($M_{1626}=-25.9$) optical transient at $z=2.9$ discovered by the Zwicky Transient Facility (ZTF). AT2020blt shares several features in common with afterglows to long-duration gamma-ray bursts (GRBs): (1) an optical light curve well-described by a broken power-law with a break at $t_\mathrm{j}=1$ day (observer-frame); (2) a luminous $(L_X = 10^{46}$ $\mathrm{erg}$ $\mathrm{s}^{-1})$ X-ray counterpart; and (3) luminous ($L_\nu = 4 \times 10^{31}$ $\mathrm{erg}$ $\mathrm{sec}^{-1}$ $\mathrm{Hz}^{-1}$ at 10 GHz) radio emission. However, no GRB was detected in the 0.74d between the last ZTF non-detection ($r > 20.64$) and the first ZTF detection ($r = 19.57$), with an upper limit on the isotropic-equivalent gamma-ray energy release of $E_{\gamma,\mathrm{iso}} < 7 \times 10^{52}$ erg. AT2020blt is thus the third afterglow-like transient discovered without a detected GRB counterpart (after PTF11agg and ZTF19abvizsw) and the second (after ZTF19abvizsw) with a redshift measurement. We conclude that the properties of AT2020blt are consistent with a classical (initial Lorentz factor $\Gamma_0 \gtrsim 100$) on-axis GRB that was missed by high-energy satellites. Furthermore, by estimating the rate of transients with light curves similar to that of AT2020blt in ZTF high-cadence data, we agree with previous results that there is no evidence for an afterglow-like phenomenon that is significantly more common than classical GRBs. We conclude by discussing the status and future of fast-transient searches in wide-field high-cadence optical surveys.

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)

We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture.

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.

Colors of Type Ia supernovae in the first few days after explosion provide a potential discriminant between different models. In this paper, we present $g-r$ colors of 65 Type Ia supernovae discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that $g-r$ colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties ($\sigma_\mathrm{noise}\sim\sigma_\mathrm{int}\sim$ 0.18 mag). Colors are nearly constant starting from 6 days after first light ($g-r\sim-0.15$ mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of $\sim-0.25$ mag day$^{-1}$) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of $^{56}$Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers ($M_\mathrm{He}=0.01\,M_\odot$) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by "double-detonation" models with larger helium masses. We finally identify a significant correlation between the early-time $g-r$ slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value=0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on $B-V$ colors.

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations.

We present ZTF18abvkwla (the "Koala"), a fast blue optical transient discovered in the Zwicky Transient Facility (ZTF) One-Day Cadence (1DC) Survey. ZTF18abvkwla has a number of features in common with the groundbreaking transient AT2018cow: blue colors at peak ($g-r\approx-0.5$ mag), a short rise time from half-max of under two days, a decay time to half-max of only three days, a high optical luminosity ($M_{g,\mathrm{peak}}\approx-20.6$mag), a hot ($\gtrsim 40,000$K) featureless spectrum at peak light, and a luminous radio counterpart. At late times ($\Delta t>80$d) the radio luminosity of ZTF18abvkwla ($\nu L_\nu \gtrsim 10^{40}$erg/s at 10 GHz, observer-frame) is most similar to that of long-duration gamma-ray bursts (GRBs). The host galaxy is a dwarf starburst galaxy ($M\approx5\times10^{8}M_\odot$, $\mathrm{SFR}\approx7 M_\odot$/yr) that is moderately metal-enriched ($\log\mathrm{[O/H]} \approx 8.5$), similar to the hosts of GRBs and superluminous supernovae. As in AT2018cow, the radio and optical emission in ZTF18abvkwla likely arise from two separate components: the radio from fast-moving ejecta ($\Gamma \beta c >0.38c$) and the optical from shock-interaction with confined dense material ($<0.07M_\odot$ in $\sim 10^{15}$cm). Compiling transients in the literature with $t_\mathrm{rise} <5$d and $M_\mathrm{peak}<-20$mag, we find that a significant number are engine-powered, and suggest that the high peak optical luminosity is directly related to the presence of this engine. From 18 months of the 1DC survey, we find that transients in this rise-luminosity phase space are at least two to three orders of magnitude less common than CC SNe. Finally, we discuss strategies for identifying such events with future facilities like the Large Synoptic Survey Telescope, and prospects for detecting accompanying X-ray and radio emission.

We present ZTF18aaqjovh (SN 2018bvw), a high-velocity ("broad-lined") stripped-envelope (Type Ic) supernova (Ic-BL SN) discovered in the Zwicky Transient Facility one-day cadence survey. ZTF18aaqjovh shares a number of features in common with engine-driven explosions: the photospheric velocity and the shape of the optical light curve are very similar to that of the Type Ic-BL SN 1998bw, which was associated with a low-luminosity gamma-ray burst (LLGRB) and had relativistic ejecta. However, the radio luminosity of ZTF18aaqjovh is almost two orders of magnitude fainter than that of SN 1998bw at the same velocity phase, and the shock velocity is at most mildly relativistic (v=0.06-0.4c). A search of high-energy catalogs reveals no compelling GRB counterpart to ZTF18aaqjovh, and the limit on the prompt GRB luminosity of $L_{\gamma,\mathrm{iso}} \approx 1.6 \times 10^{48}$ erg/sec excludes a classical GRB but not an LLGRB. Altogether, ZTF18aaqjovh represents another transition event between engine-driven SNe associated with GRBs and "ordinary" Ic-BL SNe.

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.

While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress towards this promise. Here present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, as well as longer rise times and higher disruption rates compared to the rest of the sample. The Bowen fluorescence mechanism requires a high density which can be reached at smaller radii, which in turn yields longer diffusion timescales. Thus, the difference in rise times suggests the pre-peak TDE light curves are governed not by the fallback timescale, but instead by the diffusion of photons through the tidal debris. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities and the lowest rates. We also report, for the first time, the detection of soft X-ray flares from a TDE on day timescales. Based on the fact the flares peak at a luminosity similar to the optical/UV blackbody luminosity, we attribute them to brief glimpses through a reprocessing layer that otherwise obscures the inner accretion flow.

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 field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

The Zwicky Transient Facility (ZTF) is performing a three-day cadence survey of the visible Northern sky (~3$\pi$). The transient candidates found in this survey are announced via public alerts. As a supplementary product ZTF is also conducting a large spectroscopic campaign: the ZTF Bright Transient Survey (BTS). The goal of the BTS is to spectroscopically classify all extragalactic transients brighter than 18.5 mag at peak brightness and immediately announce those classifications to the public. Extragalactic discoveries from ZTF are predominantly Supernovae (SNe). The BTS is the largest flux-limited SN survey to date. Here we present a catalog of the761 SNe that were classified during the first nine months of the survey (2018 Apr. 1 to 2018 Dec. 31). The BTS SN catalog contains redshifts based on SN template matching and spectroscopic host galaxy redshifts when available. Based on this data we perform an analysis of the redshift completeness of local galaxy catalogs, dubbed as the Redshift Completeness Fraction (RCF; the number of SN host galaxies with known spectroscopic redshift prior to SN discovery divided by the total number of SN hosts). In total, we identify the host galaxies of 512 Type Ia supernovae, 227 of which have known spectroscopic redshifts, yielding an RCF estimate of $44\% \pm1\%$. We find a steady decrease in the RCF with increasing distance in the local universe. For z<0.05, or ~200 Mpc, we find RCF=0.6, which has important ramifications when searching for multimessenger astronomical events. Prospects for dramatically increasing the RCF are limited to new multi-fiber spectroscopic instruments, or wide-field narrowband surveys. We find that existing galaxy redshift catalogs are only $50\%$ complete at $r\approx16.9$ mag. Pushing this limit several magnitudes deeper will pay huge dividends when searching for electromagnetic counterparts to gravitational wave events.

We present photometry and spectroscopy of four hydrogen-poor luminous supernovae discovered during the two-month science commissioning and early operations of the Zwicky Transient Facility (ZTF) survey. Three of these objects, SN2018bym (ZTF18aapgrxo), SN2018avk (ZTF18aaisyyp) and SN2018bgv (ZTF18aavrmcg) resemble typical SLSN-I spectroscopically, while SN2018don (ZTF18aajqcue) may be an object similar to SN2007bi experiencing considerable host galaxy reddening, or an intrinsically long-lived, luminous and red SN Ic. We analyze the light curves, spectra, and host galaxy properties of these four objects and put them in context of the population of SLSN-I. SN2018bgv stands out as the fastest-rising SLSN-I observed to date, with a rest-frame g-band rise time of just 10 days from explosion to peak -- if it is powered by magnetar spin-down, the implied ejecta mass is only ~1 M$_{\odot}$. SN2018don also displays unusual properties -- in addition to its red colors and comparatively massive host galaxy, the light curve undergoes some of the strongest light curve undulations post-peak seen in a SLSN-I, which we speculate may be due to interaction with circumstellar material. We discuss the promises and challenges of finding SLSNe in large-scale surveys like ZTF given the observed diversity in the population.

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.

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 (2017). Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model. We find that the progenitor of SN2018fif was a large red supergiant, with a radius of R=1174_{-81}^{+208} solar radii, and an ejected mass of M=5.6_{-1.0}^{+9.1} solar masses. Our model also gives information on the explosion epoch, the progenitor inner structure, the shock velocity and the extinction. The large radius differs from previously modeled objects, and the difference could be either intrinsic or due to the relatively small amount of CSM around SN2018fif, perhaps making it a "cleaner" candidate for applying shock-cooling analytical models.

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.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

We introduce AXS (Astronomy eXtensions for Spark), a scalable open-source astronomical data analysis framework built on Apache Spark, a widely used industry-standard engine for big data processing. Building on capabilities present in Spark, AXS aims to enable querying and analyzing almost arbitrarily large astronomical catalogs using familiar Python/AstroPy concepts, DataFrame APIs, and SQL statements. We achieve this by i) adding support to Spark for efficient on-line positional cross-matching and ii) supplying a Python library supporting commonly-used operations for astronomical data analysis. To support scalable cross-matching, we developed a variant of the ZONES algorithm (Gray et al. 2004) capable of operating in distributed, shared-nothing architecture. We couple this to a data partitioning scheme that enables fast catalog cross-matching and handles the data skew often present in deep all-sky data sets. The cross-match and other often-used functionalities are exposed to the end users through an easy-to-use Python API. We demonstrate AXS' technical and scientific performance on SDSS, ZTF, Gaia DR2, and AllWise catalogs. Using AXS we were able to perform on-the-fly cross-match of Gaia DR2 (1.8 billion rows) and AllWise (900 million rows) data sets in ~ 30 seconds. We discuss how cloud-ready distributed systems like AXS provide a natural way to enable comprehensive end-user analyses of large datasets such as LSST.

We present a cadence optimization strategy to unveil a large population of kilonovae using optical imaging alone. These transients are generated during binary neutron star and potentially neutron star-black hole mergers and are electromagnetic counterparts to gravitational-wave signals detectable in nearby events with Advanced LIGO, Advanced Virgo, and other interferometers that will come online in the near future. Discovering a large population of kilonovae will allow us to determine how heavy element production varies with the intrinsic parameters of the merger and across cosmic time. The rate of binary neutron star mergers is still uncertain, but only few (less than 15) events with associated kilonovae may be detectable per year within the horizon of next-generation ground-based interferometers. The rapid evolution (hours to days) at optical/infrared wavelengths, relatively low luminosity, and the low volumetric rate of kilonovae makes their discovery difficult, especially during blind surveys of the sky. We propose future large surveys to adopt a rolling cadence in which g-i observations are taken nightly for blocks of 10 consecutive nights. With the current baseline2018a cadence designed for the Large Synoptic Survey Telescope (LSST), less than 7.5 poorly-sampled kilonovae are expected to be detected in both the Wide Fast Deep (WFD) and Deep Drilling Fields (DDF) surveys per year, under optimistic assumptions on their rate, duration, and luminosity. We estimate the proposed strategy to return up to about 272 GW170817-like kilonovae throughout the LSST WFD survey, discovered independently from gravitational-wave triggers.

The Large Synoptic Survey Telescope (LSST) will open a discovery frontier for faint and fast transients with its ability to detect variable flux components down to $\sim$24.5 mag in a $\sim$30 second exposure. Spectroscopic follow-up of such phenomena - which are necessary for understanding the physics of stellar explosions - can require a rapid response and several hours with a 8-10m telescope, making it both expensive and difficult to acquire. The future Extremely Large Telescopes (ELTs) would be able to provide not only spectroscopy but capabilities such as spectropolarimetry and high-resolution diffraction-limited imaging that would contribute to future advances in our physical understanding of stellar explosions. In this white paper we focus on several specific scientific impacts in the field of explosive transient astrophysics that will be generated by the combination of LSST's discovery abilities and ELTs' follow-up capacities. First, we map the uncharted frontier of discovery phase-space in terms of intrinsic luminosity and timescales for explosive transients, where we expect the unexpected. We then focus on six areas with open science questions for known transients: the progenitors of thermonuclear supernovae (SNe), mass loss prior to core collapse, asymmetry in stellar explosions, light echoes, high-$z$ transients, and strongly lensed SNe. We conclude with a brief discussion of the practical aspects of ELT & LSST synergy.

Supernovae are among the most powerful and influential explosions in the universe. They are also ideal multi-messenger laboratories to study extreme astrophysics. However, many fundamental properties of supernovae related to their diverse progenitor systems and explosion mechanisms remain poorly constrained. Here we outline how late-time spectroscopic observations obtained during the nebular phase (several months to years after explosion), made possible with the next generation of Extremely Large Telescopes, will facilitate transformational science opportunities and rapidly accelerate the community towards our goal of achieving a complete understanding of supernova explosions. We highlight specific examples of how complementary GMT and TMT instrumentation will enable high fidelity spectroscopy from which the line profiles and luminosities of elements tracing mass loss and ejecta can be used to extract kinematic and chemical information with unprecedented detail, for hundreds of objects. This will provide uniquely powerful constraints on the evolutionary phases stars may experience approaching a supernova explosion; the subsequent explosion dynamics; their nucleosynthesis yields; and the formation of compact objects that may act as central engines.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST.

The Zwicky Transient Facility (ZTF) survey generates real-time alerts for optical transients, variables, and moving objects discovered in its wide-field survey. We describe the ZTF alert stream distribution and processing (filtering) system. The system uses existing open-source technologies developed in industry: Kafka, a real-time streaming platform, and Avro, a binary serialization format. The technologies used in this system provide a number of advantages for the ZTF use case, including (1) built-in replication, scalability, and stream rewind for the distribution mechanism; (2) structured messages with strictly enforced schemas and dynamic typing for fast parsing; and (3) a Python-based stream processing interface that is similar to batch for a familiar and user-friendly plug-in filter system, all in a modular, primarily containerized system. The production deployment has successfully supported streaming up to 1.2 million alerts or roughly 70 GB of data per night, with each alert available to a consumer within about 10 s of alert candidate production. Data transfer rates of about 80,000 alerts/minute have been observed. In this paper, we discuss this alert distribution and processing system, the design motivations for the technology choices for the framework, performance in production, and how this system may be generally suitable for other alert stream use cases, including the upcoming Large Synoptic Survey Telescope.

The Zwicky Transient Facility is a large optical survey in multiple filters producing hundreds of thousands of transient alerts per night. We describe here various machine learning (ML) implementations and plans to make the maximal use of the large data set by taking advantage of the temporal nature of the data, and further combining it with other data sets. We start with the initial steps of separating bogus candidates from real ones, separating stars and galaxies, and go on to the classification of real objects into various classes. Besides the usual methods (e.g., based on features extracted from light curves) we also describe early plans for alternate methods including the use of domain adaptation, and deep learning. In a similar fashion we describe efforts to detect fast moving asteroids. We also describe the use of the Zooniverse platform for helping with classifications through the creation of training samples, and active learning. Finally we mention the synergistic aspects of ZTF and LSST from the ML perspective.

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 ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma 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 realtime 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. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. 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.

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 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 which 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 $\sim$ 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 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 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). 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.

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.

We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over 10^10 possible asteroid trajectories in stacks of the order 10-15 4K x 4K images in under a minute using a single consumer grade GPU. We apply this algorithm to data from the 2015 campaign of the High Cadence Transient Survey (HiTS) obtained with the Dark Energy Camera (DECam). We find 39 previously unknown Kuiper Belt Objects in the 150 square degrees of the survey. Comparing these asteroids to an existing model for the inclination distribution of the Kuiper Belt we demonstrate that we recover a KBO population above our detection limit consistent with previous studies. Software used in this analysis is made available as an open source package.

Recently we have witnessed the first multi-messenger detection of colliding neutron stars through Gravitational Waves (GWs) and Electromagnetic (EM) waves (GW170817), thanks to the joint efforts of LIGO/Virgo and Space/Ground-based telescopes. In this paper, we report on the RATIR followup observation strategies and show the results for the trigger G194575. This trigger is not of astrophysical interest; however, is of great interests to the robust design of a followup engine to explore large sky error regions. We discuss the development of an image-subtraction pipeline for the 6-color, optical/NIR imaging camera RATIR. Considering a two band ($i$ and $r$) campaign in the Fall of 2015, we find that the requirement of simultaneous detection in both bands leads to a factor $\sim$10 reduction in false alarm rate, which can be further reduced using additional bands. We also show that the performance of our proposed algorithm is robust to fluctuating observing conditions, maintaining a low false alarm rate with a modest decrease in system efficiency that can be overcome utilizing repeat visits. Expanding our pipeline to search for either optical or NIR detections (3 or more bands), considering separately the optical $riZ$ and NIR $YJH$ bands, should result in a false alarm rate $\approx 1\%$ and an efficiency $\approx 90\%$. RATIR's simultaneous optical/NIR observations are expected to yield about one candidate transient in the vast 100 $\mathrm{deg^2}$ LIGO error region for prioritized followup with larger aperture telescopes.

We constrain the minimum variability timescales for 938 GRBs observed by the Fermi/GBM instrument prior to July 11, 2012. The tightest constraints on progenitor radii derived from these timescales are obtained from light curves in the hardest energy channel. In the softer bands -- or from measurements of the same GRBs in the hard X-rays from Swift -- we show that variability timescales tend to be a factor 2--3 longer. Applying a survival analysis to account for detections and upper limits, we find median minimum timescale in the rest frame for long-duration and short-duration GRBs of 45 ms and 10 ms, respectively. Fewer than 10% of GRBs show evidence for variability on timescales below 2 ms. These shortest timescales require Lorentz factors $\gtrsim 400$ and imply typical emission radii $R \approx 1 {\times} 10^{14}$ cm for long-duration GRBs and $R \approx 3 {\times} 10^{13}$ cm for short-duration GRBs. We discuss implications for the GRB fireball model and investigate whether GRB minimum timescales evolve with cosmic time.

We develop a robust technique to determine the minimum variability timescale for Gamma-ray Burst light curves, utilizing Haar wavelets. Our approach averages over the data for a given GRB, providing an aggregate measure of signal variation while also retaining sensitivity to narrow pulses within complicated time-series. In contrast to previous studies using wavelets, which simply define the minimum timescale in reference to the measurement noise floor, our approach identifies the signature of temporally-smooth features in the wavelet scaleogram and then additionally identifies a break in the scaleogram on longer timescales as signature of a true, temporally-unsmooth light curve feature or features. We apply our technique to the large sample of Swift GRB Gamma-ray light curves and for the first time -- due to the presence of a large number of GRBs with measured redshift -- determine the distribution of minimum variability timescales in the source frame. We find a median minimum timescale for long-duration GRBs in the source frame of $\Delta t_{\rm min}=0.5$ s, with the shortest timescale found being on the order of 10 ms. This short timescale suggests a compact central engine ($3 \times 10^3$ km). We discuss further implications for the GRB fireball model and present a tantalizing correlation between minimum timescale and redshift, which may in part be due to cosmological time-dilation.