[en] 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² 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.

[en] The optical light curve of the quasar PG 1302-102 at $z=0.278$ shows a strong, smooth 5.2 year periodic signal, detectable over a period of ∼20 years. Although 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. At this close separation, the nuclear black holes in PG 1302-102 will likely merge within $\sim <!--\; ???\; -->{10}^5$ years due to gravitational wave emission alone. Here, we report the rest-frame near-infrared time lags for PG 1302-102. Compiling data from WISE and Akari, we confirm that the periodic behavior reported in the optical light curve from Graham et al. is reproduced at infrared wavelengths, with best-fit observed-frame 3.4 and $4.6\mathrm{\mu <!--\; ??\; -->}\mathrm{m}$ time lags of (2219 ± 153, 2408 ± 148) days for a near face-on orientation of the torus, or (4103 ± 153, 4292 ± 148) days for an inclined system with relativistic Doppler boosting in effect. The periodicity in the infrared light curves and the light-travel time of the accretion disk photons to reach the dust glowing regions support that a source within the accretion disk is responsible for the optical variability of PG 1302-102, echoed at the farther out dusty regions. The implied distance of this dusty, assumed toroidal region is ∼1.5 pc for a near face-on geometry or ∼1.1 pc for the relativistic Doppler-boosted case.

[en] Quasars have long been known as intrinsically variable sources, but the physical mechanism underlying the temporal optical/UV variability is still not well understood. We propose a novel nonparametric method for modeling and forecasting the optical variability of quasars utilizing an AE neural network to gain insight into the underlying processes. The AE is trained with ∼15,000 decade-long quasar light curves obtained by the Catalina Real-time Transient Survey selected with negligible flux contamination from the host galaxy. The AE’s performance in forecasting the temporal flux variation of quasars is superior to that of the damped random walk process. We find a temporal asymmetry in the optical variability and a novel relation—the amplitude of the variability asymmetry decreases as luminosity and/or black hole mass increases—is suggested with the help of autoencoded features. The characteristics of the variability asymmetry are in agreement with those from the self-organized disk instability model, which predicts that the magnitude of the variability asymmetry decreases as the ratio of the diffusion mass to inflow mass in the accretion disk increases.

[en] Presented here is 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 self-consistently 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.

[en] CRTS J084133.15+200525.8 is an optically bright quasar at z = 2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V ∼ 17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V ∼ 16.2. Only ∼1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H α in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km s⁻¹ in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.

[en] Using a filter in the GROWTH Marshal based on color and the amplitude and timescale of variability, we have identified 372 objects as known or candidate cataclysmic variables (CVs) during the second year of the operation of the Zwicky Transient Facility. From the available difference imaging data, we found that 93 are previously confirmed CVs and 279 are strong candidates. Spectra of four of the candidates confirm them as CVs by the presence of Balmer emission lines, while one of the four has prominent He ii lines indicative of containing a magnetic white dwarf. Gaia EDR3 parallaxes are available for 154 of these systems, resulting in distances from 108–2096 pc and absolute magnitudes in the range of 7.5–15.0, with the largest number of candidates between 10.5 and 12.5. The total numbers are 21% higher than from the previous year of the survey with a greater number of distances available but a smaller percentage of systems close to the Galactic plane. Comparison of these findings with a machine-learning method of searching all the light curves reveals large differences in each data set related to the parameters involved in the search process.

[en] In this work, we aimed to derive the gri-band period–luminosity (PL) and period–luminosity–color (PLC) relations for late-type contact binaries, for the first time, located in globular clusters, using the homogeneous light curves collected by the Zwicky Transient Factory (ZTF). We started with 79 contact binaries in 15 globular clusters, and retained 30 contact binaries in 10 globular clusters that have adequate numbers of data points in the ZTF light curves and are unaffected by blending. Magnitudes at mean and maximum light of these contact binaries were determined using a fourth-order Fourier expansion, while extinction corrections were done using the Bayerstar2019 3D reddening map together with adopting the homogeneous distances to their host globular clusters. After removing early-type and “anomaly” contact binaries, our derived gri-band PL and period–Wesenheit (PW) relations exhibited a much larger dispersion with large errors on the fitted coefficients. Nevertheless, the gr-band PL and PW relations based on this small sample of contact binaries in globular clusters were consistent with those based on a larger sample of nearby contact binaries. Good agreements of the PL and PW relations suggested both samples of contact binaries in the local Solar neighborhood and in the distant globular clusters can be combined and used to derive and calibrate the PL, PW, and PLC relations. The final derived gr-band PL, PW, and PLC relations were much improved over those based on the limited sample of contact binaries in the globular clusters.

