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.

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

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 report the discovery and panchromatic followup observations of the young Type Ic supernova, SN2020oi, in M100, a grand design spiral galaxy at a mere distance of $14$ Mpc. We followed up with observations at radio, X-ray and optical wavelengths from only a few days to several months after explosion. The optical behaviour of the supernova is similar to those of other normal Type Ic supernovae. The event was not detected in the X-ray band but our radio observation revealed a bright mJy source ($L_{\nu} \approx 1.2 \times 10^{27} {\rm erg\,s}^{-1} {\rm Hz}^{-1}$). Given, the relatively small number of stripped envelope SNe for which radio emission is detectable, we used this opportunity to perform a detailed analysis of the comprehensive radio dataset we obtained. The radio emitting electrons initially experience a phase of inverse Compton cooling which leads to steepening of the spectral index of the radio emission. Our analysis of the cooling frequency points to a large deviation from equipartition at the level of $\epsilon_e/\epsilon_B \gtrsim 200$, similar to a few other cases of stripped envelope SNe. Our modeling of the radio data suggests that the shockwave driven by the SN ejecta into the circumstellar matter (CSM) is moving at $\sim 3\times 10^{4}\,{\rm km\,s}^{-1}$. Assuming a constant mass-loss from the stellar progenitor, we find that the mass-loss rate is $\dot{M} \approx 1.4\times 10^{-4}\,{M}_{\odot}\,{\rm yr}^{-1}$, for an assumed wind velocity of $1000\,{\rm km\,s}^{-1}$. The temporal evolution of the radio emission suggests a radial CSM density structure steeper than the standard $r^{-2}$.

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

The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital motion in a binary. However, this conclusion was based previously on sparse UV data which were not taken simultaneously with the optical data. Here we report simultaneous follow-up observations of PG1302-102 with the Ultraviolet Optical Telescope on the Neil Gehrels Swift Observatory in six optical + UV bands. The additional nine Swift observations produce light curves roughly consistent with the trend under the Doppler boost hypothesis, which predicts that UV variability should track the optical, but with a ~2.2 times higher amplitude. We perform a statistical analysis to quantitatively test this hypothesis. We find that the data are consistent with the Doppler boost hypothesis when we compare the the amplitudes in optical B-band and UV light curves. However, the ratio of UV to V-band variability is larger than expected and is consistent with the Doppler model, only if either the UV/optical spectral slopes vary, the stochastic variability makes a large contribution in the UV, or the sparse new optical data underestimate the true optical variability. We have evidence for the latter from comparison with the optical light curve from ASAS-SN. Additionally, the simultaneous analysis of all four bands strongly disfavors the Doppler boost model whenever Swift V-band is involved. Additional, simultaneous optical + UV observations tracing out another cycle of the 5.2-year proposed periodicity should lead to a definitive conclusion.

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 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 ASASSN-19dj, a nearby tidal disruption event (TDE) discovered in the post-starburst galaxy KUG 0810+227 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of d $\simeq$ 98 Mpc. We observed ASASSN-19dj from $-$21 to 392 days relative to peak UV/optical emission using high-cadence, multi-wavelength spectroscopy and photometry. From the ASAS-SN $g$-band data, we determine that the TDE began to brighten on 2019 February 6.8 and for the first 25 days the rise was consistent with a flux $\propto$ $t^2$ power-law. ASASSN-19dj peaked in the UV/optical on 2019 March 6.5 (MJD = 58548.5) at a bolometric luminosity of $L = (6.2 \pm 0.2) \times 10^{44} \text{ erg s}^{-1}$. Initially remaining roughly constant in X-rays and slowly fading in the UV/optical, the X-ray flux increased by over an order of magnitude $\sim$225 days after peak, resulting from the expansion of the X-ray emitting surface. The late-time X-ray emission is well-fit by a blackbody with an effective radius of $\sim 1 \times 10^{12} \text{ cm}$ and a temperature of $\sim 6 \times 10^{5} \text{ K}$. Analysis of Catalina Real-Time Transient Survey images reveals a nuclear outburst roughly 14.5 years earlier with a smooth decline and a luminosity of $L_V$ $\geq$ $1.4 \times 10^{43}$ erg s$^{-1}$, although the nature of the flare is unknown. ASASSN-19dj occurred in the most extreme post-starburst galaxy yet to host a TDE, with Lick H$\delta_{A}$ = $7.67 \pm 0.17$ \AA.

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.

The discovery rate of fast radio bursts (FRBs) is increasing dramatically thanks to new radio facilities. Meanwhile, wide-field instruments such as the 47 deg$^2$ Zwicky Transient Facility (ZTF) survey the optical sky to study transient and variable sources. We present serendipitous ZTF observations of the CHIME repeating source FRB 180916.J0158+65, that was localized to a spiral galaxy 149 Mpc away and is the first FRB suggesting periodic modulation in its activity. While 147 ZTF exposures corresponded to expected high-activity periods of this FRB, no single ZTF exposure was at the same time as a CHIME detection. No $>3\sigma$ optical source was found at the FRB location in 683 ZTF exposures, totalling 5.69 hours of integration time. We combined ZTF upper limits and expected repetitions from FRB 180916.J0158+65 in a statistical framework using a Weibull distribution, agnostic of periodic modulation priors. The analysis yielded a constraint on the ratio between the optical and radio fluences of $\eta \lesssim 200$, corresponding to an optical energy $E_{\rm opt} \lesssim 3 \times 10^{46}$ erg for a fiducial 10 Jy ms FRB (90% confidence). A deeper (but less statistically robust) constraint of $\eta \lesssim 3$ can be placed assuming a rate of $r(>5$ Jy ms)= hr$^{-1}$ and $1.2\pm 1.1$ FRB occurring during exposures taken in high-activity windows. The constraint can be improved with shorter per-image exposures and longer integration time, or observing FRBs at higher Galactic latitudes. This work demonstrated how current surveys can statistically constrain multi-wavelength counterparts to FRBs even without deliberately scheduled simultaneous radio observation.

We study the link between the kinematic-morphology of galaxies, as inferred from integral-field stellar kinematics, and their relation between mass and star formation rate (SFR). Our sample consists of $\sim 3200$ galaxies with integral-field spectroscopic data from the MaNGA survey with available determinations of their effective stellar angular momentum within the half-light radius $\lambda_{R_e}$. We find that for star-forming galaxies, namely along the star formation main sequence (SFMS), the $\lambda_{R_e}$ values remain large and almost unchanged over about two orders of magnitude in stellar mass, with the exception of the lowest masses $\mathcal{M}_{\star}\lesssim2\times10^{9} \mathcal{M}_{\odot}$, where $\lambda_{R_e}$ slightly decreases. The SFMS is dominated by spiral galaxies with small bulges. Below the SFMS, but above the characteristic stellar mass $\mathcal{M}_{\rm crit}\approx2\times10^{11} \mathcal{M}_{\odot}$, there is a sharp decrease in $\lambda_{R_e}$ with decreasing star formation rate: massive galaxies well below the SFMS are mainly slow-rotator early-type galaxies, namely genuinely spheroidal galaxies without disks. Below the SFMS and below $\mathcal{M}_{\rm crit}$ the decrease of $\lambda_{R_e}$ with decreasing SFR becomes modest or nearly absent: low-mass galaxies well below the SFMS, are fast-rotator early-type galaxies, and contain fast-rotating stellar disks like their star-forming counterparts. We also find a small but clear environmental dependence for the massive galaxies: in the mass range $10^{10.9}-10^{11.5} \mathcal{M}_{\odot}$, galaxies in rich groups or denser regions or classified as central galaxies have lower values of $\lambda_{R_e}$. While no environmental dependence is found for galaxies of lower mass. We discuss how our results can be understood as due to the different star formation and mass assembly histories of galaxies with varying mass.

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.

Recently, Kundu et al (2019) reported that the globular cluster NGC 5024 (M53) possesses five extra-tidal RR Lyrae. In fact, four of them were instead known members of a nearby globular cluster NGC 5053. The status of the remaining extra-tidal RR Lyrae is controversial depending on the adopted tidal radius of NGC 5024. We have also searched for additional extra-tidal RR Lyrae within an area of $\sim8$~deg$^2$ covering both globular clusters. This includes other known RR Lyrae within the search area, as well as stars that fall within the expected range of magnitudes and colors for RR Lyrae (and yet outside the cutoff of 2/3 of the tidal radii of each globular clusters for something to be called "extra-tidal") if they were extra-tidal RR Lyrae candidates for NGC 5024 or NGC 5053. Based on the the proper motion information and their locations on the color-magnitude diagram, none of the known RR Lyrae belong to the extra-tidal RR Lyrae of either globular clusters. In the cases where the stars satisfied the magnitude and color ranges of RR Lyrae, analysis of time series data taken from the Zwicky Transient Facility do not reveal periodicities, suggesting that none of these stars are RR Lyrae. We conclude that there are no extra-tidal RR Lyrae associated with either NGC 5024 or NGC 5053 located within our search area.

Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently associated one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the association of a radio-emitting tidal disruption event (AT2019dsg) with another high-energy neutrino, identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility (ZTF). The probability of finding any radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. The association suggests that tidal disruption events contribute to the cosmic neutrino flux. Unlike previous work which considered the rare subset of tidal disruption events with relativistic jets, our observations of AT2019dsg suggest an empirical model with a mildly-relativistic outflow.

Colors of Type Ia supernovae in the first few days after explosion provide a potential discriminant between different models. In this paper, we present $g-r$ colors of 65 Type Ia supernovae discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that $g-r$ colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties ($\sigma_\mathrm{noise}\sim\sigma_\mathrm{int}\sim$ 0.18 mag). Colors are nearly constant starting from 6 days after first light ($g-r\sim-0.15$ mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of $\sim-0.25$ mag day$^{-1}$) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of $^{56}$Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers ($M_\mathrm{He}=0.01\,M_\odot$) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by "double-detonation" models with larger helium masses. We finally identify a significant correlation between the early-time $g-r$ slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value=0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on $B-V$ colors.

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.

Accurate photometric redshift (photo-$z$) estimates are essential to the cosmological science goals of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). In this work we use simulated photometry for mock galaxy catalogs to explore how LSST photo-$z$ estimates can be improved by the addition of near-infrared (NIR) and/or ultraviolet (UV) photometry from the Euclid, WFIRST, and/or CASTOR space telescopes. Generally, we find that deeper optical photometry can reduce the standard deviation of the photo-$z$ estimates more than adding NIR or UV filters, but that additional filters are the only way to significantly lower the fraction of galaxies with catastrophically under- or over-estimated photo-$z$. For Euclid, we find that the addition of ${JH}$ $5{\sigma}$ photometric detections can reduce the standard deviation for galaxies with $z>1$ ($z>0.3$) by ${\sim}20\%$ (${\sim}10\%$), and the fraction of outliers by ${\sim}40\%$ (${\sim}25\%$). For WFIRST, we show how the addition of deep ${YJHK}$ photometry could reduce the standard deviation by ${\gtrsim}50\%$ at $z>1.5$ and drastically reduce the fraction of outliers to just ${\sim}2\%$ overall. For CASTOR, we find that the addition of its ${UV}$ and $u$-band photometry could reduce the standard deviation by ${\sim}30\%$ and the fraction of outliers by ${\sim}50\%$ for galaxies with $z<0.5$. We also evaluate the photo-$z$ results within sky areas that overlap with both the NIR and UV surveys, and when spectroscopic training sets built from the surveys' small-area deep fields are used.

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.

Quasars have long been known as intrinsically variable sources, but the physical mechanism underlying the temporal optical/UV variability is still not well understood. We propose a novel nonparametric method for modeling and forecasting the optical variability of quasars utilizing an autoencoder neural network to gain insight into the underlying processes. The autoencoder is trained with ~15,000 decade-long quasar light curves obtained by the Catalina Real-time Transient Survey selected with negligible flux contamination from the host galaxy. The autoencoder's performance in forecasting the temporal flux variation of quasars is superior to that of the damped random walk process. We find a temporal asymmetry in the optical variability and a novel relation - the amplitude of the variability asymmetry decreases as luminosity and/or black hole mass increases - is suggested with the help of autoencoded features. The characteristics of the variability asymmetry are in agreement with those from the self-organized disk instability model, which predicts that the magnitude of the variability asymmetry decreases as the ratio of the diffusion mass to inflow mass in the accretion disk increases.

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.

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 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 present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$\alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $\Gamma$, and the X-ray-to-optical power slopes, $\alpha_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

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

We present a 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.

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced for the Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. We report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

We present 637 low-redshift optical spectra collected by the Berkeley Supernova Ia Program (BSNIP) between 2009 and 2018, almost entirely with the Kast double spectrograph on the Shane 3~m telescope at Lick Observatory. We describe our automated spectral classification scheme and arrive at a final set of 626 spectra (of 242 objects) that are unambiguously classified as belonging to Type Ia supernovae (SNe~Ia). Of these, 70 spectra of 30 objects are classified as spectroscopically peculiar (i.e., not matching the spectral signatures of "normal" SNe~Ia) and 79 SNe~Ia (covered by 328 spectra) have complementary photometric coverage. The median SN in our final set has one epoch of spectroscopy, has a redshift of 0.0208 (with a low of 0.0007 and high of 0.1921), and is first observed spectroscopically 1.1 days after maximum light. The constituent spectra are of high quality, with a median signal-to-noise ratio of 31.8 pixel$^{-1}$, and have broad wavelength coverage, with $\sim 95\%$ covering at least 3700--9800~\AA. We analyze our dataset, focusing on quantitative measurements (e.g., velocities, pseudo-equivalent widths) of the evolution of prominent spectral features in the available early-time and late-time spectra. The data are available to the community, and we encourage future studies to incorporate our spectra in their analyses.

While it is clear that Type Ia supernovae (SNe) are the result of thermonuclear explosions in C/O white dwarfs (WDs), a great deal remains uncertain about the binary companion that facilitates the explosive disruption of the WD. Here, we present a comprehensive analysis of a unique, and large, data set of 127 SNe Ia with exquisite coverage by the Zwicky Transient Facility (ZTF). High-cadence (6 observations per night) ZTF observations allow us to measure the SN rise time and examine its initial evolution. We develop a Bayesian framework to model the early rise as a power-law in time, which enables the inclusion of priors in our model. For a volume-limited subset of normal SNe Ia, we find the mean power-law index is consistent with 2 in the $r_\mathrm{ztf}$-band ($\alpha_r = 2.01\pm0.02$), as expected in the expanding fireball model. There are, however, individual SNe that are clearly inconsistent with $\alpha_r=2$. We estimate a mean rise time of 18.5$\,$d (with a range extending from $\sim$15$-$22$\,$d), though this is subject to the adopted prior. We identify an important, previously unknown, bias whereby the rise times for higher redshift SNe within a flux-limited survey are systematically underestimated. This effect can be partially alleviated if the power-law index is fixed to $\alpha=2$, in which case we estimate a mean rise time of 21.0$\,$d (with a range from $\sim$18$-$23$\,$d). The sample includes a handful or rare and peculiar SNe Ia. Finally, we conclude with a discussion of lessons learned from the ZTF sample that can eventually be applied to Large Synoptic Survey Telescope observations.

Comet 2I/Borisov, the first unambiguous interstellar comet ever found, was discovered in August 2019 at $\sim3$ au from the Sun on its inbound leg. No pre-discovery detection beyond 3 au has yet been reported, mostly due to the comet's proximity to the Sun as seen from the Earth. Here we present a search for pre-discovery detections of comet Borisov using images taken by the Catalina Sky Survey (CSS), Pan-STARRS and Zwicky Transient Facility (ZTF), with a further comprehensive follow-up campaign being presented in \citet{Bolin2019}. We identified comet Borisov in ZTF images taken in May 2019 and use these data to update its orbit. This allowed us to identify the comet in images acquired as far back as December 2018, when it was 7.8 au from the Sun. The comet was not detected in November 2018 when it was 8.6 au from the Sun, possibly implying an onset of activity around this time. This suggests that the activity of the comet is either driven by a more volatile species other than H$_2$O, such as CO or CO$_2$, or by exothermic crystallization of amorphous ice. We derive the radius of the nucleus to be $<7$ km using the non-detection in November 2018, and estimate an area of $\sim0.5$---$10 \mathrm{km^2}$ has been active between December 2018 and September 2019, though this number is model-dependent and is highly uncertain. The behavior of comet Borisov during its inbound leg is observationally consistent with dynamically new comets observed in our solar system, suggesting some similarities between the two.

The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

Near-Earth Objects (NEOs) that orbit the Sun on or within Earth's orbit are tricky to detect for Earth-based observers due to their proximity to the Sun in the sky. These small bodies hold clues to the dynamical history of the inner solar system as well as the physical evolution of planetesimals in extreme environments. Populations in this region include the Atira and Vatira asteroids, as well as Venus and Earth co-orbital asteroids. Here we present a twilight search for these small bodies, conducted using the 1.2-m Oschin Schmidt and the Zwicky Transient Facility (ZTF) camera at Palomar Observatory. The ZTF twilight survey operates at solar elongations down to $35^\circ$ with limiting magnitude of $r=19.5$. During a total of 40 evening sessions and 62 morning sessions conducted between 2018 November 15 and 2019 June 23, we detected 6 Atiras, including 2 new discoveries 2019 AQ$_3$ and 2019 LF$_6$, but no Vatiras or Earth/Venus co-orbital asteroids. NEO population models show that these new discoveries are likely only the tip of the iceberg, with the bulk of the population yet to be found. The population models also suggest that we have only detected 5--$7\%$ of the $H<20$ Atira population over the 7-month survey. Co-orbital asteroids are smaller in diameters and require deeper surveys. A systematic and efficient survey of the near-Sun region will require deeper searches and/or facilities that can operate at small solar elongations.

We report on three redshift $z>2$ quasars with dramatic changes in their C IV emission lines, the first sample of changing-look quasars (CLQs) at high redshift. This is also the first time the changing-look behaviour has been seen in a high-ionisation emission line. SDSS J1205+3422, J1638+2827, and J2228+2201 show interesting behaviour in their observed optical light curves, and subsequent spectroscopy shows significant changes in the C IV broad emission line, with both line collapse and emergence being displayed on rest-frame timescales of $\sim$240-1640 days. These are rapid changes, especially when considering virial black hole mass estimates of $M_{\rm BH} > 10^{9} M_{\odot}$ for all three quasars. Continuum and emission line measurements from the three quasars show changes in the continuum-equivalent width plane with the CLQs seen to be on the edge of the full population distribution, and showing indications of an intrinsic Baldwin effect. We put these observations in context with recent state-change models, and note that even in their observed low-state, the C IV CLQs are generally above $\sim$5\% in Eddington luminosity.

