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

We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $\theta \sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $\rho \sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.

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)

Many inspiraling and merging stellar remnants emit both gravitational and electromagnetic radiation as they orbit or collide. These gravitational wave events together with their associated electromagnetic counterparts provide insight about the nature of the merger, allowing us to further constrain parameters about the properties of the binary. With the future launch of the Laser Interferometer Space Antenna (LISA), follow up observations and models are needed of ultracompact binary (UCB) systems. Current and upcoming long baseline time domain surveys will observe many of these UCBs. We present a new fast periodic object search tool based on the Conditional Entropy algorithm. This new implementation of Conditional Entropy allows for a grid search over both the period ($P$) and the time derivative of the period ($\dot{P}$). To demonstrate the performance and usage of this tool, we use a galactic population of UCBs generated from the population synthesis code COSMIC, as well as a Curated catalog for varying periods at fixed intrinsic parameters. We simulate light curves as likely to be observed by future time domain surveys by using an existing eclipsing binary light curve model accounting for the change in orbital period due to gravitational radiation. We find that a search with $\dot{P}$ values is necessary for detecting binaries at orbital periods less than $\sim$10 min. We also show it is useful in finding and characterizing binaries with longer periods, but at a higher computational cost. Our code is called gce (GPU-Accelerated Conditional Entropy). It is available at https://github.com/mikekatz04/gce.

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

The Zwicky Transient Facility has begun to discover binary systems with orbital periods that are less than 1\,hr. Combined with dedicated follow-up systems, which allow for high-cadence photometry of these sources, systematic confirmation and characterization of these sources are now possible. Here, we report the discovery of ZTF J190125.42+530929.5, a 40.6\,min orbital period, eclipsing double white-dwarf binary. Both photometric and spectroscopic modeling confirm its nature, yielding an estimated inclination of $i = 86.2^{+0.6}_{-0.2}\,\rm degrees$ and primary and secondary effective temperatures of $\textrm{T}_\textrm{eff} = 28,000^{+500}_{-500}\,K$ and $\textrm{T}_\textrm{eff} = 17,600^{+400}_{-400}\,K$ respectively. This system adds to a growing list of sources for future gravitational-wave detectors and contributes to the demographic analysis of double degenerates.

4FGL J0935.3+0901 is a gamma-ray source detected by the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope. We have conducted detailed analysis of the LAT data for this source and multi-wavelength studies of the source field. Its gamma-ray emission can be described with a power law ($\Gamma=2.0\pm0.2$) with an exponential cutoff ($E_c=2.9\pm1.6$ GeV), while the flux shows significant long-term variations. From analysis of archival Neil Gehrels Swift X-ray Telescope data, we find only one X-ray source in the LAT's 2$\sigma$ error region. Within a 3.7" radius error circle of the X-ray source, there is only one optical object down to $r'\sim$23 mag. Time-resolved photometry of the optical object indicates a likely 2.5~hr periodic modulation, while its spectrum shows double-peaked hydrogen and helium emission lines (similar to those seen in accretion discs in low-mass X-ray binaries). Combining these results, we conclude that we have discovered a compact X-ray emitting binary in likely association with 4FGL J0935.3+0901, i.e., a millisecond pulsar (MSP) binary. We discuss the implication of the optical spectral features: this binary could be a transitional MSP system at a sub-luminous disc state, although the other possibility, the binary in a rotation-powered state showing the optical emission lines due to intrabinary interaction processes, can not be excluded. Further observational studies will help determine detailed properties of this candidate MSP binary and thus clarify its current state.

