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 fringez (https://github.com/MichaelMedford/fringez).

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a 30 deg2 field at ΔT = 2, 9, and 33 days post-merger at a frequency of 944 MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers 89% of the LIGO/Virgo localization region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, AT2019osy, we performed multiwavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from NSBH mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational-wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019 May 10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a binary neutron star (BNS) merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 92 Mpc. We triggered Target-of-Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4 m Blanco Telescope at Cerro Tololo Inter-American Observatory in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hr of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was reclassified 1.7 days after the merger, after our observations were completed, with a "BNS merger" probability reduced from 98% to 42% in favor of a "terrestrial classification.

The discovery of a transient kilonova following the gravitational-wave (GW) event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up across the electromagnetic spectrum. In the southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4 m telescope is well suited to this task, as it is able to cover wide fields quickly while still achieving the depths required to find kilonovae like the one accompanying GW170817 to ∼500 Mpc, the binary neutron star (NS) horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility follow-up by the Global Relay of Observatories Watching Transients Happen collaboration, we describe the observations and automated data movement, data reduction, candidate discovery, and vetting pipeline of our target-of-opportunity DECam observations of S190426c, the first possible NS-black hole merger detected in GWs. Starting 7.5 hr after S190426c, over 11.28 hr of observations, we imaged an area of 525 deg2 (r band) and 437 deg2 (z band); this was 16.3% of the total original localization probability, and nearly all of the probability visible from the southern hemisphere. The machine-learning-based pipeline was optimized for fast turnaround, delivering transients for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hr before the end of our observations. One hour after our data-taking ended (roughly 20 hr after the announcement of S190426c), LVC released a refined skymap that reduced the probability coverage of our observations to 8.0%, demonstrating a critical need for localization updates on shorter (∼hour) timescales. Our observations yielded no detection of a bona fide counterpart to m z = 21.7 and m r = 22.2 at the 5σ level of significance, consistent with the refined LVC positioning. We view these observations and rapid inferencing as an important real-world test for this novel end-to-end wide-field pipeline.

Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. Here we describe how high-cadence optical observations with the Zwicky Transient Facility, with its unparalleled large field of view, led to the detection of a multiply imaged type Ia supernova, SN Zwicky, also known as SN 2022qmx. Magnified nearly 25-fold, the system was found thanks to the standard candle nature of type Ia supernovae. High-spatial-resolution imaging with the Keck telescope resolved four images of the supernova with very small angular separation, corresponding to an Einstein radius of only θE = 0.167″ and almost identical arrival times. The small θE and faintness of the lensing galaxy are very unusual, highlighting the importance of supernovae to fully characterize the properties of galaxy-scale gravitational lenses, including the impact of galaxy substructures.

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 day < t 1/2 < 12 days, 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 predominantly featureless spectra and luminous radio emission: AT2018lug (The Koala) and AT2020xnd (The Camel). Seven (five) did not have a definitive classification: AT 2020bdh showed tentative broad Hα in emission, and AT 2020bot showed unidentified broad features and was 10 kpc offset from the center of an early-type galaxy. Ten (eight) 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 24 of the transients. All X-ray and radio observations resulted in nondetections; 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.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star-black hole (NSBH) merger, or potentially a low-mass binary BH merger. Due to the low false-alarm rate and the precise localization (23 deg2 at 90%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to an NSBH merger. To search for potential counterparts, the GROWTH Collaboration performed real-time image subtraction on six nights of public Dark Energy Camera images acquired in the 3 weeks following the merger, covering >98% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be M ej < 0.04 M o˙ at polar viewing angles, or M ej < 0.03 M o˙ if the opacity is κ < 2 cm2g-1. Assuming a tidal deformability for the NS at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be χ < 0.7 for mass ratios Q < 6, with weaker constraints for more compact NSs.

We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF 19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of Mr = -13 (νLν = 9 × 106 Lo), in the luminosity gap between novae and typical supernovae (SNe). Spectra during the outburst show a red continuum, Balmer emission with a velocity width of ≈400 km s-1, Ca ii and [Ca ii] emission, and absorption features characteristic of an F-type supergiant. The spectra and multiband light curves are similar to the so-called "SN impostors" and intermediate-luminosity red transients (ILRTs). We directly identify the likely progenitor in archival Spitzer Space Telescope imaging with a 4.5 μm luminosity of M[4.5] ≈ -12.2 mag and a [3.6]-[4.5] color redder than 0.74 mag, similar to those of the prototype ILRTs SN 2008S and NGC 300 OT2008-1. Intensive monitoring of M51 with Spitzer further reveals evidence for variability of the progenitor candidate at [4.5] in the years before the OT. The progenitor is not detected in pre-outburst Hubble Space Telescope optical and near-IR images. The optical colors during outburst combined with spectroscopic temperature constraints imply a higher reddening of E(B - V) ≈ 0.7 mag and higher intrinsic luminosity of Mr ≈ -14.9 mag (νLν = 5.3 × 107 Lo) near peak than seen in previous ILRT candidates. Moreover, the extinction estimate is higher on the rise than on the plateau, suggestive of an extended phase of circumstellar dust destruction. These results, enabled by the early discovery of M51 OT2019-1 and extensive pre-outburst archival coverage, offer new clues about the debated origins of ILRTs and may challenge the hypothesis that they arise from the electron-capture induced collapse of extreme asymptotic giant branch stars.

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising (1.4 ± 0.1 mag hr-1) and luminous (Mg,peak = -20 mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity (Lbol ≳ 3 × 1044 erg s-1), the short rise time (trise = 3 days in g band), and the blue colors at peak (g-r ∼ -0.4) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature (Teff ≳ 40,000 K) spectra of a stripped-envelope SN. A retrospective search revealed luminous (Mg ∼ Mr ≈ mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release E γ,iso < 4.9 × 10 48 erg, a limit on X-ray emission LX < 1040 erg s-1, and a limit on radio emission ν Lν ≲ 1037 erg s-1. Taken together, we find that the early (< 10 days) data are best explained by shock breakout in a massive shell of dense circumstellar material (0.02 M⊙) at large radii (3 × 1014 cm) that was ejected in eruptive pre-explosion mass-loss episodes. The late-time (> 10 days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56.

Cosmic neutrinos provide a unique window into the otherwise hidden mechanism of particle acceleration in astrophysical objects. The IceCube Collaboration recently reported the likely association of one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the likely association of a radio-emitting tidal disruption event, AT2019dsg, with a second high-energy neutrino. AT2019dsg was identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility. The probability of finding any coincident radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multizone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for petaelectronvolt neutrino production. Assuming that the association is genuine, our observations suggest that tidal disruption events with mildly relativistic outflows contribute to the cosmic neutrino flux.