Boosted dark matter (BDM) is a well-motivated class of dark matter candidates in which a small subcomponent is relativistic at the present time. In this work, we lay the foundation for BDM searches via its hadronic interactions at DUNE. We investigate BDM-nucleus scattering in detail by developing new dedicated event generation techniques with a parameterized detector simulation. We demonstrate that the excellent event reconstruction capability of DUNE enables a significantly broader discovery reach.

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising ($1.4\pm0.1$ mag/hr) and luminous ($M_{g,\mathrm{peak}}=-20$ mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity ($L_{\mathrm{bol}} \gtrsim 3 \times 10^{44}$ erg $\mathrm{sec}^{-1}$), the short rise time ($t_{\mathrm{rise}}= 3$ days in $g$-band), and the blue colors at peak ($g-r\sim-0.4$) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature ($T_{\mathrm{eff}}\gtrsim40,000$ K) spectra of a stripped-envelope SN. A retrospective search revealed luminous ($M_g \sim M_r \approx -14\,$mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release $E_\mathrm{\gamma,iso}<4.9 \times 10^{48}$ erg, a limit on X-ray emission $L_{\mathrm{X}} < 10^{40}\,$erg sec$^{-1}$, and a limit on radio emission $\nu L_\nu \lesssim 10^{37}\,$erg sec$^{-1}$. Taken together, we find that the early ($<10\,$days) data are best explained by shock breakout in a massive shell of dense circumstellar material ($0.02\,M_\odot$) at large radii ($3 \times 10^{14}\,$cm) that was ejected in eruptive pre-explosion mass-loss episodes. The late-time ($>10$ days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56.

Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

Galaxy interactions leave imprints in the motions of their stars, and so observing the two-dimensional stellar kinematics allows us to uncover their formation process. Slow rotators, which have stellar orbits dominated by random motions, are thought to be the fossil relics of a sequence of multiple gas-poor mergers, in an environment where the cold gas required to form new stars is nearly absent. Indeed, observations of a handful of nearby galaxy clusters have indicated that slow rotators are preferentially found in the gas-poor, dense cores of clusters, which themselves must form by merging of smaller groups. However, the generality of this result and connection between kinematics and environment is currently unclear, as recent studies have suggested that, at given stellar mass, the environment does not influence the formation of slow rotators. Here we address this issue by combining a careful quality-assessed sample selection with two-dimensional stellar kinematics from a large galaxy survey and a novel photometric classification approach where kinematics are unavailable. We obtain the first complete census of the location of massive slow rotators in ten large clusters: in all cases, slow rotators are extremely rare and generally trace the clusters density peaks. This result unambiguously establishes that massive slow rotators are the relics of violent hierarchical cluster formation.

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

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

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

The majority of galaxy group catalogues available in the literature use the popular friends-of-friends algorithm which links galaxies using a linking length. One potential drawback to this approach is that clusters of point can be link with thin bridges which may not be desirable. Furthermore, these algorithms are designed with large-scales galaxy surveys in mind rather than small-scale, local galaxy environments, where attention to detail is important. Here we present a new simple group finder algorithm, TD-ENCLOSER, that finds the group that encloses a target galaxy of interest. TD-ENCLOSER is based on the kernel density estimation method which treats each galaxy, represented by a zero-dimensional particle, as a two-dimensional circular Gaussian. The algorithm assigns galaxies to peaks in the density field in order of density in descending order ("Top Down") so that galaxy groups "grow" around the density peaks. Outliers in under-dense regions are prevented from joining groups by a specified hard threshold, while outliers at the group edges are clipped below a soft (blurred) interior density level. The group assignments are largely insensitive to all free parameters apart from the hard density threshold and the kernel standard deviation, although this is a known feature of density-based group finder algorithms, and operates with a computing speed that increases linearly with the size of the input sample. In preparation for a companion paper, we also present a simple algorithm to select unique representative groups when duplicates occur. TD-ENCLOSER produces results comparable to those from a widely used catalogue, as shown in a companion paper. A smoothing scale of 0.3 Mpc provides the most realistic group structure.

It has been observed that low redshift early-type galaxies can be separated into slow and fast rotators according to a proxy of specific stellar angular momentum, $\lambda_{R_e}$. Detailed studies of a handful of nearby clusters have shown that slow rotators are generally found at the centres of clusters where the number density is highest, whereas the fast rotators trace the trend followed by early-type galaxies of increasing in number with local density. In this paper, we study the environmental distribution of slow and fast rotators using the stellar kinematics of about 3900 galaxies from the Sloan Digital Sky Survey's Mapping Nearby Galaxies at Apache Point Observatory survey. For galaxies in groups closer than $z=0.08$ that are not observed with MaNGA but satisfy the necessary conditions for slow rotators, we visually assign slow/fast rotator classifications to obtain a complete sample. Our final catalogue contains about 900 groups of five or more members. We observe the kinematic morphology-density (kT-$\Sigma$) relation for each group and find an increasing fraction of massive slow rotators with increasing number density. We provide evidence suggesting that the observed lack of trends in angular momentum with environment at fixed stellar mass is partly because the maximum density varies between clusters, and that the locations of massive slow rotators are strongly correlated with peak densities in galaxy groups and clusters. We conclude that the (projected) number density relative to the cluster peak is more fundamental than the absolute number density in influencing the abundance of slow rotators. We find that the kT-$\Sigma$ relation does exist at fixed stellar mass, and we rule out the hypothesis that the kT-$\Sigma$ relation is a result of dynamical friction alone, instead arguing that massive slow rotators grow hierarchically in tandem with their host clusters.

Galaxy morphology is inextricably linked to environment. The morphology-density relation quantifies this relationship. However, optical morphology is only loosely related to the kinematic structure of galaxies, and about two thirds of elliptical ("spheroidal") galaxies are actually misclassified face-on disks, and would appear flattened if view edge on. A more robust classification is the slow/fast rotator classification which describes the kinematic structure. Slow and fast rotators form a bimodality in galaxy properties and are thought to follow distinct evolutionary paths, and so a kinematic morphology-density (kT-$\Sigma$) relation is more meaningful. To date the kT-$\Sigma$ relation has only been studied for a handful of nearby clusters,or across large numbers of clusters but with incomplete coverage. Here, we combine stellar kinematics obtained with the Sloan Digital Sky Survey's (SDSS) Mapping Nearby Galaxies at Apache Point Observatory survey with classifications using a novel visual method to obtain the largest complete census of slow and fast rotators in nearby galaxy groups and clusters. To account for incompleteness in the SDSS spectroscopic catalogue, we combine the catalogue with the photometric catalogue, which we clean using empirical criteria. We test our visual classification method and find our false negative rate to be only $\sim$6% (0% for $M \gtrsim 10^{11.7}M_{\odot}$). In contrast, our false positive rate is about 50% implying half of all slow rotator candidates will be confirmed as fast rotators if stellar kinematics become available. Hence, our misclassification is essentially random with only a weak dependence on stellar mass, and as slow rotators are intrinsically rare, the absolute number of misclassifications will be small compared to the sample size. The result is the largest complete census of massive slow rotator candidates and their environments to date.

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.

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

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

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

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

Accretion disks around supermassive black holes (SMBH) are promising sites for stellar mass black hole (BH) mergers due to mass segregation and merger acceleration by disk gas torques. Here we show that a GW-kick at BH merger causes ram-pressure stripping of gas within the BH Hill sphere. If $R_{H}\geq H$, the disk height, an off-center UV flare at $a_{\rm BH} \sim 10^{3}r_{g}$ emerges within $t_{\rm UV} \sim \rm{O}(2 {\rm days})(a_{\rm BH}/10^{3}r_{g})(M_{\rm SMBH}/10^{8}M_{\odot})(v_{\rm kick}/10^{2}\rm{km/s})$ post-merger and lasts O$(R_{H}/v_{\rm kick}) \sim \rm{O}(5 t_{\rm UV}$). The flare emerges with luminosity O($10^{42}{\rm erg/s})(t_{\rm UV}/2{\rm days})^{-1}(M_{\rm Hill}/1M_{\odot})(v_{\rm kick}/10^{2}{\rm km/s})^{2}$. AGN optical/UV photometry alters and asymmetric broad emission line profiles can develop after weeks. If $R_{H}<H$, detectability depends on disk optical depth. Follow-up by large optical sky surveys is optimized for small GW error volumes and for LIGO/Virgo triggers $>50M_{\odot}$.

This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy.

In this work, BV RI light curves of 55 Type II supernovae (SNe II) from the Lick Observatory Supernova Search program obtained with the Katzman Automatic Imaging Telescope and the 1 m Nickel telescope from 2006 to 2018 are presented. Additionally, more than 150 spectra gathered with the 3 m Shane telescope are published. We conduct an analyse of the peak absolute magnitudes, decline rates, and time durations of different phases of the light and colour curves. Typically, our light curves are sampled with a median cadence of 5.5 days for a total of 5093 photometric points. In average V-band plateau declines with a rate of 1.29 mag (100 days)-1, which is consistent with previously published samples. For each band, the plateau slope correlates with the plateau length and the absolute peak magnitude: SNe II with steeper decline have shorter plateau duration and are brighter. A time-evolution analysis of spectral lines in term of velocities and pseudoequivalent widths is also presented in this paper. Our spectroscopic sample ranges between 1 and 200 days post-explosion and has a median ejecta expansion velocity at 50 days post-explosion of 6500 km/s (Halpha line) and a standard dispersion of 2000 km/s. Nebular spectra are in good agreement with theoretical models using a progenitor star having a mass <16 Msol. All the data are available to the community and will help to understand SN II diversity better, and therefore to improve their utility as cosmological distance indicators.

In this white paper, we present the scientific cases for adding narrowband optical filters to the Large Synoptic Survey Telescope (LSST). LSST is currently planning to observe the southern sky in 6 broadband optical filters. Three of the four LSST science themes would benefit from adding narrowband filter observations. We discuss the technical considerations of using narrowband filters with the LSST and lay out the scientific impact that would result on the study of AGB stars, emission line nebula (e.g., supernova remnants and planetary nebulae), photometric redshifts of galaxies, and the determination of stellar parameters.

Time domain and multi-messenger astrophysics are growing and important modes of observational astronomy that will help define astrophysics in the 2020s. Significant effort is being put into developing the components of a follow-up system for dynamically turning survey alerts into data. This system consists of: 1) brokers that will aggregate, classify, and filter alerts; 2) Target Observation Managers (TOMs) for prioritizing targets and managing observations and data; and 3) observatory interfaces, schedulers, and facilities along with data reduction software and science archives. These efforts need continued community support and funding in order to complete and maintain them. Many of the efforts can be community open-source software projects but they will benefit from the leadership of professional software developers. The coordination should be done by institutions that are involved in the follow-up system such as the national observatories (e.g. LSST/Gemini/NOAO Mid-scale/Community Science and Data Center) or a new MMA institute. These tools will help the community to produce the most science from new facilities and will provide new capabilities for all users of the facilities that adopt them.

Time domain and multi-messenger astrophysics are growing and important modes of observational astronomy that will help define astrophysics in the 2020s. Significant effort is being put into developing the components of a follow-up system for dynamically turning survey alerts into data. This system consists of: 1) brokers that will aggregate, classify, and filter alerts; 2) Target Observation Managers (TOMs) for prioritizing targets and managing observations and data; and 3) observatory interfaces, schedulers, and facilities along with data reduction software and science archives. These efforts need continued community support and funding in order to complete and maintain them. Many of the efforts can be community open-source software projects but they will benefit from the leadership of professional software developers. The coordination should be done by institutions that are involved in the follow-up system such as the national observatories (e.g. LSST/Gemini/NOAO Mid-scale/Community Science and Data Center) or a new MMA institute. These tools will help the community to produce the most science from new facilities and will provide new capabilities for all users of the facilities that adopt them.

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

We report the discovery of a mid-infrared variable AGN which is hosted by an ultraluminous infrared galaxy (ULIRG) candidate in the Sloan Stripe 82 field. \textit{WISE} \textit{J030654.88+010833.6} is a red, extended galaxy, which we estimate to be at a photometric redshift of 0.28 $\leq$ z $\leq$ 0.31, based on its optical and near-infrared spectral energy distribution (SED). The factor of two variability over 8 years seen in the \textit{WISE} 3.4 and 4.6 $\mu$m wavelength channels is not clearly correlated with optical variability in archival data. Based on our estimation of the physical parameters of the host galaxy, \textit{J030654.88+010833.6} is possibly a composite AGN/starburst ULIRG in a phase where high star formation $\sim$ 70 M$_{\odot}$ year$^{-1}$ is occurring. Our estimate of the black hole mass to stellar mass ratio also appears to be consistent with that of broad-line AGN in the local universe. The long-term variability of \textit{J030654.88+010833.6} as seen in the \textit{WISE} \textit{W1} and \textit{W2} light curves is likely due to variations in the accretion rate, with the energy being reprocessed by dust in the vicinity of the AGN.

We study the sudden optical and ultraviolet (UV) brightening of 1ES 1927+654, which until now was known as a narrow-line active galactic nucleus (AGN). 1ES 1927+654 was part of the small and peculiar class of "true Type-2" AGN, which lack broad emission lines and line-of-sight obscuration. Our high-cadence spectroscopic monitoring captures the appearance of a blue, featureless continuum, followed several weeks later by the appearance of broad Balmer emission lines. This timescale is generally consistent with the expected light travel time between the central engine and the broad-line emission region in (persistent) broad-line AGN. Hubble Space Telescope spectroscopy reveals no evidence for broad UV emission lines (e.g., CIV1549, CIII]1909, MgII2798), probably owing to dust in the broad-line emission region. To the best of our knowledge, this is the first case where the lag between the change in continuum and in broad-line emission of a "changing-look" AGN has been temporally resolved. The nature and timescales of the photometric and spectral evolution disfavor both a change in line-of-sight obscuration and a change of the overall rate of gas inflow as driving the drastic spectral transformations seen in this AGN. Although the peak luminosity and timescales are consistent with those of tidal disruption events seen in inactive galaxies, the spectral properties are not. The X-ray emission displays a markedly different behavior, with frequent flares on timescales of hours to days, and will be presented in a companion publication.

We report the discovery of an exceptional MIR flare in a Type 2 AGN, SDSS J165726.81+234528.1, at $z=0.059$. This object brightened by 3 mag in the Wide-field Infrared Survey Explorer (WISE) $W1$ and $W2$ bands between 2015 and 2017 (and is fading since 2018), without significant changes ($\lesssim$ 0.2 mag) in the optical over the same period of time. Based on the WISE light curves and near-IR imaging, the flare is more significant at longer wavelengths, suggesting an origin of hot dust emission. The estimated black hole mass ($\sim 10^{6.5}\,M_\odot$) from different methods places its peak bolometric luminosity around the Eddington limit. The high luminosity of the MIR flare and its multi-year timescale suggest that it most likely originated from reprocessed dust radiation in an extended torus surrounding the AGN, instead of from stellar explosions. The MIR color variability is consistent with known changing-look AGN and tidal disruption events (TDEs), but inconsistent with normal supernovae. We suggest that it is a turning-on Type 2 AGN or TDE, where the optical variability is obscured by the dust torus during the transition. This MIR flare event reveals a population of dramatic nuclear transients that are missed in the optical.

We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of $M_r=-13$ (${\nu}L_{\nu}=9\times10^6~L_{\odot}$), in the luminosity gap between novae and typical supernovae (SNe). Spectra during the outburst show a red continuum, Balmer emission with a velocity width of $\approx400$ km s$^{-1}$, Ca II and [Ca II] emission, and absorption features characteristic of an F-type supergiant. The spectra and multiband light curves are similar to the so-called "SN impostors" and intermediate-luminosity red transients (ILRTs). We directly identify the likely progenitor in archival Spitzer Space Telescope imaging with a $4.5~\mu$m luminosity of $M_{[4.5]}\approx-12.2$ and a $[3.6]-[4.5]$ color redder than 0.74 mag, similar to those of the prototype ILRTs SN 2008S and NGC 300 OT2008-1. Intensive monitoring of M51 with Spitzer further reveals evidence for variability of the progenitor candidate at [4.5] in the years before the OT. The progenitor is not detected in pre-outburst Hubble Space Telescope optical and near-IR images. The optical colors during outburst combined with spectroscopic temperature constraints imply a higher reddening of $E(B-V)\approx0.7$ mag and higher intrinsic luminosity of $M_r\approx-14.9$ (${\nu}L_{\nu}=5.3\times10^7~L_{\odot}$) near peak than seen in previous ILRT candidates. Moreover, the extinction estimate is higher on the rise than on the plateau, suggestive of an extended phase of circumstellar dust destruction. These results, enabled by the early discovery of M51 OT2019-1 and extensive pre-outburst archival coverage, offer new clues about the debated origins of ILRTs and may challenge the hypothesis that they arise from the electron-capture induced collapse of extreme asymptotic giant branch stars.

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.

We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

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 present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement existing sets of "Changing-Look" AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$\beta$/[OIII] ratio. We find 36 sources with declining H$\beta$ flux, 37 sources with increasing H$\beta$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$\beta$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk.

We present a novel algorithm for scheduling the observations of time-domain imaging surveys. Our Integer Linear Programming approach optimizes an observing plan for an entire night by assigning targets to temporal blocks, enabling strict control of the number of exposures obtained per field and minimizing filter changes. A subsequent optimization step minimizes slew times between each observation. Our optimization metric self-consistently weights contributions from time-varying airmass, seeing, and sky brightness to maximize the transient discovery rate. We describe the implementation of this algorithm on the surveys of the Zwicky Transient Facility and present its on-sky performance.

One clear observational prediction of the single degenerate progenitor scenario as the origin of type Ia supernovae (SNe) is the presence of relatively narrow ($\approx$1000 km s$^{-1}$) H$\alpha$ emission at nebular phases, although this feature is rarely seen. We present a compilation of nebular phase H$\alpha$ limits for SN Ia in the literature and demonstrate that this heterogenous sample has been biased towards SN Ia with relatively high luminosities and slow decline rates, as parameterized by $\Delta$m$_{15}(B)$, the difference in $B$-band magnitude between maximum light and fifteen days afterward. Motivated by the need to explore the full parameter space of SN~Ia and their subtypes, we present two new and six previously published nebular spectra of SN Ia with $\Delta$m$_{15}(B)$$ > $1.3 mag (including members of the transitional and SN1991bg-like subclasses) and measure nondetection limits of $L_{H\alpha}$$~<~$0.85--9.9$\times$10$^{36}$ ergs s$^{-1}$, which we confirmed by implanting simulated H$\alpha$ emission into our data. Based on the lastest models of swept-up material stripped from a nondegenerate companion star, these $L_{H\alpha}$ values correspond to hydrogen mass limits of $M_H$$~\lesssim~$1-3$\times$10$^{-4}$ $M_{\odot}$, roughly three orders of magnitude below that expected for the systems modeled, although we note that no simulations of H$\alpha$ nebular emission in such weak explosions have yet been performed. Despite the recent detection of strong H$\alpha$ in ASASSN-18tb (SN 2018fhw; $\Delta$m$_{15}(B)$ = 2.0 mag), we see no evidence that fast declining systems are more likely to have late time H$\alpha$ emission, although a larger sample is needed to confirm this result.

We identify minimal observing cadence requirements that enable photometric astronomical surveys to detect and recognize fast and explosive transients and fast transient features. Observations in two different filters within a short time window (e.g., g-and-i, or r-and-z, within < 0.5 hr) and a repeat of one of those filters with a longer time window (e.g., > 1.5 hr) are desirable for this purpose. Such an observing strategy delivers both the color and light curve evolution of transients on the same night. This allows the identification and initial characterization of fast transient -- or fast features of longer timescale transients -- such as rapidly declining supernovae, kilonovae, and the signatures of SN ejecta interacting with binary companion stars or circumstellar material. Some of these extragalactic transients are intrinsically rare and generally all hard to find, thus upcoming surveys like the Large Synoptic Survey Telescope (LSST) could dramatically improve our understanding of their origin and properties. We colloquially refer to such a strategy implementation for the LSST as the Presto-Color strategy (rapid-color). This cadence's minimal requirements allow for overall optimization of a survey for other science goals.

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.

Whether supernovae are a significant source of dust has been a long-standing debate. The large quantities of dust observed in high-redshift galaxies raise a fundamental question as to the origin of dust in the Universe since stars cannot have evolved to the AGB dust-producing phase in high-redshift galaxies. In contrast, supernovae occur within several millions of years after the onset of star formation. This white paper focuses on dust formation in supernova ejecta with US-Extremely Large Telescope (ELT) perspective during the era of JWST and LSST.

This White Pape motivates the time domain extragalactic science case for the NASA Near-Earth Object Camera (NEOCam). NEOCam is a NASA Planetary mission whose goal is to discover and characterize asteroids and comets, to assess the hazard to Earth from near-Earth objects, and to study the origin, evolution, and fate of asteroids and comets. NEOCam will, however, cover 68% of the extragalactic sky and as the NEOWISE-R mission has recently proved, infrared information is now vital for identifying and characterizing the $\gtrsim$10 million IR bright Active Galactic Nuclei, as well as using the IR light curve to provide deep insights into accretion disk astrophysics. NEOWISE-R data has also been used to discover Super-luminous Supernovae, dust echos in Tidal Disruption Events and detects all of the known $z\geq7$ quasars (and over 80% of the known $z\geq6.70$ quasars). As such, for relatively little additional cost, adding the capacity for additional NEOCam data processing (and/or alerting) would have a massive scientific and legacy impact on extragalactic time domain science.

