Search

The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that $-21$ can contribute no more than 87% of the total, while transients brighter than $-22$ can contribute no more than 58% of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26% of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory.

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

Tumer et al. 2022, ATel #15171, have recently reported the discovery of an X-ray source, NuSTAR J053449+2126.0, during a calibration observation which took place on 25 April 2020. We scan the Zwicky Transient Facility (ZTF) alerts and archival photometry to determine the nature of the source. Palomar Gattini-IR is searched as well. We identify no obvious counterpart candidate. Follow-up X-ray and optical studies are needed to determine the true counterpart., Comment: Accepted and Published in RNAAS

The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community., Comment: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observations, ii) the first broadband infrared spectrum ($\approx1-5\,\mu$m; $R\approx 2000$) obtained with SpeX on the NASA Infrared Telescope Facility and iii) $J$-band light curve from the Palomar Gattini-IR survey. The infrared spectrum is characterized by i) deep water absorption features (H$_2$O index $\approx 40$%), ii) strong TiO, VO and CO features, and iii) weak/absent CN lines. We show that these features are inconsistent with known red supergiants, but suggest a M8-9 III type O-rich Mira donor star. We report the discovery of large amplitude ($\Delta J\approx3.5$ mag) periodic photometric variability suggesting a pulsation period of $621\pm36\,{\rm(systematic)}\pm8\,{\rm(statistical)}$ days, which we use to constrain the donor to be a relatively luminous Mira ($M_K=-8.6\pm0.3$ mag) at a distance of $3.6^{+0.8}_{-0.7}$ kpc. Comparing these characteristics to recent models, we find the donor to be consistent with a $\approx 3-5$ M$_\odot$ star at an age of $\approx 0.1-0.3$ Gyr. Together, we show that Sct X-1 was previously mis-classified as an evolved High Mass X-ray Binary; instead it is an intermediate mass system with the first confirmed Mira donor in an X-ray binary. We discuss the implications of Mira donors in symbiotic X-ray binaries, and highlight the potential of wide-field infrared time domain surveys and broadband infrared spectroscopy to unveil their demographics., Comment: 10 pages, 4 figures. Accepted for publication in ApJ Letters

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program., Comment: 69 pages, 43 figures

The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources., Comment: 20 pages, 6 figures, 6 tables

We report second-timescale infrared photometry of the nova V1674 Her using Palomar Gattini-IR. These observations constitute the first infrared and highest temporal resolution data (resolution of ~ 0.84 s) of the nova reported to date. PGIR observed in this fast readout mode for more than an hour on three nights between 3 and 6 days after discovery. We searched for periodic variability using a Lomb-Scargle periodogram and did not detect anything down to a three sigma upper limit of 0.093 mag. This suggests that the periodic variability detected in the nova by Patterson et al. 2021 was lower by at least a factor of about 1.65 in the first week of the eruption., Comment: 3 pages, 1 figure

We present observations of SN 2020fqv, a Virgo-cluster Type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 days post-explosion; ground-based spectroscopic observations starting 1.1~days post-explosion; along with extensive photometric observations. Massive stars have complicated mass-loss histories leading up to their death as CCSNe, creating circumstellar medium (CSM) with which the SNe interact. Observations during the first few days post-explosion can provide important information about the mass-loss rate during the late stages of stellar evolution. Model fits to the quasi-bolometric light curve of SN 2020fqv reveal ~0.23 $M_{\odot}$ of CSM confined within ~1450 $R_{\odot}$ ($10^{14}$ cm) from its progenitor star. Early spectra (<4 days post-explosion), both from HST and ground-based observatories, show emission features from high-ionization metal species from the outer, optically thin part of this CSM. We find that the CSM is consistent with an eruption caused by the injection of $\sim$$5\times 10^{46}$ erg into the stellar envelope $\sim$300 days pre-explosion, potentially from a nuclear burning instability at the onset of oxygen burning. Light-curve fitting, nebular spectroscopy, and pre-explosion \textit{HST} imaging consistently point to a red supergiant (RSG) progenitor with $M_{\rm ZAMS}$$\approx$$13.5$--$15 \, M_{\odot}$, typical for SN~II progenitor stars. This finding demonstrates that a typical RSG, like the progenitor of SN 2020fqv, has a complicated mass-loss history immediately before core collapse., Comment: Published in MNRAS

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

