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

Polarized foreground emission is a potential contaminant of attempts to measure the fluctuation power spectrum of highly redshifted 21 cm HI emission from the epoch of reionization. Using the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), we present limits on the observed power spectra of all four Stokes parameters in two frequency bands, centered at 126 MHz ($z=10.3$) and 164 MHz ($z=7.66$) for a three-month observing campaign of a 32-antenna deployment, for which unpolarized power spectrum results have been reported at $z=7.7$ (Parsons et al 2014) and $7.5 < z < 10.5$ (Jacobs et al 2014). The power spectra in this paper are processed in the same way as in those works, and show no definitive detection of polarized power. This non-detection appears to be largely due to the suppression of polarized power by ionospheric rotation measure fluctuations, which strongly affect Stokes Q and U. We are able to show that the net effect of polarized leakage is a negligible contribution at the levels of of the limits reported in Parsons et al 2014 and Jacobs et al 2014.

Current-generation low frequency interferometers constructed with the objective of detecting the high-redshift 21 cm background, aim to generate power spectra of the brightness-temperature contrast of neutral hydrogen in primordial intergalactic medium. Two-dimensional power spectra (power in Fourier modes parallel and perpendicular to the line of sight) formed from interferometric visibilities have been shown to delineate a boundary between spectrally-smooth foregrounds (known as the wedge) and spectrally-structured 21 cm background emission (the EoR-window). However, polarized foregrounds are known to possess spectral structure due to Faraday rotation, which can leak into the EoR window. In this work, we create and analyze 2D power spectra from the PAPER-32 imaging array in Stokes I, Q, U and V. These allow us to observe and diagnose systematic effects in our calibration at high signal-to-noise within the Fourier space most relevant to EoR experiments. We observe well-defined windows in the Stokes visibilities, with Stokes Q, U and V power spectra sharing a similar wedge shape to that seen in Stokes I. With modest polarization calibration, we see no evidence that polarization calibration errors move power outside the wedge in any Stokes visibility, to the noise levels attained. Deeper integrations will be required to confirm that this behavior persists to the depth required for EoR detection.

The epoch of reionization power spectrum is expected to evolve strongly with redshift, and it is this variation with cosmic history that will allow us to begin to place constraints on the physics of reionization. The primary obstacle to the measurement of the EoR power spectrum is bright foreground emission. We present an analysis of observations from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) telescope which place new limits on the HI power spectrum over the redshift range of $7.5<z<10.5$, extending previously published single redshift results to cover the full range accessible to the instrument. To suppress foregrounds, we use filtering techniques that take advantage of the large instrumental bandwidth to isolate and suppress foreground leakage into the interesting regions of $k$-space. Our 500 hour integration is the longest such yet recorded and demonstrates this method to a dynamic range of $10^4$. Power spectra at different points across the redshift range reveal the variable efficacy of the foreground isolation. Noise limited measurements of $\Delta^2$ at $k=$0.2hMpc$^{-1}$ and z$=7.55$ reach as low as (48mK)$^2$ ($1\sigma$). We demonstrate that the size of the error bars in our power spectrum measurement as generated by a bootstrap method is consistent with the fluctuations due to thermal noise. Relative to this thermal noise, most spectra exhibit an excess of power at a few sigma. The likely sources of this excess include residual foreground leakage, particularly at the highest redshift, and unflagged RFI. We conclude by discussing data reduction improvements that promise to remove much of this excess.

We present new constraints on the 21cm Epoch of Reionization (EoR) power spectrum derived from 3 months of observing with a 32-antenna, dual-polarization deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. In this paper, we demonstrate the efficacy of the delay-spectrum approach to avoiding foregrounds, achieving over 8 orders of magnitude of foreground suppression (in $\textrm{mK}^2$). Combining this approach with a procedure for removing off-diagonal covariances arising from instrumental systematics, we achieve a best 2-sigma upper limit of $(41\,\textrm{mK})^2$ for $k=0.27 h\textrm{Mpc}^{-1}$ at $z=7.7$. This limit falls within an order of magnitude of the brighter predictions of the expected 21cm EoR signal level. Using the upper limits set by these measurements, we generate new constraints on the brightness temperature of 21cm emission in neutral regions for various reionization models. We show that for several ionization scenarios, our measurements are inconsistent with cold reionization. That is, heating of the neutral intergalactic medium (IGM) is necessary to remain consistent with the constraints we report. Hence, we have suggestive evidence that by $z=7.7$, the HI has been warmed from its cold primordial state, probably by X-rays from high-mass X-ray binaries or mini-quasars. The strength of this evidence depends on the ionization state of the IGM, which we are not yet able to constrain. This result is consistent with standard predictions for how reionization might have proceeded.

