It is believed that the gas accretion onto the supermassive black holes (SMBHs) is the main process of powering its luminous emission, which occurs in optical, UV and X-ray regimes and less frequently in radio waves. The observational fact that only a few percent of quasars are radio-loud is still an unresolved issue concerning the understanding of the active galactic nucleus (AGN) population. Here we present a detection of a rapid transition from the radio-quiet to the radio-loud mode in quasar 013815+00 (z=0.94) which coincides with changes of its UV-optical continuum and the low ionization MgII broadline. We interpret this as an enhancement of accretion onto a central black hole of mass about 10^9 solar masses. As a consequence a new radio-loud AGN was born. Its spectral and morphological properties indicate that it went through the short gigahertz-peaked spectrum (GPS) phase at the beginning of its activity and has now stabilized its flux density at the level of a few mJy. The radio morphology of 013815+00 is very compact and we predict that with such short-term jet activity its development will be very slow. The observed luminosity changes of the accretion disk are shorter than the lifetime of the new radio phase in 013815+00.

We investigate the imaging performance of an interferometric array in the case of wide field, high resolution, narrow band, snapshot imaging. We find that, when uv-cell sizes are sufficiently small (ie. image sizes are sufficiently large), each instantaneous visibility record is gridded into its own uv-cell. This holds even for dense arrays, like the core of the next generation VLA. In this particular, application, Uniform weighting of the gridded visibilities approaches Natural weighting, with its often deleterious consequences on the resulting synthesized beam. For a core-dominated array, we show that the resulting image noise is highly correlated on scales comparable to the spatial frequencies of the core baselines. In general, this study accentuates the fact that, for imaging applications that require high resolution (Plains array and greater), many of the core antennas can be employed as a separate subarray for low resolution science, without sacrificing the quality of the high resolution science.

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 the results of MeerKAT radio observations of eleven nearby novalike cataclysmic variables. We have detected radio emission from IM Eri, RW Sex, V3885 Sgr and V603 Aql. While RW Sex, V3885 Sgr and V603 Aql had been previously detected, this is the first reported radio detection of IM Eri. Our observations have doubled the sample of non-magnetic CVs with sensitive radio data. We observe that at our radio detection limits, a specific optical luminosity $\gtrsim 2.2\times 10^{18}$ erg/s/Hz (corresponding to $M_V \lesssim 6.0$) is required to produce a radio detection. We also observe that the X-ray and radio luminosities of our detected novalikes are on an extension of the $L_{X} \propto L_{R}^{\sim 0.7}$ power law originally proposed for non-pulsating neutron star low-mass X-ray binaries. We find no other correlations between the radio emission and emission in other wavebands or any other system parameters for the existing sample of radio-detected non-magnetic CVs. We measure in-band (0.9-1.7 GHz) radio spectral indices that are consistent with reports from earlier work. Finally, we constructed broad spectral energy distributions for our sample from published multi-wavelength data, and use them to place constraints on the mass transfer rates of these systems.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

The radio:X-ray correlation that characterises accreting black holes at all mass scales - from stellar mass black holes in binary systems to super-massive black holes powering Active Galactic Nuclei - is one of the most important pieces of observational evidence supporting the existence of a connection between the accretion process and the generation of collimated outflows - or jets - in accreting systems. Although recent studies suggest that the correlation extends down to low luminosities, only a handful of stellar mass black holes have been clearly detected, and in general only upper limits (especially at radio wavelengths) can be obtained during quiescence. We recently obtained detections of the black hole X-ray binary GX 339--4 in quiescence using the MeerKAT radio telescope and Swift X-ray Telescope instrument onboard the Neil Gehrels Swift Observatory, probing the lower end of the radio:X-ray correlation. We present the properties of accretion and of the connected generation of jets in the poorly studied low-accretion rate regime for this canonical black hole XRB system.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We report the discovery of the first transient with MeerKAT, MKT J170456.2$-$482100, discovered in ThunderKAT images of the low mass X-ray binary GX339$-$4. MKT J170456.2$-$482100 is variable in the radio, reaching a maximum flux density of $0.71\pm0.11\,\mathrm{mJy}$ on 2019 Oct 12, and is undetected in 15 out of 48 ThunderKAT epochs. MKT J170456.2$-$482100 is coincident with the chromospherically active K-type sub-giant TYC 8332-2529-1, and $\sim18\,\mathrm{years}$ of archival optical photometry of the star shows that it varies with a period of $21.25\pm0.04\,\mathrm{days}$. The shape and phase of the optical light curve changes over time, and we detect both X-ray and UV emission at the position of MKT J170456.2$-$482100, which may indicate that TYC 8332-2529-1 has large star spots. Spectroscopic analysis shows that TYC 8332-2529-1 is in a binary, and has a line-of-sight radial velocity amplitude of $43\,\mathrm{km\,s^{-1}}$. We also observe a spectral feature in anti-phase with the K-type sub-giant, with a line-of-sight radial velocity amplitude of $\sim12\pm10\,\mathrm{km\,s^{-1}}$, whose origins cannot currently be explained. Further observations and investigation are required to determine the nature of the MKT J170456.2$-$482100 system.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

We present the discovery of a nuclear transient with the Caltech-NRAO Stripe 82 Survey (CNSS), a dedicated radio transient survey carried out with the Karl G. Jansky Very Large Array (VLA). This transient, CNSS J001947.3+003527, exhibited a turn-on over a timescale of $\lesssim$1 yr, increasing in flux density at 3 GHz from $<0.14$ mJy in February 2014 to $4.4\pm0.1$ mJy in March 2015, reaching a peak luminosity of $5\times10^{28} \text{erg s}^{-1} \text{Hz}^{-1}$ around October 2015. The association of CNSS J0019+00 with the nucleus (Gaia and our VLBI positions are consistent to within 1 pc) of a nearby S0 Seyfert galaxy at 77 Mpc, together with the radio spectral evolution, implies that this transient is most likely a tidal disruption event (TDE). Our equipartition analysis indicates the presence of a $\sim$15,000 km s$^{-1}$ outflow, having energy $\sim$10$^{49}$ erg. We derive the radial density profile for the circum-nuclear material in the host galaxy to be proportional to R$^{-2.5}$. All of these properties suggest resemblance with radio-detected thermal TDEs like ASASSN-14li and XMMSL1 J0740-85. No significant X-ray or optical emission is detected from CNSS J0019+00, although this may simply be due to the thermal emission being weak during our late-time follow up observations. From the CNSS survey we have obtained the first unbiased measurement of the rate of radio TDEs, $R(>500 \mu{\rm Jy})$ of about $2\times10^{-3}$ deg$^{-2}$, or equivalently a volumetric rate of about 10 Gpc$^{-3}$ yr$^{-1}$. This rate implies that all-sky radio surveys such as the VLA Sky Survey and those planned with ASKAP, will find many tens of radio TDEs over the next few years.

We have carried out a dedicated transient survey of 300 deg$^2$ of the SDSS Stripe 82 region using the Giant Meterwavelength Radio Telescope (GMRT) at 150 MHz. Our multi-epoch observations, together with the TGSS survey, allow us to probe variability and transient activity on four different timescales, beginning with 4 hours, and up to 4 years. Data calibration, RFI flagging, source finding and transient search were carried out in a semi-automated pipeline incorporating the SPAM recipe. This has enabled us to produce superior-quality images and carry out reliable transient search over the entire survey region in under 48 hours post-observation. Among the few thousand unique point sources found in our 5$\sigma$ single-epoch catalogs (flux density thresholds of about 24 mJy, 20 mJy, 16 mJy and 18 mJy on the respective timescales), we find $<$0.08%, 0.01%, $<$0.06% and 0.05% to be variable (beyond a significance of 4$\sigma$ and fractional variability of 30%) on timescales of 4 hours, 1 day, 1 month and 4 years respectively. This is substantially lower than that in the GHz sky, where $\sim$1% of the persistent point sources are found to be variable. Although our survey was designed to probe a superior part of the transient phase space, our transient sources did not yield any significant candidates. The transient (preferentially extragalactic) rate at 150 MHz is therefore $<$0.005 on timescales of 1 month and 4 years, and $<$0.002 on timescales of 1 day and 4 hours, beyond 7$\sigma$ detection threshold. We put these results in the perspective with the previous studies and give recommendations for future low-frequency transient surveys.

