Search

Planetary Radio Interferometry and Doppler Experiment (PRIDE) is a multi-purpose experimental technique aimed at enhancing the science return of planetary missions. The technique exploits the science payload and spacecraft service systems without requiring a dedicated onboard instrumentation or imposing on the existing instrumentation any special for PRIDE requirements. PRIDE is based on the near-field phase-referencing Very Long Baseline Interferometry (VLBI) and evaluation of the Doppler shift of the radio signal transmitted by spacecraft by observing it with multiple Earth-based radio telescopes. The methodology of PRIDE has been developed initially at the Joint Institute for VLBI ERIC (JIVE) for tracking the ESA's Huygens Probe during its descent in the atmosphere of Titan in 2005. From that point on, the technique has been demonstrated for various planetary and other space science missions. The estimates of lateral position of the target spacecraft are done using the phase-referencing VLBI technique. Together with radial Doppler estimates, these observables can be used for a variety of applications, including improving the knowledge of the spacecraft state vector. The PRIDE measurements can be applied to a broad scope of research fields including studies of atmospheres through the use of radio occultations, the improvement of planetary and satellite ephemerides, as well as gravity field parameters and other geodetic properties of interest, and estimations of interplanetary plasma properties. This paper presents the implementation of PRIDE as a component of the ESA's Jupiter Icy Moons Explorer (JUICE) mission., Comment: Accepted for publication in Space Science Reviews, v. 219 (special issue on the ESA's JUICE mission)

SkyPortal is an open-source software package designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, a SkyPortal-based platform has been operating at scale for 2 yr for the Zwicky Transient Facility Phase II community, with hundreds of users, containing tens of millions of time-domain sources, interacting with dozens of telescopes, and enabling community reporting. While SkyPortal emphasizes rich user experiences (UX) across common frontend workflows, recognizing that scientific inquiry is increasingly performed programmatically, SkyPortal also surfaces an extensive and well-documented API system. From backend and frontend software to data science analysis tools and visualization frameworks, the SkyPortal design emphasizes the re-use and leveraging of best-in-class approaches, with a strong extensibility ethos. For instance, SkyPortal now leverages ChatGPT large-language models (LLMs) to automatically generate and surface source-level human-readable summaries. With the imminent re-start of the next-generation of gravitational wave detectors, SkyPortal now also includes dedicated multi-messenger features addressing the requirements of rapid multi-messenger follow-up: multi-telescope management, team/group organizing interfaces, and cross-matching of multi-messenger data streams with time-domain optical surveys, with interfaces sufficiently intuitive for the newcomers to the field. (abridged), Comment: Accepted to ApJS

Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to explosion. We performed a systematic survey of H-rich (Type II) SNe discovered within less than two days from explosion during the first phase of the Zwicky Transient Facility (ZTF) survey (2018-2020), finding thirty events for which a first spectrum was obtained within $< 2$ days from explosion. The measured fraction of events showing flash ionisation features ($>36\%$ at $95\%$ confidence level) confirms that elevated mass loss in massive stars prior to SN explosion is common. We find that SNe II showing flash ionisation features are not significantly brighter, nor bluer, nor more slowly rising than those without. This implies that CSM interaction does not contribute significantly to their early continuum emission, and that the CSM is likely optically thin. We measured the persistence duration of flash ionisation emission and find that most SNe show flash features for $\approx 5 $ days. Rarer events, with persistence timescales $>10$ days, are brighter and rise longer, suggesting these may be intermediate between regular SNe II and strongly-interacting SNe IIn.

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution., Comment: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022

Based on time-series observations collected from Zwicky Transient Facility (ZTF), we derived period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations for RR Lyrae located in globular clusters. We have applied various selection criteria to exclude RR Lyrae with problematic or spurious light curves. These selection criteria utilized information on the number of data points per light curve, amplitudes, colors, and residuals on the period-luminosity and/or period-Wesenheit relations. Due to blending, a number of RR Lyrae in globular clusters were found to be anomalously bright and have small amplitudes of their ZTF light curves. We used our final sample of ~750 RR Lyrae in 46 globular clusters covering a wide metallicity range (-2.36 < [Fe/H] < -0.54 dex) to derive PLZ and PWZ relations in gri bands. In addition, we have also derived the period-color-metallicity (PCZ) and for the first time, the PQZ relations where the Q-index is extinction-free by construction. We have compared our various relations to empirical and theoretical relations available in literature, and found a good agreement with most studies. Finally, we applied our derived PLZ relation to a dwarf galaxy, Crater II, and found its true distance modulus should be larger than the most recent determination., Comment: 28 pages, 7 Tables and 24 Figures. AJ accepted

There is a growing concern about an impact of low-Earth-orbit (LEO) satellite constellations on ground-based astronomical observations, in particular, on wide-field surveys in the optical and infrared. The Zwicky Transient Facility (ZTF), thanks to the large field of view of its camera, provides an ideal setup to study the effects of LEO megaconstellations - such as SpaceX's Starlink - on astronomical surveys. Here, we analyze the archival ZTF observations collected between 2019 November and 2021 September and find 5301 satellite streaks that can be attributed to Starlink satellites. We find that the number of affected images is increasing with time as SpaceX deploys more and more satellites. Twilight observations are particularly affected - a fraction of streaked images taken during twilight has increased from less than 0.5% in late 2019 to 18% in 2021 August. We estimate that once the size of the Starlink constellation reaches 10,000, essentially all ZTF images taken during twilight may be affected. However, despite the increase in satellite streaks observed during the analyzed period, the current science operations of ZTF are not yet strongly affected. We also find that redesigning Starlink satellites (by installing visors intended to block sunlight from reaching the satellite antennas to prevent reflection) reduces their brightness by a factor of 4.6 +/- 0.1 with respect to the original design in g, r, and i bands., Comment: accepted for publication in ApJ Letters, full table 1 included as ancillary file

Microlensing is a powerful technique to study the Galactic population of "dark" objects such as exoplanets both bound and unbound, brown dwarfs, low-luminosity stars, old white dwarfs, neutron stars, and almost the only way to study isolated stellar-mass black holes. The majority of previous efforts to search for gravitational microlensing events have concentrated towards high-density fields such as the Galactic bulge. Microlensing events in the Galactic plane have the advantage of closer proximity and better constrained relative proper motions, leading to better constrained lens mass estimates at the expense of a lower optical depth compared to events towards the Galactic bulge. We use the Zwicky Transient Facility (ZTF) Data Release 5 (DR5) compiled from 2018--2021 to survey the Galactic plane in the region of $|b| < 20^\circ$. We find a total of 60 candidate microlensing events including three that show a strong microlensing parallax effect. The rate of events traces Galactic structure, decreasing exponentially as a function Galactic longitude with scale length $\ell_0 \sim 37^\circ$. On average, we find Einstein timescales of our microlensing events to be about three times as long ($\sim60$ days) compared to those towards the Galactic bulge ($\sim20$ days). This pilot project demonstrates that microlensing towards the Galactic plane shows strong promise for characterization of dark objects within the Galatic disk., Comment: submitted to ApJ

The Zwicky Transient Facility (ZTF), a state-of-the-art optical robotic sky survey, registers on the order of a million transient events - such as supernova explosions, changes in brightness of variable sources, or moving object detections - every clear night, and generates associated real-time alerts. We present Alert-Classifying Artificial Intelligence (ACAI), an open-source deep-learning framework for the phenomenological classification of ZTF alerts. ACAI uses a set of five binary classifiers to characterize objects which, in combination with the auxiliary/contextual event information available from alert brokers, provides a powerful tool for alert stream filtering tailored to different science cases, including early identification of supernova-like and anomalous transient events. We report on the performance of ACAI during the first months of deployment in a production setting., Comment: Fourth Workshop on Machine Learning and the Physical Sciences (NeurIPS 2021)

We present observations of SN 2021csp, the second example of a newly-identified type of supernova (Type Icn) hallmarked by strong, narrow, P Cygni carbon features at early times. The SN appears as a fast and luminous blue transient at early times, reaching a peak absolute magnitude of -20 within 3 days due to strong interaction between fast SN ejecta (v ~ 30000 km/s) and a massive, dense, fast-moving C/O wind shed by the WC-like progenitor months before explosion. The narrow line features disappear from the spectrum 10-20 days after explosion and are replaced by a blue continuum dominated by broad Fe features, reminiscent of Type Ibn and IIn supernovae and indicative of weaker interaction with more extended H/He-poor material. The transient then abruptly fades ~60 days post-explosion when interaction ceases. Deep limits at later phases suggest minimal heavy-element nucleosynthesis, a low ejecta mass, or both, and imply an origin distinct from that of classical Type Ic supernovae. We place SN 2021csp in context with other fast-evolving interacting transients, and discuss various progenitor scenarios: an ultrastripped progenitor star, a pulsational pair-instability eruption, or a jet-driven fallback supernova from a Wolf-Rayet star. The fallback scenario would naturally explain the similarity between these events and radio-loud fast transients, and suggests a picture in which most stars massive enough to undergo a WR phase collapse directly to black holes at the end of their lives., Comment: Accepted to ApJ 29th December 2021. Originally submitted to ApJ on 6th August 2021

We present a software package for single-dish data processing of spacecraft signals observed with VLBI-equipped radio telescopes. The Spacecraft Doppler tracking (SDtracker) software allows one to obtain topocentric frequency detections with a sub-Hz precision, and reconstructed and residual phases of the carrier signal of any spacecraft or landing vehicle at any location in the Solar System. These data products are estimated using the ground-based telescope's highly stable oscillator as a reference, without requiring an a priori model of the spacecraft dynamics nor the downlink transmission carrier frequency. The software has been extensively validated in multiple observing campaigns of various deep space missions and is compatible with the raw sample data acquired by any standard VLBI radio telescope worldwide. In this paper, we report the numerical methodology of SDtracker, the technical operations for deployment and usage, and a summary of use cases and scientific results produced since its initial release., Comment: The paper has been accepted for publication in PASA

While it is difficult to observe the first black hole seeds in the early Universe, we can study intermediate mass black holes (IMBHs) in local dwarf galaxies for clues about their origins. In this paper we present a sample of variability--selected AGN in dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward--modeled mid--IR photometry of time--resolved Wide--field Infrared Survey Explorer ({\it WISE}) coadded images. We found that 44 out of 25,714 dwarf galaxies had optically variable AGN candidates, and 148 out of 79,879 dwarf galaxies had mid--IR variable AGN candidates, corresponding to active fractions of $0.17\pm0.03$\% and $0.19\pm0.02$\% respectively. We found that spectroscopic approaches to AGN identification would have missed 81\% of our ZTF IMBH candidates and 69\% of our {\it WISE} IMBH candidates. Only $9$ candidates have been detected previously in radio, X-ray, and variability searches for dwarf galaxy AGN. The ZTF and {\it WISE} dwarf galaxy AGN with broad Balmer lines have virial masses down to $10^{5.5}M_\odot$ and for the rest of the sample, BH masses predicted from host galaxy mass range between $10^{5.2}M_\odot

Type Ia supernovae (SNe~Ia) in the nearby Hubble flow are excellent distance indicators in cosmology. The Zwicky Transient Facility (ZTF) has observed a large sample of supernovae from an untargeted, rolling survey, reaching $20.8, 20.6, 20.3$ mag in $g$ $r$, and $i$-band, respectively. With a FoV of 47 sq.deg, ZTF discovered $>$ 3000 SNe~Ia in a little over 2.5 years. Here, we report on the sample of 761 spectroscopically classified SNe~Ia from the first year of operations (DR1). The sample has a median redshift $\bar z =$ 0.057, nearly a factor of two higher than the current low-$z$ sample. Our sample has a total of 934 spectra, of which 632 were obtained with the robotic SEDm on Palomar P60. We assess the potential for precision cosmology for a total of 305 SNe with redshifts from host galaxy spectra. The sample is already comparable in size to the entire combined literature low-$z$ anchor sample. The median first detection is 13.5 days before maximum light, about 10 days earlier than the median in the literature. Furthermore, six SNe from our sample are at $D_L < 80$ Mpc, for which host galaxy distances can be obtained in the JWST era, such that we have calibrator and Hubble flow SNe observed with the same instrument. In the entire duration of ZTF-I, we have observed nearly fifty SNe for which we can obtain calibrator distances, key for percent level distance scale measurements., Comment: 16 pages, 9 figures, accepted for publication in MNRAS. Data and associated information for the sample can be found at https://github.com/ZwickyTransientFacility/ztfcosmodr

We present the discovery of ZTF 21aaoryiz/SN 2021fcg -- an extremely low-luminosity Type Iax supernova. SN 2021fcg was discovered by the Zwicky Transient Facility in the star-forming galaxy IC0512 at a distance of $\approx$ 27 Mpc. It reached a peak absolute magnitude of $M_{r} =$ $-12.66\pm0.20$ mag, making it the least luminous thermonuclear supernova discovered to date. The E(B-V) contribution from the underlying host galaxy is unconstrained. However, even if it were as large as 0.5 mag, the peak absolute magnitude would be $M_{r} = -13.78\pm0.20$ mag -- still consistent with being the lowest luminosity SN. Optical spectra of SN 2021fcg taken at 37 and 65 days post maximum show strong [Ca II], Ca II and Na I D emission and several weak [Fe II] emission lines. The [Ca II] emission in the two spectra has extremely low velocities of $\approx 1300$ and $1000$ km s$^{-1}$ respectively. The spectra very closely resemble those of the very low luminosity Type Iax supernovae SN 2008ha, SN 2010ae and SN 2019gsc taken at similar phases. The peak bolometric luminosity of SN 2021fcg is $\approx$ $2.5^{+1.5}_{-0.3}\times10^{40}$ erg s$^{-1}$ which is a factor of three lower than that for SN 2008ha. The bolometric lightcurve of SN 2021fcg is consistent with a very low ejected nickel mass (M$_{\rm{Ni}} \approx 0.8^{+0.4}_{-0.5}\times10^{-3}$ M$_{\odot}$). The low luminosity and nickel mass of SN 2021fcg pose a challenge to the picture that low luminosity SNe Iax originate from deflagrations of near M$_{\rm{ch}}$ hybrid carbon-oxygen-neon white dwarfs. Instead, the merger of a carbon-oxygen and oxygen-neon white dwarf is a promising model to explain SN 2021fcg., Comment: accepted for publication in ApJL. 10 pages, 5 figures

White dwarfs represent the last stage of evolution of stars with mass less than about eight times that of the Sun and, like other stars, are often found in binaries. If the orbital period of the binary is short enough, energy losses from gravitational-wave radiation can shrink the orbit until the two white dwarfs come into contact and merge. Depending on the component masses, the merger can lead to a supernova of type Ia or result in a massive white dwarf. In the latter case, the white dwarf remnant is expected to be highly magnetised because of the strong magnetic dynamo that should arise during the merger, and be rapidly spinning from the conservation of the orbital angular momentum. Here we report observations of a white dwarf, ZTF J190132.9+145808.7, that exhibits these properties, but to an extreme: a rotation period of 6.94 minutes, a magnetic field ranging between 600 megagauss and 900 megagauss over its surface, and a stellar radius of about 2,100 km, slightly larger than the radius of the Moon. Such a small radius implies that the star's mass is close to the maximum white-dwarf mass, or Chandrasekhar mass. ZTF J190132.9+145808.7 is likely to be cooling through the Urca processes (neutrino emission from electron capture on sodium) because of the high densities reached in its core., Comment: 37 pages, 9 figures, published in Nature on June 30, 2021

Using a filter in the GROWTH Marshal based on color and the amplitude and the timescale of variability, we have identified 372 objects as known or candidate cataclysmic variables (CVs) during the second year of operation of the Zwicky Transient Facility (ZTF). 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 HeII 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-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 dataset related to the parameters involved in the search process., Comment: Accepted in AJ, 24 pages, 2 tables, 7 figures

We present results from the Large Adaptive optics Survey for Substellar Objects (LASSO), where the goal is to directly image new substellar companions (<70 M$_{Jup}$) at wide orbital separations ($\gtrsim$50 AU) around young ($\lesssim$300 Myrs), nearby (<100 pc), low-mass ($\approx$0.1-0.8 M$_{\odot}$) stars. We report on 427 young stars imaged in the visible (i') and near-infrared (J or H) simultaneously with Robo-AO on the Kitt Peak 2.1-m telescope and later the Maunakea University of Hawaii 2.2-m telescope. To undertake the observations, we commissioned a new infrared camera for Robo-AO that uses a low-noise high-speed SAPHIRA avalanche photodiode detector. We detected 121 companion candidates around 111 stars, of which 62 companions are physically associated based on Gaia DR2 parallaxes and proper motions, another 45 require follow-up observations to confirm physical association, and 14 are background objects. The companion separations range from 2-1101 AU and reach contrast ratios of 7.7 magnitudes in the near infrared compared to the primary. The majority of confirmed and pending candidates are stellar companions, with ~5 being potentially substellar and requiring follow-up observations for confirmation. We also detected a 43$\pm$9 M$_{Jup}$ and an 81$\pm$5 M$_{Jup}$ companion that were previously reported. We found 34 of our targets have acceleration measurements detected using Hipparcos-Gaia proper motions. Of those, 58$^{+12}_{-14}$% of the 12 stars with imaged companion candidates have significant accelerations ($\chi^2 >11.8$), while only 23$^{+11}_{-6}$% of the remaining 22 stars with no detected companion have significant accelerations. The significance of the acceleration decreases with increasing companion separation. These young accelerating low-mass stars with companions will eventually yield dynamical masses with future orbit monitoring., Comment: 44 pages: 20 pages of text with 14 figures, 24 pages of tables in the appendix. Accepted for publication in the Astronomical Journal

We present a search for extragalactic fast blue optical transients (FBOTs) during Phase I of the Zwicky Transient Facility (ZTF). We identify 38 candidates with durations above half-maximum light 1 d < t1/2 < 12 d, of which 28 have blue (g-r<-0.2 mag) colors at peak light. Of the 38 transients (28 FBOTs), 19 (13) can be spectroscopically classified as core-collapse supernovae (SNe): 11 (8) H- or He-rich (Type II/IIb/Ib) SNe, 6 (4) interacting (Type IIn/Ibn) SNe, and 2 (1) H&He-poor (Type Ic/Ic-BL) SNe. Two FBOTs (published previously) had high-S/N predominantly featureless spectra and luminous radio emission: AT2018lug and AT2020xnd. Seven (five) did not have a definitive classification: AT 2020bdh showed tentative broad H$\alpha$ in emission, and AT 2020bot showed unidentified broad features and was 10 kpc offset from the center of an early-type galaxy. Ten (six) have no spectroscopic observations or redshift measurements. We present multiwavelength (radio, millimeter, and/or X-ray) observations for five FBOTs (three Type Ibn, one Type IIn/Ibn, one Type IIb). Additionally, we search radio-survey (VLA and ASKAP) data to set limits on the presence of radio emission for 22 of the transients. All X-ray and radio observations resulted in non-detections; we rule out AT2018cow-like X-ray and radio behavior for five FBOTs and more luminous emission (such as that seen in the Camel) for four additional FBOTs. We conclude that exotic transients similar to AT2018cow, the Koala, and the Camel represent a rare subset of FBOTs, and use ZTF's SN classification experiments to measure the rate to be at most 0.1% of the local core-collapse SN rate., Comment: Replaced following peer-review process. 46 pages, 20 figures. Accepted for publication in ApJ

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 the 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 number of data points in the ZTF light curves and 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 {\tt 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 exhibit 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 than those based on the limited sample of contact binaries in the globular clusters., Comment: 26 pages, 4 Tables and 16 Figures; AJ Accepted

Cometary activity is a manifestation of sublimation-driven processes at the surface of nuclei. However, cometary outbursts may arise from other processes that are not necessarily driven by volatiles. In order to fully understand nuclear surfaces and their evolution, we must identify the causes of cometary outbursts. In that context, we present a study of mini-outbursts of comet 46P/Wirtanen. Six events are found in our long-term lightcurve of the comet around its perihelion passage in 2018. The apparent strengths range from $-0.2$ to $-1.6$ mag in a 5" radius aperture, and correspond to dust masses between $\sim10^4$ to $10^6$ kg, but with large uncertainties due to the unknown grain size distributions. However, the nominal mass estimates are the same order of magnitude as the mini-outbursts at comet 9P/Tempel 1 and 67P/Churyumov-Gerasimenko, events which were notably lacking at comet 103P/Hartley 2. We compare the frequency of outbursts at the four comets, and suggest that the surface of 46P has large-scale ($\sim$10-100 m) roughness that is intermediate to that of 67P and 103P, if not similar to the latter. The strength of the outbursts appear to be correlated with time since the last event, but a physical interpretation with respect to solar insolation is lacking. We also examine Hubble Space Telescope images taken about 2 days following a near-perihelion outburst. No evidence for macroscopic ejecta was found in the image, with a limiting radius of about 2-m., Comment: Accepted for publication in The Planetary Science Journal. 33 pages, 11 figures, 3 tables

Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsars). The discovery of SN 1998bw/GRB 980425 marked the first association of a LGRB with a collapsar and AT 2017gfo/GRB 170817A/GW170817 marked the first association of a SGRB with a binary neutron star merger, producing also gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with being the supernova (SN). Despite the GRB duration being short (rest-frame T90 of 0.65 s), our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets., Comment: submitted to Nature Astronomy

We present the goals, strategy and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short period binaries and stellar pulsators with periods less than a few hours. Between June 2018 and January 2019, we observed 64 ZTF fields resulting in 2990 deg$^2$ of high stellar density in ZTF-$r$ band along the Galactic Plane. Each field was observed continuously for 1.5 to 6 hrs with a cadence of 40 sec. Most fields have between 200 and 400 observations obtained over 2-3 continuous nights. As part of this survey we extract a total of $\approx$230 million individual objects with at least 80 epochs obtained during the high-cadence Galactic Plane survey reaching an average depth of ZTF-$r$ $\approx$20.5 mag. For four selected fields with 2 million to 10 million individual objects per field we calculate different variability statistics and find that $\approx$1-2% of the objects are astrophysically variable over the observed period. We present a progress report on recent discoveries, including a new class of compact pulsators, the first members of a new class of Roche Lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring stars. Finally we present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between ZTF-$r$ = 20-76 mmag and pulsation periods between $P$ = 5.8-16 min with a strong cluster of systems with periods $\approx$ 6 min. All of the data have now been released in either ZTF Data Release 3 or data release 4., Comment: accepted for publication in MNRAS, 16 pages, 10 figures and 5 tables

AM CVn systems are a rare type of accreting binary that consists of a white dwarf and a helium-rich, degenerate donor star. Using the Zwicky Transient Facility (ZTF), we searched for new AM CVn systems by focusing on blue, outbursting stars. We first selected outbursting stars using the ZTF alerts. We cross-matched the candidates with $Gaia$ and Pan-STARRS catalogs. The initial selection of candidates based on the $Gaia$ $BP$-$RP$ contains 1751 unknown objects. We used the Pan-STARRS $g$-$r$ and $r$-$i$ color in combination with the $Gaia$ color to identify 59 high-priority candidates. We obtained identification spectra of 35 sources, of which 18 are high priority candidates, and discovered 9 new AM CVn systems and one magnetic CV which shows only He-II lines. Using the outburst recurrence time, we estimate the orbital periods which are in the range of 29 to 50 minutes. We conclude that targeted followup of blue, outbursting sources is an efficient method to find new AM CVn systems, and we plan to followup all candidates we identified to systematically study the population of outbursting AM CVn systems., Comment: Submitted, comments and feedback welcome

We present SNIascore, a deep-learning based method for spectroscopic classification of thermonuclear supernovae (SNe Ia) based on very low-resolution (R $\sim100$) data. The goal of SNIascore is fully automated classification of SNe Ia with a very low false-positive rate (FPR) so that human intervention can be greatly reduced in large-scale SN classification efforts, such as that undertaken by the public Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS). We utilize a recurrent neural network (RNN) architecture with a combination of bidirectional long short-term memory and gated recurrent unit layers. SNIascore achieves a $<0.6\%$ FPR while classifying up to $90\%$ of the low-resolution SN Ia spectra obtained by the BTS. SNIascore simultaneously performs binary classification and predicts the redshifts of secure SNe Ia via regression (with a typical uncertainty of $<0.005$ in the range from $z = 0.01$ to $z = 0.12$). For the magnitude-limited ZTF BTS survey ($\approx70\%$ SNe Ia), deploying SNIascore reduces the amount of spectra in need of human classification or confirmation by $\approx60\%$. Furthermore, SNIascore allows SN Ia classifications to be automatically announced in real-time to the public immediately following a finished observation during the night., Comment: 12 pages, 5 figures, 2 tables, accepted for publication in ApJL

While optical surveys regularly discover slow transients like supernovae on their own, the most common way to discover extragalactic fast transients, fading away in a few nights, is via follow-up observations of gamma-ray burst and gravitational-wave triggers. However, wide-field surveys have the potential to also identify rapidly fading transients independently of such external triggers. The volumetric survey speed of the Zwicky Transient Facility (ZTF) makes it sensitive to faint and fast-fading objects as kilonovae, the optical counterparts to binary neutron stars and neutron star-black hole mergers, out to almost 200Mpc. We introduce an open-source software infrastructure, the ZTF REaltime Search and Triggering, ZTFReST, designed to identify kilonovae and fast optical transients in ZTF data. Using the ZTF alert stream combined with forced photometry, we have implemented automated candidate ranking based on their photometric evolution and fitting to kilonova models. Automated triggering of follow-up systems, such as Las Cumbres Observatory, has also been implemented. In 13 months of science validation, we found several extragalactic fast transients independent of any external trigger (though some counterparts were identified later), including at least one supernova with post-shock cooling emission, two known afterglows with an associated gamma-ray burst, two known afterglows without any known gamma-ray counterpart, and three new fast-declining sources (ZTF20abtxwfx, ZTF20acozryr, and ZTF21aagwbjr) that are likely associated with GRB200817A, GRB201103B, and GRB210204A. However, we have not found any objects which appear to be kilonovae; therefore, we constrain the rate of GW170817-like kilonovae to $R < 900$Gpc$^{-3}$yr$^{-1}$. A framework such as ZTFReST could become a prime tool for kilonova and fast transient discovery with the Vera C. Rubin Observatory., Comment: Submitted for publication

Using Zwicky Transient Facility (ZTF) observations, we identify a pair of "sibling" Type Ia supernovae (SNe Ia), i.e., hosted by the same galaxy at z = 0.0541. They exploded within 200 days from each other at a separation of $0.6^{"} $ corresponding to a projected distance of only 0.6 kpc. Performing SALT2 light curve fits to the gri ZTF photometry, we show that for these equally distant "standardizable candles", there is a difference of 2 magnitudes in their rest frame B-band peaks, and the fainter SN has a significantly red SALT2 colour $c = 0.57 \pm$ 0.04, while the stretch values $x_1$ of the two SNe are similar, suggesting that the fainter SN is attenuated by dust in the interstellar medium of the host galaxy. We use these measurements to infer the SALT2 colour standardization parameter, $\beta$ = 3.5 $\pm$ 0.3, independent of the underlying cosmology and Malmquist bias. Assuming the colour excess is entirely due to dust, the result differs by $2\sigma$ from the average Milky-Way total-to-selective extinction ratio, but is in good agreement with the colour-brightness corrections empirically derived from the most recent SN Ia Hubble-Lemaitre diagram fits. Thus we suggest that SN "siblings", which will increasingly be discovered in the coming years, can be used to probe the validity of the colour and lightcurve shape corrections using in SN Ia cosmology while avoiding important systematic effects in their inference from global multi-parameter fits to inhomogeneous data-sets, and also help constrain the role of interstellar dust in SN Ia cosmology., Comment: Replaced with version accepted in MNRAS

The many unusual properties of the enigmatic AT2018cow suggested that at least some subset of the empirical class of fast blue optical transients (FBOTs) represents a genuinely new astrophysical phenomenon. Unfortunately, the intrinsic rarity and fleeting nature of these events have made it difficult to identify additional examples early enough to acquire the observations necessary to constrain theoretical models. We present here the Zwicky Transient Facility discovery of AT2020xnd (ZTF20acigmel, the "Camel") at z=0.243, the first unambiguous AT2018cow analog to be found and confirmed in real time. AT2018cow and AT2020xnd share all key observational properties: a fast optical rise, sustained high photospheric temperature, absence of a second peak attributable to ejection of a radioactively-heated stellar envelope, extremely luminous radio, millimetre, and X-ray emission, and a dwarf-galaxy host. This supports the argument that AT2018cow-like events represent a distinct phenomenon from slower-evolving radio-quiet supernovae, likely requiring a different progenitor or a different central engine. The sample properties of the four known members of this class to date disfavour tidal disruption models but are consistent with the alternative model of an accretion powered jet following the direct collapse of a massive star to a black hole. Contextual filtering of alert streams combined with rapid photometric verification using multi-band imaging provides an efficient way to identify future members of this class, even at high redshift., Comment: Accepted to MNRAS

We present Tails, an open-source deep-learning framework for the identification and localization of comets in the image data of the Zwicky Transient Facility (ZTF), a robotic optical time-domain survey currently in operation at the Palomar Observatory in California, USA. Tails employs a custom EfficientDet-based architecture and is capable of finding comets in single images in near real time, rather than requiring multiple epochs as with traditional methods. The system achieves state-of-the-art performance with 99% recall, 0.01% false positive rate, and 1-2 pixel root mean square error in the predicted position. We report the initial results of the Tails efficiency evaluation in a production setting on the data of the ZTF Twilight survey, including the first AI-assisted discovery of a comet (C/2020 T2) and the recovery of a comet (P/2016 J3 = P/2021 A3).

We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hours on multiple days in 2020 August, September and October. The photometric lightcurves have a repeating $\lesssim$0.05 magnitude feature suggesting that (6478) Gault has a rotation period of ~2.5 hours and may have a semi-spherical or top-like shape, much like Near-Earth Asteroids Ryugu and Bennu. The rotation period of ~2.5 hours is near to the expected critical rotation period for an asteroid with the physical properties of (6478) Gault suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spun up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect., Comment: 23 pages, 18 figures; Accepted by ApJ Letters

The Zwicky Transient Facility (ZTF) has been observing the entire northern sky since the start of 2018 down to a magnitude of 20.5 ($5 \sigma$ for 30s exposure) in $g$, $r$, and $i$ filters. Over the course of two years, ZTF has obtained light curves of more than a billion sources, each with 50-1000 epochs per light curve in $g$ and $r$, and fewer in $i$. To be able to use the information contained in the light curves of variable sources for new scientific discoveries, an efficient and flexible framework is needed to classify them. In this paper, we introduce the methods and infrastructure which will be used to classify all ZTF light curves. Our approach aims to be flexible and modular and allows the use of a dynamical classification scheme and labels, continuously evolving training sets, and the use of different machine learning classifier types and architectures. With this setup, we are able to continuously update and improve the classification of ZTF light curves as new data becomes available, training samples are updated, and new classes need to be incorporated.

The Zwicky Transient Facility is a time-domain optical survey that has substantially increased our ability to observe and construct massive catalogs of astronomical objects by use of its 47 square degree camera that can observe in multiple filters. However the telescope's i-band filter suffers from significant atmospheric fringes that reduce photometric precision, especially for faint sources and in multi-epoch co-additions. Here we present a method for constructing models of these atmospheric fringes using Principal Component Analysis that can be used to identify and remove these artifacts from contaminated images. In addition, we present the Uniform Background Indicator as a quantitative measurement of the reduced correlated background noise and photometric error present after removing fringes. We conclude by evaluating the effect of our method on measuring faint sources through the injection and recovery of artificial stars in both single-image epochs and co-additions. Our method for constructing atmospheric fringe models and applying those models to produce cleaned images is available for public download in the open source python package \href{https://github.com/MichaelMedford/fringez}{fringez}., Comment: 13 pages, 10 figures, accepted to PASP (6/24/2021)

HO Puppis (HO Pup) was considered as a Be-star candidate based on its gamma-Cassiopeiae-type light curve, but lacked spectroscopic confirmation. Using distance measured from Gaia Data Release 2 and the spectral-energy-distribution (SED) fit on broadband photometry, the Be-star nature of HO Pup is ruled out. Furthermore, based on the 28,700 photometric data points collected from various time-domain surveys and dedicated intensive-monitoring observations, the light curves of HO Pup closely resemble IW And-type stars (as pointed out in Kimura et al. 2020a), exhibiting characteristics such as quasi-standstill phase, brightening, and dips. The light curve of HO Pup displays various variability timescales, including brightening cycles ranging from 23 to 61 days, variations with periods between 3.9 days and 50 minutes during the quasi-standstill phase, and a semi-regular ~14-day period for the dip events. We have also collected time-series spectra (with various spectral resolutions), in which Balmer emission lines and other expected spectral lines for an IW And-type star were detected (even though some of these lines were also expected to be present for Be stars). We detect Bowen fluorescence near the brightening phase, and that can be used to discriminate between IW And-type stars and Be stars. Finally, despite only observing for four nights, the polarization variation was detected, indicating that HO Pup has significant intrinsic polarization., Comment: 24 pages, 7 tables and 12 figures; ApJ accepted

AT2019wey (SRGA J043520.9+552226, SRGE J043523.3+552234) is a transient first reported by the ATLAS optical survey in 2019 December. It rose to prominence upon detection, three months later, by the Spektrum-Roentgen-Gamma (SRG) mission in its first all-sky survey. X-ray observations reported in Yao et al. suggest that AT2019wey is a Galactic low-mass X-ray binary (LMXB) with a black hole (BH) or neutron star (NS) accretor. Here we present ultraviolet, optical, near-infrared, and radio observations of this object. We show that the companion is a short-period (P < 16 hr) low-mass (< 1 Msun) star. We consider AT2019wey to be a candidate BH system since its locations on the L_radio--L_X and L_opt--L_X diagrams are closer to BH binaries than NS binaries. We demonstrate that from 2020 June to August, despite the more than 10 times brightening at radio and X-ray wavelengths, the optical luminosity of AT2019wey only increased by 1.3--1.4 times. We interpret the UV/optical emission before the brightening as thermal emission from a truncated disk in a hot accretion flow and the UV/optical emission after the brightening as reprocessing of the X-ray emission in the outer accretion disk. AT2019wey demonstrates that combining current wide-field optical surveys and SRG provides a way to discover the emerging population of short-period BH LMXB systems with faint X-ray outbursts., Comment: Accepted by ApJ

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 \texttt{The Tractor} to search for such offset AGN in a sample of 5493 optically variable AGN detected with the Zwicky Transient Facility (ZTF). We present the discovery of 9 AGN which may be spatially offset from their host galaxies and are candidates for recoiling SMBHs. Five of these offset AGN exhibit double-peaked broad Balmer lines which may arise 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 AGN overall. In our sample of variable AGN we also identified 52 merging galaxies, including a new spectroscopically confirmed dual AGN. Finally, we detected the dramatic rebrightening of SDSS1133, a previously discovered variable object and recoiling SMBH candidate, in ZTF. The flare was accompanied by the re-emergence of strong P-Cygni line features indicating that it may be an outbursting luminous blue variable star., Comment: Accepted by ApJ. 31 pages, 16 figures

We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$\mu$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1'' jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$\rho$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$\rho$/cross-section remained constant, it implies LD$_2$'s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $\mu$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years., Comment: 28 pages, 8 figures, 2 Tables, accepted for publication in AJ

Interaction-powered supernovae (SNe) explode within an optically-thick circumstellar medium (CSM) that could be ejected during eruptive events. To identify and characterize such pre-explosion outbursts we produce forced-photometry light curves for 196 interacting SNe, mostly of Type IIn, detected by the Zwicky Transient Facility between early 2018 and June 2020. Extensive tests demonstrate that we only expect a few false detections among the 70,000 analyzed pre-explosion images after applying quality cuts and bias corrections. We detect precursor eruptions prior to 18 Type IIn SNe and prior to the Type Ibn SN2019uo. Precursors become brighter and more frequent in the last months before the SN and month-long outbursts brighter than magnitude -13 occur prior to 25% (5 - 69%, 95% confidence range) of all Type IIn SNe within the final three months before the explosion. With radiative energies of up to $10^{49}\,\text{erg}$, precursors could eject $\sim1\,\text{M}_\odot$ of material. Nevertheless, SNe with detected precursors are not significantly more luminous than other SNe IIn and the characteristic narrow hydrogen lines in their spectra typically originate from earlier, undetected mass-loss events. The long precursor durations require ongoing energy injection and they could, for example, be powered by interaction or by a continuum-driven wind. Instabilities during the neon and oxygen burning phases are predicted to launch precursors in the final years to months before the explosion; however, the brightest precursor is 100 times more energetic than anticipated.

The Zwicky Transient Facility (ZTF) has discovered five new events belonging to an emerging class of AGN undergoing smooth flares with large amplitudes and rapid rise times. This sample consists of several transients that were initially classified as supernovae with narrow spectral lines. However, upon closer inspection, all of the host galaxies display resolved Balmer lines characteristic of a narrow-line Seyfert 1 (NLSy1) galaxy. The transient events are long-lived, over 400 days on average. We report UV and X-ray follow-up of the flares and observe persistent UV-bright emission, with two of the five transients detected with luminous X-ray emission, ruling out a supernova interpretation. We compare the properties of this sample to previously reported flaring NLSy1 galaxies, and find that they fall into three spectroscopic categories: Transients with 1) Balmer line profiles and Fe II complexes typical of NLSy1s, 2) strong He II profiles, and 3) He II profiles including Bowen fluorescence features. The latter are members of the growing class of AGN flares attributed to enhanced accretion reported by Trakhtenbrot et al. (2019). We consider physical interpretations in the context of related transients from the literature. For example, two of the sources show high amplitude rebrightening in the optical, ruling out a simple tidal disruption event scenario for those transients. We conclude that three of the sample belong to the Trakhtenbrot et al. (2019) class, and two are TDEs in NLSy1s. We also aim to understand why NLSy1s are preferentially the sites of such rapid enhanced flaring activity., Comment: Submitted to ApJ, 27 pages, 11 figures