We present a catalog of spectral measurements covering a 100-200 MHz band for 32 sources, derived from observations with a 64-antenna deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. For transit telescopes such as PAPER, calibration of the primary beam is a difficult endeavor, and errors in this calibration are a major source of error in the determination of source spectra. In order to decrease reliance on accurate beam calibration, we focus on calibrating sources in a narrow declination range from -46d to -40d. Since sources at similar declinations follow nearly identical paths through the primary beam, this restriction greatly reduces errors associated with beam calibration, yielding a dramatic improvement in the accuracy of derived source spectra. Extrapolating from higher frequency catalogs, we derive the flux scale using a Monte-Carlo fit across multiple sources that includes uncertainty from both catalog and measurement errors. Fitting spectral models to catalog data and these new PAPER measurements, we derive new flux models for Pictor A and 31 other sources at nearby declinations. 90% of these confirm and refine a power-law model for flux density. Of note is the new Pictor A flux model, which is accurate to 1.4% and shows, in contrast to previous models, that between 100 MHz and 2 GHz, the spectrum of Pictor A is consistent with a single power law given by a flux at 150 MHz of 382+/-5.4 Jy, and a spectral index of -0.76+/-0.01. This accuracy represents an order of magnitude improvement over previous measurements in this band, and is limited by the uncertainty in the catalog measurements used to estimate the absolute flux scale. The simplicity and improved accuracy of Pictor A's spectrum make it an excellent calibrator for experiments seeking to measure 21cm emission from the Epoch of Reionization.

We present new observations with the Precision Array for Probing the Epoch of Reionization (PAPER) with the aim of measuring the properties of foreground emission for 21cm Epoch of Reionization experiments at 150 MHz. We focus on the footprint of the foregrounds in cosmological Fourier space to understand which modes of the 21cm power spectrum will most likely be compromised by foreground emission. These observations confirm predictions that foregrounds can be isolated to a "wedge"-like region of 2D (k-perpendicular, k-parallel)-space, creating a window for cosmological studies at higher k-parallel values. We also find that the emission extends past the nominal edge of this wedge due to spectral structure in the foregrounds, with this feature most prominent on the shortest baselines. Finally, we filter the data to retain only this "unsmooth" emission and image specific k-parallel modes of it. The resultant images show an excess of power at the lowest modes, but no emission can be clearly localized to any one region of the sky. This image is highly suggestive that the most problematic foregrounds for 21cm EoR studies will not be easily identifiable bright sources, but rather an aggregate of fainter emission.

We present observations taken with the Precision Array for Probing the Epoch of Reionization (PAPER) of the Centaurus A field in the frequency range 114 to 188 MHz. The resulting image has a 25' resolution, a dynamic range of 3500 and an r.m.s. of 0.5 Jy\beam (for a beam size of 25' x 23'). A spectral index map of Cen A is produced across the full band. The spectral index distribution is qualitatively consistent with electron reacceleration in regions of excess turbulence in the radio lobes, as previously identified morphologically. Hence, there appears to be an association of 'severe weather' in radio lobes with energy input into the relativistic electron population. We compare the PAPER large scale radio image with the X-ray image from the ROSAT All Sky Survey. There is a tentative correlation between radio and X-ray features at the end of the southern lobe, some 200 kpc from the nucleus, as might be expected from inverse Compton scattering of the CMB by the relativistic electrons also responsible for the radio synchrotron emission. The magnetic fields derived from the (possible) IC and radio emission are of similar magnitude to fields derived under the minimum pressure assumptions, ~ 1 {\mu}G. However, the X-ray field is complex, with large scale gradients and features possibly unrelated to Cen A. If these X-ray features are unrelated to Cen A, then these fields are lower limits.

We present a new technique for calibrating the primary beam of a wide-field, drift-scanning antenna element. Drift-scan observing is not compatible with standard beam calibration routines, and the situation is further complicated by difficult-to-parametrize beam shapes and, at low frequencies, the sparsity of accurate source spectra to use as calibrators. We overcome these challenges by building up an interrelated network of source "crossing points" -- locations where the primary beam is sampled by multiple sources. Using the single assumption that a beam has 180 degree rotational symmetry, we can achieve significant beam coverage with only a few tens of sources. The resulting network of crossing points allows us to solve for both a beam model and source flux densities referenced to a single calibrator source, circumventing the need for a large sample of well-characterized calibrators. We illustrate the method with actual and simulated observations from the Precision Array for Probing the Epoch of Reionization (PAPER).

We present observations from the Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa, observed in May and September 2010. Using two nights of drift scanning with PAPER's 60\arcdeg\ FWHM beam we have made a map covering the entire sky below +10 degrees declination with an effective center frequency of 145 MHz, a 70-MHz bandwidth, and a resolution of 26\arcmin. A 4800 square-degree region of this large map with the lowest Galactic emission reaches an RMS of 0.7 Jy. We establish an absolute flux scale using sources from the 160-MHz Culgoora catalog. Using the 408-MHz Molonglo Reference Catalog (MRC) as a finding survey, we identify counterparts to 480 sources in our maps, and compare our fluxes to the MRC and to 332 sources in the Culgoora catalog. For both catalogs, the ratio of PAPER to catalog flux averages to 1, with a standard deviation of 50%. This measured variation is consistent with comparisons between independent catalogs observed at different bands. The PAPER data represent new 145-MHz flux measurements for a large number of sources in the band expected to encompass cosmic reionization, and represents a significant step toward establishing a model for removing foregrounds to the reionization signal.

