We present MeerKAT 1.28 GHz total-intensity, polarization, and spectral-index images covering the giant (projected length $l \approx 1.57$~Mpc) X-shaped radio source PKS~2014$-$55 with an unprecedented combination of brightness sensitivity and angular resolution. They show the clear "double boomerang" morphology of hydrodynamical backflows from the straight main jets deflected by the large and oblique hot-gas halo of the host galaxy PGC~064440. The magnetic field orientation in PKS~2014$-$55 follows the flow lines from the jets through the secondary wings. The radio source is embedded in faint ($T_\mathrm{b} \approx 0.5 \mathrm{\,K}$) cocoons having the uniform brightness temperature and sharp outer edges characteristic of subsonic expansion into the ambient intra-group medium. The position angle of the much smaller ($l \sim 25$~kpc) restarted central source is within $5^\circ$ of the main jets, ruling out models that invoke jet re-orientation or two independent jets. Compression and turbulence in the backflows probably produce the irregular and low polarization bright region behind the apex of each boomerang as well as several features in the flow with bright heads and dark tails.

We present a measurement of the Hubble constant made using geometric distance measurements to megamaser-hosting galaxies. We have applied an improved approach for fitting maser data and obtained better distance estimates for four galaxies previously published by the Megamaser Cosmology Project: UGC 3789, NGC 6264, NGC 6323, and NGC 5765b. Combining these updated distance measurements with those for the maser galaxies CGCG 074-064 and NGC 4258, and assuming a fixed velocity uncertainty of 250 km s$^{-1}$ associated with peculiar motions, we constrain the Hubble constant to be $H_0 = 73.9 \pm 3.0$ km s$^{-1}$ Mpc$^{-1}$ independent of distance ladders and the cosmic microwave background. This best value relies solely on maser-based distance and velocity measurements, and it does not use any peculiar velocity corrections. Different approaches for correcting peculiar velocities do not modify $H_0$ by more than ${\pm}1{\sigma}$, with the full range of best-fit Hubble constant values spanning 71.8-76.9 km s$^{-1}$ Mpc$^{-1}$. We corroborate prior indications that the local value of $H_0$ exceeds the early-Universe value, with a confidence level varying from 95-99% for different treatments of the peculiar velocities.

As part of the survey component of the Megamaser Cosmology Project, we have discovered a disk megamaser system in the galaxy CGCG 074-064. Using the GBT and the VLA, we have obtained spectral monitoring observations of this maser system at a monthly cadence over the course of two years. We find that the systemic maser features display line-of-sight accelerations of ~4.4 km s$^{-1}$ yr$^{-1}$ that are nearly constant with velocity, while the high-velocity maser features show accelerations that are consistent with zero. We have also used the HSA to make a high-sensitivity VLBI map of the maser system in CGCG 074-064, which reveals that the masers reside in a thin, edge-on disk with a diameter of ~1.5 mas (0.6 pc). Fitting a three-dimensional warped disk model to the data, we measure a black hole mass of $2.42^{+0.22}_{-0.20} \times 10^7$ M$_{\odot}$ and a geometric distance to the system of $87.6^{+7.9}_{-7.2}$ Mpc. Assuming a CMB-frame recession velocity of $7308 \pm 150$ km s$^{-1}$, we constrain the Hubble constant to $H_0 = 81.0^{+7.4}_{-6.9}$ (stat.) $\pm 1.4$ (sys.) km s$^{-1}$ Mpc$^{-1}$.

We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

We present VLBI images and kinematics of water maser emission in three active galaxies: NGC 5728, Mrk 1, and IRAS 08452-0011. IRAS 08452-0011 is a triple-peaked H2O megamaser, consistent with a Keplerian rotating disk, indicating a black hole mass of (3.2+/-0.2)x10^7 M_sun. NGC 5728 and Mrk 1 display double-peaked spectra and VLBI imaging reveal complicated gas kinematics, which do not allow for a robust determination of black hole mass. The two systems are either gas disks perturbed by AGN winds or part of outflows. We find that disturbed morphology and kinematics are a ubiquitous feature of all double-peaked maser systems, implying that these maser sources may reside in environments where AGN winds are prominent at ~1 pc scale and have significant impact on the masing gas. Such AGN tend to have black hole mass M_BH < 8x10^6 M_sun and Eddington ratios lambda_Edd >~ 0.1, while the triple-peaked megamasers show an opposite trend.

We present a sample of nearby dwarf galaxies with radio-selected accreting massive black holes (BHs), the majority of which are non-nuclear. We observed 111 galaxies using sensitive, high-resolution observations from the Karl G. Jansky Very Large Array (VLA) in its most extended A-configuration at X-band (~8-12 GHz), yielding a typical angular resolution of ~0.25" and rms noise of ~15 uJy. Our targets were selected by cross matching galaxies with stellar masses M_stellar < 3 x 10^9 M_sun and redshifts z<0.055 in the NASA-Sloan Atlas with the VLA Faint Images of the Radio Sky at Twenty centimeters (FIRST) Survey. With our new high-resolution VLA observations, we detect compact radio sources towards 39 galaxies and carefully evaluate possible origins for the radio emission including thermal HII regions, supernova remnants, younger radio supernovae, background interlopers, and AGNs in the target galaxies. We find that 13 dwarf galaxies almost certainly host active massive BHs despite the fact that only one object was previously identified as having optical signatures of an AGN. We also identify a candidate dual radio AGN in a more massive galaxy system. The majority of the radio-detected BHs are offset from the center of the host galaxies with some systems showing signs of interactions/mergers. Our results indicate that massive BHs need not always live in the nuclei of dwarf galaxies, confirming predictions from simulations. Moreover, searches attempting to constrain BH seed formation using observations of dwarf galaxies need to account for such a population of "wandering" BHs.

Nearby blue compact dwarf (BCD) galaxies are arguably our best local analogues of galaxies in the earlier Universe that may host relics of black hole (BH) seeds. Here we present high-resolution Chandra X-ray Observatory and Karl G. Jansky Very Large Array (VLA) observations of five nearby BCDs with stellar masses of less than the Small Magellanic Cloud ($M_\star \sim 10^{7} - 10^{8.4}$ $M_\odot$). We search for signatures of accreting massive BHs at X-ray and radio wavelengths, which are more sensitive to lower BH accretion rates than optical searches. We detect a total of 10 hard X-ray sources and 10 compact radio sources at luminosities consistent with star-formation-related emission. We find one case of a spatially-coincident X-ray and radio source within the astrometric uncertainties. If the X-ray and radio emission are indeed coming from the same source, the origin of the radiation is plausibly from an active massive BH with log $(M_{\rm BH}/M_{\odot}) \sim 4.8 \pm 1.1$. However, given that the X-ray and radio emission are also coincident with a young star cluster complex, we consider the combination of an X-ray binary and a supernova remnant (or \HII\ region) a viable alternative explanation. Overall, we do not find compelling evidence for active massive BHs in our target BCDs, which on average have stellar masses more than an order of magnitude lower than previous samples of dwarf galaxies found to host massive BHs. Our results suggest that moderately accreting massive BHs in BCDs are not so common as to permit unambiguous detection in a small sample.

The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a universal backend and (d) augment the VLBI facility by integrating the 12m telescope for phase referencing. These upgrades to the Arecibo telescope are critical to keep the national facility in the forefront of research in radio astronomy while maintaining its dominance in radar studies of near-Earth asteroids, planets and satellites. In the next decade, the Arecibo telescope will play a synergistic role with the upcoming facilities such as ngVLA, SKA and the now commissioned FAST telescope. Further, the observatory will be actively engaged in mentoring and training programs for students from a diverse background.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

An essential part of the paradigm describing active galactic nuclei is the alignment between the radio jet and the associated rotation axis of the sub-pc sized accretion disks. Because of the small linear and angular scales involved, this alignment has not yet been checked in a sufficient number of Low Luminosity Active Galactic Nuclei (LLAGNs). The project intends to examine the validity of this paradigm by measuring the radio continuum on the same physical scale as the accretion disks, to investigate any possible connection between these disks and the radio continuum. We observed a sample of 18 LLAGNs in the 4.8 GHz (6 cm) radio continuum using the Very Long Baseline Array (VLBA) with 3.3 to 6.5 milliarcseconds resolution. The sources were selected to show both an edge-on accretion disk revealed by 22 GHz H$_{2}$O megamaser emission and signatures of a radio jet. Furthermore, the sources were previously detected in 33GHz radio continuum observations made with the Very Large Array. Five out of 18 galaxies observed were detected at 8 sigma or higher levels (Mrk0001, Mrk1210, Mrk1419, NGC2273 and UGC3193). While all these sources are known to have maser disks, four of these five sources exhibit a maser disk with known orientation. For all four sources, the radio continuum is misaligned relative to the rotation axis of the maser disk, but with a 99.1% confidence level, the orientations are not random and are confined to a cone within 32 degree of the maser disk's normal. Among the four sources the misalignment of the radio continuum with respect to the normal vector to the maser disk is smaller when the inner radius of the maser disk is larger. Furthermore, a correlation is observed between the 5 GHz VLBA radio continuum and the [OIII] luminosity and also with the H$_{2}$O maser disk's inner radius.

We identified 15,658 NVSS radio sources among the 55,288 2MASX galaxies brighter than $k_\mathrm{20fe} = 12.25$ at $\lambda = 2.16\,\mu\mathrm{m}$ and covering the $\Omega =7.016$ sr of sky defined by J2000 $\delta > -40^\circ$ and $\vert b \vert > 20^\circ$. The complete sample of 15,043 galaxies with 1.4 GHz flux densities $S \geq 2.45 \mathrm{~mJy}$ contains a 99.9% spectroscopically complete subsample of 9,517 galaxies with $k_\mathrm{20fe} \leq 11.75$. We used only radio and infrared data to quantitatively distinguish radio sources powered primarily by recent star formation from those powered by active galactic nuclei. The radio sources with $\log[L(\mathrm{W~Hz}^{-1})] > 19.3$ that we used to derive the local spectral luminosity and power-density functions account for $>99$% of the total 1.4~GHz spectral power densities $U_\mathrm{SF} = (1.54 \pm 0.20) \times 10^{19} \mathrm{~W~Hz}^{-1} \mathrm{~Mpc}^{-3}$ and $U_\mathrm{AGN} = (4.23 \pm 0.78) \times 10^{19} \mathrm{~W~Hz}^{-1} \mathrm{~Mpc}^{-3}$ in the universe today, and the spectroscopic subsample is large enough that the quoted errors are dominated cosmic variance. The recent comoving star-formation rate density indicated by $U_\mathrm{SF}$ is $\psi \approx 0.015~ M_\odot \mathrm{~yr}^{-1} \mathrm{~Mpc}^{-3}$.