We present observations of three Core-collapse supernovae (CCSNe) in elliptical hosts, detected by the Zwicky Transient Facility Bright Transient Survey (BTS). SN 2019ape is a SN Ic that exploded in the main body of a typical elliptical galaxy. Its properties are consistent with an explosion of a regular SN Ic progenitor. A secondary g-band light curve peak could indicate interaction of the ejecta with circumstellar material (CSM). An H$\alpha$-emitting source at the explosion site suggests a residual local star formation origin. SN 2018fsh and SN 2020uik are SNe II which exploded in the outskirts of elliptical galaxies. SN 2020uik shows typical spectra for SNe II, while SN 2018fsh shows a boxy nebular H$\alpha$ profile, a signature of CSM interaction. We combine these 3 SNe with 7 events from the literature and analyze their hosts as a sample. We present multi-wavelength photometry of the hosts, and compare this to archival photometry of all BTS hosts. Using the spectroscopically complete BTS we conclude that $0.3\%^{+0.3}_{-0.1}$ of all CCSNe occur in elliptical galaxies. We derive star-formation rates and stellar masses for the host-galaxies and compare them to the properties of other SN hosts. We show that CCSNe in ellipticals have larger physical separations from their hosts compared to SNe Ia in elliptical galaxies, and discuss implications for star-forming activity in elliptical galaxies., Comment: Comments are welcome. Submitted to ApJ

A core-collapse supernova occurs when exothermic fusion ceases in the core of a massive star, typically due to exhaustion of nuclear fuel. Theory predicts that fusion could be interrupted earlier, by merging of the star with a compact binary companion. We report a luminous radio transient, VT J121001+495647, found in the Very Large Array Sky Survey. The radio emission is consistent with supernova ejecta colliding with a dense shell of material, potentially ejected by binary interaction in the centuries prior to explosion. We associate the supernova with an archival X-ray transient, which implies a relativistic jet was launched during the explosion. The combination of an early relativistic jet and late-time dense interaction is consistent with expectations for a merger-driven explosion., Comment: 33 pages, 4 figures. v2: Updated fig 1 to include panel D

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

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

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

We present optical and infrared (IR) light curves of NaSt1, also known as Wolf-Rayet (WR) 122, with observations from Palomar Gattini-IR (PGIR), the Zwicky Transient Facility (ZTF), the Katzman Automatic Imaging Telescope (KAIT), the Asteroid Terrestrial-impact Last Alert System (ATLAS), and the All-Sky Automated Survey for Supernovae (ASAS-SN). We identify a $P=309.7\pm0.7$ d photometric period from the optical and IR light curves that reveal periodic, sinusoidal variability between 2014 July and 2021 July. We also present historical IR light curves taken between 1983 July and 1989 May, which show variability consistent with the period of the present-day light curves. In the past, NaSt1 was brighter in the $J$ band with larger variability amplitudes than the present-day PGIR values, suggesting that NaSt1 exhibits variability on longer ($\gtrsim$ decade) timescales. Sinusoidal fits to the recent optical and IR light curves show that the amplitude of NaSt1's variability differs at various wavelengths and also reveal significant phase offsets of $17.0\pm2.5$ d between the ZTF $r$ and PGIR $J$ light curves. We interpret the $310$ d photometric period from NaSt1 as the orbital period of an enshrouded massive binary. We suggest that the photometric variability of NaSt1 may arise from variations in the line-of-sight optical depth toward circumstellar optical/IR emitting regions throughout its orbit due to colliding-wind dust formation. We speculate that past mass transfer in NaSt1 may have been triggered by Roche-lobe overflow (RLOF) during an eruptive phase of an Ofpe/WN9 star. Lastly, we argue that NaSt1 is no longer undergoing RLOF mass transfer., Comment: 20 Pages, 7 Figures, 5 Tables, 1 Table Attachment, Accepted to ApJ

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

Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting $\lesssim 1$ day, reaching a peak absolute magnitude of $M_V = -9.2$ mag, and a very fast decline, fading by 2 mag over 3.5 days. A steep dropoff in the light curves after 71 days and the rapid decline timescale suggest a low-mass ejection from a massive WD with $M_{\rm WD} \gtrsim 1.2~M_{\odot}$. We present an unprecedented view of the early spectroscopic evolution of such an event. Three spectra prior to the peak reveal a complex, multicomponent outflow giving rise to internal collisions and shocks in the ejecta of an He/N-class nova. We identify a coincident IR-variable counterpart in the extensive preeruption coverage of the transient location and infer the presence of a symbiotic progenitor system with an O-rich asymptotic-giant-branch donor star, as well as evidence for an earlier UV-bright outburst in 2014. We suggest that AT 2019qyl is analogous to the subset of Galactic recurrent novae with red-giant companions such as RS Oph and other embedded nova systems like V407 Cyg. Our observations provide new evidence that internal shocks between multiple, distinct outflow components likely contribute to the generation of the shock-powered emission from such systems., Comment: 32 pages, 10 figures, accepted in ApJ

