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.

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

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

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.

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

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.

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 Zwicky Transient Facility (ZTF) is currently surveying the entire northern sky, including dense Galactic plane fields. Here, we present preliminary results of the search for gravitational microlensing events in the ZTF data collected from the beginning of the survey (March 20, 2018) through June 30, 2019.

We analyze observations from Robo-AO's field M dwarf survey taken on the 2.1m Kitt Peak telescope and perform a multiplicity comparison with Gaia DR2. Through its laser-guided, automated system, the Robo-AO instrument has yielded the largest adaptive optics M dwarf multiplicity survey to date. After developing an interface to visually identify and locate stellar companions, we selected eleven lowsignificance Robo-AO detections for follow-up on the Keck II telescope using NIRC2. In the Robo-AO survey we find 553 candidate companions within 4" around 534 stars out of 5566 unique targets, most of which are new discoveries. Using a position cross match with DR2 on all targets, we assess the binary recoverability of Gaia DR2 and compare the properties of multiples resolved by both Robo-AO and Gaia. The catalog of nearby M dwarf systems and their basic properties presented here can assist other surveys which observe these stars, such as the NASA TESS mission.

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.

Comets are primitive objects that formed in the protoplanetary disk, and have been largely preserved over the history of the Solar System. However, they are not pristine, and surfaces of cometary nuclei do evolve. In order to understand the extent of their primitive nature, we must define the mechanisms that affect their surfaces and comae. We examine the lightcurve of comet 240P/NEAT over three consecutive orbits, and investigate three events of significant brightening ($\Delta m \sim -2$ mag). Unlike typical cometary outbursts, each of the three events are long-lived, with enhanced activity for at least 3 to 6 months. The third event, observed by the Zwicky Transient Facility, occurred in at least two stages. The anomalous behavior appears to have started after the comet was perturbed by Jupiter in 2007, reducing its perihelion distance from 2.53 to 2.12 au. We suggest that the brightening events are temporary transitions to a higher baseline activity level, brought on by the increased insolation, which has warmed previously insulated sub-surface layers. The new activity is isolated to one or two locations on the nucleus, indicating that the surface or immediate sub-surface is heterogeneous. Further study of this phenomenon may provide insight into cometary outbursts, the structure of the near-surface nucleus, and cometary nucleus mantling.

The Kepler space telescope observed over 15,000 stars for asteroseismic studies. Of these, 75% of dwarfs (and 8% of giants) were found to show anomalous behavior: such as suppressed oscillations (low amplitude) or no oscillations at all. The lack of solar-like oscillations may be a consequence of multiplicity, due to physical interactions with spectroscopic companions or due to the dilution of oscillation amplitudes from "wide'" (AO detected; visual) or spectroscopic companions introducing contaminating flux. We present a search for stellar companions to 327 of the Kepler asteroseismic sample, which were expected to display solar-like oscillations. We used direct imaging with Robo-AO, which can resolve secondary sources at ~0."15, and followed up detected companions with Keck AO. Directly imaged companion systems with both separations of $\leq$0."5 and amplitude dilutions >10% all have anomalous primaries, suggesting these oscillation signals are diluted by a sufficient amount of excess flux. We also used the high-resolution spectrometer ESPaDOnS at CFHT to search for spectroscopic binaries. We find tentative evidence for a higher fraction of spectroscopic binaries with high radial velocity scatter in anomalous systems, which would be consistent with previous results suggesting that oscillations are suppressed by tidal interactions in close eclipsing binaries.

The Zwicky Transient Facility (ZTF) is conducting a nightly public survey of all 13 TESS northern sectors in 2019-2020. ZTF will observe the portions of the current TESS sectors visible from Palomar Observatory each night. Each ZTF pointing will have one exposure each with $g$ and $r$ filters, totaling two images per night. ZTF is also making additional nightly $g$- and $r$-band observations of denser stellar regions (e.g. near the Galactic Plane) to better facilitate variability studies of Galactic objects. The limiting magnitude of the ZTF observations is $\approx$20.6 and ZTF saturates at magnitude $\approx$13. ZTF will release data from TESS fields in three forms: nightly alerts distributed by established ZTF brokers, nightly alerts converted to JSON format are distributed via ZTF's bucket on Google Cloud as a tarball, and monthly photometric light curves also distributed via Google Cloud.