We report the discovery of ZTF J2243+5242, an eclipsing double white dwarf binary with an orbital period of just $8.8$ minutes, the second known eclipsing binary with an orbital period less than ten minutes. The system likely consists of two low-mass white dwarfs, and will merge in approximately 400,000 years to form either an isolated hot subdwarf or an R Coronae Borealis star. Like its $6.91\, \rm min$ counterpart, ZTF J1539+5027, ZTF J2243+5242 will be among the strongest gravitational wave sources detectable by the space-based gravitational-wave detector The Laser Space Interferometer Antenna (LISA) because its gravitational-wave frequency falls near the peak of LISA's sensitivity. Based on its estimated distance of $d=2120^{+131}_{-115}\,\rm pc$, LISA should detect the source within its first few months of operation, and should achieve a signal-to-noise ratio of $87\pm5$ after four years. We find component masses of $M_A= 0.349^{+0.093}_{-0.074}\,M_\odot$ and $M_B=0.384^{+0.114}_{-0.074}\,M_\odot$, radii of $R_A=0.0308^{+0.0026}_{-0.0025}\,R_\odot$ and $R_B = 0.0291^{+0.0032}_{-0.0024}\,R_\odot$, and effective temperatures of $T_A=22200^{+1800}_{-1600}\,\rm K$ and $T_B=16200^{+1200}_{-1000}\,\rm K$. We determined all of these properties, and the distance to this system, using only photometric measurements, demonstrating a feasible way to estimate parameters for the large population of optically faint ($r>21 \, m_{\rm AB}$) gravitational-wave sources which the Vera Rubin Observatory (VRO) and LISA should identify., Comment: 12 pages, 7 figures, 3 tables, submitted

