We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

We present a classification of galaxies in the Pan-STARRS1 (PS1) 3$\pi$ survey based on their recent star formation history and morphology. Specifically, we train and test two Random Forest (RF) classifiers using photometric features (colors and moments) from the PS1 data release 2. Labels for the morphological classification are taken from Huertas-Company+2011, while labels for the star formation fraction (SFF) are from the Blanton+2005 catalog. We find that colors provide more predictive accuracy than photometric moments. We morphologically classify galaxies as either early- or late-type, and our RF model achieves a 78\% classification accuracy. Our second model classifies galaxies as having either a low-to-moderate or high SFF. This model achieves an 89\% classification accuracy. We apply both RF classifiers to the entire PS1 $3\pi$ dataset, allowing us to assign two scores to each PS1 source: $P_\mathrm{HSFF}$, which quantifies the probability of having a high SFF, and $P_\mathrm{spiral}$, which quantifies the probability of having a late-type morphology. Finally, as a proof of concept, we apply our classification framework to supernova (SN) host-galaxies from the Zwicky Transient Factory and the Lick Observatory Supernova Search samples. We show that by selecting on $P_\mathrm{HSFF}$ or $P_\mathrm{spiral}$ it is possible to significantly enhance or suppress the fraction of core-collapse SNe (or thermonuclear SNe) in the sample with respect to random guessing. This result demonstrates how contextual information can aid transient classifications at the time of first detection. In the current era of spectroscopically-starved time-domain astronomy, prompt automated classification is paramount.

For decades cataclysmic variables (CVs) were thought to be one of the few classes of accreting compact objects to not launch jets, and have consequently been used to constrain jet launching models. However, recent theoretical and observational advances indicate that CVs do in fact launch jets. Specifically, it was demonstrated that their accretion-outflow cycle is analogous to that of their higher mass cousins -- the X-ray Binaries (XRBs). Subsequent observations of the CV SS Cygni confirmed this and have consistently shown radio flaring equivalent to that in the XRBs that marks a transient jet. Based on this finding and the emission properties, several studies have concluded that the radio emission is most likely from a transient jet. Observations of other CVs, while not conclusive, are consistent with this interpretation. However, the issue is not yet settled. Later observations have raised a number of questions about this model, as well as about potential alternative radio emission mechanisms. CVs are non-relativistic and many have well-determined distances; these properties would make them ideal candidates with which to address many of our outstanding questions about fundamental jet physics. Here we review the case for jets in CVs, discuss the outstanding questions and issues, and outline the future work necessary to conclusively answer the question of whether CVs launch jets.

We present comprehensive multiwavelength radio to X-ray observations of GRB 181201A spanning from $\approx150$ s to $\approx163$ days after the burst, comprising the first joint ALMA-VLA-GMRT observations of a gamma-ray burst (GRB) afterglow. The radio and mm-band data reveal a distinct signature at $\approx3.9$ days, which we interpret as reverse shock (RS) emission. Our observations present the first time that a single radio-frequency spectral energy distribution can be decomposed directly into RS and forward shock (FS) components. We perform detailed modeling of the full multiwavelength data set, using Markov Chain Monte Carlo sampling to construct the joint posterior density function of the underlying physical parameters describing the RS and FS synchrotron emission. We uncover and account for all degeneracies in the model parameters. The joint RS-FS modeling reveals a weakly magnetized ($\sigma\approx3\times10^{-3}$), mildly relativistic RS, from which we derive an initial bulk Lorentz factor of $\Gamma_0\approx103$ for the GRB jet. Our results support the hypothesis that low-density environments are conducive to the observability of RS emission. We compare our observations to other events with strong RS detections, and find a likely observational bias selecting for longer lasting, non-relativistic reverse shocks. We present and begin to address new challenges in modeling posed by the present generation of comprehensive, multi-frequency data sets.