We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2\sigma$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3\sigma$ ($\rm5\sigma$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN.

Deep continuum surveys at radio wavelengths are able to cover large areas, yield high angular resolution, and do not suffer from dust extinction, thus providing a robust way to measure the star formation history of the universe. However, with the current sensitivities of existing radio telescopes, it remains challenging to detect galaxies that dominate the cosmic star formation history even with extremely long integrations. With the ngVLA, a new portion of parameter space will be opened up for radio continuum surveys: deep (~200 nJy/bm), large-area (~1 deg^2), sub-arcsecond surveys at high frequencies (~8 GHz), where the observed radio emission from high-z galaxies should be dominated by free-free emission, providing a robust measurement for the star formation history of the universe. By being able to image the star formation activity, unbiased by dust, for a large, homogeneous sample of galaxies with a wide range of luminosities into the epoch of reionization, such surveys with the ngVLA will be highly complementary to those conducted by JWST, which will only be sensitive to un-obscured star formation.

In combination with observations of the Cosmic Microwave Background, a measurement of the Hubble Constant provides a direct test of the standard $\Lambda$CDM cosmological model and a powerful constraint on the equation of state of dark energy. Observations of circumnuclear water vapor megamasers at 22 GHz in nearby active galaxies can be used to measure their distances, geometrically, and thereby provide a direct, one step measurement of the Hubble Constant. The measurement is independent of distance ladders and standard candles. With present-day instrumentation, the Megamaser Cosmology Project is expected to reach a $\sim$4% measurement using the megamaser technique, based on distances to fewer than 10 megamaser systems. A goal of the observational cosmology community is to reach a percent-level measurement of H$_0$ from several independent astrophysical measurements to minimize the systematics. The ngVLA will provide the sensitivity at 22 GHz required to discover and measure the additional megamaser disks that will enable an H$_0$ measurement at the $\sim$1% level.

Radio continuum observations have proven to be a workhorse in our understanding of the star formation process (i.e., stellar birth and death) from galaxies both in the nearby universe and out to the highest redshifts. In this article we focus on how the ngVLA will transform our understanding of star formation by enabling one to map and decompose the radio continuum emission from large, heterogeneous samples of nearby galaxies on $\gtrsim 10$\,pc scales to conduct a proper accounting of the energetic processes powering it. At the discussed sensitivity and angular resolution, the ngVLA will simultaneously be able to create maps of current star formation activity at $\sim$100\,pc scales, as well as detect and characterize (e.g., size, spectral shape, density, etc.) discrete H{\sc ii} regions and supernova remnants on 10\,pc scales in galaxies out to the distance of the Virgo cluster. Their properties can then be used to see how they relate to the local and global ISM and star formation conditions. Such investigations are essential for understanding the astrophysics of high-$z$ measurements of galaxies, allowing for proper modeling of galaxy formation and evolution.

H$_2$O megamasers residing in the accretion disks of active galactic nuclei (AGN) exhibit Keplerian rotation about the central supermassive black hole (SMBH). Such disk maser systems are excellent tools for diagnosing the kinematic status of the SMBH, and they currently provide the only direct and unambiguous measure of SMBH velocities outside of the Milky Way. We have measured the galaxy recession velocities for a sample of 10 maser disk systems using a combination of spatially resolved HI disk modeling, spatially integrated HI profile fitting, and optical spectral line and continuum fitting. In comparing the SMBH velocities to those of their host galaxies, we find two (out of 10) systems -- J0437+2456 and NGC 6264 -- for which the SMBH and galaxy velocities show a statistically significant ($>$3$\sigma$) difference. For NGC 6264 the apparent velocity offset can likely be explained by ionized gas motion within the host galaxy (e.g., from AGN-driven shocks). The velocity measurements for J0437+2456, however, imply a SMBH peculiar velocity of $69.6 \pm 12.7$ km s$^{-1}$ (5.5$\sigma$). We thus consider J0437+2456 to be a promising candidate for hosting either a recoiling or binary SMBH, though additional observations are necessary to exclude the possibility of a systematic offset between the galactic recession velocity and that measured using the optical spectrum.

The homogeneous, isotropic, and flat $\Lambda$CDM universe favored by observations of the cosmic microwave background can be described using only Euclidean geometry, locally correct Newtonian mechanics, and the basic postulates of special and general relativity. We present simple derivations of the most useful equations connecting astronomical observables (redshift, flux density, angular diameter, brightness, local space density,...) with the corresponding intrinsic properties of distant sources (lookback time, distance, spectral luminosity, linear size, specific intensity, source counts,...). We also present an analytic equation for lookback time that is accurate within 0.1% for all redshifts $z$. The exact equation for comoving distance is an elliptic integral that must be evaluated numerically, but we found a simple approximation with errors $< 0.2$% for all redshifts up to $z \approx 50$.

We made two new sensitive (rms noise sigma_n ~ 1 microJy/beam) high resolution (theta = 3.0" and theta = 0.66" FWHM) S--band (2 < nu < 4 GHz) images covering a single JVLA primary beam (FWHM ~ 14') centered on J2000 RA = 10 46, Dec = 59 01 in the Lockman Hole. These images yielded a catalog of 792 radio sources, 97.7 +/- 0.8% of which have infrared counterparts stronger than S ~ 2 microJy at lambda = 4.5 micron. About 91% of the radio sources found in our previously published, comparably sensitive low resolution (theta = 8" FWHM) image covering the same area were also detected at 0.66" resolution, so most radio sources with S_3GHz >~ 5 microJy have angular structure phi <~ 0.66". The ratios of peak brightness in the 0.66" and 3" images have a distribution indicating that most microJy radio sources are quite compact, with a median Gaussian angular diameter <phi> = 0.3" +/- 0.1" FWHM and an rms scatter sigma_phi <~ 0.3" of individual sizes. Most of our microJy radio sources obey the tight far-infrared/radio correlation, indicating that they are powered by star formation. The median effective angular radius enclosing half the light emitted by an exponential disk is <rho_e> ~ <phi>/2.43 ~ 0.12", so the median effective radius of star-forming galaxies at redshifts z~1 is <r_e> ~ 1.0 kpc.

We summarize the radio synchrotron background workshop that took place July 19-21, 2017 at the University of Richmond. This first scientific meeting dedicated to the topic was convened because current measurements of the diffuse radio monopole reveal a surface brightness that is several times higher than can be straightforwardly explained by known Galactic and extragalactic sources and processes, rendering it by far the least well understood photon background at present. It was the conclusion of a majority of the participants that the radio monopole level is at or near that reported by the ARCADE 2 experiment and inferred from several absolutely calibrated zero level lower frequency radio measurements, and unanimously agreed that the production of this level of surface brightness, if confirmed, represents a major outstanding question in astrophysics. The workshop reached a consensus on the next priorities for investigations of the radio synchrotron background.

We present high resolution (sub-mas) VLBI maps of nuclear H2O megamasers for seven galaxies. In UGC6093, the well-aligned systemic masers and high-velocity masers originate in an edge-on, flat disk and we determine the mass of the central SMBH to be M_SMBH = 2.58*10^7Msun(+-7%). For J1346+5228, the distribution of masers is consistent with a disk, but the faint high-velocity masers are only marginally detected, and we constrain the mass of the SMBH to be in the range 1.5-2.0*10^7Msun. The origin of the masers in Mrk1210 is less clear, as the systemic and high-velocity masers are misaligned and show a disorganized velocity structure. We present one possible model in which the masers originate in a tilted, warped disk, but we do not rule out the possibility of other explanations including outflow masers. In NGC6926, we detect a set of redshifted masers, clustered within a pc of each other, and a single blueshifted maser about 4.4pc away, an offset that would be unusually large for a maser disk system. Nevertheless, if it is a disk system, we estimate the enclosed mass to be M_SMBH<4.8*10^7 Msun . For NGC5793, we detect redshifted masers spaced about 1.4pc from a clustered set of blueshifted features. The orientation of the structure supports a disk scenario as suggested by Hagiwara et al.(2001). We estimate the enclosed mass to be M SMBH<1.3*10^7 Msun. For NGC2824 and J0350-0127, the masers may be associated with pc or sub-pc scale jets or outflows.

We investigate the nuclear environment of galaxies with observed 22 GHz water megamaser in their subparsec edge-on accretion disks, using 33 GHz (9mm) radio continuum data from VLA, with a resolution of ~ 0.2-0.5 arcsecs, and relate the maser and host galaxy properties to those of its radio continuum emission. Eighty-seven percent (21 out of 24) galaxies in our sample show 33 GHz radio continuum emission at levels of 4.5-240 $\sigma$. Five sources show extended emission, including one source with two main components and one with three main components. The remaining detected 16 sources exhibit compact cores within the sensitivity limits. Little evidence is found for extended jets (>300 pc) in most sources. Either they do not exist, or our chosen frequency of 33 GHz is too high for a detection of these supposedly steep spectrum features. In only one source among those with known maser disk orientation, NGC4388, we found an extended jet-like feature that appears to be oriented perpendicular to the water megamaser disk. Smaller 100-300 pc sized jets might also be present, as is suggested by the beam-deconvolved morphology of our sources. Whenever possible, central positions with accuracies of 20-280 mas are provided. A correlation analysis shows that the 33 GHz luminosity weakly correlates with the infrared luminosity. The 33 GHz luminosity is anticorrelated with the circular velocity of the galaxy. The black hole masses show stronger correlations with water maser luminosity than with 1.4 GHz, 33 GHz, or hard X-ray luminosities. Furthermore, the inner radii of the disks show stronger correlations with 1.4 GHz, 33 GHz, and hard X-ray luminosities than their outer radii, suggesting that the outer radii may be affected by disk warping, star formation, or peculiar density distributions.

We present Very Large Array observations of the 33 GHz radio continuum emission from 22 local ultraluminous and luminous infrared (IR) galaxies (U/LIRGs). These observations have spatial (angular) resolutions of 30--720 pc (0.07"-0.67") in a part of the spectrum that is likely to be optically thin. This allows us to estimate the size of the energetically dominant regions. We find half-light radii from 30 pc to 1.7 kpc. The 33 GHz flux density correlates well with the IR emission, and we take these sizes as indicative of the size of the region that produces most of the energy. Combining our 33 GHz sizes with unresolved measurements, we estimate the IR luminosity and star formation rate per area, and the molecular gas surface and volume densities. These quantities span a wide range (4 dex) and include some of the highest values measured for any galaxy (e.g., $\mathrm{\Sigma_{SFR}^{33GHz} \leq 10^{4.1} M_{\odot} yr^{-1} kpc^{-2}}$). At least $13$ sources appear Compton thick ($\mathrm{N_{H}^{33GHz} \geq 10^{24} cm^{-2}}$). Consistent with previous work, contrasting these data with observations of normal disk galaxies suggests a nonlinear and likely multi-valued relation between SFR and molecular gas surface density, though this result depends on the adopted CO-to-H$_{2}$ conversion factor and the assumption that our 33 GHz sizes apply to the gas. 11 sources appear to exceed the luminosity surface density predicted for starbursts supported by radiation pressure and supernovae feedback, however we note the need for more detailed observations of the inner disk structure. U/LIRGs with higher surface brightness exhibit stronger [{\sc Cii}] 158$\mu$m deficits, consistent with the suggestion that high energy densities drive this phenomenon.

We investigate extremely luminous dusty galaxies in the environments around WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7, respectively. Previous observations have detected overdensities of companion submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with overdensities of ~ 2 - 3 and ~ 5 - 6 , respectively. We find that the space densities in both samples to be overdense compared to normal star-forming galaxies and submillimetre galaxies (SMGs) in the SCUBA-2 Cosmology Legacy Survey (S2CLS). Both samples of companion sources have consistent mid-IR colours and mid-IR to submm ratios as SMGs. The brighter population around WISE/radio AGNs could be responsible for the higher overdensity reported. We also find the star formation rate density (SFRDs) are higher than the field, but consistent with clusters of dusty galaxies. WISE-selected AGNs appear to be good signposts for protoclusters at high redshift on arcmin scales. The results reported here provide an upper limit to the strength of angular clustering using the two-point correlation function. Monte Carlo simulations show no angular correlation, which could indicate protoclusters on scales larger than the SCUBA-2 1.5arcmin scale maps.

(Abridged) Our sensitive ($\sigma_{\rm n} = 572\,{\rm nJy\,beam}^{-1}$), high-resolution (FWHM $\theta_{1/2} = 220\,{\rm mas} \approx2\mathrm{\,kpc~at~}z\gtrsim1$) 10$\,$GHz image covering a single Karl G.~Jansky Very Large Array (VLA) primary beam (FWHM $\Theta_{1/2} = 4.25'$) in the GOODS-N field contains 32 sources with $S_{\rm p}\gtrsim2\,\mu{\rm Jy~beam}^{-1}$ and optical and/or near-infrared (OIR) counterparts. Most are about as large as the star-forming regions that power them. Their median FWHM major axis is $\langle\theta_{\rm M} \rangle=167\pm32\,{\rm mas} \approx 1.2\pm0.28\,{\rm kpc}$ with rms scatter 91 was $\approx$ 0.79 kpc. In units of the effective radius $r_{\rm e}$ that encloses half their flux, these radio sizes are $\langle r_{\rm e}\rangle = 69\pm13{\rm mas} \approx \pm114\mathrm{\,pc}$ and have rms scatter $38\mathrm{\,mas}\approx324\mathrm{\,pc}$. These sizes are smaller than those measured at lower radio frequencies, but agree with dust emission sizes measured at mm/sub-mm wavelengths and extinction-corrected H$\alpha$ sizes. We made a low-resolution ($\theta_{1/2}=1.0"$) image with $\approx10\times$ better brightness sensitivity to detect extended sources and measure matched-resolution spectral indices $\alpha_{1.4}^{10}$. It contains 6 new sources with $S_{\rm p}\gtrsim3.9\,\mu{\rm Jy~beam}^{-1}$ and OIR counterparts. The median redshift of all 38 sources is $1.24\pm0.15$. The 19 sources with 1.4$\,$GHz counterparts have median spectral index $-0.74\pm0.10$ with rms scatter $0.35$. Including upper limits on $\alpha$ for sources not detected at 1.4$\,$GHz flattens the median to $\gtrsim-0.61$, suggesting that the $\mu$Jy radio sources at higher redshifts, and hence selected at higher rest-frame frequencies, may have flatter spectra. If the non-thermal spectral index is -0.85, the median thermal fraction at rest-frame frequency 20$\,$GHz is $\gtrsim$48%.

During a systematic search for supermassive black holes (SMBHs) not in galactic nuclei, we identified the compact symmetric radio source B3 1715+425 with an emission-line galaxy offset ~ 8.5 kpc from the nucleus of the brightest cluster galaxy (BCG) in the redshift $z = 0.1754$ cluster ZwCl 8193. B3 1715+425 is too bright (brightness temperature $\sim 3 \times 10^{10}$ K at observing frequency 7.6 GHz) and too luminous (1.4 GHz luminosity $\sim 10^{25}$ W/Hz) to be powered by anything but a SMBH, but its host galaxy is much smaller ($\sim 0.9$ kpc $\times$ 0.6 kpc full width between half-maximum points) and optically fainter (R-band absolute magnitude $\sim -18.2$) than any other radio galaxy. Its high radial velocity $\sim 1860$ km/s relative to the BCG, continuous ionized wake extending back to the BCG nucleus, and surrounding debris indicate that the radio galaxy was tidally shredded passing through the BCG core, leaving a nearly naked supermassive black hole fleeing from the BCG with space velocity $> 2000$ km/s. The radio galaxy has mass $< 6 \times 10^9$ solar masses and infrared luminosity $\sim 3 \times 10^{11}$ solar luminosities close to its dust Eddington limit, so it is vulnerable to further mass loss from radiative feedback.

As part of the Megamaser Cosmology Project (MCP), we present VLBI maps of nuclear water masers toward five galaxies. The masers originate in sub-parsec circumnuclear disks. For three of the galaxies, we fit Keplerian rotation curves to estimate their supermassive black hole (SMBH) masses, and determine (2.9 $\pm$ 0.3) $\times~10^{6}M_\odot$ for J0437+2456, (1.7 $\pm$ 0.1) $\times~10^{7}M_\odot$ for ESO 558$-$G009, and (1.1 $\pm$ 0.2) $\times~10^{7}M_\odot$ for NGC 5495. In the other two galaxies, Mrk 1029 and NGC 1320, the geometry and dynamics are more complicated and preclude robust black hole mass estimates. Including our new results, we compiled a list of 15 VLBI-confirmed disk maser galaxies with robust SMBH mass measurements. With this sample, we confirm the empirical relation of $R_{out} \propto 0.3 M_{SMBH}$ reported in Wardle & Yusef-Zadeh (2012). We also find a tentative correlation between maser disk outer radii and WISE luminosity. We find no correlations of maser disk size with X-ray 2-10 keV luminosity or [O III] luminosity.

We use new precision measurements of black hole masses from water megamaser disks to investigate scaling relations between macroscopic galaxy properties and supermassive black hole (BH) mass. The megamaser-derived BH masses span 10^6-10^8 M_sun, while all the galaxy properties that we examine (including stellar mass, central mass density, central velocity dispersion) lie within a narrow range. Thus, no galaxy property correlates tightly with M_BH in ~L* spiral galaxies. Of them all, stellar velocity dispersion provides the tightest relation, but at fixed sigma* the mean megamaser M_BH are offset by -0.6+/-0.1 dex relative to early-type galaxies. Spiral galaxies with non-maser dynamical BH masses do not show this offset. At low mass, we do not yet know the full distribution of BH mass at fixed galaxy property; the non-maser dynamical measurements may miss the low-mass end of the BH distribution due to inability to resolve the spheres of influence and/or megamasers may preferentially occur in lower-mass BHs.

We discuss 6 GHz JVLA observations covering a volume-limited sample of 178 low redshift ($0.2 < z < 0.3$) optically selected QSOs. Our 176 radio detections fall into two clear categories: (1) About $20$\% are radio-loud QSOs (RLQs) having spectral luminosities $L_6 \gtrsim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ primarily generated in the active galactic nucleus (AGN) responsible for the excess optical luminosity that defines a \emph{bona fide} QSO. (2) The radio-quiet QSOs (RQQs) have $10^{\,21} \lesssim L_6 \lesssim 10^{\,23.2} \mathrm{~W~Hz}^{-1}$ and radio sizes $\lesssim 10 \mathrm{~kpc}$, and we suggest that the bulk of their radio emission is powered by star formation in their host galaxies. "Radio silent" QSOs ($L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$) are rare, so most RQQ host galaxies form stars faster than the Milky Way; they are not "red and dead" ellipticals. Earlier radio observations did not have the luminosity sensitivity $L_6 \lesssim 10^{\,21} \mathrm{~W~Hz}^{-1}$ needed to distinguish between such RLQs and RQQs. Strong, generally double-sided, radio emission spanning $\gg 10 \mathrm{~kpc}$ was found associated with 13 of the 18 RLQ cores having peak flux densities $S_\mathrm{p} > 5 \mathrm{~mJy~beam}^{-1}$ ($log(L) \gtrsim 24$). The radio luminosity function of optically selected QSOs and the extended radio emission associated with RLQs are both inconsistent with simple "unified" models that invoke relativistic beaming from randomly oriented QSOs to explain the difference between RLQs and RQQs. Some intrinsic property of the AGNs or their host galaxies must also determine whether or not a QSO appears radio loud.

This is the second of two papers describing the observations and source catalogues derived from sensitive 3-GHz images of the Lockman Hole North using the Karl G. Jansky Very Large Array. We describe the reduction and cataloguing process, which yielded an image with 8 arcsec resolution and instrumental noise of \sigma$_{\rm n}$=1.0 \mu Jy beam$^{-1}$ rms (before primary beam corrections) and a catalogue of 558 sources detected above 5\sigma$_{\rm n}$. We include details of how we estimate source spectral indices across the 2-GHz VLA bandwidth, finding a median index of -0.76$\pm$0.04. Stacking of source spectra reveals a flattening of spectral index with decreasing flux density. We present a source count derived from the catalogue. We show a traditional count estimate compared with a completely independent estimate made via a P(D) confusion analysis, and find very good agreement. Cross-matches of the catalogue with X-ray, optical, infrared, radio, and redshift catalogues are also presented. The X-ray, optical and infrared data, as well as AGN selection criteria allow us to classify 10 per cent as radio-loud AGN, 28 percent as radio-quiet AGN, and 58 per cent as star-forming galaxies, with only 4 per cent unclassified.

This is the first of two papers describing the observations and cataloguing of deep 3-GHz observations of the Lockman Hole North using the Karl G. Jansky Very Large Array. The aim of this paper is to investigate, through the use of simulated images, the uncertainties and accuracy of source-finding routines, as well as to quantify systematic effects due to resolution, such as source confusion and source size. While these effects are not new, this work is intended as a particular case study that can be scaled and translated to other surveys. We use the simulations to derive uncertainties in the fitted parameters, as well as bias corrections for the actual catalogue (presented in Paper 2). We compare two different source-finding routines, OBIT and AEGEAN, and two different effective resolutions, 8 and 2.75 arcsec. We find that the two routines perform comparably well, with OBIT being slightly better at de-blending sources, but slightly worse at fitting resolved sources. We show that 30 to 70 per cent of sources are missed or fit inaccurately once the source size becomes larger than the beam, possibly explaining source count errors in high-resolution surveys. We also investigate the effect of blending, finding that any sources with separations smaller than the beam size are fit as single sources. We show that the use of machine-learning techniques can correctly identify blended sources up to 90 per cent of the time, and prior-driven fitting can lead to a 70 per cent improvement in the number of de-blended sources.

As part of the Megamaser Cosmology Project (MCP), here we present a new geometric distance measurement to the megamaser galaxy NGC 5765b. Through a series of VLBI observations, we have confirmed the water masers trace a thin, sub-parsec Keplerian disk around the nucleus, implying an enclosed mass of 4.55 $\pm$ 0.40 $\times~10^{7}M_\odot$. Meanwhile, from single dish monitoring of the maser spectra over two years, we measured the secular drifts of maser features near the systemic velocity of the galaxy with rates between 0.5 and 1.2 km s$^{-1}$ yr$^{-1}$. Fitting a warped, thin disk model to these measurements, we determine a Hubble Constant $H_{0}$ of 66.0 $\pm$ 6.0 km s$^{-1}$ Mpc$^{-1}$ with the angular-diameter distance to NGC 5765b of 126.3 $\pm$ 11.6 Mpc. Apart from the distance measurement, we also investigate some physical properties related to the maser disk in NGC 5765b. The high-velocity features are spatially distributed into several clumps, which may indicate the existence of a spiral density wave associated with the accretion disk. For the red-shifted features, the envelope defined by the peak maser intensities increases with radius. The profile of the systemic masers in NGC 5765b is smooth and shows almost no structural changes over the two years of monitoring time, which differs from the more variable case of NGC 4258.

We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 microns, and we find that the sample has large mid- to far-infrared luminosity ratios consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < log P3.0 GHz (W/Hz) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log M(BH) (Msun) < 10.2. The rest frame 1-5 um SEDs are very similar to the "Hot DOGs" (Hot Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA detected sources are 9.9 < log M(ISM) (Msun) < 11.75 assuming a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) reaching several thousand Msun/yr, depending on the assumed dust temperature, however we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with ~ equal contributions from an obscured AGN and a young (10-15 Myr) compact starburst.

We present Karl G. Jansky Very Large Array (VLA) Ka band (33 GHz) and Atacama Large Millimeter Array (ALMA) Band 3 (94.5 GHz) continuum images covering the nucleus and two extranuclear star-forming regions within the nearby galaxy NGC 3627 (M 66), observed as part of the Star Formation in Radio Survey (SFRS). Both images achieve an angular resolution of $\lesssim$2\arcsec, allowing us to map the radio spectral indices and estimate thermal radio fractions at a linear resolution of $\lesssim$90 pc at the distance of NGC 3627. The thermal fraction at 33 GHz reaches unity at and around the peaks of each HII region; we additionally observed the spectral index between 33 and 94.5 GHz to become both increasingly negative and positive away from the peaks of the HII regions, indicating an increase of non-thermal extended emission from diffusing cosmic-ray electrons and the possible presence of cold dust, respectively. While the ALMA observations were optimized for collecting continuum data, they also detected line emission from the $J=1\rightarrow0$ transitions of HCN and HCO$^{+}$. The peaks of dense molecular gas traced by these two spectral lines are spatially offset from the peaks of the 33 and 94.5 GHz continuum emission for the case of the extranuclear star-forming regions, indicating that our data reach an angular resolution at which one can spatially distinguish sites of recent star formation from the sites of future star formation. Finally, we find trends of decreasing dense gas fraction and velocity dispersion with increasing star formation efficiency among the three regions observed, indicating that the dynamical state of the dense gas, rather than its abundance, plays a more significant role in the star formation process.

We use single-dish radio spectra of known 22 GHz H$_2$O megamasers, primarily gathered from the large dataset observed by the Megamaser Cosmology Project, to identify Keplerian accretion disks and to investigate several aspects of the disk physics. We test a mechanism for maser excitation proposed by Maoz & McKee (1998), whereby population inversion arises in gas behind spiral shocks traveling through the disk. Though the flux of redshifted features is larger on average than that of blueshifted features, in support of the model, the high-velocity features show none of the predicted systematic velocity drifts. We find rapid intra-day variability in the maser spectrum of ESO 558-G009 that is likely the result of interstellar scintillation, for which we favor a nearby ($D \approx 70$ pc) scattering screen. In a search for reverberation in six well-sampled sources, we find that any radially-propagating signal must be contributing $\lesssim$10% of the total variability. We also set limits on the magnetic field strengths in seven sources, using strong flaring events to check for the presence of Zeeman splitting. These limits are typically 200--300 mG ($1\sigma$), but our most stringent limits reach down to 73 mG for the galaxy NGC 1194.

We present JCMT SCUBA-2 850microns submillimetre (submm) observations of 30 mid-infrared (mid-IR) luminous AGN, detected jointly by the WISE all-sky IR survey and the NVSS/FIRST radio survey. These rare sources are selected by their extremely red mid-infrared spectral energy distributions (SEDs) and compact radio counterparts. Further investigations show that they are highly obscured, have abundant warm AGN-heated dust and are thought to be experiencing intense AGN feedback. These galaxies appear to be consistent with an AGN-dominated galaxy, and could be a transient phase of merging galaxies. When comparing the number of submm galaxies (SMGs) detected serendipitously in the surrounding 1.5-arcmin to those in blank-field submm surveys, there is a very significant overdensity, of order 5, but no sign of radial clustering centred at our primary objects. The WISE/radio-selected AGN thus reside in 10-Mpc-scale overdense environments, that could be forming in pre-viralised clusters of galaxies. WISE/radio-selected AGNs appear to be the strongest signposts of high-density regions of active, luminous and dusty galaxies. SCUBA-2 850microns observations indicate that their submm fluxes are low compared to many popular AGN SED templates, hence the WISE/radio-selected AGNs have either less cold and/or more warm dust emission than normally assumed for typical AGN. Most of the targets are not detected, only four targets are detected at SCUBA-2 850microns, and have total IR luminosities >= 10^13 L_solar, if their redshifts are consistent with the subset of the 10 SCUBA-2 undetected targets with known redshifts, z ~ 0.44 - 2.86.

The proposed VLA Sky Survey (VLASS) comprises two distinct S-band (2--4 GHz) surveys: (1) All-Sky covering the sky north of -40 deg with rms noise 69 microJy/beam = 1.5 K in the 2.5 arcsec beam and (2) Deep covering 10 square deg with rms noise 1.5 microJy/beam = 0.32 K in the 0.8 arcsec beam. This review compares the scientific goals and technical capabilities of the VLASS proposal (2015 Jan 15 version), using new calculations of performance metrics for surveys made with large fractional bandwidths.

We present observations of the H2O megamasers in the accretion disk of NGC 6323. By combining interferometric and spectral monitoring data, we estimate H$_{0} = 73^{+26}_{-22}$ km/s/Mpc, where the low strength of the systemic masers (<15 mJy) limits the accuracy of this estimate. The methods developed here for dealing with weak maser emission provide guidance for observations of similar sources, until significant increases in radio telescope sensitivity, such as anticipated from the next generation Very Large Array, are realized.

Astronomy at or below the 'survey threshold' has expanded significantly since the publication of the original 'Science with the Square Kilometer Array' in 1999 and its update in 2004. The techniques in this regime may be broadly (but far from exclusively) defined as 'confusion' or 'P(D)' analyses (analyses of one-point statistics), and 'stacking', accounting for the flux-density distribution of noise-limited images co-added at the positions of objects detected/isolated in a different waveband. Here we discuss the relevant issues, present some examples of recent analyses, and consider some of the consequences for the design and use of surveys with the SKA and its pathfinders.

We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emission surrounding both nuclei and conclude that is mostly synchrotron in nature. The spatial distributions of radio emission in both nuclei are well described by exponential profiles. These have deconvolved half-light radii of 51 and 35 pc for the eastern and western nuclei, and they match the number density profile of radio supernovae observed with very long baseline interferometry. This similarity might be due to the fast cooling of cosmic rays electrons caused by the presence of a strong (~ mG) magnetic field in this system. We estimate high luminosity surface densities of $\mathrm{\Sigma_{IR} \sim 4.2^{+1.6}_{-0.7} \times 10^{13}}$ (east) and $\mathrm{\sim 9.7^{+3.7}_{-2.4} \times 10^{13}~(west)~L_{\odot}~kpc^{-2}}$, and star formation rate surface densities of $\mathrm{\Sigma_{SFR} \sim 10^{3.7\pm0.1}}$ (east) and $\mathrm{\sim 10^{4.1\pm0.1}~(west)~M_{\odot}~yr^{-1}~kpc^{-2}}$. These values, especially for the western nucleus are, to our knowledge, the highest luminosity and star formation rate surface densities measured for any star-forming system. Despite these high values, the nuclei lie below the dusty Eddington limit in which radiation pressure is balanced only by self-gravity. The small measured sizes also imply that the nuclei of Arp 220 are only transparent in the frequency range ~ 5 to 350 GHz. Our results offer no clear evidence that an active galactic nucleus dominates the emission from either nucleus at 33 GHz.

The new broad-band capabilities of large radio interferometers such as the GMRT and JVLA allow for long-integration mosaic imaging observations to create ultra-deep full-polarization images of the sky over wide frequency ranges. Achieving rms sensitivities of order 1 $\mu$Jy, these observations explore the radio source population at flux densities well below the regime dominated by classical radio galaxies and Active Galactic Nuclei. We present initial results from radio sources revealed with deep mosaicking observations with the GMRT and JVLA at respectively 0.6 and 5 GHz, and evidence that the $\mu$Jy sensitivity level marks the transition to detection of polarized emission from a population of sources dominated by emission from magnetic fields in the disks of starburst and normal galaxies.

The incidence and properties of present-day dwarf galaxies hosting massive black holes (BHs) can provide important constraints on the origin of high-redshift BH seeds. Here we present high-resolution X-ray and radio observations of the low-metallicity, star-forming, dwarf-galaxy system Mrk 709 with the Chandra X-ray Observatory and the Karl G. Jansky Very Large Array. These data reveal spatially coincident hard X-ray and radio point sources with luminosities suggesting the presence of an accreting massive BH (M_BH ~ 10^5-10^7 Msun). Based on imaging from the Sloan Digital Sky Survey (SDSS), we find that Mrk 709 consists of a pair of compact dwarf galaxies that appear to be interacting with one another. The position of the candidate massive BH is consistent with the optical center of the southern galaxy (Mrk 709 S), while no evidence for an active BH is seen in the northern galaxy (Mrk 709 N). We derive stellar masses of M ~ 2.5 x 10^9 Msun and M ~ 1.1 x 10^9 Msun for Mrk 709 S and Mrk 709 N, respectively, and present an analysis of the SDSS spectrum of the BH-host Mrk 709 S. At a metallicity of just ~10% solar, Mrk 709 is among the most metal-poor galaxies with evidence for an active galactic nucleus. Moreover, this discovery adds to the growing body of evidence that massive BHs can form in dwarf galaxies and that deep, high-resolution X-ray and radio observations are ideally suited to reveal accreting massive BHs hidden at optical wavelengths.

Radio source counts constrain galaxy populations and evolution, as well as the global star formation history. However, there is considerable disagreement among the published 1.4-GHz source counts below 100 microJy. Here we present a statistical method for estimating the microJy and even sub-microJy source count using new deep wide-band 3-GHz data in the Lockman Hole from the Karl G. Jansky Very Large Array (VLA). We analyzed the confusion amplitude distribution P(D), which provides a fresh approach in the form of a more robust model, with a comprehensive error analysis. We tested this method on a large-scale simulation, incorporating clustering and finite source sizes. We discuss in detail our statistical methods for fitting using Monte Carlo Markov chains, handling correlations, and systematic errors from the use of wide-band radio interferometric data. We demonstrated that the source count can be constrained down to 50 nJy, a factor of 20 below the rms confusion. We found the differential source count near 10 microJy to have a slope of -1.7, decreasing to about -1.4 at fainter flux densities. At 3GHz the rms confusion in an 8arcsec FWHM beam is ~ 1.2 microJy/beam, and a radio background temperature ~ 14mK. Our counts are broadly consistent with published evolutionary models. With these results we were also able to constrain the peak of the Euclidean normalized differential source count of any possible new radio populations that would contribute to the cosmic radio background down to 50 nJy.

A community meeting on the topic of "Radio Astronomy in the LSST Era" was hosted by the National Radio Astronomy Observatory in Charlottesville, VA (2013 May 6--8). The focus of the workshop was on time domain radio astronomy and sky surveys. For the time domain, the extent to which radio and visible wavelength observations are required to understand several classes of transients was stressed, but there are also classes of radio transients for which no visible wavelength counterpart is yet known, providing an opportunity for discovery. From the LSST perspective, the LSST is expected to generate as many as 1 million alerts nightly, which will require even more selective specification and identification of the classes and characteristics of transients that can warrant follow up, at radio or any wavelength. The LSST will also conduct a deep survey of the sky, producing a catalog expected to contain over 38 billion objects in it. Deep radio wavelength sky surveys will also be conducted on a comparable time scale, and radio and visible wavelength observations are part of the multi-wavelength approach needed to classify and understand these objects. Radio wavelengths are valuable because they are unaffected by dust obscuration and, for galaxies, contain contributions both from star formation and from active galactic nuclei. The workshop touched on several other topics, on which there was consensus including the placement of other LSST "Deep Drilling Fields," inter-operability of software tools, and the challenge of filtering and exploiting the LSST data stream. There were also topics for which there was insufficient time for full discussion or for which no consensus was reached, which included the procedures for following up on LSST observations and the nature for future support of researchers desiring to use LSST data products.

The shape of the curves defined by the counts of radio sources per unit area as a function of their flux density was one of the earliest cosmological probes. Radio source counts continue to be an area of interest, used to study the relative populations of galaxy types in the Universe (as well as investigate any cosmological evolution in luminosity functions). They are a vital consideration for determining how source confusion may limit the depth of a radio interferometer observation, and are essential for characterising extragalactic foregrounds in CMB experiments. There is currently no consensus as to the relative populations of the faintest (sub-mJy) source types, where the counts turn-up. Most of the source counts in this regime are gathered from multiple observations that each use a deep, single pointing with a radio interferometer. These independent measurements show large amounts of scatter (factors ~ a few) that significantly exceeds their stated uncertainties. In this article we use a simulation of the extragalactic radio continuum emission to assess the level at which sample variance may be the cause of the scatter. We find that the scatter induced by sample variance in the simulated counts decreases towards lower flux density bins as the raw source counts increase. The field-to-field variations are significant, and could even be the sole cause at >100 {\mu}Jy. We present a method for evaluating the flux density limit that a survey must reach in order to reduce the count uncertainty induced by sample variance to a specific value. We also derive a method for correcting Poisson errors on counts in order to include the uncertainties due to the cosmological clustering of sources. An empirical constraint on the effect of sample variance at these low luminosities is unlikely to arise until the completion of new large-scale surveys with next-generation radio telescopes.

We investigate the relationship between 8.44\,GHz brightness temperatures and 1.4 to 8.44\,GHz radio spectral indices with 6.2\,$\mu$m polycyclic aromatic hydrocarbon (PAH) emission and 9.7\,$\mu$m silicate absorption features for a sample of 36 local luminous and ultra-luminous infrared galaxies. We find that galaxies having small 6.2\,$\mu$m PAH equivalent widths (EQWs), which signal the presence of weak PAH emission and/or an excess of very hot dust, also have flat spectral indices. The three active galactic nuclei (AGN) identified through their excessively large 8.44\,GHz brightness temperatures are also identified as AGN via their small 6.2\,$\mu$m PAH EQWs. We also find that the flattening of the radio spectrum increases with increasing silicate optical depth, 8.44\,GHz brightness temperature, and decreasing size of the radio source even after removing potential AGN, supporting the idea that compact starbursts show spectral flattening as the result of increased free-free absorption. These correlations additionally suggest that the dust obscuration in these galaxies must largely be coming from the vicinity of the compact starburst itself, and is not distributed throughout the (foreground) disk of the galaxy. Finally, we investigate the location of these infrared-bright systems relative to the main sequence (star formation rate vs. stellar mass) of star-forming galaxies in the local universe. We find that the radio spectral indices of galaxies flattens with increasing distance above the main sequence, or in other words, with increasing specific star formation rate. This indicates that galaxies located above the main sequence, having high specific star formation rates, are typically compact starbursts hosting deeply embedded star formation that becomes more optically thick in the radio and infrared with increased distance above the main sequence.

We used the 1.4 GHz NVSS to study radio sources in two color-selected QSO samples: a volume-limited sample of 1313 QSOs defined by M_i < -23 in the redshift range 0.2 < z < 0.45 and a magnitude-limited sample of 2471 QSOs with m_r < 18.5 and 1.8 < z < 2.5. About 10% were detected above the 2.4 mJy NVSS catalog limit and are powered primarily by AGNs. The space density of the low-redshift QSOs evolves as rho proportional to (1+z)^6. In both redshift ranges the flux-density distributions and luminosity functions of QSOs stronger than 2.4 mJy are power laws, with no features to suggest more than one kind of radio source. Extrapolating the power laws to lower luminosities predicts the remaining QSOs should be extremely radio quiet, but they are not. Most were detected statistically on the NVSS images with median peak flux densities S_p(mJy/beam) ~ 0.3 and 0.05 in the low- and high-redshift samples, corresponding to 1.4 GHz spectral luminosities log[L(W/Hz)] ~ 22.7$ and 24.1, respectively. We suggest that the faint radio sources are powered by star formation at rates ~20 M_sun per year in the moderate luminosity (median M_i ~ -23.4) low-redshift QSOs and ~500 M_sun per year in the very luminous (M_i} ~ -27.5) high-redshift QSOs. Such luminous starbursts [ log(L / L_sun) ~ 11.2 and 12.6, respectively] are consistent with "quasar mode" accretion in which cold gas flows fuel both AGN and starburst.

In Papers I and II from the Megamaser Cosmology Project (MCP), we reported initial observations of water masers in an accretion disk of a supermassive black hole at the center of the galaxy UGC 3789, which gave an angular-diameter distance to the galaxy and an estimate of Ho with 16% uncertainty. We have since conducted more VLBI observations of the spatial-velocity structure of these water masers, as well as continued monitoring of its spectrum to better measure maser accelerations. These more extensive observations, combined with improved modeling of the masers in the accretion disk of the central supermassive black hole, confirm our previous results, but with signifcantly improved accuracy. We find Ho = 68.9 +/- 7.1 km/s/Mpc; this estimate of Ho is independent of other methods and is accurate to +/-10%, including sources of systematic error. This places UGC 3789 at a distance of 49.6 +/- 5.1 Mpc, with a central supermassive black hole of (1.16 +/- 0.12) x 10^7 Msun.

We present the direct measurement of the Hubble constant, yielding the direct measurement of the angular-diameter distance to NGC 6264 using the H$_{2}$O megamaser technique. Our measurement is based on sensitive observations of the circumnuclear megamaser disk from four observations with the Very Long Baseline Array, the Green Bank Telescope and the Effelsberg Telescope. We also monitored the maser spectral profile for 2.3 years using the Green Bank Telescope to measure accelerations of maser lines by tracking their line-of-sight velocities as they change with time. The measured accelerations suggest that the systemic maser spots have a significantly wider radial distribution than in the archetypal megamaser in NGC 4258. We model the maser emission as arising from a circumnuclear disk with orbits dominated by the central black hole. The best fit of the data gives a Hubble constant of $H_{0} =$ 68$\pm$9 km s$^{-1}$ Mpc$^{-1}$, which corresponds to an angular-diameter distance of 144$\pm$19 Mpc. In addition, the fit also gives a mass of the central black hole of (3.09$\pm$0.42)$\times10^{7}$ $M_{\odot}$. The result demonstrates the feasibility of measuring distances to galaxies located well into the Hubble flow by using circumnuclear megamaser disks.

We present 33\,GHz photometry of 103 galaxy nuclei and extranuclear star-forming complexes taken with the Green Bank Telescope (GBT) as part of the Star Formation in Radio Survey (SFRS). Among the sources without evidence for an AGN, and also having lower frequency radio data, we find a median thermal fraction at 33GHz of ~76% with a dispersion of ~24%. For all sources resolved on scales <0.5kpc, the thermal fraction is even larger, being >90%. This suggests that the rest-frame 33GHz emission provides a sensitive measure of the ionizing photon rate from young star-forming regions, thus making it a robust star formation rate indicator. Taking the 33GHz star formation rates as a reference, we investigate other empirical calibrations relying on different combinations of warm 24\mu m dust, total infrared (IR; 8-1000\mu m), H\alpha\ line, and far-UV continuum emission. The recipes derived here generally agree with others found in the literature, albeit with a large dispersion that most likely stems from a combination of effects. Comparing the 33GHz to total IR flux ratios as a function of the radio spectral index, measured between 1.7 and 33GHz, we find that the ratio increases as the radio spectral index flattens which does not appear to be a distance effect. Consequently, the ratio of non-thermal to total IR emission appears relatively constant, suggesting only moderate variations in the cosmic-ray electron injection spectrum and ratio of synchrotron to total cooling processes among star-forming complexes. Assuming that this trend solely arises from an increase in the thermal fraction sets a maximum on the scatter of the non-thermal spectral indices among the star-forming regions of \sigma_\alpha^{NT} < 0.13.

We report the discovery by the Wide-field Infrared Survey Explorer of the z = 2.452 source WISE J181417.29+341224.9, the first hyperluminous source found in the WISE survey. WISE 1814+3412 is also the prototype for an all-sky sample of ~1000 extremely luminous "W1W2-dropouts" (sources faint or undetected by WISE at 3.4 and 4.6 microns and well detected at 12 or 22 microns). The WISE data and a 350 micron detection give a minimum bolometric luminosity of 3.7 x 10^13 Lsun, with ~10^14 Lsun plausible. Followup images reveal four nearby sources: a QSO and two Lyman Break Galaxies (LBGs) at z = 2.45, and an M dwarf star. The brighter LBG dominates the bolometric emission. Gravitational lensing is unlikely given the source locations and their different spectra and colors. The dominant LBG spectrum indicates a star formation rate ~300 Msun/yr, accounting for < 10% of the bolometric luminosity. Strong 22 micron emission relative to 350 microns implies that warm dust contributes significantly to the luminosity, while cooler dust normally associated with starbursts is constrained by an upper limit at 1.1 mm. Radio emission is ~10x above the far-infrared/radio correlation, indicating an active galactic nucleus is present. An obscured AGN combined with starburst and evolved stellar components can account for the observations. If the black hole mass follows the local M_BH-bulge mass relation, the implied Eddington ratio is >~4. WISE 1814+3412 may be a heavily obscured object where the peak AGN activity occurred prior to the peak era of star formation.

We used the Karl G. Jansky Very Large Array (VLA) to image one primary beam area at 3 GHz with 8 arcsec FWHM resolution and 1.0 microJy/beam rms noise near the pointing center. The P(D) distribution from the central 10 arcmin of this confusion-limited image constrains the count of discrete sources in the 1 < S(microJy/beam) < 10 range. At this level the brightness-weighted differential count S^2 n(S) is converging rapidly, as predicted by evolutionary models in which the faintest radio sources are star-forming galaxies; and ~96$% of the background originating in galaxies has been resolved into discrete sources. About 63% of the radio background is produced by AGNs, and the remaining 37% comes from star-forming galaxies that obey the far-infrared (FIR) / radio correlation and account for most of the FIR background at lambda = 160 microns. Our new data confirm that radio sources powered by AGNs and star formation evolve at about the same rate, a result consistent with AGN feedback and the rough correlation of black hole and bulge stellar masses. The confusion at centimeter wavelengths is low enough that neither the planned SKA nor its pathfinder ASKAP EMU survey should be confusion limited, and the ultimate source detection limit imposed by "natural" confusion is < 0.01 microJy at 1.4 GHz. If discrete sources dominate the bright extragalactic background reported by ARCADE2 at 3.3 GHz, they cannot be located in or near galaxies and most are < 0.03 microJy at 1.4 GHz.

Water vapor megamasers from the center of active galaxies provide a powerful tool to trace accretion disks at sub-parsec resolution and, through an entirely geometrical method, measure direct distances to galaxies up to 200 Mpc. The Megamaser Cosmology Project (MCP) is formed by a team of astronomers with the aim of identifying new maser systems, and mapping their emission at high angular resolution to determine their distance. Two types of observations are necessary to measure a distance: single-dish monitoring to measure the acceleration of gas in the disk, and sensitive VLBI imaging to measure the angular size of the disk, measure the rotation curve, and model radial displacement of the maser feature. The ultimate goal of the MCP is to make a precise measurement of H0 by measuring such distances to at least 10 maser galaxies in the Hubble flow. We present here the preliminary results from a new maser system, Mrk 1419. Through a model of the rotation from the systemic masers assuming a narrow ring, and combining these results with the acceleration measurement from the Green Bank Telescope, we determine a distance to Mrk 1419 of 81\pm10 Mpc. Given that the disk shows a significant warp that may not be entirely traced by our current observations, more sensitive observations and more sophisticated disk modeling will be essential to improve our distance estimation to this galaxy.

The Hubble constant Ho describes not only the expansion of local space at redshift z ~ 0, but is also a fundamental parameter determining the evolution of the universe. Recent measurements of Ho anchored on Cepheid observations have reached a precision of several percent. However, this problem is so important that confirmation from several methods is needed to better constrain Ho and, with it, dark energy and the curvature of space. A particularly direct method involves the determination of distances to local galaxies far enough to be part of the Hubble flow through water vapor (H2O) masers orbiting nuclear supermassive black holes. The goal of this article is to describe the relevance of Ho with respect to fundamental cosmological questions and to summarize recent progress of the the `Megamaser Cosmology Project' (MCP) related to the Hubble constant.

We are using the VLBA at 8 GHz to observe a sample of 834 nearby 2MASS galaxies that are stronger than 100 mJy in the NVSS. The goals of the project are to detect (1) supermassive black holes significantly offset from the IR positions of the host galaxy bulges and (2) binary black holes.

We use the Expanded Very Large Array to image radio continuum emission from local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in 1 GHz windows centered at 4.7, 7.2, 29, and 36 GHz. This allows us to probe the integrated radio spectral energy distribution (SED) of the most energetic galaxies in the local universe. The 4-8 GHz flux densities agree well with previous measurements. They yield spectral indices \alpha \approx -0.67 (where F_\nu \propto \nu^\alpha) with \pm 0.15 (1\sigma) scatter, typical of nonthermal (synchrotron) emission from star-forming galaxies. The contrast of our 4-8 GHz data with literature 1.5 and 8.4 GHz flux densities gives further evidence for curvature of the radio SED of U/LIRGs. The SED appears flatter near \sim 1 GHz than near \sim 6 GHz, suggesting significant optical depth effects at the lower frequencies. The high frequency (28-37 GHz) flux densities are low compared to extrapolations from the 4-8 GHz data. We confirm and extend to higher frequency a previously observed deficit of high frequency radio emission for luminous starburst galaxies.

Despite decades of study, it remains unclear whether there are distinct radio-loud and radio-quiet populations of quasi-stellar objects (QSOs). Early studies were limited by inhomogeneous QSO samples, inadequate sensitivity to probe the radio-quiet population, and degeneracy between redshift and luminosity for flux-density-limited samples. Our new 6 GHz EVLA observations allow us for the first time to obtain nearly complete (97%) radio detections in a volume-limited color-selected sample of 179 QSOs more luminous than M_i = -23 from the Sloan Digital Sky Survey (SDSS) Data Release Seven in the narrow redshift range 0.2 < z < 0.3. The dramatic improvement in radio continuum sensitivity made possible with the new EVLA allows us, in 35 minutes of integration, to detect sources as faint as 20 microJy, or log[L_6 (W/Hz)] ~ 21.5 at z = 0.25, well below the radio luminosity, log[L_6 (W/Hz)] ~ 22.5, that separates star-forming galaxies from radio-loud active galactic nuclei (AGNs) driven by accretion onto a super-massive black hole. We calculate the radio luminosity function (RLF) for these QSOs using three constraints: (a) EVLA 6 GHz observations for log[L_6 (W/Hz)] < 23.5, (b) NRAO-VLA Sky Survey (NVSS) observations for log[L_6 (W/Hz)] > 23.5, and (c) the total number of SDSS QSOs in our volume-limited sample. We show that the RLF can be explained as a superposition of two populations, dominated by AGNs at the bright end and star formation in the QSO host galaxies at the faint end.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of H\alpha+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, H\alpha+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

Observations of H$_2$O masers from circumnuclear disks in active galaxies for the Megamaser Cosmology Project allow accurate measurement of the mass of supermassive black holes (BH) in these galaxies. We present the Very Long Baseline Interferometry (VLBI) images and kinematics of water maser emission in six active galaxies: NGC~1194, NGC~2273, NGC~2960 (Mrk~1419), NGC~4388, NGC~6264 and NGC~6323. We use the Keplerian rotation curves of these six megamaser galaxies, plus a seventh previously published, to determine accurate enclosed masses within the central $\sim0.3$ pc of these galaxies, smaller than the radius of the sphere of influence of the central mass in all cases. We also set lower limits to the central mass densities of between 0.12 and 60 $\times 10^{10} M_{\odot}$~pc$^{-3}$. For six of the seven disks, the high central densities rule out clusters of stars or stellar remnants as the central objects, and this result further supports our assumption that the enclosed mass can be attributed predominantly to a supermassive black hole. The seven BHs have masses ranging between 0.76 and 6.5$\times$10$^7 M_{\odot}$. The BH mass errors are $\approx11$\%, dominated by the uncertainty of the Hubble constant. We compare the megamaser BH mass determination with other BH mass measurement techniques. The BH mass based on virial estimation in four galaxies is consistent with the megamaser BH mass given the latest empirical value of $\langle f \rangle$, but the virial mass uncertainty is much greater. MCP observations continue and we expect to obtain more maser BH masses in the future.

The black hole (BH)-bulge correlations have greatly influenced the last decade of effort to understand galaxy evolution. Current knowledge of these correlations is limited predominantly to high BH masses (M_BH> 10^8 M_sun) that can be measured using direct stellar, gas, and maser kinematics. These objects, however, do not represent the demographics of more typical L< L* galaxies. This study transcends prior limitations to probe BHs that are an order of magnitude lower in mass, using BH mass measurements derived from the dynamics of H_2O megamasers in circumnuclear disks. The masers trace the Keplerian rotation of circumnuclear molecular disks starting at radii of a few tenths of a pc from the central BH. Modeling of the rotation curves, presented by Kuo et al. (2010), yields BH masses with exquisite precision. We present stellar velocity dispersion measurements for a sample of nine megamaser disk galaxies based on long-slit observations using the B&C spectrograph on the Dupont telescope and the DIS spectrograph on the 3.5m telescope at Apache Point. We also perform bulge-to-disk decomposition of a subset of five of these galaxies with SDSS imaging. The maser galaxies as a group fall below the M_BH-sigma* relation defined by elliptical galaxies. We show, now with very precise BH mass measurements, that the low-scatter power-law relation between M_BH and sigma* seen in elliptical galaxies is not universal. The elliptical galaxy M_BH-sigma* relation cannot be used to derive the BH mass function at low mass or the zeropoint for active BH masses. The processes (perhaps BH self-regulation or minor merging) that operate at higher mass have not effectively established an M_BH-sigma* relation in this low-mass regime.

The Megamaser Cosmology Project (MCP) aims to determine H0 by measuring angular-diameter distances to galaxies in the Hubble flow using observations of water vapor megamasers in the circumnuclear accretion disks of active galaxies. The technique is based only on geometry and determines H0 in one step, independent of standard candles and the extragalactic distance ladder. In Paper I we presented a VLBI map of the maser emission from the Seyfert 2 galaxy UGC 3789. The map reveals an edge-on, sub-parsec disk in Keplerian rotation, analogous to the megamaser disk in NGC 4258. Here we present 3.2 years of monthly GBT observations of the megamaser disk in UGC 3789. We use these observations to measure the centripetal accelerations of both the systemic and high-velocity maser components. The measured accelerations suggest that maser emission lines near the systemic velocity originate on the front side of the accretion disk, primarily from segments of two narrow rings. Adopting a two-ring model for the systemic features, we determine the angular-diameter distance to UGC 3789 to be 49.9 +/- 7.0 Mpc. This is the most accurate geometric distance yet obtained to a galaxy in the Hubble flow. Based on this distance, we determine H0 = 69 +/- 11 km/s/Mpc. We also measure the mass of the central black hole to be 1.09 x 10^7 solar masses +/- 14%. With additional observations the uncertainty in the distance to this galaxy can be reduced to under 10%. Observations of megamaser disks in other galaxies will further reduce the uncertainty in H0 as measured by the MCP.

We report on the Ka-band (26-40 GHz) emission properties for 10 star-forming regions in the nearby galaxy NGC 6946. From a radio spectral decomposition, we find that the 33 GHz flux densities are typically dominated by thermal (free-free) radiation. However, we also detect excess Ka-band emission for an outer-disk star-forming region relative to what is expected given existing radio, submillimeter, and infrared data. Among the 10 targeted regions, measurable excess emission at 33 GHz is detected for half of them, but in only one region is the excess found to be statistically significant ($\approx7\sigma$). We interpret this as the first likely detection of so called `anomalous' dust emission outside of the Milky Way. We find that models explaining this feature as the result of dipole emission from rapidly rotating ultrasmall grains are able to reproduce the observations for reasonable interstellar medium conditions. While these results suggest that the use of Ka-band data as a measure of star formation activity in external galaxies may be complicated by the presence of anomalous dust, it is unclear how significant a factor this will be for globally integrated measurements as the excess emission accounts for $\la$10% of the total Ka-band flux density from all 10 regions.

Proper interpretation and understanding of astronomical data requires good knowledge of the data acquisition process. The increase in remote observing, queue observing, and the availability of large archived data products risk insulating astronomers from the telescope, potentially reducing their familiarity with the observational techniques crucial in understanding the data. Learning fundamental observing techniques can be done in at least three ways: 1) College and university operated observing facilities, 2) Student involvement in national facilities through competitive proposals, 3) Programs at national facilities to increase upper-level undergraduate and graduate student exposure to telescopes. We encourage both national organizations and universities to include programs and funding aimed at supporting hands-on experience with telescopes through the three methods mentioned.

