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.

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

The 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, 5σ 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 ∼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;

We present a novel algorithm for scheduling the observations of time-domain imaging surveys. Our Integer Linear Programming approach optimizes an observing plan for an entire night by assigning targets to temporal blocks, enabling strict control of the number of exposures obtained per field and minimizing filter changes. A subsequent optimization step minimizes slew times between each observation. Our optimization metric selfconsistently weights contributions from time-varying airmass, seeing, and sky brightness to maximize the transient discovery rate. We describe the implementation of this algorithm on the surveys of the Zwicky Transient Facility and present its on-sky performance.

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

We report the first plausible optical electromagnetic counterpart to a (candidate) binary black hole merger. Detected by the Zwicky Transient Facility, the electromagnetic flare is consistent with expectations for a kicked binary black hole merger in the accretion disk of an active galactic nucleus [B. McKernan, K. E. S. Ford, I. Bartos et al., Astrophys. J. Lett. 884, L50 (2019)] and is unlikely [< Oð0.01%Þ)] due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernova and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be < Oð0.1%Þ. If the flare is associated with S190521g, we find plausible values of total mass M BBH ∼ 100 M ⊙ , kick velocity v k ∼ 200 km s −1 at θ ∼ 60°in a disk with aspect ratio H=a ∼ 0.01 (i.e., disk height H at radius a) and gas density ρ ∼ 10 −10 g cm −3. The merger could have occurred at a disk migration trap (a ∼ 700r g ; r g ≡ GM SMBH =c 2 , where M SMBH is the mass of the active galactic nucleus supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a reencountering with the disk in ∼1.6 yrðM SMBH =10 8 M ⊙ Þða=10 3 r g Þ 3=2 .

Supermassive black hole binaries (SMBHBs) at sub-parsec separations should be common in galactic nuclei, as a result of frequent galaxy mergers. Hydrodynamical simulations of circumbinary discs predict strong periodic modulation of the mass accretion rate on time-scales comparable to the orbital period of the binary. As a result, SMB-HBs may be recognized by the periodic modulation of their brightness. We conducted a statistical search for periodic variability in a sample of 35,383 spectroscopically confirmed quasars in the photometric database of the Palomar Transient Factory (PTF). We analysed Lomb-Scargle periodograms and assessed the significance of our findings by modeling each individual quasar's variability as a damped random walk (DRW). We identified 50 quasars with significant periodicity beyond the DRW model, typically with short periods of a few hundred days. We find 33 of these to remain significant after a re-analysis of their periodograms including additional optical data from the intermediate-PTF and the Catalina Real-Time Transient Survey (CRTS). Assuming that the observed periods correspond to the redshifted orbital periods of SMBHBs, we conclude that our findings are consistent with a population of unequal-mass SMBHBs, with a typical mass ratio as low as q ≡ M 2 /M 1 ≈ 0.01.

Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). However, quantitative estimates of the CSM mass around such SNe are missing when the CSM material is optically thick. Specifically, current estimators are sensitive to uncertainties regarding the CSM density profile and the ejecta velocity. Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present new visible-light and X-ray observations of SN 2010jl (PTF 10aaxf), including the first detection of a SN in the hard X-ray band using NuSTAR. The total radiated luminosity of SN 2010jl is extreme, at least 9 × 10 50 erg. By modeling the visible-light data, we robustly show that the mass of the circumstellar material within ∼ 10 16 cm of the progenitor of SN 2010jl was in excess of 10 M ⊙ . This mass was likely ejected tens of years prior to the SN explosion. Our modeling suggests that the shock velocity during shock breakout was ∼ 6000 km s −1 , decelerating to ∼ 2600 km s −1 about two years after maximum light. Furthermore, our late-time NuSTAR and XMM spectra of the SN presumably provide the first direct measurement of SN shock velocity two years after the SN maximum light -measured to be in the range of 2000 km s −1 to 4500 km s −1 if the ions and electrons are in equilibrium, and > ∼ 2000 km s −1 if they are not in equilibrium. This measurement is in agreement with the shock velocity predicted by our modeling of the visible-light data. Our observations also show that the average radial density distribution of the CSM roughly follows an r −2 law. A possible explanation for the > ∼ 10 M ⊙ of CSM and the wind-like profile is that they are the result of multiple pulsational pair instability events prior to the SN explosion, separated from each other by years.