The Large Synoptic Survey Telescope (LSST) will open a discovery frontier for faint and fast transients with its ability to detect variable flux components down to $\sim$24.5 mag in a $\sim$30 second exposure. Spectroscopic follow-up of such phenomena - which are necessary for understanding the physics of stellar explosions - can require a rapid response and several hours with a 8-10m telescope, making it both expensive and difficult to acquire. The future Extremely Large Telescopes (ELTs) would be able to provide not only spectroscopy but capabilities such as spectropolarimetry and high-resolution diffraction-limited imaging that would contribute to future advances in our physical understanding of stellar explosions. In this white paper we focus on several specific scientific impacts in the field of explosive transient astrophysics that will be generated by the combination of LSST's discovery abilities and ELTs' follow-up capacities. First, we map the uncharted frontier of discovery phase-space in terms of intrinsic luminosity and timescales for explosive transients, where we expect the unexpected. We then focus on six areas with open science questions for known transients: the progenitors of thermonuclear supernovae (SNe), mass loss prior to core collapse, asymmetry in stellar explosions, light echoes, high-$z$ transients, and strongly lensed SNe. We conclude with a brief discussion of the practical aspects of ELT & LSST synergy.

Supernovae are among the most powerful and influential explosions in the universe. They are also ideal multi-messenger laboratories to study extreme astrophysics. However, many fundamental properties of supernovae related to their diverse progenitor systems and explosion mechanisms remain poorly constrained. Here we outline how late-time spectroscopic observations obtained during the nebular phase (several months to years after explosion), made possible with the next generation of Extremely Large Telescopes, will facilitate transformational science opportunities and rapidly accelerate the community towards our goal of achieving a complete understanding of supernova explosions. We highlight specific examples of how complementary GMT and TMT instrumentation will enable high fidelity spectroscopy from which the line profiles and luminosities of elements tracing mass loss and ejecta can be used to extract kinematic and chemical information with unprecedented detail, for hundreds of objects. This will provide uniquely powerful constraints on the evolutionary phases stars may experience approaching a supernova explosion; the subsequent explosion dynamics; their nucleosynthesis yields; and the formation of compact objects that may act as central engines.

The non-detection of companion stars in Type Ia supernova (SN) progenitor systems lends support to the notion of double-degenerate (DD) systems and explosions triggered by the merging of two white dwarfs. This very asymmetric process should lead to a conspicuous polarimetric signature. By contrast, observations consistently find very low continuum polarization as the signatures from the explosion process largely dominate over the pre-explosion configuration within several days. Critical information about the interaction of the ejecta with a companion and any circumstellar matter is encoded in the early polarization spectra. In this study, we obtain spectropolarimetry of SN\,2018gv with the ESO Very Large Telescope at $-$13.6 days relative to the $B-$band maximum light, or $\sim$5 days after the estimated explosion --- the earliest spectropolarimetric observations to date of any Type Ia SN. These early observations still show a low continuum polarization ($\lesssim$0.2\%) and moderate line polarization (0.30$\pm$0.04\% for the prominent \ion{Si}{2} $\lambda$6355 feature and 0.85$\pm$0.04\% for the high-velocity Ca component). The high degree of spherical symmetry implied by the low line and continuum polarization at this early epoch is consistent with explosion models of delayed detonations and is inconsistent with the merger-induced explosion scenario. The dense UV and optical photometry and optical spectroscopy within the first $\sim$100 days after the maximum light indicate that SN\,2018gv is a normal Type Ia SN with similar spectrophotometric behavior to SN\,2011fe.

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

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

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

We report the discovery of non-stellar hydrogen Balmer and metastable helium absorption lines accompanying a transient, high-velocity (0.05$c$) broad absorption line (BAL) system in the optical spectra of the tidal disruption event (TDE) AT2018zr ($z=0.071$). In the HST UV spectra, absorption of high- and low-ionization lines are also present at this velocity, making AT2018zr resemble a low-ionization broad absorption line (LoBAL) QSO. We conclude that these transient absorption features are more likely to arise in fast outflows produced by the TDE than absorbed by the unbound debris. In accordance with the outflow picture, we are able to reproduce the flat-topped H$\alpha$ emission in a spherically expanding medium, without invoking the typical prescription of an elliptical disk. We also report the appearance of narrow ($\sim$1000~km~s$^{-1}$) NIII$\lambda$4640, HeII$\lambda4686$, H$\alpha$, and H$\beta$, emission in the late-time optical spectra of AT2018zr, which may be a result of UV continuum hardening at late time as observed by Swift. Including AT2018zr, we find a high association rate (3 out of 4) of BALs in the UV spectra of TDEs. This suggests that outflows may be ubiquitous among TDEs and may be less sensitive to viewing angle effects compared to QSO outflows.

Over the past decade, research in resolved stellar populations has made great strides in exploring the nature of dark matter, in unraveling the star formation, chemical enrichment, and dynamical histories of the Milky Way and nearby galaxies, and in probing fundamental physics from general relativity to the structure of stars. Large surveys have been particularly important to the biggest of these discoveries. In the coming decade, current and planned surveys will push these research areas still further through a large variety of discovery spaces, giving us unprecedented views into the low surface brightness Universe, the high surface brightness Universe, the 3D motions of stars, the time domain, and the chemical abundances of stellar populations. These discovery spaces will be opened by a diverse range of facilities, including the continuing Gaia mission, imaging machines like LSST and WFIRST, massively multiplexed spectroscopic platforms like DESI, Subaru-PFS, and MSE, and telescopes with high sensitivity and spatial resolution like JWST, the ELTs, and LUVOIR. We do not know which of these facilities will prove most critical for resolved stellar populations research in the next decade. We can predict, however, that their chance of success will be maximized by granting use of the data to broad communities, that many scientific discoveries will draw on a combination of data from them, and that advances in computing will enable increasingly sophisticated analyses of the large and complex datasets that they will produce. We recommend that Astro2020 1) acknowledge the critical role that data archives will play for stellar populations and other science in the next decade, 2) recognize the opportunity that advances in computing will bring for survey data analysis, and 3) consider investments in Science Platform technology to bring these opportunities to fruition.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

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.

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST.

Active galactic nuclei (AGN) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as "changing-look" AGN, or at higher luminosities changing look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies, and open a wider interpretation of the cause of differences between broad and narrow line AGN. To better characterize extreme quasar variability, we present follow-up spectroscopy as part of a comprehensive search for CLQs across the full SDSS footprint using spectroscopically confirmed quasars from the SDSS DR7 catalog. Our primary selection requires large-amplitude (|\Delta g|>1 mag, |\Delta r|>0.5 mag) variability over any of the available time baselines probed by the SDSS and Pan-STARRS 1 surveys. We employ photometry from the Catalina Sky Survey to verify variability behavior in CLQ candidates where available, and confirm CLQs using optical spectroscopy from the William Herschel, MMT, Magellan, and Palomar telescopes. For our adopted S/N threshold on variability of broad H\beta emission, we find 17 new CLQs, yielding a confirmation rate of >~ 20%. These candidates are at lower Eddington ratio relative to the overall quasar population which supports a disk-wind model for the broad line region. Based on our sample, the CLQ fraction increases from 10% to roughly half as the continuum flux ratio between repeat spectra at 3420 Angstroms increases from 1.5 to 6. We release a catalog of over 200 highly variable candidates to facilitate future CLQ searches.

The Zwicky Transient Facility (ZTF) survey generates real-time alerts for optical transients, variables, and moving objects discovered in its wide-field survey. We describe the ZTF alert stream distribution and processing (filtering) system. The system uses existing open-source technologies developed in industry: Kafka, a real-time streaming platform, and Avro, a binary serialization format. The technologies used in this system provide a number of advantages for the ZTF use case, including (1) built-in replication, scalability, and stream rewind for the distribution mechanism; (2) structured messages with strictly enforced schemas and dynamic typing for fast parsing; and (3) a Python-based stream processing interface that is similar to batch for a familiar and user-friendly plug-in filter system, all in a modular, primarily containerized system. The production deployment has successfully supported streaming up to 1.2 million alerts or roughly 70 GB of data per night, with each alert available to a consumer within about 10 s of alert candidate production. Data transfer rates of about 80,000 alerts/minute have been observed. In this paper, we discuss this alert distribution and processing system, the design motivations for the technology choices for the framework, performance in production, and how this system may be generally suitable for other alert stream use cases, including the upcoming Large Synoptic Survey Telescope.

How an otherwise inert carbon-oxygen white dwarf can be made to explode as a Type Ia supernova remains unknown. A promising test of theoretical models is to constrain the distribution of material that is left unburned, in particular of carbon. So far, most investigations used line identification codes to detect carbon in the ejecta, a method that cannot be readily compared against model predictions because it requires assumed opacities and temperatures. Here, we instead use tomographic techniques to investigate the amount of carbon in the inner layers of SN~2011fe, starting from the previously published tomographic analysis of Mazzali et al. From the presence of the carbon feature in the optical at early epochs and its disappearance later on, we derive an average carbon mass fraction between 0.001 and 0.05 for velocities in the range $13500 \lesssim v \lesssim 16000\ \rm{km\ s^{-1}}$, and an upper limit of 0.005 inside that region. Based on our models and the assumed density profile, only small amounts of carbon should be in the neutral state, too little to be responsible for features seen in near-infrared spectra that were previously identified as due to neutral carbon; We discuss possible reasons for this discrepancy. We compare our results against a suite of explosion models, although uncertainties in both the models and our simulations make it difficult to draw definitive conclusions.

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

This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics.

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.

Aims: In this work, we aim to provide a reliable list of gravitational lens (GL) candidates based on a search performed over the entire Gaia Data Release 2 (Gaia DR2). We also show that the sole astrometric and photometric informations coming from the Gaia satellite yield sufficient insights for supervised learning methods to automatically identify GL candidates with an efficiency that is comparable to methods based on image processing. Methods: We simulated 106,623,188 lens systems composed of more than two images, based on a regular grid of parameters characterizing a non-singular isothermal ellipsoid lens model in the presence of an external shear. These simulations are used as an input for training and testing our supervised learning models consisting of Extremely Randomized Trees. The latter are finally used to assign to each of the 2,129,659 clusters of celestial objects a discriminant value that reflects the ability of our simulations to match the observed relative positions and fluxes from each cluster. Once complemented with additional constraints, these discriminant values allowed us to identify GL candidates out of the list of clusters. Results: We report the discovery of 15 new quadruply-imaged lens candidates with angular separations less than 6" and assess the performance of our approach by recovering 12 out of the 13 known quadruply-imaged systems with all their components detected in Gaia DR2 with a misclassification rate of fortuitous clusters of stars as lens systems that is below one percent. Similarly, the identification capability of our method regarding quadruply-imaged systems where three images are detected in Gaia DR2 is assessed by recovering 10 out of the 13 known quadruply-imaged systems having one of their constituting images discarded. The associated misclassification rate varying then between 5.8% and 20%, depending on the image we decided to remove.

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.

Cosmology is one of the four science pillars of LSST, which promises to be transformative for our understanding of dark energy and dark matter. The LSST Dark Energy Science Collaboration (DESC) has been tasked with deriving constraints on cosmological parameters from LSST data. Each of the cosmological probes for LSST is heavily impacted by the choice of observing strategy. This white paper is written by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy that will benefit all cosmological analyses with LSST. It is accompanied by the DESC DDF (Deep Drilling Fields) white paper (Scolnic et al.). We use a variety of metrics to understand the effects of the observing strategy on measurements of weak lensing, large-scale structure, clusters, photometric redshifts, supernovae, strong lensing and kilonovae. In order to reduce systematic uncertainties, we conclude that the current baseline observing strategy needs to be significantly modified to result in the best possible cosmological constraints. We provide some key recommendations: moving the WFD (Wide-Fast-Deep) footprint to avoid regions of high extinction, taking visit pairs in different filters, changing the 2x15s snaps to a single exposure to improve efficiency, focusing on strategies that reduce long gaps (>15 days) between observations, and prioritizing spatial uniformity at several intervals during the 10-year survey.

The nature and role of the binary companion of carbon-oxygen white dwarf stars that explode as Type Ia supernovae (SNe Ia) are not yet fully understood. Past detections of circumstellar material (CSM) that contain hydrogen for a small number of SN Ia progenitor systems suggest that at least some have a nondegenerate companion. In order to constrain the prevalence, location, and quantity of CSM in SN Ia systems, we performed a near-ultraviolet (NUV) survey with the Hubble Space Telescope (HST) to look for the high-energy signature of SN Ia ejecta interacting with CSM. Our survey revealed that SN 2015cp, a SN 1991T-like overluminous SN Ia, was experiencing late-onset interaction between its ejecta and surrounding CSM at $664$ days after its light-curve peak. We present ground- and space-based follow-up observations of SN 2015cp that reveal optical emission lines of H and Ca, typical signatures of ejecta-CSM interaction. We show how SN 2015cp was likely similar to the well-studied SN Ia-CSM event PTF11kx, making it the second case in which an unambiguously classified SN Ia was observed to interact with a distant shell of CSM that contains hydrogen ($R_{\rm CSM} \gtrsim 10^{16}\ {\rm cm}$). The remainder of our HST NUV images of SNe Ia were nondetections that we use to constrain the occurrence rate of observable late-onset CSM interaction. We apply theoretical models for the emission from ejecta-CSM interaction to our NUV nondetections, and place upper limits on the mass and radial extent of CSM in SN Ia progenitor systems.

Despite their cosmological utility, the progenitors of Type Ia supernovae (SNe Ia) are still unknown, with many efforts focused on whether accretion from a nondegenerate companion can grow a carbon-oxygen white dwarf to near the Chandrasekhar mass. The association of SNe Ia resembling SN 1991T ("91T-like") with circumstellar interaction may be evidence for this "single-degenerate" channel. However, the observed circumstellar medium (CSM) in these interacting systems is unlike a stellar wind -- of particular interest, it is sometimes detached from the stellar surface, residing at $\sim 10^{16}~{\rm cm}$. A Hubble Space Telescope (HST) program to discover detached CSM around 91T-like SNe Ia successfully discovered interaction nearly two years after explosion in SN 2015cp (Graham et al., 2018). In this work, we present radio and X-ray follow-up observations of SN 2015cp and analyze them in the framework of Harris, Nugent, & Kasen (2016) to limit the properties of a constant-density CSM shell in this system. Assuming the HST detection was shortly after the shock crossed the CSM, we constrain the total CSM mass in this system to be $< 0.5~{\rm M_\odot}$. This limit is comparable to the CSM mass of supernova PTF11kx, but does not rule out lower masses predicted for recurrent novae. From lessons learned modeling PTF11kx and SN 2015cp, we suggest a strategy for future observations of these events to increase the sample of known interacting SNe Ia.

Recently Assef et al.(2018) presented two catalogs of AGN candidates over 30,093 deg^2 selected from the Wide-field Infrared Survey Explorer (WISE) observations. From their most reliable sample, Assef et al. (2018) identified 45 AGN candidates with the highest variability levels in the AllWISE catalog, but that are not blazars. Here we present new spectroscopic observations of some of these targets to further constrain their nature. We also study their optical lightcurves using observations from CRTS, and find that only seven show significant optical variability, and that five of those seven are spectroscopically classified as AGN. In one of them, WISEA J094806.56+031801.7 (W0948+0318), we identify a transient event in the CRTS lightcurve. We present a detailed analysis of this transient, and characterize it through its CRTS lightcurve and its multi-wavelength spectral energy distribution obtained from GALEX, Pan-STARRS and WISE observations. We find that the most likely source of the transient is a super-luminous supernova (SLSN) in W0948+0318. We estimate the total radiated energy to be E=1.6\pm 0.3 x 10^52 erg, making it one of the most energetic SLSN observed. Based on the lack of change in mid-IR color throughout and after the transient event, we speculate that the location of the SLSN is within the torus of the AGN. We identify 9 possible analogs to W0948+0318 based on their WISE lightcurves. None show optically detected transients and hence suggest significant dust obscuration. Finally, we estimate a rate of >2x10^-7 yr^-1 per AGN for these transients under the conservative assumption that none of the identified analogs have a common origin with the transient in W0948+0318.

We present spectroscopy from the first three seasons of the Dark Energy Survey Supernova Program (DES-SN). We describe the supernova spectroscopic program in full: strategy, observations, data reduction, and classification. We have spectroscopically confirmed 307 supernovae, including 251 type Ia supernovae (SNe Ia) over a redshift range of $0.017 < z < 0.85$. We determine the effective spectroscopic selection function for our sample, and use it to investigate the redshift-dependent bias on the distance moduli of SNe Ia we have classified. We also provide a full overview of the strategy, observations, and data products of DES-SN, which has discovered 12,015 likely supernovae during these first three seasons. The data presented here are used for the first cosmology analysis by DES-SN ('DES-SN3YR'), the results of which are given in DES Collaboration (2018a).

We present the complete sample of stripped-envelope supernova (SN) spectra observed by the Lick Observatory Supernova Search (LOSS) collaboration over the last three decades: 888 spectra of 302 SNe, 652 published here for the first time, with 384 spectra (of 92 SNe) having photometrically-determined phases. After correcting for redshift and Milky Way dust reddening and reevaluating the spectroscopic classifications for each SN, we construct mean spectra of the three major spectral subtypes (Types IIb, Ib, and Ic) binned by phase. We compare measures of line strengths and widths made from this sample to the results of previous efforts, confirming that O I {\lambda}7774 absorption is stronger and found at higher velocity in Type Ic SNe than in Types Ib or IIb SNe in the first 30 days after peak brightness, though the widths of nebular emission lines are consistent across subtypes. We also highlight newly available observations for a few rare subpopulations of interest.

We report the spectroscopic confirmation and modeling of the quadruply imaged quasar GRAL113100-441959, the first gravitational lens (GL) to be discovered mainly from astrometric considerations. Follow-up spectra obtained with Keck/LRIS reveal the lensing nature of this quadruply-imaged quasar with redshift $z_s = 1.090 \pm 0.002$, but show no evidence of the central lens galaxy. Using the image positions and $G$-band flux ratios provided by Gaia Data Release 2 as constraints, we model the system with a singular power-law elliptical mass distribution (SPEMD) plus external shear, to different levels of complexity. We show that relaxing the isothermal constraint of the SPEMD is not statistically significant, and thus we simplify the SPEMD to a singular isothermal ellipsoid to estimate the Einstein radius of the main lens galaxy $\theta_{\text{E}} = 0.\!\!^{\prime\prime}851$, the intensity and position angle of the external shear $(\gamma,\theta_{\gamma}) = (0.044,11.\!\!^{\circ}5)$, and we predict the lensing galaxy position to be $(x_{\text{gal}},y_{\text{gal}}) = (-0.\!\!^{\prime\prime}424,-0.\!\!^{\prime\prime}744)$ with respect to image A. We provide time delay predictions for pairs of images, assuming a plausible range of lens redshift values $z_l$ between $0.5$ and $0.9$. We finally examine the impact on time delays of the so-called Source Position Transformation, a family of degeneracies existing between different lens density profiles that reproduce most of the lensing observables equally well. We show that this effect contributes significantly to the time delay error budget and cannot be ignored during the modeling. This has implications for robust cosmography applications of lensed systems. GRAL113100-441959 is the first in a series of seven new spectroscopically confirmed GLs discovered from Gaia Data Release 2.