Luminous Red Novae (LRNe) are astrophysical transients associated with the partial ejection of a binary system's common envelope (CE) shortly before its merger. Here we present the results of our photometric and spectroscopic follow-up campaign of AT2018bwo (DLT18x), a LRN discovered in NGC45, and investigate its progenitor system using binary stellar-evolution models. The transient reached a peak magnitude of $M_r=-10.97\pm0.11$ and maintained this brightness during its optical plateau of $t_p = 41\pm5$days. During this phase, it showed a rather stable photospheric temperature of ~3300K and a luminosity of ~$10^{40}$erg/s. The photosphere of AT2018bwo at early times appeared larger and cooler than other similar LRNe, likely due to an extended mass-loss episode before the merger. Towards the end of the plateau, optical spectra showed a reddened continuum with strong molecular absorption bands. The reprocessed emission by the cooling dust was also detected in the mid-infrared bands ~1.5 years after the outburst. Archival Spitzer and Hubble Space Telescope data taken 10-14 years before the transient event suggest a progenitor star with $T_{prog}\sim 6500$K, $R_{prog}\sim 100R_{\odot}$ and $L_{prog}\sim 2\times10^4L_{\odot}$, and an upper limit for optically thin warm (1000 K) dust mass of $M_d<10^{-6}M_{\odot}$. Using stellar binary-evolution models, we determined the properties of binary systems consistent with the progenitor parameter space. For AT2018bwo, we infer a primary mass of 12-16 $M_{\odot}$, which is 9-45% larger than the ~11$M_{\odot}$ obtained using single-star evolution models. The system, consistent with a yellow-supergiant primary, was likely in a stable mass-transfer regime with -2.4

Massive stars die an explosive death as a core-collapse supernova (CCSN). The exact physical processes that cause the collapsing star to rebound into an explosion are not well-understood, and the key in resolving this issue may lie in the measurement of the shape of CCSNe ejecta. Spectropolarimetry is the only way to perform this measurement for CCSNe outside of the Milky Way and Magellanic Clouds. We present an infrared (IR) spectropolarimetric detection of a CCSN, enabled by the new highly sensitive WIRC+Pol instrument at Palomar Observatory, that can observe CCSNe (M = -17 mags) out to 20 Mpc to ~0.1% polarimetric precision. IR spectropolarimetry is less affected than optical by dust scattering in the circumstellar and interstellar media, thereby providing a more unbiased probe of the intrinsic geometry of the SN ejecta. SN 2018hna, a SN 1987A-like explosion, shows 2.0+-0.3% continuum polarization in the J band oriented at ~160 degree on-sky at 182 d after the explosion. Assuming prolate geometry like in SN 1987A, we infer an ejecta axis ratio of <0.48 with the axis of symmetry pointing at 70 degree position angle. The axis ratio is similar to that of SN 1987A suggesting that they may share intrinsic geometry and inclination angle. Our data do not rule out oblate ejecta. We also observe one other core-collapse and two thermonuclear SNe in the J band. SN 2020oi, a stripped-envelope Type Ic SN in Messier 100 has p = 0.37+-0.09% at peak light, indicative of either a 10% asymmetry or host interstellar polarization. The SNe Ia, 2019ein and 2020ue have p < 0.33% and < 1.08% near peak light, indicative of asymmetries of less than 10% and 20%, respectively.