Accreting white dwarfs in binary systems known as cataclysmic variables (CVs) have in recent years been shown to produce radio flares during outbursts, qualitatively similar to those observed from neutron star and black hole X-ray binaries, but their ubiquity and energetic significance for the accretion flow has remained uncertain. We present new radio observations of the CV SS Cyg with Arcminute Microkelvin Imager Large Array, which show for the second time late-ouburst radio flaring, in April 2016. This flaring occurs during the optical flux decay phase, about ten days after the well-established early-time radio flaring. We infer that both the early- and late-outburst flares are a common feature of the radio outbursts of SS Cyg, albeit of variable amplitudes, and probably of all dwarf novae. We furthermore present new analysis of the physical conditions in the best-sampled late-outburst flare, from Feb 2016, which showed clear optical depth evolution. From this we can infer that the synchrotron-emitting plasma was expanding at about 1% of the speed of light, and at peak had a magnetic field of order 1 Gauss and total energy content > 10^{33} erg. While this result is independent of the geometry of the synchrotron-emitting region, the most likely origin is in a jet carrying away a significant amount of the available accretion power.

Powerful relativistic jets are one of the main ways in which accreting black holes provide kinetic feedback to their surroundings. Jets launched from or redirected by the accretion flow that powers them should be affected by the dynamics of the flow, which in accreting stellar-mass black holes has shown increasing evidence for precession due to frame dragging effects that occur when the black hole spin axis is misaligned with the orbital plane of its companion star. Recently, theoretical simulations have suggested that the jets can exert an additional torque on the accretion flow, although the full interplay between the dynamics of the accretion flow and the launching of the jets is not yet understood. Here we report a rapidly changing jet orientation on a timescale of minutes to hours in the black hole X-ray binary V404 Cygni, detected with very long baseline interferometry during the peak of its 2015 outburst. We show that this can be modelled as Lense-Thirring precession of a vertically-extended slim disk that arises from the super-Eddington accretion rate. Our findings suggest that the dynamics of the precessing inner accretion disk could play a role in either directly launching or redirecting the jets within the inner few hundred gravitational radii. Similar dynamics should be expected in any strongly-accreting black hole whose spin is misaligned with the inflowing gas, both affecting the observational characteristics of the jets, and distributing the black hole feedback more uniformly over the surrounding environment.

We present a cadence optimization strategy to unveil a large population of kilonovae using optical imaging alone. These transients are generated during binary neutron star and potentially neutron star-black hole mergers and are electromagnetic counterparts to gravitational-wave signals detectable in nearby events with Advanced LIGO, Advanced Virgo, and other interferometers that will come online in the near future. Discovering a large population of kilonovae will allow us to determine how heavy element production varies with the intrinsic parameters of the merger and across cosmic time. The rate of binary neutron star mergers is still uncertain, but only few (less than 15) events with associated kilonovae may be detectable per year within the horizon of next-generation ground-based interferometers. The rapid evolution (hours to days) at optical/infrared wavelengths, relatively low luminosity, and the low volumetric rate of kilonovae makes their discovery difficult, especially during blind surveys of the sky. We propose future large surveys to adopt a rolling cadence in which g-i observations are taken nightly for blocks of 10 consecutive nights. With the current baseline2018a cadence designed for the Large Synoptic Survey Telescope (LSST), less than 7.5 poorly-sampled kilonovae are expected to be detected in both the Wide Fast Deep (WFD) and Deep Drilling Fields (DDF) surveys per year, under optimistic assumptions on their rate, duration, and luminosity. We estimate the proposed strategy to return up to about 272 GW170817-like kilonovae throughout the LSST WFD survey, discovered independently from gravitational-wave triggers.

Kilonovae and radio afterglows of neutron star merger events have been identified as the two most promising counterparts, of these gravitational wave sources, that can provide arcsecond localization. While several new and existing optical search facilities have been dedicated to finding kilonovae, factors such as dust obscuration and the daytime sky may thwart these searches in a significant fraction of gavitational wave events. Radio-only searches, being almost immune to these factors, are equally capable of finding the counterparts and in fact offer a complementary discovery approach, despite the modest fields of view for many of the present-day radio interferometers. Such interferometers will be able to carry out competitive searches for the electromagnetic counterparts through the galaxy targeting approach. Adapting and improving on an existing algorithm by Rana et al. 2017. we present here a method that optimizes the placement of radio antenna pointings, integration time, and antenna slew. We simulate 3D gravitational wave localizations to find the efficacy of our algorithm; with substantial improvements in slew overhead and containment probability coverage, our algorithm performs significantly better than simple galaxy-rank-ordered observations. We propose that telescopes such as the Very Large Array, MeerKAT, Australia Telescope Compact Array and the Gaint Meterwave Radio Telescope, having fields of view $\lesssim$1 deg$^2$ and searching for the counterparts of nearby GW events over tens of square degrees or larger, will especially benefit from this optimized galaxy-targeting approach for electromagnetic counterpart searches.

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.

In this data paper we present and characterise the multi-component radio sources identified in the VLA-COSMOS Large Project at 3 GHz (0.75 arcsec resolution, 2.3 {\mu}Jy/beam rms), i.e. the radio sources which are composed of two or more radio blobs.The classification of objects into multi-components was done by visual inspection of 351 of the brightest and most extended blobs from a sample of 10,899 blobs identified by the automatic code blobcat. For that purpose we used multi-wavelength information of the field, such as the 1.4-GHz VLA-COSMOS data and the UltraVISTA stacked mosaic available for COSMOS. We have identified 67 multi-component radio sources at 3 GHz: 58 sources with AGN powered radio emission and 9 star-forming galaxies. We report 8 new detections that were not observed by the VLA-COSMOS Large Project at 1.4 GHz, due to the slightly larger area coverage at 3 GHz. The increased spatial resolution of 0.75 arcsec has allowed us to resolve (and isolate) multiple emission peaks of 28 extended radio sources not identified in the 1.4-GHz VLA-COSMOS map. We report the multi-frequency flux densities (324 MHz, 325 MHz, 1.4 GHz & 3 GHz), star-formation-rates, and stellar masses of these objects. Multi-component objects at 3-GHz VLA-COSMOS inhabit mainly massive galaxies (>10^10.5 Msun). The majority of the multi-component AGN lie below the main-sequence of star-forming galaxies (SFGs), in the green valley and the quiescent region. We provide detailed description of the objects: amongst the AGN there are 2 head-tail, 10 core-lobe, 9 wide-angle-tail (WAT), 8 double-double or Z-/X-shaped, 3 bent-tail radio sources, and 26 symmetric sources, while amongst the SFGs we find the only star-forming ring seen in radio emission in COSMOS. We report a large number (32/58) of disturbed/bent multi-component AGN, 18 of which do not lie within X-ray groups in COSMOS (0.08 < z < 1.53). [abridged]

Telescope slew and settle time markedly reduces the efficiency of wide-field multi-epoch surveys for sensitive interferometers with small fields of view. The overheads can be mitigated through the use of On-the-Fly Mosaicing (OTFM), where the the antennas are driven at a non-sidereal rate and visibilities are recorded continuously. Here we introduce the OTFM technique for the VLA, and describe its implementation for the Caltech-NRAO Stripe 82 Survey (CNSS), a dedicated 5-epoch survey for slow transients at S band (2-4 GHz). We also describe the OTFSim tool for planning dynamically-scheduled OTFM observations on the VLA, the latest imaging capabilities for OTFM in CASA, and present a comparison of OTFM observations with pointed observations. Using the subset of our observations from the CNSS pilot and final surveys, we demonstrate that the wide-band and wide-field OTFM observations with the VLA can be imaged accurately, and that this technique offers a more efficient alternative to standard mosaicing for multi-epoch shallow surveys such as the CNSS and the VLA Sky Survey (VLASS). We envisage that the new OTFM mode will facilitate new synoptic surveys and high-frequency mapping experiments on the VLA.

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

The binary neutron star merger GW170817 was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 41+/-3 Mpc. The radio and X-ray afterglows of GW170817 exhibited delayed onset, a gradual rise in the emission with time as t^0.8, a peak at about 150 days post-merger, followed by a relatively rapid decline. To date, various models have been proposed to explain the afterglow emission, including a choked-jet cocoon and a successful-jet cocoon (a.k.a. structured jet). However, the observational data have remained inconclusive as to whether GW170817 launched a successful relativistic jet. Here we show, through Very Long Baseline Interferometry, that the compact radio source associated with GW170817 exhibits superluminal motion between two epochs at 75 and 230 days post-merger. This measurement breaks the degeneracy between the models and indicates that, while the early-time radio emission was powered by a wider-angle outflow (cocoon), the late-time emission was most likely dominated by an energetic and narrowly-collimated jet, with an opening angle of <5 degrees, and observed from a viewing angle of about 20 degrees. The imaging of a collimated relativistic outflow emerging from GW170817 adds substantial weight to the growing evidence linking binary neutron star mergers and short gamma-ray bursts.

