We present the discovery and subarcsecond localization of a new Fast Radio Burst with the Karl G. Jansky Very Large Array and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc/cm3. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hours of VLA observing and 153 hours of CHIME/FRB observing. We describe a suite of statistical and data quality tests we used to verify the significance of the event and its localization precision. Follow-up optical/infrared photometry with Keck and Gemini associate the FRB to a pair of galaxies with $\rm{r}\sim23$ mag. The false-alarm rate for radio transients of this significance that are associated with a host galaxy is roughly $3\times10^{-4}\ \rm{hr}^{-1}$. The two putative host galaxies have similar photometric redshifts of $z_{\rm{phot}}\sim0.6$, but different colors and stellar masses. Comparing the host distance to that implied by the dispersion measure suggests a modest (~ 50 pc/cm3) electron column density associated with the FRB environment or host galaxy/galaxies.

We report on the detection of a burst from FRB180916 by realfast/VLA and present software for interpreting fast radio bursts (FRB) periodicity. We demonstrate a range of periodicity analyses with bursts from FRB180916, FRB121102 and FRB180814. Our results for FRB180916 and FRB121102 are consistent with published results. For FRB180814, we did not detect any significant periodic episodes. The realfast-detected and other high-frequency bursts for FRB180916 tend to lie at the beginning of the activity window, indicating a possible phase-frequency relation. The python package $\texttt{frbpa}$ can be used to reproduce and expand on this analysis to test models for repeating FRBs.

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5\sigma$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3\sigma$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

The "Space VLBI 2020: Science and Technology Futures" meeting was the second in The Future of High-Resolution Radio Interferometry in Space series. The first meeting (2018 September 5--6; Noordwijk, the Netherlands) focused on the full range of science applications possible for very long baseline interferometry (VLBI) with space-based antennas. Accordingly, the observing frequencies (wavelengths) considered ranged from below 1~MHz (> 300 m) to above 300~GHz (< 1 mm). For this second meeting, the focus was narrowed to mission concepts and the supporting technologies to enable the highest angular resolution observations at frequencies of 30~GHz and higher (< 1 cm). This narrowing of focus was driven by both scientific and technical considerations. First, results from the RadioAstron mission and the Event Horizon Telescope (EHT) have generated considerable excitement for studying the inner portions of black hole (BH) accretion disks and jets and testing elements of the General Theory of Relativity (GR). Second, the technologies and requirements involved in space-based VLBI differ considerably between 100~MHz and 100~GHz; a related consideration is that there are a number of existing instruments or mission concepts for frequencies of approximately 100~MHz and below, while it has been some time since attention has been devoted to space VLBI at frequencies above 10~GHz. This conference summary attempts to capture elements of presentations and discussions that occurred.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

PSR J1012+5307, a millisecond pulsar in orbit with a helium white dwarf (WD), has been timed with high precision for about 25 years. One of the main objectives of this long-term timing is to use the large asymmetry in gravitational binding energy between the neutron star and the WD to test gravitational theories. Such tests, however, will be eventually limited by the accuracy of the distance to the pulsar. Here, we present VLBI (very long baseline interferometry) astrometry results spanning approximately 2.5 years for PSR J1012+5307, obtained with the Very Long Baseline Array as part of the MSPSRPI project. These provide the first proper motion and absolute position for PSR J1012+5307 measured in a quasi-inertial reference frame. From the VLBI results, we measure a distance of $0.83^{+0.06}_{-0.02}$kpc (all the estimates presented in the abstract are at 68% confidence) for PSR J1012+5307, which is the most precise obtained to date. Using the new distance, we improve the uncertainty of measurements of the unmodeled contributions to orbital period decay, which, combined with three other pulsars, places new constraints on the coupling constant for dipole gravitational radiation $\kappa_D=(-1.7\pm1.7)\times 10^{-4}$ and the fractional time derivative of Newton's gravitational constant $\dot{G}/G = -1.8^{\,+5.6}_{\,-4.7}\times 10^{-13}\,{\rm yr^{-1}}$ in the local universe. As the uncertainties of the observed decays of orbital period for the four leading pulsar-WD systems become negligible in $\approx10$ years, the uncertainties for $\dot{G}/G$ and $\kappa_D$ will be improved to $\leq1.5\times10^{-13}\,{\rm yr^{-1}}$ and $\leq1.0\times10^{-4}$, respectively, predominantly limited by the distance uncertainties.

We introduce DAYENU, a linear, spectral filter for HI intensity mapping that achieves the desirable foreground mitigation and error minimization properties of inverse co-variance weighting with minimal modeling of the underlying data. Beyond 21 cm power-spectrum estimation, our filter is suitable for any analysis where high dynamic-range removal of spectrally smooth foregrounds in irregularly (or regularly) sampled data is required, something required by many other intensity mapping techniques. We show that DAYENU enables the access of large-scale line-of-sight modes that are inaccessible to tapered DFT estimators. Since these modes have the largest SNRs, DAYENU significantly increases the sensitivity of 21 cm analyses over tapered Fourier transforms. Slight modifications allow us to use DAYENU as a linear replacement for iterative delay CLEANing (DAYENUREST). An interactive jupyter tutorial on using DAYENU can be found at https://github.com/HERA-Team/uvtools/blob/master/examples/linear_clean_demo.ipynb DAYENU's source code can be found at https://github.com/HERA-Team/uvtools/blob/master/uvtools/dspec.py

Radio emission propagating over an Earth-pulsar line of sight provides a unique probe of the intervening ionized interstellar medium (ISM). Variations in the integrated electron column density along this line of sight, or dispersion measure (DM), have been observed since shortly after the discovery of pulsars. As early as 2006, frequency-dependent dispersion measures have been observed and attributed to several possible causes. Ray-path averaging over different effective light-cone volumes through the turbulent ISM contributes to this effect as will DM misestimation due to radio propagation across compact lensing structures such as those caused by "extreme scattering events". We present methods to assess the variations in frequency-dependent dispersion measures due to the turbulent ISM versus these compact lensing structures along the line of sight. We analyze recent Low-Frequency Array (LOFAR) observations of PSR J2219+4754 to test the underlying physical mechanism of the observed frequency-dependent DM. Previous analyses have indicated the presence of strong lensing due to compact overdensities halfway between the Earth and pulsar. Instead we find the frequency dependence of the DM timeseries for PSR J2219+4754 is consistent with being due solely to ISM turbulence and there is no evidence for any extreme scattering event or small-scale lensing structure. The data show possible deviations from a uniform turbulent medium, suggesting that there may be an enhanced scattering screen near one of the two ends of the line of sight. We present this analysis as an example of the power of low-frequency observations to distinguish the underlying mechanisms in frequency-dependent propagation effects.

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

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

We use a semi-analytic model to explore the potential impact of a brief and violent period of radio-loud accretion onto black-holes (The Radio Scream) during the Cosmic Dawn on the HI hyperfine 21 cm signal. We find that radio emission from super-massive black hole seeds can impact the global 21 cm signal at the level of tens to hundreds of percent provided that they were as radio loud as $z\approx1$ black holes and obscured by gas with column depths of $N_\text{H}\gtrsim 10^{23}$ cm$^{-2}$. We determine plausible sets of parameters that reproduce some of the striking features of the EDGES absorption feature including its depth, timing, and side steepness while producing radio/X-ray backgrounds and source counts that are consistent with published limits. Scenarios yielding a dramatic 21 cm signature also predict large populations of $\sim \mu$Jy point sources that will be detectable in future deep surveys from the Square Kilometer Array (SKA). Thus, 21 cm measurements, complemented by deep point source surveys, have the potential to constrain optimistic scenarios where super-massive black-hole progenitors were radio-loud.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

We present the results of the first transient survey from the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) using 31 hr of data, in which we place the most constraining limits on the instantaneous transient surface density at timescales of 13 s to a few minutes and at frequencies below 100 MHz. The OVRO-LWA is a dipole array that images the entire viewable hemisphere with 58 MHz of bandwidth from 27 to 84 MHz at 13 s cadence. No transients are detected above a 6.5$\sigma$ flux density limit of 10.5 Jy, implying an upper limit to the transient surface density of $2.5\times10^{-8}$ deg$^{-2}$ at the shortest timescales probed, which is orders of magnitude deeper than has been achieved at sub-100 MHz frequencies and comparable flux densities to date. The nondetection of transients in the OVRO-LWA survey, particularly at minutes-long timescales, allows us to place further constraints on the rate of the potential population of transients uncovered by Stewart et al. (2016). From their transient rate, we expect a detection of $8.4^{+31.8}_{-8.0}$ events, and the probability of our null detection is $1.9^{+644}_{-1.9}\times10^{-3}$, ruling out a transient rate $>1.4\times10^{-4} \text{days}^{-1} \text{deg}^{-2}$ with 95% confidence at a flux density limit of 18.1 Jy, under the assumption of a flat spectrum and wide bandwidth. We discuss the implications of our nondetection for this population and further constraints that can be made on the source spectral index, intrinsic emission bandwidth, and resulting luminosity distribution.

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

We present results of a search for late-time radio emission and Fast Radio Bursts (FRBs) from a sample of type-I superluminous supernovae (SLSNe-I). We used the Karl G. Jansky Very Large Array to observe ten SLSN-I more than 5 years old at a frequency of 3 GHz. We searched fast-sampled visibilities for FRBs and used the same data to perform a deep imaging search for late-time radio emission expected in models of magnetar-powered supernovae. No FRBs were found. One SLSN-I, PTF10hgi, is detected in deep imaging, corresponding to a luminosity of $1.2\times10^{28}$ erg s$^{-1}$. This luminosity, considered with the recent 6 GHz detection of PTF10hgi in Eftekhari et al (2019), supports the interpretation that it is powered by a young, fast-spinning ($\sim$ ms spin period) magnetar with $\sim$ 15 Msun of partially ionized ejecta. Broadly, our observations are most consistent with SLSNe-I being powered by neutron stars with fast spin periods, although most require more free-free absorption than is inferred for PTF10hgi. We predict that radio observations at higher frequencies or in the near future will detect these systems and begin constraining properties of the young pulsars and their birth environments.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

We calculate the dispersion measure (DM) contributed by the intergalactic medium (IGM) to the total measured DM for fast radio bursts (FRBs). We use the MareNostrum Instituto de Ciencias del Espacio (MICE) Onion Universe simulation (Fosalba et al. 2008) to track the evolution of the dark matter particle density over a large range of redshifts. We convert this dark matter particle number density to the corresponding free electron density and then integrate it to find the DM as a function of redshift. This approach yields an intergalactic DM of $\rm DM_{\rm IGM}(z = 1) = 800^{+7000}_{-170}$ pc cm$^{-3}$, with the large errors representative of the structure in the IGM. We place limits on the redshifts of the current population of observed FRBs. We also use our results to estimate the host galaxy contribution to the DM for the first repeater, FRB 121102, and show that the most probable host galaxy DM contribution, $\rm DM_{\rm host} \approx 310$ pc cm$^{-3}$, is consistent with the estimate made using the Balmer emission lines in the spectrum of the host galaxy, $\rm DM_{\rm Balmer} = 324$ pc cm$^{-3}$ (Tendulkar et al. 2017). We also compare our predictions for the host galaxy contribution to the DM for the observations of FRB 180924 (Bannister et al. 2019) and FRB 190523 (Ravi et al. 2019), both of which have been localized to a host galaxy.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.

We present the results of a 7-12 GHz phased-array study of the Galactic center magnetar J1745-2900 with the Karl G. Jansky Very Large Array (VLA). Using data from two 6.5 hour observations from September 2014, we find that the average profile is comprised of several distinct components at these epochs and is stable over $\sim$day timescales and $\sim$GHz frequencies. Comparison with additional phased VLA data at 8.7 GHz shows significant profile changes on longer timescales. The average profile at 7-12 GHz is dominated by the jitter of relatively narrow pulses. The pulses in each of the four main profile components seen in September 2014 are uncorrelated in phase and amplitude, though there is a small but significant correlation in the occurrence of pulses in two of the profile components. Using the brightest pulses, we measure the dispersion and scattering parameters of J1745-2900. A joint fit of 38 pulses gives a 10 GHz pulse broadening time of $\tau_{\rm sc, 10} = 0.09 \pm 0.03~\rm ms$ and a dispersion measure of ${\rm DM} = 1760^{+2.4}_{-1.3}~{\rm pc~cm}^{-3}$. Both of these results are consistent with previous measurements, which suggests that the scattering and dispersion measure of J1745-2900 may be stable on timescales of several years.

We present the results of PSR$\pi$, a large astrometric project targeting radio pulsars using the Very Long Baseline Array (VLBA). From our astrometric database of 60 pulsars, we have obtained parallax-based distance measurements for all but 3, with a parallax precision of typically 40 $\mu$as and approaching 10 $\mu$as in the best cases. Our full sample doubles the number of radio pulsars with a reliable ($\gtrsim$5$\sigma$) model-independent distance constraint. Importantly, many of the newly measured pulsars are well outside the solar neighbourhood, and so PSR$\pi$ brings a near-tenfold increase in the number of pulsars with a reliable model-independent distance at $d>2$ kpc. Using our sample along with previously published results, we show that even the most recent models of the Galactic electron density distribution model contain significant shortcomings, particularly at high Galactic latitudes. When comparing our results to pulsar timing, two of the four millisecond pulsars in our sample exhibit significant discrepancies in the estimates of proper motion obtained by at least one pulsar timing array. With additional VLBI observations to improve the absolute positional accuracy of our reference sources and an expansion of the number of millisecond pulsars, we will be able to extend the comparison of proper motion discrepancies to a larger sample of pulsar reference positions, which will provide a much more sensitive test of the applicability of the solar system ephemerides used for pulsar timing. Finally, we use our large sample to estimate the typical accuracy attainable for differential astrometry with the VLBA when observing pulsars, showing that for sufficiently bright targets observed 8 times over 18 months, a parallax uncertainty of 4 $\mu$as per arcminute of separation between the pulsar and calibrator can be expected.

We briefly review the various proposed scenarios that may lead to nonthermal radio emissions from exoplanetary systems (planetary magnetospheres, magnetosphere-ionosphere and magnetosphere-satellite coupling, and star-planet interactions), and the physical information that can be drawn from their detection. The latter scenario is especially favorable to the production of radio emission above 70\,MHz. We summarize the results of past and recent radio searches, and then discuss FAST characteristics and observation strategy, including synergies. We emphasize the importance of polarization measurements and a high duty-cycle for the very weak targets that radio-exoplanets prove to be.

Gravitational wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. Multi-wavelength observations with next-generation telescopes, including systematic large-area surveys, will provide unprecedented opportunities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales. We highlight ways in which these observations will open a new avenue for multi-messenger studies in the dawn of low-frequency (~ nHz - mHz) GW astronomy.

After a decade of great progress in understanding gas flow into, out of, and through the Milky Way, we are poised to merge observations with simulations to build a comprehensive picture of the multi-scale magnetized interstellar medium (ISM). These insights will also be crucial to four bold initiatives in the 2020s: detecting nanohertz gravitational waves with pulsar timing arrays (PTAs), decoding fast radio bursts (FRBs), cosmic B-mode detection, and imaging the Milky Way's black hole with the Event Horizon Telescope (EHT).

There is an opportunity to advance both solar system and extrasolar planetary studies that does not require the construction of new telescopes or new missions but better use and access to inter-disciplinary data sets. This approach leverages significant investment from NASA and international space agencies in exploring this solar system and using those discoveries as "ground truth" for the study of extrasolar planets. This white paper illustrates the potential, using phase curves and atmospheric modeling as specific examples. A key advance required to realize this potential is to enable seamless discovery and access within and between planetary science and astronomical data sets. Further, seamless data discovery and access also expands the availability of science, allowing researchers and students at a variety of institutions, equipped only with Internet access and a decent computer to conduct cutting-edge research.

Multimessenger observations of the binary neutron star merger GW170817 have enabled the discovery of a diverse array of electromagnetic counterparts to compact binary mergers, including an unambiguous kilonova, a short gamma-ray burst, and a late-time radio jet. Beyond these counterparts, compact binary mergers are additionally predicted to be accompanied by prompt low-frequency radio emission. The successful observation of a prompt radio counterpart would be immensely valuable, but is made difficult by the short delay between the gravitational-wave and prompt electromagnetic signals as well as the poor localization of gravitational-wave sources. Here, we present the first search for prompt radio emission accompanying a gravitational-wave event, targeting the binary black hole merger GW170104 detected by the Advanced LIGO and Virgo gravitational-wave observatories during their second (O2) observing run. Using the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA), we search a $\sim900\,\mathrm{deg}^2$ region for transient radio emission within approximately one hour of GW170104, obtaining an upper limit of $2.5\times10^{41}\,\mathrm{erg}\,\mathrm{s}^{-1}$ on its equivalent isotropic luminosity between 27-84 MHz. We additionally discuss plans to target binary neutron star mergers in Advanced LIGO and Virgo's upcoming O3 observing run.

Advancing the scientific frontier in the search for life in the universe requires support of searches for both biosignatures and technosignatures. A modest budgetary increment can expand the search for life in the universe from primitive to complex life and from the solar neighborhood to the entire Galaxy.

The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models.

Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.

Planetary radar observations offer the potential for probing the properties of characteristics of solid bodies throughout the inner solar system and at least as far as the orbit of Saturn. In addition to the direct scientific value, precise orbital determinations can be obtained from planetary radar observations, which are in turn valuable for mission planning or spacecraft navigation and planetary defense. The next-generation Very Large Array would not have to be equipped with a transmitter to be an important asset in the world's planetary radar infrastructure. Bistatic radar, in which one antenna transmits (e.g., Arecibo or Goldstone) and another receives, are used commonly today, with the Green Bank Telescope (GBT) serving as a receiver. The improved sensitivity of the ngVLA relative to the GBT would improve the signal-to-noise ratios on many targets and increase the accessible volume specifically for asteroids. Goldstone-ngVLA bistatic observations would have the potential of rivaling the sensitivity of Arecibo, but with much wider sky access.

There is a long history of radio telescopes being used to augment the radio antennas regularly used to conduct telemetry, tracking, and command of deep space spacecraft. Radio telescopes are particularly valuable during short-duration mission critical events, such as planetary landings, or when a mission lifetime itself is short, such as a probe into a giant planet's atmosphere. By virtue of its high sensitivity and frequency coverage, the next-generation Very Large Array would be a powerful addition to regular spacecraft ground systems. Further, the science focus of many of these deep-space missions provides a "ground truth" in the solar system that complements other aspects of the ngVLA's science case, such as the formation of planets in proto-planetary disks.

The science case and associated science requirements for a next-generation Very Large Array (ngVLA) are described, highlighting the five key science goals developed out of a community-driven vision of the highest scientific priorities in the next decade. Building on the superb cm observing conditions and existing infrastructure of the VLA site in the U.S. Southwest, the ngVLA is envisaged to be an interferometric array with more than 10 times the sensitivity and spatial resolution of the current VLA and ALMA, operating at frequencies spanning $\sim1.2 - 116$\,GHz with extended baselines reaching across North America. The ngVLA will be optimized for observations at wavelengths between the exquisite performance of ALMA at submm wavelengths, and the future SKA-1 at decimeter to meter wavelengths, thus lending itself to be highly complementary with these facilities. The ngVLA will be the only facility in the world that can tackle a broad range of outstanding scientific questions in modern astronomy by simultaneously delivering the capability to: (1) unveil the formation of Solar System analogues; (2) probe the initial conditions for planetary systems and life with astrochemistry; (3) characterize the assembly, structure, and evolution of galaxies from the first billion years to the present; (4) use pulsars in the Galactic center as fundamental tests of gravity; and (5) understand the formation and evolution of stellar and supermassive blackholes in the era of multi-messenger astronomy.