[en] A supermassive black hole (SMBH) ejected from the potential well of its host galaxy via gravitational wave recoil carries important information about the mass ratio and spin alignment of the pre-merger SMBH binary. Such a recoiling SMBH may be detectable as an active galactic nucleus (AGN) broad-line region offset by up to 10 kpc from a disturbed host galaxy. We describe a novel methodology using forward modeling with The Tractor to search for such offset AGNs in a sample of 5493 optically variable AGNs detected with the Zwicky Transient Facility (ZTF). We present the discovery of nine AGNs that may be spatially offset from their host galaxies and are candidates for recoiling SMBHs. Five of these offset AGNs exhibit double-peaked broad Balmer lines, which may have arisen from unobscured accretion disk emission, and four show radio emission indicative of a relativistic jet. The fraction of double-peaked emitters in our spatially offset AGN sample is significantly larger than the 16% double-peaked emitter fraction observed for ZTF AGNs overall. In our sample of variable AGNs we also identified 52 merging galaxies, including a new spectroscopically confirmed dual AGN. Finally, we detected the dramatic rebrightening of SDSS 1133, a previously discovered variable object and recoiling SMBH candidate, in ZTF. The flare was accompanied by the reemergence of strong P Cygni line features, indicating that SDSS 1133 may be an outbursting luminous blue variable star.

[en] We present ZTF20aajnksq (AT 2020blt), a fast-fading (Δr = 2.3 mag in Δt = 1.3 days) red (g − r ≈ 0.6 mag) and luminous (M _{1626 Å} = −25.9 mag) optical transient at z = 2.9 discovered by the Zwicky Transient Facility (ZTF). AT 2020blt 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 _j = 1 d (observer frame); (2) a luminous (L _{0.3–10 KeV} = 10⁴⁶ erg s⁻¹) X-ray counterpart; and (3) luminous (L _{10 GHz} = 4 × 10³¹ erg s⁻¹ Hz⁻¹) radio emission. However, no GRB was detected in the 0.74 days between the last ZTF nondetection (r > 21.36 mag) and the first ZTF detection (r = 19.60 mag), with an upper limit on the isotropic-equivalent gamma-ray energy release of E _γ,iso < 7 × 10⁵² erg. AT 2020blt 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 AT 2020blt are consistent with a classical (initial Lorentz factor Γ₀ ≳ 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 AT 2020blt 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, such as dirty fireballs. We conclude by discussing the status and future of fast-transient searches in wide-field high-cadence optical surveys.

[en] The discovery rate of fast radio bursts (FRBs) is increasing dramatically thanks to new radio facilities. Meanwhile, wide-field instruments such as the 47 deg² Zwicky Transient Facility (ZTF) survey the optical sky to study transient and variable sources. We present serendipitous ZTF observations of the Canadian Hydrogen Intensity Mapping Experiment (CHIME) repeating source FRB 180916.J0158+65 that was localized to a spiral galaxy 149 Mpc away and is the first FRB suggesting periodic modulation in its activity. While 147 ZTF exposures corresponded to expected high-activity periods of this FRB, no single ZTF exposure was at the same time as a CHIME detection. No >3σ optical source was found at the FRB location in 683 ZTF exposures, totaling 5.69 hr of integration time. We combined ZTF upper limits and expected repetitions from FRB 180916.J0158+65 in a statistical framework using a Weibull distribution, agnostic of periodic modulation priors. The analysis yielded a constraint on the ratio between the optical and radio fluences of η ≲ 200, corresponding to an optical energy E _opt ≲ 3 × 10⁴⁶ erg for a fiducial 10 Jy ms FRB (90% confidence). A deeper (but less statistically robust) constraint of η ≲ 3 can be placed assuming a rate of $r(><!--\; >\; -->5\mathrm{J}\mathrm{y}\mathrm{m}\mathrm{s})=1{\mathrm{h}\mathrm{r}}^{-<!--\; ???\; -->1}$ and 1.2 ± 1.1 FRB occurring during exposures taken in high-activity windows. The constraint can be improved with shorter per-image exposures and longer integration time, or observing FRBs at higher Galactic latitudes. This work demonstrated how current surveys can statistically constrain multiwavelength counterparts to FRBs even without deliberately scheduled simultaneous radio observation.