We present multi-frequency monitoring observations of the black hole X-ray binary V404 Cygni throughout its June 2015 outburst. Our data set includes radio and mm/sub-mm photometry, taken with the Karl G. Jansky Very Large Array, Arc-Minute MicroKelvin Imager Large Array, Sub-millimeter Array, James Clerk Maxwell Telescope, and the Northern Extended Millimetre Array, combined with publicly available infrared, optical, UV, and X-ray measurements. With these data, we report detailed diagnostics of the spectral and variability properties of the jet emission observed during different stages of this outburst. These diagnostics show that emission from discrete jet ejecta dominated the jet emission during the brightest stages of the outburst. We find that the ejecta became fainter, slower, less frequent, and less energetic, before the emission transitioned (over 1--2 days) to being dominated by a compact jet, as the outburst decayed toward quiescence. While the broad-band spectrum of this compact jet showed very little evolution throughout the outburst decay (with the optically thick to thin synchrotron jet spectral break residing in the near-infrared/optical bands; $\sim2-5\times10^{14}$ Hz), the emission still remained intermittently variable at mm/sub-mm frequencies. Additionally, we present a comparison between the radio jet emission throughout the 2015 and previous 1989 outbursts, confirming that the radio emission in the 2015 outburst decayed significantly faster than in 1989. Lastly, we detail our sub-mm observations taken during the December 2015 mini-outburst of V404 Cygni, which demonstrate that, similar to the main outburst, the source was likely launching jet ejecta during this short period of renewed activity.

The Hubble constant ($H_0$) measures the current expansion rate of the Universe, and plays a fundamental role in cosmology. Tremendous effort has been dedicated over the past decades to measure $H_0$. Notably, Planck cosmic microwave background (CMB) and the local Cepheid-supernovae distance ladder measurements determine $H_0$ with a precision of $\sim 1\%$ and $\sim 2\%$ respectively. A $3$-$\sigma$ level of discrepancy exists between the two measurements, for reasons that have yet to be understood. Gravitational wave (GW) sources accompanied by electromagnetic (EM) counterparts offer a completely independent standard siren (the GW analogue of an astronomical standard candle) measurement of $H_0$, as demonstrated following the discovery of the neutron star merger, GW170817. This measurement does not assume a cosmological model and is independent of a cosmic distance ladder. The first joint analysis of the GW signal from GW170817 and its EM localization led to a measurement of $H_0=74^{+16}_{-8}$ km/s/Mpc (median and symmetric $68\%$ credible interval). In this analysis, the degeneracy in the GW signal between the source distance and the weakly constrained viewing angle dominated the $H_0$ measurement uncertainty. Recently, Mooley et al. (2018) obtained tight constraints on the viewing angle using high angular resolution imaging of the radio counterpart of GW170817. Here we obtain a significantly improved measurement $H_0=68.9^{+4.7}_{-4.6}$ km/s/Mpc by using these new radio observations, combined with the previous GW and EM data. We estimate that 15 more localized GW170817-like events (comparable signal-to-noise ratio, favorable orientation), having radio images and light curve data, will potentially bring resolution to the tension between the Planck and Cepheid-supernova measurements, as compared to 50-100 GW events without such data.

Calcium rich gap transients represent an intriguing new class of faint and fast evolving supernovae that exhibit strong [Ca II] emission in their nebular phase spectra. In this paper, we present the discovery and follow-up observations of iPTF 16hgs -- an intermediate luminosity and fast evolving transient that exhibited a double peaked light curve. Exhibiting a typical Type Ib spectrum in the photospheric phase and an early transition to a [Ca II] dominated nebular phase, we show that iPTF 16hgs shows properties consistent with the class of Ca-rich gap transients, with two interesting exceptions. First, while the second peak of the light curve is similar to other Ca-rich gap transients (suggesting $M_{ej}$ of 0.4 M$_\odot$ and peak luminosity of $3 \times 10^{41}$ ergs s$^{-1}$), we show that the first blue and fast declining (over $2$ days) peak is unique to this source. Second, with Integral Field Unit observations of the host galaxy, we find that iPTF 16hgs occurred in the outskirts (projected offset of $6$ kpc $ = 1.9 R_{eff}$) of a low metallicity (0.4 Z$_\odot$), star forming, dwarf spiral galaxy. Using deep late-time VLA and uGMRT observations, we place stringent limits on the local environment of the source, ruling out a large parameter space of circumstellar densities and mass loss environments of the progenitor. If iPTF 16hgs shares explosion physics with the class of Ca-rich gap transients, the presence of the first peak can be explained by enhanced mixing of 0.01 M$_\odot$ of $^{56}$Ni into the outer layers the ejecta, reminiscent of some models of He-shell detonations on WDs. On the other hand, if iPTF 16hgs is physically unrelated to the class, the first peak is consistent with shock cooling emission (of an envelope with a mass of 0.08 M$_\odot$ and radius of 13 R$_\odot$) associated with a core-collapse explosion of a highly stripped massive star in a close binary system.

We present late-time radio observations of GW170817, the first binary neutron star merger discovered through gravitational waves by the advanced LIGO and Virgo detectors. Our observations, carried out with the Karl G. Jansky Very Large Array, were optimized to detect polarized radio emission, and thus to constrain the linear polarization fraction of GW170817. At an epoch of ~244 days after the merger, we rule out linearly polarized emission above a fraction of ~12% at a frequency of 2.8 GHz (99% confidence). Within the structured jet scenario (a.k.a. successful jet plus cocoon system) for GW170817, the derived upper-limit on the radio continuum linear polarization fraction strongly constrains the magnetic field configuration in the shocked ejecta. We show that our results for GW170817 are compatible with the low level of linear polarization found in afterglows of cosmological long gamma-ray bursts. Finally, we discuss our findings in the context of future expectations for the study of radio counterparts of binary neutron star mergers identified by ground-based gravitational-wave detectors.

We present 2-9 GHz radio observations of GW170817 covering the period 125-200 days post-merger, taken with the Australia Telescope Compact Array and the Karl G. Jansky Very Large Array. Our observations demonstrate that the radio afterglow peaked at $149\pm2$ days post-merger and is now declining in flux density. We see no evidence for evolution in the radio-only spectral index, which remains consistent with optically-thin synchrotron emission connecting the radio, optical, and X-ray regimes. The peak implies a total energy in the synchrotron-emitting component of a ${\rm few}\times 10^{50}\,$erg. The temporal decay rate is most consistent with mildly- or non-relativistic material and we do not see evidence for a very energetic off-axis jet, but we cannot distinguish between a lower-energy jet and more isotropic emission.

We apply a newly-developed On-the-Fly mosaicing technique on the NSF's Karl G. Jansky Very Large Array (VLA) at 3 GHz in order to carry out a sensitive search for an afterglow from the Advanced LIGO binary black hole merger event GW 151226. In three epochs between 1.5 and 6 months post-merger we observed a 100 sq. deg region, with more than 80% of the survey region having a RMS sensitivity of better than 150 uJy/beam, in the northern hemisphere having a merger containment probability of 10%. The data were processed in near-real-time, and analyzed to search for transients and variables. No transients were found but we have demonstrated the ability to conduct blind searches in a time-frequency phase space where the predicted afterglow signals are strongest. If the gravitational wave event is contained within our survey region, the upper limit on any late-time radio afterglow from the merger event at an assumed mean distance of 440 Mpc is about 1e29 erg/s/Hz. Approximately 1.5% of the radio sources in the field showed variability at a level of 30%, and can be attributed to normal activity from active galactic nuclei. The low rate of false positives in the radio sky suggests that wide-field imaging searches at a few Gigahertz can be an efficient and competitive search strategy. We discuss our search method in the context of the recent afterglow detection from GW 170817 and radio follow-up in future gravitational wave observing runs.

We present the discovery by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) of a likely supernova (SN) in NGC 3556 at only 8.8 Mpc, which, despite its proximity, was not detected by optical searches. A luminous infrared (IR) transient at $M_{[4.5]} = -16.7$ mag (Vega), SPIRITS 16tn is coincident with a dust lane in the inclined, star-forming disk of the host. Using IR, optical, and radio observations, we attempt to determine the nature of this event. We estimate $A_V \approx$ 8 - 9 mag of extinction, placing it among the three most highly obscured IR-discovered SNe to date. The [4.5] light curve declined at a rate of 0.013 mag day$^{-1}$, and the $[3.6] - [4.5]$ color grew redder from 0.7 to $\gtrsim$ 1.0 mag by 184.7 days post discovery. Optical/IR spectroscopy shows a red continuum, but no clearly discernible features, preventing a definitive spectroscopic classification. Deep radio observations constrain the radio luminosity of SPIRITS 16tn to $L_{\nu} \lesssim 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ between 3 - 15 GHz, excluding many varieties of radio core-collapse SNe. A type Ia SN is ruled out by the observed red IR color, and lack of features normally attributed to Fe-peak elements in the optical and IR spectra. SPIRITS 16tn was fainter at [4.5] than typical stripped-envelope SNe by $\approx$ 1 mag. Comparison of the spectral energy distribution to SNe II suggests SPIRITS 16tn was both highly obscured, and intrinsically dim, possibly akin to the low-luminosity SN 2005cs. We infer the presence of an IR dust echo powered by a peak luminosity of the transient of $5 \times 10^{40}$ erg s$^{-1} < L_{\mathrm{peak}} < 4\times10^{43}$ erg s$^{-1}$, consistent with the observed range for SNe II. This discovery illustrates the power of IR surveys to overcome the compounding effects of visible extinction and optically sub-luminous events in completing the inventory of nearby SNe.