We report on VLBI observations of the fast and blue optical transient (FBOT), AT 2018cow. At ~62 Mpc, AT 2018cow is the first relatively nearby FBOT. The nature of AT 2018cow is not clear, although various hypotheses from a tidal disruption event to different kinds of supernovae have been suggested. It had a very fast rise time (3.5 d) and an almost featureless blue spectrum although high photospheric velocities (40,000 km s$^{-1}$) were suggested early on. The X-ray luminosity was very high, ~$1.4 \times 10^{43}$ erg s$^{-1}$, larger than those of ordinary SNe, and more consistent with those of SNe associated with gamma-ray bursts. Variable hard X-ray emission hints at a long-lived "central engine." It was also fairly radio luminous, with a peak 8.4-GHz spectral luminosity of ~$4 \times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$, allowing us to make VLBI observations at ages between 22 and 287 d. We do not resolve AT 2018cow. Assuming a circularly symmetric source, our observations constrain the average apparent expansion velocity to be <0.49$c$ by t = 98 d (3$\sigma$ limit). We also constrain the proper motion of AT 2018cow to be <0.51$c$. Since the radio emission generally traces the fastest ejecta, our observations make the presence of a long-lived relativistic jet with a lifetime of more than one month very unlikely.

We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_{\nu}\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_{\nu}\propto \nu^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $\theta_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $\Gamma_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $\theta_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $\epsilon_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $\epsilon_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5\sigma$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (\Gamma\beta)^{-\alpha}$ into the environment, finding that shallow stratification indexes $\alpha\le6$ are disfavored.

We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 110-362$ days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously indicated by radio and X-ray data. Combined with contemporaneous radio and X-ray observations, we find no spectral evolution, with a weighted average spectral index of $\langle \beta \rangle = -0.583 \pm 0.013$, demonstrating that no synchrotron break frequencies evolve between the radio and X-ray bands over these timescales. We find that an extrapolation of the post-peak temporal slope of GW170817 to the luminosities of cosmological short GRBs matches their observed jet break times, suggesting that their explosion properties are similar, and that the primary difference in GW170817 is viewing angle. Additionally, we place a deep limit on the luminosity and mass of an underlying globular cluster of $L \lesssim 6.7 \times 10^{3}\,L_{\odot}$, or $M \lesssim 1.3 \times 10^{4}\,M_{\odot}$, at least 4 standard deviations below the peak of the globular cluster mass function of the host galaxy, NGC4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a globular cluster. As highlighted here, HST (and soon JWST) enables critical observations of the optical emission from neutron star merger jets and outflows.

We present comprehensive observations and analysis of the energetic H-stripped SN 2016coi (a.k.a. ASASSN-16fp), spanning the $\gamma$-ray through optical and radio wavelengths, acquired within the first hours to $\sim$420 days post explosion. Our campaign confirms the identification of He in the SN ejecta, which we interpret to be caused by a larger mixing of Ni into the outer ejecta layers. From the modeling of the broad bolometric light curve we derive a large ejecta mass to kinetic energy ratio ($M_{\rm{ej}}\sim 4-7\,\rm{M_{\odot}}$, $E_{\rm{k}}\sim 7-8\times 10^{51}\,\rm{erg}$). The small [\ion{Ca}{ii}] \lam\lam7291,7324 to [\ion{O}{i}] \lam\lam6300,6364 ratio ($\sim$0.2) observed in our late-time optical spectra is suggestive of a large progenitor core mass at the time of collapse. We find that SN 2016coi is a luminous source of X-rays ($L_{X}>10^{39}\,\rm{erg\,s^{-1}}$ in the first $\sim100$ days post explosion) and radio emission ($L_{8.5\,GHz}\sim7\times 10^{27}\,\rm{erg\,s^{-1}Hz^{-1}}$ at peak). These values are in line with those of relativistic SNe (2009bb, 2012ap). However, for SN 2016coi we infer substantial pre-explosion progenitor mass-loss with rate $\dot M \sim (1-2)\times 10^{-4}\,\rm{M_{\odot}yr^{-1}}$ and a sub-relativistic shock velocity $v_{sh}\sim0.15c$, in stark contrast with relativistic SNe and similar to normal SNe. Finally, we find no evidence for a SN-associated shock breakout $\gamma$-ray pulse with energy $E_{\gamma}>2\times 10^{46}\,\rm{erg}$. While we cannot exclude the presence of a companion in a binary system, taken together, our findings are consistent with a massive single star progenitor that experienced large mass loss in the years leading up to core-collapse, but was unable to achieve complete stripping of its outer layers before explosion.