Wide-angle surveys have been an engine for new discoveries throughout the modern history of astronomy, and have been among the most highly cited and scientifically productive observing facilities in recent years. This trend is likely to continue over the next decade, as many of the most important questions in astrophysics are best tackled with massive surveys, often in synergy with each other and in tandem with the more traditional observatories. We argue that these surveys are most productive and have the greatest impact when the data from the surveys are made public in a timely manner. The rise of the "survey astronomer" is a substantial change in the demographics of our field; one of the most important challenges of the next decade is to find ways to recognize the intellectual contributions of those who work on the infrastructure of surveys (hardware, software, survey planning and operations, and databases/data distribution), and to make career paths to allow them to thrive.

The Megamaser Cosmology Project (MCP) seeks to measure the Hubble Constant (Ho) in order to improve the extragalactic distance scale and constrain the nature of dark energy. We are searching for sources of water maser emission from AGN with sub-pc accretion disks, as in NGC 4258, and following up these discoveries with Very Long Baseline Interferometric (VLBI) imaging and spectral monitoring. Here we present a VLBI map of the water masers toward UGC 3789, a galaxy well into the Hubble Flow. We have observed masers moving at rotational speeds up to 800 km/s at radii as small as 0.08 pc. Our map reveals masers in a nearly edge-on disk in Keplerian rotation about a 10^7 Msun supermassive black hole. When combined with centripetal accelerations, obtained by observing spectral drifts of maser features (to be presented in Paper II), the UGC 3789 masers may provide an accurate determination of Ho, independent of luminosities and metallicity and extinction corrections.