We present DeepStreaks, a convolutional-neural-network, deep-learning system designed to efficiently identify streaking fast-moving near-Earth objects that are detected in the data of the Zwicky Transient Facility (ZTF), a wide-field, time-domain survey using a dedicated 47 sq. deg camera attached to the Samuel Oschin 48-inch Telescope at the Palomar Observatory in California, United States. The system demonstrates a 96-98% true positive rate, depending on the night, while keeping the false positive rate below 1%. The sensitivity of DeepStreaks is quantified by the performance on the test data sets as well as using known near-Earth objects observed by ZTF. The system is deployed and adapted for usage within the ZTF Solar-System framework and has significantly reduced human involvement in the streak identification process, from several hours to typically under 10 minutes per day.

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.

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.

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.

Efficient automated detection of flux-transient, reoccurring flux-variable, and moving objects is increasingly important for large-scale astronomical surveys. We present braai, a convolutional-neural-network, deep-learning real/bogus classifier designed to separate genuine astrophysical events and objects from false positive, or bogus, detections in the data of the Zwicky Transient Facility (ZTF), a new robotic time-domain survey currently in operation at the Palomar Observatory in California, USA. Braai demonstrates a state-of-the-art performance as quantified by its low false negative and false positive rates. We describe the open-source software tools used internally at Caltech to archive and access ZTF's alerts and light curves (Kowalski), and to label the data (Zwickyverse). We also report the initial results of the classifier deployment on the Edge Tensor Processing Units (TPUs) that show comparable performance in terms of accuracy, but in a much more (cost-) efficient manner, which has significant implications for current and future surveys.

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

We report on the discovery of a hot Saturn-sized planet (9.916 +/- 0.985 R_earth) around a late F star, EPIC 246865365, observed in Campaign 13 of the K2 mission. We began studying this planet candidate because prior to the release of Gaia DR2, the host star was thought to have been a member (> 90% membership probability) of the approximately 1 Gyr open cluster NGC 1817 based on its kinematics and photometric distance. We identify the host star (among three stars within the K2 photometric aperture) using seeing-limited photometry and rule out false positive scenarios using adaptive optics imaging and radial velocity observations. We statistically validate EPIC 246865365b by calculating a false positive probability rate of 0.01%. However, we also show using new kinematic measurements provided by Gaia DR2 and our measured radial velocity of the system that EPIC 246865365 is unassociated with the cluster NGC 1817. Therefore, the long-running search for a giant transiting planet in an open cluster remains fruitless. Finally, we note that our use of seeing-limited photometry is a good demonstration of similar techniques that are already being used to follow up TESS planet candidates, especially in crowded regions.

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.

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 describe ZStreak, a semi-real-time pipeline specialized in detecting small, fast-moving near-Earth asteroids (NEAs) that is currently operating on the data from the newly-commissioned Zwicky Transient Facility (ZTF) survey. Based on a prototype originally developed by Waszczak et al. (2017) for the Palomar Transient Factory (PTF), the predecessor of ZTF, ZStreak features an improved machine-learning model that can cope with the $10\times$ data rate increment between PTF and ZTF. Since its first discovery on 2018 February 5 (2018 CL), ZTF/ZStreak has discovered $45$ confirmed new NEAs over a total of 232 observable nights until 2018 December 31. Most of the discoveries are small NEAs, with diameters less than $\sim100$ m. By analyzing the discovery circumstances, we find that objects having the first to last detection time interval under 2 hr are at risk of being lost. We will further improve real-time follow-up capabilities, and work on suppressing false positives using deep learning.