Pulsar radio emission undergoes dispersion due to the presence of free electrons in the interstellar medium (ISM). The dispersive delay in the arrival time of pulsar signal changes over time due to the varying ISM electron column density along the line of sight. Correcting for this delay accurately is crucial for the detection of nanohertz gravitational waves using Pulsar Timing Arrays. In this work, we present in-band and inter-band DM estimates of four pulsars observed with uGMRT over the timescale of a year using two different template alignment methods. The DMs obtained using both these methods show only subtle differences for PSR 1713+0747 and J1909$-$3744. A considerable offset is seen in the DM of PSR J1939+2134 and J2145$-$0750 between the two methods. This could be due to the presence of scattering in the former and profile evolution in the latter. We find that both methods are useful but could have a systematic offset between the DMs obtained. Irrespective of the template alignment methods followed, the precision on the DMs obtained is about $10^{-3}$ pc cm$^{-3}$ using only BAND3 and $10^{-4}$ pc cm$^{-3}$ after combining data from BAND3 and BAND5 of the uGMRT. In a particular result, we have detected a DM excess of about $5\times10^{-3}$ pc cm$^{-3}$ on 24 February 2019 for PSR J2145$-$0750. This excess appears to be due to the interaction region created by fast solar wind from a coronal hole and a coronal mass ejection (CME) observed from the Sun on that epoch. A detailed analysis of this interesting event is presented., Comment: 11 pages, 6 figures, 2 tables. Accepted by A&A

The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for Type Ia supernovae. The rate has been previously estimated in the range of $\approx10-300$ yr$^{-1}$, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galaxies; both methods are subject to debatable assumptions. The total discovery rate of optical novae remains much smaller ($\approx5-10$ yr$^{-1}$) than these estimates, even with the advent of all-sky optical time domain surveys. Here, we present a systematic sample of 12 spectroscopically confirmed Galactic novae detected in the first 17 months of Palomar Gattini-IR (PGIR), a wide-field near-infrared time domain survey. Operating in $J$-band ($\approx1.2$ $\mu$m) that is relatively immune to dust extinction, the extinction distribution of the PGIR sample is highly skewed to large extinction values ($> 50$% of events obscured by $A_V\gtrsim5$ mag). Using recent estimates for the distribution of mass and dust in the Galaxy, we show that the observed extinction distribution of the PGIR sample is commensurate with that expected from dust models. The PGIR extinction distribution is inconsistent with that reported in previous optical searches (null hypothesis probability $< 0.01$%), suggesting that a large population of highly obscured novae have been systematically missed in previous optical searches. We perform the first quantitative simulation of a $3\pi$ time domain survey to estimate the Galactic nova rate using PGIR, and derive a rate of $\approx 46.0^{+12.5}_{-12.4}$ yr$^{-1}$. Our results suggest that all-sky near-infrared time-domain surveys are well poised to uncover the Galactic nova population., Comment: 25 pages, 12 figures. Submitted to ApJ

