Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

We present the discovery of the optical afterglow and host galaxy of the {\it Swift} short-duration gamma-ray burst, GRB\,181123B. Observations with Gemini-North starting at $\approx 9.1$~hr after the burst reveal a faint optical afterglow with $i\approx25.1$~mag, at an angular offset of 0.59 $\pm$ 0.16$''$ from its host galaxy. Using $grizYJHK$ observations, we measure a photometric redshift of the host galaxy of $z = 1.77^{+0.30}_{-0.17}$. From a combination of Gemini and Keck spectroscopy of the host galaxy spanning 4500-18000~\AA , we detect a single emission line at 13390~\AA, inferred as H$\beta$ at $z = 1.754 \pm 0.001$ and corroborating the photometric redshift. The host galaxy properties of GRB\,181123B are typical to those of other SGRB hosts, with an inferred stellar mass of $\approx 1.7 \times 10^{10}\,M_{\odot}$, mass-weighted age of $\approx 0.9$~Gyr and optical luminosity of $\approx 0.9L^{*}$. At $z=1.754$, GRB\,181123B is the most distant secure SGRB with an optical afterglow detection, and one of only three at $z>1.5$. Motivated by a growing number of high-$z$ SGRBs, we explore the effects of a missing $z>1.5$ SGRB population among the current {\it Swift} sample on delay time distribution models. We find that log-normal models with mean delay times of $\approx 4-6$~Gyr are consistent with the observed distribution, but can be ruled out to $95\%$ confidence with an additional $\approx1-5$~{\it Swift} SGRBs recovered at $z>1.5$. In contrast, power-law models with $\propto$ $t^{-1}$ are consistent with the redshift distribution and can accommodate up to $\approx30$ SGRBs at these redshifts. Under this model, we predict that $\approx 1/3$ of the current {\it Swift} population of SGRBs is at $z>1$. The future discovery or recovery of existing high-$z$ SGRBs will provide significant discriminating power on their delay time distributions, and thus their formation channels.

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$.

We introduce pinta, a pipeline for reducing the 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. Furthermore, the results of a calibration experiment carried out to determine the relative time offsets between different back-end modes and the correct interpretation of the observatory frequency settings at the uGMRT, which are crucial for performing precision pulsar timing, are also discussed.

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). At $z=0.0866$, ZTF19aawfbtg (SN2019hge) has more than 10 optical spectra at phases from $-41$ to $+103$ days. Most of these spectra match well with that of PTF10hgi, suggesting SN2019hge as the second case of a He-rich SLSN-I. Confirmation comes from its near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features, yielding a lower limit to the He-rich fraction of $\sim9$\%. It is suggested that the excitation of He I atoms requires non-thermal radiation. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. However, how the magnetar energy is transported to the outer layers is unclear. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities, around $-20.06\pm0.16$ mag ($g$-band, rest-frame), in comparison to the full range of $-20$ to $-23$ mag covered by SLSN-I.

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.

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)

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.

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.

We present optical, radio, and X-ray observations of SN2020bvc (=ASASSN20bs; ZTF20aalxlis), a nearby ($z=0.0252$; $d$=114 Mpc) broad-lined (BL) Type Ic supernova (SN). Our observations show that SN2020bvc shares several properties in common with the Ic-BL SN2006aj, which was associated with the low-luminosity gamma-ray burst (LLGRB) 060218. First, the 10 GHz radio light curve is on the faint end of LLGRB-SNe ($L_\mathrm{radio} \approx 10^{37}$erg/s): we model our VLA observations (spanning 13-43 d) as synchrotron emission from a mildly relativistic ($v \gtrsim 0.3c$) forward shock. Second, with Swift and Chandra we detect X-ray emission ($L_X \approx 10^{41}$erg/s) that is not naturally explained as inverse Compton emission or as part of the same synchrotron spectrum as the radio emission. Third, high-cadence ($6\times$/night) data from the Zwicky Transient Facility (ZTF) shows a double-peaked optical light curve, the first peak from shock-cooling emission from extended low-mass material (mass $M<10^{-2} M_\odot$ at radius $R>10^{12}$cm) and the second peak from the radioactive decay of Ni-56. SN2020bvc is the first confirmed double-peaked Ic-BL SN discovered without a GRB trigger, and shows X-ray and radio emission similar to LLGRB-SNe: this is consistent with models in which the same mechanism produces both the LLGRB and the shock-cooling emission. For four of the five other nearby ($z\lesssim0.05$) Ic-BL SNe with ZTF high-cadence data, we rule out a first peak like that seen in SN2006aj and SN2020bvc, i.e. that lasts $\approx 1$d and reaches a peak luminosity $M \approx -18$. X-ray and radio follow-up observations of future such events will establish whether double-peaked optical light curves are indeed predictive of LLGRB-like X-ray and radio emission.

