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)

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

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

Stars aligned in thin stream-like features (feathers), with widths of $\delta\sim1-10^{\circ}$ and lengths as large as $\Delta l\sim180^{\circ}$, have been observed towards the Anticenter of our Galaxy and their properties mapped in abundances and phase-space. We study their origin by analysing similar features arising in an N-body simulation of a Galactic disc interacting with a Sagittarius-like dwarf spheroidal galaxy (Sgr). By following the orbits of the particles identified as contributing to feathers backwards in time, we trace their excitation to one of Sgr's previous pericentric passages. These particles initially span a large range of phase-angles but a tight range of radii, suggesting they provide a probe of populations in distinct annuli in the outer Galactic disc. The structures are long lived and persist after multiple passages on timescales of $\sim4 \,\rm{Gyrs}$. On the sky, they exhibit oscillatory motion that can be traced with a single orbit mapped over much of their full length and with amplitudes and gradients similar to those observed. We demonstrate how these properties of feathers may be exploited to measure the potential, its flattening, as well as infer the strength of recent potential perturbations.

The Monoceros Ring (also known as the Galactic Anticenter Stellar Structure) and A13 are stellar overdensities at estimated heliocentric distances of $d \sim 11$ kpc and 15 kpc observed at low Galactic latitudes towards the anticenter of our Galaxy. While these overdensities were initially thought to be remnants of a tidally-disrupted satellite galaxy, an alternate scenario is that they are composed of stars from the Milky Way (MW) disk kicked out to their current location due to interactions between a satellite galaxy and the disk. To test this scenario, we study the stellar populations of the Monoceros Ring and A13 by measuring the number of RR Lyrae and M giant stars associated with these overdensities. We obtain low-resolution spectroscopy for RR Lyrae stars in the two structures and measure radial velocities to compare with previously measured velocities for M giant stars in the regions of the Monoceros Ring and A13, to assess the fraction of RR Lyrae to M giant stars ($f_{RR:MG}$) in A13 and Mon/GASS. We perform velocity modeling on 153 RR Lyrae stars (116 in the Monoceros Ring and 37 in A13) and find that both structures have very low $f_{RR:MG}$. The results support a scenario in which stars in A13 and Mon/GASS formed in the MW disk. We discuss a possible association between Mon/GASS, A13, and the Triangulum-Andromeda overdensity based on their similar velocity distributions and $f_{RR:MG}$.