We present an analysis of nine years of Catalina Surveys optical photometry for 1.4 million spectroscopically confirmed SDSS galaxies. We find 717 outburst events that were not reported by ongoing transient surveys. These events have timescales ranging from weeks to years. More than two thirds of these new events are found in starforming galaxies, while such galaxies only constitute ~20% of our sample. Based on the properties of the hosts and events, we find that almost all of the new events are likely to be associated with regular supernovae. However, a small number of long-timescale events are found among the galaxies containing AGN. These events have similar properties to those recently found in the analyses of light curves of large samples of AGN. Given the lack of such events among the more than a million passive galaxies in the sample, we suggest that the long outbursts are associated with super-massive black holes or their environments.

Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear more likely than changes in obscuration. Here we report on four epochs of spectroscopy of SDSS J110057.70-005304.5, a quasar at a redshift of $z=0.378$ whose UV continuum and broad hydrogen emission lines have faded, and then returned over the past $\approx$20 years. The change in this quasar was initially identified in the infrared, and an archival spectrum from 2010 shows an intermediate phase of the transition during which the flux below rest-frame $\approx$3400\AA\ has decreased by close to an order of magnitude. This combination is unique compared to previously published examples of changing-look quasars, and is best explained by dramatic changes in the innermost regions of the accretion disk. The optical continuum has been rising since mid-2016, leading to a prediction of a rise in hydrogen emission line flux in the next year. Increases in the infrared flux are beginning to follow, delayed by a $\sim$3 year observed timescale. If our model is confirmed, the physics of changing-look quasars are governed by processes at the innermost stable circular orbit (ISCO) around the black hole, and the structure of the innermost disk. The easily identifiable and monitored changing-look quasars would then provide a new probe and laboratory of the nuclear central engine.

By inverting the distributions of galaxies' apparent ellipticities and misalignment angles (measured around the projected half-light radius $R_{\rm e}$) between their photometric and kinematic axes, we study the intrinsic shape distribution of 189 slow rotator early-type galaxies with stellar masses $2\times 10^{11} M_{\odot}<M_\ast<2\times 10^{12} M_{\odot}$, extracted from a sample of about 2200 galaxies with integral-field stellar kinematics from the DR14 of the SDSS-IV MaNGA IFU survey. Thanks to the large sample of slow rotators, Graham+18 showed that there is clear structure in the misalignment angle distribution, with two peaks at both $0^{\circ}$ and $90^{\circ}$ misalignment (characteristic of oblate and prolate rotation respectively). Here we invert the observed distribution from Graham+18. The large sample allows us to go beyond the known fact that slow rotators are weakly triaxial and to place useful constraints on their intrinsic triaxiality distribution (around $1R_{\rm e}$) for the first time. The shape inversion is generally non-unique. However, we find that, for a wide set of model assumptions, the observed distribution clearly requires a dominant triaxial-oblate population. For some of our models, the data suggest a hint for a minor triaxial-prolate population, but a dominant prolate population is ruled out.

We present the discovery of a gravitationally lensed dust-reddened QSO at $z=2.517$ discovered in a survey for red QSOs by infrared selection. $Hubble~Space~Telescope$ imaging in the WFC3/IR F160W and F125W filters reveals a quadruply lensed system in a cusp configuration. We find that compared to the central image of the cusp, the nearby, brightest image is anomalous by a factor of $\sim7-11$. Although the source is extremely bright in the mid-infrared, a magnification by a factor of $\sim50-120$ places an upper limit of 1.35 mJy on the intrinsic mid-infrared brightness, well below the $WISE~W4$ detection limit of 6 mJy. We find that this QSO is moderately reddened, with $E(B-V)=0.7$ and that $\sim1\%$ of the intrinsic spectrum is leaked back into the line of sight resulting in an upturn in its UV spectrum. We conclude that the QSO's reddening is intrinsic and not due to the lens. Consistent with previous red quasar samples, this source exhibits outflows in its spectrum as well as morphological properties suggestive of it being in a merger-driven transitional phase. Depending on how $L_{\rm bol}$ is computed, the quasar's accretion rate may be as high as $0.26~L_{\rm Edd}$. We detect two Lyman limit systems, at $z=2.102$ and $z=2.431$, with absorption by metal lines likely at small impact parameter to the QSO, and a putative lens redshift of $z=0.599$. Given the rarity of quad lenses, the discovery of this source allows detailed study of a less luminous, more typical infrared-selected quasar at high redshift.

We present an analysis of 507 spectra of 173 stripped-envelope (SE) supernovae (SNe) discovered by the untargeted Palomar Transient Factory (PTF) and intermediate PTF (iPTF) surveys. Our sample contains 55 Type IIb SNe (SNe IIb), 45 Type Ib SNe (SNe Ib), 56 Type Ic SNe (SNe Ic), and 17 Type Ib/c SNe (SNe Ib/c). We compare the SE SN subtypes via measurements of the pseudo-equivalent widths (pEWs) and velocities of the He I $\lambda\lambda5876, 7065$ and O I $\lambda7774$ absorption lines. Consistent with previous work, we find that SNe Ic show higher pEWs and velocities in O I $\lambda7774$ compared to SNe IIb and Ib. The pEWs of the He I $\lambda\lambda5876, 7065$ lines are similar in SNe Ib and IIb after maximum light. The He I $\lambda\lambda5876, 7065$ velocities at maximum light are higher in SNe Ib compared to SNe IIb. We have identified an anticorrelation between the He I $\lambda7065$ pEW and O I $\lambda7774$ velocity among SNe IIb and Ib. This can be interpreted as a continuum in the amount of He present at the time of explosion. It has been suggested that SNe Ib and Ic have similar amounts of He, and that lower mixing could be responsible for hiding He in SNe Ic. However, our data contradict this mixing hypothesis. The observed difference in the expansion rate of the ejecta around maximum light of SNe Ic ($V_{\mathrm{m}}=\sqrt{2E_{\mathrm{k}}/M_{\mathrm{ej}}}\approx15,000$ km s$^{-1}$) and SNe Ib ($V_{\mathrm{m}}\approx9000$ km s$^{-1}$) would imply an average He mass difference of $\sim1.4$ $M_{\odot}$, if the other explosion parameters are assumed to be unchanged between the SE SN subtypes. We conclude that SNe Ic do not hide He but lose He due to envelope stripping.

The Type Ia supernovae (SNe Ia) 2011by, hosted in NGC 3972, and 2011fe, hosted in M101, are optical "twins," having almost identical optical light-curve shapes, colours, and near-maximum-brightness spectra. However, SN 2011fe had significantly more ultraviolet (UV; 1600 < lambda < 2500 A) flux than SN 2011by before and at peak luminosity. Theory suggests that SNe Ia with higher progenitor metallicity should (1) have additional UV opacity near peak and thus lower UV flux; (2) have an essentially unchanged optical spectral-energy distribution; (3) have a similar optical light-curve shape; and (4) because of the excess neutrons, produce more stable Fe-group elements at the expense of radioactive 56Ni and thus have a lower peak luminosity. Foley & Kirshner (2013) suggested that the difference in UV flux between SNe 2011by and 2011fe was the result of their progenitors having significantly different metallicities. The SNe also had a large, but insignificant, difference between their peak absolute magnitudes (Delta M_V, peak = 0.60 +/- 0.36 mag), with SN 2011fe being more luminous. We present a new Cepheid-based distance to NGC 3972, significantly improving the precision of the distance measurement for SN 2011by. With these new data, we determine that the SNe have significantly different peak luminosities (Delta M_V, peak = 0.335 +/- 0.069 mag), corresponding to SN 2011fe having produced 38% more 56Ni than SN 2011by, and providing additional evidence for progenitor metallicity differences for these SNe. We discuss how progenitor metallicity differences can contribute to the intrinsic scatter for light-curve-shape-corrected SN luminosities, the use of "twin" SNe for measuring distances, and implications for using SNe Ia for constraining cosmological parameters.

We present nebular phase optical and near-infrared spectroscopy of the Type Ia supernova (SN) 2017cbv. The early light curves of SN~2017cbv showed a prominent blue bump in the $U$, $B$ and $g$ bands lasting for $\sim$5 d. One interpretation of the early light curve was that the excess blue light was due to shocking of the SN ejecta against a nondegenerate companion star -- a signature of the single degenerate scenario. If this is the correct interpretation, the interaction between the SN ejecta and the companion star could result in significant H$\alpha$ (or helium) emission at late times, possibly along with other species, depending on the companion star and its orbital separation. A search for H$\alpha$ emission in our +302 d spectrum yields a nondetection, with a $L_{H\alpha}$$<$8.0$\times$10$^{35}$ erg/s (given an assumed distance of $D$=12.3 Mpc), which we have verified by implanting simulated H$\alpha$ emission into our data. We make a quantitative comparison to models of swept-up material stripped from a nondegenerate companion star, and limit the mass of hydrogen that might remain undetected to $M_{\rm H} < 1 \times 10^{-4}$ $\rm M_{\odot}$. A similar analysis of helium star related lines yields a $M_{\rm He} < 5 \times 10^{-4}$ $\rm M_{\odot}$. Taken at face value, these results argue against a nondegenerate H or He-rich companion in Roche lobe overflow as the progenitor of SN 2017cbv. Alternatively, there could be weaknesses in the envelope-stripping and radiative transfer models necessary to interpret the strong H and He flux limits.

Modern galaxy surveys produce redshift probability density functions (PDFs) in addition to traditional photometric redshift (photo-$z$) point estimates. However, the storage of photo-$z$ PDFs may present a challenge with increasingly large catalogs, as we face a trade-off between the accuracy of subsequent science measurements and the limitation of finite storage resources. This paper presents $\texttt{qp}$, a Python package for manipulating parametrizations of 1-dimensional PDFs, as suitable for photo-$z$ PDF compression. We use $\texttt{qp}$ to investigate the performance of three simple PDF storage formats (quantiles, samples, and step functions) as a function of the number of stored parameters on two realistic mock datasets, representative of upcoming surveys with different data qualities. We propose some best practices for choosing a photo-$z$ PDF approximation scheme and demonstrate the approach on a science case using performance metrics on both ensembles of individual photo-$z$ PDFs and an estimator of the overall redshift distribution function. We show that both the properties of the set of PDFs we wish to approximate and the chosen fidelity metric(s) affect the optimal parametrization. Additionally, we find that quantiles and samples outperform step functions, and we encourage further consideration of these formats for PDF approximation.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

Thanks to its spatial resolution the ESA/Gaia space mission offers a unique opportunity to discover new multiply-imaged quasars and to study the already known lensed systems at sub-milliarcsecond astrometric precisions. In this paper, we address the detection of the known multiply-imaged quasars from the Gaia Data Release 2 and determine the astrometric and photometric properties of the individually detected images found in the Gaia DR2 catalogue. We have compiled an exhaustive list of quasar gravitational lenses from the literature to search for counterparts in the Gaia Data Release 2. We then analyze the astrometric and photometric properties of these Gaia's detections. To highlight the tremendous potential of Gaia at the sub-milliarcsecond level we finally perform a simple Bayesian modeling of the well-known gravitational lens system HE0435-1223, using Gaia Data Release 2 and HST astrometry. From 478 known multiply imaged quasars, 200 have at least one image found in the Gaia Data Release 2. Among the 41 known quadruply-imaged quasars of the list, 26 have at least one image in the Gaia Data Release 2, 12 of which are fully detected, 6 have three counterparts, 7 have two and 1 has only one. As expected, the modeling of HE0435-1223 shows that the model parameters are significantly better constrained when using Gaia astrometry compared to HST astrometry, in particular the relative positions of the background quasar source and the centroid of the deflector. The Gaia sub-milliarcsecond astrometry also significantly reduces the parameter correlations. Besides providing an up-to-date list of multiply imaged quasars and their detection in the Gaia DR2, this paper shows that more complex modeling scenarios will certainly benefit from Gaia sub-milliarcsecond astrometry.

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios.

We present a spectroscopically complete sample of 147 infrared-color-selected AGN down to a 22 $\mu$m flux limit of 20 mJy over the $\sim$270 deg$^2$ of the SDSS Stripe 82 region. Most of these sources are in the QSO luminosity regime ($L_{\rm bol} \gtrsim 10^{12} L_\odot$) and are found out to $z\simeq3$. We classify the AGN into three types, finding: 57 blue, unobscured Type-1 (broad-lined) sources; 69 obscured, Type-2 (narrow-lined) sources; and 21 moderately-reddened Type-1 sources (broad-lined and $E(B-V) > 0.25$). We study a subset of this sample in X-rays and analyze their obscuration to find that our spectroscopic classifications are in broad agreement with low, moderate, and large amounts of absorption for Type-1, red Type-1 and Type-2 AGN, respectively. We also investigate how their X-ray luminosities correlate with other known bolometric luminosity indicators such as [O III] line luminosity ($L_{\rm [OIII]}$) and infrared luminosity ($L_{6 \mu{\rm m}}$). While the X-ray correlation with $L_{\rm [OIII]}$ is consistent with previous findings, the most infrared-luminous sources appear to deviate from established relations such that they are either under-luminous in X-rays or over-luminous in the infrared. Finally, we examine the luminosity function (LF) evolution of our sample, and by AGN type, in combination with the complementary, infrared-selected, AGN sample of Lacy et al. (2013), spanning over two orders of magnitude in luminosity. We find that the two obscured populations evolve differently, with reddened Type-1 AGN dominating the obscured AGN fraction ($\sim$30%) for $L_{5 \mu{\rm m}} > 10^{45}$ erg s$^{-1}$, while the fraction of Type-2 AGN with $L_{5 \mu{\rm m}} < 10^{45}$ erg s$^{-1}$ rises sharply from 40% to 80% of the overall AGN population.

Context. Strong gravitationally lensed quasars are among the most interesting and useful observable extragalactic phenomena. Because their study constitutes a unique tool in various fields of astronomy, they are highly sought, not without difficulty. Indeed, even in this era of all-sky surveys, their recognition remains a great challenge, with barely a few hundred currently known systems. Aims. In this work we aim to detect new strongly lensed quasar candidates in the recently published Gaia Data Release 2 (DR2), which is the highest spatial resolution astrometric and photometric all-sky survey, attaining effective resolutions from 0.4" to 2.2". Methods. We cross-matched a merged list of quasars and candidates with the Gaia DR2 and found 1,839,143 counterparts within 0.5". We then searched matches with more than two Gaia DR2 counterparts within 6". We further narrowed the resulting list using astrometry and photometry compatibility criteria between the Gaia DR2 counterparts. A supervised machine learning method, Extremely Randomized Trees, is finally adopted to assign to each remaining system a probability of being lensed. Results. We report the discovery of three quadruply-imaged quasar candidates that are fully detected in Gaia DR2. These are the most promising new quasar lens candidates from Gaia DR2 and a simple singular isothermal ellipsoid lens model is able to reproduce their image positions to within $\sim$1 mas. This letter demonstrates the gravitational lens discovery potential of Gaia.

The dominant mechanism and time scales over which stripped-envelope supernovae (SNe) progenitor stars shed their hydrogen envelopes are uncertain. Observations of Type Ib and Ic SNe at late phases could reveal the optical signatures of interaction with distant circumstellar material (CSM) providing important clues on the origin of the necessary pre-SN mass loss. We report deep late-time optical spectroscopy of the Type Ib explosion SN 2004dk 4684 days (13 years) after discovery. Prominent intermediate-width H-alpha emission is detected, signaling that the SN blast wave has caught up with the hydrogen-rich CSM lost by the progenitor system. The line luminosity is the highest ever reported for a SN at this late stage. Prominent emission features of He, Fe, and Ca are also detected. The spectral characteristics are consistent with CSM energized by the forward shock, and resemble the late-time spectra of the persistently interacting Type IIn SNe 2005ip and 1988Z. We suggest that the onset of interaction with H-rich CSM was associated with a previously reported radio rebrightening at ~1700 days. The data indicate that the mode of pre-SN mass loss was a relatively slow dense wind that persisted millennia before the SN, followed by a short-lived Wolf-Rayet phase that preceded core-collapse and created a cavity within an extended distribution of CSM. We also present new spectra of SNe 2014C, PTF11iqb, and 2009ip, all of which also exhibit continued interaction with extended CSM distributions.

We measure $\lambda_{R_e}$, a proxy for galaxy specific stellar angular momentum within one effective radius, and the ellipticity, $\epsilon$, for about 2300 galaxies of all morphological types observed with integral field spectroscopy as part of the MaNGA survey, the largest such sample to date. We use the $(\lambda_{R_e}, \epsilon)$ diagram to separate early-type galaxies into fast and slow rotators. We also visually classify each galaxy according to its optical morphology and two-dimensional stellar velocity field. Comparing these classifications to quantitative $\lambda_{R_e}$ measurements reveals tight relationships between angular momentum and galaxy structure. In order to account for atmospheric seeing, we use realistic models of galaxy kinematics to derive a general approximate analytic correction for $\lambda_{R_e}$. Thanks to the size of the sample and the large number of massive galaxies, we unambiguously detect a clear bimodality in the $(\lambda_{R_e}, \epsilon)$ diagram which may result from fundamental differences in galaxy assembly history. There is a sharp secondary density peak inside the region of the diagram with low $\lambda_{R_e}$ and $\epsilon < 0.4$, previously suggested as the definition for slow rotators. Most of these galaxies are visually classified as non-regular rotators and have high velocity dispersion. The intrinsic bimodality must be stronger, as it tends to be smoothed by noise and inclination. The large sample of slow rotators allows us for the first time to unveil a secondary peak at +/-90 degrees in their distribution of the misalignments between the photometric and kinematic position angles. We confirm that genuine slow rotators start appearing above a stellar mass of 2\times10^{11} M_{\odot}$ where a significant number of high-mass fast rotators also exist.

We present photometry, spectra, and spectropolarimetry of supernova (SN) 2012ab, mostly obtained over the course of $\sim 300$ days after discovery. SN 2012ab was a Type IIn (SN IIn) event discovered near the nucleus of spiral galaxy 2MASXJ12224762+0536247. While its light curve resembles that of SN 1998S, its spectral evolution does not. We see indications of CSM interaction in the strong intermediate-width emission features, the high luminosity (peak at absolute magnitude $M=-19.5$), and the lack of broad absorption features in the spectrum. The H$\alpha$ emission undergoes a peculiar transition. At early times it shows a broad blue emission wing out to $-14{,}000$ km $\mathrm{s^{-1}}$ and a truncated red wing. Then at late times ($>$ 100$\,$days) it shows a truncated blue wing and a very broad red emission wing out to roughly $+20{,}000$ km $\mathrm{s^{-1}}$. This late-time broad red wing probably arises in the reverse shock. Spectra also show an asymmetric intermediate-width H$\alpha$ component with stronger emission on the red side at late times. The evolution of the asymmetric profiles requires a density structure in the distant CSM that is highly aspherical. Our spectropolarimetric data also suggest asphericity with a strong continuum polarization of $\sim 1-3$% and depolarization in the H$\alpha$ line, indicating asphericity in the CSM at a level comparable to that in other SNe IIn. We estimate a mass-loss rate of $\dot{M} = 0.050\, {\rm M}_{\odot}\,\mathrm{yr^{-1}}$ for $v_{\rm pre} = 100$$\,$km$\,$$\mathrm{s^{-1}}$ extending back at least 75$\,$yr prior to the SN. The strong departure from axisymmetry in the CSM of SN 2012ab may suggest that the progenitor was an eccentric binary system undergoing eruptive mass loss.

In this paper we present and characterize a nearest-neighbors color-matching photometric redshift estimator that features a direct relationship between the precision and accuracy of the input magnitudes and the output photometric redshifts. This aspect makes our estimator an ideal tool for evaluating the impact of changes to LSST survey parameters that affect the measurement errors of the photometry, which is the main motivation of our work (i.e., it is not intended to provide the "best" photometric redshifts for LSST data). We show how the photometric redshifts will improve with time over the 10-year LSST survey and confirm that the nominal distribution of visits per filter provides the most accurate photo-$z$ results. The LSST survey strategy naturally produces observations over a range of airmass, which offers the opportunity of using an SED- and $z$-dependent atmospheric affect on the observed photometry as a color-independent redshift indicator. We show that measuring this airmass effect and including it as a prior has the potential to improve the photometric redshifts and can ameliorate extreme outliers, but that it will only be adequately measured for the brightest galaxies, which limits its overall impact on LSST photometric redshifts. We furthermore demonstrate how this airmass effect can induce a bias in the photo-$z$ results, and caution against survey strategies that prioritize high-airmass observations for the purpose of improving this prior. Ultimately, we intend for this work to serve as a guide for the expectations and preparations of the LSST science community with regards to the minimum quality of photo-$z$ as the survey progresses.