We report the discovery of a likely outbursting Class I young stellar object, associated with the star-forming region NGC 281-W (distance $\sim 2.8$ kpc). The source is currently seen only at infrared wavelengths, appearing in both the Palomar Gattini InfraRed ($1.2~\mu$m) and the Near Earth Object Widefield Infrared Survey Explorer ($3.4$ and $4.6~\mu$m) photometric time-domain surveys. Recent near-infrared imaging reveals a new, extended scattered light nebula. Recent near-infrared spectroscopy confirms the similarity of PGIR 20dci to FU Ori type sources, based on strong molecular absorption in CO, H$_2$O, and OH, weak absorption in several atomic lines, and a warm wind/outflow as indicated by a P Cygni profile in the HeI 10830 A line. This is a rare case of an FU Ori star with a well-measured long term photometric rise before a sharper outburst, and the second instance of an FU Ori star with a documented two-step brightening in the mid-infrared., Comment: accepted to AAS Journals on 2021, February 5

We are undertaking the first systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, beginning with IR light curves from the Palomar Gattini IR (PGIR) survey. PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that is surveying 18000 deg$^{2}$ of the northern sky ($\delta>-28^{o}$) at a cadence of 2 days. We present PGIR light curves for 922 RCB candidates selected from a mid-IR color-based catalog (Tisserand et al. 2020). Of these 922, 149 are promising RCB candidates as they show pulsations or declines similar to RCB stars. Majority of the candidates that are not RCB stars are either long period variables (LPVs) or RV-Tauri stars. We identify IR color-based criteria to better distinguish between RCB stars and LPVs. As part of a pilot spectroscopic run, we obtained NIR spectra for 26 out of the 149 promising candidates and spectroscopically confirm 11 new RCB stars. We detect strong He I $\lambda 10830$ features in spectra of all RCB stars, likely originating within high velocity (200-400 km-s$^{-1}$) winds in their atmospheres. 9 of these RCB stars show $^{12}$C$^{16}$O and $^{12}$C$^{18}$O molecular absorption features, suggesting that they are formed through a white dwarf merger. We detect quasi-periodic pulsations in the light curves of 5 RCB stars. The periods range between 30-125 days and likely originate from the strange-mode instability in these stars. Our pilot run results motivate a dedicated IR spectroscopic campaign to classify all RCB candidates., Comment: 20 pages, 14 figures, Submitted to ApJ

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

We present Spitzer/IRAC observations of dust formation from six extragalactic carbon-rich Wolf-Rayet (WC) binary candidates in low-metallicity (Z $\lesssim0.65$ Z$_\odot$) environments using multi-epoch mid-infrared (IR) imaging data from the SPitzer InfraRed Intensive Transients Survey (SPIRITS). Optical follow-up spectroscopy of SPIRITS~16ln, 19q, 16df, 18hb, and 14apu reveals emission features from C IV $\lambda5801\text{-}12$~and/or the C III-IV $\lambda4650$ He II $\lambda4686$~blend that are consistent with early-type WC stars. We identify SPIRITS~16ln as the variable mid-IR counterpart of the recently discovered colliding-wind WC4+O binary candidate, N604-WRXc, located in the sub-solar metallicity NGC 604 H II~region in M33. We interpret the mid-IR variability from SPIRITS~16ln as a dust-formation episode in an eccentric colliding-wind WC binary. SPIRITS~19q, 16df, 14apu, and 18hb exhibit absolute [3.6] magnitudes exceeding one of most IR-luminous dust-forming WC systems known, WR~104 (M$_\mathrm{[3.6]}\lesssim-12.3$). An analysis of dust formation in the mid-IR outburst from SPIRITS~19q reveals a high dust production rate of $\dot{M}_d\gtrsim2\times10^{-6}$ M$_\odot$ yr$^{-1}$, which may therefore exceed that of the most efficient dust-forming WC systems known. We demonstrate that efficient dust-formation is feasible from early-type WC binaries in the theoretical framework of colliding-wind binary dust formation if the systems host an O-type companion with a high mass-loss rate ($\dot{M}\gtrsim1.6\times10^{-6}$ M$_\odot$ yr$^{-1}$). This efficient dust-formation from early-type WC binaries highlights their potential role as significant sources of dust in low-metallicity environments., Comment: 28 pages, 10 Figures, and 5 Tables. Accepted to ApJ on 2021 Jan 4 (Replacement of version submitted to arXiv on 2020 Nov 19)