The Very Large Array (VLA) Low-frequency Sky Survey (VLSS) has imaged 95% of the 3*pi sr of sky north of declination = -30 degrees at a frequency of 74 MHz (4 meter wavelength). The resolution is 80" (FWHM) throughout, and the typical RMS noise level is ~0.1 Jy/beam. The typical point-source detection limit is 0.7 Jy/beam and so far nearly 70,000 sources have been catalogued. This survey used the 74 MHz system added to the VLA in 1998. It required new imaging algorithms to remove the large ionospheric distortions at this very low frequency throughout the entire ~11.9 degree field of view. This paper describes the observation and data reduction methods used for the VLSS and presents the survey images and source catalog. All of the calibrated images and the source catalog are available online (http://lwa.nrl.navy.mil/VLSS) for use by the astronomical community.

The Spitzer Space Telescope (formerly SIRTF) extragalactic First-Look Survey covered about 5 square degrees centered on J2000 17:18 +59:30 in order to characterize the infrared sky with high sensitivity. We used the 100-m Green Bank Telescope to image the 21cm Galactic HI emission over a 3x3 degree field covering this position with an effective angular resolution of 9.8 arcmin and a velocity resolution of 0.62 km/s. In the central square degree of the image the average column density is N(HI) = 2.5 x 10^{20} cm-2 with an rms fluctuation of 0.3 x 10^{20}. The Galactic HI in this region has a very interesting structure. There is a high-velocity cloud, several intermediate-velocity clouds (one of which is probably part of the Draco nebula), and narrow-line low velocity filaments. The HI emission shows a strong and detailed correlation with dust. Except for the high-velocity cloud, all features in the HI map have counterparts in an E(B-V) map derived from infrared data. Relatively high E(B-V)/N(HI) ratios in some directions suggest the presence of molecular gas. The best diagnostic of such regions is the peak HI line brightness temperature, not the total N(HI): directions where Tb > 12 K have E(B-V)/N(HI) significantly above the average value. The data corrected for stray radiation have been released via the Web.