The current generation of all-sky surveys is rapidly expanding our ability to study variable and transient sources. These surveys, with a variety of sensitivities, cadences, and fields of view, probe many ranges of timescale and magnitude. Data from the Zwicky Transient Facility (ZTF) yields an opportunity to find variables on timescales from minutes to months. In this paper, we present the codebase, ztfperiodic, and the computational metrics employed for the catalogue based on ZTF's Second Data Release. We describe the publicly available, graphical-process-unit optimized period-finding algorithms employed, and highlight the benefit of existing and future graphical-process-unit clusters. We show how generating metrics as input to catalogues of this scale is possible for future ZTF data releases. Further work will be needed for future data from the Vera C. Rubin Observatory's Legacy Survey of Space and Time.

LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility (ZTF). ZTF observed $\sim$\,48\% of S200105ae and $\sim$\,22\% of S200115j's localization probabilities, with observations sensitive to KNe brighter than $-$17.5\,mag fading at 0.5\,mag/day in g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art KN models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with depths of $\rm m_{\rm AB}\approx 22$ mag, attainable in meter-class, wide field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high BH spins, and large neutron star radii.

Using photometry collected with the Zwicky Transient Facility (ZTF), we are conducting an ongoing survey for binary systems with short orbital periods ($P_{\rm b}<1\rm \,hr)$ with the goal of identifying new gravitational-wave sources detectable by the upcoming Laser Interferometer Space Antenna (LISA). Here, we present a sample of fifteen binary systems discovered thus far, with orbital periods ranging from $6.91\rm\,min$ to $56.35\rm\,min$. Of the fifteen systems, seven are eclipsing systems which do not show signs of significant mass transfer. Additionally, we have discovered two AM Canum Venaticorum (AM CVn) systems and six systems exhibiting primarily ellipsoidal variations in their light curves. We present follow-up spectroscopy and high-speed photometry confirming the nature of these systems, estimates of their LISA signal-to-noise ratios (SNR), and a discussion of their physical characteristics., Comment: 32 pages, 12 figures, 7 tables, submitted to ApJ