We present the discovery of a luminous X-ray transient, serendipitously detected by Swift's X-ray Telescope (XRT) on 2020 February 5, located in the nucleus of the galaxy SDSS J143359.16+400636.0 at z=0.099 (luminosity distance $D_{\rm L}=456$ Mpc). The transient was observed to reach a peak luminosity of $\sim10^{44}$ erg s$^{-1}$ in the 0.3--10 keV X-ray band, which was $\sim20$ times more than the peak optical/UV luminosity. Optical, UV, and X-ray lightcurves from the Zwicky Transient Facility (ZTF) and Swift show a decline in flux from the source consistent with $t^{-5/3}$, and observations with NuSTAR and Chandra show a soft X-ray spectrum with photon index $\Gamma=2.9\pm0.1$. The X-ray/UV properties are inconsistent with well known AGN properties and have more in common with known X-ray tidal disruption events (TDE), leading us to conclude that it was likely a TDE. The broadband spectral energy distribution (SED) can be described well by a disk blackbody model with an inner disk temperature of $7.3^{+0.3}_{-0.8}\times10^{5}$ K, with a large fraction ($>40$%) of the disk emission up-scattered into the X-ray band. An optical spectrum taken with Keck/LRIS after the X-ray detection reveals LINER line ratios in the host galaxy, suggesting low-level accretion on to the supermassive black hole prior to the event, but no broad lines or other indications of a TDE were seen. The stellar velocity dispersion implies the mass of the supermassive black hole powering the event is log($M_{\rm BH}$/$M_{\odot}$)$=7.41\pm0.41$, and we estimate that at peak the Eddington fraction of this event was $\sim$50%. This likely TDE was not identified by wide-field optical surveys, nor optical spectroscopy, indicating that more events like this would be missed without wide-field UV or X-ray surveys., Comment: Accepted for publication by ApJ. Accepted version now replaces initial submission

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.

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.

The nature of the peculiar `Ca-rich' SN 2019ehk in the nearby galaxy M100 remains unclear. Its origin has been debated as either a stripped core-collapse supernova or a thermonuclear helium detonation event. Here, we present very late-time photometry of the transient obtained with the Keck I telescope at $\approx280$ days from peak light. Using the photometry to perform accurate flux calibration of a contemporaneous nebular phase spectrum, we measure an [O I] luminosity of $(0.19-1.08)\times10^{38}$ erg s$^{-1}$ and [Ca II] luminosity of $(2.7-15.6)\times10^{38}$ erg s$^{-1}$ over the range of the uncertain extinction along the line of sight. We use these measurements to derive lower limits on the synthesized oxygen mass of $\approx0.004-0.069$ M$_\odot$. The oxygen mass is a sensitive tracer of the progenitor mass for core-collapse supernovae, and our estimate is consistent with explosions of very low mass CO cores of $1.45-1.5$ M$_\odot$, corresponding to He core masses of $\approx1.8-2.0$ M$_\odot$. We present high quality peak light optical spectra of the transient and highlight features of hydrogen in both the early (`flash') and photospheric phase spectra, that suggest the presence of $\gtrsim0.02$ M$_\odot$ of hydrogen in the progenitor at the time of explosion. The presence of H, together with the large [Ca II]/[O I] ratio ($\approx10-15$) in the nebular phase is consistent with SN 2019ehk being a Type IIb core-collapse supernova from a stripped low mass ($\approx9-9.5$ M$_\odot$) progenitor, similar to the Ca-rich SN IIb iPTF 15eqv. These results provide evidence for a likely class of `Ca-rich' core-collapse supernovae from stripped low mass progenitors in star forming environments, distinct from the thermonuclear Ca-rich gap transients in old environments., Comment: 12 pages, 5 figures, Accepted for publication in ApJL

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.