Using selection criteria based on amplitude, time and color, we have identified 329 objects as known or candidate cataclysmic variable (CVs) during the first year of testing and operation of the Zwicky Transient Facility (ZTF). Of these, 90 are previously confirmed CVs, 218 are strong candidates based on the shape and color of their light curves obtained during 3-562 days of observations, and the remaining 21 are possible CVs but with too few data points to be listed as good candidates. Almost half the strong candidates are within 10 deg of the galactic plane, in contrast to most other large surveys which have avoided crowded fields. The available Gaia parallaxes are consistent with sampling the low mass transfer CVs, as predicted by population models. Our followup spectra have confirmed Balmer/helium emission lines in 27 objects, with four showing high excitation HeII emission, including candidates for an AM CVn, a polar and an intermediate polar. Our results demonstrate that a complete survey of the galactic plane is needed to accomplish an accurate determination of the number of CVs existing in the Milky Way.

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star.

The main objective of the present work is to extend these investigations by computing the gravity and limb-darkening coefficients for white dwarf atmosphere models with hydrogen, helium, or mixed compositions (types DA, DB, and DBA). We computed gravity and limb-darkening coefficients for DA, DB, and DBA white dwarfs atmosphere models, covering the transmission curves of the Sloan, UBVRI, Kepler, TESS, and Gaia photometric systems. Specific calculations for the HiPERCAM instrument were also carried out. For all calculations of the limb-darkening coefficients we used the least-squares method. Concerning the effects of tidal and rotational distortions, we also computed for the first time the gravity-darkening coefficients $y(\lambda)$ for white dwarfs using the same models of stellar atmospheres as in the case of limb-darkening. A more general differential equation was introduced to derive these quantities, including the partial derivative $\left(\partial{\ln I_o(\lambda)}/{\partial{\ln g}}\right)_{T_{\rm eff}}$. Six laws were adopted to describe the specific intensity distribution: linear, quadratic, square root, logarithmic, power-2, and a more general one with four coefficients. The computations are presented for the chemical compositions log[H/He] = $-$10.0 (DB), $-$2.0 (DBA) and He/H = 0 (DA), with log g varying between 5.0 and 9.5 and effective temperatures between 3750 K-100\,000 K. For effective temperatures higher than 40\,000 K, the models were also computed adopting nonlocal thermal equilibirum (DA). The adopted mixing-length parameters are ML2/$\alpha = 0.8$ (DA case) and 1.25 (DB and DBA). The results are presented in the form of 112 tables. Additional calculations, such as for other photometric systems and/or different values of log[H/He], $\log g,$ and T$_{\rm eff}$ can be performed upon request.

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 report the discovery of a detached double white dwarf binary with an orbital period of $\approx20.6\,\rm minutes$, PTF J053332.05+020911.6. The visible object in this binary, PTF J0533+0209B, is a $\approx 0.17\,M_\odot$ mass white dwarf with a helium-dominated atmosphere containing traces of hydrogen (DBA). This object exhibits ellipsoidal variations due to tidal deformation, and is the visible component in a single-lined spectroscopic binary with a velocity semi-amplitude of $K_B=618.7\pm6.9 \, \rm km \, s^{-1}$. We have detected significant orbital decay due to the emission of gravitational radiation, and we expect that the Laser Interferometer Space Antenna (LISA) will detect this system with a signal to noise of $8.4^{+4.2}_{-3.0}$ after four years of operation. Because this system already has a well determined orbital period, radial velocity semi-amplitude, temperature, atmospheric composition, surface gravity, and orbital decay rate, a LISA signal will help fully constrain the properties of this system by providing a direct measurement of its inclination. Thus, this binary demonstrates the synergy between electromagnetic and gravitational radiation for constraining the physical properties of an astrophysical object.

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

General relativity predicts that short orbital period binaries emit significant gravitational radiation, and the upcoming Laser Interferometer Space Antenna (LISA) is expected to detect tens of thousands of such systems; however, few have been identified, and only one is eclipsing--the double white dwarf binary SDSS J065133.338+284423.37, which has an orbital period of 12.75 minutes. Here, we report the discovery of an eclipsing double white dwarf binary system with an orbital period of only 6.91 minutes, ZTF J153932.16+502738.8. This system has an orbital period close to half that of SDSS J065133.338+284423.37 and an orbit so compact that the entire binary could fit within the diameter of the planet Saturn. The system exhibits a deep eclipse, and a double-lined spectroscopic nature. We observe rapid orbital decay, consistent with that expected from general relativity. ZTF J153932.16+502738.8 is a significant source of gravitational radiation close to the peak of LISA's sensitivity, and should be detected within the first week of LISA observations.