Main-belt asteroid (6478) Gault unexpectedly sprouted two tails in late 2018 and early 2019, identifying it as a new active asteroid. Here we present observations obtained by the 1.2-m Zwicky Transient Facility survey telescope that provide detailed time-series coverage of the onset and evolution of Gault's activity. Gault exhibited two brightening events, with the first one starting on 2018 Oct. 18$\pm5$ days and a second one starting on 2018 Dec. 24$\pm1$ days. The amounts of mass released are $2\times10^7$ kg and $1\times10^6$ kg, respectively. Based on photometric measurements, each event persisted for about a month. Gault's color has not changed appreciably over time, with a pre-outburst color of $g_\mathrm{PS1}-r_\mathrm{PS1}=0.50\pm0.04$ and $g_\mathrm{PS1}-r_\mathrm{PS1}=0.46\pm0.04$ during the two outbursts. Simulations of dust dynamics shows that the ejecta consists of dust grains of up to 10 $\mu$m in size that are ejected at low velocities below $1~\mathrm{m~s^{-1}}$ regardless of particle sizes. This is consistent with non-sublimation-driven ejection events. The size distribution of the dust exhibits a broken power-law, with particles at 10--20 $\mu$m following a power-law of $-2.5$ to $-3.0$, while larger particles follow a steeper slope of $-4.0$. The derived properties can be explained by either rotational excitation of the nucleus or a merger of a near-contact binary, with the latter scenario to be statistically more likely.

Context. Radio occultation is a technique used to study planetary atmospheres by means of the refraction and absorption of a spacecraft carrier signal through the atmosphere of the celestial body of interest, as detected from a ground station on Earth. This technique is usually employed by the deep space tracking and communication facilities (e.g., NASA's Deep Space Network (DSN), ESA's Estrack). Aims. We want to characterize the capabilities of the Planetary Radio Interferometry and Doppler Experiment (PRIDE) technique for radio occultation experiments, using radio telescopes equipped with Very Long Baseline Interferometry (VLBI) instrumentation. Methods. We conducted a test with ESA's Venus Express (VEX), to evaluate the performance of the PRIDE technique for this particular application. We explain in detail the data processing pipeline of radio occultation experiments with PRIDE, based on the collection of so-called open-loop Doppler data with VLBI stations, and perform an error propagation analysis of the technique. Results. With the VEX test case and the corresponding error analysis, we have demonstrated that the PRIDE setup and processing pipeline is suited for radio occultation experiments of planetary bodies. The noise budget of the open-loop Doppler data collected with PRIDE indicated that the uncertainties in the derived density and temperature profiles remain within the range of uncertainties reported in previous Venus' studies. Open-loop Doppler data can probe deeper layers of thick atmospheres, such as that of Venus, when compared to closed-loop Doppler data. Furthermore, PRIDE through the VLBI networks around the world, provides a wide coverage and range of large antenna dishes, that can be used for this type of experiments.

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.

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.

We describe a dynamic science portal called the GROWTH Marshal that allows time-domain astronomers to define science programs, program filters to save sources from different discovery streams, co-ordinate follow-up with various robotic or classical telescopes, analyze the panchromatic follow-up data and generate summary tables for publication. The GROWTH marshal currently serves 137 scientists, 38 science programs and 67 telescopes. Every night, in real-time, several science programs apply various customized filters to the 10^5 nightly alerts from the Zwicky Transient Facility. Here, we describe the schematic and explain the functionality of the various components of this international collaborative platform.

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.

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.

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.

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 report on initial results from the largest infrared AO direct imaging survey searching for wide orbit (>100 AU) massive exoplanets and brown dwarfs as companions around young nearby stars using Robo-AO at the 2.1-m telescope on Kitt Peak, Arizona. The occurrence rates of these rare substellar companions are critical to furthering our understanding of the origin of planetary-mass companions on wide orbits. The observing efficiency of Robo-AO allows us to conduct a survey an order of magnitude larger than previously possible. We commissioned a low-noise high-speed SAPHIRA near-infrared camera to conduct this survey and report on its sensitivity, performance, and data reduction process.

The Kepler light curves used to detect thousands of planetary candidates are susceptible to dilution due to blending with previously unknown nearby stars. With the automated laser adaptive optics instrument, Robo-AO, we have observed 620 nearby stars around 3857 planetary candidates host stars. Many of the nearby stars, however, are not bound to the KOI. In this paper, we quantify the association probability between each KOI and detected nearby stars through several methods. Galactic stellar models and the observed stellar density are used to estimate the number and properties of unbound stars. We estimate the spectral type and distance to 145 KOIs with nearby stars using multi-band observations from Robo-AO and Keck-AO. We find most nearby stars within 1" of a Kepler planetary candidate are likely bound, in agreement with past studies. We use likely bound stars as well as the precise stellar parameters from the California Kepler Survey to search for correlations between stellar binarity and planetary properties. No significant difference between the binarity fraction of single and multiple planet systems is found, and planet hosting stars follow similar binarity trends as field stars, many of which likely host their own non-aligned planets. We find that hot Jupiters are ~4x more likely than other planets to reside in a binary star system. We correct the radius estimates of the planet candidates in characterized systems and find that for likely bound systems, the estimated planetary candidate radii will increase on average by a factor of 1.77, if either star is equally likely to host the planet. We find that the planetary radius gap is robust to the impact of dilution, and find an intriguing 95%-confidence discrepancy between the radius distribution of small planets in single and binary systems.