We report on late time radio and X-ray observations of the tidal disruption event candidate ASASSN-14li, covering the first 1000 days of the decay phase. For the first $\sim200$ days the radio and X-ray emission fade in concert. This phase is better fit by an exponential decay at X-ray wavelengths, while the radio emission is well described by either an exponential or the canonical $t^{-5/3}$ decay assumed for tidal disruption events. The correlation between radio and X-ray emission during this period can be fit as $L_{R}\propto L_{X}^{1.9\pm0.2}$. After 400 days the radio emission at $15.5\,\textrm{GHz}$ has reached a plateau level of $244\pm8\,\mu\textrm{Jy}$ which it maintains for at least the next 600 days, while the X-ray emission continues to fade exponentially. This steady level of radio emission is likely due to relic radio lobes from the weak AGN-like activity implied by historical radio observations. We note that while most existing models are based upon the evolution of ejecta which are decoupled from the central black hole, the radio : X-ray correlation during the declining phase is also consistent with core jet emission coupled to a radiatively efficient accretion flow.

Observing nearby galaxies would facilitate the search for artificial radio signals by sampling many billions of stars simultaneously, but few efforts have been made to exploit this opportunity. An added attraction is that the Milky Way is the second-largest member of the Local Group, so our galaxy might be a probable target for hypothetical broadcasters in nearby galaxies. We present the first relatively high spectral resolution (<1 kHz) 21 cm band search for intelligent radio signals of complete galaxies in the Local Group with the Jansky VLA, observing the galaxies M31 (Andromeda) and M33 (Triangulum) - the first and third largest members of the group respectively - sampling more stars than any prior search of this kind. We used 122 Hz channels over a 1 MHz spectral window in the target galaxy velocity frame of reference, and 15 Hz channels over a 125 kHz window in our local standard of rest. No narrowband signals were detected above a signal-to-noise ratio of 7, suggesting the absence of continuous narrowband flux greater than approximately 0.24 Jy and 1.33 Jy in the respective spectral windows illuminating our part of the Milky Way during our observations in December 2014 and January 2015. This is also the first study in which the upgraded VLA has been used for SETI.

We describe an image-based method that uses two radio criteria, compactness and spectral index, to identify promising pulsar candidates among Fermi Large Area Telescope (LAT) unassociated sources. These criteria are applied to those radio sources from the Giant Metrewave Radio Telescope all-sky survey at 150 MHz (TGSS ADR1) found within the error ellipses of unassociated sources from the 3FGL catalog and a preliminary source list based on 7 years of LAT data. After follow-up interferometric observations to identify extended or variable sources, a list of 16 compact, steep-spectrum candidates is generated. An ongoing search for pulsations in these candidates, in gamma rays and radio, has found six millisecond pulsars and one normal pulsar. A comparison of this method with existing selection criteria based on gamma-ray spectral and variability properties suggests that the pulsar discovery space using Fermi may be larger than previously thought. Radio imaging is a hitherto underutilized source selection method that can be used, as with other multi-wavelength techniques, in the search for Fermi pulsars.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

The evolution of extragalactic radio sources has been a fundamental problem in the study of active galactic nuclei for many years. A standard evolutionary model has been created based on observations of a wide range of radio sources. In the general scenario of the evolution, the younger and smaller Gigahertz-Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sources become large-scale FRI and FRII objects. However, a growing number of observations of low power radio sources suggests that the model cannot explain all their properties and there are still some aspects of the evolutionary path that remain unclear. There are indications, that some sources may be short-lived objects on timescales of $10^4$ - $10^5$ years. Those objects represent a new population of active galaxies. Here, we present the discovery of several radio transient sources on timescales of 5-20 years, largely associated with renewed AGN (Active Galactic Nucleus) activity. These changing-look AGNs possibly represent behaviour typical for many active galaxies.

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.

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.

Relativistic plasma jets are observed in many accreting black holes. According to theory, coiled magnetic fields close to the black hole accelerate and collimate the plasma, leading to a jet being launched. Isolating emission from this acceleration and collimation zone is key to measuring its size and understanding jet formation physics. But this is challenging because emission from the jet base cannot be easily disentangled from other accreting components. Here, we show that rapid optical flux variations from a Galactic black-hole binary are delayed with respect to X-rays radiated from close to the black hole by ~0.1 seconds, and that this delayed signal appears together with a brightening radio jet. The origin of these sub-second optical variations has hitherto been controversial. Not only does our work strongly support a jet origin for the optical variations, it also sets a characteristic elevation of <~10$^3$ Schwarzschild radii for the main inner optical emission zone above the black hole, constraining both internal shock and magnetohydrodynamic models. Similarities with blazars suggest that jet structure and launching physics could potentially be unified under mass-invariant models. Two of the best-studied jetted black hole binaries show very similar optical lags, so this size scale may be a defining feature of such systems.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

Gravitational waves have been detected from a binary neutron star merger event, GW170817. The detection of electromagnetic radiation from the same source has shown that the merger occurred in the outskirts of the galaxy NGC 4993, at a distance of 40 megaparsecs from Earth. We report the detection of a counterpart radio source that appears 16 days after the event, allowing us to diagnose the energetics and environment of the merger. The observed radio emission can be explained by either a collimated ultra-relativistic jet viewed off-axis, or a cocoon of mildly relativistic ejecta. Within 100 days of the merger, the radio light curves will distinguish between these models and very long baseline interferometry will have the capability to directly measure the angular velocity and geometry of the debris.

We present the Arcminute Microkelvin Imager (AMI) Large Array catalogue of 139 gamma-ray bursts (GRBs). AMI observes at a central frequency of 15.7 GHz and is equipped with a fully automated rapid-response mode, which enables the telescope to respond to high-energy transients detected by Swift. On receiving a transient alert, AMI can be on-target within two minutes, scheduling later start times if the source is below the horizon. Further AMI observations are manually scheduled for several days following the trigger. The AMI GRB programme probes the early-time (< 1 day) radio properties of GRBs, and has obtained some of the earliest radio detections (GRB 130427A at 0.36 and GRB 130907A at 0.51 days post-burst). As all Swift GRBs visible to AMI are observed, this catalogue provides the first representative sample of GRB radio properties, unbiased by multi-wavelength selection criteria. We report the detection of six GRB radio afterglows that were not previously detected by other radio telescopes, increasing the rate of radio detections by 50% over an 18-month period. The AMI catalogue implies a Swift GRB radio detection rate of >15%, down to ~0.2 mJy/beam. However, scaling this by the fraction of GRBs AMI would have detected in the Chandra & Frail sample (all radio-observed GRBs between 1997 - 2011), it is possible ~44 - 56% of Swift GRBs are radio-bright, down to ~0.1 - 0.15 mJy/beam. This increase from the Chandra & Frail rate (~30%) is likely due to the AMI rapid-response mode, which allows observations to begin while the reverse-shock is contributing to the radio afterglow.

We present multi-wavelength follow-up campaigns by the AstroSat-CZTI and GROWTH collaborations to search for an electromagnetic counterpart to the gravitational wave event GW170104. At the time of the GW170104 trigger, the AstroSat CZTI field-of-view covered 50.3\% of the sky localization. We do not detect any hard X-ray (>100 keV) signal at this time, and place an upper limit of $\approx 4.5 \times 10^{-7}~{\rm erg~cm}^{-2}{\rm~s}^{-1}$ for a 1\,s timescale. Separately, the ATLAS survey reported a rapidly fading optical source dubbed ATLAS17aeu in the error circle of GW170104. Our panchromatic investigation of ATLAS17aeu shows that it is the afterglow of an unrelated long, soft GRB~170105A, with only a fortuitous spatial coincidence with GW170104. We then discuss the properties of this transient in the context of standard long GRB afterglow models.