We present a public catalog of transients from the Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS), a magnitude-limited (m<19 mag in either the g or r filter) survey for extragalactic transients in the ZTF public stream. We introduce cuts on survey coverage, sky visibility around peak light, and other properties unconnected to the nature of the transient, and show that the resulting statistical sample is spectroscopically 97% complete at <18 mag, 93% complete at <18.5 mag, and 75% complete at <19 mag. We summarize the fundamental properties of this population, identifying distinct duration-luminosity correlations in a variety of supernova (SN) classes and associating the majority of fast optical transients with well-established spectroscopic SN types (primarily SN Ibn and II/IIb). We measure the Type Ia SN and core-collapse (CC) SN rates and luminosity functions, which show good consistency with recent work. About 7% of CC SNe explode in very low-luminosity galaxies (M_i > -16 mag), 10% in red-sequence galaxies, and 1% in massive ellipticals. We find no significant difference in the luminosity or color distributions between the host galaxies of Type II and Type Ib/c supernovae, suggesting that line-driven wind stripping does not play a major role in the loss of the hydrogen envelope from their progenitors. Future large-scale classification efforts with ZTF and other wide-area surveys will provide high-quality measurements of the rates, properties, and environments of all known types of optical transients and limits on the existence of theoretically predicted but as of yet unobserved explosions., Comment: Accepted to ApJ. Live BTS statistics are available online at https://sites.astro.caltech.edu/ztf/bts/bts.php and an interactive catalog of our sample is available at https://sites.astro.caltech.edu/ztf/bts/explorer.php