We present the overall statistical results from the Robo-AO Kepler planetary candidate survey, comprising of 3857 high-angular resolution observations of planetary candidate systems with Robo-AO, an automated laser adaptive optics system. These observations reveal previously unknown nearby stars blended with the planetary candidate host star which alter the derived planetary radii or may be the source of an astrophysical false positive transit signal. In the first three papers in the survey, we detected 440 nearby stars around 3313 planetary candidate host stars. In this paper, we present observations of 532 planetary candidate host stars, detecting 94 companions around 88 stars; 84 of these companions have not previously been observed in high-resolution. We also report 50 more-widely-separated companions near 715 targets previously observed by Robo-AO. We derive corrected planetary radius estimates for the 814 planetary candidates in systems with a detected nearby star. If planetary candidates are equally likely to orbit the primary or secondary star, the radius estimates for planetary candidates in systems with likely bound nearby stars increase by a factor of 1.54, on average. We find that 35 previously-believed rocky planet candidates are likely not rocky due to the presence of nearby stars. From the combined data sets from the complete Robo-AO KOI survey, we find that 14.5\pm0.5% of planetary candidate hosts have a nearby star with 4", while 1.2% have two nearby stars and 0.08% have three. We find that 16% of Earth-sized, 13% of Neptune-sized, 14% of Saturn-sized, and 19% of Jupiter-sized planet candidates have detected nearby stars.

We present an approach to testing the gravitational redshift effect using the RadioAstron satellite. The experiment is based on a modification of the Gravity Probe A scheme of nonrelativistic Doppler compensation and benefits from the highly eccentric orbit and ultra-stable atomic hydrogen maser frequency standard of the RadioAstron satellite. Using the presented techniques we expect to reach an accuracy of the gravitational redshift test of order $10^{-5}$, a magnitude better than that of Gravity Probe A. Data processing is ongoing, our preliminary results agree with the validity of the Einstein Equivalence Principle.

Context. Closed-loop Doppler data obtained by deep space tracking networks (e.g., NASA's DSN and ESA's Estrack) are routinely used for navigation and science applications. By "shadow tracking" the spacecraft signal, Earth-based radio telescopes involved in Planetary Radio Interferometry and Doppler Experiment (PRIDE) can provide open-loop Doppler tracking data when the dedicated deep space tracking facilities are operating in closed-loop mode only. Aims. We explain in detail the data processing pipeline, discuss the capabilities of the technique and its potential applications in planetary science. Methods. We provide the formulation of the observed and computed values of the Doppler data in PRIDE tracking of spacecraft, and demonstrate the quality of the results using as a test case an experiment with ESA's Mars Express spacecraft. Results. We find that the Doppler residuals and the corresponding noise budget of the open-loop Doppler detections obtained with the PRIDE stations are comparable to the closed-loop Doppler detections obtained with the dedicated deep space tracking facilities.

We present the discovery of EPIC 228735255b, a P= 6.57 days Jupiter-mass (M$_P$=1.019$\pm$0.070 M$_{Jup}$) planet transiting a V=12.5 (G5-spectral type) star in an eccentric orbit (e=$0.120^{+0.056}_{-0.046}$) detected using a combination of K2 photometry and ground-based observations. With a radius of 1.095$\pm$0.018R$_{Jup}$ the planet has a bulk density of 0.726$\pm$0.062$\rho_{Jup}$. The host star has a [Fe/H] of 0.12$\pm$0.045, and from the K2 light curve we find a rotation period for the star of 16.3$\pm$0.1 days. This discovery is the 9th hot Jupiter from K2 and highlights K2's ability to detect transiting giant planets at periods slightly longer than traditional, ground-based surveys. This planet is slightly inflated, but much less than others with similar incident fluxes. These are of interest for investigating the inflation mechanism of hot Jupiters.