We are developing the Precision Array for Probing the Epoch of Reionization (PAPER) to detect 21cm emission from the early Universe, when the first stars and galaxies were forming. We describe the overall experiment strategy and architecture and summarize two PAPER deployments: a 4-antenna array in the low-RFI environment of Western Australia and an 8-antenna array at our prototyping site in Green Bank, WV. From these activities we report on system performance, including primary beam model verification, dependence of system gain on ambient temperature, measurements of receiver and overall system temperatures, and characterization of the RFI environment at each deployment site. We present an all-sky map synthesized between 139 MHz and 174 MHz using data from both arrays that reaches down to 80 mJy (4.9 K, for a beam size of 2.15e-5 steradians at 154 MHz), with a 10 mJy (620 mK) thermal noise level that indicates what would be achievable with better foreground subtraction. We calculate angular power spectra ($C_\ell$) in a cold patch and determine them to be dominated by point sources, but with contributions from galactic synchrotron emission at lower radio frequencies and angular wavemodes. Although the cosmic variance of foregrounds dominates errors in these power spectra, we measure a thermal noise level of 310 mK at $\ell=100$ for a 1.46-MHz band centered at 164.5 MHz. This sensitivity level is approximately three orders of magnitude in temperature above the level of the fluctuations in 21cm emission associated with reionization.

Using the Green Bank Telescope, we conducted a ``snapshot'' survey for water maser emission toward the nuclei of 611 galaxies and detected eight new sources. The sample consisted of nearby (v < 5000 km/s) and luminous (M_B < -19.5) galaxies, some with known nuclear activity but most not previously known to host AGNs. Our detections include both megamasers associated with AGNs and relatively low luminosity masers probably associated with star formation. The detection in UGC 3789 is particularly intriguing because the spectrum shows both systemic and high-velocity lines indicative of emission from an AGN accretion disk seen edge-on. Based on six months of monitoring, we detected accelerations among the systemic features ranging from 2 to 8 km/s/yr, the larger values belonging to the most redshifted systemic components. High-velocity maser lines in UGC 3789 show no detectable drift over the same period. Although UGC 3789 was not known to be an AGN prior to this survey, the presence of a disk maser is strong evidence for nuclear activity, and an optical spectrum obtained later has confirmed it. With follow up observations, it may be possible to measure a geometric distance to UGC 3789.

We present multi-frequency VLBA observations of two polarized Compact Symmetric Objects (CSOs), J0000+4054 and J1826+1831, and a polarized CSO candidate, J1915+6548. Using the wavelength-squared dependence of Faraday rotation, we obtained rotation measures (RMs) of -180 \pm 10 rad m^-2 and 1540 \pm 7 rad m^-2 for the latter two sources. These are lower than what is expected of CSOs (several 1000 rad m^-2) and, depending on the path length of the Faraday screens, require magnetic fields from 0.03 to 6 \mu G. These CSOs may be more heavily affected by Doppler boosting than their unpolarized counterparts, suggesting that a jet-axis orientation more inclined towards the line of sight is necessary to detect any polarization. This allows for low RMs if the polarized components are oriented away from the depolarizing circumnuclear torus. These observations also add a fourth epoch to the proper motion studies of J0000+4054 and J1826+1831, constraining their kinematic age estimates to >610 yrs and 2600 \pm 490 yrs, respectively. The morphology, spectrum, and component motions of J1915+6548 are discussed in light of its new classification as a CSO candidate, and its angle to the line of sight (~50\deg) is determined from relativistic beaming arguments.

We present Very Long Baseline Array (VLBA) observations of the nucleus of NGC 1275, the central, dominant galaxy in the Perseus cluster of galaxies. These are the first observations to resolve the linearly polarized emission from 3C84, and from them we determine a Faraday rotation measure (RM) ranging from 6500 to 7500 rad/m^2 across the tip of the bright southern jet component. At 22 GHz some polarization is also detected from the central parsec of 3C84, indicating the presence of even more extreme RMs that depolarize the core at lower frequencies. The nature of the Faraday screen is most consistent with being produced by magnetic fields associated with the optical filaments of ionized gas in the Perseus Cluster.

We present multi-epoch VLBA observations of Compact Symmetric Objects (CSOs) from the COINS sample (CSOs Observed In the Northern Sky). These observations allow us to make estimates of, or place limits on, the kinematic ages for those sources with well-identified hot spots. This study significantly increases the number of CSOs with well-determined ages or limits. The age distribution is found to be sharply peaked under 500 years, suggesting that many CSOs die young, or are episodic in nature, and very few survive to evolve into FR II sources like Cygnus A. Jet components are found to have higher velocities than hot spots which is consistent with their movement down cleared channels. We also report on the first detections of significant polarization in two CSOs, J0000+4054 (2.1%) and J1826+1831 (8.8%). In both cases the polarized emission is found in jet components on the stronger side of the center of activity.