On 2019 April 25.346 and 26.640 UT the LIGO and Virgo gravitational wave (GW) observatories announced the detection of the first candidate events in Observing Run 3 that contain at least one neutron star. S190425z is a likely binary neutron star (BNS) merger at $d_L = 156 \pm 41$ Mpc, while S190426c is possibly the first NS-BH merger ever detected, at $d_L = 377 \pm 100$ Mpc, although with marginal statistical significance. Here we report our optical follow-up observations for both events using the MMT 6.5-m telescope, as well as our spectroscopic follow-up of candidate counterparts (which turned out to be unrelated) with the 4.1-m SOAR telescope. We compare to publicly reported searches, explore the overall areal coverage and depth, and evaluate those in relation to the optical/NIR kilonova emission from the BNS merger GW170817, to theoretical kilonova models, and to short GRB afterglows. We find that for a GW170817-like kilonova, the partial volume covered spans up to about 40% for S190425z and 60% for S190426c. For an on-axis jet typical of short GRBs, the search effective volume is larger, but such a configuration is expected in at most a few percent of mergers. We further find that wide-field $\gamma$-ray and X-ray limits rule out luminous on-axis SGRBs, for a large fraction of the localization regions, although these searches are not sufficiently deep in the context of the $\gamma$-ray emission from GW170817 or off-axis SGRB afterglows. The results indicate that some optical follow-up searches are sufficiently deep for counterpart identification to about 300 Mpc, but that localizations better than 1000 deg$^2$ are likely essential.

Broad-line type Ic Supernovae (BL-Ic SNe) are characterized by high ejecta velocity ($\gtrsim 10^4$ km s$^{-1}$) and are sometimes associated with the relativistic jets typical of long duration ($\gtrsim 2$ s) Gamma-Ray Bursts (L-GRBs). The reason why a small fraction of BL-Ic SNe harbor relativistic jets is not known. Here we present deep X-ray and radio observations of the BL-Ic SN2014ad extending from $13$ to $930$ days post explosion. SN2014ad was not detected at either frequency and has no observational evidence of a GRB counterpart. The proximity of SN2014ad ($d\sim 26$ Mpc) enables very deep constraints on the progenitor mass-loss rate $\dot{M}$ and on the total energy of the fast ejecta $E$. We consider two synchrotron emission scenarios for a wind-like circumstellar medium (CSM): (i) uncollimated non-relativistic ejecta, and (ii) off-axis relativistic jet. Within the first scenario our observations are consistent with GRB-less BL-Ic SNe characterized by a modest energy budget of their fast ejecta ($E \lesssim 10^{45}$ erg), like SNe 2002ap and 2010ay. For jetted explosions, we cannot rule out a GRB with $E \lesssim 10^{51}$ erg (beam-corrected) with a narrow opening angle ($\theta_j \sim 5^{\circ}$) observed moderately off-axis ($\theta_{\rm obs} \gtrsim 30^{\circ}$) and expanding in a very low CSM density ($\dot{M}$ $\lesssim 10^{-6}$ M$_{\odot}$ yr$^{-1}$). Our study shows that off-axis low-energy jets expanding in a low-density medium cannot be ruled out even in the most nearby BL-Ic SNe with extensive deep observations, and might be a common feature of BL-Ic SNe.

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.

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.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

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.

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources.

We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations.

Despite the prevalence of jets in accreting systems and their impact on the surrounding medium, the fundamental physics of how they are launched and collimated is not fully understood. Radio observations of local compact objects, including accreting stellar mass black holes, neutron stars and white dwarfs, probe their jet emission. Coupled with multi-wavelength observations, this allows us to test the underlying accretion-outflow connection and to establish the relationship between the accretor properties and the jet power, which is necessary to accurately model jets. Compact accretors are nearby, numerous and come in a range of accretor properties, and hence are ideal probes for the underlying jet physics. Despite this there are a number of key outstanding questions regarding accretion-driven outflows in these objects that cannot be answered with current radio observations. The vastly improved sensitivity, polarization capabilities, spatial resolution and high-frequency coverage of the ngVLA will be crucial to answering these, and subsequently determining the fundamental physics behind accretion and jets at all physical scales.