Robo-AO is an autonomous laser guide star adaptive optics system recently commissioned at the Kitt Peak 2.1-m telescope. Now operating every clear night, Robo-AO at the 2.1-m telescope is the first dedicated adaptive optics observatory. This paper presents the imaging performance of the adaptive optics system in its first eighteen months of operations. For a median seeing value of $1.31^{\prime\prime}$, the average Strehl ratio is 4\% in the $i^\prime$ band and 29\% in the J band. After post-processing, the contrast ratio under sub-arcsecond seeing for a $2\leq i^{\prime} \leq 16$ primary star is five and seven magnitudes at radial offsets of $0.5^{\prime\prime}$ and $1.0^{\prime\prime}$, respectively. The data processing and archiving pipelines run automatically at the end of each night. The first stage of the processing pipeline shifts and adds the data using techniques alternately optimized for stars with high and low SNRs. The second "high contrast" stage of the pipeline is eponymously well suited to finding faint stellar companions.

Two transiting planets have been identified orbiting K2 target EPIC 220674823. One object is an ultra-short-period planet (USP) with a period of just 0.57 days (13.7 hours), while the other has a period of 13.3 days. Both planets are small, with the former having a radius of R_p1=1.5 R_E and the latter R_p2=2.5 R_E. Follow-up observations, including radial velocity (with uncertainties of 110 m/s) and high-resolution adaptive optics imagery, show no signs of stellar companions. EPIC 220674823 is the 12th confirmed or validated planetary system in which an ultra-short-period planet (i.e., having an orbital period less than 1 day) is accompanied by at least one additional planet, suggesting that such systems may be common and must be accounted for in models for the formation and evolution of such extreme systems.

We report the discovery of two super-Earth-sized planets transiting the bright (V = 8.94, K = 7.07) nearby late G-dwarf HD 3167, using data collected by the K2 mission. The inner planet, HD 3167 b, has a radius of 1.6 R_e and an ultra-short orbital period of only 0.96 days. The outer planet, HD 3167 c, has a radius of 2.9 R_e and orbits its host star every 29.85 days. At a distance of just 45.8 +/- 2.2 pc, HD 3167 is one of the closest and brightest stars hosting multiple transiting planets, making HD 3167 b and c well suited for follow-up observations. The star is chromospherically inactive with low rotational line-broadening, ideal for radial velocity observations to measure the planets' masses. The outer planet is large enough that it likely has a thick gaseous envelope which could be studied via transmission spectroscopy. Planets transiting bright, nearby stars like HD 3167 are valuable objects to study leading up to the launch of the James Webb Space Telescope.

We have started an initial three-year deployment of Robo-AO at the 2.1-m telescope at Kitt Peak, Arizona as of November 2015. We report here on the project status and two new developments with the Robo-AO KP system: the commissioning of a sub-electron readnoise SAPHIRA near-infrared camera, which will allow us to widen the scope of possible targets to low-mass stellar and substellar objects; and, performance analysis and tuning of the adaptive optics system, which will improve the sensitivity to these objects. Commissioning of the near-infrared camera and optimizing the AO performance occur in parallel with ongoing visible-light science programs.

The Kepler mission has revealed a great diversity of planetary systems and architectures, but most of the planets discovered by Kepler orbit faint stars. Using new data from the K2 mission, we present the discovery of a five planet system transiting a bright (V = 8.9, K = 7.7) star called HIP 41378. HIP 41378 is a slightly metal-poor late F-type star with moderate rotation (v sin(i) = 7 km/s) and lies at a distance of 116 +/- 18 from Earth. We find that HIP 41378 hosts two sub-Neptune sized planets orbiting 3.5% outside a 2:1 period commensurability in 15.6 and 31.7 day orbits. In addition, we detect three planets which each transit once during the 75 days spanned by K2 observations. One planet is Neptune sized in a likely ~160 day orbit, one is sub-Saturn sized likely in a ~130 day orbit, and one is a Jupiter sized planet in a likely ~1 year orbit. We show that these estimates for the orbital periods can be made more precise by taking into account dynamical stability considerations. We also calculate the distribution of stellar reflex velocities expected for this system, and show that it provides a good target for future radial velocity observations. If a precise orbital period can be determined for the outer Jovian planet through future observations, it will be an excellent candidate for follow-up transit observations to study its atmosphere and measure its oblateness.