Using the SPITZER SPACE TELESCOPE and the VLA, we present the first DIRECT evidence that the well-known far-infrared/radio correlation is valid to cosmologically significant redshift. We also confirm, with improved statistics compared with previous surveys, a similar result for the Mid-IR/radio correlation. We explore the dependence of monochromatic q_{24} and q_{70} on z. The results were obtained by matching Spitzer sources at 24 and 70 microns with VLA 1.4 GHz micro-Jy radio sources obtained for the SPITZER FLS. Spectroscopic redshifts have been obtained for over 500 matched IR/radio sources using observations at WIYN, Keck and archival SDSS data extending out to z > 2. We find that q_{24} shows significantly more dispersion than q_{70}. By comparing the observed fluxes at 70, 24 and 4.5 microns with a library of SED templates, we find that the larger dispersion in q_{24} is predictable in terms of systematic variations in SED shape throughout the population. Although the models are not able to encompass the full range of observed behavior (both the presence of either extremely flat or extremely steep IR SEDs), the fitting parameters were used to `k-correct' the higher-z galaxies which resulted in a reduced scatter in q. For comparison, we also corrected these data using the SED for M82. The results for 24 and 70 microns provide strong consistent evidence for the universality of the mid-IR/radio and far-IR/radio correlations out to redshifts of at least z = 1.

The NRAO Very Long Baseline Array (VLBA) was used at 1.4 GHz to image 20 compact sources in the inner 1-degree radius of the Spitzer VLA First-Look Survey. Twelve sources were detected with the VLBA, with peak flux densities above 2 mJy at 9-mas resolution. The positions of these VLBA detections have been measured in the International Celestial Reference Frame (ICRF) to 1-sigma accuracies of 8 mas East-West and 4 mas North-South. These positions can help anchor ancillary data for the First-Look Survey to the ICRF and are being reported quickly to the community.

We have searched the Faint Images of the Radio Sky at Twenty centimeters (FIRST) and the NRAO VLA Sky Survey (NVSS) 1.4 GHz radio surveys for sources that are coincident with emission-line galaxy (ELG) candidates from the KPNO International Spectroscopic Survey (KISS). A total of 207 of the 2157 KISS ELGs (~10%) in the first two H-alpha-selected survey lists were found to possess radio detections in FIRST and/or NVSS. Follow-up spectra exist for all of the radio detections, allowing us to determine the activity type (star-forming vs. AGN) for the entire sample. We explore the properties of the radio-detected KISS galaxies in order to gain a better insight into the nature of radio-emitting galaxies in the local universe (z < 0.1). No dwarf galaxies were detected, despite the large numbers of low-luminosity galaxies present in KISS, suggesting that lower mass, lower luminosity objects do not possess strong galaxian-scale magnetic fields. Due to the selection technique used for KISS, our radio ELGs represent a quasi-volume-limited sample, which allows us to develop a clearer picture of the radio galaxy population at low redshift. Nearly 2/3rds of the KISS radio galaxies are starburst/star-forming galaxies, which is in stark contrast to the results of flux-limited radio surveys that are dominated by AGNs and elliptical galaxies (i.e., classic radio galaxies). While there are many AGNs among the KISS radio galaxies, there are no objects with large radio powers in our local volume. We derive a radio luminosity function (RLF) for the KISS ELGs that agrees very well with previous RLFs that adequately sample the lower-luminosity radio population.

New multi-band observations of the famous compact group of galaxies Stephan's Quintet (SQ) are presented and analyzed. These include far infrared (FIR) images at 60$\mu m$ and 100$\mu m$ (ISOPHOT C-100 camera), radio continuum images at 1.4 GHz (VLA B-array) and 4.86 GHz (VLA C-array), and long-slit optical spectrographs (Palomar $200"$ telescope). With these new data, we aim to learn more about the X-ray/radio ridge in the middle of the intragroup medium (IGM) and the IGM starburst SQ-A, both are likely to be caused by the high speed collision ($\sim 900$ km s$^{-1}$) between the intruder galaxy NGC 7318b ($v = 5700$ km s$^{-1}$) and the IGM ($v = 6600$ km s$^{-1}$).

Three nearby galaxies which have abnormally high infrared to radio continuum ratios, NGC 1377, IC 1953 and NGC 4491, are investigated with a view to understanding the physical origin of their peculiarity. We review the existing data and present new radio continuum measurements along with near-infrared integral-field spectroscopy and molecular gas observations. The three galaxies have low luminosities but starburst-like infrared colors; in NGC 1377, no synchrotron emission is detected at any wavelength; in IC 1953, the observed synchrotron component is attributable to the spiral disk alone, and is lacking in the central regions; the radio spectrum of NGC 4491 is unusually flat. We also compare and contrast them with NGC 4418, a heavily extinguished galaxy which shares some attributes with them. After examining various scenarios, we conclude that these galaxies are most likely observed within a few Myr of the onset of an intense star formation episode after being quiescent for at least ~ 100 Myr. This starburst, while heating the dust, has not produced optical signatures nor a normal amount of cosmic rays yet. We briefly discuss the statistics of such galaxies and what they imply for star formation surveys.

New infrared and millimeter observations from Keck, Palomar, ISO, and OVRO and archival data from the NRAO VLA and IRAS are presented for the low metallicity blue compact dwarf galaxy SBS 0335-052. Mid-infrared imaging shows this young star-forming system is compact (0.31"; 80 pc) at 12.5 microns. The large Br-gamma equivalent width (235 Angstroms) measured from integral field spectroscopy is indicative of a ~5 Myr starburst. The central source appears to be optically thin in emission, containing both a warm (~80 K) and a hot (~210 K) dust component, and the overall interstellar radiation field is quite intense, about 10,000 times the intensity in the solar neighborhood. CO emission is not detected, though the galaxy shows an extremely high global H I gas-to-dust mass ratio, high even for blue compact dwarfs. Finally, the galaxy's mid-infrared-to-optical and mid-to-near-infrared luminosity ratios are quite high, whereas its far-infrared-to-radio and far-infrared-to-optical flux ratios are surprisingly similar to what is seen in normal star-forming galaxies. The relatively high bolometric infrared-to-radio ratio is more easily understood in the context of such a young system with negligible nonthermal radio continuum emission. These new lines of evidence may outline features common to primordial galaxies found at high redshift.

Observations for seven infrared luminous starburst galaxies are reported in the mid-infrared from 8 - 18 microns using the Keck Telescopes with spatial resolution approaching the diffraction limit. All of the galaxies observed show evidence of strong interactions based on optical morphologies. For these galaxies, a substantial fraction, usually more than 50%, of the infrared luminosity is generated in regions ranging in sizes from 100pc -- 1 Kpc. Nuclear starbursts often dominate the infrared luminosity, but this is not always true. In some galaxies, most notably NGC 6090, substantial infrared luminosity greatly in excess of the nuclear luminosity is generated in regions associated with the physical interaction between two galaxies. The radio emission is a good tracer of the location of high luminosity young stars. The visual/ultraviolet radiation output of the nearby star forming galaxies is dominated by emission from regions that are generally not producing the copious infrared luminosity of the systems. The regions of high infrared luminosity in local starburst galaxies are significantly smaller than the galaxies as a whole. The integrated spectral energy distributions (SEDs) of these galaxies are very different from the SEDs of the regions of star formation. If the SEDs of star-forming regions in these galaxies reflect the SEDs found in forming galaxies at high redshift, the distant galaxies should be dominated by the mid and far-infrared luminosity output far more than the integrated luminous output of nearby starburst galaxies would suggest.

The radio counterparts to the IRAS Redshift Survey galaxies are identified in the NRAO VLA Sky Survey (NVSS) catalog. Our new catalog of the infrared flux-limited (S60 >= 2 Jy) complete sample of 1809 galaxies lists accurate radio positions, redshifts, and 1.4 GHz radio and IRAS fluxes. This sample is six times larger in size and five times deeper in redshift coverage (to z ~ 0.15) compared with those used in earlier studies of the radio and far-infrared (FIR) properties of galaxies in the local volume. The well known radio-FIR correlation is obeyed by the overwhelming majority (>98%) of the infrared-selected galaxies, and the radio AGNs identified by their excess radio emission constitute only about 1% of the sample, independent of the IR luminosity. These FIR-selected galaxies can account for the entire population of late-type field galaxies in the local volume, and their radio continuum may be used directly to infer the extinction-free star formation rate in most cases. Both the 1.4 GHz radio and 60 micron infrared luminosity functions are reasonably well described by linear sums of two Schechter functions, one representing normal, late-type field galaxies and the second representing starbursts and other luminous infrared galaxies. The integrated FIR luminosity density for the local volume is 4.8+/-0.5 x 10^7 L(sun) per Mpc^3, less than 10% of which is contributed by the luminous infrared galaxies with L(FIR) > 10^{11} L(sun) . The inferred extinction-free star formation density for the local volume is 0.015 +/- 0.005 solar mass per year per Mpc^3.

The Very Large Array (VLA) has been used in C-configuration to map an area \~0.3 square degrees at 1.4 GHz with 5sigma sensitivities of 0.305, 0.325, 0.380 and 0.450 mJy/beam over four equal subareas. Radio properties are presented for 62 detected sources. Deep optical imaging to Gunn r~25 mag using the Hale 5-m telescope covering ~0.21 square degrees is reported for a subset of 43 sources. This optical follow-up is much deeper than that of existing larger area radio surveys of similar radio sensitivity. Archival J, H and K-band photometry from the Two-Micron All Sky Survey (2MASS) is also presented. Using a robust likelihood-ratio technique, we optically identified 26 radio sources with probabilities >80%. Nine are uncertain/ambiguous detections and eight are empty fields. Comparisons with a stellar synthesis model that includes radio emission and dust reddening suggests that the near-infrared--to--optical emission in a small, bright sub-sample is reddened by `optically-thin' dust with absorption A_V ~ 2-2.5 mag, regardless of morphological type. This is consistent with other, more direct determinations of absorption. The radio--to--optical(--near-infrared) flux ratios of early-type galaxies require significant contamination in the radio by an active galactic nucleus (AGN), consistent with the current paradigm. Using our simple modeling approach, we also discuss a potential diagnostic for selecting Ultraluminous Infrared Galaxies (ULIGS) to z~1.6 from microJansky radio surveys.