We present an improved photometric error analysis for the 7,100 CRTS (Catalina Real-Time Transient Survey) optical light curves for quasars from the SDSS (Sloan Digital Sky Survey) Stripe 82 catalogue. The SDSS imaging survey has provided a time-resolved photometric data set which greatly improved our understanding of the quasar optical continuum variability: Data for monthly and longer time-scales are consistent with a damped random walk (DRW). Recently, newer data obtained by CRTS provided puzzling evidence for enhanced variability, compared to SDSS results, on monthly time-scales. Quantitatively, SDSS results predict about 0.06 mag root-mean-square (rms) variability for monthly time-scales, while CRTS data show about a factor of 2 larger rms, for spectroscopically confirmed SDSS quasars. Our analysis has successfully resolved this discrepancy as due to slightly underestimated photometric uncertainties from the CRTS image processing pipelines. As a result, the correction for observational noise is too small and the implied quasar variability is too large. The CRTS photometric error correction factors, derived from detailed analysis of non-variable SDSS standard stars that were re-observed by CRTS, are about 20-30%, and result in reconciling quasar variability behaviour implied by the CRTS data with earlier SDSS results. An additional analysis based on independent light curve data for the same objects obtained by the Palomar Transient Factory provides further support for this conclusion. In summary, the quasar variability constraints on weekly and monthly time-scales from SDSS, CRTS and PTF surveys are mutually compatible, as well as consistent with DRW model.

Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required.

The Zwicky Transient Facility is a new robotic-observing program, in which a newly engineered 600-MP digital camera with a pioneeringly large field of view, 47~square degrees, will be installed into the 48-inch Samuel Oschin Telescope at the Palomar Observatory. The camera will generate $\sim 1$~petabyte of raw image data over three years of operations. In parallel related work, new hardware and software systems are being developed to process these data in real time and build a long-term archive for the processed products. The first public release of archived products is planned for early 2019, which will include processed images and astronomical-source catalogs of the northern sky in the $g$ and $r$ bands. Source catalogs based on two different methods will be generated for the archive: aperture photometry and point-spread-function fitting.

SN Hun248 was classified as a nonterminal eruption (a SN "impostor") from a directly identified and highly variable cool hypergiant star. The 2014 outburst achieved peak luminosity equivalent to that of the historic eruption of luminous blue variable (LBV) Eta Car, and exhibited a multipeaked optical light curve that rapidly faded after ~100 days. We report ultraviolet (UV) through optical observations of SN Hunt248 with the Hubble Space Telescope (HST) about 1 yr after the outburst, and mid-infrared observations with the Spitzer Space Telescope before the burst and in decline. The HST data reveal a source that is a factor of ~10 dimmer in apparent brightness than the faintest available measurement of the precursor star. The UV-optical spectral energy distribution (SED) requires a strong Balmer continuum, consistent with a hot B4-B5 photosphere attenuated by grey circumstellar extinction. Substantial mid-infrared excess of the source is consistent with thermal emission from hot dust with a mass of 10^{-6} - 10^{-5} MSun and a geometric extent that is comparable to the expansion radius of the ejecta from the 2014 event. SED modeling indicates that the dust consists of relatively large grains (>0.3 um), which could be related to the grey circumstellar extinction that we infer for the UV-optical counterpart. Revised analysis of the precursor photometry is also consistent with grey extinction by circumstellar dust, and suggests that the initial mass of the star could be twice as large as previously estimated (nearly ~60 MSun). Reanalysis of the earlier outburst data shows that the peak luminosity and outflow velocity of the eruption are consistent with a trend exhibited by stellar merger candidates, prompting speculation that SN Hunt248 may also have stemmed from a massive stellar merger or common-envelope ejection.

Astronomy light curves are sparse, gappy, and heteroscedastic. As a result standard time series methods regularly used for financial and similar datasets are of little help and astronomers are usually left to their own instruments and techniques to classify light curves. A common approach is to derive statistical features from the time series and to use machine learning methods, generally supervised, to separate objects into a few of the standard classes. In this work, we transform the time series to two-dimensional light curve representations in order to classify them using modern deep learning techniques. In particular, we show that convolutional neural networks based classifiers work well for broad characterization and classification. We use labeled datasets of periodic variables from CRTS survey and show how this opens doors for a quick classification of diverse classes with several possible exciting extensions.

We present the results of an extensive observational campaign on the nearby Type Ibn SN 2015G, including data from radio through ultraviolet wavelengths. SN 2015G was asymmetric, showing late-time nebular lines redshifted by ~1000 km/s. It shared many features with the prototypical SN In 2006jc, including extremely strong He I emssion lines and a late-time blue pseudocontinuum. The young SN 2015G showed narrow P-Cygni profiles of He I, but never in its evolution did it show any signature of hydrogen - arguing for a dense, ionized, and hydrogen-free circumstellar medium moving outward with a velocity of ~1000 km/s and created by relatively recent mass loss from the progenitor star. Ultraviolet through infrared observations show that the fading SN 2015G (which was probably discovered some 20 days post-peak) had a spectral energy distribution that was well described by a simple, single-component blackbody. Archival HST images provide upper limits on the luminosity of SN 2015G's progenitor, while nondetections of any luminous radio afterglow and optical nondetections of outbursts over the past two decades provide constraints upon its mass-loss history.

We present observations of two new hydrogen-poor superluminous supernovae (SLSN-I), iPTF15esb and iPTF16bad, showing late-time H-alpha emission with line luminosities of (1-3)e+41 erg/s and velocity widths of (4000-6000) km/s. Including the previously published iPTF13ehe, this makes up a total of three such events to date. iPTF13ehe is one of the most luminous and the slowest evolving SLSNe-I, whereas the other two are less luminous and fast decliners. We interpret this as a result of the ejecta running into a neutral H-shell located at a radius of ~ 1.0e+16cm. This implies that violent mass loss must have occurred several decades before the supernova explosion. Such a short time interval suggests that eruptive mass loss could be common shortly prior to the death of a massive star as a SLSN. And more importantly, helium is unlikely to be completely stripped off the progenitor stars and could be present in the ejecta. It is a mystery why helium features are not detected, even though non-thermal energy sources, capable of ionizing He atoms, may exist as suggested by the O II absorption series in the early time spectra. At late times (+240d), our spectra appear to have intrinsically lower [O I]6300A luminosities than that of SN2015bn and SN2007bi, possibly an indication of smaller oxygen masses (<10-30Msun). The blue-shifted H-alpha emission relative to the hosts for all three events may be in tension with the binary star model proposed for iPTF13ehe. Finally, iPTF15esb has a peculiar light curve with three peaks separated from one another by ~ 22 days. The LC undulation is higher in bluer bands. One possible explanation is eject-CSM interaction.

We present late-time spectra of eight Type Ia supernovae (SNe Ia) obtained at $>200$ days after peak brightness using the Gemini South and Keck telescopes. All of the SNe Ia in our sample were nearby, well separated from their host galaxy's light, and have early-time photometry and spectroscopy from the Las Cumbres Observatory (LCO). Parameters are derived from the light curves and spectra such as peak brightness, decline rate, photospheric velocity, and the widths and velocities of the forbidden nebular emission lines. We discuss the physical interpretations of these parameters for the individual SNe Ia and the sample in general, including comparisons to well-observed SNe Ia from the literature. There are possible correlations between early-time and late-time spectral features that may indicate an asymmetric explosion, so we discuss our sample of SNe within the context of models for an offset ignition and/or white dwarf collisions. A subset of our late-time spectra are uncontaminated by host emission, and we statistically evaluate our nondetections of H$\alpha$ emission to limit the amount of hydrogen in these systems. Finally, we consider the late-time evolution of the iron emission lines, finding that not all of our SNe follow the established trend of a redward migration at $>200$ days after maximum brightness.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

The jet compositions, central engines, and progenitors of gamma-ray bursts (GRBs) remain open questions in GRB physics. Applying broadband observations, including GRB prompt emission and afterglow properties derived from {\em Fermi} and {\em Swift} data, as well as from Keck host-galaxy observations, we address these questions for the peculiar, bright GRB 110731A. By using the pair-opacity method, we derive $\Gamma_{0}>190$ during the prompt emission phase. Alternatively, we derive $\Gamma_{0} \approx 580$ and $\Gamma_{0} \approx 154$ by invoking the early-afterglow phase within the homogeneous density and wind cases, respectively. On the other hand, nondetection of a thermal component in the spectra suggests that the prompt emission is likely powered by dissipation of a Poynting-flux-dominated jet leading to synchrotron radiation in an optically thin region. The nondetection of a jet break in the X-ray and optical bands allows us to place a lower limit on the jet opening angle $\theta_{\rm j}>5.5^{\circ}$. Within a millisecond magnetar central engine scenario, we derive the period $P_{0}$ and polar magnetic field strength $B_{\rm p}$, which have extreme (but still allowed) values. The moderately short observed duration (7.3\,s) and relatively large redshift ($z=2.83$) places the burst as a "rest-frame short" GRB, so the progenitor of the burst is subject to debate. Its relatively large $f_{{\rm eff}, z}$ parameter (ratio of the 1\,s peak flux of a pseudo-GRB and the background flux) and a large physical offset from a potential host galaxy suggest that the progenitor of GRB 110731A may be a compact-star merger.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

A systematic study on the long-term periodicities of known Galactic cataclysmic variables (CVs) was conducted. Among 1580 known CVs, 344 sources were matched and extracted from the Palomar Transient Factory (PTF) data repository. The PTF light curves were combined with the Catalina Real-Time Transient Survey (CRTS) light curves and analyzed. Ten targets were found to exhibit long-term periodic variability, which is not frequently observed in the CV systems. These long-term variations are possibly caused by various mechanisms, such as the precession of the accretion disk, hierarchical triple star system, magnetic field change of the companion star, and other possible mechanisms. We discuss the possible mechanisms in this study. If the long-term period is less than several tens of days, the disk precession period scenario is favored. However, the hierarchical triple star system or the variations in magnetic field strengths are most likely the predominant mechanisms for longer periods.

There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being atop of the normal, stochastic variability of quasars. We have identified 51 events from over 900,000 known quasars and high probability quasar candidates, typically lasting 900 days and with a median peak amplitude of $\Delta m = 1.25$ mag. Characterizing the flare profile with a Weibull distribution, we find that nine of the sources are well described by a single-point single-lens model. This supports the proposal by Lawrence et al. (2016) that microlensing is a plausible physical mechanism for extreme variability. However, we attribute the majority of our events to explosive stellar-related activity in the accretion disk: superluminous supernovae, tidal disruption events, and mergers of stellar mass black holes.

The optical transient PTF11kx exhibited both the characteristic spectral features of Type Ia supernovae (SNe Ia) and the signature of ejecta interacting with circumstellar material (CSM) containing hydrogen, indicating the presence of a nondegenerate companion. We present an optical spectrum at $1342$ days after peak from Keck Observatory, in which the broad component of H$\alpha$ emission persists with a similar profile as in early-time observations. We also present $Spitzer$ IRAC detections obtained $1237$ and $1818$ days after peak, and an upper limit from $HST$ ultraviolet imaging at $2133$ days. We interpret our late-time observations in context with published results - and reinterpret the early-time observations - in order to constrain the CSM's physical parameters and compare to theoretical predictions for recurrent nova systems. We find that the CSM's radial extent may be several times the distance between the star and the CSM's inner edge, and that the CSM column density may be two orders of magnitude lower than previous estimates. We show that the H$\alpha$ luminosity decline is similar to other SNe with CSM interaction, and demonstrate how our infrared photometry is evidence for newly formed, collisionally heated dust. We create a model for PTF11kx's late-time CSM interaction and find that X-ray reprocessing by photoionization and recombination cannot reproduce the observed H$\alpha$ luminosity, suggesting that the X-rays are thermalized and that H$\alpha$ radiates from collisional excitation. Finally, we discuss the implications of our results regarding the progenitor scenario and the geometric properties of the CSM for the PTF11kx system.

The 2012 brightening of SN 2009ip was dominated by emission from the interaction of ejecta with the surrounding circumstellar material (CSM) produced by episodic mass loss from the progenitor, complicating the diagnosis of whether the underlying explosion was a true supernova or a nonterminal eruption of a massive star. In this paper, we contribute a time series of optical photometric and spectroscopic observations for SN 2009ip from 1 to 3 years after the 2012 outburst, collected at the Las Cumbres Observatory and the Keck Observatory. We find that the brightness of SN 2009ip continues to decline with no deviations from a linear slope of $0.0030\pm0.0005$ $\rm mag\ day^{-1}$ in the $r^{\prime}$ band, and demonstrate that this is similar to both observations and models of CSM-ejecta interaction. We show that the late-time spectra continue to be dominated by the signature features of CSM interaction, and that the large ratio of $L_{\rm H\alpha}/L_{\rm H\beta}\approx40$ implies that the material remains optically thick to Balmer photons ("Case C" recombination). We combine our late-time photometry and spectra with early-time data for SN 2009ip and provide a comprehensive discussion that incorporates recently published models and observations for transient phenomena dominated by CSM-ejecta interaction, and conclude that the presence of broad H$\alpha$ at early times remains among the best evidence that a terminal supernova has occurred. Finally, we compare our late-time spectra to those of Type IIn SNe and SN impostors at late phases and find that although SN 2009ip has some similarities with both types, it has more differences with late-time impostor spectra.

CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.

We present the most precise estimate to date of the clustering of quasars on very small scales, based on a sample of 47 binary quasars with magnitudes of $g<20.85$ and proper transverse separations of $\sim 25\,h^{-1}$\,kpc. Our sample of binary quasars, which is about 6 times larger than any previous spectroscopically confirmed sample on these scales, is targeted using a Kernel Density Estimation technique (KDE) applied to Sloan Digital Sky Survey (SDSS) imaging over most of the SDSS area. Our sample is "complete" in that all of the KDE target pairs with $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc in our area of interest have been spectroscopically confirmed from a combination of previous surveys and our own long-slit observational campaign. We catalogue 230 candidate quasar pairs with angular separations of $<8\arcsec$, from which our binary quasars were identified. We determine the projected correlation function of quasars ($\bar W_{\rm p}$) in four bins of proper transverse scale over the range $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc. The implied small-scale quasar clustering amplitude from the projected correlation function, integrated across our entire redshift range, is $A=24.1\pm3.6$ at $\sim 26.6 ~h^{-1}$\,kpc. Our sample is the first spectroscopically confirmed sample of quasar pairs that is sufficiently large to study how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ kpc. We find that empirical descriptions of how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ Mpc also adequately describe the evolution of quasar clustering at $\sim 25 ~h^{-1}$ kpc.

The time domain is the emerging forefront of astronomical research with new facilities and instruments providing unprecedented amounts of data on the temporal behavior of astrophysical populations. Dealing with the size and complexity of this requires new techniques and methodologies. Quasars are an ideal work set for developing and applying these: they vary in a detectable but not easily quantifiable manner whose physical origins are poorly understood. In this paper, we will review how quasars are identified by their variability and how these techniques can be improved, what physical insights into their variability can be gained from studying extreme examples of variability, and what approaches can be taken to increase the number of quasars known. These will demonstrate how astroinformatics is essential to discovering and understanding this important population.

SN 2013ej is a well-studied core-collapse supernova (SN) that stemmed from a directly identified red supergiant (RSG) progenitor in galaxy M74. The source exhibits signs of substantial geometric asphericity, X-rays from persistent interaction with circumstellar material (CSM), thermal emission from warm dust, and a light curve that appears intermediate between supernovae of Types II-P and II-L. The proximity of this source motivates a close inspection of these physical characteristics and their potential interconnection. We present multi-epoch spectropolarimetry of SN 2013ej during the first 107 days, and deep optical spectroscopy and ultraviolet through infrared photometry past ~800 days. SN 2013ej exhibits the strongest and most persistent continuum and line polarization ever observed for a SN of its class during the recombination phase. Modeling indicates that the data are consistent with an oblate ellipsoidal photosphere, viewed nearly edge-on, and probably augmented by optical scattering from circumstellar dust. We suggest that interaction with an equatorial distribution of CSM, perhaps the result of binary evolution, is responsible for generating the photospheric asphericity. Relatedly, our late-time optical imaging and spectroscopy shows that asymmetric CSM interaction is ongoing, and the morphology of broad H-alpha emission from shock-excited ejecta provides additional evidence that the geometry of the interaction region is ellipsoidal. Alternatively, a prolate ellipsoidal geometry from an intrinsically bipolar explosion is also a plausible interpretation of the data, but would probably require a ballistic jet of radioactive material capable of penetrating the hydrogen envelope early in the recombination phase (abridged).

Herein we analyse late-time (post-plateau; 103 < t < 1229 d) optical spectra of low-redshift (z < 0.016), hydrogen-rich Type IIP supernovae (SNe IIP). Our newly constructed sample contains 91 nebular spectra of 38 SNe IIP, which is the largest dataset of its kind ever analysed in one study, and many of the objects have complementary photometric data. We determined the peak and total luminosity, velocity of the peak, HWHM intensity, and profile shape for many emission lines. Temporal evolution of these values and various flux ratios are studied. We also investigate the correlations between these measurements and photometric observables, such as the peak and plateau absolute magnitudes and the late-time light curve decline rates in various optical bands. The strongest and most robust result we find is that the luminosities of all spectral features (except those of helium) tend to be higher in objects with steeper late-time V-band decline rates. A steep late-time V-band slope likely arises from less efficient trapping of gamma-rays and positrons, which could be caused by multidimensional effects such as clumping of the ejecta or asphericity of the explosion itself. Furthermore, if gamma-rays and positrons can escape more easily, then so can photons via the observed emission lines, leading to more luminous spectral features. It is also shown that SNe IIP with larger progenitor stars have ejecta with a more physically extended oxygen layer that is well-mixed with the hydrogen layer. In addition, we find a subset of objects with evidence for asymmetric Ni-56 ejection, likely bipolar in shape. We also compare our observations to theoretical late-time spectral models of SNe IIP from two separate groups and find moderate-to-good agreement with both sets of models. Our SNe IIP spectra are consistent with models of 12-15 M_Sun progenitor stars having relatively low metallicity (Z $\le$ 0.01).

SN2005ip was a TypeIIn event notable for its sustained strong interaction with circumstellar material (CSM), coronal emission lines, and IR excess, interpreted as shock interaction with the very dense and clumpy wind of an extreme red supergiant. We present a series of late-time spectra of SN2005ip and a first radio detection of this SN, plus late-time X-rays, all of which indicate that its CSM interaction is still strong a decade post-explosion. We also present and discuss new spectra of geriatric SNe with continued CSM interaction: SN1988Z, SN1993J, and SN1998S. From 3-10 yr post-explosion, SN2005ip's H-alpha luminosity and other observed characteristics were nearly identical to those of the radio-luminous SN1988Z, and much more luminous than SNe1993J and 1998S. At 10 yr after explosion, SN2005ip showed a drop in H$\alpha$ luminosity, followed by a quick resurgence over several months. We interpret this variability as ejecta crashing into a dense shell located at around 0.05 pc from the star, which may be the same shell that caused the IR echo at earlier epochs. The extreme H-alpha luminosities in SN2005ip and SN1988Z are still dominated by the forward shock at 10 yr post-explosion, whereas SN1993J and SN1998S are dominated by the reverse shock at a similar age. Continuous strong CSM interaction in SNe~2005ip and 1988Z is indicative of enhanced mass loss for about 1e3 yr before core collapse, longer than Ne, O, or Si burning phases. Instead, the episodic mass loss must extend back through C burning and perhaps even part of He burning.