The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler measurements of planetary spacecraft using the Mars Express spacecraft. The technique will be used in the framework of Planetary Radio Interferometry and Doppler Experiments in various planetary missions, in particular in fly-by mode. We advanced a novel approach to spacecraft data processing using the techniques of Doppler and phase-referenced very long baseline interferometry spacecraft tracking. We achieved, on average, mHz precision (30 {\mu}m/s at a 10 seconds integration time) for radial three-way Doppler estimates and sub-nanoradian precision for lateral position measurements, which in a linear measure (at a distance of 1.4 AU) corresponds to ~50 m.

A test of a cornerstone of general relativity, the gravitational redshift effect, is currently being conducted with the RadioAstron spacecraft, which is on a highly eccentric orbit around Earth. Using ground radio telescopes to record the spacecraft signal, synchronized to its ultra-stable on-board H-maser, we can probe the varying flow of time on board with unprecedented accuracy. The observations performed so far, currently being analyzed, have already allowed us to measure the effect with a relative accuracy of $4\times10^{-4}$. We expect to reach $2.5\times10^{-5}$ with additional observations in 2016, an improvement of almost a magnitude over the 40-year old result of the GP-A mission.

In this paper a description is given of the SFXC software correlator, developed and maintained at the Joint Institute for VLBI in Europe (JIVE). The software is designed to run on generic Linux-based computing clusters. The correlation algorithm is explained in detail, as are some of the novel modes that software correlation has enabled, such as wide-field VLBI imaging through the use of multiple phase centres and pulsar gating and binning. This is followed by an overview of the software architecture. Finally, the performance of the correlator as a function of number of CPU cores, telescopes and spectral channels is shown.

The accuracy of orbit determination has a strong impact on the scientific output of the Space VLBI mission RadioAstron. The aim of this work is to improve the RadioAstron orbit reconstruction by means of sophisticated dynamical modelling of its motion in combination with multi-station Doppler tracking of the RadioAstron spacecraft. The improved orbital solution is demonstrated using Doppler measurements of the RadioAstron downlink signal and by correlating VLBI observations made by RadioAstron with ground-based telescopes using the enhanced orbit determination data. Orbit determination accuracy has been significantly improved from ~ 600 m in 3D position and ~ 2 cm/s in 3D velocity to several tens of metres and mm/s, respectively.

Aims. The phase scintillation of the European Space Agency's (ESA) Venus Express (VEX) spacecraft telemetry signal was observed at X-band (\lambda = 3.6 cm) with a number of radio telescopes of the European VLBI Network (EVN) in the period 2009-2013. Methods. We found a phase fluctuation spectrum along the Venus orbit with a nearly constant spectral index of -2.42 +/-0.25 over the full range of solar elongation angles from 0{\deg} to 45{\deg}, which is consistent with Kolmogorov turbulence. Radio astronomical observations of spacecraft signals within the solar system give a unique opportunity to study the temporal behaviour of the signal's phase fluctuations caused by its propagation through the interplanetary plasma and the Earth's ionosphere. This gives complementary data to the classical interplanetary scintillation (IPS) study based on observations of the flux variability of distant natural radio sources. Results. We present here our technique and the results on IPS. We compare these with the total electron content (TEC) for the line of sight through the solar wind. Finally, we evaluate the applicability of the presented technique to phase-referencing Very Long Baseline Interferometry (VLBI) and Doppler observations of currently operational and prospective space missions.

We present the results of several multi-station Very Long Baseline Interferometry (VLBI) experiments conducted with the ESA spacecraft Venus Express as a target. To determine the true capabilities of VLBI tracking for future planetary missions in the solar system, it is necessary to demonstrate the accuracy of the method for existing operational spacecraft. We describe the software pipeline for the processing of phase referencing near-field VLBI observations and present results of the ESA Venus Express spacecraft observing campaign conducted in 2010-2011. We show that a highly accurate determination of spacecraft state-vectors is achievable with our method. The consistency of the positions indicates that an internal rms accuracy of 0.1 mas has been achieved. However, systematic effects produce offsets up to 1 mas, but can be reduced by better modelling of the troposphere and ionosphere and closer target-calibrator configurations.