We present ground-based and \textit{Swift} observations of iPTF16fnl, a likely tidal disruption event (TDE) discovered by the intermediate Palomar Transient Factory (iPTF) survey at 66.6 Mpc. The lightcurve of the object peaked at absolute $M_g=-17.2$ mag. The maximum bolometric luminosity (from optical and UV) was $L_p~\simeq~(1.0\,\pm\,0.15) \times 10^{43}$ erg/s, an order of magnitude fainter than any other optical TDE discovered so far. The luminosity in the first 60 days is consistent with an exponential decay, with $L \propto e^{-(t-t_0)/\tau}$, where $t_0$=~57631.0 (MJD) and $\tau\simeq 15$ days. The X-ray shows a marginal detection at $L_X=2.4^{1.9}_{-1.1}\times 10^{39}$ erg/s (\textit{Swift} X-ray Telescope). No radio counterpart was detected down to 3$\sigma$, providing upper limits for monochromatic radio luminosity of $\nu L_{\nu} < 2.3\times10^{36}$ erg/s and $\nu L_{\nu}<1.7\times 10^{37}$ erg/s (VLA, 6.1 and 22 GHz). The blackbody temperature, obtained from combined \textit{Swift} UV and optical photometry, shows a constant value of 19,000 K. The transient spectrum at peak is characterized by broad He II and H$\alpha$ emission lines, with an FWHM of about 14,000 km/s and 10,000 km/s respectively. He I lines are also detected at $\lambda\lambda$ 5875 and 6678. The spectrum of the host is dominated by strong Balmer absorption lines, which are consistent with a post-starburst (E+A) galaxy with an age of $\sim$650 Myr and solar metallicity. The characteristics of iPTF16fnl make it an outlier on both luminosity and decay timescales, as compared to other optically selected TDEs. The discovery of such a faint optical event suggests a higher rate of tidal disruptions, as low luminosity events may have gone unnoticed in previous searches.

We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

Supermassive binary black holes (BBHs) on sub-parsec scales are prime targets for gravitational wave experiments. They also provide insights on close binary evolution and hierarchical structure formation. Sub-parsec BBHs cannot be spatially resolved but indirect methods can identify candidates. In 2015 Liu et al. reported an optical-continuum periodicity in the quasar PSO J334.2028+01.4075, with the estimated mass and rest-frame period suggesting an orbital separation of about 0.006 pc (0.7 micro-arcsec). The persistence of the quasar's optical periodicity has recently been disfavored over an extended baseline. However, if a radio jet is launched from a sub-parsec BBH, the binary's properties can influence the radio structure on larger scales. Here, we use the Very Long Baseline Array (VLBA) and Karl G. Jansky Very Large Array (VLA) to study the parsec- and kiloparsec-scale emission energized by the quasar's putative BBH. We find two VLBA components separated by 3.6 mas (30 pc), tentatively identifying one as the VLBA "core" from which the other was ejected. The VLBA components contribute to a point-like, time-variable VLA source that is straddled by lobes spanning 8 arcsec (66 kpc). We classify PSO J334.2028+01.4075 as a lobe-dominated quasar, albeit with an atypically large twist of 39 deg between its elongation position angles on parsec- and kiloparsec-scales. By analogy with 3C 207, a well-studied lobe-dominated quasar with a similarly-rare twist, we speculate that PSO J334.2028+01.4075 could be ejecting jet components over an inner cone that traces a precessing jet in a BBH system.

We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.

We present first results from a series of $NuSTAR$ observations of the black hole X-ray binary V404 Cyg obtained during its summer 2015 outburst, primarily focusing on observations during the height of this outburst activity. The $NuSTAR$ data show extreme variability in both the flux and spectral properties of the source. This is partly driven by strong and variable line-of-sight absorption, similar to previous outbursts. The latter stages of this observation are dominated by strong flares, reaching luminosities close to Eddington. During these flares, the central source appears to be relatively unobscured and the data show clear evidence for a strong contribution from relativistic reflection, providing a means to probe the geometry of the innermost accretion flow. Based on the flare properties, analogy with other Galactic black hole binaries, and also the simultaneous onset of radio activity, we argue that this intense X-ray flaring is related to transient jet activity during which the ejected plasma is the primary source of illumination for the accretion disk. If this is the case, then our reflection modelling implies that these jets are launched in close proximity to the black hole (as close as a few gravitational radii), consistent with expectations for jet launching models that tap either the spin of the central black hole, or the very innermost accretion disk. Our analysis also allows us to place the first constraints on the black hole spin for this source, which we find to be $a^* > 0.92$ (99% statistical uncertainty, based on an idealized lamppost geometry).

The Fermi Large Area Telescope has detected an extended region of GeV emission toward the Galactic Center that is currently thought to be powered by dark matter annihilation or a population of young and/or millisecond pulsars. In a test of the pulsar hypothesis, we have carried out an initial search of a 20 deg**2 area centered on the peak of the galactic center GeV excess. Candidate pulsars were identified as a compact, steep spectrum continuum radio source on interferometric images and followed with targeted single-dish pulsation searches. We report the discovery of the recycled pulsar PSR 1751-2737 with a spin period of 2.23 ms. PSR 1751-2737 appears to be an isolated recycled pulsar located within the disk of our Galaxy, and it is not part of the putative bulge population of pulsars that are thought to be responsible for the excess GeV emission. However, our initial success in this small pilot survey suggests that this hybrid method (i.e. wide-field interferometric imaging followed up with single dish pulsation searches) may be an efficient alternative strategy for testing whether a putative bulge population of pulsars is responsible for the GeV excess.

We present a radio light curve of supernova (SN) 2014C taken with the Arcminute Microkelvin Imager (AMI) Large Array at 15.7 GHz. Optical observations presented by Milisavljevic et al. demonstrated that SN 2014C metamorphosed from a stripped-envelope Type Ib SN into a strongly interacting Type IIn SN within 1 year. The AMI light curve clearly shows two distinct radio peaks, the second being a factor of 4 times more luminous than the first peak. This double bump morphology indicates two distinct phases of mass-loss from the progenitor star with the transition between density regimes occurring at 100-200 days. This reinforces the interpretation that SN 2014C exploded in a low density region before encountering a dense Hydrogen-rich shell of circumstellar material that was likely ejected by the progenitor prior to the explosion. The AMI flux measurements of the first light curve bump are the only reported observations taken within ~50 to ~125 days post-explosion, before the blast-wave encountered the Hydrogen shell. Simplistic synchrotron self-absorption (SSA) and free-free absorption (FFA) modelling suggest that some physical properties of SN 2014C, such as the mass-loss rate, are consistent with the properties of other Type Ibc and IIn SNe. However, our single frequency data does not allow us to distinguish between these two models, which implies they are likely too simplistic to describe the complex environment surrounding this event. Lastly, we present the precise radio location of SN 2014C obtained with eMERLIN, which will be useful for future VLBI observations of the SN.

We present a radio-quiet quasar at z=0.237 discovered "turning on" by the intermediate Palomar Transient Factory (iPTF). The transient, iPTF 16bco, was detected by iPTF in the nucleus of a galaxy with an archival SDSS spectrum with weak narrow-line emission characteristic of a low-ionization emission line region (LINER). Our follow-up spectra show the dramatic appearance of broad Balmer lines and a power-law continuum characteristic of a luminous (L_bol~10^45 erg/s) type 1 quasar 12 years later. Our photometric monitoring with PTF from 2009-2012, and serendipitous X-ray observations from the XMM-Newton Slew Survey in 2011 and 2015, constrain the change of state to have occurred less than 500 days before the iPTF detection. An enhanced broad Halpha to [OIII]5007 line ratio in the type 1 state relative to other changing-look quasars also is suggestive of the most rapid change of state yet observed in a quasar. We argue that the >10 increase in Eddington ratio inferred from the brightening in UV and X-ray continuum flux is more likely due to an intrinsic change in the accretion rate of a pre-existing accretion disk, than an external mechanism such as variable obscuration, microlensing, or the tidal disruption of a star. However, further monitoring will be helpful in better constraining the mechanism driving this change of state. The rapid "turn on" of the quasar is much shorter than the viscous infall timescale of an accretion disk, and requires a disk instability that can develop around a ~10^8 M_sun black hole on timescales less than a year.

The connection between accretion and jet production in accreting white dwarf binary systems, especially dwarf novae, is not well understood. Radio wavelengths provide key insights into the mechanisms responsible for accelerating electrons, including jets and outflows. Here we present densely-sampled radio coverage, obtained with the Arcminute MicroKelvin Imager Large Array, of the dwarf nova SS Cyg during its February 2016 anomalous outburst. The outburst displayed a slower rise (3 days per mag) in the optical than typical ones, and lasted for more than 3 weeks. Rapid radio flaring on timescales <1 hour was seen throughout the outburst. The most intriguing behavior in the radio was towards the end of the outburst where a fast, luminous ("giant"), flare peaking at ~20 mJy and lasting for 15 minutes was observed. This is the first time that such a flare has been observed in SS Cyg, and insufficient coverage could explain its non-detection in previous outbursts. These data, together with past radio observations, are consistent with synchrotron emission from plasma ejection events as being the origin of the radio flares. However, the production of the giant flare during the declining accretion rate phase remains unexplained within the standard accretion-jet framework and appears to be markedly different to similar patterns of behavior in X-ray binaries.