We present the first observations of a Type I superluminous supernova (SLSN) at $\gtrsim 1000$ days after maximum light. We observed SN 2015bn using the Hubble Space Telescope Advanced Camera for Surveys in the F475W, F625W and F775W filters at 721 days and 1068 days. SN 2015bn is clearly detected and resolved from its compact host, allowing reliable photometry. A galaxy template constructed from these data further enables us to isolate the SLSN flux in deep ground-based imaging. We measure a light curve decline rate at $>700$ days of $0.19 \pm 0.03$ mag (100 d)$^{-1}$, much shallower than the earlier evolution, and slower than previous SLSNe (at any phase) or the decay rate of $^{56}$Co. Neither additional radioactive isotopes nor a light echo can consistently account for the slow decline. A spectrum at 1083 days shows the same [O I] and [Ca II] lines as seen at $\sim300-400$ days, with no new features to indicate strong circumstellar interaction. Radio limits with the Very Large Array rule out an extended wind for mass-loss rates $10^{-2.7} \lesssim \dot{M}/v_{10} \lesssim 10^{-1.1}$ M$_\odot$ yr$^{-1}$ (where $v_{10}$ is the wind velocity in units of 10 km s$^{-1}$). The optical light curve is consistent with $L \propto t^{-4}$, which we show is expected for magnetar spin-down with inefficient trapping; furthermore, the evolution matches predictions from earlier magnetar model fits. The opacity to magnetar radiation is constrained at $\sim 0.01$ cm$^2$ g$^{-1}$, consistent with photon-matter pair-production over a broad $\sim$GeV-TeV range. This suggests the magnetar spectral energy distribution, and hence the 'missing energy' leaking from the ejecta, may peak in this range.

To understand the physical processes governing accretion discs we can study active galactic nuclei (AGN), X-ray binary systems (XRBs) and cataclysmic variables (CVs). It has been shown that XRBs and CVs show similar observational properties such as recurrent outbursts and aperiodic variability. The latter has been extensively studied for XRBs, but only recently have direct phenomenological analogies been found between XRBs and CVs, including the discovery of the rms--flux relation and the optical detection of Fourier-dependent time-lags. We present a Fourier analysis of the well-known CV SS Cyg in quiescence based on data collected at the 4.2--m William Herschel Telescope using ULTRACAM. Light curves in SDSS filters $u'$, $g'$ and $r'$ were taken simultaneously with sub-second cadence. The high cadence and sensitivity of the camera allow us to study the broad-band noise component of the source in the time range $\sim 10000 - 0.5$ s ($10^{-4}-2$ Hz). Soft/negative lags with an amplitude ~ 5 s at a time-scale of ~ 250 s were observed, indicating that the emission in the redder bands lags the emission in the bluer bands. This effect could be explained by thermal reprocessing of hard photons in the innermost region of the accretion disc, assuming a high viscosity parameter $\alpha > 0.3$, and high irradiation of the disc. Alternatively, it could be associated with the recombination time-scale on the upper layer of the accretions disc.

We present a sub-arcsec resolution radio imaging survey of a sample of 76 low-luminosity active galactic nuclei (LLAGN) that were previously not detected with the Very Large Array at 15 GHz. Compact, parsec-scale radio emission has been detected above a flux density of 40 $\mu$Jy in 60% (45 of 76) of the LLAGN sample. We detect 20 out of 31 (64%) low-ionization nuclear emission-line region (LINER) nuclei, ten out of 14 (71%) low-luminosity Seyfert galaxies, and 15 out of 31 (48%) transition objects. We use this sample to explore correlations between different emission lines and the radio luminosity. We also populate the X-ray and the optical fundamental plane of black hole activity and further refine its parameters. We obtain a fundamental plane relation of log L$_\textrm{R}$ = $0.48\,(\pm0.04$) log L$_\textrm{X}$ + $0.79\,(\pm0.03$) log $\textrm{M}$ and an optical fundamental plane relation of log L$_\textrm{R}$ = $0.63\,(\pm 0.05)$ log L$_{[\rm O~III]}$ + $0.67\,(\pm 0.03)$ log $\textrm{M}$ after including all the LLAGN detected at high resolution at 15 GHz, and the best-studied hard-state X-ray binaries (luminosities are given in erg s$^{-1}$ while the masses are in units of solar mass). Finally, we find conclusive evidence that the nuclear 15 GHz radio luminosity function (RLF) of all the detected Palomar Sample LLAGN has a turnover at the low-luminosity end, and is best-fitted with a broken power law. The break in the power law occurs at a critical mass accretion rate of 1.2$\times$10$^{-3}$ M$_{\odot}$/yr, which translates to an Eddington ratio of $\rm \dot m_{Edd} \sim 5.1 \times 10^{-5}$, assuming a black hole mass of 10$^9 M_{\odot}$. The local group stands closer to the extrapolation of the higher-luminosity sources, and the classical Seyferts agree with the nuclear RLF of the LLAGN in the local universe.