We report the results from spectroscopic observations of the multiple images of the strongly lensed Type Ia supernova (SN Ia), iPTF16geu, obtained with ground based telescopes and the Hubble Space Telescope (HST). From a single epoch of slitless spectroscopy with HST, we can resolve spectra of individual lensed supernova images for the first time. This allows us to perform an independent measurement of the time-delay between the two brightest images, $\Delta t = 1.4 \pm 5.0$ days, which is consistent with the time-delay measured from the light-curves. We also present measurements of narrow emission and absorption lines characterizing the interstellar medium in the host galaxy at z=0.4087, as well as in the foreground lensing galaxy at z=0.2163. We detect strong Na ID absorption in the host galaxy, indicating that iPTF16geu belongs to a subclass of SNe Ia displaying "anomalously" large Na ID column densities in comparison to the amount of dust extinction derived from their light curves. For the deflecting galaxy, we refine the measurement of the velocity dispersion, $\sigma = 129 \pm 4$ km/s, which significantly constrains the lens model. Since the time-delay between the SN images is negligible, we can use unresolved ground based spectroscopy, boosted by a factor ~70 from lensing magnification, to study the properties of a high-z SN Ia with unprecedented signal-to-noise ratio. The spectral properties of the supernova, such as pseudo-Equivalent widths of several absorption features and velocities of the Si II-line indicate that iPTF16geu, besides being lensed, is a normal SN Ia, indistinguishable from well-studied ones in the local universe, providing support for the use of SNe Ia in precision cosmology. We do not detect any significant deviations of the SN spectral energy distribution from microlensing of the SN photosphere by stars and compact objects in the lensing galaxy.