After more than 26 years in quiescence, the black hole transient V404 Cyg went into a luminous outburst in June 2015, and additional activity was detected in late December of the same year. Here, we present an optical spectroscopic follow-up of the December mini-outburst, together with X-ray, optical and radio monitoring that spanned more than a month. Strong flares with gradually increasing intensity are detected in the three spectral ranges during the ~10 days following the Swift trigger. Our optical spectra reveal the presence of a fast outflowing wind, as implied by the detection of a P-Cyg profile (He I - 5876 A) with a terminal velocity of ~2500 km/s. Nebular-like spectra -- with an H_alpha equivalent width of ~500 A -- are also observed. All these features are similar to those seen during the main June 2015 outburst. Thus, the fast optical wind simultaneous with the radio jet is most likely present in every V404 Cyg outburst. Finally, we report on the detection of a strong radio flare in late January 2016, when X-ray and optical monitoring had stopped due to Sun constraints.

We present the first full release of a survey of the 150 MHz radio sky, observed with the Giant Metrewave Radio Telescope between April 2010 and March 2012 as part of the TGSS project. Aimed at producing a reliable compact source survey, our automated data reduction pipeline efficiently processed more than 2000 hours of observations with minimal human interaction. Through application of innovative techniques such as image-based flagging, direction-dependent calibration of ionospheric phase errors, correcting for systematic offsets in antenna pointing, and improving the primary beam model, we created good quality images for over 95 percent of the 5336 pointings. Our data release covers 36,900 square degrees (or 3.6 pi steradians) of the sky between -53 deg and +90 deg DEC, which is 90 percent of the total sky. The majority of pointing images have a background RMS noise below 5 mJy/beam with an approximate resolution of 25" x 25" (or 25" x 25" / cos (DEC - 19 deg) for pointings south of 19 deg DEC). We have produced a catalog of 0.62 Million radio sources derived from an initial, high reliability source extraction at the 7 sigma level. For the bulk of the survey, the measured overall astrometric accuracy is better than 2" in RA and DEC, while the flux density accuracy is estimated at ~10 percent. Within the scope of the TGSS ADR project, the source catalog as well as 5336 mosaic images (5 deg x 5 deg) and an image cutout service, are made publicly available online as a service to the astronomical community. Next to enabling a wide range of different scientific investigations, we anticipate that these survey products provide a solid reference for various new low-frequency radio aperture array telescopes (LOFAR, LWA, MWA, SKA-low), and can play an important role in characterizing the EoR foreground. The TGSS ADR project aims at continuously improving the quality of the survey data products.

We present optical and near-IR linear polarimetry of V404 Cyg during its 2015 outburst and in quiescence. We obtained time resolved r'-band polarimetry when the source was in outburst, near-IR polarimetry when the source was near quiescence and multiple wave-band optical polarimetry later in quiescence. The optical to near-IR linear polarization spectrum can be described by interstellar dust and an intrinsic variable component. The intrinsic optical polarization, detected during the rise of one of the brightest flares of the outburst, is variable, peaking at 4.5 per cent and decaying to 3.5 per cent. We present several arguments that favour a synchrotron jet origin to this variable polarization, with the optical emission originating close to the jet base. The polarization flare occurs during the initial rise of a major radio flare event that peaks later, and is consistent with a classically evolving synchrotron flare from an ejection event. We conclude that the optical polarization flare represents a jet launching event; the birth of a major ejection. For this event we measure a rather stable polarization position angle of -9 degrees E of N, implying that the magnetic field near the base of the jet is approximately perpendicular to the jet axis. This may be due to the compression of magnetic field lines in shocks in the accelerated plasma, resulting in a partially ordered transverse field that have now been seen during the 2015 outburst. We also find that this ejection occurred at a similar stage in the repetitive cycles of flares.

We present the first results from an on-going survey to characterize the circumgalactic medium (CGM) of the massive high-redshift galaxies detected as submillimeter galaxies (SMGs). We constructed a parent sample of 163 SMG-QSO pairs with separations less than $\sim$36 arcsec by cross-matching far-infrared-selected galaxies from Herschel with spectroscopically confirmed QSOs. The Herschel sources were selected to match the properties of SMGs. We determined the sub-arcsecond positions of six Herschel sources with the Very Large Array and obtained secure redshift identification for three of those with near-infrared spectroscopy. The QSO sightlines probe transverse proper distances of 112, 157, and 198 kpc at foreground redshifts of 2.043, 2.515, and 2.184, respectively, which are comparable to the virial radius of the $\sim10^{13}$ Msun halos expected to host SMGs. High-quality absorption-line spectroscopy of the QSOs reveals systematically strong HI Lyman-alpha absorption around all three SMGs, with rest-frame equivalent widths of $\sim2-3$ \AA. However, none of the three absorbers exhibits compelling evidence for optically thick HI gas or metal absorption, in contrast to the dominance of strong neutral absorbers in the CGM of luminous $z \sim 2$ QSOs. The low covering factor of optically thick HI gas around SMGs tentatively indicates that SMGs may not have as prominent cool gas reservoirs in their halos as the co-eval QSOs and that they may inhabit less massive halos than previously thought.

We report on Fermi/Large Area Telescope observations of the accreting black hole low-mass X-ray binary V404 Cygni during its outburst in June-July 2015. Detailed analyses reveal a possible excess of $\gamma$-ray emission on 26 June 2015, with a very soft spectrum above $100$ MeV, at a position consistent with the direction of V404 Cyg (within the $95\%$ confidence region and a chance probability of $4 \times 10^{-4}$). This emission cannot be associated with any previously-known Fermi source. Its temporal coincidence with the brightest radio and hard X-ray flare in the lightcurve of V404 Cyg, at the end of the main active phase of its outburst, strengthens the association with V404 Cyg. If the $\gamma$-ray emission is associated with V404 Cyg, the simultaneous detection of $511\,$keV annihilation emission by INTEGRAL requires that the high-energy $\gamma$ rays originate away from the corona, possibly in a Blandford-Znajek jet. The data give support to models involving a magnetically-arrested disk where a bright $\gamma$-ray jet can re-form after the occurrence of a major transient ejection seen in the radio.

We present observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz toward a sub-field of the XXL-North 25 deg$^{2}$ field targeting the first supercluster discovered in the XXL Survey. The structure has been found at a spectroscopic redshift of 0.43 and extending over $0.35\times0.1$ deg$^{2}$ on the sky. We present the 3 GHz VLA radio continuum observations, the radio mosaic and radio source catalogue, and, perform an analysis of the supercluster in the optical and radio regimes using photometric redshifts from the CFHTLS survey and our new VLA data. Our 3 GHz radio mosaic has a resolution of $3.2\times1.9$ arcsec$^2$, and covers an area of $41\times41$ arcmin$^2$ with rms noise level below $\sim20 \mu$Jy beam$^{-1}$ ($\sim11\mu$Jy beam$^{-1}$ in central $15\times15$ arcmin$^2$). We extract a catalogue of 155 radio sources $\geq6\sigma$, 8 of which are multicomponent sources, and 123 ($79\%$) can be associated with optical sources in the CFHTLS W1 catalogue. Applying Voronoi tessellation analysis (VTA) in the area around the X-ray identified supercluster we identify a total of 17 overdensities at $z=0.35-0.50$, 7 of which are associated with clusters detected in the $XMM-Newton$ XXL data. We find a mean photometric redshift of 0.43 for our overdensities, consistent with the spectroscopic redshifts of 7 X-ray detected clusters. The full VTA-identified structure extends over $\sim0.6\times0.2$ deg$^2$ on the sky, ($\sim12\times4$ Mpc$^{2}$ at $z=0.43$). We associate 8 radio sources with potential group/cluster member galaxies none of which are large radio galaxies. The spatial distribution of the red and blue potential group member galaxies suggests that the clusters are dynamically young and not virialised yet, as expected for hierarchical structure growth in a $\Lambda$CDM universe. Further spectroscopic data are required to analyze the dynamical state of the groups.

We report on a search for steep spectrum radio sources within the 95% confidence error ellipses of the Fermi unassociated sources from the Large Array Telescope (LAT). Using existing catalogs and the newly released GMRT all-sky survey at 150 MHz we identify compact radio sources that are bright at MHz frequencies but faint or absent at GHz frequencies. Such steep spectrum radio sources are rare and constitute a sample of pulsar candidates, selected independently of period, dispersion measure, interstellar scattering and orbital parameters. We find point-like, steep spectrum candidates toward 11 Fermi sources. Based on the gamma-ray/radio positional coincidence, the rarity of such radio sources, and the properties of the 3FGL sources themselves, we argue that many of these sources could be pulsars. They may have been missed by previous radio periodicity searches due to interstellar propagation effects or because they lie in an unusually tight binary. If this hypothesis is correct, then renewed gamma-ray and radio periodicity searches at the positions of the steep spectrum radio sources may reveal pulsations.