The energy source powering the extreme optical luminosity of hydrogen-stripped Superluminous Supernovae (SLSNe-I) is not known, but recent studies have highlighted the case for a central engine. Radio and/or X-ray observations are best placed to track the fastest ejecta and probe the presence of outflows from a central engine. We compile all the published radio observations of SLSNe-I to date and present three new observations of two new SLSNe-I. None were detected. Through modeling the radio emission, we constrain the sub-parsec environments and possible outflows in SLSNe-I. In this sample we rule out on-axis collimated relativistic jets of the kind detected in Gamma-Ray Bursts (GRBs). We constrain off-axis jets with opening angles of 5\arcdeg\ (30\arcdeg) to energies of $\rm{E_k<4\times10^{50}\,erg}$ ($\rm{E_k<10^{50}\,erg}$) in environments shaped by progenitors with mass-loss rates of $\dot{M}<10^{-4}\,M_{\odot}\,{\rm yr}^{-1}$ ($\dot{M}<10^{-5}\,M_{\odot}\,{\rm yr}^{-1}$) for all off-axis angles, assuming fiducial values $\epsilon_e=0.1$ and $\epsilon_B=0.01$. The deepest limits rule out emission of the kind seen in faint un-collimated GRBs (with the exception of GRB\,060218), and from relativistic supernovae. Finally, for the closest SLSN-I SN 2017egm we constrained the energy of an uncollimated non-relativistic outflow like those observed in normal SNe to $E_{\rm k}\lesssim10^{48}$ erg.

ThunderKAT is the image-plane transients programme for MeerKAT. The goal as outlined in 2010, and still today, is to find, identify and understand high-energy astrophysical processes via their radio emission (often in concert with observations at other wavelengths). Through a comprehensive and complementary programme of surveying and monitoring Galactic synchrotron transients (across a range of compact accretors and a range of other explosive phenomena) and exploring distinct populations of extragalactic synchrotron transients (microquasars, supernovae and possibly yet unknown transient phenomena) - both from direct surveys and commensal observations - we will revolutionise our understanding of the dynamic and explosive transient radio sky. As well as performing targeted programmes of our own, we have made agreements with the other MeerKAT large survey projects (LSPs) that we will also search their data for transients. This commensal use of the other surveys, which remains one of our key programme goals in 2016, means that the combined MeerKAT LSPs will produce by far the largest GHz-frequency radio transient programme to date.

The luminosity distance measurement of GW170817 derived from GW analysis in Abbott et al. 2017 (here, A17:H0) is highly correlated with the measured inclination of the NS-NS system. To improve the precision of the distance measurement, we attempt to constrain the inclination by modeling the broad-band X-ray-to-radio emission from GW170817, which is dominated by the interaction of the jet with the environment. We update our previous analysis and we consider the radio and X-ray data obtained at $t<40$ days since merger. We find that the afterglow emission from GW170817 is consistent with an off-axis relativistic jet with energy $10^{48}\,\rm{erg}<E_{k}\le 3\times 10^{50} \,\rm{erg}$ propagating into an environment with density $n\sim10^{-2}-10^{-4} \,\rm{cm^{-3}}$, with preference for wider jets (opening angle $\theta_j=15$ deg). For these jets, our modeling indicates an off-axis angle $\theta_{\rm obs}\sim25-50$ deg. We combine our constraints on $\theta_{\rm obs}$ with the joint distance-inclination constraint from LIGO. Using the same $\sim 170$ km/sec peculiar velocity uncertainty assumed in A17:H0 but with an inclination constraint from the afterglow data, we get a value of $H_0=$$74.0 \pm \frac{11.5}{7.5}$ $\mbox{km/s/Mpc}$, which is higher than the value of $H_0=$$70.0 \pm \frac{12.0}{8.0}$ $\mbox{km/s/Mpc}$ found in A17:H0. Further, using a more realistic peculiar velocity uncertainty of 250 km/sec derived from previous work, we find $H_0=$$75.5 \pm \frac{11.6}{9.6}$ km/s/Mpc for H0 from this system. We note that this is in modestly better agreement with the local distance ladder than the Planck CMB, though a significant such discrimination will require $\sim 50$ such events. Future measurements at $t>100$ days of the X-ray and radio emission will lead to tighter constraints.