The Type~Ia supernova (SN~Ia) 2016coj in NGC 4125 (redshift $z=0.004523$) was discovered by the Lick Observatory Supernova Search 4.9 days after the fitted first-light time (FFLT; 11.1 days before $B$-band maximum). Our first detection (pre-discovery) is merely $0.6\pm0.5$ day after the FFLT, making SN 2016coj one of the earliest known detections of a SN Ia. A spectrum was taken only 3.7 hr after discovery (5.0 days after the FFLT) and classified as a normal SN Ia. We performed high-quality photometry, low- and high-resolution spectroscopy, and spectropolarimetry, finding that SN 2016coj is a spectroscopically normal SN Ia, but with a high velocity of \ion{Si}{2} $\lambda$6355 ($\sim 12,600$\,\kms\ around peak brightness). The \ion{Si}{2} $\lambda$6355 velocity evolution can be well fit by a broken-power-law function for up to a month after the FFLT. SN 2016coj has a normal peak luminosity ($M_B \approx -18.9 \pm 0.2$ mag), and it reaches a $B$-band maximum \about16.0~d after the FFLT. We estimate there to be low host-galaxy extinction based on the absence of Na~I~D absorption lines in our low- and high-resolution spectra. The spectropolarimetric data exhibit weak polarization in the continuum, but the \ion{Si}{2} line polarization is quite strong ($\sim 0.9\% \pm 0.1\%$) at peak brightness.

We present an analysis of the joint $XMM$-$Newton$ and $NuSTAR$ observations of the radio-quiet quasar QSO B2202$-$209. Using an optical observation from the Hale Telescope at the Palomar Observatory, we revise the redshift of the source from the previously reported $z=1.77$ to $z=0.532$, and we estimate the mass of the central black hole, $\log (M_{\rm BH}/M_{\odot}) = 9.08 \pm 0.18$. The X-ray spectrum of this source can be well described by a power-law of photon index $\Gamma = 1.82 \pm 0.05$ with $E_{\rm cut} = 152_{-54}^{+103}\,{\rm keV}$, in the rest frame of the source. Assuming a Comptonisation model, we estimate the coronal temperature to be $kT_{\rm e}=42\pm 3 \,{\rm keV}$ and $kT_{\rm e}= 56 \pm 3\,{\rm keV}$ for a spherical and a slab geometry, respectively. The coronal properties are comparable to the ones derived for local AGN, despite a difference of around one order of magnitude in black hole mass and X-ray luminosity ($L_{2-10} = 1.93\times 10^{45}\,{\rm erg\,s^{-1}}$). The quasar is X-ray loud, with an unusually flat observed optical-to-X-ray spectral slope $\alpha_{\rm OX} = 1.00 \pm 0.02$, and has an exceptionally strong optical [O III] line. Assuming that both the X-ray emission and the [O III] line are isotropic, these two extreme properties can be explained by a nearly edge-on disk, leading to a reduction in the observed UV continuum light.

The rate of type Ia supernovae (SNe Ia) in a galaxy depends not only on stellar mass, but also on star formation history. Here we show that two simple observational quantities ($g-r$ or $u-r$ host galaxy color, and $r$-band luminosity), coupled with an assumed delay time distribution (the rate of SNe Ia as a function of time for an instantaneous burst of star formation), are sufficient to accurately determine a galaxy's SN Ia rate, with very little sensitivity to the precise details of the star formation history. Using this result, we compare observed and predicted color distributions of SN Ia hosts for the MENeaCS cluster supernova survey, and for the SDSS Stripe 82 supernova survey. The observations are consistent with a continuous delay time distribution (DTD), without any cutoff. For old progenitor systems the power-law slope for the DTD is found to be $-1.50 ^{+0.19} _{-0.15}$. This result favours the double degenerate scenario for SN Ia, though other interpretations are possible. We find that the late-time slopes of the delay time distribution are different at the 1$\sigma$ level for low and high stretch supernova, which suggest a single degenerate scenario for the latter. However, due to ambiguity in the current models' DTD predictions, single degenerate progenitors can neither be confirmed as causing high stretch supernovae nor ruled out from contributing to the overall sample.

We present ultraviolet through near-infrared photometry and spectroscopy of the host galaxies of all superluminous supernovae (SLSNe) discovered by the Palomar Transient Factory prior to 2013, and derive measurements of their luminosities, star-formation rates, stellar masses, and gas-phase metallicities. We find that Type I (hydrogen-poor) SLSNe are found almost exclusively in low-mass (M < 2x10^9 M_sun) and metal-poor (12+log[O/H] < 8.4) galaxies. We compare the mass and metallicity distributions of our sample to nearby galaxy catalogs in detail and conclude that the rate of SLSNe-I as a fraction of all SNe is heavily suppressed in galaxies with metallicities >0.5 Z_sun. Extremely low metallicities are not required, and indeed provide no further increase in the relative SLSN rate. Several SLSN-I hosts are undergoing vigorous starbursts, but this may simply be a side effect of metallicity dependence: dwarf galaxies tend to have bursty star-formation histories. Type-II (hydrogen-rich) SLSNe are found over the entire range of galaxy masses and metallicities, and their integrated properties do not suggest a strong preference for (or against) low-mass/low-metallicity galaxies. Two hosts exhibit unusual properties: PTF 10uhf is a Type I SLSN in a massive, luminous infrared galaxy at redshift z=0.29, while PTF 10tpz is a Type II SLSN located in the nucleus of an early-type host at z=0.04.

We present the catalogue of the Mg II absorption systems detected at a high significance level using an automated search algorithm in the spectra of quasars from the twelfth data release of the Sloan Digital Sky Survey. A total of 266,433 background quasars were searched for the presence of absorption systems in their spectra. The continuum modelling for the quasar spectra was performed using a mean filter. A pseudo-continuum derived using a median filter was used to trace the emission lines. The absorption system catalogue contains 39,694 Mg II systems detected at a 6.0, 3.0$\sigma$ level respectively for the two lines of the doublet. The catalogue was constrained to an absorption line redshift of 0.35 $\le$ z$_{2796}$ $\le$ 2.3. The rest-frame equivalent width of the $\lambda$2796 line ranges between 0.2 $\le$ W$_r$ $\le$ 6.2 $\AA$. Using Gaussian-noise only simulations we estimate a false positive rate of 7.7 per cent in the catalogue. We measured the number density $\partial N^{2796}/\partial z$ of Mg II absorbers and find evidence for steeper evolution of the systems with W$_r \ge$ 1.2 $\AA$ at low redshifts (z$_{2796}$ $\le$ 1.0), consistent with other earlier studies. A suite of null tests over the redshift range 0.5 $\le$ z$_{2796}$ $\le$ 1.5 was used to study the presence of systematics and selection effects like the dependence of the number density evolution of the absorption systems on the properties of the background quasar spectra. The null tests do not indicate the presence of any selection effects in the absorption catalogue if the quasars with spectral signal-to-noise level less than 5.0 are removed. The resultant catalogue contains 36,981 absorption systems. The Mg II absorption catalogue is publicly available and can be downloaded from the link http://srini.ph.unimelb.edu.au/mgii.php

Extensive photometric and spectroscopic observations are presented for SN 2014cx, a type IIP supernova (SN) exploding in the nearby galaxy NGC 337. The observations are performed in optical and ultraviolet bands, covering from -20 to +400 days from the peak light. The stringent detection limit from prediscovery images suggests that this supernova was actually detected within about 1 day after explosion. Evolution of the very early-time light curve of SN 2014cx is similar to that predicted from a shock breakout and post-shock cooling decline before reaching the optical peak. Our photometric observations show that SN 2014cx has a plateau duration of ~ 100 days, an absolute V-band magnitude of ~ -16.5 mag at t~50 days, and a nickel mass of 0.056+-0.008 Msun. The spectral evolution of SN 2014cx resembles that of normal SNe IIP like SN 1999em and SN 2004et, except that it has a slightly higher expansion velocity (~ 4200 km/s at 50 days). From the cooling curve of photospheric temperature, we derive that the progenitor has a pre-explosion radius of ~ 640 Rsun, consistent with those obtained from SNEC modeling (~ 620 Rsun) and hydrodynamical modeling of the observables (~ 570 Rsun). Moreover, the hydrodynamical simulations yield a total explosion energy of ~ 0.4*10e51 erg, and an ejected mass of ~ 8 Msun. These results indicate that the immediate progenitor of SN 2014cx is likely a red supergiant star with a mass of ~ 10 Msun.

This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

Supermassive black hole binaries (SMBHBs) at sub-parsec separations should be common in galactic nuclei, as a result of frequent galaxy mergers. Hydrodynamical simulations of circumbinary discs predict strong periodic modulation of the mass accretion rate on time-scales comparable to the orbital period of the binary. As a result, SMBHBs may be recognized by the periodic modulation of their brightness. We conducted a statistical search for periodic variability in a sample of 35,383 spectroscopically confirmed quasars in the photometric database of the Palomar Transient Factory (PTF). We analysed Lomb-Scargle periodograms and assessed the significance of our findings by modeling each individual quasar's variability as a damped random walk (DRW). We identified 50 quasars with significant periodicity beyond the DRW model, typically with short periods of a few hundred days. We find 33 of these to remain significant after a re-analysis of their periodograms including additional optical data from the intermediate-PTF and the Catalina Real-Time Transient Survey (CRTS). Assuming that the observed periods correspond to the redshifted orbital periods of SMBHBs, we conclude that our findings are consistent with a population of unequal-mass SMBHBs, with a typical mass ratio as low as q = M2/M1 ~ 0.01.

Supernova (SN) 2015U (also known as PSN J07285387+3349106) was discovered in NGC 2388 on 2015 Feb. 11. A rapidly evolving and luminous event, it showed effectively hydrogen-free spectra dominated by relatively narrow helium P-Cygni spectral features and it was classified as a SN Ibn. In this paper we present photometric, spectroscopic, and spectropolarimetric observations of SN 2015U, including a Keck/DEIMOS spectrum (resolution $\approx$ 5000) which fully resolves the optical emission and absorption features. We find that SN 2015U is best understood via models of shock breakout from extended and dense circumstellar material (CSM), likely created by a history of mass loss from the progenitor with an extreme outburst within $\sim$1-2 yr of core collapse (but we do not detect any outburst in our archival imaging of NGC 2388). We argue that the high luminosity of SN 2015U was powered not through $^{56}$Ni decay but via the deposition of kinetic energy into the ejecta/CSM shock interface. Though our analysis is hampered by strong host-galaxy dust obscuration (which likely exhibits multiple components), our dataset makes SN 2015U one of the best-studied Type Ibn supernovae and provides a bridge of understanding to other rapidly fading transients, both luminous and relatively faint.

We present the second multi-frequency radio detection of a reverse shock in a $\gamma$-ray burst. By combining our extensive radio observations of the Fermi-LAT GRB 160509A at $z = 1.17$ up to $20$ days after the burst with Swift X-ray observations and ground-based optical and near-infrared data, we show that the afterglow emission comprises distinct reverse shock and forward shock contributions: the reverse shock emission dominates in the radio band at $\lesssim10~$days, while the forward shock emission dominates in the X-ray, optical, and near-infrared bands. Through multi-wavelength modeling, we determine a circumburst density of $n_0\approx10^{-3}~$cm$^{-3}$, supporting our previous suggestion that a low-density circumburst environment is conducive to the production of long-lasting reverse shock radiation in the radio band. We infer the presence of a large excess X-ray absorption column, $N_{\rm H} \approx 1.5\times10^{22}~$cm$^{-2}$, and a high rest-frame optical extinction, $A_{\rm V}\approx3.4~$mag. We identify a jet break in the X-ray light curve at $t_{\rm jet}\approx6~$d, and thus derive a jet opening angle of $\theta_{\rm jet}\approx4~$deg, yielding a beaming-corrected kinetic energy and radiated $\gamma$-ray energy of $E_{\rm K}\approx4\times10^{50}~$erg and $E_{\gamma}\approx1.3\times10^{51}~$erg ($1$-$10^4~$keV, rest frame), respectively. Consistency arguments connecting the forward and reverse shocks suggest a deceleration time of $t_{\rm dec} \approx 460~$s$~\approx T_{90}$, a Lorentz factor of $\Gamma(t_{\rm dec})\approx330$, and a reverse shock to forward shock fractional magnetic energy density ratio of $R_{\rm B}\equiv\epsilon_{\rm B,RS}/\epsilon_{\rm B,FS}\approx8$.

We study the compatibility of large quasar groups with the concordance cosmological model. Large Quasar Groups are very large spatial associations of quasars in the cosmic web, with sizes of 50-250h^-1 Mpc. In particular, the largest large quasar group known, named Huge-LQG, has a longest axis of ~860h^-1 Mpc, larger than the scale of homogeneity (~260 Mpc), which has been noted as a possible violation of the cosmological principle. Using mock catalogues constructed from the Horizon Run 2 cosmological simulation, we found that large quasar groups size, quasar member number and mean overdensity distributions in the mocks agree with observations. The Huge-LQG is found to be a rare group with a probability of 0.3 per cent of finding a group as large or larger than the observed, but an extreme value analysis shows that it is an expected maximum in the sample volume with a probability of 19 per cent of observing a largest quasar group as large or larger than Huge-LQG. The Huge-LQG is expected to be the largest structure in a volume at least 5.3 +- 1 times larger than the one currently studied.

High-quality collections of Type II supernova (SN) light curves are scarce because they evolve for hundreds of days, making follow-up observations time consuming and often extending over multiple observing seasons. In light of these difficulties, the diversity of SNe II is not fully understood. Here we present ultraviolet and optical photometry of 12 SNe II monitored by the Las Cumbres Observatory Global Telescope Network (LCOGT) during 2013-2014, and compare them with previously studied SNe having well-sampled light curves. We explore SN II diversity by searching for correlations between the slope of the linear light-curve decay after maximum light (historically used to divide SNe II into IIL and IIP) and other measured physical properties. While SNe IIL are found to be on average more luminous than SNe IIP, SNe IIL do not appear to synthesize more 56Ni than SNe IIP. Finally, optical nebular spectra obtained for several SNe in our sample are found to be consistent with models of red supergiant progenitors in the 12-16 Msun range. Consequently, SNe IIL appear not to account for the deficit of massive red supergiants as SN II progenitors.

We investigate the strength of ultraviolet Fe II emission in fainter quasars compared with brighter quasars for 1.0 <= z <= 1.8, using the SDSS (Sloan Digital Sky Survey) DR7QSO catalogue and spectra of Schneider et al., and the SFQS (SDSS Faint Quasar Survey) catalogue and spectra of Jiang et al. We quantify the strength of the UV Fe II emission using the W2400 equivalent width of Weymann et al., which is defined between two rest-frame continuum windows at 2240-2255 and 2665-2695 Ang. The main results are the following. (1) We find that for W2400 >~ 25 Ang. there is a universal (i.e. for quasars in general) strengthening of W2400 with decreasing intrinsic luminosity, L3000. (2) In conjunction with previous work by Clowes et al., we find that there is a further, differential, strengthening of W2400 with decreasing L3000 for those quasars that are members of Large Quasar Groups (LQGs). (3) We find that increasingly strong W2400 tends to be associated with decreasing FWHM of the neighbouring Mg II {\lambda}2798 broad emission line. (4) We suggest that the dependence of W2400 on L3000 arises from Ly{\alpha} fluorescence. (5) We find that stronger W2400 tends to be associated with smaller virial estimates from Shen et al. of the mass of the central black hole, by a factor ~ 2 between the ultrastrong emitters and the weak. Stronger W2400 emission would correspond to smaller black holes that are still growing. The differential effect for LQG members might then arise from preferentially younger quasars in the LQG environments.

We present extensive optical ($UBVRI$, $g'r'i'z'$, and open CCD) and near-infrared ($ZYJH$) photometry for the very nearby Type IIP SN ~2013ej extending from +1 to +461 days after shock breakout, estimated to be MJD $56496.9\pm0.3$. Substantial time series ultraviolet and optical spectroscopy obtained from +8 to +135 days are also presented. Considering well-observed SNe IIP from the literature, we derive $UBVRIJHK$ bolometric calibrations from $UBVRI$ and unfiltered measurements that potentially reach 2\% precision with a $B-V$ color-dependent correction. We observe moderately strong Si II $\lambda6355$ as early as +8 days. The photospheric velocity ($v_{\rm ph}$) is determined by modeling the spectra in the vicinity of Fe II $\lambda5169$ whenever observed, and interpolating at photometric epochs based on a semianalytic method. This gives $v_{\rm ph} = 4500\pm500$ km s$^{-1}$ at +50 days. We also observe spectral homogeneity of ultraviolet spectra at +10--12 days for SNe IIP, while variations are evident a week after explosion. Using the expanding photosphere method, from combined analysis of SN 2013ej and SN 2002ap, we estimate the distance to the host galaxy to be $9.0_{-0.6}^{+0.4}$ Mpc, consistent with distance estimates from other methods. Photometric and spectroscopic analysis during the plateau phase, which we estimated to be $94\pm7$ days long, yields an explosion energy of $0.9\pm0.3\times10^{51}$ ergs, a final pre-explosion progenitor mass of $15.2\pm4.2$~M$_\odot$ and a radius of $250\pm70$~R$_\odot$. We observe a broken exponential profile beyond +120 days, with a break point at +$183\pm16$ days. Measurements beyond this break time yield a $^{56}$Ni mass of $0.013\pm0.001$~M$_\odot$.

In Hubble Space Telescope (HST) imaging taken on 10 November 2014, four images of supernova (SN) "Refsdal" (redshift z=1.49) appeared in an Einstein-cross--like configuration (images S1-S4) around an early-type galaxy in the cluster MACS J1149.5+2223 (z=0.54). Almost all lens models of the cluster have predicted that the SN should reappear within a year in a second host-galaxy image created by the cluster's potential. In HST observations taken on 11 December 2015, we find a new source at the predicted position of the new image of SN Refsdal approximately 8" from the previous images S1-S4. This marks the first time the appearance of a SN at a particular time and location in the sky was successfully predicted in advance! We use these data and the light curve from the first four observed images of SN Refsdal to place constraints on the relative time delay and magnification of the new image (SX), compared to images S1-S4. This enables us, for the first time, to test "blind" lens model predictions of both magnifications and time delays for a lensed SN. We find that the timing and brightness of the new image are consistent with the blind predictions of a fraction of the models. The reappearance illustrates the discriminatory power of this blind test and its utility to uncover sources of systematic uncertainty. From planned HST photometry, we expect to reach a precision of 1-2% on the time delay between S1-S4 and SX.

The LaSilla/QUEST Variability Survey (LSQ) and the Carnegie Supernova Project (CSP II) are collaborating to discover and obtain photometric light curves for a large sample of low redshift (z < 0.1) Type Ia supernovae. The supernovae are discovered in the LSQ survey using the 1 m ESO Schmidt telescope at the La Silla Observatory with the 10 square degree QUEST camera. The follow-up photometric observations are carried out using the 1 m Swope telescope and the 2.5 m du Pont telescopes at the Las Campanas Observatory. This paper describes the survey, discusses the methods of analyzing the data and presents the light curves for the first 31 Type Ia supernovae obtained in the survey. The SALT 2.4 supernova light curve fitter was used to analyze the photometric data, and the Hubble diagram for this first sample is presented. The measurement errors for these supernovae averaged 4%, and their intrinsic spread was 14%.

