We report on the results from the first six months of the Catalina Real-Time Transient Survey (CRTS). In order to search for optical transients (OTs) with timescales of minutes to years, the CRTS analyses data from the Catalina Sky Survey which repeatedly covers 26,000 of square degrees on the sky. The CRTS provides a public stream of transients that are bright enough to be followed up using small telescopes. Since the beginning of the survey, all CRTS transients have been made available to astronomers around the world in real time using HTML tables,RSS feeds, and VOEvents. As part of our public outreach program, the detections are now also available in Keyhole Markup Language through Google Sky. The initial discoveries include over 350 unique OTs rising more than 2 mag from past measurements. Sixty two of these are classified as supernovae (SNe), based on light curves, prior deep imaging and spectroscopic data. Seventy seven are due to cataclysmic variables (CVs; only 13 previously known), while an additional 100 transients were too infrequently sampled to distinguish between faint CVs and SNe. The remaining OTs include active galactic nucleus, blazars, high-proper-motions stars, highly variable stars (such as UV Ceti stars), and transients of an unknown nature. Our results suggest that there is a large population of SNe missed by many current SN surveys because of selection biases. These objects appear to be associated with faint host galaxies. We also discuss the unexpected discovery of white dwarf binary systems through dramatic eclipses.

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.

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;

With the discovery of the first transiting extrasolar planetary system back in 1999, a great number of projects started to hunt for other similar systems. Because the incidence rate of such systems was unknown and the length of the shallow transit events is only a few percent of the orbital period, the goal was to monitor continuously as many stars as possible for at least a period of a few months. Small aperture, large field of view automated telescope systems have been installed with a parallel development of new data reduction and analysis methods, leading to better than 1% per data point precision for thousands of stars. With the successful launch of the photometric satellites CoRoT and Kepler, the precision increased further by one-two orders of magnitude. Millions of stars have been analyzed and searched for transits. In the history of variable star astronomy this is the biggest undertaking so far, resulting in photometric time series inventories immensely valuable for the whole field. In this review we briefly discuss the methods of data analysis that were inspired by the main science driver of these surveys and highlight some of the most interesting variable star results that impact the field of variable star astronomy.

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.

Quasars have long been known to be variable sources at all wavelengths. Their optical variability is stochastic, can be due to a variety of physical mechanisms, and is well-described statistically in terms of a damped random walk model 1 . The recent availability of large collections of astronomical time series of flux measurements (light curves) 2-5 offers new data sets for a systematic exploration of quasar variability. Here we report on the detection of a strong, smooth periodic signal in the optical variability of the quasar PG 1302-102 with a mean observed period of 1,884 ± 88 days. It was identified in a search for periodic variability in a data set of light curves for 247,000 known, spectroscopically confirmed quasars with a temporal baseline of ∼ 9 years. While the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. Such systems are an expected consequence of galaxy mergers and can provide important constraints on models of galaxy formation and evolution.Subparsec supermassive black-hole (SMBH) binary systems are not resolvable except possibly with long baseline radio interferometry. An alternative approach to their detection is through a modulated variability -caused by, for example, perturbations in their accretion disks or precession of relativistic jets, if they are present (see Fig. 1). The best known candidate, OJ 287 6 , has shown a pair of outburst peaks every 12.2 years for at least the past century: this object can be interpreted as a secondary SMBH perturbing the accretion disk of the primary SMBH at regular intervals 7 . Systematic searches for equivalent systems to date 8-10 have attempted to identify them from broad-line velocity offsets in their optical and near-infrared spectra but cannot detect the closest pairs (with 0.1 pc separation).We applied a novel joint wavelet and autocorrelation function (ACF) based technique that identifies objects exhibiting strongly periodic behavior in their light curves (M.J.G. et al, manuscript in preparation) to the largest set of quasar time series currently available. These are drawn from the Catalina Real-time Transient Survey (CRTS) [11][12][13] . PG 1302-102 (see Fig. 2) is the strongest 1 arXiv:1501.01375v1 [astro-ph.GA] 7 Jan 2015 periodic candidate out of twenty objects meeting the selection criteria: strong constant wavelet peak, strong ACF detection of periodic behaviour, sufficient temporal coverage for 1.5 or greater cycles at the detected period, and a phased light curve well-described by a sinusoid. For statistical comparison, we have also generated a simulated light curve for each known quasar based on a damped random walk (DRW) model, a standard statistical description of the optical variability of quasars 1 , using the CRTS time sampling. We find that only one object from the simulated data of 247,000 quasars satisfies the same selection criteria, showing that the number of quasars selected is statistically significant and...

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 .