We present the results from a sensitive X-ray survey of 26 nearby hydrogen-poor superluminous supernovae (SLSNe-I) with Swift, Chandra and XMM. This dataset constrains the SLSN evolution from a few days until ~2000 days after explosion, reaching a luminosity L_x~10^40 erg/s and revealing the presence of significant X-ray emission at the location of PTF12dam. No SLSN-I is detected above L_x~10^41 erg/s, suggesting that the luminous X-ray emission L_x~10^45 erg/s associated with SCP60F6 is not common among SLSNe-I. We constrain the presence of off-axis GRB jets, ionization breakouts from magnetar central engines and the density in the sub-pc environments of SLSNe-I through Inverse Compton emission. The deepest limits rule out the weakest uncollimated GRB outflows, suggesting that IF the similarity of SLSNe-I with GRB/SNe extends to their fastest ejecta, then SLSNe-I are either powered by energetic jets pointed far away from our line of sight theta>30 deg, or harbor failed jets that do not successfully break through the stellar envelope. Furthermore, IF a magnetar central engine is responsible for the exceptional luminosity of SLSNe-I, our X-ray analysis favors large magnetic fields B>2x10^(14) G and ejecta masses M_ej>3 Msun in agreement with optical/UV studies. Finally, we constrain the pre-explosion mass-loss rate of stellar progenitors of SLSNe-I. For PTF12dam we infer Mdot<2x10^(-5) Msun/yr, suggesting that the SN shock interaction with the CSM is unlikely to supply the main source of energy powering the optical transient and that some SLSN-I progenitors end their life as compact stars surrounded by a low-density medium similar to long GRBs and Type Ib/c SNe.

High-redshift quasars are important to study galaxy and active galactic nuclei (AGN) evolution, test cosmological models, and study supermassive black hole growth. Optical searches for high-redshift sources have been very successful, but radio searches are not hampered by dust obscuration and should be more effective at finding sources at even higher redshifts. Identifying high-redshift sources based on radio data is, however, not trivial. Here we report on new multi-frequency Giant Metrewave Radio Telescope (GMRT) observations of eight z>4.5 sources previously studied at high angular resolution with very long baseline interferometry (VLBI). Combining these observations with those from the literature, we construct broad-band radio spectra of all 30 z>4.5 sources that have been observed with VLBI. In the sample we found flat, steep and peaked spectra in approximately equal proportions. Despite several selection effects, we conclude that the z>4.5 VLBI (and likely also non-VLBI) sources have diverse spectra and that only about a quarter of the sources in the sample have flat spectra. Previously, the majority of high-redshift radio sources were identified based on their ultra-steep spectra (USS). Recently a new method has been proposed to identify these objects based on their megahertz-peaked spectra (MPS). Neither method would have identified more than 18% of the high-redshift sources in this sample. More effective methods are necessary to reliably identify complete samples of high-redshift sources based on radio data.

We present the detection of persistent soft X-ray radiation with L_x ~ 10^41-10^42 erg/s at the location of the extremely luminous, double-humped transient ASASSN-15lh as revealed by Chandra and Swift. We interpret this finding in the context of observations from our multiwavelength campaign, which revealed the presence of weak narrow nebular emission features from the host-galaxy nucleus and clear differences with respect to superluminous supernova optical spectra. Significant UV flux variability on short time-scales detected at the time of the re-brightening disfavors the shock interaction scenario as the source of energy powering the long-lived UV emission, while deep radio limits exclude the presence of relativistic jets propagating into a low-density environment. We propose a model where the extreme luminosity and double-peaked temporal structure of ASASSN-15lh is powered by a central source of ionizing radiation that produces a sudden change of the ejecta opacity at later times. As a result, UV radiation can more easily escape, producing the second bump in the light-curve. We discuss different interpretations for the intrinsic nature of the ionizing source. We conclude that, IF the X-ray source is physically associated with the optical-UV transient, ASASSN-15lh most likely represents the tidal disruption of a main-sequence star by the most massive spinning black hole detected to date. In this case, ASASSN-15lh and similar events discovered in the future would constitute the most direct probes of very massive, dormant, spinning, supermassive black holes in galaxies. Future monitoring of the X-rays may allow us to distinguish between the supernova and TDE hypothesis.