We report evidence for excess blue light from the Type Ia supernova SN 2012cg at fifteen and sixteen days before maximum B-band brightness. The emission is consistent with predictions for the impact of the supernova on a non-degenerate binary companion. This is the first evidence for emission from a companion to a SN Ia. Sixteen days before maximum light, the B-V color of SN 2012cg is 0.2 mag bluer than for other normal SN~Ia. At later times, this supernova has a typical SN Ia light curve, with extinction-corrected M_B = -19.62 +/- 0.02 mag and Delta m_{15}(B) = 0.86 +/- 0.02. Our data set is extensive, with photometry in 7 filters from 5 independent sources. Early spectra also show the effects of blue light, and high-velocity features are observed at early times. Near maximum, the spectra are normal with a silicon velocity v_{Si} = -10,500$ km s^{-1}. Comparing the early data with models by Kasen (2010) favors a main-sequence companion of about 6 solar masses. It is possible that many other SN Ia have main-sequence companions that have eluded detection because the emission from the impact is fleeting and faint.

We present extensive observations of the Type Ib/c SN2013ge from -13 to +457 days, including spectra and Swift UV-optical photometry beginning 2-4 days post-explosion. This data set makes SN2013ge one of the best observed normal Type Ib/c SN at early times---when the light curve is particularly sensitive to the progenitor configuration and mixing of radioactive elements---and reveals two distinct light curve components in the UV bands. The first component rises over 4-5 days and is visible for the first week post-explosion. Spectra of the first component have blue continua and show a plethora of high velocity (~15,000 km/s) but narrow (~3500 km/s) features, indicating that the line-forming region is restricted. The explosion parameters estimated for the bulk explosion are standard for Type Ib/c SN, and there is evidence for weak He features at early times. In addition, SN2013ge exploded in a low metallicity environment and we have obtained some of the deepest radio and X-ray limits for a Type Ib/c SN to date, which constrain the progenitor mass-loss rate. We are left with two distinct progenitor scenarios for SN2013ge, depending on our interpretation of the early emission. If the first component is cooling envelope emission, then the progenitor of SN2013ge either possessed a low-mass extended envelope or ejected a portion of its envelope in the final <1 year before core-collapse. Alternatively, if the first component is due to outwardly mixed Ni-56, then our observations are consistent with the asymmetric ejection of a distinct clump of nickel-rich material at high velocities. Current models for the collision of a SN shock with a binary companion cannot reproduce both the timescale and luminosity of the early emission in SN2013ge. Finally, the spectra of the first component of SN2013ge are similar to those of the rapidly-declining SN2002bj.

Type IIP supernovae (SNe IIP) have recently been proposed as metallicity ($Z$) probes. The spectral models of Dessart et al. (2014) showed that the pseudo-equivalent width of Fe II $\lambda$5018 (pEW$_{5018}$) during the plateau phase depends on the primordial $Z$, but there was a paucity of SNe IIP exhibiting pEW$_{5018}$ that were compatible with $Z < 0.4$ ${\rm Z}_{\odot}$. This lack might be due to some physical property of the SN II population or to the fact that those SNe have been discovered in luminous, metal-rich targeted galaxies. Here we use SN II observations from the untargeted (intermediate) Palomar Transient Factory [(i)PTF] survey, aiming to investigate the pEW$_{5018}$ distribution of this SN population and, in particular, to look for the presence of SNe II at lower $Z$. We perform pEW$_{5018}$ measurements on the spectra of a sample of 39 (i)PTF SNe II, selected to have well-constrained explosion epochs and light-curve properties. Based on the comparison with the pEW$_{5018}$ spectral models, we subgrouped our SNe into four $Z$ bins from $Z \approx 0.1$ $Z_{\odot}$ up to $Z \approx 2$ $Z_{\odot}$. We also independently investigated the $Z$ of the hosts by using their absolute magnitudes and colors and, in a few cases, using strong-line diagnostics from spectra. We searched for possible correlations between SN observables, such as their peak magnitudes and the $Z$ inferred from pEW$_{5018}$. We found 11 events with pEW$_{5018}$ that were small enough to indicate $Z \approx 0.1$ $Z_{\odot}$. The trend of pEW$_{5018}$ with $Z$ matches the $Z$ estimates obtained from the host-galaxy photometry, although the significance of the correlation is weak. We also found that SNe with brighter peak magnitudes have smaller pEW$_{5018}$ and occur at lower $Z$.

The Galactic transient V1309 Sco was the result of a merger in a low-mass star system, while V838 Mon was thought to be a similar merger event from a more massive B-type progenitor. In this paper we study an optical/IR transient discovered in the nearby galaxy NGC4490, which appeared similar to these merger events (unobscured progenitor, irregular multi-peaked light curve, increasingly red color, similar optical spectrum, IR excess at late times), but which had a higher peak luminosity and longer duration in outburst. NGC4490-OT has less in common with the class of SN~2008S-like transients. A progenitor detected in pre-eruption HST images, combined with upper limits in the IR, requires a luminous and blue progenitor that has faded in late-time HST images. The same source was detected by Spitzer and ground-based data as a luminous IR transient, indicating a transition to a self-obscured state qualitatively similar to the evolution seen in other stellar mergers and in LBVs. The post-outburst dust-obscured source is too luminous and too warm at late times to be explained with an IR echo, suggesting that the object survived the event. The luminosity of the enshrouded IR source is similar to that of the progenitor. Compared to proposed merger events, the more massive progenitor of NGC4490-OT seems to extend a correlation between stellar mass and peak luminosity, and may suggest that both of these correlate with duration. We show that spectra of NGC4490-OT and V838 Mon also resemble light-echo spectra of eta Car, prompting us to speculate that eta Car may be an extreme extension of this phenomenon.

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

We present a high-precision measurement of the parallax for the 12-day Cepheid SS Canis Majoris, obtained via spatial scanning with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). Spatial scanning enables astrometric measurements with a precision of 20-40 muas, an order of magnitude better than pointed observations. SS CMa is the second Cepheid targeted for parallax measurement with HST, and is the first of a sample of eighteen long-period >~ 10 days) Cepheids selected in order to improve the calibration of their period-luminosity relation and eventually permit a determination of the Hubble constant H_0 to better than 2%. The parallax of SS CMa is found to be 348 +/- 38 muas, corresponding to a distance of 2.9 +/- 0.3 kpc. We also present a refinement of the static geometric distortion of WFC3 obtained using spatial scanning observations of calibration fields, with a typical magnitude <~0.01 pixels on scales of 100 pixels.

We have acquired Hubble Space Telescope (HST) and Very Large Telescope near-infrared spectra and images of supernova (SN) Refsdal after its discovery as an Einstein cross in Fall 2014. The HST light curve of SN Refsdal matches the distinctive, slowly rising light curves of SN 1987A-like supernovae (SNe), and we find strong evidence for a broad H-alpha P-Cygni profile in the HST grism spectrum at the redshift (z = 1.49) of the spiral host galaxy. SNe IIn, powered by circumstellar interaction, could provide a good match to the light curve of SN Refsdal, but the spectrum of a SN IIn would not show broad and strong H-alpha absorption. From the grism spectrum, we measure an H-alpha expansion velocity consistent with those of SN 1987A-like SNe at a similar phase. The luminosity, evolution, and Gaussian profile of the H-alpha emission of the WFC3 and X-shooter spectra, separated by ~2.5 months in the rest frame, provide additional evidence that supports the SN 1987A-like classification. In comparison with other examples of SN 1987A-like SNe, SN Refsdal has a blue B-V color and a high luminosity for the assumed range of potential magnifications. If SN Refsdal can be modeled as a scaled version of SN 1987A, we estimate it would have an ejecta mass of 20+-5 solar masses. The evolution of the light curve at late times will provide additional evidence about the potential existence of any substantial circumstellar material (CSM). Using MOSFIRE and X-shooter spectra, we estimate a subsolar host-galaxy metallicity (8.3+-0.1 dex and <8.4 dex, respectively) near the explosion site.

Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EBL), or the radiation near the supermassive black hole lying at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS 1441+25 (z=0.939) during April 2015, a period of high activity across all wavelengths. This observation of PKS 1441+25 suggests that the emission region is located thousands of Schwarzschild radii away from the black hole. The gamma-ray detection also sets a stringent upper limit on the near-ultraviolet to near-infrared EBL intensity, suggesting that galaxy surveys have resolved most, if not all, of the sources of the EBL at these wavelengths.

We present an optical nebular spectrum of the nearby Type Ia supernova 2011fe, obtained 981 days after explosion. SN 2011fe exhibits little evolution since the +593 day optical spectrum, but there are several curious aspects in this new extremely late-time regime. We suggest that the persistence of the $\sim5800$~\AA\ feature is due to Na I D, and that a new emission feature at $\sim7300$~\AA\ may be [Ca II]. Also, we discuss whether the new emission feature at $\sim6400$~\AA\ might be [Fe I] or the high-velocity hydrogen predicted by Mazzali et al. The nebular feature at 5200~\AA\ exhibits a linear velocity evolution of $\sim350$ $\rm km\ s^{-1}$ per 100 days from at least +220 to +980 days, but the line's shape also changes in this time, suggesting that line blending contributes to the evolution. At $\sim 1000$ days after explosion, flux from the SN has declined to a point where contribution from a luminous secondary could be detected. In this work we make the first observational tests for a post-impact remnant star and constrain its temperature and luminosity to $T \gtrsim 10^4$ $\rm K$ and $L \lesssim 10^4$ $\rm L_{\odot}$. Additionally, we do not see any evidence for narrow H$\alpha$ emission in our spectrum. We conclude that observations continue to strongly exclude many single-degenerate scenarios for SN 2011fe.

The optical light curve of the quasar PG 1302-102 at $z = 0.278$ shows a strong, smooth 5.2 yr periodic signal, detectable over a period of $\sim 20$ yr. Although the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. At this close separation, the nuclear black holes in PG 1302-102 will likely merge within $\sim 10^{5}$ yr due to gravitational wave emission alone. Here we report the rest-frame near-infrared time lags for PG 1302-102. Compiling data from {\it WISE} and {\it Akari}, we confirm that the periodic behavior reported in the optical light curve from Graham et al. (2015) is reproduced at infrared wavelengths, with best-fit observed-frame 3.4 and $4.6 \mu$m time lags of $(2219 \pm 153, 2408 \pm 148)$ days for a near face-on orientation of the torus, or $(4103\pm 153, 4292 \pm 148)$ days for an inclined system with relativistic Doppler boosting in effect. The periodicity in the infrared light curves and the light-travel time of the accretion disk photons to reach the dust glowing regions support that a source within the accretion disk is responsible for the optical variability of PG 1302-102, echoed at the further out dusty regions. The implied distance of this dusty, assumed toroidal region is $\sim$ 1.5 pc for a near face-on geometry, or $\sim$1.1 pc for the relativistic Doppler boosted case.

While SN impostors resemble the Great Eruption of eta Car in the sense that their spectra show narrow H lines and they have typical peak absolute magnitudes of -13 to -14 mag, most extragalactic events observed so far are quite different from eta Car in duration. Their bright phases typically last for 100~d or less, rather than persisting for several years. The transient object UGC2773-OT had a similar peak absolute magnitude to other SN impostors, but with a gradual 5-yr prediscovery rise. In the 6 yr since discovery, it has faded very slowly (0.26 mag/yr). Overall, we suggest that its decade-long eruption is so far the best known analog of eta Car's 19th century eruption. We discuss extensive spectroscopy of the ongoing eruption. The spectra show interesting changes in velocity and line shape that we discuss in detail, including an asymmetric Halpha emission line that we show is consistent with the ejection of a bipolar nebula that could be very much like the Homunculus of eta Car. Moreover, changes in the line width, line profile, blue excess emission resembling that of Type IIn supernovae, and the intensity of Halpha suggest the presence of strong circumstellar interaction in the eruption at late times. This supports the hypothesis that the extended plateau of eta Car's eruption may have been powered by shock interaction as well. One interesting difference compared to eta Car, however, is that UGC2773-OT so far does not exhibit the repeated brief spikes in luminosity that have been associated with binary periastron events.

VOSpace is the IVOA interface to distributed storage. This specification presents the first RESTful version of the interface, which is functionally equivalent to the SOAP-based VOSpace 1.1 specification. Note that all prior VOSpace clients will not work with this new version of the interface.

SN 2013dy is a Type Ia supernova for which we have compiled an extraordinary dataset spanning from 0.1 to ~ 500 days after explosion. We present 10 epochs of ultraviolet (UV) through near-infrared (NIR) spectra with HST/STIS, 47 epochs of optical spectra (15 of them having high resolution), and more than 500 photometric observations in the BVrRiIZYJH bands. SN 2013dy has a broad and slowly declining light curve (delta m(B) = 0.92 mag), shallow Si II 6355 absorption, and a low velocity gradient. We detect strong C II in our earliest spectra, probing unburned progenitor material in the outermost layers of the SN ejecta, but this feature fades within a few days. The UV continuum of SN 2013dy, which is strongly affected by the metal abundance of the progenitor star, suggests that SN 2013dy had a relatively high-metallicity progenitor. Examining one of the largest single set of high-resolution spectra for a SN Ia, we find no evidence of variable absorption from circumstellar material. Combining our UV spectra, NIR photometry, and high-cadence optical photometry, we construct a bolometric light curve, showing that SN 2013dy had a maximum luminosity of 10.0^{+4.8}_{-3.8} * 10^{42} erg/s. We compare the synthetic light curves and spectra of several models to SN 2013dy, finding that SN 2013dy is in good agreement with a solar-metallicity W7 model.

Hierarchical assembly models predict a population of supermassive black hole (SMBH) binaries. These are not resolvable by direct imaging but may be detectable via periodic variability (or nanohertz frequency gravitational waves). Following our detection of a 5.2 year periodic signal in the quasar PG 1302-102 (Graham et al. 2015), we present a novel analysis of the optical variability of 243,500 known spectroscopically confirmed quasars using data from the Catalina Real-time Transient Survey (CRTS) to look for close (< 0.1 pc) SMBH systems. Looking for a strong Keplerian periodic signal with at least 1.5 cycles over a baseline of nine years, we find a sample of 111 candidate objects. This is in conservative agreement with theoretical predictions from models of binary SMBH populations. Simulated data sets, assuming stochastic variability, also produce no equivalent candidates implying a low likelihood of spurious detections. The periodicity seen is likely attributable to either jet precession, warped accretion disks or periodic accretion associated with a close SMBH binary system. We also consider how other SMBH binary candidates in the literature appear in CRTS data and show that none of these are equivalent to the identified objects. Finally, the distribution of objects found is consistent with that expected from a gravitational wave-driven population. This implies that circumbinary gas is present at small orbital radii and is being perturbed by the black holes. None of the sources is expected to merge within at least the next century. This study opens a new unique window to study a population of close SMBH binaries that must exist according to our current understanding of galaxy and SMBH evolution.

The light curves of Type Ia supernovae (SNe Ia) are powered by the radioactive decay of $^{56}$Ni to $^{56}$Co at early times, and the decay of $^{56}$Co to $^{56}$Fe from ~60 days after explosion. We examine the evolution of the [Co III] 5892 A emission complex during the nebular phase for SNe Ia with multiple nebular spectra and show that the line flux follows the square of the mass of $^{56}$Co as a function of time. This result indicates both efficient local energy deposition from positrons produced in $^{56}$Co decay, and long-term stability of the ionization state of the nebula. We compile 77 nebular spectra of 25 SN Ia from the literature and present 17 new nebular spectra of 7 SNe Ia, including SN2014J. From these we measure the flux in the [Co III] 5892 A line and remove its well-behaved time dependence to infer the initial mass of $^{56}$Ni ($M_{Ni}$) produced in the explosion. We then examine $^{56}$Ni yields for different SN Ia ejected masses ($M_{ej}$ - calculated using the relation between light curve width and ejected mass) and find the $^{56}$Ni masses of SNe Ia fall into two regimes: for narrow light curves (low stretch s~0.7-0.9), $M_{Ni}$ is clustered near $M_{Ni}$ ~ 0.4$M_\odot$ and shows a shallow increase as $M_{ej}$ increases from ~1-1.4$M_\odot$; at high stretch, $M_{ej}$ clusters at the Chandrasekhar mass (1.4$M_\odot$) while $M_{Ni}$ spans a broad range from 0.6-1.2$M_\odot$. This could constitute evidence for two distinct SN Ia explosion mechanisms.

We present Keck spectroscopic measurements of the millisecond pulsar binary J2215+5135. These data indicate a neutron-star (NS) mass M_NS=1.6Mo, much less than previously estimated. The pulsar heats the companion face to T_D~9000K; the large heating efficiency may be mediated by the intrabinary shock dominating the X-ray light curve. At the best-fit inclination i=88.8deg, the pulsar should be eclipsed. We find weak evidence for such eclipses in the pulsed gamma-rays; an improved radio ephemeris allows use of up to 5 times more Fermi-LAT gamma-ray photons for a definitive test of this picture. If confirmed, the gamma-ray eclipse provides a novel probe of the dense companion wind and the pulsar magnetosphere.

We present deep Hubble Space Telescope imaging at the locations of four, potentially hostless, long-faded Type Ia supernovae (SNe Ia) in low-redshift, rich galaxy clusters that were identified in the Multi-Epoch Nearby Cluster Survey. Assuming a steep faint-end slope for the galaxy cluster luminosity function ($\alpha_d=-1.5$), our data includes all but $\lesssim0.2\%$ percent of the stellar mass in cluster galaxies ($\lesssim0.005\%$ with $\alpha_d=-1.0$), a factor of 10 better than our ground-based imaging. Two of the four SNe Ia still have no possible host galaxy associated with them ($M_R>-9.2$), confirming that their progenitors belong to the intracluster stellar population. The third SNe Ia appears near a faint disk galaxy ($M_V=-12.2$) which has a relatively high probability of being a chance alignment. A faint, red, point source coincident with the fourth SN Ia's explosion position ($M_V=-8.4$) may be either a globular cluster (GC) or faint dwarf galaxy. We estimate the local surface densities of GCs and dwarfs to show that a GC is more likely, due to the proximity of an elliptical galaxy, but neither can be ruled out. This faint host implies that the SN Ia rate in dwarfs or GCs may be enhanced, but remains within previous observational constraints. We demonstrate that our results do not preclude the use of SNe Ia as bright tracers of intracluster light at higher redshifts, but that it will be necessary to first refine the constraints on their rate in dwarfs and GCs with deep imaging for a larger sample of low-redshift, apparently hostless SNe Ia.

We present densely-sampled ultraviolet/optical photometric and low-resolution optical spectroscopic observations of the type IIP supernova 2013ab in the nearby ($\sim$24 Mpc) galaxy NGC 5669, from 2 to 190d after explosion. Continuous photometric observations, with the cadence of typically a day to one week, were acquired with the 1-2m class telescopes in the LCOGT network, ARIES telescopes in India and various other telescopes around the globe. The light curve and spectra suggest that the SN is a normal type IIP event with a plateau duration of $ \sim80 $ days with mid plateau absolute visual magnitude of -16.7, although with a steeper decline during the plateau (0.92 mag 100 d$ ^{-1} $ in $ V $ band) relative to other archetypal SNe of similar brightness. The velocity profile of SN 2013ab shows striking resemblance with those of SNe 1999em and 2012aw. Following the Rabinak & Waxman (2011) prescription, the initial temperature evolution of the SN emission allows us to estimate the progenitor radius to be $ \sim $ 800 R$_{\odot}$, indicating that the SN originated from a red supergiant star. The distance to the SN host galaxy is estimated to be 24.3 Mpc from expanding photosphere method (EPM). From our observations, we estimate that 0.064 M$_{\odot}$ of $^{56}$Ni was synthesized in the explosion. General relativistic, radiation hydrodynamical modeling of the SN infers an explosion energy of $ 0.35\times10^{51} $ erg, a progenitor mass (at the time of explosion) of $ \sim9 $ M$_{\odot}$ and an initial radius of $ \sim600 $ R$_{\odot}$.

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