Several thousand core-collapse supernovae (CCSNe) of different flavors have been discovered so far. However, identifying their progenitors has remained an outstanding open question in astrophysics. Studies of SN host galaxies have proven to be powerful in providing constraints on the progenitor populations. In this paper, we present all CCSNe detected between 2009 and 2017 by the Palomar Transient Factory. This sample includes 888 SNe of 12 distinct classes out to redshift $z\approx1$. We present the photometric properties of their host galaxies from the far-ultraviolet to the mid-infrared and model the host-galaxy spectral energy distributions to derive physical properties. The galaxy mass functions of Type Ic, Ib, IIb, II, and IIn SNe ranges from $10^{5}$ to $10^{11.5}~M_\odot$, probing the entire mass range of star-forming galaxies down to the least-massive star-forming galaxies known. Moreover, the galaxy mass distributions are consistent with models of star-formation-weighted mass functions. Regular CCSNe are hence direct tracers of star formation. Small but notable differences exist between some of the SN classes. Type Ib/c SNe prefer galaxies with slightly higher masses (i.e., higher metallicities) and star-formation rates than Type IIb and II SNe. These differences are less pronounced than previously thought. H-poor SLSNe and SNe~Ic-BL are scarce in galaxies above $10^{10}~M_\odot$. Their progenitors require environments with metallicities of $<0.4$ and $<1$ solar, respectively. In addition, the hosts of H-poor SLSNe are dominated by a younger stellar population than all other classes of CCSNe. Our findings corroborate the notion that low-metallicity \textit{and} young age play an important role in the formation of SLSN progenitors., Comment: 24 pages main text, 14 figures, 9 Tables, catalogue available at http://www.github.com/steveschulze/PTF

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

The Galactic magnetar SGR1935+2154 has been reported to produce the first known example of a bright millisecond duration radio burst (FRB 200428) similar to the cosmological population of fast radio bursts (FRBs), bolstering the association of FRBs to active magnetars. The detection of a coincident bright X-ray burst has revealed the first observed multi-wavelength counterpart of a FRB. However, the search for similar emission at optical wavelengths has been hampered by the high inferred extinction on the line of sight. Here, we present results from the first search for second-timescale emission from the source at near-infrared wavelengths using the Palomar Gattini-IR observing system in J-band, made possible by a recently implemented detector read-out mode that allowed for short exposure times of 0.84 s with 99.9% observing efficiency. With a total observing time of 12 hours (47728 images) on source, we place median $3\,\sigma$ limits on the second-timescale emission of $< 20$ mJy (13.1 AB mag). We present non-detection limits from epochs of four simultaneous X-ray bursts detected by the Insight-{\it HXMT} and {\it NuSTAR} telescopes during our observing campaign. The limits translate to an extinction corrected fluence limit of $< 125$ Jy ms for an estimated extinction of $A_J = 2.0$ mag. These limits provide the most stringent constraints to date on the fluence of flares at frequencies of $\sim 10^{14}$ Hz, and constrain the ratio of the near-infrared (NIR) fluence to that of coincident X-ray bursts to $R_{\rm NIR} < 2.5 \times 10^{-2}$. Our observations were sensitive enough to easily detect a near-infrared counterpart of FRB 200428 if the NIR emission falls on the same power law as that observed across its radio to X-ray spectrum. The non-detection of NIR emission around the coincident X-ray bursts constrains the fluence index of the brightest burst to be steeper than $0.35$., Comment: 10 pages, 4 figures, submitted to ApJL. Comments welcome

We introduce pinta, a pipeline for reducing the upgraded Giant Metre-wave Radio Telescope (uGMRT) raw pulsar timing data, developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. We also discuss a novel and independent determination of the relative time offsets between the different back-end modes of uGMRT and the interpretation of the uGMRT observation frequency settings, and their agreement with results obtained from engineering tests. Further, we demonstrate the capability of pinta to generate data products which can produce high-precision TOAs using PSR J1909-3744 as an example. These results are crucial for performing precision pulsar timing with the uGMRT.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged), Comment: Submitted to ApJ