Gravitational-wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. The ngVLA offers unparalleled capabilities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales, opening a new avenue for multi-messenger studies in the dawn of low-frequency GW astronomy.

Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied.

The frequency dependence of radio pulse arrival times provides a probe of structures in the intervening media. Demorest et al. 2013 was the first to show a short-term (~100-200 days) reduction in the electron content along the line of sight to PSR J1713+0747 in data from 2008 (approximately MJD 54750) based on an apparent dip in the dispersion measure of the pulsar. We report on a similar event in 2016 (approximately MJD 57510), with average residual pulse-arrival times of approximately 3.0,-1.3, and -0.7 microseconds at 820, 1400, and 2300 MHz, respectively. Timing analyses indicate possible departures from the standard nu^-2 dispersive-delay dependence. We discuss and rule out a wide variety of potential interpretations. We find the likeliest scenario to be lensing of the radio emission by some structure in the interstellar medium, which causes multiple frequency-dependent pulse arrival-time delays.

Precision pulsar timing requires optimization against measurement errors and astrophysical variance from the neutron stars themselves and the interstellar medium. We investigate optimization of arrival time precision as a function of radio frequency and bandwidth. We find that increases in bandwidth that reduce the contribution from receiver noise are countered by the strong chromatic dependence of interstellar effects and intrinsic pulse-profile evolution. The resulting optimal frequency range is therefore telescope and pulsar dependent. We demonstrate the results for five pulsars included in current pulsar timing arrays and determine that they are not optimally observed at current center frequencies. For those objects, we find that better choices of total bandwidth as well as center frequency can improve the arrival-time precision. Wideband receivers centered at somewhat higher frequencies with respect to the currently adopted receivers can reduce required overall integration times and provide significant improvements in arrival time uncertainty by a factor of ~sqrt(2) in most cases, assuming a fixed integration time. We also discuss how timing programs can be extended to pulsars with larger dispersion measures through the use of higher-frequency observations.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.

The Large-Aperture Experiment to Detect the Dark Age (LEDA) was designed to detect the predicted O(100)mK sky-averaged absorption of the Cosmic Microwave Background by Hydrogen in the neutral pre- and intergalactic medium just after the cosmological Dark Age. The spectral signature would be associated with emergence of a diffuse Ly$\alpha$ background from starlight during 'Cosmic Dawn'. Recently, Bowman et al. (2018) have reported detection of this predicted absorption feature, with an unexpectedly large amplitude of 530 mK, centered at 78 MHz. Verification of this result by an independent experiment, such as LEDA, is pressing. In this paper, we detail design and characterization of the LEDA radiometer systems, and a first-generation pipeline that instantiates a signal path model. Sited at the Owens Valley Radio Observatory Long Wavelength Array, LEDA systems include the station correlator, five well-separated redundant dual polarization radiometers and backend electronics. The radiometers deliver a 30-85MHz band (16<z<34) and operate as part of the larger interferometric array, for purposes ultimately of in situ calibration. Here, we report on the LEDA system design, calibration approach, and progress in characterization as of January 2016. The LEDA systems are currently being modified to improve performance near 78 MHz in order to verify the purported absorption feature.