(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.

We present ZTF18abvkwla (the "Koala"), a fast blue optical transient discovered in the Zwicky Transient Facility (ZTF) One-Day Cadence (1DC) Survey. ZTF18abvkwla has a number of features in common with the groundbreaking transient AT2018cow: blue colors at peak ($g-r\approx-0.5$ mag), a short rise time from half-max of under two days, a decay time to half-max of only three days, a high optical luminosity ($M_{g,\mathrm{peak}}\approx-20.6$mag), a hot ($\gtrsim 40,000$K) featureless spectrum at peak light, and a luminous radio counterpart. At late times ($\Delta t>80$d) the radio luminosity of ZTF18abvkwla ($\nu L_\nu \gtrsim 10^{40}$erg/s at 10 GHz, observer-frame) is most similar to that of long-duration gamma-ray bursts (GRBs). The host galaxy is a dwarf starburst galaxy ($M\approx5\times10^{8}M_\odot$, $\mathrm{SFR}\approx7 M_\odot$/yr) that is moderately metal-enriched ($\log\mathrm{[O/H]} \approx 8.5$), similar to the hosts of GRBs and superluminous supernovae. As in AT2018cow, the radio and optical emission in ZTF18abvkwla likely arise from two separate components: the radio from fast-moving ejecta ($\Gamma \beta c >0.38c$) and the optical from shock-interaction with confined dense material ($<0.07M_\odot$ in $\sim 10^{15}$cm). Compiling transients in the literature with $t_\mathrm{rise} <5$d and $M_\mathrm{peak}<-20$mag, we find that a significant number are engine-powered, and suggest that the high peak optical luminosity is directly related to the presence of this engine. From 18 months of the 1DC survey, we find that transients in this rise-luminosity phase space are at least two to three orders of magnitude less common than CC SNe. Finally, we discuss strategies for identifying such events with future facilities like the Large Synoptic Survey Telescope, and prospects for detecting accompanying X-ray and radio emission.

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.

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.

We present a systematic study of the most luminous ($M_{\mathrm{IR}}$ [Vega magnitudes] brighter than $-14$) infrared (IR) transients discovered by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) between 2014 and 2018 in nearby galaxies ($D < 35$ Mpc). The sample consists of nine events that span peak IR luminosities of $M_{[4.5],\mathrm{peak}}$ between $-14$ and $-18.2$, show IR colors between $0.2 < ([3.6]{-}[4.5]) < 3.0$, and fade on timescales between $55$ days $< t_{\mathrm{fade}} < 480$ days. The two reddest events ($A_V > 12$) show multiple, luminous IR outbursts over several years and have directly detected, massive progenitors in archival imaging. With analyses of extensive, multiwavelength follow-up, we suggest the following possible classifications: five obscured core-collapse supernovae (CCSNe), two erupting massive stars, one luminous red nova, and one intermediate-luminosity red transient. We define a control sample of all optically discovered transients recovered in SPIRITS galaxies and satisfying the same selection criteria. The control sample consists of eight CCSNe and one Type Iax SN. We find that 7 of the 13 CCSNe in the SPIRITS sample have lower bounds on their extinction of $2 < A_V < 8$. We estimate a nominal fraction of CCSNe in nearby galaxies that are missed by optical surveys as high as $38.5^{+26.0}_{-21.9}$% (90% confidence). This study suggests that a significant fraction of CCSNe may be heavily obscured by dust and therefore undercounted in the census of nearby CCSNe from optical searches.

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.

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 black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three 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 we 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_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star 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 $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising ($1.4\pm0.1$ mag/hr) and luminous ($M_{g,\mathrm{peak}}=-20$ mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity ($L_{\mathrm{bol}} \gtrsim 3 \times 10^{44}$ erg $\mathrm{sec}^{-1}$), the short rise time ($t_{\mathrm{rise}}= 3$ days in $g$-band), and the blue colors at peak ($g-r\sim-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 ($T_{\mathrm{eff}}\gtrsim40,000$ K) spectra of a stripped-envelope SN. A retrospective search revealed luminous ($M_g \sim M_r \approx -14\,$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_\mathrm{\gamma,iso}<4.9 \times 10^{48}$ erg, a limit on X-ray emission $L_{\mathrm{X}} < 10^{40}\,$erg sec$^{-1}$, and a limit on radio emission $\nu L_\nu \lesssim 10^{37}\,$erg sec$^{-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_\odot$) at large radii ($3 \times 10^{14}\,$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.

We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of $M_r = -17.98$, peak bolometric luminosity of ${\sim} 1.4 \times 10^{43}$erg s$^{\mathrm{-1}}$ and a rapid decline rate of $0.48 \pm 0.03$ mag day$^{\textrm{-1}}$ in the $r$ band. The bolometric luminosity evolves too quickly to be explained by pure $^{\mathrm{56}}$Ni heating, necessitating the inclusion of an alternative powering source. Incorporating the spin-down of a magnetized neutron star adequately describes the lightcurve and we estimate a small ejecta mass of $M_\mathrm{ej} = 0.10 \pm 0.05$ $\textrm{M}_{\odot}$. Our spectral modelling suggests the ejecta is composed of intermediate mass elements including O, Si and Mg and trace amounts of Fe-peak elements, which disfavours a binary neutron star merger. We discuss three explosion scenarios for SN2018kzr, given the low ejecta mass, intermediate mass element composition and the high likelihood of additional powering - core collapse of an ultra-stripped progenitor, the accretion induced collapse of a white dwarf and the merger of a white dwarf and neutron star. The requirement for an alternative input energy source favours either the accretion induced collapse with magnetar powering or a white dwarf - neutron star merger with energy from disk wind shocks.

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\,{\rm deg}^2$ field at $\Delta 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 localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength 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 neutron star-black hole 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.

(Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans $\approx 7500$ square degrees of the sky every night to a median 5$\sigma$ depth of $15.7$ AB mag outside the Galactic plane. The survey covers $\approx 15000$ square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with PSF-fit source catalogs and transient candidates identified from subtractions within a median delay of $\approx 4$ hours from the time of observation. The calibrated data products achieve an astrometric accuracy (RMS) of $\approx 0.7$" with respect to Gaia DR2 for sources with S/N $> 10$, and better than $\approx 0.35$" for sources brighter than $\approx 12$ Vega mag. The photometric accuracy (RMS) achieved in the PSF-fit source catalogs is better than $\approx 3$% for sources brighter than $\approx 12$ Vega mag, as calibrated against the 2MASS catalog. With a field of view $\approx 40\times$ larger than any other existing near infrared imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened micro-lensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period.

The Zwicky Transient Facility (ZTF) is performing a three-day cadence survey of the visible Northern sky (~3$\pi$). The transient candidates found in this survey are announced via public alerts. As a supplementary product ZTF is also conducting a large spectroscopic campaign: the ZTF Bright Transient Survey (BTS). The goal of the BTS is to spectroscopically classify all extragalactic transients brighter than 18.5 mag at peak brightness and immediately announce those classifications to the public. Extragalactic discoveries from ZTF are predominantly Supernovae (SNe). The BTS is the largest flux-limited SN survey to date. Here we present a catalog of the761 SNe that were classified during the first nine months of the survey (2018 Apr. 1 to 2018 Dec. 31). The BTS SN catalog contains redshifts based on SN template matching and spectroscopic host galaxy redshifts when available. Based on this data we perform an analysis of the redshift completeness of local galaxy catalogs, dubbed as the Redshift Completeness Fraction (RCF; the number of SN host galaxies with known spectroscopic redshift prior to SN discovery divided by the total number of SN hosts). In total, we identify the host galaxies of 512 Type Ia supernovae, 227 of which have known spectroscopic redshifts, yielding an RCF estimate of $44\% \pm1\%$. We find a steady decrease in the RCF with increasing distance in the local universe. For z<0.05, or ~200 Mpc, we find RCF=0.6, which has important ramifications when searching for multimessenger astronomical events. Prospects for dramatically increasing the RCF are limited to new multi-fiber spectroscopic instruments, or wide-field narrowband surveys. We find that existing galaxy redshift catalogs are only $50\%$ complete at $r\approx16.9$ mag. Pushing this limit several magnitudes deeper will pay huge dividends when searching for electromagnetic counterparts to gravitational wave events.

Context: Type Ibn supernovae (SNe) are rapidly-evolving and luminous (M$_\text{R,peak}$ $\sim-19$) transients interacting with He-rich circumstellar material (CSM). SN 2018bcc, detected by the ZTF shortly after explosion, provides the best constraints on the shape of the rising lightcurve of a fast Type Ibn. Aims: We use the high-quality data set of SN 2018bcc to study observational signatures of the class. Additionally, the powering mechanism of SN 2018bcc offers insights into the debated progenitor connection of Type Ibn SNe. Methods: Well-constrained lightcurve properties obtained from empirical models are compared with the literature. We fit the pseudo-bolometric lightcurve with semi-analytical models powered by radioactive decay and CSM interaction. The line profiles and emissivity of the prominent He I lines are modeled to study the formation of P-Cygni profiles and to estimate CSM properties. Results: SN 2018bcc has a risetime to peak of $5.6^{+0.2}_{-0.1}$ days in the rest frame with a rising shape powerlaw index near 2, and seems to be a typical rapidly-evolving Type Ibn SN. The spectrum lacks signatures of SN-like ejecta and is dominated by over 15 He emission features at 20 days past peak, alongside Ca and Mg, all with V$_{\text{FWHM}} \sim 2000~\text{km}~\text{s}^{-1}$. The luminous and rapidly-evolving lightcurve can be powered by CSM interaction but not by the decay of radioactive $^{56}$Ni. Modeling of the He I lines indicates a dense and optically thick CSM that can explain the P-Cygni profiles. Conclusions: Like other rapidly-evolving Type Ibn SNe, SN 2018bcc is a luminous transient with a rapid rise to peak powered by shock-interaction inside a dense and He-rich CSM. Its spectra do not support the existence of two Type Ibn spectral classes. We also note the remarkable observational match to pulsational pair instability (PPI) SN models.

We present photometry and spectroscopy of four hydrogen-poor luminous supernovae discovered during the two-month science commissioning and early operations of the Zwicky Transient Facility (ZTF) survey. Three of these objects, SN2018bym (ZTF18aapgrxo), SN2018avk (ZTF18aaisyyp) and SN2018bgv (ZTF18aavrmcg) resemble typical SLSN-I spectroscopically, while SN2018don (ZTF18aajqcue) may be an object similar to SN2007bi experiencing considerable host galaxy reddening, or an intrinsically long-lived, luminous and red SN Ic. We analyze the light curves, spectra, and host galaxy properties of these four objects and put them in context of the population of SLSN-I. SN2018bgv stands out as the fastest-rising SLSN-I observed to date, with a rest-frame g-band rise time of just 10 days from explosion to peak -- if it is powered by magnetar spin-down, the implied ejecta mass is only ~1 M$_{\odot}$. SN2018don also displays unusual properties -- in addition to its red colors and comparatively massive host galaxy, the light curve undergoes some of the strongest light curve undulations post-peak seen in a SLSN-I, which we speculate may be due to interaction with circumstellar material. We discuss the promises and challenges of finding SLSNe in large-scale surveys like ZTF given the observed diversity in the population.

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain their progenitor properties. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in 2018. We describe our method to perform forced point spread function (PSF) photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night ($3g+3r$), all of the 127 SNe Ia are detected in both $g$ and $r$ band more than 10\,d (in the rest frame) prior to the epoch of $g$-band maximum light. The redshifts of these objects range from $z=0.0181$ to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14\,d prior to maximum light (in the rest frame), with a subset of 9 objects being detected more than 17\,d before $g$-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample, including one 02cx-like event ZTF18abclfee (SN\,2018crl), one Ia-CSM SN ZTF18aaykjei (SN\,2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN\,2018eul), ZTF18abdpvnd (SN\,2018dvf), ZTF18aawpcel (SN\,2018cir) and ZTF18abddmrf (SN\,2018dsx).

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

Supernova (SN) 2014C is a unique explosion where a seemingly typical hydrogen-poor stripped envelope SN started to interact with a dense, hydrogen-rich circumstellar medium (CSM) a few months after the explosion. The delayed interaction suggests a detached CSM shell, unlike in a typical SN IIn where the CSM is much closer and the interaction commences earlier post-explosion; indicating a different mass loss history. We present near- to mid-infrared observations of SN 2014C from 1-5 years after the explosion, including uncommon 9.7 $\mu$m imaging with COMICS on the Subaru telescope. Spectroscopy shows that the interaction is still ongoing, with the intermediate-width He I 1.083 $\mu$m emission present out to our latest epoch 1639 days post-explosion. The last Spitzer/IRAC photometry at 1920 days post-explosion further confirms ongoing CSM interaction. The 1-10 $\mu$m spectral energy distributions (SEDs) can be explained by a dust model with a mixture of 69% carbonaceous and 31% silicate dust, pointing to a chemically inhomogeneous CSM. The inference of silicate dust is the first among interacting SNe. An SED model with purely carbonaceous CSM dust is possible, but would require more than 0.22 $M_{\odot}$ of dust, which is an order of magnitude larger than what observed in any other SNe, measured in the same way, at this epoch. The light curve beyond 500 days is well fit by an interaction model with a wind-driven CSM and a mass loss rate of $\sim 10^{-3} \, M_{\odot}\,\rm yr^{-1}$, which presents an additional CSM density component exterior to the constant density shell reported previously in the literature. SN 2014C could originate in a binary system, similar to RY Scuti, which would explain the observed chemical and density profile inhomogeneity in the CSM.

We report the discovery of six active galactic nuclei (AGN) caught "turning on" during the first nine months of the Zwicky Transient Facility (ZTF) survey. The host galaxies were classified as LINERs by weak narrow forbidden line emission in their archival SDSS spectra, and detected by ZTF as nuclear transients. In five of the cases, we found via follow-up spectroscopy that they had transformed into broad-line AGN, reminiscent of the changing-look LINER iPTF 16bco. In one case, ZTF18aajupnt/AT2018dyk, follow-up HST UV and ground-based optical spectra revealed the transformation into a narrow-line Seyfert 1 (NLS1) with strong [Fe VII, X, XIV] and He II 4686 coronal lines. Swift monitoring observations of this source reveal bright UV emission that tracks the optical flare, accompanied by a luminous soft X-ray flare that peaks ~60 days later. Spitzer follow-up observations also detect a luminous mid-infrared flare implying a large covering fraction of dust. Archival light curves of the entire sample from CRTS, ATLAS, and ASAS-SN constrain the onset of the optical nuclear flaring from a prolonged quiescent state. Here we present the systematic selection and follow-up of this new class of changing-look LINERs, compare their properties to previously reported changing-look Seyfert galaxies, and conclude that they are a unique class of transients well-suited to test the uncertain physical processes associated with the LINER accretion state.

We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of $M_r=-13$ (${\nu}L_{\nu}=9\times10^6~L_{\odot}$), 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 $\approx400$ 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~\mu$m luminosity of $M_{[4.5]}\approx-12.2$ 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)\approx0.7$ mag and higher intrinsic luminosity of $M_r\approx-14.9$ (${\nu}L_{\nu}=5.3\times10^7~L_{\odot}$) 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.

The discovery of a transient kilonova following the gravitational-wave 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 $\sim$500 Mpc, the binary neutron star horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility followup by the Global Relay of Observatories Watching Transients Happen (GROWTH) 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 neutron star--black hole merger detected via gravitational waves. Starting 7.5hr after S190426c, over 11.28\,hr of observations, we imaged an area of 525\,deg$^2$ ($r$-band) and 437\,deg$^2$ ($z$-band); this was 16.3\% of the total original localization probability and nearly all of the probability density visible from the Southern hemisphere. The machine-learning based pipeline was optimized for fast turnaround, delivering transient candidates for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hours before the end of our observations. Our observations yielded no detection of a bona fide counterpart to $m_z = 22.5$ and $m_r = 22.9$ at the 5$\sigma$ 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.

We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei (AGN). The supergiant candidate counterparts exhibit a bi-modal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "blue" ULXs fall in IRAC colors $[3.6] - [4.5]\sim0.7$ and $[3.6] - [4.5]\sim0.0$, respectively. The mid-IR colors and absolute magnitudes of 4 "red" and 5 "blue" ULXs are consistent with that of supergiant B[e] (sgB[e]) and red supergiant (RSG) stars, respectively. While "blue", RSG-like mid-IR ULX counterparts likely host RSG mass donors, we propose the "red" counterparts are ULXs exhibiting the "B[e] phenomenon'' rather than hosts of sgB[e] mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of $\dot{M}\sim1\times10^{-4}$ M$_\odot$ yr$^{-1}$ in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, $\alpha=-0.19\pm0.1$, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet.

We investigate ZTF18aalrxas, a double-peaked Type IIb core-collapse supernova (SN) discovered during science validation of the Zwicky Transient Facility (ZTF). ZTF18aalrxas was discovered while the optical emission was still rising towards the initial cooling peak (0.7 mag over 2 days). Our observations consist of multi-band (UV, optical) light-curves, and optical spectra spanning from $\approx0.7$ d to $\approx180$ d past the explosion. We use a Monte-Carlo based non-local thermodynamic equilibrium (NLTE) model, that simultanously reproduces both the $\rm ^{56}Ni$ powered bolometric light curve and our nebular spectrum. This model is used to constrain the synthesized radioactive nickel mass (0.17 $\mathrm{M}_{\odot}$) and the total ejecta mass (1.7 $\mathrm{M}_{\odot}$) of the SN. The cooling emission is modeled using semi-analytical extended envelope models to constrain the progenitor radius ($790-1050$ $\mathrm{R}_{\odot}$) at the time of explosion. Our nebular spectrum shows signs of interaction with a dense circumstellar medium (CSM), and this spetrum is modeled and analysed to constrain the amount of ejected oxygen ($0.3-0.5$ $\mathrm{M}_{\odot}$) and the total hydrogen mass ($\approx0.15$ $\mathrm{M}_{\odot}$) in the envelope of the progenitor. The oxygen mass of ZTF18aalrxas is consistent with a low ($12-13$ $\mathrm{M}_{\odot}$) Zero Age Main Sequence mass progenitor. The light curves and spectra of ZTF18aalrxas are not consistent with massive single star SN Type IIb progenitor models. The presence of an extended hydrogen envelope of low mass, the presence of a dense CSM, the derived ejecta mass, and the late-time oxygen emission can all be explained in a binary model scenario.

Opening up the dynamic infrared sky for systematic time-domain exploration would yield many scientific advances. Multi-messenger pursuits such as localizing gravitational waves from neutron star mergers and quantifying the nucleosynthetic yields require the infrared. Another multi-messenger endeavor that needs infrared surveyors is the study of the much-awaited supernova in our own Milky Way. Understanding shocks in novae, true rates of supernovae and stellar mergers are some other examples of stellar evolution and high energy physics wherein the answers are buried in the infrared. We discuss some of the challenges in the infrared and pathfinders to overcome them. We conclude with recommendations on both infrared discovery engines and infrared follow-up machines that would enable this field to flourish in the next decade.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

There is significant interest in the models for production of short gamma-ray bursts. Until now, the number of known short gamma-ray bursts with multi-wavelength afterglows has been small. While the {\it Fermi} Gamma-Ray Burst Monitor detects many gamma-ray bursts relative to the Neil Gehrels {\it Swift} Observatory, the large localization regions makes the search for counterparts difficult. With the Zwicky Transient Facility recently achieving first light, it is now fruitful to use its combination of depth ($m_\textrm{AB} \sim 20.6$), field of view ($\approx$ 47 square degrees), and survey cadence (every $\sim 3$ days) to perform Target of Opportunity observations. We demonstrate this capability on GRB 180523B, which was recently announced by the {\it Fermi} Gamma-Ray Burst Monitor as a short gamma-ray burst. ZTF imaged $\approx$ 2900\,square degrees of the localization region, resulting in the coverage of 61.6\,\% of the enclosed probability over 2 nights to a depth of $m_\textrm{AB} \sim 20.5$. We characterized 14 previously unidentified transients, and none were found to be consistent with a short gamma-ray burst counterpart. This search with the Zwicky Transient Facility shows it is an efficient camera for searching for coarsely-localized short gamma-ray burst and gravitational-wave counterparts, allowing for a sensitive search with minimal interruption to its nominal cadence.

We present infrared (IR) photometry and spectroscopy of the Type II-P SN 2017eaw and its progenitor in the nearby galaxy NGC 6946. Progenitor observations in the Ks band in 4 epochs from 1 year to 1 day before the explosion reveal no significant variability in the progenitor star greater than 6% that last longer than 200 days. SN 2017eaw is a typical SN II-P with near-IR and mid-IR photometric evolution similar to those of SNe 2002hh and 2004et, other normal SNe II-P in the same galaxy. Spectroscopic monitoring between 389 and 480 days post explosion reveals strong CO first overtone emission at 389 d, with a line profile matching that of SN 1987A from the same epoch, indicating $\sim 10^{-3} \, M_{\odot}$ of CO at 1,800 K. From the 389 d epoch until the most recent observation at 566 d, the first overtone feature fades while the 4.5 $\mu$m excess, likely from the CO fundamental band, remains. This behavior indicates that the CO has not been destroyed, but that the gas has cooled enough that the levels responsible for first overtone emissions are no longer populated. Finally, the evolution of Spitzer 3.6 $\mu$m photometry shows evidence for dust formation in SN 2017eaw, with a dust mass of $10^{-6}$ or $10^{-4}\,M_{\odot}$ assuming carbonaceous or silicate grains respectively.

The detonation of a helium shell on a white dwarf has been proposed as a possible explosion triggering mechanism for Type Ia supernovae. Here, we report ZTF 18aaqeasu (SN 2018byg/ATLAS 18pqq), a peculiar Type I supernova, consistent with being a helium-shell double-detonation. With a rise time of $\approx 18$ days from explosion, the transient reached a peak absolute magnitude of $M_R \approx -18.2$ mag, exhibiting a light curve akin to sub-luminous SN 1991bg-like Type Ia supernovae, albeit with an unusually steep increase in brightness within a week from explosion. Spectra taken near peak light exhibit prominent Si absorption features together with an unusually red color ($g-r \approx 2$ mag) arising from nearly complete line blanketing of flux blue-wards of 5000 \AA. This behavior is unlike any previously observed thermonuclear transient. Nebular phase spectra taken at and after $\approx 30$ days from peak light reveal evidence of a thermonuclear detonation event dominated by Fe-group nucleosynthesis. We show that the peculiar properties of ZTF 18aaqeasu are consistent with the detonation of a massive ($\approx 0.15$ M$_\odot$) helium shell on a sub-Chandrasekhar mass ($\approx 0.75$ M$_\odot$) white dwarf after including mixing of $\approx 0.2$ M$_\odot$ of material in the outer ejecta. These observations provide evidence of a likely rare class of thermonuclear supernovae arising from detonations of massive helium shells.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

Compact neutron star binary systems are produced from binary massive stars through stellar evolution involving up to two supernova explosions. The final stages in the formation of these systems have not been directly observed. We report the discovery of iPTF 14gqr (SN 2014ft), a Type Ic supernova with a fast evolving light curve indicating an extremely low ejecta mass ($\approx 0.2$ solar masses) and low kinetic energy ($\approx 2 \times 10^{50}$ ergs). Early photometry and spectroscopy reveal evidence of shock cooling of an extended He-rich envelope, likely ejected in an intense pre-explosion mass loss episode of the progenitor. Taken together, we interpret iPTF 14gqr as evidence for ultra-stripped supernovae that form neutron stars in compact binary systems.

Calcium rich gap transients represent an intriguing new class of faint and fast evolving supernovae that exhibit strong [Ca II] emission in their nebular phase spectra. In this paper, we present the discovery and follow-up observations of iPTF 16hgs -- an intermediate luminosity and fast evolving transient that exhibited a double peaked light curve. Exhibiting a typical Type Ib spectrum in the photospheric phase and an early transition to a [Ca II] dominated nebular phase, we show that iPTF 16hgs shows properties consistent with the class of Ca-rich gap transients, with two interesting exceptions. First, while the second peak of the light curve is similar to other Ca-rich gap transients (suggesting $M_{ej}$ of 0.4 M$_\odot$ and peak luminosity of $3 \times 10^{41}$ ergs s$^{-1}$), we show that the first blue and fast declining (over $2$ days) peak is unique to this source. Second, with Integral Field Unit observations of the host galaxy, we find that iPTF 16hgs occurred in the outskirts (projected offset of $6$ kpc $ = 1.9 R_{eff}$) of a low metallicity (0.4 Z$_\odot$), star forming, dwarf spiral galaxy. Using deep late-time VLA and uGMRT observations, we place stringent limits on the local environment of the source, ruling out a large parameter space of circumstellar densities and mass loss environments of the progenitor. If iPTF 16hgs shares explosion physics with the class of Ca-rich gap transients, the presence of the first peak can be explained by enhanced mixing of 0.01 M$_\odot$ of $^{56}$Ni into the outer layers the ejecta, reminiscent of some models of He-shell detonations on WDs. On the other hand, if iPTF 16hgs is physically unrelated to the class, the first peak is consistent with shock cooling emission (of an envelope with a mass of 0.08 M$_\odot$ and radius of 13 R$_\odot$) associated with a core-collapse explosion of a highly stripped massive star in a close binary system.

We present the discovery by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) of a likely supernova (SN) in NGC 3556 at only 8.8 Mpc, which, despite its proximity, was not detected by optical searches. A luminous infrared (IR) transient at $M_{[4.5]} = -16.7$ mag (Vega), SPIRITS 16tn is coincident with a dust lane in the inclined, star-forming disk of the host. Using IR, optical, and radio observations, we attempt to determine the nature of this event. We estimate $A_V \approx$ 8 - 9 mag of extinction, placing it among the three most highly obscured IR-discovered SNe to date. The [4.5] light curve declined at a rate of 0.013 mag day$^{-1}$, and the $[3.6] - [4.5]$ color grew redder from 0.7 to $\gtrsim$ 1.0 mag by 184.7 days post discovery. Optical/IR spectroscopy shows a red continuum, but no clearly discernible features, preventing a definitive spectroscopic classification. Deep radio observations constrain the radio luminosity of SPIRITS 16tn to $L_{\nu} \lesssim 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ between 3 - 15 GHz, excluding many varieties of radio core-collapse SNe. A type Ia SN is ruled out by the observed red IR color, and lack of features normally attributed to Fe-peak elements in the optical and IR spectra. SPIRITS 16tn was fainter at [4.5] than typical stripped-envelope SNe by $\approx$ 1 mag. Comparison of the spectral energy distribution to SNe II suggests SPIRITS 16tn was both highly obscured, and intrinsically dim, possibly akin to the low-luminosity SN 2005cs. We infer the presence of an IR dust echo powered by a peak luminosity of the transient of $5 \times 10^{40}$ erg s$^{-1} < L_{\mathrm{peak}} < 4\times10^{43}$ erg s$^{-1}$, consistent with the observed range for SNe II. This discovery illustrates the power of IR surveys to overcome the compounding effects of visible extinction and optically sub-luminous events in completing the inventory of nearby SNe.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

Gravitational waves have been detected from a binary neutron star merger event, GW170817. The detection of electromagnetic radiation from the same source has shown that the merger occurred in the outskirts of the galaxy NGC 4993, at a distance of 40 megaparsecs from Earth. We report the detection of a counterpart radio source that appears 16 days after the event, allowing us to diagnose the energetics and environment of the merger. The observed radio emission can be explained by either a collimated ultra-relativistic jet viewed off-axis, or a cocoon of mildly relativistic ejecta. Within 100 days of the merger, the radio light curves will distinguish between these models and very long baseline interferometry will have the capability to directly measure the angular velocity and geometry of the debris.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

AX J1745.6-2901 is a high-inclination (eclipsing) transient neutron star (NS) Low Mass X-ray Binary (LMXB) showcasing intense ionised Fe K absorption. We present here the analysis of 11 XMM-Newton and 15 NuSTAR new data-sets (obtained between 2013-2016), therefore tripling the number of observations of AX J1745.6-2901 in outburst. Thanks to simultaneous XMM-Newton and NuSTAR spectra, we greatly improve on the fitting of the X-ray continuum. During the soft state the emission can be described by a disk black body ($kT\sim1.1-1.2$ keV and inner disc radius $r_{DBB}\sim14$ km), plus hot ($kT\sim2.2-3.0$ keV) black body radiation with a small emitting radius ($r_{BB}\sim0.5-0.8$ km) likely associated with the boundary layer or NS surface, plus a faint Comptonisation component. Imprinted on the spectra are clear absorption features created by both neutral and ionised matter. Additionally, positive residuals suggestive of an emission Fe K$\alpha$ disc line and consistent with relativistic ionised reflection are present during the soft state, while such residuals are not significant during the hard state. The hard state spectra are characterised by a hard ($\Gamma\sim1.9-2.1$) power law, showing no evidence for a high energy cut off ($kT_e>60-140$ keV) and implying a small optical depth ($\tau<1.6$). The new observations confirm the previously witnessed trend of exhibiting strong Fe K absorption in the soft state, that significantly weakens during the hard state. Optical (GROND) and radio (GMRT) observations suggest for AX J1745.6-2901 a standard broad band SED as typically observed in accreting neutron stars.

Microstructure emission, involving short time scale, often quasi-periodic, intensity fluctuations in subpulse emission, is well known in normal period pulsars. In this letter, we present the first detections of quasi-periodic microstructure emission from millisecond pulsars (MSPs), from Giant Metrewave Radio Telescope (GMRT) observations of two MSPs at 325 and 610 MHz. Similar to the characteristics of microstructure observed in normal period pulsars, we find that these features are often highly polarized, and exhibit quasi-periodic behavior on top of broader subpulse emission, with periods of the order of a few $\mu$s. By measuring their widths and periodicities from single pulse intensity profiles and their autocorrelation functions, we extend the microstructure timescale - rotation period relationship by more than an order of magnitude down to rotation periods $\sim$ 5 ms, and find it to be consistent with the relationship derived earlier for normal pulsars. The similarity of behavior is remarkable, given the significantly different physical properties of MSPs and normal period pulsars, and rules out several previous speculations about the possible different characteristics of microstructure in MSP radio emission. We discuss the possible reasons for the non-detection of these features in previous high time resolution MSP studies along with the physical implications of our results, both in terms of a geometric beam sweeping model and temporal modulation model for micropulse production.

The orbital period evolution of X-ray binaries provides fundamental clues to understanding mechanisms of angular momentum loss from these systems. We present an X-ray eclipse timing analysis of the transient low mass X-ray binary AX J1745.6-2901. This system shows full eclipses and thus is one of the few objects for which accurate orbital evolution studies using this method can be carried out. We report on XMM-Newton and ASCA observations covering 30 complete X-ray eclipses spanning an interval of more than 20 years. We improve the determination of the orbital period to a relative precision of $2\times10^{-8}$, two orders of magnitudes better than previous estimates. We determine, for the first time, a highly significant rate of decrease of the orbital period $\dot{P}_{orb}=-4.03\pm0.32\times10^{-11}$~s/s. This is at least one order of magnitude larger than expected from conservative mass transfer and angular momentum losses due to gravitational waves and magnetic breaking, and might result from non-conservative mass transfer. Imprinted on the long term evolution of the orbit, we observe highly significant eclipse leads-delays of ~10-20 s, characterised by a clear state dependence in which, on average, eclipses occur earlier during the hard state.

The advent of the new generation of X-ray telescopes yielded a significant step forward in our understanding of ionised absorption generated in the accretion discs of X-ray binaries. It has become evident that these relatively weak and narrow absorption features, sporadically present in the X-ray spectra of some systems, are actually the signature of equatorial outflows, which might carry away more matter than that being accreted. Therefore, they play a major role in the accretion phenomenon. These outflows (or ionised atmospheres) are ubiquitous during the softer states but absent during the power-law dominated, hard states, suggesting a strong link with the state of the inner accretion disc, presence of the radio-jet and the properties of the central source. Here, we discuss the current understanding of this field.

A fully real-time coherent dedispersion system has been developed for the pulsar back-end at the Giant Metrewave Radio Telescope (GMRT). The dedispersion pipeline uses the single phased array voltage beam produced by the existing GMRT software back-end (GSB) to produce coherently dedispersed intensity output in real time, for the currently operational bandwidths of 16 MHz and 32 MHz. Provision has also been made to coherently dedisperse voltage beam data from observations recorded on disk. We discuss the design and implementation of the real-time coherent dedispersion system, describing the steps carried out to optimise the performance of the pipeline. Presently functioning on an Intel Xeon X5550 CPU equipped with a NVIDIA Tesla C2075 GPU, the pipeline allows dispersion free, high time resolution data to be obtained in real-time. We illustrate the significant improvements over the existing incoherent dedispersion system at the GMRT, and present some preliminary results obtained from studies of pulsars using this system, demonstrating its potential as a useful tool for low frequency pulsar observations. We describe the salient features of our implementation, comparing it with other recently developed real-time coherent dedispersion systems. This implementation of a real-time coherent dedispersion pipeline for a large, low frequency array instrument like the GMRT, will enable long-term observing programs using coherent dedispersion to be carried out routinely at the observatory. We also outline the possible improvements for such a pipeline, including prospects for the upgraded GMRT which will have bandwidths about ten times larger than at present.