We observed 17 nearby low luminosity FR I radio galaxies using the NRAO Very Long Baseline Array (VLBA) at 1.67 GHz, as part of a multi-wavelength study of a complete sample of 21 sources selected by radio flux density from the Uppsala General Catalogue of Galaxies. We detected radio emission from all 17 galaxies. At a FWHM resolution of ~ 10x4 mas, five galaxies show only an unresolved radio core, ten galaxies show core-jet structures, and two galaxies show twin-jet structures. Comparing these VLBA images with images previously obtained with the NRAO VLA, we find that all detected VLBA jets are well-aligned on parsec scales with the VLA jets on kilo parsec scales, and that the jet to counter-jet surface brightness ratios, or the sidedness, decreases systematically with increasing distance along the jet. We attribute the sidedness to the Doppler boosting effect and its decline to the deceleration of the jets. We show that a distribution of Lorentz factor centered near Gamma = 5 can reproduce our VLBA detection statistics for core, core-jet and twin-jet sources. We also note that the luminosity per unit length, L_j, of the VLBA jets drops quickly with distance, r, along the jet, approximately as L_j being proportional to r^(-2.0). We discuss three different mechanisms to explain this jet fading: (1) the decrease of Doppler boosting due to jet deceleration, (2) synchrotron losses, and (3) expansion losses in constant velocity but adiabatically spreading jets. Mechanisms (1) and (2) are inconsistent with the observations, while mechanism (3) is consistent with the observations provided the magnetic field lines in the jets are aligned perpendicular to the jet axis. This implies that the deceleration of the jets required by the Unified Scheme does not occur on the tens of parsec scales, but must occur on larger scales.