We present ZTF20aajnksq (AT2020blt), a fast-fading ($\Delta r=2.4$ mag in $\Delta t=1.3$ days) red ($g-r\approx0.6$ mag) and luminous ($M_{1626}=-25.9$) optical transient at $z=2.9$ discovered by the Zwicky Transient Facility (ZTF). AT2020blt 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_\mathrm{j}=1$ day (observer-frame); (2) a luminous $(L_X = 10^{46}$ $\mathrm{erg}$ $\mathrm{s}^{-1})$ X-ray counterpart; and (3) luminous ($L_\nu = 4 \times 10^{31}$ $\mathrm{erg}$ $\mathrm{sec}^{-1}$ $\mathrm{Hz}^{-1}$ at 10 GHz) radio emission. However, no GRB was detected in the 0.74d between the last ZTF non-detection ($r > 20.64$) and the first ZTF detection ($r = 19.57$), with an upper limit on the isotropic-equivalent gamma-ray energy release of $E_{\gamma,\mathrm{iso}} < 7 \times 10^{52}$ erg. AT2020blt 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 AT2020blt are consistent with a classical (initial Lorentz factor $\Gamma_0 \gtrsim 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 AT2020blt 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. We conclude by discussing the status and future of fast-transient searches in wide-field high-cadence optical surveys., Comment: 19 pages, 4 figures, 1 table. Version accepted by ApJ on 17 Oct 2020

We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture., Comment: 29 pages, 23 figures, submitted to ApJ

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations., Comment: 30 figures, 6 tables. Submitted to ApJ

We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% on the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modelling of these events., Comment: Submitted to ApJ. Comments are welcome

While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress towards this promise. Here present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, as well as longer rise times and higher disruption rates compared to the rest of the sample. The Bowen fluorescence mechanism requires a high density which can be reached at smaller radii, which in turn yields longer diffusion timescales. Thus, the difference in rise times suggests the pre-peak TDE light curves are governed not by the fallback timescale, but instead by the diffusion of photons through the tidal debris. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities and the lowest rates. We also report, for the first time, the detection of soft X-ray flares from a TDE on day timescales. Based on the fact the flares peak at a luminosity similar to the optical/UV blackbody luminosity, we attribute them to brief glimpses through a reprocessing layer that otherwise obscures the inner accretion flow., Comment: 30 pages, 18 figures, 8 tables, to be submitted to ApJ, comments welcome

We present ZTF18aaqjovh (SN 2018bvw), a high-velocity ("broad-lined") stripped-envelope (Type Ic) supernova (Ic-BL SN) discovered in the Zwicky Transient Facility one-day cadence survey. ZTF18aaqjovh shares a number of features in common with engine-driven explosions: the photospheric velocity and the shape of the optical light curve are very similar to that of the Type Ic-BL SN 1998bw, which was associated with a low-luminosity gamma-ray burst (LLGRB) and had relativistic ejecta. However, the radio luminosity of ZTF18aaqjovh is almost two orders of magnitude fainter than that of SN 1998bw at the same velocity phase, and the shock velocity is at most mildly relativistic (v=0.06-0.4c). A search of high-energy catalogs reveals no compelling GRB counterpart to ZTF18aaqjovh, and the limit on the prompt GRB luminosity of $L_{\gamma,\mathrm{iso}} \approx 1.6 \times 10^{48}$ erg/sec excludes a classical GRB but not an LLGRB. Altogether, ZTF18aaqjovh represents another transition event between engine-driven SNe associated with GRBs and "ordinary" Ic-BL SNe., Comment: 18 pages, 8 figures. Accepted to the Astrophysical Journal on 2 March 2020. Revised from proofs on 4 Apr 2020

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium., Comment: 22 pages, 7 figures, 1 table. Accepted for publication in AJ on 12 May 2020