The Zwicky Transient Facility ZTF is a new optical timedomain survey that uses the Palomar 48 inch Schmidt telescope. A custombuilt widefield camera provides a 47 degSUP2SUP 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 nearrealtime reduction to identify moving and varying objects. We outline the analysis pipelines data products and associated archive. Finally we present onsky performance analysis and first scientific results from commissioning and the early survey. ZTFs public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.
['Astrophysics - Instrumentation and Methods for Astrophysics']
The Zwicky Transient Facility System Overview Performance and First Results
{'The Zwicky Transient Facility: System Overview, Performance, and First Results': "ERIC C. BELLM, 1 SHRINIVAS R. KULKARNI, 2 MATTHEW J. GRAHAM, 2 RICHARD DEKANY, 3 ROGER M. SMITH, 3 REED RIDDLE, 3 FRANK J. MASCI, 4 GEORGE HELOU, 4 THOMAS A. PRINCE, 2 SCOTT M. ADAMS, 2 C. BARBARINO, 5 TOM BARLOW, 2 JAMES BAUER, 6 RON BECK, 4 JUSTIN BELICKI, 3 RAHUL BISWAS, 7 NADEJDA BLAGORODNOVA, 2 DENNIS BODEWITS, 6 BRYCE BOLIN, 1, 8 , ∗ VALERY BRINNEL, 9 TIM BROOKE, 4 BRIAN BUE, 10 MATTIA BULLA, 7 RICK BURRUSS, 3 S. BRADLEY CENKO, 11, 12 CHAN-KAO CHANG ( GLYPH<224> U ¥ ), 13 ANDREW CONNOLLY, 1 MICHAEL COUGHLIN, 14 JOHN CROMER, 3 VIRGINIA CUNNINGHAM, 6 KISHALAY DE, 2 ALEX DELACROIX, 3 VANDANA DESAI, 4 DMITRY A. DUEV, 2 GWENDOLYN EADIE, 15, 16, 1 , † TONY L. FARNHAM, 6 MICHAEL FEENEY, 3 ULRICH FEINDT, 7 DAVID FLYNN, 4 ANNA FRANCKOWIAK, 17 S. FREDERICK, 6 C. FREMLING, 2 AVISHAY GAL-YAM, 18 SUVI GEZARI, 6, 12 MATTEO GIOMI, 19 DANIEL A. GOLDSTEIN, 2, ‡ V. ZACH GOLKHOU, 1, 16, † ARIEL GOOBAR, 7 STEVEN GROOM, 4 EUGEAN HACOPIANS, 4 DAVID HALE, 3 JOHN HENNING, 3 ANNA Y. Q. HO, 2 DAVID HOVER, 3 JUSTIN HOWELL, 4 TIARA HUNG, 6 DANIELA HUPPENKOTHEN, 1 DAVID IMEL, 4 WING-HUEN IP ( I 8 GLYPH<220> ), 13, 20 ŽELJKO IVEZI ' C, 1 EDWARD JACKSON, 4 LYNNE JONES, 1 MARIO JURIC, 1 MANSI M. KASLIWAL, 2 S. KASPI, 21 STEPHEN KAYE, 3 MICHAEL S. P. KELLEY, 6 MAREK KOWALSKI, 9, 17 EMILY KRAMER, 22 THOMAS KUPFER, 23, 24, 2 WALTER LANDRY, 4 RUSS R. LAHER, 4 CHIEN-DE LEE, 13 HSING WEN LIN ( GLYPH<151> GLYPH<1> GLYPH<135> ), 25, 13 ZHONG-YI LIN ( GLYPH<151> GLYPH<224> ' ), 13 RAGNHILD LUNNAN, 5 MATTEO GIOMI, 19 ASHISH MAHABAL, 2, 26 PETER MAO, 3 ADAM A. MILLER, 27, 28 SERGE MONKEWITZ, 4 PATRICK MURPHY, 29 CHOW-CHOONG NGEOW, 13 JAKOB NORDIN, 19 PETER NUGENT, 30, 31 ERAN OFEK, 32 MARIA T. PATTERSON, 1 BRYAN PENPRASE, 33 MICHAEL PORTER, 3 LUDWIG RAUCH, 17 UMAA REBBAPRAGADA, 10 DAN REILEY, 3 MICKAEL RIGAULT, 34 HECTOR RODRIGUEZ, 3 JAN VAN ROESTEL, 35 BEN RUSHOLME, 4 JAKOB VAN SANTEN, 17 S. SCHULZE, 18 DAVID L. SHUPE, 4 LEO P. SINGER, 11, 12 MAAYANE T. SOUMAGNAC, 18 ROBERT STEIN, 17 JASON SURACE, 4 JESPER SOLLERMAN, 5 PAULA SZKODY, 15 F. TADDIA, 5 SCOTT TEREK, 4 ANGELA VAN SISTINE, 36 SJOERT VAN VELZEN, 6 W. THOMAS VESTRAND, 37 RICHARD WALTERS, 3 CHARLOTTE WARD, 6 QUAN-ZHI YE ( I GLYPH<201> GLYPH<215> ), 4, 2 PO-CHIEH YU ( GLYPH<222>/ GLYPH<145> ), 13 LIN YAN, 2 AND JEFFRY ZOLKOWER 3  \n1 DIRAC Institute, Department of Astronomy, University of Washington, 3910 15th Avenue NE, Seattle, WA 98195, USA 2 Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 3 Caltech Optical Observatories, California Institute of Technology, Pasadena, CA 91125, USA 4 Infrared Processing and Analysis Center, California Institute of Technology, MS 100-22, Pasadena, CA 91125, USA 5 The Oskar Klein Centre & Department of Astronomy, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden 6 Department of Astronomy, University of Maryland, College Park, MD 20742, USA 7 The Oskar Klein Centre, Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden 8 B612 Asteroid Institute, 20 Sunnyside Ave, Suite 427, Mill Valley, Ca 94941 9 Institute of Physics, Humboldt-Universität zu Berlin, Newtonstr. 15, 124 89 Berlin, Germany 10 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA 11 Astrophysics Science Division, NASA Goddard Space Flight Center, MC 661, Greenbelt, MD 20771, USA 12 Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA 13 Institute of Astronomy, National Central University, 32001, Taiwan 14 Division of Physics, Math, and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 15 Department of Astronomy, University of Washington, Seattle, WA 98195, USA 16 The eScience Institute, University of Washington, Seattle, WA 98195, USA 17 Deutsches Elektronensynchrotron, Platanenallee 6, D-15738, Zeuthen, Germany 18 Department of Particle Physics and Astrophysics, Weizmann Institute of Science 234 Herzl St., Rehovot, 76100, Israel 19  \nInstitute of Physics, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany  \n20 Space Science Institute, Macau University of Science and Technology, Macau  \n21 School of Physics & Astronomy and Wise Observatory, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel-Aviv 6997801,  \nIsrael  \n22 Jet Propulsion Laboratory, Pasadena, CA 91109, USA  \n23 Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA  \n24 Department of Physics, University of California, Santa Barbara, CA 93106, USA  \n25 Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA  \n26 Center for Data Driven Discovery, California Institute of Technology, Pasadena, CA 91125, USA  \n27 Center for Interdisciplinary Exploration and Research in Astrophysics and Department of Physics and Astronomy, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA  \nCorresponding author: Eric Bellm  \necbellm@uw.edu", 'BELLM ET AL.': "<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nMag  \nFigure 10. ZTF r -band light curves of asteroid (11014) Svatopluk, (34771) 2001 QO252, and (182312) 2001 OT59.Figure 11. ZTF observations of the outburst of comet Echeclus. The top panel shows the 1.5 arcmin FOV images of Echeclus, oriented North up and East to the left, on 12 Nov. 2017, 16 Nov. 2017, 04 Dec. 2017, 10 Dec. 2017, 15 Dec. 2017, 20 Dec. 2017, 27 Dec. 2017, 30 Dec. 2017, 31 Dec. 2017, 01 Jan. 2018, 02 Jan. 2018, and 06 Jan. 2018 (UT dates) left to right. The plot shows results from aperture photometry (7 arcsecond radius aperture) from the Echeclus data spanning dates between 12 November, 2017 through 16 January, 2018. We converted these to the equivalent Af ρ values (a proxy for dust production) in log-cm units for the corresponding dates. The magenta dashed line indicates the derived Af ρ value for a magnitude value corresponding to a bare nucleus. The images were primarily taken in the ZTF r -band, while those taken on the 12-16 Nov. 2017, or on or after 27 Dec. 2017 were ZTF g -band images.  \n<!-- image -->  \nWith the P48 focal plane now filled with CCDs, future sky surveys with the P48 will require substantial effort to achieve further performance improvements relative to ZTF, although further gains in angular resolution, wavelength coverage, and/or time sampling may be contemplated. Instead, most third-generation sky surveys will look to naturally scalable networks of small and medium telescopes distributed geographically, following the example of ASAS-SN, Las Cumbres Observatory, ATLAS, KMTNet, and BlackGEM. The  \nalternative model is large new monolithic facilities purposebuilt for time-domain observations, with LSST serving as the exemplar.  \nIndeed, while the survey characteristics are quite different, ZTF will serve as a useful precursor for LSST. ZTF will stream one million time-domain detections nightly using a prototype of the LSST alert distribution system, providing several years of community experience ahead of LSST's flood of ten million nightly alerts.  \nFigure 12. ZTF images of comet C/2016 R2 (PanSTARRS) in the g -band. Four epochs from the commissioning phase are shown: (a) 2017 Nov 11; (b) 2017 Dec 23; (c) 2018 Jan 13; and (d) 2018 Feb 19. The field of view is 37 ' × 21 ' and the projected sunward vector is along the x-axis. The images are logarithmically scaled, except near the background where they are linearly scaled, in order to enhance details in the tail.  \n<!-- image -->  \nBecause of its larger field of view, ZTF obtains on average about four times more observations of any area of the sky than LSST, and these visits are split among a smaller set of filters. The resulting finer time-sampling will enable earlier discovery of transients and better classification of events based on their lightcurves. Moreover, ZTF's smaller aperture means that all of the ZTF-discovered events are accessible for spectroscopic followup with moderate-aperture telescopes. In fact, ZTF's discovery rate of transients brighter than 21 st magnitude is greater than LSST's (Bellm 2016). ZTF should thus provide large samples of bright transients and variables that will be crucial for interpreting LSST's deeper and more challenging survey.  \nBased on observations obtained with the Samuel Oschin Telescope 48-inch and the 60-inch Telescope at the Palomar Observatory as part of the Zwicky Transient Facility project. Major funding has been provided by the U.S. National Science Foundation under Grant No. AST-1440341 and by the ZTF partner institutions: the California Institute of Technology, the Oskar Klein Centre, the Weizmann Institute of Sci-  \nence, the University of Maryland, the University of Washington, Deutsches Elektronen-Synchrotron, the University of Wisconsin-Milwaukee, and the TANGO Program of the University System of Taiwan.  \nThis work is partly based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. The data presented here were obtained in part with ALFOSC, which is provided by the Instituto de Astrofisica de Andalucia (IAA) under a joint agreement with the University of Copenhagen and NOTSA. This work is partly based on observations made with DOLoRes@TNG.  \nJ. Bauer, T. Farnham, and M. Kelley gratefully acknowledge the NASA/University of Maryland/MPC Augmentation through the NASA Planetary Data System Cooperative Agreement NNX16AB16A.  \nE. Bellm, B. Bolin, A. Connolly, V. Z. Golkhou, D. Huppenkothen, Z. Ivezi'c, L. Jones, M. Juric, and M. Patterson acknowledge support from the University of Washington Col-  \nege of Arts and Sciences, Department of Astronomy, and the DIRAC Institute. University of Washington's DIRAC Institute is supported through generous gifts from the Charles and Lisa Simonyi Fund for Arts and Sciences, and the Washington Research Foundation. M. Juric and A. Connolly acknowledge the support of the Washington Research Foundation Data Science Term Chair fund, and the UW Provost's Initiative in Data-Intensive Discovery.  \n- E. Bellm, A. Connolly, Z. Ivezi'c, L. Jones, M. Juric, and M. Patterson acknowledge support from the Large Synoptic Survey Telescope, which is supported in part by the National Science Foundation through Cooperative Agreement 1258333 managed by the Association of Universities for Research in Astronomy (AURA), and the Department of Energy under Contract No. DE-AC02-76SF00515 with the SLAC National Accelerator Laboratory. Additional LSST funding comes from private donations, grants to universities, and inkind support from LSSTC Institutional Members.\n- E. Bellm is supported in part by the NSF AAG grant 1812779 and grant #2018-0908 from the Heising-Simons Foundation.\n- B.T. Bolin acknowledges funding for the Asteroid Institute program provided by B612 Foundation, W.K. Bowes Jr. Foundation, P. Rawls Family Fund and two anonymous donors in addition to general support from the B612 Founding Circle.\n- M. Bulla acknowledges support from the Swedish Research Council (Vetenskapsradet) and the Swedish National Space Board.\n- C.-K. Chang, W.-H. Ip, C.-D. Lee, Z.-Y. Lin, C.-C. Ngeow and P.-C. Yu thank the funding from Ministry of Science and Technology (Taiwan) under grant 104-2923-M-008-004MY5, 104-2112-M-008-014-MY3, 105-2112-M-008-002MY3, 106-2811-M-008-081 and 106-2112-M-008-007.\n- A. Gal-Yam is supported by the EU via ERC grant No. 725161, the Quantum Universe I-Core program, the ISF, the BSF Transformative program and by a Kimmel award.\n- S. Gezari is supported in part by NSF CAREER grant 1454816 and NSF AAG grant 1616566.\n- D. A. Goldstein acknowledges support from Hubble Fellowship grant HST-HF2-51408.001-A. Support for Program number HST-HF2-51408.001-A is provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555.\n- M. M. Kasliwal and Q.-Z. Ye acknowledge support by the GROWTH (Global Relay of Observatories Watching Transients Happen) project funded by the National Science Foundation PIRE (Partnership in International Research and Education) program under Grant No 1545949.\n- A. A. Mahabal acknowledges support from the following grants: NSF AST-1749235, NSF-1640818 and NASA 16ADAP16-0232.\n- A. A. Miller is funded by the Large Synoptic Survey Telescope Corporation in support of the Data Science Fellowship Program.\n- E. Ofek is grateful for support by a grant from the Israeli Ministry of Science, ISF, Minerva, BSF, BSF transformative program, and the I-CORE Program of the Planning and Budgeting Committee and The Israel Science Foundation (grant No 1829/12).\n- M. Rigault is supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement no. 759194 USNAC).\n- J. Sollerman acknowledges support from the Knut and Alice Wallenberg Foundation.\n- M. T. Soumagnac acknowledges support by a grant from IMOS/ISA, the Ilan Ramon fellowship from the Israel Ministry of Science and Technology and the Benoziyo center for Astrophysics at the Weizmann Institute of Science.  \nPart of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.  \nFacilities: PO:1.2m, PO:1.5m, ARC, DCT, NOT", 'ABSTRACT': "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.", '1. INTRODUCTION': "Large optical sky surveys have served as engines of discovery throughout the history of astronomy. By cataloging large samples of astrophysical objects, these surveys provide literal and metaphorical finder charts for detailed followup observations with larger and more expensive telescopes.  \nIn the last century, among the most influential sky surveys were the Palomar Observatory Sky Surveys. Conducted with photographic plates using the wide-field Palomar 48inch Schmidt telescope (Harrington 1952), the first and second sky surveys (POSS-I, Minkowski & Abell 1963; POSSII, Reid et al. 1991) mapped the Northern Hemisphere sky and enabled fifty years of discovery. A digitized version 1 (Lasker 1994; Djorgovski et al. 1998) is still widely used today.  \nThe advent of solid-state charge coupled devices (CCDs) provided a huge leap in the quantum efficiency (QE) of astronomical cameras, enabling existing telescopes to reach greater depths with shorter exposures. Contemporaneous improvements in CCD controller readout time and computer processing speed have increased data volumes while allowing data processing to keep pace. Beginning especially with the Sloan Digital Sky Survey (SDSS; York et al. 2000) but also including the Optical Gravitational Lensing Experiment (OGLE; Udalski et al. 1992), the All-Sky Automated Survey (ASAS; Pojmanski 1997), the Lincoln Near-Earth Asteroid  \nResearch survey (LINEAR; Stokes et al. 2000), the Supernova Legacy Survey (SNLS; Astier et al. 2006), PalomarQuest (Djorgovski et al. 2008), the Catalina Sky Survey (CSS; Larson et al. 2003) and associated Catalina Real-Time Transient Survey (CRTS; Drake et al. 2009), Skymapper (Keller et al. 2007), PanSTARRS (Kaiser et al. 2010), the Palomar Transient Factory (PTF; Law et al. 2009), the AllSky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014), the Asteroid Terrestrial-impact Last Alert System (ATLAS; Tonry et al. 2018), the Korea Microlensing Telescope Network (KMTNet; Kim et al. 2016), the Dark Energy Survey (Dark Energy Survey Collaboration 2005; Dark Energy Survey Collaboration et al. 2016) and other surveys using the Dark Energy Camera (Flaugher et al. 2015), Hyper Suprime-Cam (Miyazaki et al. 2018; Aihara et al. 2018), and the Evryscope (Law et al. 2015), surveys exploited these new capabilities to improve a subset of depth, areal coverage, filter selection, and/or time-domain sampling.  \nHere we present a new sky survey, the Zwicky Transient Facility (ZTF) 2 . ZTF's new CCD camera for the first time fills the focal plane of the Palomar 48-inch Schmidt, providing three orders of magnitude improvement in survey speed relative to the photographic surveys, by virtue of higher QE and substantial reduction in time between exposures. If it could ignore daylight, ZTF could repeat the entire POSS survey in one day.  \nThis paper provides a general overview of the ZTF observing and data systems, describes the on-sky performance and public surveys, and presents initial results in transient, vari-  \nable, and solar system science. Additional papers discuss specific ZTF aspects in greater detail: Graham et al. (2019) describe the scientific objectives of ZTF. Dekany et al. (2019) provide an in-depth description of the design of the observing system. Bellm et al. (2019) discuss the ZTF surveys and scheduler. Masci et al. (2019) detail the ZTF data system. Patterson et al. (2019) present the alert distribution system employed by ZTF. Mahabal et al. (2019) discuss several applications of machine learning used by ZTF. Tachibana & Miller (2018) presents a new star/galaxy classifier developed for ZTF from the PanSTARRS DR1 catalog (Chambers et al. 2016; Flewelling et al. 2016). Kasliwal et al. (2019) describe a web-based interface used by the ZTF collaboration to identify, track, and follow up transients of interest.", '2. OBSERVING SYSTEM': 'The capability of a survey camera to discover astrophysical transients can be quantified by its volumetric survey speed: the spatial volume within which it can detect an object of given absolute magnitude, divided by the total time per exposure (Bellm 2016). This quantity combines limiting magnitude, field of view, and exposure and overhead times into a single metric capturing how quickly a survey can probe physical space for new events.  \nThe ZTF concept assumed reuse of the Palomar 48-inch Samuel Oschin Schmidt Telescope. The subsequent design of the ZTF observing system-the camera, telescope, and associated subsystems-then attempted to maximize the volumetric survey speed of the system within a fixed cost envelope subject to this constraint. This goal required maximizing the field of view of the camera while maintaining image quality, minimizing beam obstruction, and minimizing readout and slew overheads. The final design achieves more than an order of magnitude improvement in survey speed relative to PTF 3 .  \nDekany et al. (2019) describes the as-built observing system in greater detail.', '2.1. CCD Mosaic': "The P48 was designed to use 14-inch square photographic plates, providing a field of view of 43.56 deg 2 (Harrington 1952). Large-format 'wafer-scale' CCDs proved the most cost-effective means of filling this large area and had the additional advantage of minimizing losses due to chip gaps. Our goal of maximizing throughput while minimizing cost motivated our decision to survey primarily in filters near the peak quantum efficiency of standard silicon. We selected  \n3 ZTF's median overhead time is about 10.2 s compared to 42.0 s for PTF, which had median R -band limiting magnitudes of 20.7 mag in 60 s exposures. For a fiducial object with Mr = -19 mag, then, ˙ V -19 = 3 . 5 × 10 4 Mpc 3 s -1 for ZTF as built, a factor of 14.9 larger than for PTF (Bellm 2016).  \nFigure 1. Image of the ZTF focal plane. The top and bottom rows of 6k × 6k science CCDs have a single-layer anti-reflective coating, while the middle rows have a dual-layer coating. Four 2k × 2k CCDs are located on the perimeter of the mosaic; one serves as a guider while the remaining three control tip, tilt, and focus. North is up and East is left.  \n<!-- image -->  \nbackside-illuminated standard silicon CCD231-C6 devices from e2v, Inc. The 15 µ m pixels provided critical sampling of the expected 2.0 '' FWHM point spread function (§2.4) at a plate scale of 1.01 '' pixel -1 while moderating data volume. (This pixel scale also matched that of the PTF camera.) Half of ZTF's CCDs have a single-layer anti-reflective coating, while the other half has a dual-layer coating that provides improved quantum efficiency in the g and r bands (Figures 1 &2).  \nThe CCDs are nearly perfect cosmetically having only a few blocked columns. QE is uniform to a few percent on large scales. Response non-uniformity on short scales is 0.55% at 400 nm falling linearly to 0.3% at 650nm. Dark current is negligible in maximum exposure times contemplated (300 s). Well capacity is typically 350,000 e-, and charge transfer inefficiency is < 5ppm per pixel shift.  \nFour 2k × 2k CCDs located around the perimeter of the mosaic serve as guide and focus sensors. These are STAdesigned fully depleted thick CCDs that were delta-doped and multi-layer anti-reflection coated by the JPL Micro Devices Laboratory. Three are offset 1.45 mm beyond the plane of science CCDs to allow determination of tip, tilt, and focus by computing the square root of the 2nd moment of the out-of-focus images. The fourth in-focus CCD is used for guiding.  \nTable 1. Specifications of the ZTF Observing System  \nFigure 2. On-axis filter transmission for the ZTF g , r , and i -band filters (blue, orange, and red lines). Grey and green points are measurements of the quantum efficiencies of the CCDs with single- and double-layer anti-reflective coatings, respectively. Shaded regions show the range of these measurements, while grey and green lines show a model of the quantum efficiency for each configuration.  \n<!-- image -->  \nBecause the focus of a Schmidt telescope is located within the telescope tube itself, maximizing throughput requires minimizing the beam obstruction caused by the ZTF camera  \nFigure 3. Cutaway view of the ZTF cryostat.  \n<!-- image -->  \nand related components. We located the readout electronics (§2.3), shutter (§2.5), and filter exchanger (§2.7) outside the telescope tube.  \nThe cryostat can achieve its extraordinary compactness by a signal routing strategy based on a vacuum interface board, a printed circuit board having embedded traces and sandwiched between two O-rings that maintain vacuum (Figures 3 & 4). On the outer edge of the interface board, commercial off-the-shelf connectors mount flush against the side of the cryostat, behind the beam footprint of the front window. Low obscuration (22.4% including spiders) is achieved at the Schmidt prime focus despite the mosaic having compara-  \nFigure 4. Rear view of the vacuum interface board during cryostat assembly. The vacuum gasket and connectors are visible around the perimeter. Holes in the interior provide space for the CCD flex cables as well as the control rods used during assembly.  \n<!-- image -->  \nble diameter to other major CCD cameras: The ZTF mosaic measures 560 mm from corner to corner, similar in size to the Dark Energy Camera (525 mm diameter) and about half of the area of the Large Synoptic Survey Telescope camera (640 mm diameter).", '2.3. Readout Electronics': "Each four-CCD quadrant of the mosaic is operated completely independently by a sixteen-channel CCD controller, the 'Archon' made by Semiconductor Technology Associates (STA) with 100 MHz video processor employing digital correlated double sampling. A fifth Archon reads the three focus CCDs and guider though two channels each. All controllers share a 100 MHz master clock and trigger to support the pixel-synchronous readout required to eliminate patterns that would be caused by crosstalk from clocks on one controller to the video signal of another. The guide/focus CCDs cannot operate at the same speed, so one pixel is read for every three science pixels, to retain synchronization. True differential outputs of the science CCDs provide clock feedthrough attenuation and crosstalk suppression, which in combination with clock slew rate minimization allows pixel time to be reduced to 830 ns (Smith et al. 2018). Novel concurrent parallel clocking conceals line shift overhead so that readout time is only 8.2 sec, while delivering 10.3 electrons median read noise, well below the minimum shot noise in the sky (27 electrons).", '2.4. Optics': "Maintaining PTF's moderate image quality (2.0 '' FWHM in r -band) over the larger ZTF focal plane required novel optics. The focal surface of the Schmidt telescope is curved; the glass planes used in the photographic surveys conducted with P48 were physically bent on a mandrel to conform to this shape (Harrington 1952). For PTF, an optically powered dewar window was sufficient to provide good image quality  \nover the flat CCD array. However, this approach alone was insufficient to correct the much larger field of view of ZTF.  \nThe final ZTF optical design consists of four opticallypowered elements as well as the flat filters (§2.7). In front of the existing achromatic doublet Schmidt corrector that was installed for the Second Palomar Sky Surveys (Reid et al. 1991), we installed a 'trim plate' to modify the aspheric coefficient of the Schmidt corrector at the telescope pupil by about 10%. The trim plate was figured by the Nanjing Institute for Astronomical and Optical Technologies (NIAOT) in China, from a Corning fused silica blank. To handle the extreme field curvature of the Schmidt focus, the dewar vacuum window is a meniscus (with a conductive ITO coating on the inner surface providing resistive heating.) The CCDs themselves are mounted to a faceted cold plate, where each facet is a chord of the residual field curvature. Finally, to compensate for remaining curvature over each large science CCD, we mount ∼ 2 mm thick fused silica field flattener lenses 2 mmabove each detector (Figure 3).  \nUltimately, the useful field of the view of ZTF is limited by the Schmidt telescope design. At extreme field angles part of the beam from falls beyond the edge of the telescope primary mirror, with vignetting reaching 30% in the corners.", '2.5. Shutter': 'To minimize beam obscuration within the telescope tube, we placed the shutter at the entrance pupil of the telescope. This shutter was developed through a collaboration of Deutsches Elektronen-Synchrotron (DESY) and Bonn Shutter, who successfully delivered a bi-parting shutter with 1.2meter aperture that opens and closes within 430 msec while imparting less than 3 grams unbalanced force to the top of the telescope with negligible impact on image stability.', '2.6. Telescope': 'In order to optimize survey productivity as a system, and increase reliability, we invested in numerous upgrades to the Samuel Oschin Telescope. To minimize slew overheads, we upgraded both telescope drive axes as well as the dome drive system to enable faster slews. After upgrades, the (hour angle, declination, dome) drive accelerates to and decelerates from a top speed of (2.5, 3.0, 3.0) · s -1 at (0.4, 0.5, 0.5) · s -2 . With these upgrades, the telescope can slew and settle between adjacent fields, separated by 7 degrees, entirely during the CCD readout time.  \nOther upgrades (Dekany et al. 2019) included a new threevane instrument spider (to reduce beam obscuration), improved baffling of the telescope tube (to reduce scattered light), facility electrical improvements and lightning protection, a new dry air system (to inhibit condensation on window and refrigerant lines), refurbishment of the wind screen, and better thermal management in the dome.', '2.7. Filters and Filter Exchanger': 'ZTF has a complement of three custom filters, ZTFg , ZTFr , and ZTFi . Given the differences of the ZTF system relative to potential calibrators (SDSS, PS1, Gaia ), we did not attempt to match any existing filter bandpasses exactly. Instead, we sought to maximize the signal-to-noise ratio achieved by avoiding major lines in the Palomar sky background and to control costs on the large filters. Materion Precision Optics manufactured the g and r band filters and Asahi Spectra produced the i -band filter.  \nOur desire to minimize beam obstruction motivated an unusual design for the filter exchanger. We store all three filters in slots in a filter cabinet mounted in an access hatch of the telescope tube. A commercial robotic arm stows in a similar position. During the filter exchange, the arm uses a set of solenoid-deactivated magnets and redundant latches on its manipulator to dock with the frame holding the desired filter in the storage cabinet, move it to the camera, and secure it there. The arm then disconnects from the filter frame and stows against the wall of the telescope tube. At present, for safety, the exchanges are only conducted when the telescope is in a quasi-horizontal stow position where none of the moving elements physically pass over the telescope primary mirror. Including the slew time to and from the stow position, the additional overhead to change filters is ∼ 100 seconds in typical operations. Additional optimization of the arm motion profile and exchanging closer to zenith is expected to reduce this further.', '2.8. Robotic Observing System': 'The ZTF Robotic Observing Software (ROS) is based on the Robo-AO observing system (Riddle et al. 2012; Dekany et al. 2019) though many of the underlying tools were retained and upgraded. ROS is based on a modular, fail-safe, multi-threaded, multi-daemon software architecture. It has been designed to be able to run continuously for an extended period, while allowing human operators to monitor the system to determine its performance, track nightly errors, and reconfigure parameters if necessary. Configuration files support engineering and science operation modes. Extensive telemetry is aggregated from all ZTF hardware and telescope control subsystems.  \nROS is hosted on a single supervisory computer (which also controls the guide and focus CCDs) and four Archon camera control computers, each of which is responsible for readout of one quadrant of four science CCDs. Sufficient onsite data storage exists for at least two weeks of regular observing, in the unlikely event the microwave link from Palomar (§3.1) were to suffer an outage.', '2.9. Scheduler': 'The ZTF scheduler determines which fields to observe and in what order. Integer Linear Programming techniques inspired by Lampoudi et al. (2015) maximize the volumetric survey speed using slot-based lookahead throughout the night. Bellm et al. (2019) describes the scheduling system in detail.  \nDue to the desire to simplify the data processing for image subtraction (§3.4), all ZTF images are obtained on a fixed grid of fields with minimal dithering. The primary grid covers the entire sky with an average overlap between fields of about 0.26 · in Declination and about 0.29 · in Right Ascension. The fields are aligned to cover the Galactic Plane region with the fewest pointings, improving the efficiency of both Galactic and extragalactic surveys. We also took care to ensure that large local galaxies were placed effectively. A secondary grid of pointings, offset diagonally by about 5 · , fills in most of the CCD gaps and improves the fill factor within the survey footprint from 87.5% to 99.2% assuming no dithering.', '2.10. Flat Field Illuminator': 'PTF constructed its flat fields from science images taken each night. In addition to preventing final reduction of the images until the end of the night (Laher et al. 2014), this scheme was negatively affected by fringing of sky lines and scattered light from the moon and other bright sources and proved to be among the factors limiting PTF\'s photometric precision. For ZTF, a new Flat Field Illuminator system enables stable calibration frames to be taken before the night\'s observing.  \nThe Flat Field Illuminator consists of a reflective screen, LED illuminators, and a baffled enclosure. It is mounted on a tower close to the P48 dome. Twenty four narrowspectrum LEDs in each of 15 wavelengths spanning the ZTF filter bandpasses are spaced around a ring pointing towards a screen. The screen is constructed from aluminum honeycomb paneling which makes it lightweight, stiff, and flat. Many coats of Avian-D white paint provide a very uniform lambertian scattering surface. The forward baffle mounted on the telescope docks to a similar baffle surrounding the flat field system to fully enclose the path between flat field screen and detectors. The enclosure walls are heavily baffled and covered with 2% reflective Avian-D black paint facing the screen and black flocking facing the telescope. Similar baffles have been installed along the entire length of the enclosed telescope tube at sufficiently close spacing to block all single-bounce paths between flat field screen and primary mirror.  \nThis screen provides smooth and stable illumination for removing mid- to small-scale spatial frequencies in the sensitivity pattern. The 7% radial intensity variation at the screen integrates to <2% flat fielding error at the focal plane. This  \nresidual error occurs on large spatial scales that are easily corrected by calibrations derived from observing the relative photometry of stars as they are moved across the field (\'starflats\').  \nLEDs are driven by constant current sources, and their forward voltage is monitored to sense junction temperature so that temperature compensation can be applied if required. Flats are acquired separately in each LED wavelength and then combined with a relative weighting which minimizes the manifestation of CCD QE patterns in the "star-flats" which should only show mosaic-scale patterns.  \nThe principal error observed in flats is a 6% increase close to East and West edges of the CCDs where light scatters off the frames holding the field flattener lenses. This additive background must be removed from flats since it does not represent enhanced sensitivity. Fortunately it rises rapidly close to the edge of the CCD and can be fitted with sub-percent accuracy.', '3. DATA SYSTEM': 'The ZTF data processing system is housed at the Infrared Processing and Analysis Center (IPAC) and builds on the lessons learned processing data from PTF and iPTF (Ofek et al. 2012; Laher et al. 2014; Masci et al. 2017). Masci et al. (2019) provides a complete description of the ZTF pipelines.', '3.1. Data Transfer': 'The CCD controllers sample the video signal at 100 MHz and 16 bit resolution, averaging multiple samples to produce a floating point output with about 18 bits of dynamic range. We use the fpack program (Pence et al. 2010) to compress each quadrant and each overscan separately, allowing the compression to be optimized for the measured width of the core of the histogram in each image. In practice this noise root variance σ is dominated by sky noise ( ≥ 25 e -s -1 pixel -1 ). fpack converts the floating point data to integers applying a normalization factor q = 2, which results in σ = 2 for the histogram of integers. Lossless Rice compression is then applied. We apply a pseudo-random dither prior to normalization to avoid biases produced by rounding. The same dither values are subtracted after decompression (using the funpack program) so that the dither does not add noise. The result is that number of bits per pixel is reduced to an average of 5 during data transport at a cost of a 1% increase in sky noise, due to quantization by the normalization step. Despite the slight increase in noise, our tests confirm that this procedure does not appreciably bias image coaddition or photometry (cf. Price-Whelan & Hogg 2010; Pence et al. 2010, and references therein).  \nThe observing system transfers the images to IPAC via the NSF-funded High Performance Wireless Research and Education Network (HPWREN) administered by the University of California San Diego. Typical transfer times are  \n<25seconds, sufficient to keep up with the fastest observing cadence (38.3 seconds) throughout the night.', '3.2. Image Processing': 'Upon arrival, each multi-extension FITS image 4 is split into four readout quadrants per CCD and farmed out in parallel to the processing cluster. All subsequent processing is conducted independently on each CCD readout quadrant. Each image is tagged with the observing program that obtained it (public, collaboration, or Caltech), and the access permissions for all of the downstream data products are propagated accordingly.  \nThe image processing pipeline first subtracts bias frames and applies the flat field correction. The pipeline then calls the SCAMP package (Bertin 2006) to determine a World Coordinate System using Gaia DR1 (Gaia Collaboration et al. 2016) as the reference catalog. Subsequently the pipeline fits a zero point and color term to photometrically calibrate the quadrant to PanSTARRS 1 (Chambers et al. 2016). The pipeline sets appropriate mask bits for saturation, bad pixels, ghosts, and other instrumental artifacts.  \nThe pipeline produces both point-spread function (PSF) fit (DAOPHOT, Stetson 1987) and aperture (SExtractor, Bertin & Arnouts 1996a) photometry catalogs from the processed direct image, and the raw and processed images and catalogs are archived (§3.8).', '3.3. Reference Image Generation': 'Coadded reference images are required for image differencing (§3.4) as well as lightcurve source association (§3.7). Weconstruct reference images for each field, filter, and quadrant combination. Typical stacks have at least 15 images. We use Swarp (Bertin et al. 2002) to map the images to a common footprint and then compute an outlier-rejected average. Reference building pipelines are triggered automatically at the end of the night.', '3.4. Image Differencing': "The image differencing pipeline identifies moving and changing sources. It first prepares the processed science image and reference image by matching their photometric throughputs, warps the reference image onto the science image, and matches their backgrounds at low spatial frequencies. PSF-matching, image differencing, and the creation of an accompanying match-filtered image optimized for detecting point sources on the difference are then performed using the ZOGY algorithm (Zackay et al. 2016). The pipeline then detects both positive and negative 'candidate' sources at a signal-to-noise ratio greater than 5. The pipeline also measures a variety of pixel-based features for each candidate  \n(e.g., the number of positive and negative pixels in a region around the candidate) to provide to the Real-Bogus machine learning algorithm (Mahabal et al. 2019). Each candidate is loaded into a database and then packaged with other contextual information into an alert packet (§3.5) for distribution.  \nThe realtime pipeline runs from raw images to transient alerts in about four minutes.", '3.5. Transient Alert Stream': "The ZTF alert distribution system provides near-real-time access to transient and variable events identified by the image differencing pipelines. To aid the user in filtering the full alert stream for sources of interest, the ZTF alert stream provides rich alert packets containing not only the measurement that triggered the alert, but also a wide variety of contextual information. These include a Real-Bogus score (Mahabal et al. 2019) assessing the probability the candidate is astrophysical, a lightcurve of previous detections (or upper limits) from the last 30 days, a summary of prior detections before the 30-day window, cross-matches to the Pan-STARRS1 catalog along with a probabilistic star-galaxy score (Tachibana & Miller 2018), and FITS cutouts of the science, reference, and difference images.  \nThe alert packets themselves are serialized in the opensource Apache Avro format 5 . Schemas, example packets, and complete documentation of the packet fields are available 6 .  \nThe alert packets are distributed using the open-source queue system Apache Kafka 7 . Kafka provides a distributed queue that is scalable to the alert volumes expected by LSST. Patterson et al. (2019) describes the architecture and implementation of the alert distribution system more fully.  \nAlerts from ZTF's public survey stream in near-real time to community brokers, including the Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES; Narayan et al. 2018), ALeRCE 8 , Lasair, and Las Cumbres Observatory 9 which will provide public access. While the community brokers come online, we are also providing a bulk nightly release of public alerts 10 .", '3.6. Solar System Processing': 'Solar System Processing is divided between searches for streaked Near-Earth Objects and point-like moving objects. Both are detected in the difference image processing.  \nscience/  \nStreaked objects are identified by a dedicated pipeline originally developed for PTF (Waszczak et al. 2017).  \nPoint-like moving object candidates are identified at the end of the night by the ZTF Moving-Object Discovery Engine (ZMODE). ZMODE attempts to link tracklets from the last three observing nights and then fit orbits to them. Highquality objects are forwarded to human scanners for vetting and then reported to the Minor Planet Center.', '3.7. Direct Imaging Lightcurves': "For archival studies of variable stars and AGN in uncrowded fields, lightcurves built from direct (un-subtracted) images provide a higher-fidelity data product because they avoid the subtraction artifacts and additional noise produced by difference imaging. We build lightcurves every few months from the calibrated PSF photometry catalogs produced from the unsubtracted epochal images (§3.2). Starting from the catalogs built from the deep reference images, we associate the sources in each epochal PSF photometry catalog with the nearest source in the reference catalog. The resulting lightcurves are stored in HDF5 'matchfiles' on a field, quadrant basis 11 . To further improve the photometric precision, we solve for a small per-epoch shift in the absolute calibration zeropoint by minimizing the scatter of non-varying stars (Ofek et al. 2011), achieving better than 10 mmag repeatability for bright, unsaturated sources. Additionally, we store a variety of timeseries features (cf. Richards et al. 2011) computed on the lightcurve to aid in identification of variable sources.", '3.8. Archive and Data Releases': 'The Infrared Science Archive (IRSA) at IPAC provides archival access to ZTF images, catalogs, lightcurves, and archived alert packets 12 . Both interactive web-based and programmatic queries are supported. The first release of data products (other than alerts) from the public surveys is planned for one year after the start of the survey, in the second quarter of 2019 13 .', '4. ON-SKY PERFORMANCE': "ZTF achieved first light in October 2017. Commissioning activities continued through March 2018 and combined technical activities to verify the performance of the observing and data systems with science validation experiments.  \nFigure 5. Normalized histogram of point-source full-width at halfmaximum (FWHM) for all images in g (blue), r (orange), and i (red) bands during June 2018.  \n<!-- image -->  \nFormal survey operations began on March 20, 2018, although routine operations of the filter exchanger and guide and focus CCDs occurred only in April and June 2018 respectively.  \nFigure 5 shows the delivered image quality for all three ZTF filters. Median image quality for the subset of observations above airmass 1.2 was 2.1 '' FWHM in g -band, 2.0 '' FWHMin r -band, and 2.1 '' FWHMin i -band.  \nFigure 6 shows the limiting magnitudes obtained in all three filters over one lunation. Median five-sigma model limiting magnitudes are 20.8 mag in g -band, 20.6 mag in r -band, and 19.9 in i -band. Restricting to ± 3 days around new moon, the dark-time median limiting magnitudes are 21.1 mag in g -band, 20.9 mag in r -band, and 20.2 in i -band.", '5. SURVEY STRATEGY': "ZTF divides its observing time between three high-level programs: public surveys (40%), ZTF collaboration surveys (40%), and Caltech surveys (20%). Each program in turn divides its time between multiple sub-surveys. All of the available surveys are interleaved simultaneously by the survey scheduler (Bellm et al. 2019), which optimizes each night's schedule for volumetric survey speed while maintaining balance among the programs. Private surveys are not allowed to use the observation history of the public surveys in the scheduling process. In addition to performing the regularly scheduled surveys, the scheduler can perform Target of Opportunity (TOO) observations in response to external triggers. Each image is taken for one and only one owner in order to simplify access to derived data products (images, catalogs, lightcurve points, alerts). As the public surveys cover the entire available sky, some 'duplicate' observations are unavoidable. Here we give an overview of the public surveys; a detailed discussion of the surveys and on-sky scheduler per-  \nformance will appear in Bellm et al. (2019). Graham et al. (2019) provides an overview of some of the expected scientific returns.  \nDuring its public time, ZTF is conducting the two surveys of broad scientific utility that we proposed to the NSF MidScale Innovations Program (MSIP): a Northern Sky Survey and a Galactic Plane Survey. Motivated by the LSST baseline cadence (e.g., Ivezi'c et al. 2008), the Northern Sky Survey is a three-day cadence survey of all fields with centers north of δ = -31 · , except those in the Galactic Plane Survey 14 . The Galactic Plane Survey is a nightly survey of all visible fields in the region | b | < 7 · , δ > -31 · . For both surveys, each night a field is observed, it is visited twice, once in g -band and once in r -band, with at least 30 minutes separation between the two visits (cf. Miller et al. 2017). We expect to run these public surveys for at least the first eighteen months of the ZTF survey.  \nWe will attempt to obtain low-resolution spectra for all likely extragalactic transients brighter than 18.5 mag using the SED Machine (Blagorodnova et al. 2018) on the Palomar 60-inch and will publicly report these classifications (Fremling et al. 2018).", '6. FIRST RESULTS': "ZTF will enable new discoveries of many classes of astrophysical objects, including explosive extragalactic transients, optical counterparts of multi-wavelength and multimessenger phenomena, variable stars, Tidal Disruption Events, Active Galactic Nuclei, and solar system objects. Graham et al. (2019) presents ZTF's science goals in detail. In this section we present initial results in these areas from the early ZTF survey.", '6.1. Transient Science': "During commissioning of the alerts system, we searched the incoming alerts for astrophysical transients, both providing feedback for the machine learning by marking 'bogus' sources, and flagging potential supernovae for follow-up. Transient alerts were filtered and vetted via the GROWTH marshal system (Kasliwal et al. 2019) and using a machinelearning based classifier (Mahabal et al. 2019). In two months of commissioning data, we classified a total of 38 supernovae. Of these, 15 were only discovered by ZTF, while another 13 were first discovered by ZTF and later picked up by other surveys. The relatively modest yield is expected due to the limited set of reference images available, the need to maintain high thresholds while training the Real/Bogus system, and poor winter weather. All classified supernovae from commissioning data have been made public on the Transient  \nFigure 6. Left: Histogram of five-sigma limiting magnitudes in 30 second exposures for g (blue), r (orange), and i (red) bands over one lunation. Right: Limiting magnitudes for observations obtained within ± 3 days of new moon.  \n<!-- image -->  \nName Server (see Kulkarni 2018 and Lunnan et al. 2018 for details). The classification spectra are available on WISeREP (Yaron & Gal-Yam 2012).  \nAs an illustrative example from the commissioning data, we present ZTF18aaayemw (SN 2018yt), one of the first objects found, and one not identified by any other surveys. ZTF18aaayemw was discovered as a rising transient on 2018-02-07.26 (UT). Because this object was discovered so early in the survey, flux is also seen in the reference image that was built from data taken over the previous nights, so we cannot constrain the explosion date exactly. The light curve is shown in the left panel of Figure 7. An initial spectrum taken with the Nordic Optical Telescope on 2018-02-14 shows a featureless, blue continuum indicating a blackbody temperature of ∼ 12 , 000 K; narrow emission lines from the host galaxy sets the redshift at z = 0 . 0512. We continued to follow ZTF18aaayemw, and the sequence of spectra obtained is shown in the right panel of Figure 7. The spectrum remained blue and featureless for at least two weeks after discovery; spectra taken a month later show broad H α emission, classifying ZTF18aaayemw as a SN II. Details of the data collection and reduction are found in Appendix A.  \nThe early spectral evolution of ZTF18aaayemw is similar to that of other SNe II such as SN IIb iPTF13ast (Gal-Yam et al. 2014) and SN IIn iPTF11iqb (Smith et al. 2015), which also did not show broad features until later than 15 days postexplosion. These two supernovae also showed flash spectroscopy features (i.e., features from the stellar envelope or circumstellar material ionized by the supernova shock breakout), which we do not observe in ZTF18aaayemw. This could be because no such features were present, or because they have faded by our earliest spectrum at > 7 days. In the compilation of Khazov et al. (2016), only 1/13 SNe II where the first spectrum was taken 7-9 days after explosion showed  \nflash features, while 3/13 showed blue, featureless continua like we see in ZTF18aaayemw.", '6.2. Target of Opportunity and Multi-Messenger Science': "We performed Target-of-Opportunity (ToO) follow-up observations in response to IceCube-171106A (IceCube Collaboration 2017), a neutrino of likely astrophysical origin with an estimated energy in excess of 1 PeV. The neutrino was detected by the IceCube Neutrino Observatory, and was distributed as part of the IceCube Realtime Program (Aartsen et al. 2017). It was well-localized, with a sub-degree angular resolution, and was followed-up by ZTF with singlepointing observations. With ZTF's large field-of-view, such events will typically be covered by observations in a single field. Though the field was observed in this commissioning phase multiple times over a period of days, comparisons to reference images did not reveal any optical counterpart. Nonetheless, this example illustrates the potential of the ZTF ToO program to undertake multi-messenger observations of neutrino and gravitational-wave events.  \nZTF also observed the localization region of the short gamma-ray burst GRB180523A (trigger 548793993) detected by Fermi -GBM. ZTF obtained a series of r and g -band images covering 2900 square degrees beginning at 3:51 UT on 2018 May 24 (9.1 hours after the burst trigger time), corresponding to approximately 70% of the probability enclosed in the localization region. Images in r and g bands were again taken the following night. More than 100 highsignificance transient and variable candidates were identified by our pipeline in this area, all of which had previous detections with ZTF in the days and weeks prior to the GRB trigger time. No viable optical counterparts were thus identified. The median 5 sigma upper limit for an isolated point source in our images was r > 20 . 3 and g > 20 . 6 mag.  \n<!-- image -->  \nFigure 7. Left: Light curve of ZTF18aaayemw. The rise is well captured in ZTF data. 20 days after discovery, the supernova is still detectable in the UV. The spectroscopic epochs are marked along the bottom axis. Right: Sequence of spectra of ZTF18aaayemw. The spectrum stays featureless and blue for at least the first ∼ 20 days, before finally developing broad H α classifying ZTF18aaayemw as a SN II.  \n<!-- image -->", '6.3. Variable Science': "During commissioning we also validated ZTF's utility for studying variable stars using direct (non-difference) imaging.", '6.3.1. Variability of Be Stars': 'A fraction of Be stars are known to exhibit photometric variability due to the non-radial pulsation, ejected material, stellar winds, or instability of the decretion disk (see review in Rivinius et al. 2013, and references therein). A variety of kinds of variability with different time scales have been reported, including outbursts, long-term variation, and periodic variations (Okazaki 1997; Hubert & Floquet 1998; LabadieBartz et al. 2017). Using the ZTF commissioning data, we explored the variability on timescales of days to months of 83 Be star candidates in open clusters selected from Yu et al. (2018). In our preliminary examination of these data, we found that less than ∼ 10% of our Be star candidates show qualitative variability. Figure 8 gives one example of a Be star candidate that exhibits variability (upper panel) and another one that does not (lower panel). We expect that a longer time baseline as well as further refinements of the lightcurve pipeline will provide valuable constraints on the variability  \nof Be stars (such as variable fraction, amplitude of variation, outburst activity, and so on).', '6.3.2. RR Lyrae': 'The homogeneous gri -band light curves for RR Lyrae provided by the ZTF are also a useful tool to investigate their pulsational properties. For example, the period-color and amplitude-color relations of RR Lyrae can be used to probe the interaction of photosphere with the hydrogen ionization front in these type of pulsating stars (e.g., Ngeow et al. 2017, and references therein). To check the light curve quality for large-amplitude variable stars such as RR Lyrae, we constructed the light curves of known RR Lyrae in one ZTF field based on the ZTF commissioning data. Figure 9 shows the saw-tooth shape light curves for one bright and one faint RR Lyrae located in the selected ZTF field, demonstrating the expected light curve quality when ZTF is in full science operation. The finding of faint ( ∼ 20 . 5 mag), and hence distant, RR Lyrae will be useful for the study of the Galactic halo (e.g., Cohen et al. 2017, and references therein).', '6.4. Small Bodies in the Solar System': 'Figure 8. ZTF light curves of two Be star candidates selected from Yu et al. (2018) in g (filled blue triangles) and r (open red circles) bands. The magnitudes are based on the PSF photometry but have not had relative photometry corrections applied (see Masci et al. 2019), leading to larger observed scatter on a handful of nights.  \n<!-- image -->  \nFigure 9. ZTF r -band light curves for two known RR Lyrae based on PSF photometry but without relative photometric correction (see Masci et al. 2019). The pulsation periods P are taken from literature and not derived from the ZTF light curves.  \n<!-- image -->  \nSmall solar system bodies encompass comets and asteroids, Trojans, Centaurs, near-Earth objects (NEOs), and trans-Neptunian objects. ZTF will provide extensive observations of thousands of small bodies, allowing long-duration measurements of their positions, motions and brightnesses as a function of time. Below we discuss the potential science return from the ZTF observations of solar system objects, and highlight four examples of early results from the first months of operation.', '6.4.1. Near Earth Objects': 'The NEO search activity of ZTF comprises two components: detection of point-like NEOs, and detection of natural fast-moving objects that are moving more than a few degrees per day and hence appear as streaks. The ZTF Data System (Masci et al. 2019) scans all ZTF difference images for these two types of objects and releases candidate detections in near real-time. Screening of new detections and submission to the Minor Planet Center (MPC) has been done on a best effort basis since February 2018 for those fields for which good reference images exist. On a clear night with cadence and fields suitable for asteroid detection, ZTF can produce ∼ 100 , 000 detections of ∼ 25 , 000 asteroids.  \nBy 2018 May 4, after three months of operation, ZTF had submitted ∼ 600 , 000 measurements to the MPC and been assigned designations for about 320 new objects. The new discoveries include seven Near-Earth Asteroids (Table 2), of which one (2018 CL) is a Potentially Hazardous Asteroidan object with a minimum orbit intersection distance with Earth of less than 0.05 A.U. and MH < 22. Five of these seven new NEOs were detected by the dedicated streak-detection pipeline (Waszczak et al. 2017). Current efforts are aimed towards optimizing this pipeline for better rejection of false positives as we better characterize the new camera and detectors, and using citizen science through Zooniverse to increase the size of the training sample (for details, see Mahabal et al. 2019). Efficient new algorithms are also under development (Nir et al. 2018).', '6.4.2. Asteroid Light Curves': 'Asteroid light curves obtained from high-cadence observations can secure the measurements of their rotation periods and, moreover, facilitate the discovery of super-fast rotating asteroids (cf. Chang et al. 2017, and references therein). Wide-field facilities such as ZTF are particularly powerful for this type of science because of the efficiency of collecting numerous light curves within a short period of time (e.g., Masiero et al. 2009; Polishook & Brosch 2009; Dermawan et al. 2011; Polishook et al. 2012; Chang et al. 2014, 2015, 2016; Waszczak et al. 2015). To demonstrate the ability of the ZTF for this task, we conducted a pilot campaign on December 15, 2017, in which we repeatedly scanned between two ZTF fields on the ecliptic plane at opposition for ∼ 3', 'THE ZWICKY TRANSIENT FACILITY': 'Table 3. Log of ZTF18aaayemw Spectroscopic Observations', '6.4.3. Activity of Comets and Centaurs': "By covering the entire Northern sky approximately every three days (Sec. 5), ZTF acquires serendipitous observations of a large number of comets and centaurs. Through ZTF's high cadence and sensitivity, it is well suited to monitor the activity development of comets and to look for temporal variability, including both secular changes and rotational modulation of the activity, as well as transient events such as outbursts.  \nWeidentified comets and Centaurs in the ZTF data by comparing the telescope's observing logs to the ephemeris positions of all comets with predicted brightness V < 22 mag. This brightness limit is below ZTF's detection limit, but it is used not only because comet brightness predictions are notoriously poor, but also because an outburst could make a normally faint comet detectable. As of 24 April 2018, we estimate that ZTF had made 15000 observations of 186 comets brighter than 22 mag, and 3300 observations of 41 comets brighter than 18 mag.  \nZTF imaging of C/2016 R2 (PanSTARRS) acquired between November 11, 2017 and February 19, 2018 is presented in Fig. 12. The images show the comet before perihelion, approaching the Sun from 3.2 to 2.7 AU. At such  \nheliocentric distances, water sublimation rates are low, yet the comet had an impressive ion tail spanning over 0 . 5 · . This emission is fluorescence by CO + ions within the g band (Cochran & McKay 2018). No other volatiles have been detected and this comet appears to have an extremely high chemical abundance of CO (Cochran & McKay 2018), suggesting that CO sublimation drives the activity of this comet. Changes in the morphology of the ion tail reflect temporal variations in the comet's activity and in the local solar wind conditions (cf. Jones et al. 2018). ZTF monitoring will allow us to follow the comet's activity evolution until it falls below V > 21, anticipated around 8 AU from the Sun (JPL/Horizons).  \nThe first outburst observed by ZTF was seen when the Centaur Echeclus (q=5.8 AU, e=0.46, i=4 . 3 · ) exhibited an increase in activity at 7.3 AU from the Sun. The outburst, originally discovered by Brian Skiff at Lowell Observatory, occurred on 7 December 2017 UT, and was first observed by ZTF on 10 December. It produced a dust coma with peak Af ρ (a proxy for dust production; A'Hearn et al. 1984) of 20 , 000 ± 2500 cm (Fig. 11), similar to previously observed outbursts of this object (Bauer et al. 2008). Assuming a dust ejection velocity near 50 m/s for ∼ 1 µ m grain radii, we find a dust production rate ∼ 300 kg/s (cf. Bauer et al. 2008). The August/September 2016 outburst produced brightening that lasted just over a month, while the late 2017 outburst also lasted roughly 30 days, as shown in the ZTF data.", '7. SUMMARY': "ZTF will survey the Northern Hemisphere sky hundreds of times in three bands, with observations taken on timescales from minutes to years. We expect the resulting datasets to enable discovery of young supernovae and rare relativistic transients; construction of systematic samples of Tidal Disruption Events, Active Galactic Nuclei, and variable stars; and detailed measurements of a variety of solar system objects. Thanks to ZTF's moderate depth, ZTF discoveries will be readily amenable to follow-up observations with 1-5 m-class telescopes. Graham et al. (2019) provides a more thorough overview of the ZTF science case.", 'REFERENCES': "Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS,  \nAartsen, M. G., Ackermann, M., Adams, J., et al. 2017,  \nAstroparticle Physics, 92, 30  \nAfricano, B. M., Christensen, E. J., Fuls, D. C., et al. 2018, Minor  \nPlanet Electronic Circulars, 2018-H70  \nA'Hearn, M. F., Schleicher, D. G., Millis, R. L., Feldman, p. D., &  \nThompson, D. T. 1984, Astronomical Journal, 89, 579  \nAihara, H., Arimoto, N., Armstrong, R., et al. 2018, PASJ, 70, S4  \n219, 12 Astier, P., Guy, J., Regnault, N., et al. 2006, A&A, 447, 31 Bacci, P., Maestripieri, M., Tesi, L., et al. 2018, Minor Planet Electronic Circulars, 2018-C73 Bauer, J. M., Choi, Y.-J., Weissman, P. R., et al. 2008, PASP, 120, 393  \nBellm, E. C. 2016, PASP, 128, 084501  \nMinkowski, R. L., & Abell, G. O. 1963, The National Geographic Society-Palomar Observatory Sky Survey, ed. K. A. Strand (the University of Chicago Press), 481  \n- Miyazaki, S., Komiyama, Y., Kawanomoto, S., et al. 2018, PASJ, 70, S1  \nNarayan, G., Zaidi, T., Soraisam, M. D., et al. 2018, ApJS, 236, 9  \n- Ngeow, C.-C., Kanbur, S. M., Bhardwaj, A., Schrecengost, Z., & Singh, H. P. 2017, ApJ, 834, 160\n- Nir, G., Zackay, B., & Ofek, E. O. 2018, ArXiv e-prints, arXiv:1806.04204  \nOfek, E. O., Frail, D. A., Breslauer, B., et al. 2011, ApJ, 740, 65  \nOfek, E. O., Laher, R., Law, N., et al. 2012, PASP, 124, 62  \nOkazaki, A. T. 1997, A&A, 318, 548  \n- Oke, J. B., & Gunn, J. E. 1982, PASP, 94, 586\n- Patterson, M. T., Bellm, E. C., Rusholme, B., et al. 2019, PASP, 131, 018001  \nPence, W. D., White, R. L., & Seaman, R. 2010, PASP, 122, 1065 Pojmanski, G. 1997, AcA, 47, 467  \nPolishook, D., & Brosch, N. 2009, Icarus, 199, 319  \n- Polishook, D., Ofek, E. O., Waszczak, A., et al. 2012, MNRAS, 421, 2094  \nPrice-Whelan, A. M., & Hogg, D. W. 2010, PASP, 122, 207 Reid, I. N., Brewer, C., Brucato, R. J., et al. 1991, PASP, 103, 661 Richards, J. W., Starr, D. L., Butler, N. R., et al. 2011, ApJ, 733, 10 Riddle, R. L., Burse, M. P., Law, N. M., et al. 2012, in Proc. SPIE, Vol. 8447, Adaptive Optics Systems III, 84472O", 'A. OBSERVATIONS OF ZTF18AAAYEMW': 'A.1. Light curves  \nHost-subtracted PSF photometry of ZTF18aaayemw was produced from our P48 observations using a Pan-STARRS1 r -band stack as the reference image, since our P48 references contain SN light. The image subtraction and photometry methods follow Fremling et al. (2016).  \nIn addition to the P48 observations, we observed ZTF18aaayemw with the Centurion 28-inch telescope (C28) and 1-m telescope at the WISE observatory (Israel). In addition, we also obtained several epochs of UV photometry with Ultraviolet/Optical Telescope (UVOT; Roming et al. 2005) of the Neil Gehrels Swift Observatory (Gehrels et al. 2004). The data were reduced with routines in IRAF (Tody 1986) version 2.16. The world-coordinate system was calibrated with the software package astrometry.net (Lang et al. 2010) version 0.69. We measured the brightness of the transients using circular apertures in Source Extractor version 2.19.5 (Bertin & Arnouts 1996b) where the aperture diameter had a size of 1 × FWHMof the stellar point-spread function. The absolute flux calibration was secured with the Sloan Digital Sky Survey DR12 (Alam et al. 2015). The UVOT data were retrieved from the Swift Data Archive 15 . We used the standard UVOT data analysis software distributed with HEASOFT version 6.19, along with the standard calibration data. All photometry is summarized in Table 4, and has not been corrected for foreground or host reddening.', 'A.2. Spectroscopy': 'Spectra of ZTF18aaayemw were obtained with the AndaLucia Faint Object Spectrograph and Camera (ALFOSC) on the Nordic Optical Telescope, with the Double-Beamed Spectrograph (DBSP; Oke & Gunn 1982) on the 200-in Hale Telescope at Palomar Observatory, the Device Optimized for the LOw RESolution (DOLORES) on the Telescopio Nazionale Galileo, the DeVeny Spectrograph on the Discovery Channel Telescope, and the Dual Imager Spectrograph on the 3.5 ARC telescope at Apache Point Observatory. Details of the observations are listed in Table 3.  \nNOT and TNG spectra were reduced using a combination of IRAF and MATLAB scripts, which included bias and flat-field corrections; extraction of the 1D spectrum; wavelength calibration of the spectrum by comparison with the spectrum of an arc lamp; flux calibrations using the sensitivity function built from the spectra of a spectral standard star observed the same night. The TNG spectra from the two different grisms (see Table 3) were combined together.  \nThe APO+DIS spectrum was reduced using pydis 16 . A spectrophotometric standard star observed on the same night in the same instrumental configuration was used for flux calibration.  \nThe DCT DeVeny spectrum was reduced using standard IRAF routines. We first corrected for bias and flat-field then extracted the 1D spectrum. Wavelength and flux calibration were done by using a comparing with spectra of an arc lamp and the flux standard Feige34.  \nTable 4  \nBELLM ET AL. . ZTF18aaayemw Light Curve  \nTable 4 continued  \nTHE ZWICKY TRANSIENT FACILITY Table 4 (continued)'}