Spectroscopic detection of narrow emission lines traces the presence of circumstellar mass distributions around massive stars exploding as core-collapse supernovae. Transient emission lines disappearing shortly after the supernova explosion suggest that the spatial extent of such material is compact, and hence imply an increased mass loss shortly prior to explosion. Here, we present a systematic survey for such transient emission lines (Flash Spectroscopy) among Type II supernovae detected in the first year of the Zwicky Transient Facility (ZTF) survey. We find that at least six out of ten events for which a spectrum was obtained within two days of estimated explosion time show evidence for such transient flash lines. Our measured flash event fraction ($>30\%$ at $95\%$ confidence level) indicates that elevated mass loss is a common process occurring in massive stars that are about to explode as supernovae.

We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections., Comment: 26 pages, 8 figures, 3 Tables, accepted for publication in ApJL

The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded $>11,200$ candidates, 24 of which passed quality checks and strict selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting, catalog cross-matching, and study of their color evolution. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. In addition, we identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB190106A, and the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat and linearly decaying light curves and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be $R < 1775$ Gpc$^{-3}$ yr$^{-1}$ at 95% confidence level by requiring at least 2 high-significance detections. By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of $R < 4029$ Gpc$^{-3}$ yr$^{-1}$., Comment: Submitted for publication in ApJ

We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD., Comment: accepted for publication in ApJL, 11 pages, 5 figures and 1 tables. arXiv admin note: text overlap with arXiv:2002.01485

This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end., Comment: Editors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist. "EVN Vision Document (2020)", 7 chapters, 3 appendices, 191 pages, 64 figures. Chapter coordinators and all the contributors are listed in the document