We report on an infrared 1 to 17 micron study of the nearby (cz = 4600 km/s) interacting spiral galaxy system, UGC 12914/12915, using the ground-based Palomar 200" telescope and PFIRCAM near-infrared detector and space-based mid-infrared imaging and spectral observations using ISOCAM and Phot-S on the Infrared Space Observatory (ISO). Strong line emission from aromatic band features at 6.2, 7.7 and 11.3 microns are seen in both the mid-IR imaging and Phot-S spectro-photometry centered on the nuclei. In the connecting bridge, the mid-infrared imaging reveals dust grains intermixed with the atomic hydrogen gas. The heating mechanism for the hot dust is likely to be UV photons diffusing out from the galactic disks and the H II complex located along the extreme northeastern portion of the bridge: The dust emission, or mid-infrared intensity per atomic hydrogen column density ratio, is consistent with heating from the local (bridge) interstellar radiation field.

We cross-identified the ROSAT Bright Source Catalog (RBSC) and the NRAO VLA Sky Survey (NVSS) to construct the RBSC-NVSS sample of the brightest X-ray sources (>= 0.1 counts/s or ~1E-12 ergs/cm/cm/s in the 0.1-2.4 keV band) that are also radio sources (S >= 2.5 mJy at 1.4 GHz) in the 7.8 sr of extragalactic sky with |b| > 15 degrees. and delta > -40 degrees. The sky density of NVSS sources is low enough that they can be reliably identified with RBSC sources having average rms positional uncertainties = 10 arcsec. We used the more accurate radio positions to make reliable X-ray/radio/optical identifications down to the POSS plate limits. We obtained optical spectra for many of the bright identifications lacking published redshifts. The resulting X-ray/radio sample is unique in its size (N ~ 1500 objects), composition (a mixture of nearly normal galaxies, Seyfert galaxies, quasars, and clusters), and low average redshift (<z> ~ 0.1).