We present optical observations of the peculiar stripped-envelope supernovae (SNe) LSQ12btw and LSQ13ccw discovered by the La Silla-QUEST survey. LSQ12btw reaches an absolute peak magnitude of M(g) = -19.3 +- 0.2, and shows an asymmetric light curve. Stringent prediscovery limits constrain its rise time to maximum light to less than 4 days, with a slower post-peak luminosity decline, similar to that experienced by the prototypical SN~Ibn 2006jc. LSQ13ccw is somewhat different: while it also exhibits a very fast rise to maximum, it reaches a fainter absolute peak magnitude (M(g) = -18.4 +- 0.2), and experiences an extremely rapid post-peak decline similar to that observed in the peculiar SN~Ib 2002bj. A stringent prediscovery limit and an early marginal detection of LSQ13ccw allow us to determine the explosion time with an uncertainty of 1 day. The spectra of LSQ12btw show the typical narrow He~I emission lines characterising Type Ibn SNe, suggesting that the SN ejecta are interacting with He-rich circumstellar material. The He I lines in the spectra of LSQ13ccw exhibit weak narrow emissions superposed on broad components. An unresolved Halpha line is also detected, suggesting a tentative Type Ibn/IIn classification. As for other SNe~Ibn, we argue that LSQ12btw and LSQ13ccw likely result from the explosions of Wolf-Rayet stars that experienced instability phases prior to core collapse. We inspect the host galaxies of SNe Ibn, and we show that all of them but one are hosted in spiral galaxies, likely in environments spanning a wide metallicity range.

We present optical observations of the peculiar Type Ibn supernova (SN Ibn) OGLE-2012-SN-006, discovered and monitored by the OGLE-IV survey, and spectroscopically followed by PESSTO at late phases. Stringent pre-discovery limits constrain the explosion epoch with fair precision to JD = 2456203.8 +- 4.0. The rise time to the I-band light curve maximum is about two weeks. The object reaches the peak absolute magnitude M(I) = -19.65 +- 0.19 on JD = 2456218.1 +- 1.8. After maximum, the light curve declines for about 25 days with a rate of 4 mag per 100d. The symmetric I-band peak resembles that of canonical Type Ib/c supernovae (SNe), whereas SNe Ibn usually exhibit asymmetric and narrower early-time light curves. Since 25 days past maximum, the light curve flattens with a decline rate slower than that of the 56Co to 56Fe decay, although at very late phases it steepens to approach that rate. An early-time spectrum is dominated by a blue continuum, with only a marginal evidence for the presence of He I lines marking this SN Type. This spectrum shows broad absorptions bluewards than 5000A, likely O II lines, which are similar to spectral features observed in super-luminous SNe at early epochs. The object has been spectroscopically monitored by PESSTO from 90 to 180 days after peak, and these spectra show the typical features observed in a number of SN 2006jc-like events, including a blue spectral energy distribution and prominent and narrow (v(FWHM) ~ 1900 km/s) He I emission lines. This suggests that the ejecta are interacting with He-rich circumstellar material. The detection of broad (10000 km/s) O I and Ca II features likely produced in the SN ejecta (including the [O I] 6300A,6364A doublet in the latest spectra) lends support to the interpretation of OGLE-2012-SN-006 as a core-collapse event.

PTF11iqb was initially classified as a TypeIIn event caught very early after explosion. It showed narrow Wolf-Rayet (WR) spectral features on day 2, but the narrow emission weakened quickly and the spectrum morphed to resemble those of Types II-L and II-P. At late times, Halpha emission exhibited a complex, multipeaked profile reminiscent of SN1998S. In terms of spectroscopic evolution, we find that PTF11iqb was a near twin of SN~1998S, although with weaker interaction with circumstellar material (CSM) at early times, and stronger CSM interaction at late times. We interpret the spectral changes as caused by early interaction with asymmetric CSM that is quickly (by day 20) enveloped by the expanding SN ejecta photosphere, but then revealed again after the end of the plateau when the photosphere recedes. The light curve can be matched with a simple model for weak CSM interaction added to the light curve of a normal SN~II-P. This plateau requires that the progenitor had an extended H envelope like a red supergiant, consistent with the slow progenitor wind speed indicated by narrow emission. The cool supergiant progenitor is significant because PTF11iqb showed WR features in its early spectrum --- meaning that the presence of such WR features in an early SN spectrum does not necessarily indicate a WR-like progenitor. [abridged] Overall, PTF11iqb bridges SNe~IIn with weaker pre-SN mass loss seen in SNe II-L and II-P, implying a continuum between these types.

Quasars have long been known to be variable sources at all wavelengths. Their optical variability is stochastic, can be due to a variety of physical mechanisms, and is well-described statistically in terms of a damped random walk model. The recent availability of large collections of astronomical time series of flux measurements (light curves) offers new data sets for a systematic exploration of quasar variability. Here we report on the detection of a strong, smooth periodic signal in the optical variability of the quasar PG 1302-102 with a mean observed period of 1,884 $\pm$ 88 days. It was identified in a search for periodic variability in a data set of light curves for 247,000 known, spectroscopically confirmed quasars with a temporal baseline of $\sim9$ years. While the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. Such systems are an expected consequence of galaxy mergers and can provide important constraints on models of galaxy formation and evolution.

We use seven year's worth of observations from the Catalina Sky Survey and the Siding Spring Survey covering most of the northern and southern hemisphere at galactic latitudes higher than 20 degrees to search for serendipitously imaged moving objects in the outer solar system. These slowly moving objects would appear as stationary transients in these fast cadence asteroids surveys, so we develop methods to discover objects in the outer solar system using individual observations spaced by months, rather than spaced by hours, as is typically done. While we independently discover 8 known bright objects in the outer solar system, the faintest having $V=19.8\pm0.1$, no new objects are discovered. We find that the survey is nearly 100% efficient at detecting objects beyond 25 AU for $V\lesssim 19.1$ ($V\lesssim18.6$ in the southern hemisphere) and that the probability that there is one or more remaining outer solar system object of this brightness left to be discovered in the unsurveyed regions of the galactic plane is approximately 32%.

We present photometric and spectroscopic observations of SN Hunt 248, a new supernova (SN) impostor in NGC 5806, which began a multi-stage outburst in 2014 May. The initial '2014a' discovery brightening exhibited an absolute magnitude of M~-12 and the spectral characteristics of a cool dense outflow, including P-Cygni lines of Fe II, H I, Na I, and strong line blanketing from metals. The source rapidly climbed and peaked near M~-15 mag after two additional weeks. During this bright '2014b' phase the spectrum became dominated by Balmer emission and a stronger blue continuum, similar to the SN impostor SN 1997bs. Archival images from the Hubble Space Telescope between 1997 and 2005 reveal a luminous (4e5 Lsun) variable precursor star. Its location on the Hertzsprung-Russell diagram is consistent with a massive (M_init~30 Msun) cool hypergiant having an extremely dense wind and an Eddington ratio just below unity. At the onset of the recent 2014a brightening, however, the object became super-Eddington. The subsequent boost in luminosity during the 2014b phase probably resulted from circumstellar interaction. SN Hunt 248 provides the first case of a cool hypergiant undergoing a giant eruption reminiscent of outbursts from luminous blue variable (LBV) stars. This lends support to the hypothesis that some cool hypergiants, such as Rho Cas, could be LBVs masquerading under a pseudo-photosphere created by their extremely dense winds. Moreover, SN Hunt 248 demonstrates that eruptions stemming from such stars can rival in peak luminosity the giant outbursts of much more massive systems like Eta Car.

Torques acting on galaxies lead to physical alignments, but the resulting ellipticity correlations are difficult to predict. As they constitute a major contaminant for cosmic shear studies, it is important to constrain the intrinsic alignment signal observationally. We measured the alignments of satellite galaxies within 90 massive galaxy clusters in the redshift range 0.05<z<0.55 and quantified their impact on the cosmic shear signal. We combined a sample of 38,104 galaxies with spectroscopic redshifts with high-quality data from the Canada-France-Hawaii Telescope. We used phase-space information to select 14,576 cluster members, 14,250 of which have shape measurements and measured three different types of alignment: the radial alignment of satellite galaxies toward the brightest cluster galaxies (BCGs), the common orientations of satellite galaxies and BCGs, and the radial alignments of satellites with each other. Residual systematic effects are much smaller than the statistical uncertainties. We detect no galaxy alignment of any kind out to at least 3 r200. The signal is consistent with zero for both blue and red galaxies, bright and faint ones, and also for subsamples of clusters based on redshift, dynamical mass, and dynamical state. These conclusions are unchanged if we expand the sample with bright cluster members from the red sequence. We augment our constraints with those from the literature to estimate the importance of the intrinsic alignments of satellites compared to those of central galaxies, for which the alignments are described by the linear alignment model. Comparison of the alignment signals to the expected uncertainties of current surveys such as the Kilo-Degree Survey suggests that the linear alignment model is an adequate treatment of intrinsic alignments, but it is not clear whether this will be the case for larger surveys.

We present a time series of the highest resolution spectra yet published for the nearby Type Ia supernova (SN) 2014J in M82. They were obtained at 11 epochs over 33 days around peak brightness with the Levy Spectrograph (resolution R~110,000) on the 2.4m Automated Planet Finder telescope at Lick Observatory. We identify multiple Na I D and K I absorption features, as well as absorption by Ca I H & K and several of the more common diffuse interstellar bands (DIBs). We see no evolution in any component of Na I D, Ca I, or in the DIBs, but do establish the dissipation/weakening of the two most blueshifted components of K I. We present several potential physical explanations, finding the most plausible to be photoionization of circumstellar material, and discuss the implications of our results with respect to the progenitor scenario of SN 2014J.

We present the results of a deep, wide-area variability survey in the Southern hemisphere, the first of its kind. As part of the Catalina Sky Surveys, the Siding Spring Survey (SSS) has covered $14,800$ square degrees in the declination range of $-75^{\circ}\leq\delta\leq-15^{\circ}$. To mine the enormous SSS dataset efficiently we have developed two algorithms: Automatic Period Selection (APS) and Automatic Fourier Decomposition (AFD), which aim to sharpen the period estimation and produce robust lightcurve models. Armed with the APS and AFD outputs we classify $10,540$ ab-type RR Lyrae (RRab) stars ($\sim$90% of which are new) across the Southern sky. As well as the positional information we supply photometric metallicities, and unreddened distances. For the RRab stars in the halo, a study of the photometric metallicity distribution reveals a nearly Gaussian shape with a mean metallicity of ${\rm [Fe/H]}=-1.4$ dex and a dispersion of $0.3$ dex. A spatial study of the RRab metallicities shows no significant radial gradient in the first $\sim7$ kpc from the Galaxy center. However, further out, a small negative gradient is clearly present. This is complemented by a very obvious correlation of the mean RR Lyrae metallicity with distance above the Galactic plane, $z$. We have also carried out an initial substructure search using the discovered RRab, and present the properties of the candidates with significance greater than $2 \sigma$. Most prominent among these is a southern extension of the Sagittarius dwarf galaxy's stream system, reaching down to declinations $\sim -40\deg$.

Effective data visualization is a key part of the discovery process in the era of big data. It is the bridge between the quantitative content of the data and human intuition, and thus an essential component of the scientific path from data into knowledge and understanding. Visualization is also essential in the data mining process, directing the choice of the applicable algorithms, and in helping to identify and remove bad data from the analysis. However, a high complexity or a high dimensionality of modern data sets represents a critical obstacle. How do we visualize interesting structures and patterns that may exist in hyper-dimensional data spaces? A better understanding of how we can perceive and interact with multi dimensional information poses some deep questions in the field of cognition technology and human computer interaction. To this effect, we are exploring the use of immersive virtual reality platforms for scientific data visualization, both as software and inexpensive commodity hardware. These potentially powerful and innovative tools for multi dimensional data visualization can also provide an easy and natural path to a collaborative data visualization and exploration, where scientists can interact with their data and their colleagues in the same visual space. Immersion provides benefits beyond the traditional desktop visualization tools: it leads to a demonstrably better perception of a datascape geometry, more intuitive data understanding, and a better retention of the perceived relationships in the data.

We report on the results of a search for a Weakly Interacting Massive Particle (WIMP) signal in low-energy data of the Cryogenic Dark Matter Search (CDMS~II) experiment using a maximum likelihood analysis. A background model is constructed using GEANT4 to simulate the surface-event background from $^{210}$Pb decay-chain events, while using independent calibration data to model the gamma background. Fitting this background model to the data results in no statistically significant WIMP component. In addition, we perform fits using an analytic ad hoc background model proposed by Collar and Fields, who claimed to find a large excess of signal-like events in our data. We confirm the strong preference for a signal hypothesis in their analysis under these assumptions, but excesses are observed in both single- and multiple-scatter events, which implies the signal is not caused by WIMPs, but rather reflects the inadequacy of their background model.

A growing subset of Type Ia supernovae (SNe Ia) show evidence for unexpected interaction with a dense circumstellar medium (SNe Ia-CSM). The precise nature of the progenitor, however, remains debated owing to spectral ambiguities arising from a strong contribution from the CSM interaction. Late-time spectra offer potential insight if the post-shock cold, dense shell becomes sufficiently thin and/or the ejecta begin to cross the reverse shock. To date, few high-quality spectra of this kind exist. Here we report on the late-time optical and infrared spectra of the SNe~Ia-CSM 2012ca and 2013dn. These SNe Ia-CSM spectra exhibit low [Fe III]/[Fe II] ratios and strong [Ca II] at late epochs. Such characteristics are reminiscent of the super-Chandrasekhar-mass (SC) candidate SN 2009dc, for which these features suggested a low-ionisation state due to high densities, although the broad Fe features admittedly show similarities to the blue "quasi-continuum" observed in some core-collapse SNe Ibn and IIn. Neither SN 2012ca nor any of the other SNe Ia-CSM show evidence for broad oxygen, carbon, or magnesium in their spectra. Similar to the interacting Type IIn SN 2005ip, a number of high-ionisation lines are identified in SN 2012ca, including [S III], [Ar III], [Ar X], [Fe VIII], [Fe X], and possibly [Fe XI]. The total bolometric energy output does not exceed 10^51 erg, but does require a large kinetic-to-radiative conversion efficiency. All of these observations taken together suggest that SNe Ia-CSM are more consistent with a thermonuclear explosion than a core-collapse event, although detailed radiative transfer models are certainly necessary to confirm these results.

We present photospheric-phase observations of LSQ12gdj, a slowly-declining, UV-bright Type Ia supernova. Classified well before maximum light, LSQ12gdj has extinction-corrected absolute magnitude $M_B = -19.8$, and pre-maximum spectroscopic evolution similar to SN 1991T and the super-Chandrasekhar-mass SN 2007if. We use ultraviolet photometry from Swift, ground-based optical photometry, and corrections from a near-infrared photometric template to construct the bolometric (1600-23800 \AA) light curve out to 45 days past $B$-band maximum light. We estimate that LSQ12gdj produced $0.96 \pm 0.07$ $M_\odot$ of $^{56}$Ni, with an ejected mass near or slightly above the Chandrasekhar mass. As much as 27% of the flux at the earliest observed phases, and 17% at maximum light, is emitted bluewards of 3300 \AA. The absence of excess luminosity at late times, the cutoff of the spectral energy distribution bluewards of 3000 \AA, and the absence of narrow line emission and strong Na I D absorption all argue against a significant contribution from ongoing shock interaction. However, up to 10% of LSQ12gdj's luminosity near maximum light could be produced by the release of trapped radiation, including kinetic energy thermalized during a brief interaction with a compact, hydrogen-poor envelope (radius $< 10^{13}$ cm) shortly after explosion; such an envelope arises generically in double-degenerate merger scenarios.

The nearby Type Ia supernovae (SNe Ia) 2011fe and 2011by had nearly identical photospheric phase optical spectra, light-curve widths, and photometric colours, but at peak brightness SN 2011by reached a fainter absolute magnitude in all optical bands and exhibited lower flux in the near-ultraviolet (NUV). Based on those data, Foley & Kirshner (2013) argue that the progenitors of SNe 2011by and 2011fe were supersolar and subsolar, respectively, and that SN 2011fe generated 1.7 times the amount of 56Ni as SN 2011by. With this work, we extend the comparison of these SNe Ia to 10 days before and 300 days after maximum brightness with new spectra and photometry. We show that the nebular phase spectra of SNe 2011fe and 2011by are almost identical, and do not support a factor of 1.7 difference in 56Ni mass. Instead, we find it plausible that the Tully-Fisher distance for SN 2011by is an underestimate, in which case these SNe Ia may have reached similar peak luminosity, formed similar amounts of 56Ni, and had lower metallicity progenitors than previously estimated. Regardless of the true distance to SN 2011by, we find that the relative progenitor metallicity difference remains well supported by their disparity in NUV flux, which we show to be even stronger at pre-maximum epochs --- although contributions from differences in total ejecta mass, viewing angle, or progenitor density cannot be ruled out. We also demonstrate that, independent of distance modulus, SN 2011by exhibits a late-time luminosity excess that cannot be explained by a light echo, but is more likely to be the result of greater energy trapping by the nucleosynthetic products of SN 2011by.

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

We investigate the properties of 367 ultra-short period binary candidates selected from 31,000 sources recently identified from Catalina Surveys data. Based on light curve morphology, along with WISE, SDSS and GALEX multi-colour photometry, we identify two distinct groups of binaries with periods below the 0.22 day contact binary minimum. In contrast to most recent work, we spectroscopically confirm the existence of M-dwarf+M-dwarf contact binary systems. By measuring the radial velocity variations for five of the shortest-period systems, we find examples of rare cool-white dwarf+M-dwarf binaries. Only a few such systems are currently known. Unlike warmer white dwarf systems, their UV flux and their optical colours and spectra are dominated by the M-dwarf companion. We contrast our discoveries with previous photometrically-selected ultra-short period contact binary candidates, and highlight the ongoing need for confirmation using spectra and associated radial velocity measurements. Overall, our analysis increases the number of ultra-short period contact binary candidates by more than an order of magnitude.

Astronomy is undergoing through a methodological revolution triggered by an unprecedented wealth of complex and accurate data. The new panchromatic, synoptic sky surveys require advanced tools for discovering patterns and trends hidden behind data which are both complex and of high dimensionality. We present DAMEWARE (DAta Mining & Exploration Web Application REsource): a general purpose, web-based, distributed data mining environment developed for the exploration of large datasets, and finely tuned for astronomical applications. By means of graphical user interfaces, it allows the user to perform classification, regression or clustering tasks with machine learning methods. Salient features of DAMEWARE include its capability to work on large datasets with minimal human intervention, and to deal with a wide variety of real problems such as the classification of globular clusters in the galaxy NGC1399, the evaluation of photometric redshifts and, finally, the identification of candidate Active Galactic Nuclei in multiband photometric surveys. In all these applications, DAMEWARE allowed to achieve better results than those attained with more traditional methods. With the aim of providing potential users with all needed information, in this paper we briefly describe the technological background of DAMEWARE, give a short introduction to some relevant aspects of data mining, followed by a summary of some science cases and, finally, we provide a detailed description of a template use case.

We constrain the properties of the progenitor system of the highly reddened Type Ia supernova (SN) 2014J in Messier 82 (M82; d ~ 3.5 Mpc). We determine the SN location using Keck-II K-band adaptive optics images, and we find no evidence for flux from a progenitor system in pre-explosion near-ultraviolet through near-infrared Hubble Space Telescope (HST) images. Our upper limits exclude systems having a bright red giant companion, including symbiotic novae with luminosities comparable to that of RS Ophiuchi. While the flux constraints are also inconsistent with predictions for comparatively cool He-donor systems (T < ~35,000 K), we cannot preclude a system similar to V445 Puppis. The progenitor constraints are robust across a wide range of R_V and A_V values, but significantly greater values than those inferred from the SN light curve and spectrum would yield proportionally brighter luminosity limits. The comparatively faint flux expected from a binary progenitor system consisting of white dwarf stars would not have been detected in the pre-explosion HST imaging. Infrared HST exposures yield more stringent constraints on the luminosities of very cool (T < 3000 K) companion stars than was possible in the case of SN Ia 2011fe.