Using high-cadence observations from the Zwicky Transient Facility at low Galactic latitudes, we have discovered a new class of pulsating, hot, compact stars. We have found four candidates, exhibiting blue colors ($g-r\leq-0.1$ mag), pulsation amplitudes of $>5\%$, and pulsation periods of $200 - 475$ sec. Fourier transforms of the lightcurves show only one dominant frequency. Phase-resolved spectroscopy for three objects reveals significant radial velocity, T$_{\rm eff}$ and log(g) variations over the pulsation cycle, consistent with large amplitude radial oscillations. The mean T$_{\rm eff}$ and log(g) for these stars are consistent with hot subdwarf B (sdB) effective temperatures and surface gravities. We calculate evolutionary tracks using MESA and adiabatic pulsations using GYRE for low-mass helium-core pre-white dwarfs and low mass helium-burning stars. Comparison of low-order radial oscillation mode periods with the observed pulsation periods show better agreement with the pre-white dwarf models. Therefore, we suggest that these new pulsators and Blue Large-Amplitude Pulsators (BLAPs) could be members of the same class of pulsators, composed of young $\approx0.25-0.35$ M$_\odot$ helium-core pre-white dwarfs.

We present a detailed study of a transient in the center of SDSS1115+0544 based on the extensive UV, optical, mid-IR light curves (LC) and spectra over 1200 days. The host galaxy is a quiescent early type galaxy at $z$ = 0.0899 with a blackhole mass of $2\times10^7M_\odot$. The transient underwent a 2.5 magnitude brightening over $\sim120$ days, reaching a peak $V$-band luminosity (extinction corrected) of $-20.9$ magnitude, then fading 0.5 magnitude over 200 days, settling into a plateau of $>600$ days. Following the optical brightening are the significant mid-IR flares at $3.4$ and $4.5\mu$m, with a peak time delay of $\sim180$ days. The mid-IR LCs are explained as the echo of UV photons by a dust medium with a radius of $5\times10^{17}$ cm, consistent with $\rm E(B-V)$ of 0.58 inferred from the spectra. This event is very energetic with an extinction corrected $L_{bol} \sim 4\times10^{44}$ erg s$^{-1}$. Optical spectra over 400 days in the plateau phase revealed newly formed broad H$\alpha, \beta$ emission with a FWHM of $\sim3750$ km s$^{-1}$ and narrow coronal lines such as [Fe VII], [Ne V]. This flare also has a steeply rising UV continuum, detected by multi-epoch $Swift$ data at $+700$ to $+900$ days post optical peak. The broad Balmer lines and the UV continuum do not show significant temporal variations. The slow evolving LCs over 1200 days, the constant Balmer lines and UV continuum at late-times rule out TDE and SN IIn as the physical model for this event. We propose that this event is a `turn-on' AGN, transitioning from a quiescent state to a type 1 AGN with a sub-Eddington accretion rate of $0.017M_\odot$/yr. This change occurred on a very short time scale of $\sim 120- 200$ days. The discovery of such a rapid `turn-on' AGN poses challenges to accretion disk theories and may indicate such event is not extremely rare.

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The realtime pipeline employs a novel image-differencing algorithm, optimized for the detection of point source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae, rare flux transients, variable stars, eclipsing binaries, variability from active galactic nuclei, counterparts to gravitational wave sources, a more complete census of Type Ia supernovae, and Solar System objects.

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 is a large optical survey in multiple filters producing hundreds of thousands of transient alerts per night. We describe here various machine learning (ML) implementations and plans to make the maximal use of the large data set by taking advantage of the temporal nature of the data, and further combining it with other data sets. We start with the initial steps of separating bogus candidates from real ones, separating stars and galaxies, and go on to the classification of real objects into various classes. Besides the usual methods (e.g., based on features extracted from light curves) we also describe early plans for alternate methods including the use of domain adaptation, and deep learning. In a similar fashion we describe efforts to detect fast moving asteroids. We also describe the use of the Zooniverse platform for helping with classifications through the creation of training samples, and active learning. Finally we mention the synergistic aspects of ZTF and LSST from the ML perspective.