We present 8--12\,GHz radio light curves of five dwarf nova (DN) type Cataclysmic Variable stars (CVs) in outburst (RX And, U Gem and Z Cam), or superoutburst (SU UMa and YZ Cnc), increasing the number of radio-detected DN by a factor of two. The observed radio emission was variable on time-scales of minutes to days, and we argue that it is likely to be synchrotron emission. This sample shows no correlation between the radio luminosity and optical luminosity, orbital period, CV class, or outburst type; however higher-cadence observations are necessary to test this, as the measured luminosity is dependent on the timing of the observations in these variable objects. The observations show that the previously detected radio emission from SS Cyg is not unique in type, luminosity (in the plateau phase of the outburst), or variability time-scales. Our results prove that DN, as a class, are radio emitters in outburst.

The Outburst Catalogue contains a wide variety of observational properties for 722 dwarf nova-type (DN) cataclysmic variables (CVs) and 309 CVs of other types from the Catalina Real-time Transient Survey. In particular, it includes the apparent outburst and quiescent V-band magnitudes, duty cycles, limits on the recurrence time, upper- and lower-limits on the distance and absolute quiescent magnitudes, colour information, orbital parameters, and X-ray counterparts. These properties were determined by means of a classification script presented in this paper. The DN in the catalogue show a correlation between the outburst duty cycle and the orbital period (and outburst recurrence time), as well as between the quiescent absolute magnitude and the orbital period (and duty cycle). This is the largest sample of dwarf nova properties collected to date. Besides serving as a useful reference for individual systems and a means of selecting objects for targeted studies, it will prove valuable for statistical studies that aim to shed light on the formation and evolution of cataclysmic variables.

Radio emission from non-magnetic cataclysmic variables (CVs, accreting white dwarfs) could allow detailed studies of outflows and possibly accretion flows in these nearby, numerous and non-relativistic compact accretors. Up to now, however, very few CVs have been detected in the radio. We have conducted a VLA pilot survey of four close and optically-bright novalike CVs at 6 GHz, detecting three, and thereby doubling the number of radio detections of these systems. RW Sex, V603 Aql and the old nova TT Ari were detected in both of the epochs, while V1084 Her was not detected (to a $3\sigma$ upper-limit of 7.8 $\mu\rm{Jy}\,\rm{beam}^{-1}$). These observations clearly show that the sensitivity of previous surveys was typically too low to detect these objects and that non-magnetic CVs can indeed be significant radio emitters. The three detected sources show a range of properties, including flaring and variability on both short ($\sim$200 s) and longer-term (days) time-scales, as well as circular polarization levels of up to 100\%. The spectral indices range from steep to inverted; TT Ari shows a spectral turnover at $\sim$6.5 GHz, while the spectral index of V603 Aql flattened from $\alpha=0.54\pm0.05$ to $0.16\pm0.08$ ($F_{\nu}\propto \nu^{\alpha}$) in the week between observations. This range of properties suggests that more than one emission process can be responsible for the radio emission in non-magnetic CVs. In this sample we find that individual systems are consistent with optically thick synchrotron emission, gyrosynchrotron emission or cyclotron maser emission.

Time series photometry of 20 Cataclysmic Variables detected by the Catalina Real-time Transient Survey is presented. 14 of these systems have not been observed previously and only two have been examined in-depth. From the observations we determined 12 new orbital periods and independently found a further two. Eight of the CVs are eclipsing systems, five of which have eclipse depths of more than 0.9 mag. Included in the sample are six SU UMa systems (three of which show superhumps in our photometry), a polar (SSS1944-42) and one system (CSS1417-18) that displays an abnormally fast decline from outburst.