A radio source that faded over six days, with a redshift of $z\approx0.5$ host, has been identified by Keane et al. (2016) as the transient afterglow to a fast radio burst (FRB 150418). We report follow-up radio and optical observations of the afterglow candidate and find a source that is consistent with an active galactic nucleus. If the afterglow candidate is nonetheless a prototypical FRB afterglow, existing slow-transient surveys limit the fraction of FRBs that produce afterglows to 0.25 for afterglows with fractional variation, $m=2|S_1-S_2|/(S_1+S_2)\geq0.7$, and 0.07 for $m\geq1$, at 95% confidence. In anticipation of a barrage of bursts expected from future FRB surveys, we provide a simple framework for statistical association of FRBs with afterglows. Our framework properly accounts for statistical uncertainties, and ensures consistency with limits set by slow-transient surveys.

We present the first dedicated radio continuum survey of a Kepler K2 mission field, Field 1 covering the North Galactic Cap. The survey is wide field, contemporaneous, multi-epoch, and multi-resolution in nature and was conducted at low radio frequencies between 140 and 200 MHz. The multi-epoch and ultra wide field (but relatively low resolution) part of the survey was provided by 15 nights of observation with the Murchison Widefield Array (MWA) over a period of approximately a month, contemporaneous with K2 observations of the field. The multi-resolution aspect of the survey was provided by the low resolution (4') MWA imaging, complemented by non-contemporaneous but much higher resolution (20") observations using the Giant Metrewave Radio Telescope (GMRT). The survey is therefore sensitive to the details of radio structures across a wide range of angular scales. Consistent with other recent low radio frequency surveys, no significant radio transients or variables were detected in the survey. The resulting source catalogs consist of 1,085 and 1,468 detections in the two MWA observation bands (centered at 154 and 185 MHz, respectively) and 7,445 detections in the GMRT observation band (centered at 148 MHz), over 314 square degrees. The survey is presented as a significant resource for multi-wavelength investigations of the more than 21,000 target objects in the K2 field. We briefly examine our survey data against K2 target lists for dwarf star types (stellar types M and L) that have been known to produce radio flares.

We have used the 150 MHz radio continuum survey (TGSS ADR) from the Giant Metrewave Radio Telescope (GMRT) to search for phase-averaged emission toward all well-localized radio pulsars north of -53deg Declination. We detect emission toward 200 pulsars with high confidence (>=5-sigma) and another 88 pulsars at fainter levels. We show that most of our identifications are likely from pulsars, except for a small number where the measured flux density is confused by an associated supernova or pulsar-wind nebula, or a globular cluster. We investigate the radio properties of the 150 MHz sample and we find an unusually high number of gamma-ray binary millisecond pulsars with very steep spectral indices. We also note a discrepancy in the measured flux densities between GMRT and LOFAR pulsar samples, suggesting that the flux density scale for the LOFAR pulsar sample may be in error by approximately a factor two. We carry out a separate search of 30 well-localized gamma-ray, radio-quiet pulsars in an effort to detect a widening of the radio beam into the line-of-sight at lower frequencies. No steep spectrum emission was detected either toward individual pulsars or in a weighted stack of all 30 images.

Accretion of matter onto black holes is universally associated with strong radiative feedback and powerful outflows. In particular, black hole transients show outflows whose properties are strongly coupled to those of the accretion flow. This includes X-ray winds of ionized material, expelled from the accretion disc encircling the black hole, and collimated radio jets. Very recently, a distinct optical variability pattern has been reported in the transient black hole transient V404 Cyg, and interpreted as disrupted mass flow into the inner regions of its large accretion disc. Here, we report on the discovery of a sustained outer accretion disc wind in V404 Cyg, which is unlike any seen previously. We find that the outflowing wind is neutral, has a large covering factor, expands at 1% of the speed of light and triggers a nebular phase once accretion sharply drops and the ejecta become optically thin. The large expelled mass (> 10^-8 Msun) indicates that the outburst was prematurely ended when a sizeable fraction of the outer disc was depleted by the wind, detaching the inner regions from the rest of the disc. The luminous, but brief, accretion phases shown by transients with large accretion discs imply that this outflow is most likely a new fundamental ingredient regulating mass accretion onto black holes.

Recently, studies of brown dwarfs have demonstrated that they possess strong magnetic fields and have the potential to produce radio and optical auroral emissions powered by magnetospheric currents. This emission provides the only window on magnetic fields in the coolest brown dwarfs and identifying additional benchmark objects is key to constraining dynamo theory in this regime. To this end, we conducted a new red optical (6300 - 9700 Angstrom) survey with the Keck telescopes looking for H-alpha emission from a sample of late L dwarfs and T dwarfs. Our survey gathered optical spectra for 29 targets, 18 of which did not have previous optical spectra in the literature, greatly expanding the number of moderate resolution (R~2000) spectra available at these spectral types. Combining our sample with previous surveys, we confirm an H-alpha detection rate of 9.2 (+3.5/-2.1) % for L and T dwarfs in the optical spectral range of L4 - T8. This detection rate is consistent with the recently measured detection rate for auroral radio emission from Kao et al. (2016), suggesting that geometrical selection effects due to the beaming of the radio emission are small or absent. We also provide the first detection of H-alpha emission from 2MASS 0036+1821, previously notable as the only electron cyclotron maser radio source without a confirmed detection of H-alpha emission. Finally, we also establish optical standards for spectral types T3 and T4, filling in the previous gap between T2 and T5.

We present observations of rapid (sub-second) optical flux variability in V404 Cyg during its 2015 June outburst. Simultaneous three-band observations with the ULTRACAM fast imager on four nights show steep power spectra dominated by slow variations on ~100-1000s timescales. Near the peak of the outburst on June 26, a dramatic change occurs and additional, persistent sub-second optical flaring appears close in time to giant radio and X-ray flaring. The flares reach peak optical luminosities of ~few x 10^{36} erg/s. Some are unresolved down to a time resolution of 24 milliseconds. Whereas the fast flares are stronger in the red, the slow variations are bluer when brighter. The redder slopes, emitted power, and characteristic timescales of the fast flares can be explained as optically-thin synchrotron emission from a compact jet arising on size scales ~140-500 Gravitational radii (with a possible additional contribution by a thermal particle distribution). The origin of the slower variations is unclear. The optical continuum spectral slopes are strongly affected by dereddening uncertainties and contamination by strong H{alpha} emission, but the variations of these slopes follow relatively stable loci as a function of flux. Cross-correlating the slow variations between the different bands shows asymmetries on all nights consistent with a small red skew (i.e., red lag). X-ray reprocessing and non-thermal emission could both contribute to these. These data reveal a complex mix of components over five decades in timescale during the outburst.

We have commenced a multi-year program, the Caltech-NRAO Stripe 82 Survey (CNSS), to search for radio transients with the Jansky VLA in the SDSS Stripe 82 region. The CNSS will deliver five epochs over the entire $\sim$270 deg$^2$ of Stripe 82, an eventual deep combined map with a rms noise of $\sim$40 $\mu$Jy and catalogs at a frequency of 3 GHz, and having a spatial resolution of 3". This first paper presents the results from an initial pilot survey of a 50 deg$^2$ region of Stripe 82, involving four epochs spanning logarithmic timescales between one week and 1.5 years, with the combined map having a median rms noise of 35 $\mu$Jy. This pilot survey enabled the development of the hardware and software for rapid data processing, as well as transient detection and follow-up, necessary for the full 270 deg$^2$ survey. Classification of variable and transient sources relied heavily on the wealth of multi-wavelength data in the Stripe 82 region, supplemented by repeated mapping of the region by the Palomar Transient Factory. $3.9^{+0.5}_{-0.9}$% of the detected point sources were found to vary by greater than 30%, consistent with similar studies at 1.4 GHz and 5 GHz. Multi-wavelength photometric data and light curves suggest that the variability is mostly due to shock-induced flaring in the jets of AGN. Although this was only a pilot survey, we detected two bona fide transients, associated with an RS CVn binary and a dKe star. Comparison with existing radio survey data revealed additional highly variable and transient sources on timescales between 5-20 years, largely associated with renewed AGN activity. The rates of such AGN possibly imply episodes of enhanced accretion and jet activity occurring once every $\sim$40,000 years in these galaxies. We compile the revised radio transient rates and make recommendations for future transient surveys and joint radio-optical experiments. (Abridged)

We present an X-ray analysis of a 4 Gyr old open cluster, M67, using archival XMM-Newton data. The aim of this study was to find new X-ray members of M67, and to use the updated member list for studying X-ray variability, and derive the X-ray luminosity functions (XLFs) of different stellar types and compare them with other star clusters of similar age. We report the detection of X-ray emission from 25 members of M67, with membership based primarily on their proper motion, of which one X-ray source is a new member. Supplementing this study with previous ROSAT and Chandra studies of M67, and using the most recent proper motion study by Vereshchagin et al., we have compiled a revised list of X-ray emitting members of M67 consisting of 43 stars. Sixteen of these are known RS CVn type binaries having orbital periods $<$ 10 days, and near-circular orbits, 5 are contact binaries with orbital periods $<$ 6 hours, 5 are yellow and blue stragglers, 2 are Algol-type binaries, and one source is a cataclysmic variable. Fourteen members do not have any orbital information and cannot be classified. Fourteen of the X-ray sources detected do not have any optical counterpart down to a magnitude of $V\simeq22$, and their membership is uncertain. Finally, we report the X-ray luminosity functions of RS CVn type and other types of stars in M67 and compare them with other open clusters of intermediate-to-old age.

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.

We present deep 3 GHz VLA observations of the potentially recoiling black hole CID-42 in the COSMOS field. This galaxy shows two optical nuclei in the HST/ACS image and a large velocity offset of ~ 1300 km/s between the broad and narrow H beta emission line although the spectrum is not spacially resolved (Civano et al. 2010). The new 3 GHz VLA data has a bandwidth of 2 GHz and to correctly interpret the flux densities imaging was done with two different methods: multi-scale multi-frequency synthesis and spectral windows stacking. The final resolutions and sensitivities of these maps are 0.7" with rms = 4.6 muJy/beam and 0.9" with rms = 4.8 muJy/beam respectively. With a 7 sigma detection we find that the entire observed 3 GHz radio emission can be associated with the South-Eastern component of CID-42, coincident with the detected X-ray emission. We use our 3 GHz data combined with other radio data from the literature ranging from 320 MHz to 9 GHz, which include the VLA, VLBA and GMRT data, to construct a radio synchrotron spectrum of CID-42. The radio spectrum suggests a type I unobscured radio-quiet flat-spectrum AGN in the South-Eastern component which may be surrounded by a more extended region of old synchrotron electron population or shocks generated by the outflow from the supermassive black hole. Our data are consistent with the recoiling black hole picture but cannot rule out the presence of an obscured and radio-quiet SMBH in the North-Western component.

We present the discovery, classification, and extensive panchromatic (from radio to X-ray) follow-up observations of PTF11qcj, a supernova discovered by the Palomar Transient Factory. PTF11qcj is located at a distance of dL ~ 124 Mpc. Our observations with the Karl G. Jansky Very Large Array show that this event is radio-loud: PTF11qcj reached a radio peak luminosity comparable to that of the famous gamma-ray-burst-associated supernova 1998bw (L_{5GHz} ~ 10^{29} erg/s/Hz). PTF11qcj is also detected in X-rays with the Chandra observatory, and in the infrared band with Spitzer. Our multi-wavelength analysis probes the supernova interaction with circumstellar material. The radio observations suggest a progenitor mass-loss rate of ~10^{-4} Msun/yr x (v_w/1000 km/s), and a velocity of ~(0.3-0.5)c for the fastest moving ejecta (at ~10d after explosion). However, these estimates are derived assuming the simplest model of supernova ejecta interacting with a smooth circumstellar material characterized by radial power-law density profile, and do not account for possible inhomogeneities in the medium and asphericity of the explosion. The radio light curve shows deviations from such a simple model, as well as a re-brightening at late times. The X-ray flux from PTF11qcj is compatible with the high-frequency extrapolation of the radio synchrotron emission (within the large uncertainties). An IR light echo from pre-existing dust is in agreement with our infrared data. Our analysis of pre-explosion data from the Palomar Transient Factory suggests that a precursor eruption of absolute magnitude M_r ~ -13 mag may have occurred ~ 2.5 yr prior to the supernova explosion. Based on our panchromatic follow-up campaign, we conclude that PTF11qcj fits the expectations from the explosion of a Wolf-Rayet star. Precursor eruptions may be a feature characterizing the final pre-explosion evolution of such stars.

Advanced LIGO and Advanced Virgo will be all-sky monitors for merging compact objects within a few hundred Mpc. Finding the electromagnetic counterparts to these events will require an understanding of the transient sky at low red-shift (z<0.1). We performed a systematic search for extragalactic, low red-shift, transient events in the XMM-Newton Slew Survey. In a flux limited sample, we found that highly-variable objects comprised 10% of the sample, and that of these, 10% were spatially coincident with cataloged optical galaxies. This led to 4x10^-4 transients per square degree above a flux threshold of 3x10^-12 erg cm-2 s-1 [0.2-2 keV] which might be confused with LIGO/Virgo counterparts. This represents the first extragalactic measurement of the soft X-ray transient rate within the Advanced LIGO/Virgo horizon. Our search revealed six objects that were spatially coincident with previously cataloged galaxies, lacked evidence for optical AGNs, displayed high luminosities around 10^43 erg s-1, and varied in flux by more than a factor of ten when compared with the ROSAT All-Sky Survey. At least four of these displayed properties consistent with previously observed tidal disruption events.

We describe the results of a search for early-type stars associated with the Taurus-Auriga molecular cloud complex, a diffuse nearby star-forming region noted as lacking young stars of intermediate and high mass. We investigate several sets of possible O, B and early A spectral class members. The first is a group of stars for which mid-infrared images show bright nebulae, all of which can be associated with stars of spectral type B. The second group consists of early-type stars compiled from (i) literature listings in SIMBAD; (ii) B stars with infrared excesses selected from the Spitzer Space Telescope survey of the Taurus cloud; (iii) magnitude- and color-selected point sources from the 2MASS; and (iv) spectroscopically identified early-type stars from the SDSS coverage of the Taurus region. We evaluated stars for membership in the Taurus-Auriga star formation region based on criteria involving: spectroscopic and parallactic distances, proper motions and radial velocities, and infrared excesses or line emission indicative of stellar youth. For selected objects, we also model the scattered and emitted radiation from reflection nebulosity and compare the results with the observed spectral energy distributions to further test the plausibility of physical association of the B stars with the Taurus cloud. This investigation newly identifies as probable Taurus members three B-type stars: HR 1445 (HD 28929), tau Tau (HD 29763), 72 Tau (HD 28149), and two A-type stars: HD 31305 and HD 26212, thus doubling the number of stars A5 or earlier associated with the Taurus clouds. Several additional early-type sources including HD 29659 and HD 283815 meet some, but not all, of the membership criteria and therefore are plausible, though not secure, members.

We report on an analysis of the Extended Chandra Deep Field South (E-CDFS) region using archival data from the Very Large Array, with the goal of studying radio variability and transients at the sub-mJy level. The 49 epochs of E-CDFS observations at 1.4 GHz sample timescales from one day to 3 months. We find that only a fraction (1%) of unresolved radio sources above 40 uJy are variable at the 4-sigma level. There is no evidence that the fractional variability changes along with the known transition of radio source populations below one milliJansky. Optical identifications of the sources show that the variable radio emission is associated with the central regions of an active galactic nucleus or a star-forming galaxy. After a detailed comparison of the efficacy of various source-finding algorithms, we use the best to carry out a transient search. No transients were found. This implies that the areal density of transients with peak flux density greater than 0.21 mJy is less than 0.37 deg^-2 (at a confidence level of 95%). This result is approximately an order of magnitude below the transient rate measured at 5 GHz by Bower et al. (2007) but it is consistent with more recent upper limits from Frail et al. (2012). Our findings suggest that the radio sky at 1.4 GHz is relatively quiet. For multi-wavelength transient searches, such as the electromagnetic counterparts to gravitational waves, this frequency may be optimal for reducing the high background of false positives.

Active galactic nuclei (AGN), powered by long-term accretion onto central supermassive black holes, produce relativistic jets with lifetimes of greater than one million yr that preclude observations at birth. Transient accretion onto a supermassive black hole, for example through the tidal disruption of a stray star, may therefore offer a unique opportunity to observe and study the birth of a relativistic jet. On 2011 March 25, the Swift {\gamma}-ray satellite discovered an unusual transient source (Swift J164449.3+573451) potentially representing such an event. Here we present the discovery of a luminous radio transient associated with Swift J164449.3+573451, and an extensive set of observations spanning centimeter to millimeter wavelengths and covering the first month of evolution. These observations lead to a positional coincidence with the nucleus of an inactive galaxy, and provide direct evidence for a newly-formed relativistic outflow, launched by transient accretion onto a million solar mass black hole. While a relativistic outflow was not predicted in this scenario, we show that the tidal disruption of a star naturally explains the high-energy properties, radio luminosity, and the inferred rate of such events. The weaker beaming in the radio compared to {\gamma}-rays/X-rays, suggests that radio searches may uncover similar events out to redshifts of z ~ 6.