We report the first detections of the repeating fast radio burst source FRB 121102 above 5.2 GHz. Observations were performed using the 4$-$8 GHz receiver of the Robert C. Byrd Green Bank Telescope with the Breakthrough Listen digital backend. We present the spectral, temporal and polarization properties of 21 bursts detected within the first 60 minutes of a total 6-hour observations. These observations comprise the highest burst density yet reported in the literature, with 18 bursts being detected in the first 30 minutes. A few bursts clearly show temporal sub-structures with distinct spectral properties. These sub-structures superimpose to provide enhanced peak signal-to-noise ratio at higher trial dispersion measures. Broad features occur in $\sim 1$ GHz wide subbands that typically differ in peak frequency between bursts within the band. Finer-scale structures ($\sim 10-50$ MHz) within these bursts are consistent with that expected from Galactic diffractive interstellar scintillation. The bursts exhibit nearly 100% linear polarization, and a large average rotation measure of 9.359$\pm$0.012 $\times$ 10$^{\rm 4}$ rad m$^{\rm -2}$ (in the observer's frame). No circular polarization was found for any burst. We measure an approximately constant polarization position angle in the 13 brightest bursts. The peak flux densities of the reported bursts have average values (0.2$\pm$0.1 Jy), similar to those seen at lower frequencies ($<3$ GHz), while the average burst widths (0.64$\pm$0.46 ms) are relatively narrower.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.

This paper reports the first OH 18-cm line observation of the first detected interstellar object 1I/2017 U1 (`Oumuamua) using the Green Bank Telescope. We have observed the OH lines at 1665.402 MHz, 1667.359, and 1720.53 MHz frequencies with a spectral resolution of 357 Hz (approximately 0.06 km-s^{-1}). At the time of the observation, `Oumuamua was at topocentric distance and velocity of 1.07 au and 63.4 km-s^{-1}, respectively, or at heliocentric distance and velocity of 1.8 au and 39 km-s^{-1}, respectively. Based on a detailed data reduction and an analogy-based inversion, our final results confirm the asteroidal origin of `Oumuamua (as discussed in Meech et al., 2017) with an upper bound of OH production of Q[OH] < 0.17 x 10^{28} s^{-1}.

This paper presents Very Long Baseline Array (VLBA) observations of 13 double-peaked [O III] emission-line type-2 Active Galactic Nuclei (AGNs) at redshifts 0.06 < z < 0.41 (with a median redshift of z~0.15) identified in the Sloan Digital Sky Survey. Such double-peaked emission-line objects may result from jets or outflows from the central engine or from a dual AGN. The VLBA provides an angular resolution of <~10 pc at the distance of many of these galaxies, sufficient to resolve the radio emission from extremely close dual AGNs and to contribute to understanding the origin of double-peaked [O III] emission lines. Of the 13 galaxies observed at 3.6 cm (8.4 GHz), we detect six at a 1\sigma\ sensitivity level of ~0.15 mJy/beam, two of which show clear jet structures on scales ranging from a few milliarcseconds to tens of milliarcseconds (corresponding to a few pc to tens of pc at a median redshift of 0.15). We suggest that radio-loud double-peaked emission-line type-2 AGNs may be indicative of jet produced structures, but a larger sample of double-peaked [O III] AGNs with high angular resolution radio observations will be required to confirm this suggestion.

Radio interferometers have the ability to precisely localize and better characterize the properties of sources. This ability is having a powerful impact on the study of fast radio transients, where a few milliseconds of data is enough to pinpoint a source at cosmological distances. However, recording interferometric data at millisecond cadence produces a terabyte-per-hour data stream that strains networks, computing systems, and archives. This challenge mirrors that of other domains of science, where the science scope is limited by the computational architecture as much as the physical processes at play. Here, we present a solution to this problem in the context of radio transients: realfast, a commensal, fast transient search system at the Jansky Very Large Array. Realfast uses a novel architecture to distribute fast-sampled interferometric data to a 32-node, 64-GPU cluster for real-time imaging and transient detection. By detecting transients in situ, we can trigger the recording of data for those rare, brief instants when the event occurs and reduce the recorded data volume by a factor of 1000. This makes it possible to commensally search a data stream that would otherwise be impossible to record. This system will search for millisecond transients in more than 1000 hours of data per year, potentially localizing several Fast Radio Bursts, pulsars, and other sources of impulsive radio emission. We describe the science scope for realfast, the system design, expected outcomes, and ways real-time analysis can help in other fields of astrophysics.

Fast radio bursts (FRBs) are millisecond-duration, extragalactic radio flashes of unknown physical origin. FRB 121102, the only known repeating FRB source, has been localized to a star-forming region in a dwarf galaxy at redshift z = 0.193, and is spatially coincident with a compact, persistent radio source. The origin of the bursts, the nature of the persistent source, and the properties of the local environment are still debated. Here we present bursts that show ~100% linearly polarized emission at a very high and variable Faraday rotation measure in the source frame: RM_src = +1.46 x 10^5 rad m^-2 and +1.33 x 10^5 rad m^-2 at epochs separated by 7 months, in addition to narrow (< 30 mus) temporal structure. The large and variable rotation measure demonstrates that FRB 121102 is in an extreme and dynamic magneto-ionic environment, while the short burst durations argue for a neutron star origin. Such large rotation measures have, until now, only been observed in the vicinities of massive black holes (M_BH > 10^4 MSun). Indeed, the properties of the persistent radio source are compatible with those of a low-luminosity, accreting massive black hole. The bursts may thus come from a neutron star in such an environment. However, the observed properties may also be explainable in other models, such as a highly magnetized wind nebula or supernova remnant surrounding a young neutron star.

We have conducted a search for radio emission consistent with an artificial source targeting 1I/'Oumuamua with the Robert C. Byrd Green Bank Telescope (GBT) between 1.1 and 11.6 GHz. We searched the data for narrowband signals and found none. Given the close proximity to this interstellar object, we can place limits to putative transmitters with extremely low power (0.08 W).

In-spiraling supermassive black holes should emit gravitational waves, which would produce characteristic distortions in the time of arrival residuals from millisecond pulsars. Multiple national and regional consortia have constructed pulsar timing arrays by precise timing of different sets of millisecond pulsars. An essential aspect of precision timing is the transfer of the times of arrival to a (quasi-)inertial frame, conventionally the solar system barycenter. The barycenter is determined from the knowledge of the planetary masses and orbits, which has been refined over the past 50 years by multiple spacecraft. Within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), uncertainties on the solar system barycenter are emerging as an important element of the NANOGrav noise budget. We describe what is known about the solar system barycenter, touch upon how uncertainties in it affect gravitational wave studies with pulsar timing arrays, and consider future trends in spacecraft navigation.

Gravitational wave (GW) detection with pulsar timing arrays (PTAs) requires accurate noise characterization. The noise of our Galactic-scale GW detector has been systematically evaluated by the Noise Budget and Interstellar Medium Mitigation working groups within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Intrinsically, individual radio millisecond pulsars (MSPs) used by NANOGrav can have some degree of achromatic red spin noise, as well as white noise due to pulse phase jitter. Along any given line-of-sight, the ionized interstellar medium contributes chromatic noise through dispersion measure (DM) variations, interstellar scintillation, and scattering. These effects contain both red and white components. In the future, with wideband receivers, the effects of frequency-dependent DM will become important. Having anticipated and measured these diverse sources of detector noise, the NANOGrav PTA remains well-poised to detect low-frequency GWs.

This document describes some of the fundamental astrophysical problems that require observing capabilities at millimeter- and centimeter wavelengths well beyond those of existing, or already planned, telescopes. The results summarized in this report follow a solicitation from the National Radio Astronomy Observatory to develop key science cases for a future U. S.-led radio telescope, the "next generation Very Large Array" (ngVLA). The ngVLA will have roughly 10 times the collecting area of the Jansky VLA, operate at frequencies from 1 GHz to 116 GHz with up to 20 GHz of bandwidth, possess a compact core for high surface-brightness sensitivity, and extended baselines of at least hundreds of kilometers and ultimately across the continent to provide high-resolution imaging. The ngVLA builds on the scientific and technical legacy of the Jansky VLA and ALMA, and will be designed to provide the next leap forward in our understanding of planets, galaxies, and black holes.

We have conducted the most sensitive low frequency (below 100 MHz) search to date for prompt, low-frequency radio emission associated with short-duration gamma-ray bursts (GRBs), using the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA). The OVRO-LWA's nearly full-hemisphere field-of-view ($\sim20$,$000$ square degrees) allows us to search for low-frequency (sub-$100$ MHz) counterparts for a large sample of the subset of GRB events for which prompt radio emission has been predicted. Following the detection of short GRB 170112A by Swift, we used all-sky OVRO-LWA images spanning one hour prior to and two hours following the GRB event to search for a transient source coincident with the position of GRB 170112A. We detect no transient source, with our most constraining $1\sigma$ flux density limit of $650~\text{mJy}$ for frequencies spanning $27~\text{MHz}-84~\text{MHz}$. We place constraints on a number of models predicting prompt, low-frequency radio emission accompanying short GRBs and their potential binary neutron star merger progenitors, and place an upper limit of $L_\text{radio}/L_\gamma \lesssim 7\times10^{-16}$ on the fraction of energy released in the prompt radio emission. These observations serve as a pilot effort for a program targeting a wider sample of both short and long GRBs with the OVRO-LWA, including bursts with confirmed redshift measurements which are critical to placing the most constraining limits on prompt radio emission models, as well as a program for the follow-up of gravitational wave compact binary coalescence events detected by advanced LIGO and Virgo.

A host of new low-frequency radio telescopes seek to measure the 21-cm transition of neutral hydrogen from the early universe. These telescopes have the potential to directly probe star and galaxy formation at redshifts $20 \gtrsim z \gtrsim 7$, but are limited by the dynamic range they can achieve against foreground sources of low-frequency radio emission. Consequently, there is a growing demand for modern, high-fidelity maps of the sky at frequencies below 200 MHz for use in foreground modeling and removal. We describe a new widefield imaging technique for drift-scanning interferometers, Tikhonov-regularized $m$-mode analysis imaging. This technique constructs images of the entire sky in a single synthesis imaging step with exact treatment of widefield effects. We describe how the CLEAN algorithm can be adapted to deconvolve maps generated by $m$-mode analysis imaging. We demonstrate Tikhonov-regularized $m$-mode analysis imaging using the Owens Valley Long Wavelength Array (OVRO-LWA) by generating 8 new maps of the sky north of $\delta=-30^\circ$ with 15 arcmin angular resolution, at frequencies evenly spaced between 36.528 MHz and 73.152 MHz, and $\sim$800 mJy/beam thermal noise. These maps are a 10-fold improvement in angular resolution over existing full-sky maps at comparable frequencies, which have angular resolutions $\ge 2^\circ$. Each map is constructed exclusively from interferometric observations and does not represent the globally averaged sky brightness. Future improvements will incorporate total power radiometry, improved thermal noise, and improved angular resolution -- due to the planned expansion of the OVRO-LWA to 2.6 km baselines. These maps serve as a first step on the path to the use of more sophisticated foreground filters in 21-cm cosmology incorporating the measured angular and frequency structure of all foreground contaminants.

We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done.

Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion $M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large and poorly constrained concentrations of gas and stars. There are currently no observations of merging SMBHBs--- it is in fact possible that they stall at their final parsec of separation and never merge. While LIGO has detected high frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is inaccessible to ground-based interferometers, but possible with Pulsar Timing Arrays (PTAs). Using data from local galaxies in the 2 Micron All-Sky Survey, together with galaxy merger rates from Illustris, we find that there are on average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries which will never merge. Local unresolved SMBHBs can contribute to GW background anisotropy at a level of $\sim20\%$, and if the GW background can be successfully isolated, GWs from at least one local SMBHB can be detected in 10 years.

Radio Exploration of Planetary Habitability was the fifth in the series of American Astronomical Society's Topical Conference Series. Notable aspects of the conference included the interdisciplinary nature of both the topics and the intellectual breadth of the participants, the diversity of approaches to studying this topic presented by recent discoveries and of the participants themselves, the expanding meaning of the topic of "star-planet interactions," and the expectation of an increasingly statistical approach to the topic. Potential areas of future research include the actual extent to which planetary magnetic fields shield planetary atmospheres; the planetary dynamo process itself, particularly once multiple extrasolar planetary magnetic fields are confirmed; and "planet-star interactions." A major major topic of the conference concerned observational opportunities, highlighted by a number of new or upcoming, specialized observatories to observe exoplanets especially at radio wavelengths. This article summarizes these main points of the conference and expands briefly upon these potential avenues for future investigation. A future meeting on this topic, given the variety of data sets being generated over the next few years, is warranted.

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

Plasma lenses in the host galaxies of fast radio bursts (FRBs) can strongly modulate FRB amplitudes for a wide range of distances, including the $\sim $ Gpc distance of the repeater FRB121102. To produce caustics, the lens' dispersion-measure depth (${\rm DM}_{\ell}$), scale size ($a$), and distance from the source ($d_{\rm sl}$) must satisfy ${\rm DM}_{\ell} d_{\rm sl} / a^2 \gtrsim 0.65~ {\rm pc^2 \ AU^{-2} \ cm^{-3}}$. Caustics produce strong magnifications ($\lesssim 10^2$) on short time scales ($\sim$ hours to days and perhaps shorter) along with narrow, epoch dependent spectral peaks (0.1 to 1~GHz). However, strong suppression also occurs in long-duration ($\sim$ months) troughs. For geometries that produce multiple images, the resulting burst components will arrive differentially by $< 1~\mu$s to tens of ms and they will show different apparent dispersion measures, $\delta{\rm DM}_{\rm apparent} \sim 1$ pc cm$^{-3}$. Arrival time perturbations may mask any underlying periodicity with period $\lesssim 1$ s. When arrival times differ by less than the burst width, interference effects in dynamic spectra are expected. Strong lensing requires source sizes smaller than $({\rm Fresnel~scale)^2} / a$, which can be satisfied by compact objects such as neutron star magnetospheres but not by AGNs. Much of the phenomenology of the repeating fast radio burst source FRB121102 is similar to lensing effects. The overall picture can be tested by obtaining wideband spectra of bursts (from $<1$ to 10 GHz and possibly higher), which can also be used to characterize the plasma environment near FRB sources. A rich variety of phenomena is expected from an ensemble of lenses near the FRB source. We discuss constraints on densities, magnetic fields, and locations of plasma lenses related to requirements for lensing to occur.

The precise localization of the repeating fast radio burst (FRB 121102) has provided the first unambiguous association (chance coincidence probability $p\lesssim3\times10^{-4}$) of an FRB with an optical and persistent radio counterpart. We report on optical imaging and spectroscopy of the counterpart and find that it is an extended ($0.6^{\prime\prime}-0.8^{\prime\prime}$) object displaying prominent Balmer and [OIII] emission lines. Based on the spectrum and emission line ratios, we classify the counterpart as a low-metallicity, star-forming, $m_{r^\prime} = 25.1$ AB mag dwarf galaxy at a redshift of $z=0.19273(8)$, corresponding to a luminosity distance of 972 Mpc. From the angular size, the redshift, and luminosity, we estimate the host galaxy to have a diameter $\lesssim4$ kpc and a stellar mass of $M_*\sim4-7\times 10^{7}\,M_\odot$, assuming a mass-to-light ratio between 2 to 3$\,M_\odot\,L_\odot^{-1}$. Based on the H$\alpha$ flux, we estimate the star formation rate of the host to be $0.4\,M_\odot\,\mathrm{yr^{-1}}$ and a substantial host dispersion measure depth $\lesssim 324\,\mathrm{pc\,cm^{-3}}$. The net dispersion measure contribution of the host galaxy to FRB 121102 is likely to be lower than this value depending on geometrical factors. We show that the persistent radio source at FRB 121102's location reported by Marcote et al (2017) is offset from the galaxy's center of light by $\sim$200 mas and the host galaxy does not show optical signatures for AGN activity. If FRB 121102 is typical of the wider FRB population and if future interferometric localizations preferentially find them in dwarf galaxies with low metallicities and prominent emission lines, they would share such a preference with long gamma ray bursts and superluminous supernovae.

Fast radio bursts are astronomical radio flashes of unknown physical nature with durations of milliseconds. Their dispersive arrival times suggest an extragalactic origin and imply radio luminosities orders of magnitude larger than any other kind of known short-duration radio transient. Thus far, all FRBs have been detected with large single-dish telescopes with arcminute localizations, and attempts to identify their counterparts (source or host galaxy) have relied on contemporaneous variability of field sources or the presence of peculiar field stars or galaxies. These attempts have not resulted in an unambiguous association with a host or multi-wavelength counterpart. Here we report the sub-arcsecond localization of FRB 121102, the only known repeating burst source, using high-time-resolution radio interferometric observations that directly image the bursts themselves. Our precise localization reveals that FRB 121102 originates within 100 mas of a faint 180 uJy persistent radio source with a continuum spectrum that is consistent with non-thermal emission, and a faint (25th magnitude) optical counterpart. The flux density of the persistent radio source varies by tens of percent on day timescales, and very long baseline radio interferometry yields an angular size less than 1.7 mas. Our observations are inconsistent with the fast radio burst having a Galactic origin or its source being located within a prominent star-forming galaxy. Instead, the source appears to be co-located with a low-luminosity active galactic nucleus or a previously unknown type of extragalactic source. [Truncated] If other fast radio bursts have similarly faint radio and optical counterparts, our findings imply that direct sub-arcsecond localizations of FRBs may be the only way to provide reliable associations.

The millisecond-duration radio flashes known as Fast Radio Bursts (FRBs) represent an enigmatic astrophysical phenomenon. Recently, the sub-arcsecond localization (~ 100mas precision) of FRB121102 using the VLA has led to its unambiguous association with persistent radio and optical counterparts, and to the identification of its host galaxy. However, an even more precise localization is needed in order to probe the direct physical relationship between the millisecond bursts themselves and the associated persistent emission. Here we report very-long-baseline radio interferometric observations using the European VLBI Network and the 305-m Arecibo telescope, which simultaneously detect both the bursts and the persistent radio emission at milliarcsecond angular scales and show that they are co-located to within a projected linear separation of < 40pc (< 12mas angular separation, at 95% confidence). We detect consistent angular broadening of the bursts and persistent radio source (~ 2-4mas at 1.7GHz), which are both similar to the expected Milky Way scattering contribution. The persistent radio source has a projected size constrained to be < 0.7pc (< 0.2mas angular extent at 5.0GHz) and a lower limit for the brightness temperature of T_b > 5 x 10^7K. Together, these observations provide strong evidence for a direct physical link between FRB121102 and the compact persistent radio source. We argue that a burst source associated with a low-luminosity active galactic nucleus or a young neutron star energizing a supernova remnant are the two scenarios for FRB121102 that best match the observed data.

We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the 9-year NANOGrav data release and constrain the sources of these variations. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 13 pulsars show both effects. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized ISM features and, in one case, an upper limit to the size of the dispersing region of 13.2 AU. Finally, five pulsars show very nearly quadratic structure functions, which could be indicative of an underlying Kolmogorov medium. Four pulsars show roughly Kolmogorov structure functions and another four show structure functions less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and show that the presence of other trends in the data is the most likely cause.

The Galactic center has some of the highest stellar densities in the Galaxy and a range of interstellar scattering properties that may aid in the detection of new radio-selected transient events. Here we describe a search for radio transients in the Galactic center using over 200 hours of archival data from the Very Large Array (VLA) at 5 and 8.4 GHz. Every observation of SgrA* from 1985$-$2005 has been searched using an automated processing and detection pipeline sensitive to transients with timescales between 30 seconds and five minutes with a typical detection threshold of $\sim$100 mJy. Eight possible candidates pass tests to filter false-positives from radio-frequency interference, calibration errors, and imaging artifacts. Two events are identified as promising candidates based on the smoothness of their light curves. Despite the high quality of their light curves, these detections remain suspect due to evidence of incomplete subtraction of the complex structure in the Galactic center, and apparent contingency of one detection on reduction routines. Events of this intensity ($\sim$100 mJy) and duration ($\sim$100 s) are not obviously associated with known astrophysical sources, and no counterparts are found in data at other wavelengths. We consider potential sources, including Galactic center pulsars, dwarf stars, sources like GCRT J1745-3009, and bursts from X-ray binaries. None can fully explain the observed transients, suggesting either a new astrophysical source or a subtle imaging artifact. More sensitive multiwavelength studies are necessary to characterize these events which, if real, occur with a rate of $14^{+32}_{-12}~{\rm hr}^{-1}\,{\rm deg}^{-2}$ in the Galactic center.

Mergers of binary neutron stars and black hole-neutron star binaries produce gravitational-wave (GW) emission and outflows with significant kinetic energies. These outflows result in radio emissions through synchrotron radiation. We explore the detectability of these synchrotron generated radio signals by follow-up observations of GW merger events lacking a detection of electromagnetic counterparts in other wavelengths. We model radio light curves arising from (i) sub-relativistic merger ejecta and (ii) ultra-relativistic jets. The former produces radio remnants on timescales of a few years and the latter produces $\gamma$-ray bursts in the direction of the jet and orphan-radio afterglows extending over wider angles on timescales of weeks. Based on the derived light curves, we suggest an optimized survey at $1.4$ GHz with five epochs separated by a logarithmic time interval. We estimate the detectability of the radio counterparts of simulated GW-merger events to be detected by advanced LIGO and Virgo by current and future radio facilities. The detectable distances for these GW merger events could be as high as 1 Gpc. $20$--$60\%$ of the long-lasting radio remnants will be detectable in the case of the moderate kinetic energy of $3\cdot 10^{50}$ erg and a circum-merger density of $0.1 {\rm cm^{-3}}$ or larger, while $5$--$20\%$ of the orphan radio afterglows with kinetic energy of $10^{48}$ erg will be detectable. The detection likelihood increases if one focuses on the well-localizable GW events. We discuss the background noise due to radio fluxes of host galaxies and false positives arising from extragalactic radio transients and variable Active Galactic Nuclei and we show that the quiet radio transient sky is of great advantage when searching for the radio counterparts.

Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way.

Model-independent distance constraints to binary millisecond pulsars (MSPs) are of great value to both the timing observations of the radio pulsars, and multiwavelength observations of their companion stars. Very Long Baseline Interferometry (VLBI) astrometry can be employed to provide these model-independent distances with very high precision via the detection of annual geometric parallax. Using the Very Long Baseline Array, we have observed two binary millisecond pulsars, PSR J1022+1001 and J2145-0750, over a two-year period and measured their distances to be 700 +14 -10 pc and 613 +16 -14 pc respectively. We use the well-calibrated distance in conjunction with revised analysis of optical photometry to tightly constrain the nature of their massive (M ~ 0.85 Msun) white dwarf companions. Finally, we show that several measurements of their parallax and proper motion of PSR J1022+1001 and PSR J2145-0750 obtained by pulsar timing array projects are incorrect, differing from the more precise VLBI values by up to 5 sigma. We investigate possible causes for the discrepancy, and find that imperfect modeling of the solar wind is a likely candidate for the timing model errors given the low ecliptic latitude of these two pulsars.

Detection of the global HI 21 cm signal from Cosmic Dawn and Epoch of Reionization is the key science driver for several ongoing ground-based and future ground/space-based experiments. The crucial spectral features in the global 21 cm signal (turning points) occur at low radio frequencies <100 MHz. In addition to the human-generated RFI, Earth's ionosphere drastically corrupts low-frequency radio observations from the ground. In this paper, we examine the effects of time-varying ionospheric refraction, absorption and thermal emission at these low radio frequencies and their combined effect on any ground-based global 21 cm experiment. It should be noted that this is the first study of the effect of a dynamic ionosphere on global 21 cm experiments. The fluctuations in the ionosphere are influenced by solar activity with flicker noise characteristics. The same characteristics are reflected in the ionospheric corruption to any radio signal passing through the ionosphere. As a result, any ground based observations of the faint global 21 cm signal are corrupted by flicker noise (or "$1/f$" noise, where "$f$" is the dynamical frequency) which scales as $\nu^{-2}$ (where $\nu$ is the frequency of observation) in the presence of a bright galactic foreground ($\propto \nu^{-s}$, where $s$ is radio spectral index). Hence, the calibration of the ionosphere for any such experiment is critical. Any attempt to calibrate the ionospheric effects will be subject to the inaccuracies in the current ionospheric measurements using GPS ionospheric measurements, riometer measurements, ionospheric soundings, etc. Even considering an optimistic improvement in the accuracy of GPS-TEC (Total Electron Content) measurements, we conclude that the detection of the global 21 cm signal below 100 MHz is best done from above the Earth's atmosphere in orbit of the Moon.

Decade-long timing observations of arrays of millisecond pulsars have placed highly constraining upper limits on the amplitude of the nanohertz gravitational-wave stochastic signal from the mergers of supermassive black-hole binaries ($\sim 10^{-15}$ strain at $f = 1/\mathrm{yr}$). These limits suggest that binary merger rates have been overestimated, or that environmental influences from nuclear gas or stars accelerate orbital decay, reducing the gravitational-wave signal at the lowest, most sensitive frequencies. This prompts the question whether nanohertz gravitational waves are likely to be detected in the near future. In this letter, we answer this question quantitatively using simple statistical estimates, deriving the range of true signal amplitudes that are compatible with current upper limits, and computing expected detection probabilities as a function of observation time. We conclude that small arrays consisting of the pulsars with the least timing noise, which yield the tightest upper limits, have discouraging prospects of making a detection in the next two decades. By contrast, we find large arrays are crucial to detection because the quadrupolar spatial correlations induced by gravitational waves can be well sampled by many pulsar pairs. Indeed, timing programs which monitor a large and expanding set of pulsars have an $\sim 80\%$ probability of detecting gravitational waves within the next ten years, under assumptions on merger rates and environmental influences ranging from optimistic to conservative. Even in the extreme case where $90\%$ of binaries stall before merger and environmental coupling effects diminish low-frequency gravitational-wave power, detection is delayed by at most a few years.

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.

We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the NANOGrav pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ~800-MHz wide frequency bands imply dramatic changes in scintillation bandwidth across the bandpass, and a stretching routine has been included to account for this scaling. For most of the 10 pulsars for which the scaling has been measured, the bandwidths scale with frequency less steeply than expected for a Kolmogorov medium. We find estimated scattering delay values that vary with time by up to an order of magnitude. The mean measured scattering delays are similar to previously published values and slightly higher than predicted by interstellar medium models. We investigate the possibility of increasing the timing precision by mitigating timing errors introduced by the scattering delays. For most of the pulsars, the uncertainty in the time of arrival of a single timing point is much larger than the maximum variation of the scattering delay, suggesting that diffractive scintillation remains only a negligible part of their noise budget.

We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.

Most black holes (BHs) will absorb a neutron star (NS) companion fully intact, without tidal disruption, suggesting the pair will remain dark to telescopes. Even without tidal disruption, electromagnetic luminosity is generated from the battery phase of the binary when the BH interacts with the NS magnetic field. Originally the luminosity was expected in high-energy X-rays or gamma-rays, however we conjecture that some of the battery power is emitted in the radio bandwidth. While the luminosity and timescale are suggestive of fast radio bursts (FRBs; millisecond-scale radio transients) NS--BH coalescence rates are too low to make these a primary FRB source. Instead, we propose the transients form a FRB sub-population, distinguishable by a double peak with a precursor. The rapid ramp-up in luminosity manifests as a precursor to the burst which is $20\%-80\%$ as luminous, given 0.5~ms timing resolution. The main burst is from the peak luminosity before merger. The post-merger burst follows from the NS magnetic field migration to the BH, causing a shock. NS--BH pairs are especially desirable for ground-based gravitational wave (GW) observatories since the pair might not otherwise be detected, with electromagnetic counterparts greatly augmenting the scientific leverage beyond the GW signal. Valuably, the electromagnetic signal can break degeneracies in the parameters encoded in the GW as well as probe the NS magnetic field strength, yielding insights into open problems in NS magnetic field decay.

We present initial pulsar results from the first station of the Long Wavelength Array (LWA1) obtained during the commissioning period of LWA1 and early science results. We present detections of periodic emission from 44 previously known pulsars, including 3 millisecond pulsars (MSPs). The effects of the interstellar medium on pulsar emission are significantly enhanced at the low frequencies of the LWA1 band (10--88 MHz), making LWA1 a very sensitive instrument for characterizing changes in dispersion measures (DM) and other effects from the interstellar medium. Pulsars also often have significant evolution in their pulse profile at low frequency and a break in their spectral index. We report DM measurements for 44 pulsars, mean flux density measurements for 36 pulsars, and multi-frequency component spacing and widths for 15 pulsars with more than one profile component. For 27 pulsars, we report spectral index measurements within our frequency range. We also introduce the LWA1 Pulsar Data Archive, which stores reduced data products from LWA1 pulsar observations. Reduced data products for the observations presented here can be found on the archive. Reduced data products from future LWA1 pulsar observations will also be made available through the archive.

The U.S. Virtual Astronomical Observatory was a software infrastructure and development project designed both to begin the establishment of an operational Virtual Observatory (VO) and to provide the U.S. coordination with the international VO effort. The concept of the VO is to provide the means by which an astronomer is able to discover, access, and process data seamlessly, regardless of its physical location. This paper describes the origins of the VAO, including the predecessor efforts within the U.S. National Virtual Observatory, and summarizes its main accomplishments. These accomplishments include the development of both scripting toolkits that allow scientists to incorporate VO data directly into their reduction and analysis environments and high-level science applications for data discovery, integration, analysis, and catalog cross-comparison. Working with the international community, and based on the experience from the software development, the VAO was a major contributor to international standards within the International Virtual Observatory Alliance. The VAO also demonstrated how an operational virtual observatory could be deployed, providing a robust operational environment in which VO services worldwide were routinely checked for aliveness and compliance with international standards. Finally, the VAO engaged in community outreach, developing a comprehensive web site with on-line tutorials, announcements, links to both U.S. and internationally developed tools and services, and exhibits and hands-on training .... All digital products of the VAO Project, including software, documentation, and tutorials, are stored in a repository for community access. The enduring legacy of the VAO is an increasing expectation that new telescopes and facilities incorporate VO capabilities during the design of their data management systems.

We use MERLIN, VLA and VLBA observations of Galactic \HI absorption towards 3C~138 to estimate the structure function of the \HI opacity fluctuations at AU scales. Using Monte Carlo simulations, we show that there is likely to be a significant bias in the estimated structure function at signal-to-noise ratios characteristic of our observations, if the structure function is constructed in the manner most commonly used in the literature. We develop a new estimator that is free from this bias and use it to estimate the true underlying structure function slope on length scales ranging $5$ to $40$~AU. From a power law fit to the structure function, we derive a slope of $0.81^{+0.14}_{-0.13}$, i.e. similar to the value observed at parsec scales. The estimated upper limit for the amplitude of the structure function is also consistent with the measurements carried out at parsec scales. Our measurements are hence consistent with the \HI opacity fluctuation in the Galaxy being characterized by a power law structure function over length scales that span six orders of magnitude. This result implies that the dissipation scale has to be smaller than a few AU if the fluctuations are produced by turbulence. This inferred smaller dissipation scale implies that the dissipation occurs either in (i) regions with densities $\gtrsim 10^3 $cm$^-3$ (i.e. similar to that inferred for "tiny scale" atomic clouds or (ii) regions with a mix of ionized and atomic gas (i.e. the observed structure in the atomic gas has a magneto-hydrodynamic origin).

Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.

The discovery and timing of radio pulsars within the Galactic centre is a fundamental aspect of the SKA Science Case, responding to the topic of "Strong Field Tests of Gravity with Pulsars and Black Holes" (Kramer et al. 2004; Cordes et al. 2004). Pulsars have in many ways proven to be excellent tools for testing the General theory of Relativity and alternative gravity theories (see Wex (2014) for a recent review). Timing a pulsar in orbit around a companion, provides a unique way of probing the relativistic dynamics and spacetime of such a system. The strictest tests of gravity, in strong field conditions, are expected to come from a pulsar orbiting a black hole. In this sense, a pulsar in a close orbit ($P_{\rm orb}$ < 1 yr) around our nearest supermassive black hole candidate, Sagittarius A* - at a distance of ~8.3 kpc in the Galactic centre (Gillessen et al. 2009a) - would be the ideal tool. Given the size of the orbit and the relativistic effects associated with it, even a slowly spinning pulsar would allow the black hole spacetime to be explored in great detail (Liu et al. 2012). For example, measurement of the frame dragging caused by the rotation of the supermassive black hole, would allow a test of the "cosmic censorship conjecture." The "no-hair theorem" can be tested by measuring the quadrupole moment of the black hole. These are two of the prime examples for the fundamental studies of gravity one could do with a pulsar around Sagittarius A*. As will be shown here, SKA1-MID and ultimately the SKA will provide the opportunity to begin to find and time the pulsars in this extreme environment.

We report on the first millisecond timescale radio interferometric search for the new class of transient known as fast radio bursts (FRBs). We used the Very Large Array (VLA) for a 166-hour, millisecond imaging campaign to detect and precisely localize an FRB. We observed at 1.4 GHz and produced visibilities with 5 ms time resolution over 256 MHz of bandwidth. Dedispersed images were searched for transients with dispersion measures from 0 to 3000 pc/cm3. No transients were detected in observations of high Galactic latitude fields taken from September 2013 though October 2014. Observations of a known pulsar show that images typically had a thermal-noise limited sensitivity of 120 mJy/beam (8 sigma; Stokes I) in 5 ms and could detect and localize transients over a wide field of view. Our nondetection limits the FRB rate to less than 7e4/sky/day (95% confidence) above a fluence limit of 1.2 Jy-ms. Assuming a Euclidean flux distribution, the VLA rate limit is inconsistent with the published rate of Thornton et al. We recalculate previously published rates with a homogeneous consideration of the effects of primary beam attenuation, dispersion, pulse width, and sky brightness. This revises the FRB rate downward and shows that the VLA observations had a roughly 60% chance of detecting a typical FRB and that a 95% confidence constraint would require roughly 500 hours of similar VLA observing. Our survey also limits the repetition rate of an FRB to 2 times less than any known repeating millisecond radio transient.

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

The period between the creation of the cosmic microwave background at a redshift of ~1000 and the formation of the first stars and black holes that re-ionize the intergalactic medium at redshifts of 10-20 is currently unobservable. The baryonic component of the universe during this period is almost entirely neutral hydrogen, which falls into local regions of higher dark matter density. This seeds the formation of large-scale structures including the cosmic web that we see today in the filamentary distribution of galaxies and clusters of galaxies. The only detectable signal from these dark ages is the 21-cm spectral line of hydrogen, redshifted down to frequencies of approximately 10-100 MHz. Space-based observations of this signal will allow us to determine the formation epoch and physics of the first sources of ionizing radiation, and potentially detect evidence for the decay of dark matter particles. JPL is developing deployable low frequency antenna and receiver prototypes to enable both all-sky spectral measurements of neutral hydrogen and ultimately to map the spatial distribution of the signal as a function of redshift. Such observations must be done from space because of Earth's ionosphere and ubiquitous radio interference. A specific application of these technologies is the Dark Ages Radio Explorer (DARE) mission. This small Explorer class mission is designed to measure the sky-averaged hydrogen signal from the shielded region above the far side of the Moon. These data will complement ground-based radio observations of the final stages of intergalactic re-ionization at higher frequencies. DARE will also serve as a scientific percursor for space-based interferometry missions to image the distribution of hydrogen during the cosmic dark ages.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

We present the results of a recent re-reduction of the data from the Very Large Array (VLA) Low-frequency Sky Survey (VLSS). We used the VLSS catalog as a sky model to correct the ionospheric distortions in the data and create a new set of sky maps and corresponding catalog at 73.8 MHz. The VLSS Redux (VLSSr) has a resolution of 75 arcsec, and an average map RMS noise level of $\sigma\sim0.1$ Jy beam$^{-1}$. The clean bias is $0.66\times\sigma$, and the theoretical largest angular size is 36 arcmin. Six previously un-imaged fields are included in the VLSSr, which has an unbroken sky coverage over 9.3 sr above an irregular southern boundary. The final catalog includes 92,964 sources. The VLSSr improves upon the original VLSS in a number of areas including imaging of large sources, image sensitivity, and clean bias; however the most critical improvement is the replacement of an inaccurate primary beam correction which caused source flux errors which vary as a function of radius to nearest pointing center in the VLSS.

We present the Chasing the Identification of ASCA Galactic Objects (ChIcAGO) survey, which is designed to identify the unknown X-ray sources discovered during the ASCA Galactic Plane Survey (AGPS). Little is known about most of the AGPS sources, especially those that emit primarily in hard X-rays (2-10 keV) within the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 X-ray flux range. In ChIcAGO, the subarcsecond localization capabilities of Chandra have been combined with a detailed multi-wavelength follow-up program, with the ultimate goal of classifying the >100 unidentified sources in the AGPS. Overall to date, 93 unidentified AGPS sources have been observed with Chandra as part of the ChIcAGO survey. A total of 253 X-ray point sources have been detected in these Chandra observations within 3' of the original ASCA positions. We have identified infrared and optical counterparts to the majority of these sources, using both new observations and catalogs from existing Galactic plane surveys. X-ray and infrared population statistics for the X-ray point sources detected in the Chandra observations reveal that the primary populations of Galactic plane X-ray sources that emit in the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 flux range are active stellar coronae, massive stars with strong stellar winds that are possibly in colliding-wind binaries, X-ray binaries, and magnetars. There is also a fifth population that is still unidentified but, based on its X-ray and infrared properties, likely comprise partly of Galactic sources and partly of active galactic nuclei.

This whitepaper describes how the VLASS could be designed in a manner to allow the identification of candidate dual active galactic nuclei (AGN) at separations <7 kpc. Dual AGN represent a clear marker of two supermassive black holes within an ongoing merger. A dual AGN survey will provide a wealth of studies in structure growth and gravitational-wave science. Radio wavelengths are ideal for identifying close pairs, as disturbed stellar and gaseous material can obscure their presence in optical and shorter wavelengths. With sufficiently high resolution and sensitivity, a large-scale radio imaging survey like the VLASS will uncover many of these systems and provide the means to broadly study the radio properties of candidate dual systems revealed at other wavelengths. We determine that the ideal survey for our purposes will be at as high a resolution as possible, with significantly more science return in A array at L-band or higher, or B array at C-band or higher. We describe a range of potential survey parameters within this document. Based on the analysis outlined in this whitepaper, our ideal survey would create a catalogue of $\gtrsim$100 dual AGN in either: 1) a medium-sensitivity (~1 mJy detection threshold), wide-field (few thousand square degree) survey, or 2) a high-sensitivity (~10 uJy threshold) survey of several hundred square degrees.

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.

We present flux density measurements and pulse profiles for the millisecond pulsar PSR J2145-0750 spanning 37 to 81 MHz using data obtained from the first station of the Long Wavelength Array. These measurements represent the lowest frequency detection of pulsed emission from a millisecond pulsar to date. We find that the pulse profile is similar to that observed at 102 MHz. We also find that the flux density spectrum between ~40 MHz to 5 GHz is suggestive of a break and may be better fit by a model that includes spectral curvature with a rollover around 730 MHz rather than a single power law.

We report the detection and observed characteristics of giant pulses from the Crab Nebula pulsar (B0531+21) in four frequency bands covering 20-84 MHz using the recently-completed Long Wavelength Array Station 1 (LWA1) radio telescope. In 10 hours of observations distributed over a 72-day period in Fall of 2012, 33 giant pulses having peak flux densities between 400 Jy and 2000 Jy were detected. Twenty-two of these pulses were detected simultaneously in channels of 16 MHz bandwidth centered at 44 MHz, 60 MHz, and 76 MHz, including one pulse which was also detected in a channel centered at 28 MHz. We quantify statistics of pulse amplitude and pulse shape characteristics, including pulse broadening. Amplitude statistics are consistent with expectations based on extrapolations from previous work at higher and lower frequencies. Pulse broadening is found to be relatively high, but not significantly greater than expected. We present procedures that have been found to be effective for observing giant pulses in this frequency range.

Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.

A novel concept is presented in this paper for a human mission to the lunar L2 (Lagrange) point that would be a proving ground for future exploration missions to deep space while also overseeing scientifically important investigations. In an L2 halo orbit above the lunar farside, the astronauts aboard the Orion Crew Vehicle would travel 15% farther from Earth than did the Apollo astronauts and spend almost three times longer in deep space. Such a mission would serve as a first step beyond low Earth orbit and prove out operational spaceflight capabilities such as life support, communication, high speed re-entry, and radiation protection prior to more difficult human exploration missions. On this proposed mission, the crew would teleoperate landers and rovers on the unexplored lunar farside, which would obtain samples from the geologically interesting farside and deploy a low radio frequency telescope. Sampling the South Pole-Aitken basin, one of the oldest impact basins in the solar system, is a key science objective of the 2011 Planetary Science Decadal Survey. Observations at low radio frequencies to track the effects of the Universe's first stars/galaxies on the intergalactic medium are a priority of the 2010 Astronomy and Astrophysics Decadal Survey. Such telerobotic oversight would also demonstrate capability for human and robotic cooperation on future, more complex deep space missions such as exploring Mars.

The discovery of a pulsar (PSR) in orbit around a black hole (BH) is expected to provide a superb new probe of relativistic gravity and BH properties. Apart from a precise mass measurement for the BH, one could expect a clean verification of the dragging of space-time caused by the BH spin. In order to measure the quadrupole moment of the BH for testing the no-hair theorem of general relativity (GR), one has to hope for a sufficiently massive BH. In this respect, a PSR orbiting the super-massive BH in the center of our Galaxy would be the ultimate laboratory for gravity tests with PSRs. But even for gravity theories that predict the same properties for BHs as GR, a PSR-BH system would constitute an excellent test system, due to the high grade of asymmetry in the strong field properties of these two components. Here we highlight some of the potential gravity tests that one could expect from different PSR-BH systems, utilizing present and future radio telescopes, like FAST and SKA.

[Abridged] The Moon is a unique platform from and on which to conduct astrophysical measurements. The Lunar University Network for Astrophysics Research (LUNAR) and the Center for Lunar Origins and Evolution (CLOE) teams within the NASA Lunar Science Institute (NLSI) are illustrating how the Moon can be used as a platform to advance important goals in astrophysics. Of relevance to Astrophysics and aligned with NASA strategic goals, all three of the primary research themes articulated by New Worlds, New Horizons in Astronomy & Astrophysics are being addressed by LUNAR and CLOE, namely Probing Cosmic Dawn, Understanding New Worlds, and Physics of the Universe

The detection of radio pulsars within the central few parsecs of the Galaxy would provide a unique probe of the gravitational and magneto-ionic environments in the Galactic Center (GC) and, if close enough to Sgr A*, precise tests of general relativity in the strong-field regime. While it is difficult to find pulsars at radio wavelengths because of interstellar scattering, the payoff from detailed timing of pulsars in the GC warrants a concerted effort. To motivate pulsar surveys and help define search parameters for them, we constrain the pulsar number and spatial distribution using a wide range of multiwavelength measurements. These include the five known radio pulsars within 15 arcmin of Sgr A*, radio and gamma-ray measurements of diffuse emission, non-detections in high frequency pulsar surveys of the central parsec, a catalog of radio point sources from an imaging survey, infrared observations of massive star populations in the central few parsecs, candidate pulsar wind nebulae in the inner 20 pc and estimates of the core-collapse supernova rate based on X-ray measurements. We find that under current observational constraints, the inner parsec of the Galaxy could harbor as many as ~10^3 active radio pulsars that are beamed towards Earth. Such a large population would distort the low-frequency measurements of both the intrinsic spectrum of Sgr A* and the free-free absorption along the line of sight of Sgr A*.

The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities of ASKAP to enable the discovery and investigation of variable and transient phenomena from the local to the cosmological, including flare stars, intermittent pulsars, X-ray binaries, magnetars, extreme scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. In addition, it will allow us to probe unexplored regions of parameter space where new classes of transient sources may be detected. In this paper we review the known radio transient and variable populations and the current results from blind radio surveys. We outline a comprehensive program based on a multi-tiered survey strategy to characterise the radio transient sky through detection and monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and greater. We also present an analysis of the expected source populations that we will be able to detect with VAST.

The Virtual Observatory (VO) is realizing global electronic integration of astronomy data. One of the long-term goals of the U.S. VO project, the Virtual Astronomical Observatory (VAO), is development of services and protocols that respond to the growing size and complexity of astronomy data sets. This paper describes how VAO staff are active in such development efforts, especially in innovative strategies and techniques that recognize the limited operating budgets likely available to astronomers even as demand increases. The project has a program of professional outreach whereby new services and protocols are evaluated.

The U.S. Virtual Astronomical Observatory (VAO; http://www.us-vao.org/) has been in operation since May 2010. Its goal is to enable new science through efficient integration of distributed multi-wavelength data. This paper describes the management and organization of the VAO, and emphasizes the techniques used to ensure efficiency in a distributed organization. Management methods include using an annual program plan as the basis for establishing contracts with member organizations, regular communication, and monitoring of processes.

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

The structure function of opacity fluctuations is a useful statistical tool to study tiny scale structures of neutral hydrogen. Here we present high resolution observation of HI absorption towards 3C 138, and estimate the structure function of opacity fluctuations from the combined VLA, MERLIN and VLBA data. The angular scales probed in this work are ~ 10-200 milliarcsec (about 5-100 AU). The structure function in this range is found to be well represented by a power law S_tau(x) ~ x^{beta} with index beta ~ 0.33 +/- 0.07 corresponding to a power spectrum P_tau(U) ~ U^{-2.33}. This is slightly shallower than the earlier reported power law index of ~ 2.5-3.0 at ~ 1000 AU to few pc scales. The amplitude of the derived structure function is a factor of ~ 20-60 times higher than the extrapolated amplitude from observation of Cas A at larger scales. On the other hand, extrapolating the AU scale structure function for 3C 138 predicts the observed structure function for Cas A at the pc scale correctly. These results clearly establish that the atomic gas has significantly more structures in AU scales than expected from earlier pc scale observations. Some plausible reasons are identified and discussed here to explain these results. The observational evidence of a shallower slope and the presence of rich small scale structures may have implications for the current understanding of the interstellar turbulence.

We report the results of a low frequency radio variability and slow transient search using archival observations from the Very Long Array. We selected six 325 MHz radio observations from the spring of 2006, each centered on the Spitzer-Space-Telescope Wide-area Infrared Extragalactic Survey (SWIRE) Deep Field: 1046+59. Observations were spaced between 1 day to 3 months, with a typical single-epoch peak flux sensitivity below 0.2 \mjb near the field pointing center. We describe the observation parameters, data post-processing, and search methodology used to identify variable and transient emission. Our search revealed multiple variable sources and the presence of one, day-scale transient event with no apparent astronomical counterpart. This detection implies a transient rate of 1$\pm$1 event per 6.5 $\deg^2$ per 72 observing hours in the direction of 1046+59 and an isotropic transient surface density $\Sigma = 0.12 \deg^{-2}$ at 95% confidence for sources with average peak flux density higher than 2.1 mJy over 12 hr.

We have used a relatively long, contiguous VHF observation of a bright cosmic radio source (Cygnus A) with the Very Large Array (VLA) to demonstrate the capability of this instrument to study the ionosphere. This interferometer, and others like it, can observe ionospheric total electron content (TEC) fluctuations on a much wider range of scales than is possible with many other instruments. We have shown that with a bright source, the VLA can measure differential TEC values between pairs of antennas (delta-TEC) with an precision of 0.0003 TECU. Here, we detail the data reduction and processing techniques used to achieve this level of precision. In addition, we demonstrate techniques for exploiting these high-precision delta-TEC measurements to compute the TEC gradient observed by the array as well as small-scale fluctuations within the TEC gradient surface. A companion paper details specialized spectral analysis techniques used to characterize the properties of wave-like fluctuations within this data.

We have used a relatively long, contiguous VHF observation of a bright cosmic radio source (Cygnus A) with the Very Large Array (VLA) through the nighttime, midlatitude ionosphere to demonstrate the phenomena observable with this instrument. In a companion paper, we showed that the VLA can detect fluctuations in total electron content (TEC) with amplitudes of <0.001 TECU and can measure TEC gradients with a precision of about 0.0002 TECU/km. We detail two complementary techniques for producing spectral analysis of these TEC gradient measurements. The first is able to track individual waves with wavelengths of about half the size of the array (~20 km) or more. This technique was successful in detecting and characterizing many medium-scale traveling ionospheric disturbances (MSTIDs) seen intermittently throughout the night and has been partially validated using concurrent GPS measurements. Smaller waves are also seen with this technique at nearly all times, many of which move in similar directions as the detected MSTIDs. The second technique allows for the detection and statistical description of the properties of groups of waves moving in similar directions with wavelengths as small as 5 km. Combining the results of both spectral techniques, we found a class of intermediate and small scale waves which are likely the quasi-periodic (QP) echoes that have been observed to occur within sporadic-E (Es) layers. We find two distinct populations of these waves. The members of one population are coincident in time with MSTIDs and are consistent with being generated within Es layers by the E-F coupling instability. The other population seems more influenced by the neutral wind, similar to the predominant types of QP echoes found by the Sporadic-E Experiments over Kyushu (Fukao et al. 1998; Yamamoto et al. 2005).

The discovery of radio pulsars in compact orbits around Sgr A* would allow an unprecedented and detailed investigation of the spacetime of the supermassive black hole. This paper shows that pulsar timing, including that of a single pulsar, has the potential to provide novel tests of general relativity, in particular its cosmic censorship conjecture and no-hair theorem for rotating black holes. These experiments can be performed by timing observations with 100 micro-second precision, achievable with the Square Kilometre Array for a normal pulsar at frequency above 15 GHz. Based on the standard pulsar timing technique, we develop a method that allows the determination of the mass, spin, and quadrupole moment of Sgr A*, and provides a consistent covariance analysis of the measurement errors. Furthermore, we test this method in detailed mock data simulations. It seems likely that only for orbital periods below ~0.3 yr is there the possibility of having negligible external perturbations. For such orbits we expect a ~10^-3 test of the frame dragging and a ~10^-2 test of the no-hair theorem within 5 years, if Sgr A* is spinning rapidly. Our method is also capable of identifying perturbations caused by distributed mass around Sgr A*, thus providing high confidence in these gravity tests. Our analysis is not affected by uncertainties in our knowledge of the distance to the Galactic center, R0. A combination of pulsar timing with the astrometric results of stellar orbits would greatly improve the measurement precision of R0.

Obtaining pulsar parallaxes via relative astrometry (also known as differential astrometry) yields distances and transverse velocities that can be used to probe properties of the pulsar population and the interstellar medium. Large programs are essential to obtain the sample sizes necessary for these population studies, but they must be efficiently conducted to avoid requiring an infeasible amount of observing time. This paper describes the PSRPI astrometric program, including the use of new features in the DiFX software correlator to efficiently locate calibrator sources, selection and observing strategies for a sample of 60 pulsars, initial results, and likely science outcomes. Potential applications of high-precision relative astrometry to measure source structure evolution in defining sources of the International Celestial Referent Frame are also discussed.

We present 325 MHz (90 cm wavelength) radio observations of ultracool dwarfs TVLM 513-46546 and 2MASS J0036+1821104 using the Very Large Array (VLA) in June 2007. Ultracool dwarfs are expected to be undetectable at radio frequencies, yet observations at 8.5 GHz (3.5 cm) and 4.9 GHz (6 cm) of have revealed sources with > 100 {\mu}Jy quiescent radio flux and > 1 mJy pulses coincident with stellar rotation. The anomalous emission is likely a combination of gyrosynchrotron and cyclotron maser processes in a long-duration, large-scale magnetic field. Since the characteristic frequency for each process scales directly with the magnetic field magnitude, emission at lower frequencies may be detectable from regions with weaker field strength. We detect no significant radio emission at 325 MHz from TVLM 513-46546 or 2MASS J0036+1821104 over multiple stellar rotations, establishing 2.5{\sigma} total flux limits of 795 {\mu}Jy and 942 {\mu}Jy respectively. Analysis of an archival VLA 1.4 GHz observation of 2MASS J0036+1821104 from January 2005 also yields a non-detection at the level of < 130 {\mu}Jy . The combined radio observation history (0.3 GHz to 8.5 GHz) for these sources suggests a continuum emission spectrum for ultracool dwarfs which is either flat or inverted below 2-3 GHz. Further, if the cyclotron maser instability is responsible for the pulsed radio emission observed on some ultracool dwarfs, our low-frequency non-detections suggest that the active region responsible for the high-frequency bursts is confined within 2 stellar radii and driven by electron beams with energies less than 5 keV.

A concept for a new space-based cosmology mission called the Dark Ages Radio Explore (DARE) is presented in this paper. DARE's science objectives include (1) When did the first stars form? (2) When did the first accreting black holes form? (3) When did Reionization begin? (4) What surprises does the end of the Dark Ages hold (e.g., Dark Matter decay)? DARE will use the highly-redshifted hyperfine 21-cm transition from neutral hydrogen to track the formation of the first luminous objects by their impact on the intergalactic medium during the end of the Dark Ages and during Cosmic Dawn (redshifts z=11-35). It will measure the sky-averaged spin temperature of neutral hydrogen at the unexplored epoch 80-420 million years after the Big Bang, providing the first evidence of the earliest stars and galaxies to illuminate the cosmos and testing our models of galaxy formation. DARE's approach is to measure the expected spectral features in the sky-averaged, redshifted 21-cm signal over a radio bandpass of 40-120 MHz. DARE orbits the Moon for a mission lifetime of 3 years and takes data above the lunar farside, the only location in the inner solar system proven to be free of human-generated radio frequency interference and any significant ionosphere. The science instrument is composed of a three-element radiometer, including electrically-short, tapered, bi-conical dipole antennas, a receiver, and a digital spectrometer. The smooth frequency response of the antennas and the differential spectral calibration approach using a Markov Chain Monte Carlo technique will be applied to detect the weak cosmic 21-cm signal in the presence of the intense solar system and Galactic foreground emissions.

By volume, more than 99% of the solar system has not been imaged at radio frequencies. Almost all of this space (the solar wind) can be traversed by fast electrons producing radio emissions at frequencies lower than the terrestrial ionospheric cutoff, which prevents observation from the ground. To date, radio astronomy-capable space missions consist of one or a few satellites, typically far from each other, which measure total power from the radio sources, but cannot produce images with useful angular resolution. To produce such images, we require arrays of antennas distributed over many wavelengths (hundreds of meters to kilometers) to permit aperture synthesis imaging. Such arrays could be free-flying arrays of microsatellites or antennas laid out on the lunar surface. In this white paper, we present the lunar option. If such an array were in place by 2020, it would provide context for observations during Solar Probe Plus perihelion passes. Studies of the lunar ionosphere's density and time variability are also important goals. This white paper applies to the Solar and Heliospheric Physics study panel.

We report on a systematic and sensitive search for pulsars and transient sources in the nearby spiral galaxy M33, conducted at 1.4 GHz with the Arecibo telescope's seven-beam receiver system, ALFA. Data were searched for both periodic and aperiodic sources, up to 1000 pc cm^{-3} in dispersion measure and on timescales from $\sim$50 \us to several seconds. The galaxy was sampled with 12 ALFA pointings, or 84 pixels in total, each of which was searched for 2-3 hr. We describe the observations, search methodologies and analysis strategies applicable to multibeam systems, and comment on the data quality and statistics of spurious events that arise due to radio frequency interference. While these searches have not led to any conclusive signals of periodic or transient nature that originate in the galaxy, they illustrate some of the underlying challenges and difficulties in such searches and the efficacy of simultaneous multiple beams in the analysis of search output. The implied limits are $\la$5 uJy Mpc^2 in luminosity (at 1400 MHz) for periodic sources in M33 with duty cycles $\la$5%. For short-duration transient signals (with pulse widths $\la$100 \us), the limiting peak flux density is 100 mJy, which would correspond to a 5-$\sigma$ detection of bright giant pulses ($\sim$20 kJy) from Crab-like pulsars if located at the distance of M33. We discuss the implications of our null results for various source populations within the galaxy and comment on the future prospects to conduct even more sensitive searches with the upcoming next-generation instruments including the Square Kilometer Array and its pathfinders.

The Australian Square Kilometre Array Pathfinder (ASKAP) is a 36-element array with a 30-square-degree field of view being built at the proposed SKA site in Western Australia. We are conducting a Design Study for pulsar observations with ASKAP, planning both timing and search observations. We provide an overview of the ASKAP telescope and an update on pulsar-related progress.

We present X-ray, infrared, optical and radio observations of four previously unidentified Galactic plane X-ray sources, AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. Detection of each source with the Chandra X-ray Observatory has provided sub-arcsecond localizations, which we use to identify bright infrared counterparts to all four objects. Infrared and optical spectroscopy of these counterparts demonstrate that all four X-ray sources are extremely massive stars, with spectral classifications Ofpe/WN9 (AX J163252-4746), WN7 (AX J184738-0156 = WR121a), WN7-8h (AX J144701-5919) and OIf+ (AX J144547-5931). AX J163252-4746 and AX J184738-0156 are both luminous, hard, X-ray emitters with strong Fe XXV emission lines in their X-ray spectra at ~6.7 keV. The multi-wavelength properties of AX J163252-4746 and AX J184738-0156 are not consistent with isolated massive stars or accretion onto a compact companion; we conclude that their X-ray emission is most likely generated in a colliding-wind binary system. For both AX J144701-5919 and AX J144547-5931, the X-ray emission is an order of magnitude less luminous and with a softer spectrum. These properties are consistent with a colliding-wind binary interpretation for these two sources also, but other mechanisms for the generation of X-rays cannot be excluded. There are many other as yet unidentified X-ray sources in the Galactic plane, with X-ray properties similar to those seen for AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. This may indicate a substantial population of X-ray-emitting massive stars and colliding-wind binaries in the Milky Way.

This paper presents a search for radio transients at a frequency of 73.8 MHz (4 m wavelength) using the all-sky imaging capabilities of the Long Wavelength Demonstrator Array (LWDA). The LWDA was a 16-dipole phased array telescope, located on the site of the Very Large Array in New Mexico. The field of view of the individual dipoles was essentially the entire sky, and the number of dipoles was sufficiently small that a simple software correlator could be used to make all-sky images. From 2006 October to 2007 February, we conducted an all-sky transient search program, acquiring a total of 106 hr of data; the time sampling varied, being 5 minutes at the start of the program and improving to 2 minutes by the end of the program. We were able to detect solar flares, and in a special-purpose mode, radio reflections from ionized meteor trails during the 2006 Leonid meteor shower. We detected no transients originating outside of the solar system above a flux density limit of 500 Jy, equivalent to a limit of no more than about 10^{-2} events/yr/deg^2, having a pulse energy density >~ 1.5 x 10^{-20} J/m^2/Hz at 73.8 MHz for pulse widths of about 300 s. This event rate is comparable to that determined from previous all-sky transient searches, but at a lower frequency than most previous all-sky searches. We believe that the LWDA illustrates how an all-sky imaging mode could be a useful operational model for low-frequency instruments such as the Low Frequency Array, the Long Wavelength Array station, the low-frequency component of the Square Kilometre Array, and potentially the Lunar Radio Array.

This paper reports Very Large Array observations at 325 and 1425 MHz (90cm and 20cm) during and near the periastron passage of HD 80606b on 2007 November 20. We obtain flux density limits (3-sigma) of 1.7 mJy and 48 microJy at 325 and 1425 MHz, respectively, equivalent to planetary luminosity limits of 2.3 x 10^{24} erg/s and 2.7 x 10^{23} erg/s. These are well above the Jovian value (at 40 MHz) of 2 x 10^{18} erg/s. The motivation for these observations was that the planetary magnetospheric emission is driven by a stellar wind-planetary magnetosphere interaction so that the planetary luminosity would be elevated. Near periastron, HD 80606b might be as much as 3000 times more luminous than Jupiter. Recent transit observations of HD 80606b provide stringent constraints on the planetary mass and radius, and, because of the planet's highly eccentric orbit, its rotation period is likely to be "pseudo-synchronized" to its orbital period, allowing a robust estimate of the former. We are able to make robust estimates of the emission frequency of the planetary magnetospheric emission and find it to be around 60--90 MHz. We compare HD 80606b to other high-eccentricity systems and assess the detection possibilities for both near-term and more distant future systems. Of the known high eccentricity planets, only HD 80606b is likely to be detectable, as HD 20782B b and HD 4113b are both likely to have weaker magnetic field strengths. Both the forthcoming "EVLA low band" system and the Low Frequency Array may be able to improve upon our limits for HD 80606b, and do so at a more optimum frequency. If the low-frequency component of the Square Kilometre Array (SKA-lo) and a future lunar radio array are able to approach their thermal noise limits, they should be able to detect an HD 80606b-like planet, unless the planet's luminosity increases by substantially less than a factor of 3000.

This paper considers the suitability of a number of emerging and future instruments for the study of radio recombination lines (RRLs) at frequencies below 200 MHz. These lines arise only in low-density regions of the ionized interstellar medium, and they may represent a frequency-dependent foreground for next-generation experiments trying to detect H I signals from the Epoch of Reionization and Dark Ages ("21-cm cosmology"). We summarize existing decametre-wavelength observations of RRLs, which have detected only carbon RRLs. We then show that, for an interferometric array, the primary instrumental factor limiting detection and study of the RRLs is the areal filling factor of the array. We consider the Long Wavelength Array (LWA-1), the LOw Frequency ARray (LOFAR), the low-frequency component of the Square Kilometre Array (SKA-lo), and a future Lunar Radio Array (LRA), all of which will operate at decametre wavelengths. These arrays offer digital signal processing, which should produce more stable and better defined spectral bandpasses; larger frequency tuning ranges; and better angular resolution than that of the previous generation of instruments that have been used in the past for RRL observations. Detecting Galactic carbon RRLs, with optical depths at the level of 10^-3, appears feasible for all of these arrays, with integration times of no more than 100 hr. The SKA-lo and LRA, and the LWA-1 and LOFAR at the lowest frequencies, should have a high enough filling factor to detect lines with much lower optical depths, of order 10^-4 in a few hundred hours. The amount of RRL-hosting gas present in the Galaxy at the high Galactic latitudes likely to be targeted in 21-cm cosmology studies is currently unknown. If present, however, the spectral fluctuations from RRLs could be comparable to or exceed the anticipated H I signals.

We report detection of strong circularly polarized emission from the transient bursting source GCRT J1745-3009 based on new analysis of 325 MHz GMRT observations conducted on 28 September 2003. We place 8 Solar radius as the upper limit on the size of the emission region. The implied high brightness temperature required for an object beyond 1 pc and the high fraction of circular polarization firmly establish the emission as coherent. Electron cyclotron or plasma emission from a highly subsolar magnetically dominated dwarf located less than 4 kpc away could have given rise to the GCRT radio emission.

We are developing a purely commensal survey experiment for fast (<5s) transient radio sources. Short-timescale transients are associated with the most energetic and brightest single events in the Universe. Our objective is to cover the enormous volume of transients parameter space made available by ASKAP, with an unprecedented combination of sensitivity and field of view. Fast timescale transients open new vistas on the physics of high brightness temperature emission, extreme states of matter and the physics of strong gravitational fields. In addition, the detection of extragalactic objects affords us an entirely new and extremely sensitive probe on the huge reservoir of baryons present in the IGM. We outline here our approach to the considerable challenge involved in detecting fast transients, particularly the development of hardware fast enough to dedisperse and search the ASKAP data stream at or near real-time rates. Through CRAFT, ASKAP will provide the testbed of many of the key technologies and survey modes proposed for high time resolution science with the SKA.

This paper reports a blind search for magnetospheric emissions from planets around nearby stars. Young stars are likely to have much stronger stellar winds than the Sun, and because planetary magnetospheric emissions are powered by stellar winds, stronger stellar winds may enhance the radio luminosity of any orbiting planets. Using various stellar catalogs, we selected nearby stars (<~ 30 pc) with relatively young age estimates (< 3 Gyr). We constructed different samples from the stellar catalogs, finding between 100 and several hundred stars. We stacked images from the 74-MHz (4-m wavelength) VLA Low-frequency Sky Survey (VLSS), obtaining 3\sigma limits on planetary emission in the stacked images of between 10 and 33 mJy. These flux density limits correspond to average planetary luminosities less than 5--10 x 10^{23} erg/s. Using recent models for the scaling of stellar wind velocity, density, and magnetic field with stellar age, we estimate scaling factors for the strength of stellar winds, relative to the Sun, in our samples. The typical kinetic energy carried by the stellar winds in our samples is 15--50 times larger than that of the Sun, and the typical magnetic energy is 5--10 times larger. If we assume that every star is orbited by a Jupiter-like planet with a luminosity larger than that of the Jovian decametric radiation by the above factors, our limits on planetary luminosities from the stacking analysis are likely to be a factor of 10--100 above what would be required to detect the planets in a statistical sense. Similar statistical analyses with observations by future instruments, such as the Low Frequency Array (LOFAR) and the Long Wavelength Array (LWA), offer the promise of improvements by factors of 10--100.

We report the discovery of three pulsars whose large dispersion measures and angular proximity to \sgr indicate the existence of a Galactic center population of neutron stars. The relatively long periods (0.98 to 1.48 s) most likely reflect strong selection against short-period pulsars from radio-wave scattering at the observation frequency of 2 GHz used in our survey with the Green Bank Telescope. One object (PSR J1746-2850I) has a characteristic spindown age of only 13 kyr along with a high surface magnetic field $\sim 4\times 10^{13}$ G. It and a second object found in the same telescope pointing, PSR J1746-2850II (which has the highest known dispersion measure among pulsars), may have originated from recent star formation in the Arches or Quintuplet clusters given their angular locations. Along with a third object, PSR J1745-2910, and two similar high-dispersion, long-period pulsars reported by Johnston et al. (2006), the five objects found so far are 10 to 15 arc min from \sgr, consistent with there being a large pulsar population in the Galactic center, most of whose members are undetectable in relatively low-frequency surveys because of pulse broadening from the same scattering volume that angularly broadens \sgr and OH/IR masers.

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames.

This paper reports dual-epoch, Very Long Baseline Array observations of H I absorption toward 3C 147. One of these epochs (2005) represents new observations while one (1998) represents the reprocessing of previous observations to obtain higher signal-to-noise results. Significant H I opacity and column density variations, both spatially and temporally, are observed with typical variations at the level of \Delta\tau ~ 0.20 and in some cases as large as \Delta\tau ~ 0.70, corresponding to column density fluctuations of order 5 x 10^{19} cm^{-2} for an assumed 50 K spin temperature. The typical angular scale is 15 mas; while the distance to the absorbing gas is highly uncertain, the equivalent linear scale is likely to be about 10 AU. Approximately 10% of the face of the source is covered by these opacity variations, probably implying a volume filling factor for the small-scale absorbing gas of no more than about 1%. Comparing our results with earlier results toward 3C 138 (Brogan et al.), we find numerous similarities, and we conclude that small-scale absorbing gas is a ubiquitous phenomenon, albeit with a low probability of intercept on any given line of sight. Further, we compare the volumes sampled by the line of sight through the Galaxy between our two epochs and conclude that, on the basis of the motion of the Sun alone, these two volumes are likely to be substantially different. In order to place more significant constraints on the various models for the origin of these small-scale structures, more frequent sampling is required in any future observations.

We report the detection of a new transient radio source, GCRT J1742-3001, located ~1 degree from the Galactic center. The source was detected ten times from late 2006 to 2007 May in our 235 MHz transient monitoring program with the Giant Metrewave Radio Telescope (GMRT). The radio emission brightened in about one month, reaching a peak observed flux density of ~100 mJy on 2007 January 28, and decaying to ~50 mJy by 2007 May when our last monitoring observation was made. Two additional faint, isolated 235 MHz detections were made in mid-2006, also with the GMRT. GCRT J1742-3001 is unresolved at each epoch, with typical resolutions of ~20 arcsec x 10 arcsec. No polarization information is available from the observations. Based on nondetections in observations obtained simultaneously at 610 MHz, we deduce that the spectrum of GCRT J1742-3001 is very steep, with a spectral index less than about -2. Follow-up radio observations in 2007 September at 330 MHz and 1.4 GHz, and in 2008 February at 235 MHz yielded no detections. No X-ray counterpart is detected in a serendipitous observation obtained with the X-ray telescope aboard the Swift satellite during the peak of the radio emission in early 2007. We consider the possibilities that GCRT J1742-3001 is either a new member of an existing class of radio transients, or is representative of a new class having no associated X-ray emission.

We present a catalog of 9017 X-ray sources identified in Chandra observations of a 2 by 0.8 degree field around the Galactic center. We increase the number of known X-ray sources in the region by a factor of 2.5. The catalog incorporates all of the ACIS-I observations as of 2007 August, which total 2.25 Msec of exposure. At the distance to the Galactic center (8 kpc), we are sensitive to sources with luminosities >4e32 erg/s (0.5-8.0 keV; 90% confidence) over an area of one square degree, and up to an order of magnitude more sensitive in the deepest exposure (1.0 Msec) around Sgr A*. The positions of 60% of our sources are accurate to <1" (95% confidence), and 20% have positions accurate to <0.5". We search for variable sources, and find that 3% exhibit flux variations within an observation, 10% exhibit variations from observation-to-observation. We also find one source, CXOUGC J174622.7-285218, with a periodic 1745 s signal (1.4% chance probability), which is probably a magnetically-accreting cataclysmic variable. We compare the spatial distribution of X-ray sources to a model for the stellar distribution, and find 2.8 sigma evidence for excesses in the numbers of X-ray sources in the region of recent star formation encompassed by the Arches, Quintuplet, and Galactic center star clusters. These excess sources are also seen in the luminosity distribution of the X-ray sources, which is flatter near the Arches and Quintuplet than elsewhere in the field. These excess point sources, along with a similar longitudinal asymmetry in the distribution of diffuse iron emission that has been reported by other authors, probably have their origin in the young stars that are prominent at l~0.1 degree.

The radio sky is poorly sampled for rapidly varying transients because of the narrow field-of-view of most imaging radio telescopes at cm and shorter wavelengths. The emergence of sensitive long wavelength observations with intrinsically larger fields-of-view are changing this situation, as partly illustrated by our ongoing meter-wavelength monitoring observations and archival studies of the Galactic Center. In this search, we discovered a transient, bursting, radio source in the direction of the Galactic Center, GCRT J1745-3009, with extremely unusual properties. Its flux and rapid variability imply a brightness temperature >10^12 K if it is at a distance >70 pc, implying that it is a coherent emitter. I will discuss the discovery of the source and the subsequent re-detections, as well as searches for counterparts at other wavelengths, and several proposed models.

We present an optical/near-infrared search for a counterpart to the perplexing radio transient GCRT J1745-3009, a source located ~1 degree from the Galactic Center. Motivated by some similarities to radio bursts from nearby ultracool dwarfs, and by a distance upper limit of 70 pc for the emission to not violate the 1e12 K brightness temperature limit for incoherent radiation, we searched for a nearby star at the position of GCRT J1745-3009. We found only a single marginal candidate, limiting the presence of any late-type star to >1 kpc (spectral types earlier than M9), >200 pc (spectral types L and T0-T4), and >100 pc (spectral types T4-T7), thus severely restricting the possible local counterparts to GCRT J1745-3009. We also exclude any white dwarf within 1 kpc or a supergiant star out to the distance of the Galactic Center as possible counterparts. This implies that GCRT J1745-3009 likely requires a coherent emission process, although whether or not it reflects a new class of sources is unclear.

Using both the Very Large Array (VLA) at 7mm wavelength, and the Australia Telescope Compact Array (ATCA) at 3mm, we have searched for microwave emission from from cool dust in the extrasolar planetary system Gliese 876 (Gl 876). Having detected no emission above our 3-sigma detection threshold of 135 microJy, we rule out any dust disk with either a mass greater than 0.0006 Earth masses or less than ~250 AU across. This result improves on previous detection aperture thresholds an order of magnitude greater, and it has some implications for the dynamical modeling of the system. It also is consistent with the Greaves et al. hypothesis that relates the presence of a debris disk to close-in planets. Due to the dust-planetesimal relationship, our null result may also provide a constraint on the population or composition of the dust and small bodies around this nearby M dwarf.

As a service to the community, we have compiled radio frequency spectra from the literature for all sources within the VLA Low Frequency Sky Survey (VLSS) that are brighter than 15 Jy at 74 MHz. Over 160 references were used to maximize the amount of spectral data used in the compilation of the spectra, while also taking care to determine the corrections needed to put the flux densities from all reference on the same absolute flux density scale. With the new VLSS data, we are able to vastly improve upon previous efforts to compile spectra of bright radio sources to frequencies below 100 MHz because (1) the VLSS flux densities are more reliable than those from some previous low frequency surveys and (2) the VLSS covers a much larger area of the sky (declination >-30 deg.) than many other low frequency surveys (e.g., the 8C survey). In this paper, we discuss how the spectra were constructed and how parameters quantifying the shapes of the spectra were derived. Both the spectra and the shape parameters are made available here to assist in the calibration of observations made with current and future low frequency radio facilities.

All of the solar system gas giants produce electron cyclotron masers, driven by the solar wind impinging on their magnetospheres. Extrapolating to the planet orbiting tau Boo, various authors have predicted that it may be within the detection limits of the 4-meter wavelength (74 MHz) system on the Very Large Array. This paper reports three epochs of observations of tau Boo. In no epoch do we detect the planet; various means of determining the upper limit to the emission yield single-epoch limits ranging from 135 to 300 mJy. We develop a likelihood method for multi-epoch observations and use it to constrain various radiation properties of the planet. Assuming that the planet does radiate at our observation wavelength, its typical luminosity must be less than about 10^{16} W, unless its radiation is highly beamed into a solid angle Omega << 1 sr. While within the range of luminosities predicted by various authors for this planet, this value is lower than recent estimates which attempt to take into account the stellar wind of tau Boo using the known properties of the star itself. Electron cyclotron maser emission from solar systems planets is beamed, but with characteristic solid angles of approximately 1 sr illuminated. Future long-wavelength instruments (e.g., the Long Wavelength Array and the Low Frequency Array) must be able to make typical flux density measurements on short time scales (~ 15 min.) of approximately 25 mJy in order to improve these constraints significantly.

We analyze a sample of 58 multi-wavelength, Very Long Baseline Array observations of active galactic nuclei (AGN) to determine their scattering properties. Approximately 75% of the sample consists of AGN that exhibit centimeter-wavelength intraday variability (interstellar scintillation) while the other 25% do not show intraday variability. We find that interstellar scattering is measurable for most of these AGN, and the typical broadening diameter is 2 mas at 1 GHz. We find that the scintillating AGN are typically at lower Galactic latitudes than the non-scintillating AGN, consistent with the scenario that intraday variability is a propagation effect from the Galactic interstellar medium. The magnitude of the inferred interstellar broadening measured toward the scintillating AGN, when scaled to higher frequencies, is comparable to the diameters inferred from analyses of the light curves for the more well-known intraday variable sources. However, we find no difference in the amount of scattering measured toward the scintillating versus non-scintillating AGN. A consistent picture is one in which the scintillation results from localized regions ("clumps") distributed throughout the Galactic disk, but which individually make little contribution to the angular broadening. Of the 58 AGN observed, 37 (64%) have measured redshifts. At best, a marginal trend is found for scintillating (non-scintillating) AGN to have smaller (larger) angular diameters at higher redshifts. We also use our observations to try to constrain the possibility of intergalactic scattering. While broadly consistent with the scenario of a highly turbulent intergalactic medium, our observations do not place significant constraints on its properties.

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 Naval Research Laboratory and the National Radio Astronomy Observatory completed implementation of a low frequency capability on the VLA at 73.8 MHz in 1998. This frequency band offers unprecedented sensitivity (~25 mJy/beam) and resolution (~25 arcsec) for low-frequency observations. We review the hardware, the calibration and imaging strategies, comparing them to those at higher frequencies, including aspects of interference excision and wide-field imaging. Ionospheric phase fluctuations pose the major difficulty in calibrating the array. Over restricted fields of view or at times of extremely quiescent ionospheric ``weather'', an angle-invariant calibration strategy can be used. In this approach a single phase correction is devised for each antenna, typically via self-calibration. Over larger fields of view or at times of more normal ionospheric ``weather'' when the ionospheric isoplanatic patch size is smaller than the field of view, we adopt a field-based strategy in which the phase correction depends upon location within the field of view. This second calibration strategy was implemented by modeling the ionosphere above the array using Zernike polynomials. Images of 3C sources of moderate strength are provided as examples of routine, angle-invariant calibration and imaging. Flux density measurements indicate that the 74 MHz flux scale at the VLA is stable to a few percent, and tied to the Baars et al. value of Cygnus A at the 5 percent level. We also present an example of a wide-field image, devoid of bright objects and containing hundreds of weaker sources, constructed from the field-based calibration. We close with a summary of lessons the 74 MHz system offers as a model for new and developing low-frequency telescopes. (Abridged)

GCRT J1745-3009 is a transient bursting radio source located in the direction of the Galactic center. It was discovered in a 330 MHz VLA observation from 2002 September 30--October 1 and subsequently rediscovered in a 330 MHz GMRT observation from 2003 September 28 by Hyman et al. Here we report a new radio detection of the source in 330 MHz GMRT data taken on 2004 March 20. The observed properties of the single burst detected differ significantly from those measured previously, suggesting that GCRT J1745-3009 was detected in a new physical state. The 2004 flux density was ~0.05 Jy, ~10x weaker than the single 2003 burst and ~30x weaker than the five bursts detected in 2002. We derive a very steep spectral index, alpha = -13.5 +/- 3.0, across the bandpass, a new result previously not detectable due to limitations in the analysis of the 2002 and 2003 observations. Also, the burst was detected for only ~2 min., in contrast to the 10 min. duration observed in the earlier bursts. Due to sparse sampling, only the single burst was detected in 2004, as in the 2003 epoch, and we cannot rule out additional undetected bursts that may have occurred with the same ~77 min. periodicity observed in 2002 or with a different periodicity. Considering our total time on source throughout both our archival and active monitoring campaigns, we estimate the source exhibits detectable bursting activity ~7% of the time.

The recently started Arecibo L-band Feed Array (ALFA) pulsar survey aims to find ~1000 new pulsars. Due to its high time and frequency resolution the survey is especially sensitive to millisecond pulsars, which have the potential to test gravitational theories, detect gravitational waves and probe the neutron-star equation of state. Here we report the results of our preliminary analysis: in the first months we have discovered 21 new pulsars. One of these, PSR J1906+0746, is a young 144-ms pulsar in a highly relativistic 3.98-hr low-eccentricity orbit. The 2.61 +- 0.02 solar-mass system is expected to coalesce in ~300 Myr and contributes significantly to the computed cosmic inspiral rate of compact binary systems.

This paper presents observations of Cygnus A at 74 and 327 MHz at angular resolutions of approximately 10" and 3", respectively. These observations are among the highest angular resolutions obtained below 1000 MHz for this object. While the angular resolution at 74 MHz is not sufficient to separate clearly the hot spots from the lobes, guided by 151 and 327 MHz images, we have estimated the 74 MHz emission from the hot spots. We confirm that the emission from both the western and eastern hot spots flattens at low frequencies and that there is a spectral asymmetry between the two. For the eastern hot spot, a low-energy cutoff in the electron energy spectrum appears to explain the flattening, which implies a cutoff Lorentz factor \gamma_min ~ 300, though we cannot exclude the possibility that there might be a moderate level of free-free absorption. For the western hot spot, the current observations are not sufficient to distinguish between a free-free absorped power-law spectrum and a synchrotron self-absorbed spectrum.

At low radio frequencies (< 100 MHz), classical HII regions may become optically thick (optical depth > 1) and can be observed as discrete absorption regions against the Galactic nonthermal background emission created by Galactic cosmic ray electrons spiraling around magnetic fields. In this work we present 74 MHz observations in the region 26>l>-15, -5<b<5, and report the detection of 92 absorption features associated with known HII regions, and derive the brightness temperature of the Galactic cosmic ray electron synchrotron emission emanating from the column behind these regions. For the 42 HII regions with known distances, the average emissivity of the column behind the HII region is derived. 74 MHz emissivity values range between 0.3 and 1.0 Kelvin per parsec for a model assuming uniform distribution of emissivity. Methods for utilizing this type of data to model the 3-dimensional distribution of cosmic ray emissivity and the possibility of using this method to break the HII region kinematic distance degeneracy are discussed.

We report the discovery of up to 35 new supernova remnants (SNRs) from a 42 arcsec resolution 90cm multi-configuration Very Large Array survey of the Galactic plane covering 4.5 deg< l <22.0 deg and |b| < 1.25 deg. Archival 20cm, 11cm, and 8 micron data have also been used to identify the SNRs and constrain their properties. The 90cm image is sensitive to SNRs with diameters 2.5 arcmin to 50 arcmin and down to a surface brightness limit of about 10^{-21} W m^{-2} Hz^{-1} sr^{-1}. This survey has nearly tripled the number of SNRs known in this part of the Galaxy, and represents an overall 15% increase in the total number of Galactic SNRs. These results suggest that further deep low frequency surveys of the inner Galaxy will solve the discrepancy between the expected number of Galactic SNRs and the significantly smaller number of currently known SNRs.

We report the discovery of PSR J1906+0746, a young 144-ms pulsar in a highly relativistic 3.98-hr orbit with an eccentricity of 0.085 and expected gravitational wave coalescence time of 300 Myr. The new pulsar was found during precursor survey observations with the Arecibo 1.4-GHz feed array system and retrospectively detected in the Parkes Multibeam plane pulsar survey data. From radio follow-up observations with Arecibo, Jodrell Bank, Green Bank, and Parkes, we have measured the spin-down and binary parameters of the pulsar and its basic spectral and polarization properties. We also present evidence for pulse profile evolution, which is likely due to geodetic precession, a relativistic effect caused by the misalignment of the pulsar spin and total angular momentum vectors. Our measurements show that PSR J1906+0746 is a young object with a characteristic age of 112 kyr. From the measured rate of orbital periastron advance 7.57+/-0.03 deg/yr, we infer a total system mass of 2.61+/-0.02 Msun. While these parameters suggest that the PSR J1906+0746 binary system might be a younger version of the double pulsar system, intensive searches for radio pulses from the companion have so far been unsuccessful. It is therefore not known whether the companion is another neutron star or a massive white dwarf. Regardless of the nature of the companion, a simple calculation suggests that the Galactic birth rate of binaries similar to PSR J1906+0746 is 60/Myr. This implies that PSR J1906+0746 will make a significant contribution to the computed cosmic inspiral rate of compact binary systems.

From 1979--1996 the Green Bank Interferometer was used by the Naval Research Laboratory to monitor the flux density from 146 compact radio sources at frequencies near 2 and 8 GHz. We filter the ``light curves'' to separate intrinsic variations on times of a year or more from more rapid interstellar scintilation (ISS) on times of 5--50 d. Whereas the intrinsic variation at 2 GHz is similar to that at 8 GHz (though diminished in amplitude), the ISS variation is much stronger at 2 than at 8 GHz. We characterize the ISS variation by an rms amplitude and a timescale and examine the statistics of these parameters for the 121 sources with significant ISS at 2 GHz. We model the scintillations using the NE2001 Galactic electron model assuming the sources are brightness-limited. We find the observed rms amplitude to be in general agreement with the model, provided that the compact components of the sources have about 50% of their flux density in a component with maximum brightness temperatures $10^{11}$--$10^{12}$K. Thus our results are consistent with cm-wavelength VLBI studies of compact AGNs, in that the maximum brightness temperatures found are consistent with the inverse synchrotron limit at $3 \times 10^{11}$ K, boosted in jet configurations by Doppler factors up to about 20. The average of the observed 2 GHz ISS timescales is in reasonable agreement with the model at Galactic latitudes above about 10\de. At lower latitudes the observed timescales are too fast, suggesting that the transverse plasma velocity increases more than expected beyond about 1 kpc.

The highest velocity neutron stars establish stringent constraints on natal kicks, asymmetries in supernova core collapse, and the evolution of close binary systems. Here we present the first results of a long-term pulsar astrometry program using the VLBA. We measure a proper motion and parallax for the pulsar B1508+55, leading to model-independent estimates of its distance (2.37+0.23-0.20 kpc) and transverse velocity (1083+103-90 km/s), the highest velocity directly measured for a neutron star. We trace the pulsar back from its present Galactic latitude of 52.3 degrees to a birth site in the Galactic plane near the Cyg OB associations, and find that it will inevitably escape the Galaxy. Binary disruption alone is insufficient to impart the required birth velocity, and a natal kick is indicated. A composite scenario including a large kick along with binary disruption can plausibly account for the high velocity.

Determination of pulsar parallaxes and proper motions addresses fundamental astrophysical questions. We have recently finished a VLBI astrometry project to determine the proper motions and parallaxes of 27 pulsars, thereby doubling the total number of pulsar parallaxes. Here we summarise our astrometric technique and present the discovery of a pulsar moving in excess of 1000 km/s. As an example of the application of high precision pulsar astrometry we also infer the identification of 2 pulsars originating from a disrupted binary in the Cygnus Superbubble.

GCRT J1745-3009 is a transient bursting radio source located in the direction of the Galactic center, discovered in 330 MHz VLA observations from 2002 September 30--October 1 by Hyman et al. We have searched for bursting activity from GCRT J1745-3009 in nearly all of the available 330 MHz VLA observations of the Galactic center since 1989 as well as in 2003 GMRT observations. We report a new radio detection of the source in 330 MHz GMRT data taken on 2003 September 28. A single ~0.5 Jy burst was detected, approximately 3x weaker than the five bursts detected in 2002. Due to the sparse sampling of the 2003 observation, only the decay portion of a single burst was detected. We present additional evidence indicating that this burst is an isolated one, but we cannot completely rule out additional undetected bursts that may have occured with the same ~77 min. periodicity observed in 2002 or with a different periodicity. Assuming the peak emission was detected, the decay time of the burst, ~2 min, is consistent with that determined for the 2002 bursts. Based on the total time for which we have observations, we estimate that the source has a duty cycle of roughly 10%.

We present the discovery of non-thermal radio and X-ray emission positionally coincident with the TeV $\gamma$-ray source HESS J1813--178. We demonstrate that the non-thermal radio emission is due to a young shell-type supernova remnant (SNR) G12.8--0.0, and constrain its distance to be greater than 4 kpc. The X-ray emission is primarily non-thermal and is consistent with either an SNR shell or unidentified pulsar/pulsar wind nebula origin; pulsed emission is not detected in archival ASCA data. A simple synchrotron+inverse Compton model for the broadband emission assuming that all of the emission arises from the SNR shell implies maximum energies of $(30-450)(B/10 microG)^{-0.5}$ TeV. Further observations are needed to confirm that the broadband emission has a common origin and better constrain the X-ray spectrum.

New low-frequency 74 and 330 MHz observations of the Galactic center (GC) region reveal the presence of a large-scale ($6\arcdeg\times 2\arcdeg$) diffuse source of nonthermal synchrotron emission. A minimum energy analysis of this emission yields a total energy of $\sim (\phi^{4/7}f^{3/7})\times 10^{52}$ ergs and a magnetic field strength of $\sim 6(\phi/f)^{2/7}$ \muG (where $\phi$ is the proton to electron energy ratio and $f$ is the filling factor of the synchrotron emitting gas). The equipartition particle energy density is $1.2(\phi/f)^{2/7}$ \evcm, a value consistent with cosmic-ray data. However, the derived magnetic field is several orders of magnitude below the 1 mG field commonly invoked for the GC. With this field the source can be maintained with the SN rate inferred from the GC star formation. Furthermore, a strong magnetic field implies an abnormally low GC cosmic-ray energy density. We conclude that the mean magnetic field in the GC region must be weak, of order 10 \muG (at least on size scales $\ga 125\arcsec$).

We have previously presented the results of a 325 MHz radio survey of M31, conducted with the A-configuration of the Very Large Array. In this survey, a total of 405 radio sources between <6" and 170" in extent were mapped with a resolution of 6" and a sensitivity of ~0.6 mJy/beam. Here, we compare the resultant source list and image with other radio, IR, optical, X-ray observations and catalogs of the region. Through this, we were able to identify five supernova remnant (SNR) candidates and three pulsar wind nebula (PWN) candidates in M31, as well as three Milky Way radio stars, a possible Milky Way Planetary Nebula, and a bevy of interesting extragalactic objects: a BL Lac, a Giant Radio Galaxy, a galaxy merger, and several high-z radio galaxy candidates. In addition, a large number of compact (<6") extremely steep spectrum sources ($\alpha \leq -1.6; S_{\nu} \propto \nu^{\alpha}$) were detected whose nature is unknown.

The sheer number of pulsars discovered by the SKA, in combination with the exceptional timing precision it can provide, will revolutionize the field of pulsar astrophysics. The SKA will provide a complete census of pulsars in both the Galaxy and in Galactic globular clusters that can be used to provide a detailed map of the electron density and magnetic fields, the dynamics of the systems, and their evolutionary histories. This complete census will provide examples of nearly every possible outcome of the evolution of massive stars, including the discovery of very exotic systems such as pulsar black-hole systems and sub-millisecond pulsars, if they exist. These exotic systems will allow unique tests of the strong field limit of relativistic gravity and the equation of state at extreme densities. Masses of pulsars and their binary companions -- planets, white dwarfs, other neutron stars, and black holes -- will be determined to $\sim 1$% for hundreds of objects. With the SKA we can discover and time highly-stable millisecond pulsars that comprise a pulsar-timing array for the detection of low-frequency gravitational waves. The SKA will also provide partial censuses of nearby galaxies through periodicity and single-pulse detections, yielding important information on the intergalactic medium.

We present three epochs of VLBA observations of Galactic HI absorption toward the quasar 3C~138 with resolutions of 20 mas (~ 10 AU). This analysis includes VLBA data from observations in 1999 and 2002 along with a reexamination of 1995 VLBA data. Improved data reduction and imaging techniques have led to an order of magnitude improvement in sensitivity compared to previous work. With these new data we confirm the previously detected milliarcsecond scale spatial variations in the HI opacity at the level of Delta(tau_{max}) =0.50 \pm 0.05. The typical size scale of the optical depth variations is ~ 50 mas or 25 AU. In addition, for the first time we see clear evidence for temporal variations in the HI opacity over the seven year time span of our three epochs of data. We also attempted to detect the magnetic field strength in the HI gas using the Zeeman effect. From this analysis we have been able to place a 3 sigma upper limit on the magnetic field strength per pixel of ~45 muG. We have also been able to calculate for the first time the plane of sky covering fraction of the small scale HI gas of ~10%. We also find that the line widths of the milliarcsecond sizescale HI features are comparable to those determined from previous single dish measurements toward 3C~138, suggesting that the opacity variations cannot be due to changes in the HI spin temperature. From these results we favor a density enhancement interpretation for the small scale HI structures, although these enhancements appear to be of short duration and are unlikely to be in equilibrium.

Transient astronomical sources are typically powered by compact objects and usually signify highly explosive or dynamic events. While radio astronomy has an impressive record of obtaining high time resolution observations, usually it is achieved in quite narrow fields-of-view. Consequently, the dynamic radio sky is poorly sampled, in contrast to the situation in the X- and gamma-ray bands in which wide-field instruments routinely detect transient sources. Here we report a new transient source, GCRT J1745-3009, detected in 2002 during a moderately wide-field radio transient monitoring program of the Galactic center (GC) region at 0.33 GHz. The characteristics of its bursts are unlike those known for any other class of radio transient. If located in or near the GC, its brightness temperature (~10^16 K) and the implied energy density within GCRT J1745-3009 vastly exceeds that observed in most other classes of radio astronomical sources, and is consistent with coherent emission processes rarely observed. We conclude that GCRT J1745-3009 is the first member of a new class of radio transient sources, the first of possibly many new classes to be identified through current and upcoming radio surveys.

We present VLA images of the supernova remnant (SNR) 3C~391 at 74, 330, and 1465 MHz. This remnant has been known for some time to exhibit a turnover in its integrated radio continuum spectrum at frequencies < 100 MHz, indicative of free-free absorption from thermal ionized gas along the line of sight. For the first time, our data reveal the spatially resolved morphology of the low frequency free-free absorption with a resolution of ~70 arcsec. Contrary to the expectation that such absorption arises from unrelated low density HII regions (or their envelopes) along the line of sight, these data suggest that in this case the absorbing medium is directly linked to the SNR itself. 3C~391 has been shown in a number of recent papers to be interacting with a molecular cloud. Indeed, it exhibits a number of signposts of SNR/molecular cloud shocks including OH (1720 MHz) masers and broad molecular emission lines. Comparison of the regions of strongest 74 MHz absorption with existing X-ray, IR, and molecular data suggests that the free-free absorption originates from the SNR/molecular cloud shock boundaries due to ionized gas created from the passage of a J-type shock with a speed of ~100 km/s. This makes only the second SNR for which such (extrinsic) spatially resolved absorption has been measured, and the only one for which the absorption is thought to arise from a SNR/molecular cloud interface region.

RFI mitigation is a critically important issue in radio astronomy using existing instruments as well as in the development of next-generation radio telescopes, such as the Square Kilometer Array (SKA). Most designs for the SKA involve multiple stations with spacings of up to a few thousands of kilometers and thus can exploit the drastically different RFI environments at different stations. As demonstrator observations and analysis for SKA-like instruments, and to develop RFI mitigation schemes that will be useful in the near term, we recently conducted simultaneous observations with Arecibo Observatory and the Green Bank Telescope (GBT). The observations were aimed at diagnosing RFI and using the mostly uncorrelated RFI between the two sites to excise RFI from several generic kinds of measurements such as giant pulses from Crab-like pulsars and weak HI emission from galaxies in bands heavily contaminated by RFI. This paper presents observations, analysis, and RFI identification and excision procedures that are effective for both time series and spectroscopy applications using multi-station data.

The Square Kilometre Array is conceived as a telescope which will both test fundamental physical laws and transform our current picture of the Universe. However, the scientific challenges outlined in this book are today's problems--will they still be the outstanding problems that will confront astronomers in the period 2020 to 2050 and beyond, when the SKA will be in its most productive years? If history is any example, the excitement of the SKA will not be in the old questions which are answered, but the new questions that will be raised by the new types of observations it will permit. The SKA is a tool for as-yet-unborn users and there is an onus on its designers to allow for the exploration of the unknown. We outline a philosophy for the design and operation of the SKA that can lead the radio astronomers in the 21st century to add to the many discoveries of new phenomena made by radio astronomers in the 20th century.

Radio waves are imprinted with propagation effects from ionized media through which they pass. Owing to electron density fluctuations, compact sources (pulsars, masers, and compact extragalactic sources) can display a wide variety of scattering effects. These scattering effects, particularly interstellar scintillation, can be exploited to provide *superresolution*, with achievable angular resolutions (<~ 1 microarcsecond) far in excess of what can be obtained by very long baseline interferometry on terrestrial baselines. Scattering effects also provide a powerful sub-AU probe of the microphysics of the interstellar medium, potentially to spatial scales smaller than 100 km, as well as a tracer of the Galactic distribution of energy input into the interstellar medium through a variety of integrated measures. Coupled with future gamma-ray observations, SKA observations also may provide a means of detecting fainter compact gamma-ray sources. Though it is not yet clear that propagation effects due to the intergalactic medium are significant, the SKA will either detect or place stringent constraints on intergalactic scattering.

Transient radio sources are necessarily compact and usually are the locations of explosive or dynamic events, therefore offering unique opportunities for probing fundamental physics and astrophysics. In addition, short-duration transients are powerful probes of intervening media owing to dispersion, scattering and Faraday rotation that modify the signals. While radio astronomy has an impressive record obtaining high time resolution, usually it is achieved in quite narrow fields of view. Consequently, the dynamic radio sky is poorly sampled, in contrast to the situation in the X-ray and gamma-ray bands. The SKA has the potential to change this situation, opening up new parameter space in the search for radio transients. We summarize the wide variety of known and hypothesized radio transients and demonstrate that the SKA offers considerable power in exploring this parameter space. Requirements on the SKA to search the parameter space include the abilities to (1) Make targeted searches using beamforming capability; (2) Conduct blind, all-sky surveys with dense sampling of the frequency-time plane in wide fields; (3) Sample the sky with multiple fields of view from spatially well-separated sites in order to discriminate celestial and terrestrial signals; (4) Utilize as much of the SKA's aggregate collecting area as possible in blind surveys, thus requiring a centrally condensed configuration; and (5) Localize repeating transient sources to high angular precision, requiring a configuration with long baselines, thus requiring collecting area in both a centrally condensed "core" array and sufficient area on long baselines.

The sensitivity of the SKA enables a number of tests of theories of gravity. A Galactic Census of pulsars will discover most of the active pulsars in the Galaxy beamed toward us. In this census will almost certainly be pulsar-black hole binaries as well as pulsars orbiting the super-massive black hole in the Galactic centre. These systems are unique in their capability to probe the ultra-strong field limit of relativistic gravity. These measurements can be used to test the Cosmic Censorship Conjecture and the No-Hair theorem. The large number of millisecond pulsars discovered with the SKA will also provide a dense array of precision clocks on the sky. These clocks will act as the multiple arms of a huge gravitational wave detector, which can be used to detect and measure the stochastic cosmological gravitational wave background that is expected from a number of sources.

We present a wide field, sub-arcminute resolution VLA image of the Galactic Center region at 330 MHz. With a resolution of ~ 7" X 12" and an RMS noise of 1.6 mJy/beam, this image represents a significant increase in resolution and sensitivity over the previously published VLA image at this frequency. The improved sensitivity has more than tripled the census of small diameter sources in the region, has resulted in the detection of two new Non Thermal Filaments (NTFs), 18 NTF candidates, 30 pulsar candidates, reveals previously known extended sources in greater detail, and has resulted in the first detection of Sagittarius A* in this frequency range. A version of this paper containing full resolution images may be found at http://lwa.nrl.navy.mil/nord/AAAB.pdf.

We present here the results of a 325 MHz radio survey of M31, conducted with the A-configuration of the Very Large Array. The survey covered an area of 7.6 deg$^2$, and a total of 405 radio sources between $\la$6\arcsec and 170\arcsec in extent were mapped with a resolution of 6\arcsec and a 1$\sigma$ sensitivity of $\sim$0.6 \mjyb. For each source, its morphological class, major axis $\theta_M$, minor axis $\theta_m$, position angle $\theta_{PA}$, peak flux $I$, integrated flux density $S$, spectral index $\alpha$ and spectral curvature parameter $\phi$ were calculated. A comparison of the flux and radial distribution -- both in the plane of the sky and in the plane of M31 -- of these sources with those of the XMM--LSS and WENSS radio surveys revealed that a vast majority of sources detected are background radio galaxies. As a result of this analysis, we expect that only a few sources are intrinsic to M31. These sources are identified and discussed in an accompanying paper.

New high-resolution, wide-field 90 cm VLA observations of the Galactic Center region by Nord et al. have revealed twenty nonthermal filament (NTF) candidates. We report 6 cm polarization observations of six of these. All of the candidates have the expected NTF morphology, and two show extended polarization confirming their identification as NTFs. One of the new NTFs appears to be part of a system of NTFs located in the Sgr B region, 64 pc in projection north of Sgr A. These filaments cross the Galactic plane with an orientation similar to the filaments in the Galactic Center Radio Arc. They extend the scale over which the NTF phenomena is known to occur to almost 300 pc along the Galactic plane. Another NTF was found in the Galactic plane south of the Sgr C filament but with an orientation of 45 degrees to the Galactic plane. This is only the second of 12 confirmed NTFs that is not oriented perpendicular to the Galactic plane. An additional candidate in the Sgr C region was resolved into multiple filamentary structures. Polarization was detected only at the brightness peak of one of the filaments. Several of these filaments run parallel to the Galactic plane and can be considered additional evidence for non-poloidal magnetic fields at the GC. Together the 90 and 6 cm observations indicate that the Galactic center magnetic field may be more complex than a simple globally ordered dipole field.

We report the detection of Sagittarius A*, the radio source associated with our Galaxy's central massive black hole, at 330 MHz with the Very Large Array. Implications for the spectrum and emission processes of Sagittarius A* are discussed and several hypothetical geometries of the central region are considered.

We present the first pulsar parallaxes measured with phase-referenced pulsar VLBI observations at 5 GHz. Due to the steep spectra of pulsars, previous astrometric measurements have been at lower frequencies. However, the strongest pulsars can be observed at 5 GHz, offering the benefit of lower combined ionospheric and tropospheric phase errors, which usually limit VLBI astrometric accuracy. The pulsars B0329+54, B0355+54 and B1929+10 were observed for 7 epochs spread evenly over 2 years. For B0329+54, large systematic errors lead to only an upper limit on the parallax (pi < 1.5 mas). A new proper motion and parallax were measured for B0355+54 (pi = 0.91 +- 0.16 mas), implying a distance of 1.04+0.21-0.16 kpc and a transverse velocity of 61+12-9 km/s. The parallax and proper motion for B1929+10 were significantly improved (pi = 2.77 +- 0.07 mas), yielding a distance of 361+10-8 pc and a transverse velocity of 177+4-5 km/s. We demonstrate that the astrometric errors are correlated with the angular separation between the phase reference calibrator and the target source, with significantly lower errors at 5 GHz compared to 1.6 GHz. Finally, based on our new distance determinations for B1929+10 and B0355+54, we derive or constrain the luminosities of each pulsar at high energies. We show that, for thermal emission models, the emitting area for X-rays from PSR B1929+10 is roughly consistent with the canonical size for a heated polar cap, and that the conversion of spin-down power to gamma-ray luminosity in B0355+54 must be low. The new proper motion for B1929+10 also implies that its progenitor is unlikely to have been the binary companion of the runaway O-star zeta-Ophiuchi.

We present a 74 MHz survey of a 165 square degree region located near the north galactic pole. This survey has an unprecedented combination of both resolution (25'' FWHM) and sensitivity (RMS noise as low as 24 mJy/beam). We detect 949 sources at the 5-sigma level in this region, enough to begin exploring the nature of the 74 MHz source population. We present differential source counts, spectral index measurements and the size distribution as determined from counterparts in the high resolution FIRST 1.4 GHz survey. We find a trend of steeper spectral indices for the brighter sources. Further, there is a clear correlation between spectral index and median source size, with the flat spectrum sources being much smaller on average. Ultra-steep spectrum objects (alpha < -1.2) are identified, and we present high resolution VLA follow-up observations of these sources which, identified at such a low frequency, are excellent candidates for high redshift radio galaxies.

We have imaged a 1 deg^2 field centered on the known Galactic supernova remnant (SNR) G11.2-0.3 at 74, 330, and 1465 MHz with the Very Large Array radio telescope (VLA) and 235 MHz with the Giant Metrewave Radio Telescope (GMRT). The 235, 330, and 1465 MHz data have a resolution of 25 arcsec, while the 74 MHz data have a resolution of 100 arcsec. The addition of this low frequency data has allowed us to confirm the previously reported low frequency turnover in the radio continuum spectra of the two known SNRs in the field: G11.2-0.3 and G11.4-0.1 with unprecedented precision. Such low frequency turnovers are believed to arise from free-free absorption in ionized thermal gas along the lines of site to the SNRs. Our data suggest that the 74 MHz optical depths of the absorbing gas is 0.56 and 1.1 for G11.2-0.3 and G11.4-0.1, respectively. In addition to adding much needed low frequency integrated flux measurements for two known SNRs, we have also detected three new SNRs: G11.15-0.71, G11.03-0.05, and G11.18+0.11. These new SNRs have integrated spectral indices between -0.44 and -0.80. Because of confusion with thermal sources, the high resolution (compared to previous Galactic radio frequency surveys) and surface brightness sensitivity of our observations have been essential to the identification of these new SNRs. With this study we have more than doubled the number of SNRs within just a 1 deg^2 field of view in the inner Galactic plane. This result suggests that future low frequency observations of the Galactic plane of similar quality may go a long way toward alleviating the long recognized incompleteness of Galactic SNR catalogs.

High resolution VLA observations of the Arches cluster near the Galactic center show evidence of continuum emission at $\lambda$3.6, 6, 20 and 90cm. The continuum emission at $\lambda$90cm is particularly striking because thermal sources generally become optically thick at longer wavelengths and fall off in brightness whereas non-thermal sources increase in brightness. It is argued that the radio emission from this unique source has compact and diffuse components produced by thermal and nonthermal processes, respectively. Compact sources within the cluster arise from stellar winds of mass-losing stars (Lang, Goss & Rodriguez 2001a) whereas diffuse emission is likely to be due to colliding wind shocks of the cluster flow generating relativistic particles due to diffuse shock acceleration. We also discuss the possibility that $\gamma$-ray emission from 3EG J1746--2851, located within 3.3$'$ of the Arches cluster, results from the inverse Compton scattering of the radiation field of the cluster.

In Paper I we present a new model for the Galactic distribution of free electrons. In this paper we describe the input data and methodology for determining the structure and parameters of the model. Tables of the input data are provided and several figures are used to demonstrate why particular Galactic structures are needed. We identify lines of sight on which discrete regions either enhance or diminish the dispersion or the scattering. Most do not coincide with known \ion{H}{2} regions or supershells, most likely because the enhancements correspond to column densities smaller than detection thresholds for the emission measure in recombination-line surveys.

We present a new model for the Galactic distribution of free electrons. It (a) describes the distribution of the free electrons responsible for pulsar dispersion measures and thus can be used for estimating the distances to pulsars; (b) describes large-scale variations in the strength of fluctuations in electron density that underly interstellar scattering; (c) can be used to interpret interstellar scattering and scintillation observations of Galactic objects and of extragalactic objects, such as intrinsically compact AGNs and Gamma-ray burst afterglows; and (d) serves as a preliminary, smooth spatial model of the warm ionized component of the interstellar gas. This work builds upon and supercedes the Taylor & Cordes (1993) model by exploiting new observations and analysis methods. For lines of sight directed out of the Galactic plane, the new model yields substantially larger values for pulsar dispersion measures, except for directions dominated by the local hot bubble. Unlike the TC93 model, the new model provides sufficient electrons to account for the dispersion measures of the vast majority of known, Galactic pulsars. The new model is described and exemplified using plots of astronomically useful quantities on Galactic-coordinate grids. Software available on the Internet is also described. Future observations and analysis techniques that will improve the Galactic model are outlined.

This paper reports new, more sensitive observations of nine of the eleven extrastatistical signals in the Megachannel Extraterrestrial Assay (META). These extrastatistical signals had all of the expected characteristics of a transmission from an extraterrestrial transmitter, except that they did not repeat. Cordes, Lazio, & Sagan showed that this lack of repeatability could be explained by the high detection thresholds used in the reobservations of these candidates combined with interstellar scintillation of intrinsically steady sources. We use the Cordes, Lazio, & Sagan methodology, correcting an error in the original presentation, and our new observations to rule out this scintillation hypothesis at a confidence level of at least 97.8% (for the case of an intrinsically weak source) to a level in excess of 99% (if the source strengths are comparable to that favored by Cordes, Lazio, & Sagan. We also demonstrate that gravitational microlensing cannot account for the initial detection of these candidate signals nor is gravitational lensing likely to play a role in future SETI programs. We conclude that the META candidates do not reflect a large population of powerful, strong beacons.

We report the detection of a new radio transient source, GCRT J1746-2757, located only 1.1 degrees north of the Galactic center. Consistent with other radio transients toward the Galactic center, this source brightened and faded on a time scale of a few months. No X-ray counterpart was detected. We also report new 0.33 GHz measurements of the radio counterpart to the X-ray transient source, XTE J1748-288, previously detected and monitored at higher radio frequencies. We show that the spectrum of XTE J1748-288 steepened considerably during a period of a few months after its peak. We also discuss the need for a more efficient means of finding additional radio transients.

We report 60 and 90 micron observations of 7 millisecond pulsars with the ISOPHOT instrument and describe our analysis procedures. No pulsars were detected, and typical (3\sigma) upper limits are 150 mJy. We combine our results with others in the literature and use them to place constraints on the existence of protoplanetary or dust disks around millisecond pulsars.

We report the discovery of a system of isolated nonthermal filaments approximately 0.5 deg. northwest (75 pc in projection) of Sgr A. Unlike other isolated nonthermal filaments which show subfilamentation, braiding of subfilaments, and flaring at their ends, these filaments are simple linear structures and more closely resemble the parallel bundled filaments in the Galactic center radio arc. However, the most unusual feature of these filaments is that the 20/90 cm spectral index uniformly decreases as a function of length, in contrast to all other nonthermal filaments in the Galactic center. This spectral gradient may not be due to simple particle aging but could be explained by a curved electron energy spectrum embedded in a diverging magnetic field. If so, the scale of the magnetic gradient is not consistent with a large scale magnetic field centered on Sgr A* suggesting that this filament system is tracing a local magnetic field.

Supernova remnants expanding into adjacent molecular clouds are believed to be sites of cosmic ray acceleration and sources of energetic gamma-rays. Under certain environmental conditions, such interactions also give rise to unusual OH masers in which the 1720 MHz satellite line dominates over the more common 1665/7 MHz emission. Motivated by the apparent coincidence of a handful of EGRET sources with OH(1720 MHz) maser-producing supernova remnants, we have carried out a search using the Very Large Array for new OH(1720 MHz) masers within the error regions of 11 unidentified EGRET sources at low Galactic latitude. While a previously known maser associated with an HII region was serendipitously detected, initial results indicate that no new masers were found down to a limiting flux of, typically, 50 mJy. We discuss the implications of this result on the nature of the unidentified Galactic EGRET sources.

We present multi-frequency observations with the Very Long Baseline Array of the compact double radio source B2 2050+36. Our observations are at 0.33, 0.61, 1.67, 2.3, and 8.4 GHz, with the 0.61 GHz observations forming the third epoch of observation of this source at that frequency. At 0.61 GHz, the structure of B2 2050+36 is dominated by two components 56 mas apart. Within the uncertainties of the various measurements, this separation has remained unchanged for the past 16 years. Any differential image wander caused by refractive interstellar scattering is less than 4 mas. Both the lack of differential image wander and the frequency dependence of the angular diameter of B2 2050+36 below 1 GHz indicate that the electron density power spectrum along this line of sight has a spectral index near the Kolmogorov value, with a value of 4 being highly unlikely. We conclude that diffractive scattering dominates along this line of sight; results in the literature indicate that this conclusion also holds true for the line of sight to the pulsar PSR B2020+28 (8.7 deg. from B2 2050+36). Comparison of our 1.67 GHz observations with those obtained 21 years previously place a limit on the projected linear separation velocity of the two components of c.

We present sub-arcminute resolution imaging of the Galactic supernova remnant W49B at 74 MHz (25") and 327 MHz (6"), the former being the lowest frequency at which the source has been resolved. While the 327 MHz image shows a shell-like morphology similar to that seen at higher frequencies, the 74 MHz image is considerably different, with the southwest region of the remnant almost completely attenuated. The implied 74 MHz optical depth (~ 1.6) is much higher than the intrinsic absorption levels seen inside two other relatively young remnants, Cas A and the Crab Nebula, nor are natural variations in the relativistic electron energy spectra expected at such levels. The geometry of the absorption is also inconsistent with intrinsic absorption. We attribute the absorption to extrinsic free-free absorption by a intervening cloud of thermal electrons. Its presence has already been inferred from the low-frequency turnover in the integrated continuum spectrum and from the detection of radio recombination lines toward the remnant. Our observations confirm the basic conclusions of those measurements, and our observations have resolved the absorber into a complex of classical HII regions surrounded either partially or fully by low-density HII gas. We identify this low-density gas as an extended HII region envelope (EHE), whose statistical properties were inferred from low resolution meter- and centimeter-wavelength recombination line observations. Comparison of our radio images with HI and H_2CO observations show that the intervening thermal gas is likely associated with neutral and molecular material as well.

We present light curves for 149 sources monitored with the Green Bank Interferometer. The light curves are at two radio frequencies (approximately 2.5 and 8.2 GHz) and range from 3 to 15 yrs in length, covering the interval 1979--1996, and have a typical sampling of one flux density measurement every 2 days. We have used these light curves to conduct various variability analysis (rms flux density variations and autoregressive, integrated, moving average modeling) of these sources. We find suggestive, though not unambiguous evidence, that these sources have a common, broadband mechanism for intrinsic variations, in agreement with previous studies of a subset of these source. We also find that the sources generally display a short-term variability (~ 10 d) that arises from radio-wave scattering in an extended medium. These conclusions extend those of Fiedler et al. (1987) who used a sub-sample of these data. The primary motivation for this monitoring program was the identification of extreme scattering events. In an effort to identify ESEs in a systematic manner, we have taken the wavelet transform of the light curves. We find 15 events in the light curves of 12 sources that we classify as probable ESEs. However, we also find that five ESEs previously identified from these data do not survive our wavelet selection criteria. Future identification of ESEs will probably continue to rely on both visual and systematic methods. Instructions for obtaining the data are also presented.

This paper considers scattering screens that have arbitrary spatial variations of scattering strength transverse to the line of sight, including screens that are spatially well confined, such as disks and filaments. We calculate the scattered image of a point source and the observed pulse shape of a scattered impulse. The consequences of screen confinement include: (1) Source image shapes that are determined by the physical extent of the screen rather than by the shapes of much-smaller diffracting microirregularities. These include image elongations and orientations that are frequency dependent. (2) Variation with frequency of angular broadening that is much weaker than the trademark \nu^{-2} scaling law (for a cold, unmagnetized plasma), including frequency-independent cases; and (3) Similar departure of the pulse broadening time from the usually expected \nu^{-4} scaling law. We briefly discuss applications that include scattering of pulses from the Crab pulsar by filaments in the Crab Nebula; image asymmetries from Galactic scattering of the sources Cyg X-3, Sgr A*, and NGC 6334B; and scattering of background active galactic nuclei by intervening galaxies. We also address the consequences for inferences about the shape of the wavenumber spectrum of electron density irregularities, which depend on scaling laws for the image size and the pulse broadening. Future low-frequency (< 100 MHz) array observations will also be strongly affected by the Galactic structure of scattering material. Our formalism is derived in the context of radio scattering by plasma density fluctuations. It is also applicable to optical, UV and X-ray scattering by grains in the interstellar medium.

Results are presented from a long-term astrometry program on PSR B0919+06 using the NRAO Very Long Baseline Array. With ten observations (seven epochs) between 1994--2000, we measure a proper motion of 18.35 +/- 0.06 mas/yr in RA, 86.56 +/- 0.12 mas/yr in Dec, and a parallax of 0.83 +/- 0.13 mas (68% confidence intervals). This yields a pulsar distance of 1.21 +/- 0.19 kpc, making PSR B0919+06 the farthest pulsar for which a trigonometric parallax has been obtained, and the implied pulsar transverse speed is 505 +/- 80 km/s. Combining the distance estimate with interstellar scintillation data spanning 20 years, we infer the existence of a patchy or clumpy scattering screen along the line of sight in addition to the distributed electron density predicted by models for the Galaxy, and constrain the location of this scattering region to within about 250 parsecs of the Sun. Comparison with the lines of sight towards other pulsars in the same quadrant of the Galaxy permits refinement of our knowledge of the local interstellar matter in this direction.

We report multi-epoch VLA H I absorption observations of the source 1741-038 (OT-068) before and during an extreme scattering event (ESE). Observations at four epochs, three during the ESE, were obtained. We find no changes in the equivalent width, maximum optical depth, or velocity of maximum optical depth during the ESE, but we do find a secular trend of decreasing maximum optical depth between our observations and ones by other observers a decade prior. The resulting limit on the H I column density change during the ESE for a structure with a spin temperature T_s is 6.4 x 10^{17} cm^{-2} (T_s/10 K). Tiny-scale atomic structures (TSAS), with a column density N_H ~ 3 x 10^{18} cm^{-2}, are ruled out marginally by this limit, though geometric arguments may allow this limit to be relaxed. Galactic halo molecular clouds, that are opaque in the H I line, cannot be excluded because the observed velocity range covers only 25% of their allowed velocity range.

We review observational constraints on the structures responsible for extreme scattering events, focussing on a series of observations of the quasar PKS 1741-038. VLA observations were conducted to search for changes in the rotation measure and H I absorption during the ESE, while VLBI observations sought ESE-induced changes in the source's image. No RM changes were found implying B_|| < 12 mG, and no H I opacity changes were found implying N(H I) < 6.4 x 10^{17} cm^{-2}. No multiple imaging was observed, but the diameter of the source increased by 0.7 mas, contrary to what is predicted by simple refractive lens modeling of ESEs. We summarize what these limits imply about the structure responsible for this ESE.

(Abridged) We report multi-epoch VLBI observations of the source PKS 1741-038 as it underwent an extreme scattering event. Observations at four epochs were obtained, and images were produced at three of these. During the event the source consisted of a dominant, compact component, essentially identical to the structure seen outside the event. However, the source's diameter increased slightly at 13 cm during the ESE. An increase in the source's diameter is inconsistent with a simple refractive model. We also see no evidence for ESE-induced substructure within the source or the formation of multiple images, as would occur in a strongly refractive lens. However, a model in which the decrease in flux density during the ESE occurs solely because of stochastic broadening within the lens requires a larger broadening diameter during the ESE than is observed. Thus, the ESE toward 1741-038 involved both stochastic broadening and refractive defocussing within the lens. If the structure responsible for the ESE has a size of order 1 AU, the level of scattering within an ESE lens may be a factor of 10^7 larger than that in the ambient medium. A filamentary structure could reduce the difference between the strength of scattering in the lens and ambient medium, but we conclude that, if ESEs arise from filamentary structures, they occur when the filamentary structures are seen lengthwise. We predict the amount of pulse broadening that would result from a comparable lens passing in front of a pulsar. The pulse broadening would be no more than 1.1 microseconds, consistent with the lack of pulse broadening detected during ESEs toward the pulsars PSR B1937+21 and PSR J1643-1224.

We report a search for radio continuum emission from the gamma-ray pulsar Geminga. We have used the VLA to image the location of the optical counterpart of Geminga at 74 and 326 MHz. We detect no radio counterpart. We derive upper limits to the pulse-averaged flux density of Geminga, taking diffractive scintillation into account. We find that diffractive scintillation is probably quenched at 74 MHz and does not influence our upper limit, S < 56 mJy (2\sigma), but that a 95% confidence level at 326 MHz is S < 5 mJy. Owing to uncertainties on the other low-frequency detections and the possibility of intrinsic variability or extrinsic variability (refractive interstellar scintillation) or both, our non-detections are nominally consistent with these previous detections.

We report the discovery of a new non-thermal filament, G358.85+0.47, the ``Pelican'', located ~225 pc in projection from SgrA, and oriented parallel to the Galactic plane. VLA continuum observations at 20 cm reveal that this 7' (17.5 pc) structure bends at its northern extension and is comprised of parallel strands, most apparent at its ends. Observations at 6 and 3.6 cm reveal that the Pelican is a synchrotron-emitting source and is strongly linearly polarized over much of its extent. The spectral index of the filament changes from alpha(20/6)=-0.8 to alpha(6/3.6)=-1.5. The rotation measures exhibit a smooth gradient, with values ranging from -1000 rad/m2 to +500 rad/m2. The intrinsic magnetic field is well-aligned along the length of the filament. Based on these properties, we classify the Pelican as one of the non-thermal filaments unique to the Galactic center. Since these filaments (most of which are oriented perpendicular to the Galactic plane) are believed to trace the overall magnetic field in the inner Galaxy, the Pelican is the first detection of a component of this field parallel to the plane. The Pelican may thus mark a transition region of the magnetic field orientation in the inner kiloparsec of the Galaxy.

We report the discovery of a new filamentary structure, G358.85+0.47, consisting of at least three mutually parallel but bent `strands', located about 1.5 degrees SW of the Sgr A complex. Unlike all the other known Galactic center filaments, G358.85+0.47 is oriented parallel to the galactic plane. This orientation of the filament may have implications for the large scale structure of the magnetic field near the Galactic Center. This feature was identified in a new wide-field image of the Galactic Center region at 90 cm. Further higher resolution observations at 20 cm revealed a filamentary structure which is similar to several other known linear features in the vicinity of the Galactic center. Based on its appearance in a 20 cm image, we give it the name 'the Pelican' for further reference.

G359.87+0.18 is an enigmatic object located 15' from Sgr A*. It has been variously classified as an extragalactic source, Galactic jet source, and young supernova remnant. We present new observations of G359.87+0.18 between 0.33 and 15 GHz and use these to argue that this source is an Faranoff-Riley II radio galaxy. We are able to place a crude limit on its redshift of z > 0.1. The source has a spectral index \alpha < -1 (S \propto \nu^\alpha), suggestive of a radio galaxy with a redshift z >~ 2. The scattering diameters of Sgr A* and several nearby OH masers (~ 1" at 1 GHz) indicate that a region of enhanced scattering is along the line of sight to the Galactic center. If the region covers the Galactic center uniformly, the implied diameter for a background source is at least 600" at 0.33 GHz, in contrast with the observed 20" diameter of G359.87+0.18. Using the scattering diameter of a nearby OH maser OH 359.762+0.120 and the widths of two, nearby, non-thermal threads, G0.08+0.15 and G359.79+0.17, we show that a uniform scattering region should cover G359.87+0.18. We therefore conclude that the Galactic center scattering region is inhomogeneous on a scale of 5' (~ 10 pc at a distance of 8.5 kpc). This scale is comparable to the size scale of molecular clouds in the Galactic center. The close agreement between these two lengths scales is an indication that the scattering region is linked intimately to the Galactic center molecular clouds.

The scattering diameters of Sgr A* and several nearby OH masers (~ 1" at 1 GHz) indicate that a region of enhanced scattering is along the line of sight to the Galactic center. The scattering diameter of an extragalactic source seen through this scattering region is larger by the ratio (Sun-GC distance)/(GC-scattering region separation). This ratio could be a factor of a few, if the scattering region is far from the GC and only a random superposition with it, to more than 100, if the scattering region is within the GC. We have used the VLA to survey 10 (11) fields at 20 cm (6 cm) ranging from 7' to 137' from Sgr A*. Our objective was to identify extragalactic sources and measure their scattering diameters so as to constrain this ratio. In order to find sources within these fields, we have employed pdfCLEAN, an algorithm in which sources are identified in an image by comparing the intensity histogram of the image to that expected from a noise-only image. We found over 100 sources, with the faintest sources being 3 mJy. The average number of sources per field is approximately 10, though fields close to Sgr A* tend to contain fewer sources. In a companion paper we use this survey in assessing the likelihood that the scattering region is within the GC. A number of Galactic sources are included in our source catalog. We discuss a double-lobed source, potentially an X-ray quiet version of 1E1740.7-2942, a shell-like structure with a central point source, and a possible radio transient.

The scattering diameters of Sgr A* and several nearby OH masers (~ 1" at 1 GHz) indicate that a region of enhanced scattering is along the line of sight to the Galactic center. We combine radio-wave scattering data and free-free emission and absorption measurements in a likelihood analysis that constrains the following parameters of the GC scattering region: The GC-scattering region separation, d; the angular extent of the region, \psi_l; the outer scale on which density fluctuations occur, l_0; and the gas temperature, T. The maximum likelihood estimates of these parameters are d = 133_{-80}^{+200} pc, 0.5 degrees <= \psi_l <~ 1 degrees, and (l_0/1 pc)^{2/3}T^{-1/2} = 10^{-7 +/- 0.8}. As host media for the scattering, we consider the photoionized surface layers of molecular clouds and the interfaces between molecular clouds and the 10^7 K ambient gas. We are unable to make an unambiguous determination, but we favor an interface model in which the scattering medium is hot (T ~ 10^6 K) and dense (n_e ~ 10 cm^{-3}). The GC scattering region produces a 1 GHz scattering diameter for an extragalactic source of 90", if the region is a single screen, or 180", if the region wraps around the GC, as appears probable. We modify the Taylor-Cordes model for the Galactic distribution of free electrons in order to include an explicit GC component. Pulsars seen through this region will have a dispersion measure of approximately 2000 pc cm^{-3}, of which 75% arises from the GC component. We stress the uniqueness of the GC scattering region, probably resulting from the high-pressure environment in the GC.

We use radio-wave scattering data to constrain the distribution of ionized gas in the outer Galaxy. Like previous models, our model for the H II disk includes parameters for the radial scale length and scale height of the H II, but we allow the H II disk to warp and flare. Our model also includes the Perseus arm. We use a likelihood analysis on 11 extragalactic sources and 7 pulsars. Scattering in the Perseus arm is no more than 60% of the level contributed by spiral arms in the inner Galaxy, equivalent to a 1 GHz scattering diameter of 1.5 mas. Our analysis favors an unwarped, nonflaring disk with a 1 kpc scale height, though this may reflect the non-uniform and coarse coverage provided by the available data. The lack of a warp indicates that VLBI observations near 1 GHz with an orbiting station having baseline lengths of a few Earth diameters will not be affected by interstellar scattering at Galactic latitudes |b| ~ 15 degrees. The radial scale length is 15--20 kpc, but the data cannot distinguish between a gradual decrease in the electron density and a truncated distribution. We favor a truncated one, because we associate the scattering with massive star formation, which is also truncated near 20 kpc. The distribution of electron density turbulence decreases more rapidly with Galactocentric distance than does the hydrogen distribution. Alternate ionizing and turbulent agents---the intergalactic ionizing flux and satellite galaxies passing through the disk---do not contribute significantly to scattering. We cannot exclude the possibility that a largely ionized, but quiescent disk extends to >~ 100 kpc, similar to that for some Ly-alpha absorbers.

We report multifrequency Very Long Baseline Array observations of twelve active galactic nuclei seen toward the Galactic anticenter. All of the sources are at |b| < 10 degrees and seven have |b| < 0.5 degrees. Our VLBA observations can detect an enhancement in the angular broadening of these sources due to an extended H II disk, if the orientation of the H II disk in the outer Galaxy is similar to that of the H I disk. Such an extended H II disk is suggested by the C IV absorption in a quasar's spectrum, the appearance of H I disks of nearby spiral galaxies, and models of Ly-alpha cloud absorbers and the Galactic fountain. We detect eleven of the twelve sources at one or more frequencies; nine of the sources are compact and suitable for an angular broadening analysis. A preliminary analysis of the observed angular diameters suggests that the H II disk does not display considerable warping or flaring and does not extend to large Galactocentric distances (R >~ 100 kpc). A companion paper (Lazio & Cordes 1997) combines these observations with those in the literature and presents a more comprehensive analysis.

We consider the geometrical optics for refraction of a distant radio source by an interstellar plasma lens, with application to a lens with a Gaussian electron column density profile. The refractive properties of the lens are specified completely by a dimensionless parameter, alpha, which is a function of the wavelength of observation, the lens' electron column density, the lens-observer distance, and the transverse diameter of the lens. Relative motion of the observer and lens produces modulations in the source's light curve. Plasma lenses are diverging so the light curve displays a minimum, when the lens is on-axis, surrounded by enhancements above the unlensed flux density. Lensing can also produce caustics, multiple imaging, and angular position wander of the background source. If caustics are formed, the separation of the outer caustics can constrain alpha, while the separation of the inner caustics can constrain the size of the lens. We apply our analysis to 0954+654, a source for which we can identify caustics in its light curve, and 1741-038, for which polarization observations were obtained during and after the scattering event. We find general agreement between modelled and observed light curves at 2.25 GHz, but poor agreement at 8.1 GHz. The discrepancies may result from a combination of lens substructure or anisotropic shape, a lens that only grazes the source, or unresolved source substructure. Our analysis places the following constraints on the lenses: Toward 0954+654 (1741-038) the lens was 0.38 AU (0.065 AU) in diameter, with a peak column density of 0.24 pc cm^{-3} (1E-4 pc cm^{-3}) and an electron density of 1E5 cm^{-3} (300 cm^{-3}). The angular wander caused by the lens was 250 mas (0.4 mas) at 2.25 GHz. For 1741-038, we place an upper limit of 100 mG on the lens' magnetic field.

We consider interstellar scintillations as a cause of intermittency in radio signals from extraterrestrial intelligence (ETI). We demonstrate that scintillations are very likely to allow initial detections of narrowband signals from distant sources (> 100 pc), while making redetections improbable. We consider three models in order to assess the non-repeating, narrowband events found in recent SETI and to analyze large surveys in general: (I) Radiometer noise; (II) A population of constant Galactic sources undergoing interstellar scintillation,; and (III) Real, transient signals (or hardware errors) of either terrestrial or ET origin. We apply likelihood and Bayesian tests of the models to The Planetary Society/Harvard META data. We find that Models II and III are both highly preferred to Model I, but that Models II and III are about equally likely. Ruling out Model II in favor of Model III requires many more reobservations than were conducted in META *or* the reobservation threshold must be much lower than was used in META. *We cannot rule out the possibility that META events are real, intrinsically steady ETI signals.* We recommend that future surveys use thresholds far below the typical false-alarm threshold to lessen the effects of intermittency. The threshold level is best defined in terms of the recording and computational technology that is available at a cost commensurate with other survey costs.

Radio-wave scattering is enhanced dramatically for Galactic center sources in a region with radius >~ 15 arc min. Using scattering from Sgr A* and other sources, we show that pulse broadening for pulsars in the Galactic center is {\em at least} 6.3 \nu^{-4} seconds (\nu = radio frequency in GHz) and is most likely 50--200 times larger because the relevant scattering screen appears to be within the Galactic center region itself. Pulsars beyond---but viewed through---the Galactic center suffer even greater pulse broadening and are angularly broadened by <~ 2 {\em arc min}. Periodicity searches at radio frequencies are likely to find only long period pulsars and, then, only if optimized by using frequencies >~ 7 GHz and by testing for small numbers of harmonics in the power spectrum. The optimal frequency is $\nu ~ 7.3 GHz (\Delta_{0.1}P\sqrt{\alpha})^{-1/4}$ where \Delta_{0.1} is the distance of the scattering region from Sgr A* in units of 0.1 kpc, P is the period (seconds), and \alpha is the spectral index. A search for compact sources using aperture synthesis may be far more successful than searches for periodicities because the angular broadening is not so large as to desensitize the survey. We estimate that the number of {\em detectable} pulsars in the Galactic center may range from <= 1 to 100, with the larger values resulting from recent, vigorous starbursts. Such pulsars provide unique opportunities for probing the ionized gas, gravitational potential, and stellar population near Sgr A*.

We have observed the hard X-ray source 1E1740.7-2942 in the HCO+ (1-0) line using the Owens Valley millimeter interferometer. Previous single dish observations have found HCO+ emission coincident with the location of the radio continuum hot spots of the radio source associated with 1E1740.7-2942. Our higher resolution observations show a 15 arcsec offset between the HCO+ emission and the location of the radio hot spots. We propose that the lack of emission results from a large ionization rate, exceeding $10^{-15}$ s$^{-1}$, in the neighborhood of 1E1740.7-2942.