Observations of ultraluminous infrared galaxies (ULIRGs) with an achieved resolution approaching the diffraction limit in the mid-infrared from 8 - 25 $\mu$m using the Keck Telescopes are reported. We find extremely compact structures, with spatial scales of $< 0.3''$ (diameter) in six of the seven ULIRGs observed. These compact sources emit between 30% and 100% of the mid-infrared energy from these galaxies. We have utilized the compact mid-infrared structures as a diagnostic of whether an AGN or a compact (100 -- 300 pc) starburst is the primary power source in these ULIRGs. In Markarian 231, the upper limit on the diameter of the 12.5 $\mu$m source, 0.13$''$, shows that the size of the infrared source must increase with increasing wavelength, consistent with AGN models. In IRAS 05189-2524 and IRAS 08572+3915 there is strong evidence that the source size increases with increasing wavelength. This suggests heating by a central source rather than an extended luminosity source, consistent with the optical classification as an AGN. The compact mid-infrared sources seen in the other galaxies cannot be used to distinguish the ultimate luminosity source. If these ULIRGs are powered by compact starbursts, the star formation rates seen in the central few hundred parsecs far exceed the global rates seen in nearby starburst galaxies, and approach the surface brightness of individual clusters in nearby starburst galaxies.

We present a new sample of Parkes half-Jansky flat-spectrum radio sources having made a particular effort to find any previously unidentified sources. The sample contains 323 sources selected according to a flux limit of 0.5 Jy at 2.7 GHz, a spectral index measured between 2.7 and 5.0 GHz of alpha > -0.5 (where S(v) is proportional to v to the power of alpha), Galactic latitude abs(b) > 20 deg and -45 deg < Declination (B1950) < +10 deg. The sample was selected from a region 3.90 steradians in area. We have obtained accurate radio positions for all the unresolved sources in this sample and combined these with accurate optical positions from digitised photographic sky survey data to check all the optical identifications. We report new identifications based on \R- and \Kn-band imaging and new spectroscopic measurements of many of the sources. We present a catalogue of the 323 sources of which 321 now have identified optical counterparts and 277 have measured spectral redshifts.