The Kitt Peak Electron Multiplying CCD (EMCCD) demonstrator is a new instrument that has been developed for use at the Kitt Peak National Observatory's 84-inch telescope. The EMCCD enables single-band optical imaging in the Sloan \textit{g} and \textit{r} bands and Johnson \textit{UVRI} filters. The EMCCD is chosen for its sub-electron effective read noise using large multiplicative gains. As a result, frame rates of greater than 1\,Hz are possible. The field-of-view is 4.4$^\prime$ $\times$ 4.4$^\prime$ and the pixel size is 0.259$^{\prime\prime}$. This camera, coupled with a fully roboticized telescope, is ideal for follow-up of short period, white dwarf binary candidates, as well as short duration transient and periodic sources identified by large field-of-view all-sky surveys such as the Zwicky Transient Facility.

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.

Although dwarf carbon (dC) stars are thought universally to be binaries to explain the presence of $C_2$ in their spectra while still near main sequence luminosity, direct observational evidence for binarity is remarkably scarce. Here we report the detection of a 2.92 d periodicity in both photometry and radial velocity of SDSS J125017.90+252427.6, an $r=16.4$ dC star. This is the first photometric binary dC, and only the second dC spectroscopic binary. The relative phase of the photometric period to the spectroscopic observations suggests that the photometric variations are a reflection effect due to heating from an unseen companion. The observed radial velocity amplitude of the dC component ($K = 98.8\pm10.7$ km/s) is consistent with a white dwarf companion, presumably the evolved star that earlier donated the carbon to the dC, although substantial orbital evolution must have occurred. Large synoptic photometric surveys such as the Palomar Transient Factory, used for this work, may prove useful for identifying binaries among the shorter period dC stars.

We report the discovery and analysis of 36 new eclipsing EL CVn-type binaries, consisting of a core helium-composition pre-white dwarf and an early-type main-sequence companion, more than doubling the known population of these systems. We have used supervised machine learning methods to search 0.8 million lightcurves from the Palomar Transient Factory, combined with SDSS, Pan-STARRS and 2MASS colours. The new systems range in orbital periods from 0.46-3.8 d and in apparent brightness from ~14-16 mag in the PTF $R$ or $g^{\prime}$ filters. For twelve of the systems, we obtained radial velocity curves with the Intermediate Dispersion Spectrograph at the Isaac Newton Telescope. We modelled the lightcurves, radial velocity curves and spectral energy distributions to determine the system parameters. The radii (0.3-0.7 $\mathrm{R_{\odot}}$) and effective temperatures (8000-17000 K) of the pre-He-WDs are consistent with stellar evolution models, but the masses (0.12-0.28 $\mathrm{M_{\odot}}$) show more variance than models predicted. This study shows that using machine learning techniques on large synoptic survey data is a powerful way to discover substantial samples of binary systems in short-lived evolutionary stages.

The measurement of the direction of WIMP-induced nuclear recoils is a compelling but technologically challenging strategy to provide an unambiguous signature of the detection of Galactic dark matter. Most directional detectors aim to reconstruct the dark-matter-induced nuclear recoil tracks, either in gas or solid targets. The main challenge with directional detection is the need for high spatial resolution over large volumes, which puts strong requirements on the readout technologies. In this paper we review the various detector readout technologies used by directional detectors. In particular, we summarize the challenges, advantages and drawbacks of each approach, and discuss future prospects for these technologies.