The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentially explained by a progenitor system that is asymmetric in terms of the stellar component masses, i.e. with a mass ratio $q$ of 0.7-0.8, the known Galactic population of merging double neutron star (DNS) systems (i.e. those that will coalesce within billions of years or less) has, until now, only consisted of nearly equal-mass ($q > 0.9$) binaries. PSR J1913+1102 is a DNS system in a 5-hour, low-eccentricity ($e = 0.09$) orbit, implying an orbital separation of 1.8 solar radii, with the two neutron stars predicted to coalesce in 470 million years due to gravitational-wave emission. Here we report that the masses of the two neutron stars, as measured by a dedicated pulsar timing campaign, are $1.62 \pm 0.03$ and $1.27 \pm 0.03$ solar masses for the pulsar and companion neutron star, respectively; with a measured mass ratio $q = 0.78 \pm 0.03$, it is the most asymmetric DNS among known merging systems. Based on this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30% (90% confidence) of the total population of merging DNS binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

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.

Pulsar dispersion measures (DMs) have been used to model the electron density of the interstellar medium (ISM) in the Galactic disk as a plane-parallel medium, despite significant scatter in the DM-distance distribution and strong evidence for inhomogeneities in the ISM. We use a sample of pulsars with independent distance measurements to evaluate a model of the local ISM in the thick disk of the Galaxy that incorporates turbulent fluctuations, clumps, and voids in the electron density. The latter two components are required because $\sim 1/3$ of the lines of sight are discrepant from a strictly plane parallel model. A likelihood analysis for smooth components of the model yields a scale height $z_0=1.57^{+0.15}_{-0.14}$ kpc and a mid-plane density $n_0=0.015 \pm 0.001$ cm$^{-3}$. The scatter in the DM-distance distribution is dominated by clumps and voids but receives significant contributions from a broad spectrum of density fluctuations, such as a Kolmogorov spectrum. The model is used to identify lines of sight with outlier values of DM. Three of these pulsars, J1614$-$2230, J1623$-$0908, and J1643$-$1224, lie behind known HII regions, and the electron density model is combined with H$\alpha$ intensity data to constrain the filling factors and other substructure properties of the HII regions (Sh 2-7 and Sh 2-27). Several pulsars also exhibit enhanced DM fluctuations that are likely caused by their lines of sight intersecting the superbubble GSH 238+00+09.

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.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

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.

Precision pulsar timing can be used for a variety of astrophysical tests from the detection of gravitational waves to probing the properties of the interstellar medium (ISM). Here we present analyses of the noise contributions to pulsar timing residuals from continuous multi-hour observations of seven millisecond pulsars (MSPs). We present scintillation bandwidth measurements for all MSPs in the sample, some for the first time, and scintillation timescale measurements and lower limits for all MSPs for the first time. In addition, we present upper limits on the contribution of pulse phase jitter to the timing residual error for all MSPs. These long observations also allow us to constrain variations in dispersion measures (DMs) on hour-long timescales for several millisecond pulsars. We find that there are no apparent DM variations in any of the MSPs studied on these timescales as expected. In light of new radio telescopes such as the Canadian Hydrogen Intensity Mapping Experiment (CHIME), which will be able to time many pulsars for a short time each day, we search for differences in timing precisions from continuous TOAs and from equivalent length time-discontinuous TOAs. We find no differences in the precision for any MSP in our sample, as expected. We conclude that the TOA variations are consistent with the expected breakdown into template-fitting, jitter, and scintillation errors.

Gravitational scattering events, in which the path of an interstellar object is deflected by a pulsar or the solar system, give rise to reflex motion which can potentially be detected using pulsar timing. We determine the form of the timing signal expected from a gravitational scattering event, which is ramp-like and resembles the signal produced by a glitch or a gravitational wave burst with memory (BWM), and investigate the prospects for detecting such a signal using a pulsar timing array. The level of timing precision currently achieved for some millisecond pulsars makes it possible to detect objects as small as $10^{-10}$ solar masses, less than the mass of the dwarf planet Ceres, at impact parameters as large as 1 AU. The signals produced by gravitational scattering could provide independent constraints on models of dark matter involving asteroid-mass objects or subhalos, and should be considered as potential false positives in searches for BWMs.

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.

Some short gamma-ray bursts (SGRBs) are thought to be caused by the mergers of binary neutron stars which may sometimes produce massive neutron star remnants capable of producing extragalactic fast radio bursts (FRBs). We conducted a deep search for FRBs from the sites of six low-redshift SGRBs. We collected high time- and frequency-resolution data from each of the sites for 10 hours using the 2 GHz receiver of the Green Bank Telescope. Two of the SGRB sites we targeted were visible with the Arecibo Radio Telescope with which we conducted an additional 10 hours of 1.4 GHz observations for each. We searched our data for FRBs using the GPU-optimized dedispersion algorithm $\texttt{heimdall}$ and the machine-learning-based package $\texttt{FETCH}$ (Fast Extragalactic Transient Candidate Hunter). We did not discover any FRBs, but would have detected any with peak flux densities in excess of 87 mJy at the Green Bank Telescope or 21 mJy at Arecibo with a signal-to-noise ratio of at least 10. The isotropic-equivalent energy of any FRBs emitted from these sites in our bands during our observations must not have exceeded a few times $10^{38}$ erg, comparable to some of the lowest energy bursts yet seen from the first known repeating FRB 121102.

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.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

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.

We report on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo PALFA survey. Phase-coherent timing solutions derived from 2.5 to 5 years of observations carried out at Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65-ms pulsar in a 39.6-day orbit for which we detect $\gamma$-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mass-function system with an orbital period of 3.3 hours that belongs to the non-eclipsing black widow population. We also present PSR J1932+1756, the longest-orbital-period (41.5 days) intermediate-mass binary pulsar known to date. In light of the numerous discoveries of binary MSPs over the past years, we characterize the Galactic distribution of known MSP binaries in terms of binary class. Our results support and strengthen previous claims that the scatter in the Galactic scale height distribution correlates inversely with the binary mass function. We provide evidence of observational biases against detecting the most recycled pulsars near the Galactic plane, which overestimates the scale height of lighter systems. A possible bimodality in the mass function of MSPs with massive white dwarfs is also reported.

We summarize our understanding of millisecond radio bursts from an extragalactic population of sources. FRBs occur at an extraordinary rate, thousands per day over the entire sky with radiation energy densities at the source about ten billion times larger than those from Galactic pulsars. We survey FRB phenomenology, source models and host galaxies, coherent radiation models, and the role of plasma propagation effects in burst detection. The FRB field is guaranteed to be exciting: new telescopes will expand the sample from the current ~80 unique burst sources (and a few secure localizations and redshifts) to thousands, with burst localizations that enable host-galaxy redshifts emerging directly from interferometric surveys. * FRBs are now established as an extragalactic phenomenon. * Only a few sources are known to repeat. Despite the failure to redetect other FRBs, they are not inconsistent with all being repeaters. * FRB sources may be new, exotic kinds of objects or known types in extreme circumstances. Many inventive models exist, ranging from alien spacecraft to cosmic strings but those concerning compact objects and supermassive black holes have gained the most attention. A rapidly rotating magnetar is a promising explanation for FRB 121102 along with the persistent source associated with it, but alternative source models are not ruled out for it or other FRBs. * FRBs are powerful tracers of circumsource environments, `missing baryons' in the IGM, and dark matter. * The relative contributions of host galaxies and the IGM to propagation effects have yet to be disentangled, so dispersion measure distances have large uncertainties.

We describe the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) and its efforts to directly detect and study gravitational waves and other synergistic physics and astrophysics using radio timing observations of millisecond pulsars.

We present the DSA-2000: a world-leading radio survey telescope and multi-messenger discovery engine for the next decade. The array will be the first true radio camera, outputting science-ready image data over the 0.7 - 2 GHz frequency range with a spatial resolution of 3.5 arcsec. With 2000 x 5 m dishes, the DSA-2000 will have an equivalent point-source sensitivity to SKA1-mid, but with ten times the survey speed. The DSA-2000 is envisaged as an all-sky survey instrument complementary to the ngVLA, and as a counterpart to the LSST (optical), SPHEREx (near-infrared) and SRG/eROSITA (X-ray) all-sky surveys. Over a five-year prime phase, the DSA-2000 will image the entire sky above declination -30 degrees every four months, detecting > 1 unique billion radio sources in a combined full-Stokes sky map with 500 nJy/beam rms noise. This all-sky survey will be complemented by intermediate and deep surveys, as well as spectral and polarization image cubes. The array will be a cornerstone for multi-messenger science, serving as the principal instrument for the US pulsar timing array community, and by searching for radio afterglows of compact object mergers detected by LIGO and Virgo. The array will simultaneously detect and localize ~10,000 fast radio bursts each year, realizing their ultimate use as a cosmological tool. The DSA-2000 will be proposed to the NSF Mid-Scale Research Infrastructure-2 program with a view to first light in 2026

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

In this paper, we present the results of timing observations of PSRs J1949+3106 and J1950+2414, two binary millisecond pulsars discovered in data from the Arecibo ALFA pulsar survey (PALFA). The timing parameters include precise measurements of the proper motions of both pulsars, which show that PSR J1949+3106 has a transversal motion very similar to that of an object in the local standard of rest. The timing also includes measurements of the Shapiro delay and the rate of advance of periastron for both systems. Assuming general relativity, these allow estimates of the masses of the components of the two systems; for PSR J1949+3106, the pulsar mass is $M_p \, = \, 1.34^{+0.17}_{-0.15} \, M_{\odot}$ and the companion mass $M_c \, = \, 0.81^{+0.06}_{-0.05}\, M_{\odot}$; for PSR J1950+2414 $M_p \, = \, 1.496 \, \pm \, 0.023\, M_{\odot}$ and $M_c \, = \, 0.280^{+0.005}_{-0.004}\, M_{\odot}$ (all values 68.3 % confidence limits). We use these masses and proper motions to investigate the evolutionary history of both systems: PSR J1949+3106 is likely the product of a low-kick supernova; PSR J1950+2414 is a member of a new class of eccentric millisecond pulsar binaries with an unknown formation mechanism. We discuss the proposed hypotheses for the formations of these systems in light of our new mass measurements.

We present 41 bursts from the first repeating fast radio burst discovered (FRB 121102). A deep search has allowed us to probe unprecedentedly low burst energies during two consecutive observations (separated by one day) using the Arecibo telescope at 1.4 GHz. The bursts are generally detected in less than a third of the 580-MHz observing bandwidth, demonstrating that narrow-band FRB signals may be more common than previously thought. We show that the bursts are likely faint versions of previously reported multi-component bursts. There is a striking lack of bursts detected below 1.35 GHz and simultaneous VLA observations at 3 GHz did not detect any of the 41 bursts, but did detect one that was not seen with Arecibo, suggesting preferred radio emission frequencies that vary with epoch. A power law approximation of the cumulative distribution of burst energies yields an index $-1.8\pm0.3$ that is much steeper than the previously reported value of $\sim-0.7$. The discrepancy may be evidence for a more complex energy distribution. We place constraints on the possibility that the associated persistent radio source is generated by the emission of many faint bursts ($\sim700$ ms$^{-1}$). We do not see a connection between burst fluence and wait time. The distribution of wait times follows a log-normal distribution centered around $\sim200$ s; however, some bursts have wait times below 1 s and as short as 26 ms, which is consistent with previous reports of a bimodal distribution. We caution against exclusively integrating over the full observing band during FRB searches, because this can lower signal-to-noise.

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.

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.

Remarkable progress has been made in understanding turbulent astrophysical plasmas in past decades including, notably, the solar wind and the interstellar medium. In the case of the solar wind, much of this progress has relied on in situ measurements from space-borne instruments. However, ground-based radio observations also have played a significant role and have the potential to play an even bigger role. In particular, using distant background sources (quasars, pulsars, satellite beacons) to transilluminate the foreground corona and solar wind, a variety of radio propagation phenomena can be used to map plasma properties of the solar corona and heliosphere, as well as the warm interstellar medium. These include angular broadening, interplanetary and interstellar scintillations, and differential Faraday rotation. These observations are highly complementary to in situ observations of the solar wind, and could be a mainstay of investigations into turbulence of the ISM. We point out that the Next Generation Very Large Array (ngVLA) fulfills all the requirements necessary to exploit radio observations of astrophysical turbulence fully.

Precision pulsar timing at the level of tens to hundreds of nanoseconds allows detection of nanohertz gravitational waves (GWs) from supermassive binary black holes (SMBBHs) at the cores of merging galaxies and, potentially, from exotic sources such as cosmic strings. The same timing data used for GW astronomy also yield precision masses of neutron stars orbiting other compact objects, constraints on the equation of state of nuclear matter, and precision tests of General Relativity, the Strong Equivalence Principle, and alternative theories of gravity. Timing can also lead to stringent constraints on the photon mass and on changes in fundamental constants and could reveal low mass objects (rogue planets, dark matter clumps) that traverse pulsar lines of sight. Data sets also allow modeling of the density, magnetic field, and turbulence in the interstellar plasma. Roughly 100 millisecond pulsars (MSPs) are currently being timed at $\sim$GHz frequencies using the largest radio telescopes in the world. These data sets currently represent ~1000 pulsar-years and will increase dramatically over the next decade. These topics are presented as a program of key science with flowdown technical requirements for achieving the science.

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).

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.

We present the detection of an unresolved radio source coincident with the position of the Type I superluminous supernova (SLSN) PTF10hgi ($z=0.098$) about 7.5 years post-explosion, with a luminosity of $L_\nu(6\,\,{\rm GHz})\approx 1.1\times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$. This represents the first detection of radio emission coincident with a SLSN on any timescale. We investigate various scenarios for the origin of the radio emission: star formation activity, an active galactic nucleus, an off-axis jet, and a non-relativistic supernova blastwave. While any of these would be quite novel if confirmed, none appear likely when taken in context of the other properties of the host galaxy, previous radio observations of SLSNe, the sample of long gamma-ray bursts (LGRBs), and the general population of hydrogen-poor SNe. Instead, the radio emission is reminiscent of the quiescent radio source associated with the repeating FRB 121102, which has been argued to be powered by a magnetar born in a SLSN or LGRB explosion several decades ago. We show that such a central engine powered nebula is consistent with the age and luminosity of the radio source. Our directed search for FRBs from the location of PTF10hgi using 40 min of VLA phased-array data reveals no detections to a limit of $22$ mJy ($7\sigma$; 10 ms duration). We outline several follow-up observations that can conclusively establish the origin of the radio emission.

Birefringence in ionized, magnetized media is usually measured as Faraday rotation of linearly polarized radiation. However, pulses propagating through regions with very large Faraday rotation measures (RMs) can split into circularly polarized components with measurable differences in arrival times $\propto \nu^{-3}\, \rm RM$, where $\nu$ is the radio frequency. Differential refraction from gradients in DM (dispersion measure) and RM can contribute a splitting time $\propto \vert \mathbf{\nabla_{\perp}} \rm DM \vert \vert \mathbf{\nabla_{\perp}}\rm RM \vert\,\nu^{-5}$. Regardless of whether the emitted pulse is unpolarized or linearly polarized, net circular polarization will be measured when splitting is a significant fraction of the pulse width. However, the initial polarization may be inferable from the noise statistics of the bursts. Extreme multipath scattering that broadens pulses can mask splitting effects. We discuss particular cases such as the Galactic center magnetar, J1749$-$2900, and the repeating fast radio burst source, FRB 121102. Both lines of sight have $\rm RM \sim 10^5~\rm rad~m^{-2}$ that yields millisecond splittings at frequencies well below $\sim 1 \text{ GHz}$. We also consider the splitting of nanosecond shot pulses in giant pulses from the Crab pulsar and the minimal effects of birefringence on precision pulsar timing. Finally, we explore the utility of two-dimensional coherent dedispersion with DM and RM as parameters.

FRB 121102 is the only known repeating fast radio burst source. Here we analyze a wide-frequency-range (1-8 GHz) sample of high-signal-to-noise, coherently dedispersed bursts detected using the Arecibo and Green Bank telescopes. These bursts reveal complex time-frequency structures that include sub-bursts with finite bandwidths. The frequency-dependent burst structure complicates the determination of a dispersion measure (DM); we argue that it is appropriate to use a DM metric that maximizes frequency-averaged pulse structure, as opposed to peak signal-to-noise, and find DM = 560.57 +/- 0.07 pc/cc at MJD 57644. After correcting for dispersive delay, we find that the sub-bursts have characteristic frequencies that typically drift lower at later times in the total burst envelope. In the 1.1-1.7 GHz band, the ~ 0.5-1-ms sub-bursts have typical bandwidths ranging from 100-400 MHz, and a characteristic drift rate of ~ 200 MHz/ms towards lower frequencies. At higher radio frequencies, the sub-burst bandwidths and drift rate are larger, on average. While these features could be intrinsic to the burst emission mechanism, they could also be imparted by propagation effects in the medium local to the source. Comparison of the burst DMs with previous values in the literature suggests an increase of Delta(DM) ~ 1-3 pc/cc in 4 years, though this could be a stochastic variation as opposed to a secular trend. This implies changes in the local medium or an additional source of frequency-dependent delay. Overall, the results are consistent with previously proposed scenarios in which FRB 121102 is embedded in a dense nebula.

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.

Plasma lensing events can have significant observational consequences, including flux density modulations and perturbations in pulse arrival times. In this paper we develop and apply a formalism that extends geometrical optics to describe the effects of two dimensional plasma lenses of arbitrary shape. We apply insights from catastrophe theory and the study of uniform asymptotic expansions of integrals to describe the lensing amplification close to fold caustics and in shadow regions, and explore the effects of image appearance and disappearance at caustics in the time of arrival (TOA) perturbations due to lensing. The enhanced geometric optics approach successfully reproduces the predictions from wave optics and can be efficiently used to simulate multifrequency TOA residuals during lensing events. Lensing will introduce perturbations both in the way the residuals change as a function of frequency and also in the magnitude and sign of the residuals averaged over a frequency band. The deviations from the expected dispersive $\nu^{-2}$ scaling will be most significant when including observations at low frequencies. We examine the consequences of lensing in the context of precision pulsar timing and touch on its potential relevance to the study of FRBs.

Radio SETI experiments aim to test the hypothesis that extraterrestrial civilizations emit detectable signals from communication, propulsion, or other technologies. The unprecedented capabilities of next generation radio telescopes, including ngVLA, will allow us to probe hitherto unexplored regions of parameter space, thereby placing meaningful limits on the prevalence of technological civilizations in the Universe (or, if we are fortunate, making one of the most significant discoveries in the history of science). ngVLA provides critical capabilities in the 10 - 100 GHz range, and will be a valuable complement to SKA in the southern hemisphere, as well as surveying the sky at frequencies underexplored by previous SETI experiments.

Pulsars in the Galactic Center (GC) are important probes of General Relativity, star formation, stellar dynamics, stellar evolution, and the interstellar medium. Despite years of searching, only a handful of pulsars in the central 0.5 deg are known. The high-frequency sensitivity of ngVLA will open a new window for discovery and characterization of pulsars in the GC. A pulsar in orbit around the GC black hole, Sgr A*, will provide an unprecedented probe of black hole physics and General Relativity.

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.

We report the detection of 72 new pulses from the repeating fast radio burst FRB 121102 in Breakthrough Listen C-band (4-8 GHz) observations at the Green Bank Telescope. The new pulses were found with a convolutional neural network in data taken on August 26, 2017, where 21 bursts have been previously detected. Our technique combines neural network detection with dedispersion verification. For the current application we demonstrate its advantage over a traditional brute-force dedis- persion algorithm in terms of higher sensitivity, lower false positive rates, and faster computational speed. Together with the 21 previously reported pulses, this observa- tion marks the highest number of FRB 121102 pulses from a single observation, total- ing 93 pulses in five hours, including 45 pulses within the first 30 minutes. The number of data points reveal trends in pulse fluence, pulse detection rate, and pulse frequency structure. We introduce a new periodicity search technique, based on the Rayleigh test, to analyze the time of arrivals, with which we exclude with 99% confidence pe- riodicity in time of arrivals with periods larger than 5.1 times the model-dependent time-stamp uncertainty. In particular, we rule out constant periods >10 ms in the barycentric arrival times, though intrinsic periodicity in the time of emission remains plausible.

The second data release from the Gaia mission (Gaia DR2) includes, among its billion entries, as- trometric parameters for binary companions to a number of known pulsars, including white dwarf companions to millisecond pulsars and the non-degenerate components of so-called "black widow" and "redback" systems. We find 22 such counterparts in DR2, of which 12 have statistically sig- nificant measurements of parallax. These DR2 optical proper motions and parallaxes provide new measurements of the distances and transverse velocities of the associated pulsars. For the most part, the results agree with existing radio interferometric and pulsar timing-based astrometry, as well as other distance estimates based on photometry or associations, and for some pulsars they provide the best known distance and velocity estimates. In particular, two of these pulsars have no previous dis- tance measurement: PSR J1227$-$4853, for which Gaia measures a parallax of $0.62\pm 0.16$ mas, and PSR J1431$-$4715, with a Gaia parallax of $0.64\pm 0.16$ mas. Using the Gaia distance measurements, we find that dispersion measure-based distance estimates calculated using the NE2001 and YMW16 Galactic electron density models are on average slightly underestimated, which may be a selection effect due to the over-representation of pulsars at high Galactic latitudes in the present Gaia sample. While the Gaia DR2 results do not quite match the precision that can be achieved by dedicated pulsar timing or radio interferometry, taken together they constitute a small but important improvement to the pulsar distance scale, and the subset of millisecond pulsars with distances measured by Gaia may help improve the sensitivity of pulsar timing arrays to nanohertz gravitational waves.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the solar wind from the full collection of DM measurements in the recently released NANOGrav 11-yr data set. We filter out long time scale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the solar wind as spherically symmetric and constant in time, we find the electron number density at 1 A.U. to be $7.9\pm0.2$ cm$^{-3}$. Our data shows little evidence of long-term variation in the density of the wind. We argue that our techniques paired with a high cadence, low radio frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind.

We present a newly implemented single-pulse pipeline for the PALFA survey to efficiently identify single radio pulses from pulsars, Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs). We have conducted a sensitivity analysis of this new pipeline in which multiple single pulses with a wide range of parameters were injected into PALFA data sets and run through the pipeline. Based on the recovered pulses, we find that for pulse widths $\rm < 5\ ms$ the sensitivity of the PALFA pipeline is at most a factor of $\rm \sim 2$ less sensitive to single pulses than our theoretical predictions. For pulse widths $\rm > 10\ ms$, as the $\rm DM$ decreases, the degradation in sensitivity gets worse and can increase up to a factor of $\rm \sim 4.5$. Using this pipeline, we have thus far discovered 7 pulsars and 2 RRATs and identified 3 candidate RRATs and 1 candidate FRB. The confirmed pulsars and RRATs have DMs ranging from 133 to 386 pc cm$^{-3}$ and flux densities ranging from 20 to 160 mJy. The pulsar periods range from 0.4 to 2.1 s. We report on candidate FRB 141113, which we argue is likely astrophysical and extragalactic, having $\rm DM \simeq 400\ pc~cm^{-3}$, which represents an excess over the Galactic maximum along this line of sight of $\rm \sim$ 100 - 200 pc cm$^{-3}$. We consider implications for the FRB population and show via simulations that if FRB 141113 is real and extragalactic, the slope $\alpha$ of the distribution of integral source counts as a function of flux density ($N (>S) \propto S^{-\alpha}$) is $1.4 \pm 0.5$ (95% confidence range). However this conclusion is dependent on several assumptions that require verification.

FRB 121102, the only repeating fast radio burst (FRB) known to date, was discovered at 1.4 GHz and shortly after the discovery of its repeating nature, detected up to 2.4 GHz. Here we present three bursts detected with the 100-m Effelsberg radio telescope at 4.85 GHz. All three bursts exhibited frequency structure on broad and narrow frequency scales. Using an autocorrelation function analysis, we measured a characteristic bandwidth of the small-scale structure of 6.4$\pm$1.6 MHz, which is consistent with the diffractive scintillation bandwidth for this line of sight through the Galactic interstellar medium (ISM) predicted by the NE2001 model. These were the only detections in a campaign totaling 22 hours in 10 observing epochs spanning five months. The observed burst detection rate within this observation was inconsistent with a Poisson process with a constant average occurrence rate; three bursts arrived in the final 0.3 hr of a 2 hr observation on 2016 August 20. We therefore observed a change in the rate of detectable bursts during this observation, and we argue that boosting by diffractive interstellar scintillations may have played a role in the detectability. Understanding whether changes in the detection rate of bursts from FRB 121102 observed at other radio frequencies and epochs are also a product of propagation effects, such as scintillation boosting by the Galactic ISM or plasma lensing in the host galaxy, or an intrinsic property of the burst emission will require further observations.

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}$.

The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fourier analyses (typically $\sim$16), or due to the strong effect of red noise at low modulation frequencies. We address this reduction in sensitivity by using a time-domain search technique: the Fast-Folding Algorithm (FFA). We designed a program that implements a FFA-based search in the PALFA processing pipeline, and tested the efficiency of the algorithm by performing tests under both ideal, white noise conditions, as well as with real PALFA observational data. In the two scenarios, we show that the time-domain algorithm has the ability to outperform the FFT-based periodicity search implemented in the survey. We perform simulations to compare the previously reported PALFA sensitivity with that obtained using our new FFA implementation. These simulations show that for a pulsar having a pulse duty cycle of roughly 3%, the performance of our FFA pipeline exceeds that of our FFT pipeline for pulses with DM $\lesssim$ 40 pc cm$^{-3}$ and for periods as short as $\sim$500 ms, and that the survey sensitivity is improved by at least a factor of two for periods $\gtrsim$ 6 sec. Discoveries from the implementation of the algorithm in PALFA are also presented in this paper.

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.

We report the discovery and initial follow-up of a double neutron star (DNS) system, PSR J1946$+$2052, with the Arecibo L-Band Feed Array pulsar (PALFA) survey. PSR J1946$+$2052 is a 17-ms pulsar in a 1.88-hour, eccentric ($e \, =\, 0.06$) orbit with a $\gtrsim 1.2 \, M_\odot$ companion. We have used the Jansky Very Large Array to localize PSR J1946$+$2052 to a precision of 0.09 arcseconds using a new phase binning mode. We have searched multiwavelength catalogs for coincident sources but did not find any counterparts. The improved position enabled a measurement of the spin period derivative of the pulsar ($\dot{P} \, = \, 9\,\pm \, 2 \,\times 10^{-19}$); the small inferred magnetic field strength at the surface ($B_S \, = \, 4 \, \times \, 10^9 \, \rm G$) indicates that this pulsar has been recycled. This and the orbital eccentricity lead to the conclusion that PSR J1946$+$2052 is in a DNS system. Among all known radio pulsars in DNS systems, PSR J1946$+$2052 has the shortest orbital period and the shortest estimated merger timescale, 46 Myr; at that time it will display the largest spin effects on gravitational wave waveforms of any such system discovered to date. We have measured the advance of periastron passage for this system, $\dot{\omega} \, = \, 25.6 \, \pm \, 0.3\, \deg \rm yr^{-1}$, implying a total system mass of only 2.50 $\pm$ 0.04 $M_\odot$, so it is among the lowest mass DNS systems. This total mass measurement combined with the minimum companion mass constrains the pulsar mass to $\lesssim 1.3 \, M_\odot$.

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.

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.

Fast radio bursts are a new class of transient radio phenomena currently detected as millisecond radio pulses with very high dispersion measures. As new radio surveys begin searching for FRBs a large population is expected to be detected in real-time, triggering a range of multi-wavelength and multi-messenger telescopes to search for repeating bursts and/or associated emission. Here we propose a method for disseminating FRB triggers using Virtual Observatory Events (VOEvents). This format was developed and is used successfully for transient alerts across the electromagnetic spectrum and for multi-messenger signals such as gravitational waves. In this paper we outline a proposed VOEvent standard for FRBs that includes the essential parameters of the event and where these parameters should be specified within the structure of the event. An additional advantage to the use of VOEvents for FRBs is that the events can automatically be ingested into the FRB Catalogue (FRBCAT) enabling real-time updates for public use. We welcome feedback from the community on the proposed standard outlined below and encourage those interested to join the nascent working group forming around this topic.

Pulsar timing arrays (PTAs) are used to search for long-wavelength gravitational waves (GWs) by monitoring a set of spin-stable millisecond pulsars. Most theoretical analyses assume that the relevant GW sources are much more distant from Earth than the pulsars comprising the array. Unlike ground- or solar system-based GW detectors, PTAs might well contain embedded GW sources. We derive the PTA response from sources at any distance, with a specific focus on GW bursts with memory (BWMs). We consider supernovae and compact binary mergers as potential Galactic BWM sources and evaluate the signature for an array with pulsars in globular clusters or in the Galactic center. Understanding the response of PTAs to nearby sources of BWM is a step towards investigating other more complex Galactic sources.

We undertook coordinated campaigns with the Green Bank, Effelsberg, and Arecibo radio telescopes during Chandra X-ray Observatory and XMM-Newton observations of the repeating fast radio burst FRB 121102 to search for simultaneous radio and X-ray bursts. We find 12 radio bursts from FRB 121102 during 70 ks total of X-ray observations. We detect no X-ray photons at the times of radio bursts from FRB 121102 and further detect no X-ray bursts above the measured background at any time. We place a 5$\sigma$ upper limit of $3\times10^{-11}$ erg cm$^{-2}$ on the 0.5--10 keV fluence for X-ray bursts at the time of radio bursts for durations $<700$ ms, which corresponds to a burst energy of $4\times10^{45}$ erg at the measured distance of FRB 121102. We also place limits on the 0.5--10 keV fluence of $5\times10^{-10}$ erg cm$^{-2}$ and $1\times10^{-9}$ erg cm$^{-2}$ for bursts emitted at any time during the XMM-Newton and Chandra observations, respectively, assuming a typical X-ray burst duration of 5 ms. We analyze data from the Fermi Gamma-ray Space Telescope Gamma-ray Burst Monitor and place a 5$\sigma$ upper limit on the 10--100 keV fluence of $4\times10^{-9}$ erg cm$^{-2}$ ($5\times10^{47}$ erg at the distance of FRB 121102) for gamma-ray bursts at the time of radio bursts. We also present a deep search for a persistent X-ray source using all of the X-ray observations taken to date and place a 5$\sigma$ upper limit on the 0.5--10 keV flux of $4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ ($3\times10^{41}$ erg~s$^{-1}$ at the distance of FRB 121102). We discuss these non-detections in the context of the host environment of FRB 121102 and of possible sources of fast radio bursts in general.

We present optical, near- and mid-infrared imaging of the host galaxy of FRB 121102 with the Gemini North telescope, the Hubble Space Telescope and the Spitzer Space Telescope. The FRB 121102 host galaxy is resolved, revealing a bright star forming region located in the outskirts of the irregular, low-metallicity dwarf galaxy. The star forming region has a half-light radius of 0.68 kpc (0.20 arcsec), encompassing the projected location of the compact (<0.7 pc), persistent radio source that is associated with FRB 121102. The half-light diameter of the dwarf galaxy is 5 to 7 kpc, and broadband spectral energy distribution fitting indicates that it has a total stellar mass of M*~10^8 Msun. The metallicity of the host galaxy is low, 12+log10 ([O/H])=8.0+-0.1. The properties of the host galaxy of FRB 121102 are comparable to those of extreme emission line galaxies, also known to host hydrogen-poor superluminous supernovae and long-duration gamma-ray bursts. The projected location of FRB 121102 within the star forming region supports the proposed connection of FRBs with newly born neutron stars or magnetars.

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.

While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here as "time-domain" observations. In this whitepaper we identify some of these policies, describe the scientific applications they impair, and suggest changes that would increase ALMA's science impact in Cycle 6 and beyond. Parties interested in time-domain science with ALMA are encouraged to join the ALMA Time-domain Special Interest Group (ATSIG) by signing up for the ATSIG mailing list at https://groups.google.com/group/alma-td-sig .

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.

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.

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.

We report on the discovery and timing observations of 29 distant long-period pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, confirmation and timing observations of these pulsars over several years at Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation and radiation properties. We have used multi-frequency data to measure the interstellar scattering properties of some of these pulsars. Most of the pulsars have properties that mirror those of the previously known pulsar population, although four show some notable characteristics. PSRs J1907+0631 and J1925+1720 are young and are associated with supernova remnants or plerionic nebulae: J1907+0631 lies close to the center of SNR G40.5-0.5, while J1925+1720 is coincident with a high-energy Fermi gamma-ray source. One pulsar, J1932+1500, is in a surprisingly eccentric, 199-day binary orbit with a companion having a minimum mass of 0.33 solar masses. Several of the sources exhibit timing noise, and two, PSRs J0611+1436 and J1907+0631, have both suffered large glitches, but with very different post-glitch rotation properties. In particular, the rotational period of PSR J0611+1436 will not recover to its pre-glitch value for about 12 years, a far greater recovery timescale than seen following any other large glitches.

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.

We report the discovery of two long-term intermittent radio pulsars in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, extended observations of these pulsars over several years at Jodrell Bank Observatory have revealed the details of their rotation and radiation properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bi-modal intermittency, switching between active (ON) and inactive (OFF) emission states and indicating the presence of a large, hitherto unrecognised, underlying population of such objects. For PSR J1929+1357, the initial duty cycle was fON=0.008, but two years later this changed, quite abruptly, to fON=0.16. This is the first time that a significant evolution in the activity of an intermittent pulsar has been seen and we show that the spin-down rate of the pulsar is proportional to the activity. The spin-down rate of PSR J1929+1357 is increased by a factor of 1.8 when it is in active mode, similar to the increase seen in the other three known long-term intermittent pulsars.

Gravitational wave astronomy using a pulsar timing array requires high-quality millisecond pulsars, correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for millisecond pulsars observed by the North American Nanohertz Observatory for Gravitational Waves. We characterize the excess noise using variance and structure function analyses. We find that 26 out of 37 pulsars show inconsistencies with a white-noise-only model based on the short timescale analysis of each pulsar and we demonstrate that the excess noise has a red power spectrum for 15 pulsars. We also decompose the excess noise into chromatic (radio-frequency-dependent) and achromatic components. Associating the achromatic red-noise component with spin noise and including additional power-spectrum-based estimates from the literature, we estimate a scaling law in terms of spin parameters (frequency and frequency derivative) and data-span length and compare it to the scaling law of Shannon \& Cordes (2010). We briefly discuss our results in terms of detection of gravitational waves at nanohertz frequencies.

In pulsar timing, timing residuals are the differences between the observed times of arrival and the predictions from the timing model. A comprehensive timing model will produce featureless residuals, which are presumably composed of dominating noise and weak physical effects excluded from the timing model (e.g. gravitational waves). In order to apply the optimal statistical methods for detecting the weak gravitational wave signals, we need to know the statistical properties of the noise components in the residuals. In this paper we utilize a variety of non-parametric statistical tests to analyze the whiteness and Gaussianity of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 5-year timing data which are obtained from the Arecibo Observatory and the Green Bank Telescope from 2005 to 2010 (Demorest et al. 2013). We find that most of the data are consistent with white noise; Many data deviate from Gaussianity at different levels, nevertheless, removing outliers in some pulsars will mitigate the deviations.

We report on radio and X-ray observations of the only known repeating Fast Radio Burst (FRB) source, FRB 121102. We have detected six additional radio bursts from this source: five with the Green Bank Telescope at 2 GHz, and one at 1.4 GHz at the Arecibo Observatory for a total of 17 bursts from this source. All have dispersion measures consistent with a single value ($\sim559$ pc cm$^{-3}$) that is three times the predicted maximum Galactic value. The 2-GHz bursts have highly variable spectra like those at 1.4 GHz, indicating that the frequency structure seen across the individual 1.4 and 2-GHz bandpasses is part of a wideband process. X-ray observations of the FRB 121102 field with the Swift and Chandra observatories show at least one possible counterpart; however, the probability of chance superposition is high. A radio imaging observation of the field with the Jansky Very Large Array at 1.6 GHz yields a 5$\sigma$ upper limit of 0.3 mJy on any point-source continuum emission. This upper limit, combined with archival WISE 22-$\mu$m and IPHAS H$\alpha$ surveys, rules out the presence of an intervening Galactic HII region. We update our estimate of the FRB detection rate in the PALFA survey to be 1.1$^{+3.7}_{-1.0} \times 10^4$ FRBs sky$^{-1}$ day$^{-1}$ (95% confidence) for peak flux density at 1.4 GHz above 300 mJy. We find that the intrinsic widths of the 12 FRB 121102 bursts from Arecibo are, on average, significantly longer than the intrinsic widths of the 13 single-component FRBs detected with the Parkes telescope.

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.

We analyze deterministic and random temporal variations in dispersion measure (DM) from the full three-dimensional velocities of pulsars with respect to the solar system, combined with electron-density variations on a wide range of length scales. Previous treatments have largely ignored the pulsar's changing distance while favoring interpretations involving the change in sky position from transverse motion. Linear trends in pulsar DMs seen over 5-10~year timescales may signify sizable DM gradients in the interstellar medium (ISM) sampled by the changing direction of the line of sight to the pulsar. We show that motions parallel to the line of sight can also account for linear trends, for the apparent excess of DM variance over that extrapolated from scintillation measurements, and for the apparent non-Kolmogorov scalings of DM structure functions inferred in some cases. Pulsar motions through atomic gas may produce bow-shock ionized gas that also contributes to DM variations. We discuss possible causes of periodic or quasi-periodic changes in DM, including seasonal changes in the ionosphere, annual variation of the solar elongation angle, structure in the heliosphere-ISM boundary, and substructure in the ISM. We assess the solar cycle's role on the amplitude of ionospheric and solar-wind variations. Interstellar refraction can produce cyclic timing variations from the error in transforming arrival times to the solar system barycenter. We apply our methods to DM time series and DM gradient measurements in the literature and assess consistency with a Kolmogorov medium. Finally, we discuss the implications of DM modeling in precision pulsar timing experiments.

The use of pulsars as astrophysical clocks for gravitational wave experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for Gravitational Waves to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of RMS jitter in our pulsars to be $\sim 1\%$ of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve gravitational wave detection pipelines.

To obtain the most accurate pulse arrival times from radio pulsars, it is necessary to correct or mitigate the effects of the propagation of radio waves through the warm and ionised interstellar medium. We examine both the strength of propagation effects associated with large-scale electron-density variations and the methodology used to estimate infinite-frequency arrival times. Using simulations of two-dimensional phase-varying screens, we assess the strength and non-stationarity of timing perturbations associated with large-scale density variations. We identify additional contributions to arrival times that are stochastic in both radio frequency and time and therefore not amenable to correction solely using times of arrival. We attribute this to the frequency dependence of the trajectories of the propagating radio waves. We find that this limits the efficacy of low-frequency (metre-wavelength) observations. Incorporating low-frequency pulsar observations into precision timing campaigns is increasingly problematic for pulsars with larger dispersion measures.

We report the discovery and timing results for five millisecond pulsars (MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726, J1938+2516, and J1957+2516. Timing observations of the 5 pulsars were conducted with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3 yr. All of the MSPs except one (PSR J1914+0659) are in binary systems with low eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ~0.1 $M_\odot$ companion and possible eclipses that last ~10% of the orbit. The position of PSR J1957+2516 is also coincident with a NIR source. All 5 MSPs are distant (>3.1 kpc) as determined from their dispersion measures, and none of them show evidence of $\gamma$-ray pulsations in a search of Fermi Gamma-Ray Space Telescope data. These 5 MSPs bring the total number of MSPs discovered by the PALFA survey to 26 and further demonstrate the power of this survey in finding distant, highly dispersed MSPs deep in the Galactic plane.

We report here the Einstein@Home discovery of PSR J1913+1102, a 27.3-ms pulsar found in data from the ongoing Arecibo PALFA pulsar survey. The pulsar is in a 4.95-hr double neutron star (DNS) system with an eccentricity of 0.089. From radio timing with the Arecibo 305-m telescope, we measure the rate of advance of periastron to be 5.632(18) deg/yr. Assuming general relativity accurately models the orbital motion, this corresponds to a total system mass of 2.875(14) solar masses, similar to the mass of the most massive DNS known to date, B1913+16, but with a much smaller eccentricity. The small eccentricity indicates that the second-formed neutron star (the companion of PSR J1913+1102) was born in a supernova with a very small associated kick and mass loss. In that case this companion is likely, by analogy with other systems, to be a light (1.2 solar mass) neutron star; the system would then be highly asymmetric. A search for radio pulsations from the companion yielded no plausible detections, so we can't yet confirm this mass asymmetry. By the end of 2016, timing observations should permit the detection of two additional post-Keplerian parameters: the Einstein delay, which will enable precise mass measurements and a verification of the possible mass asymmetry of the system, and the orbital decay due to the emission of gravitational waves, which will allow another test of the radiative properties of gravity. The latter effect will cause the system to coalesce in ~0.5 Gyr.

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.

We analyze plasma dispersion and scattering of fast radio bursts (FRBs) to identify the dominant locations of free electrons along their lines of sight and thus constrain the distances of the burst sources themselves. We establish the average $\tau$-DM relation for Galactic pulsars and use it as a benchmark for discussing FRB scattering. Though scattering times $\tau$ for FRBs are large in the majority of the 17 events we analyze, they are systematically smaller than those of Galactic pulsars that have similar dispersion measures (DMs). The lack of any correlation between $\tau$ and DM for FRBs suggests that the intergalactic medium (IGM) cannot account for both $\tau$ and DM. We therefore consider mixed models involving the IGM and host galaxies. If the IGM contributes significantly to DM while host galaxies dominate $\tau$, the scattering deficit with respect to the mean Galactic trend can be explained with a $\tau$-DM relation in the host that matches that for the Milky Way. However, it is possible that hosts dominate both $\tau$ and DM, in which case the observed scattering deficits require free electrons in the host to be less turbulent than in the Galaxy, such as if they are in hot rather than warm ionized regions. Our results imply that distances or redshifts of FRB sources can be significantly overestimated if they are based on the assumption that the extragalactic portion of DM is dominated by the IGM.

Fast Radio Bursts are millisecond-duration astronomical radio pulses of unknown physical origin that appear to come from extragalactic distances. Previous follow-up observations have failed to find additional bursts at the same dispersion measures (i.e. integrated column density of free electrons between source and telescope) and sky position as the original detections. The apparent non-repeating nature of the fast radio bursts has led several authors to hypothesise that they originate in cataclysmic astrophysical events. Here we report the detection of ten additional bursts from the direction of FRB121102, using the 305-m Arecibo telescope. These new bursts have dispersion measures and sky positions consistent with the original burst. This unambiguously identifies FRB121102 as repeating and demonstrates that its source survives the energetic events that cause the bursts. Additionally, the bursts from FRB121102 show a wide range of spectral shapes that appear to be predominantly intrinsic to the source and which vary on timescales of minutes or shorter. While there may be multiple physical origins for the population of fast radio bursts, the repeat bursts with high dispersion measure and variable spectra specifically seen from FRB121102 support models that propose an origin in a young, highly magnetised, extragalactic neutron star.

In a search of the unidentified Fermi gamma-ray source 3FGL J1417.5-4402 with the Parkes radio telescope, we discovered PSR J1417-4402, a 2.66 ms pulsar having the same 5.4 day orbital period as the optical and X-ray binary identified by Strader et al. The existence of radio pulsations implies that the neutron star is currently not accreting. Substantial outflows from the companion render the radio pulsar undetectable for more than half of the orbit, and may contribute to the observed Halpha emission. Our initial pulsar observations, together with the optically inferred orbit and inclination, imply a mass ratio of 0.171+/-0.002, a companion mass of M_2=0.33+/-0.03 Msun, and a neutron star mass in the range 1.77<M_1<2.13 Msun. However, there remains a discrepancy between the distance of 4.4 kpc inferred from the optical properties of the companion and the smaller radio dispersion measure distance of 1.6 kpc. The smaller distance would reduce the inferred Roche-lobe filling factor, increase the inferred inclination angle, and decrease the masses. As a wide binary, PSR J1417-4402 differs from the radio-eclipsing black widow and redback pulsars being discovered in large numbers by Fermi. It is probably a system that began mass transfer onto the neutron star after the companion star left the main sequence. The companion should end its evolution as a He white dwarf in a 6-20 day orbit, i.e., as a typical binary millisecond pulsar companion.

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.

We consider radio bursts that originate from extragalactic neutron stars (NSs) by addressing three questions about source distances. What are the physical limitations on coherent radiation at GHz frequencies? Do they permit detection at cosmological distances? How many bursts per NS are needed to produce the inferred burst rate $\sim 10^3$-$10^4 $sky$^{-1}$ day$^{-1}$? The burst rate is comparable to the NS formation rate in a Hubble volume, requiring only one per NS if they are bright enough. However, radiation physics causes us to favor a closer population. More bursts per NS are then required but repeats in 10 to 100 yr could still be negligible. Bursts are modeled as sub-ns, coherent shot pulses superposed incoherently to produce ms-duration $\sim 1$ Jy amplitudes; each shot-pulse can be much weaker than the burst amplitude, placing less restrictive requirements on the emission process. Nonetheless, single shot pulses are similar to the extreme, unresolved ($< 0.4$ ns) MJy shot pulse seen from the Crab pulsar, which is consistent with coherent curvature radiation emitted near the light cylinder by an almost neutral clump with net charge $\sim \pm 10^{21}e$ and total energy $\gtrsim 10^{23}$ ergs. Bursts from Gpc distances require incoherent superposition of $\sim 10^{12}d_{\rm Gpc}^2$ shot pulses or a total energy $\gtrsim 10^{35} d_{\rm Gpc}^2$ erg. The energy reservoir near the light cylinder limits the detection distance to $\lesssim {\rm few} \times 100$ Mpc for a fluence $\sim 1$ Jy ms unless conditions are more extreme than for the Crab pulsar. Similarly, extreme single pulses from ordinary pulsars and magnetars could be detectable from throughout the Local Group and perhaps farther. Contributions to dispersion measures from galaxy clusters will be significant for some of the bursts. We discuss tests for the signatures of bursts associated with extragalactic NSs.

We report here on key science topics for the Next Generation Very Large Array in the areas of time domain, fundamental physics, and cosmology. Key science cases considered are pulsars in orbit around the Galactic Center massive black hole, Sagittarius A*, electromagnetic counterparts to gravitational waves, and astrometric cosmology. These areas all have the potential for ground-breaking and transformative discovery. Numerous other topics were discussed during the preparation of this report and some of those discussions are summarized here, as well. There is no doubt that further investigation of the science case will reveal rich and compelling opportunities.

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.

We report the discovery of the millisecond pulsar (MSP) PSR J1950+2414 ($P=4.3$ ms) in a binary system with an eccentric ($e=0.08$) 22-day orbit in Pulsar ALFA survey observations with the Arecibo telescope. Its companion star has a median mass of 0.3 $M_\odot$ and is most likely a white dwarf. Fully recycled MSPs like this one are thought to be old neutron stars spun-up by mass transfer from a companion star. This process should circularize the orbit, as is observed for the vast majority of binary MSPs, which predominantly have orbital eccentricities $e < 0.001$. However, four recently discovered binary MSPs have orbits with $0.027 < e < 0.44$; PSR J1950+2414 is the fifth such system to be discovered. The upper limits for its intrinsic spin period derivative and inferred surface magnetic field strength are comparable to those of the general MSP population. The large eccentricities are incompatible with the predictions of the standard recycling scenario: something unusual happened during their evolution. Proposed scenarios are a) initial evolution of the pulsar in a triple system which became dynamically unstable, b) origin in an exchange encounter in an environment with high stellar density, c) rotationally delayed accretion-induced collapse of a super-Chandrasekhar white dwarf, and d) dynamical interaction of the binary with a circumbinary disk. We compare the properties of all five known eccentric MSPs with the predictions of these formation channels. Future measurements of the masses and proper motion might allow us to firmly exclude some of the proposed formation scenarios.

The on-going PALFA survey at the Arecibo Observatory began in 2004 and is searching for radio pulsars in the Galactic plane at 1.4 GHz. Observations since 2009 have been made with new wider-bandwidth spectrometers than were previously employed in this survey. A new data reduction pipeline has been in place since mid-2011 which consists of standard methods using dedispersion, searches for accelerated periodic sources, and search for single pulses, as well as new interference-excision strategies and candidate selection heuristics. This pipeline has been used to discover 41 pulsars, including 8 millisecond pulsars (MSPs; P < 10 ms), bringing the PALFA survey's discovery totals to 145 pulsars, including 17 MSPs, and one Fast Radio Burst (FRB). The pipeline presented here has also re-detected 188 previously known pulsars including 60 found in PALFA data by re-analyzing observations previously searched by other pipelines. A comprehensive description of the survey sensitivity, including the effect of interference and red noise, has been determined using synthetic pulsar signals with various parameters and amplitudes injected into real survey observations and subsequently recovered with the data reduction pipeline. We have confirmed that the PALFA survey achieves the sensitivity to MSPs predicted by theoretical models. However, we also find that compared to theoretical survey sensitivity models commonly used there is a degradation in sensitivity to pulsars with periods P >= 100 ms that gradually becomes up to a factor of ~10 worse for P > 4 s at DM < 150 pc/cc. This degradation of sensitivity at long periods is largely due to red noise. We find that 35 +- 3% of pulsars are missed despite being bright enough to be detected in the absence of red noise. This reduced sensitivity could have implications on estimates of the number of long-period pulsars in the Galaxy.

We analyze the frequency dependence of the dispersion measure (DM), the column density of free electrons to a pulsar, caused by multipath scattering from small scale electron-density fluctuations. The DM is slightly different along each propagation path and the transverse spread of paths varies greatly with frequency, yielding time-of-arrival (TOA) perturbations that scale differently than the inverse square of the frequency, the expected dependence for a cold, unmagnetized plasma. We quantify DM and TOA perturbations analytically for thin phase screens and extended media and verify the results with simulations of thin screens. The rms difference between DMs across an octave band near 1.5~GHz $\sim 4\times10^{-5}\,{\rm pc\ cm^{-3}}$ for pulsars at $\sim 1$~kpc distance. TOA errors from chromatic DMs are of order a few to hundreds of nanoseconds for pulsars with DM $\lesssim 30$~pc~cm$^{-3}$ observed across an octave band but increase rapidly to microseconds or larger for larger DMs and wider frequency ranges. Frequency-dependent DMs introduce correlated noise into timing residuals whose power spectrum is `low pass' in form. The correlation time is of order the geometric mean of the refraction times for the highest and lowest radio frequencies used and thus ranges from days to years, depending on the pulsar. We discuss the implications for methodologies that use large frequency separations or wide bandwidth receivers for timing measurements. Chromatic DMs are partially mitigable by using an additional chromatic term in arrival time models. Without mitigation, our results provide an additional term in the noise model for pulsar timing; they also indicate that in combination with measurement errors from radiometer noise, an arbitrary increase in total frequency range (or bandwidth) will yield diminishing benefits and may be detrimental to overall timing precision.

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}$.

Free electrons in the interstellar medium cause frequency-dependent delays in pulse arrival times due to both scattering and dispersion. Multi-frequency measurements are used to estimate and remove dispersion delays. In this paper, we focus on the effect of any non-simultaneity of multi-frequency observations on dispersive delay estimation and removal. Interstellar density variations combined with changes in the line-of-sight from pulsar and observer motions cause dispersion measure variations with an approximately power-law power spectrum, augmented in some cases by linear trends. We simulate time series, estimate the magnitude and statistical properties of timing errors that result from non-simultaneous observations, and derive prescriptions for data acquisition that are needed in order to achieve a specified timing precision. For nearby, highly stable pulsars, measurements need to be simultaneous to within about one day in order that the timing error from asynchronous DM correction is less than about 10 ns. We discuss how timing precision improves when increasing the number of dual-frequency observations used in dispersion measure estimation for a given epoch. For a Kolmogorov wavenumber spectrum, we find about a factor of two improvement in precision timing when increasing from two to three observations but diminishing returns thereafter.

We present the discovery of five millisecond pulsars (MSPs) from the PALFA Galactic plane survey using Arecibo. Four of these (PSRs J0557+1551, J1850+0244, J1902+0300, and J1943+2210) are binary pulsars whose companions are likely white dwarfs, and one (PSR J1905+0453) is isolated. Phase-coherent timing solutions, ranging from $\sim$1 to $\sim$3 years in length, and based on observations from the Jodrell Bank and Arecibo telescopes, provide precise determinations of spin, orbital, and astrometric parameters. All five pulsars have large dispersion measures ($>100$ pc cm$^{-3}$, within the top 20% of all known Galactic field MSPs) and are faint (1.4 GHz flux density < 0.1 mJy, within the faintest 5% of all known Galactic field MSPs), illustrating PALFA's ability to find increasingly faint, distant MSPs in the Galactic plane. In particular, PSR J1850+0244 has a dispersion measure of 540 pc cm$^{-3}$, the highest of all known MSPs. Such distant, faint MSPs are important input for accurately modeling the total Galactic MSP population.

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.

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 vast collecting area of the Square Kilometre Array (SKA), harnessed by sensitive receivers, flexible digital electronics and increased computational capacity, could permit the most sensitive and exhaustive search for technologically-produced radio emission from advanced extraterrestrial intelligence (SETI) ever performed. For example, SKA1-MID will be capable of detecting a source roughly analogous to terrestrial high-power radars (e.g. air route surveillance or ballistic missile warning radars, EIRP (EIRP = equivalent isotropic radiated power, ~10^17 erg sec^-1) at 10 pc in less than 15 minutes, and with a modest four beam SETI observing system could, in one minute, search every star in the primary beam out to ~100 pc for radio emission comparable to that emitted by the Arecibo Planetary Radar (EIRP ~2 x 10^20 erg sec^-1). The flexibility of the signal detection systems used for SETI searches with the SKA will allow new algorithms to be employed that will provide sensitivity to a much wider variety of signal types than previously searched for. Here we discuss the astrobiological and astrophysical motivations for radio SETI and describe how the technical capabilities of the SKA will explore the radio SETI parameter space. We detail several conceivable SETI experimental programs on all components of SKA1, including commensal, primary-user, targeted and survey programs and project the enhancements to them possible with SKA2. We also discuss target selection criteria for these programs, and in the case of commensal observing, how the varied use cases of other primary observers can be used to full advantage for SETI.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The Pulsar Arecibo L-band Feed Array (PALFA) Survey uses the ALFA 7-beam receiver to search both inner and outer Galactic sectors visible from Arecibo ($32^{\circ}\lesssim \ell \lesssim 77^{\circ}$ and $168^{\circ}\lesssim \ell \lesssim 214^{\circ}$) close to the Galactic plane ($|b|\lesssim5^{\circ}$) for pulsars. In this paper we detail a precursor survey of this region with PALFA, which observed a subset of the full region (slightly more restrictive in $\ell$ and $|b|\lesssim1^{\circ}$) and detected 45 pulsars. For both Galactic millisecond and normal pulsar populations, we compare the survey's detections with simulations to model these populations and, in particular, to estimate the number of observable pulsars in the Galaxy. We place 95\% confidence intervals of 82,000 to 143,000 on the number of detectable normal pulsars and 9,000 to 100,000 on the number of detectable millisecond pulsars in the Galactic disk. These are consistent with previous estimates. Given the most likely population size in each case (107,000 and 15,000 for normal and millisecond pulsars, respectively) we extend survey detection simulations to predict that, when complete, the full PALFA survey should have detected $1,000\substack{+330 \\ -230}$ normal pulsars and $30\substack{+200 \\ -20}$ millisecond pulsars. Identical estimation techniques predict that $490\substack{+160 \\ -115}$ normal pulsars and $12\substack{+70 \\ -5}$ millisecond pulsars would be detected by the beginning of 2014; at the time, the PALFA survey had detected 283 normal pulsars and 31 millisecond pulsars, respectively.

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.

Highly energetic astrophysical phenomena like supermassive black hole binary (SMBHB) mergers are predicted to emit prodigious amounts of gravitational waves (GWs). An anticipated component of the gravitational waveform known as "memory" is permanent and non-oscillatory. For SMBHB mergers, the memory is created primarily during the most violent moments of the inspiral immediately preceding the final plunge and ring-down when the strongest gravitational fields are at work and the non-linearities of general relativity are most pronounced. The essentially time-domain nature of memory makes it forbiddingly difficult to detect with ground based GW detectors, leaving pulsar timing array (PTA) experiments as the most promising means by which it may be detected and studied. In this paper, we discuss how GW bursts with memory (BWMs) influence pulsar timing experiments and develop methods to assess how sensitive modern timing efforts are to such GW events. We discuss how PTA searches for BWMs can be used to constrain the rate of BWMs and how these constraints relate to information regarding the population of SMBHBs.

Recent work has exploited pulsar survey data to identify temporally isolated, millisecond-duration radio bursts with large dispersion measures (DMs). These bursts have been interpreted as arising from a population of extragalactic sources, in which case they would provide unprecedented opportunities for probing the intergalactic medium; they may also be linked to new source classes. Until now, however, all so-called fast radio bursts (FRBs) have been detected with the Parkes radio telescope and its 13-beam receiver, casting some concern about the astrophysical nature of these signals. Here we present FRB 121102, the first FRB discovery from a geographic location other than Parkes. FRB 121102 was found in the Galactic anti-center region in the 1.4-GHz Pulsar ALFA survey with the Arecibo Observatory with a DM = 557.4 $\pm$ 3 pc cm$^{-3}$, pulse width of $3\; \pm 0.5$ ms, and no evidence of interstellar scattering. The observed delay of the signal arrival time with frequency agrees precisely with the expectation of dispersion through an ionized medium. Despite its low Galactic latitude ($b = -0.2^{\circ}$), the burst has three times the maximum Galactic DM expected along this particular line-of-sight, suggesting an extragalactic origin. A peculiar aspect of the signal is an inverted spectrum; we interpret this as a consequence of being detected in a sidelobe of the ALFA receiver. FRB 121102's brightness, duration, and the inferred event rate are all consistent with the properties of the previously detected Parkes bursts.

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.

In this paper, we present a novel artificial intelligence (AI) program that identifies pulsars from recent surveys using image pattern recognition with deep neural nets---the PICS (Pulsar Image-based Classification System) AI. The AI mimics human experts and distinguishes pulsars from noise and interferences by looking for patterns from candidate. The information from each pulsar candidate is synthesized in four diagnostic plots, which consist of up to thousands pixel of image data. The AI takes these data from each candidate as its input and uses thousands of such candidates to train its ~9000 neurons. Different from other pulsar selection programs which use pre-designed patterns, the PICS AI teaches itself the salient features of different pulsars from a set of human-labeled candidates through machine learning. The deep neural networks in this AI system grant it superior ability in recognizing various types of pulsars as well as their harmonic signals. The trained AI's performance has been validated with a large set of candidates different from the training set. In this completely independent test, PICS ranked 264 out of 277 pulsar-related candidates, including all 56 previously known pulsars, to the top 961 (1%) of 90008 test candidates, missing only 13 harmonics. The first non-pulsar candidate appears at rank 187, following 45 pulsars and 141 harmonics. In other words, 100% of the pulsars were ranked in the top 1% of all candidates, while 80% were ranked higher than any noise or interference. The performance of this system can be improved over time as more training data are accumulated. This AI system has been integrated into the PALFA survey pipeline and has discovered six new pulsars to date.

Precision astrometry is an integral component of successful pulsar timing campaigns. Astrometric parameters are commonly derived by fitting them as parameters of a timing model to a series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over several-year spans can yield position measurements with sub-milliarcsecond precision. However, timing-based astrometry can become biased if a pulsar displays any red spin noise, which can be compared to the red noise signal produced by the stochastic gravitational wave background. We investigate how noise of different spectral types is absorbed by timing models, leading to significant estimation errors in the astrometric parameters. We find that commonly used techniques for fitting timing models in the presence of red noise (Cholesky whitening) prevent the absorption of noise into the timing model remarkably well if the time baseline of observations exceeds several years, but are inadequate for dealing with shorter pulsar data sets. Independent of timing, pulsar-optimized very long baseline interferometry (VLBI) is capable of providing position estimates precise to the sub-milliarcsecond levels needed for high-precision timing. We compute a necessary transformation between the International Celestial Reference Frame and pulsar timing frames and quantitatively discuss how the transformation will improve in coming years. We find that incorporating VLBI astrometry into the timing models of pulsars for which only a couple of years of timing data exist will lead to more realistic assessments of red spin noise and could enhance the amplitude of gravitational wave signatures in post-fit timing residuals by factors of 20 or more.

We present temporal scattering measurements of single pulses and average profiles of PSR J1745--2900, a magnetar recently discovered only 3 arcsec away from Sagittarius A* (Sgr A*), from 1.2 - 18.95 GHz using the Effelsberg 100-m Radio Telescope, the Nan\c{c}ay Decimetric Radio Telescope, and the Jodrell Bank Lovell Telescope. Single pulse analysis shows that the integrated pulse profile above 2 GHz is dominated by pulse jitter, while below 2 GHz the pulse profile shape is dominated by scattering. The high dispersion measure and rotation measure of the magnetar suggest that it is close to Sgr A* (within ~0.1 pc). This is the first object in the GC with both pulse broadening and angular broadening measurements. We measure a pulse broadening spectral index of alpha = -3.8 +/- 0.2 and a pulse broadening time scale at 1 GHz of tau_GHz = 1.3 +/- 0.2 s, which is several orders of magnitude lower than the scattering predicted by the NE2001 model (Cordes and Lazio 2002). If this scattering timescale is representative of the GC as a whole, then previous surveys should have detected many pulsars. The lack of detections implies either our understanding of scattering in the GC is incomplete or there are fewer pulsars in the GC than previously predicted. Given that magnetars are a rare class of radio pulsar, we believe that there are many canonical and millisecond pulsars in the GC, and not surprisingly, scattering regions in the GC have complex spatial structures.

We use Fermi Gamma-ray Space Telescope detections and upper limits on non-recycled pulsars obtained from the Large Area Telescope (LAT) to constrain how the gamma-ray luminosity L depends on the period P and the period derivative \dot{P}. We use a Bayesian analysis to calculate a best-fit luminosity law, or dependence of L on P and \dot{P}, including different methods for modeling the beaming factor. An outer gap (OG) magnetosphere geometry provides the best-fit model, which is L \propto P^{-a} \dot{P}^{b} where a=1.36\pm0.03 and b=0.44\pm0.02, similar to but not identical to the commonly assumed L \propto \sqrt{\dot{E}} \propto P^{-1.5} \dot{P}^{0.5}. Given upper limits on gamma-ray fluxes of currently known radio pulsars and using the OG model, we find that about 92% of the radio-detected pulsars have gamma-ray beams that intersect our line of sight. By modeling the misalignment of radio and gamma-ray beams of these pulsars, we find an average gamma-ray beaming solid angle of about 3.7{\pi} for the OG model, assuming a uniform beam. Using LAT-measured diffuse fluxes, we place a 2{\sigma} upper limit on the average braking index and a 2{\sigma} lower limit on the average surface magnetic field strength of the pulsar population of 3.8 and 3.2 X 10^{10} G, respectively. We then predict the number of non-recycled pulsars detectable by the LAT based on our population model. Using the two-year sensitivity, we find that the LAT is capable of detecting emission from about 380 non-recycled pulsars, including 150 currently identified radio pulsars. Using the expected five-year sensitivity, about 620 non-recycled pulsars are detectable, including about 220 currently identified radio pulsars. We note that these predictions significantly depend on our model assumptions.

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 193 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using Pulsar ALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 minutes. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722's pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~ 100 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.

We have made extensive observations of 35 distant slow (non-recycled) pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Timing observations of these pulsars over several years at Arecibo Observatory and Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation properties. Despite being a relatively distant population, these pulsars have properties that mirror those of the previously known pulsar population. Many of the sources exhibit timing noise, and one underwent a small glitch. We have used multifrequency data to measure the interstellar scattering properties of these pulsars. We find scattering to be higher than predicted along some lines of sight, particularly in the Cygnus region. Lastly, we present XMM-Newton and Chandra observations of the youngest and most energetic of the pulsars, J1856+0245, which has previously been associated with the GeV-TeV pulsar wind nebula HESS J1857+026.

Modern radio pulsar surveys produce a large volume of prospective candidates, the majority of which are polluted by human-created radio frequency interference or other forms of noise. Typically, large numbers of candidates need to be visually inspected in order to determine if they are real pulsars. This process can be labor intensive. In this paper, we introduce an algorithm called PEACE (Pulsar Evaluation Algorithm for Candidate Extraction) which improves the efficiency of identifying pulsar signals. The algorithm ranks the candidates based on a score function. Unlike popular machine-learning based algorithms, no prior training data sets are required. This algorithm has been applied to data from several large-scale radio pulsar surveys. Using the human-based ranking results generated by students in the Arecibo Remote Command enter programme, the statistical performance of PEACE was evaluated. It was found that PEACE ranked 68% of the student-identified pulsars within the top 0.17% of sorted candidates, 95% within the top 0.34%, and 100% within the top 3.7%. This clearly demonstrates that PEACE significantly increases the pulsar identification rate by a factor of about 50 to 1000. To date, PEACE has been directly responsible for the discovery of 47 new pulsars, 5 of which are millisecond pulsars that may be useful for pulsar timing based gravitational-wave detection projects.

Markov processes are shown to be consistent with metastable states seen in pulsar phenomena, including intensity nulling, pulse-shape mode changes, subpulse drift rates, spindown rates, and X-ray emission, based on the typically broad and monotonic distributions of state lifetimes. Markovianity implies a nonlinear magnetospheric system in which state changes occur stochastically, corresponding to transitions between local minima in an effective potential. State durations (though not transition times) are thus largely decoupled from the characteristic time scales of various magnetospheric processes. Dyadic states are common but some objects show at least four states with some transitions forbidden. Another case is the long-term intermittent pulsar B1931+24 that has binary radio-emission and torque states with wide, but non-monotonic duration distributions. It also shows a quasi-period of $38\pm5$ days in a 13-yr time sequence, suggesting stochastic resonance in a Markov system with a forcing function that could be strictly periodic or quasi-periodic. Nonlinear phenomena are associated with time-dependent activity in the acceleration region near each magnetic polar cap. The polar-cap diode is altered by feedback from the outer magnetosphere and by return currents from an equatorial disk that may also cause the neutron star to episodically charge and discharge. Orbital perturbations in the disk provide a natural periodicity for the forcing function in the stochastic resonance interpretation of B1931+24. Disk dynamics may introduce additional time scales in observed phenomena. Future work can test the Markov interpretation, identify which pulsar types have a propensity for state changes, and clarify the role of selection effects.

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.

The on-going PALFA survey is searching the Galactic plane (|b| < 5 deg., 32 < l < 77 deg. and 168 < l < 214 deg.) for radio pulsars at 1.4 GHz using ALFA, the 7-beam receiver installed at the Arecibo Observatory. By the end of August 2012, the PALFA survey has discovered 100 pulsars, including 17 millisecond pulsars (P < 30 ms). Many of these discoveries are among the pulsars with the largest DM/P ratios, proving that the PALFA survey is capable of probing the Galactic plane for millisecond pulsars to a much greater depth than any previous survey. This is due to the survey's high sensitivity, relatively high observing frequency, and its high time and frequency resolution. Recently the rate of discoveries has increased, due to a new more sensitive spectrometer, two updated complementary search pipelines, the development of online collaborative tools, and access to new computing resources. Looking forward, focus has shifted to the application of artificial intelligence systems to identify pulsar-like candidates, and the development of an improved full-resolution pipeline incorporating more sophisticated radio interference rejection. The new pipeline will be used in a complete second analysis of data already taken, and will be applied to future survey observations. An overview of recent developments, and highlights of exciting discoveries will be presented.

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.

Most millisecond pulsars, like essentially all other radio pulsars, show timing errors well in excess of what is expected from additive radiometer noise alone. We show that changes in amplitude, shape and pulse phase for the millisecond pulsar J1713+0747 cause this excess error. These changes appear to be uncorrelated from one pulse period to the next. The resulting time of arrival variations are correlated across a wide frequency range and are observed with different backend processors on different days, confirming that they are intrinsic in origin and not an instrumental effect or caused by strongly frequency dependent interstellar scattering. Centroids of single pulses show an rms phase variation \approx 40 microsec, which dominates the timing error and is the same phase jitter phenomenon long known in slower spinning, canonical pulsars. We show that the amplitude modulations of single pulses are modestly correlated with their arrival time fluctuations. We also demonstrate that single-pulse variations are completely consistent with arrival time variations of pulse profiles obtained by integrating N pulses such that the arrival time error decreases proportional to 1/\sqrt{N}. We investigate methods for correcting times of arrival for these pulse shape changes, including multi-component TOA fitting and principal component analysis. These techniques are not found to improve the timing precision of the observations. We conclude that when pulse shape changes dominate timing errors, the timing precision of PSR J1713+0747 can be only improved by averaging over a larger number of pulses.

We present the discovery and phase-coherent timing of four highly dispersed millisecond pulsars (MSPs) from the Arecibo PALFA Galactic plane survey: PSRs J1844+0115, J1850+0124, J1900+0308, and J1944+2236. Three of the four pulsars are in binary systems with low-mass companions, which are most likely white dwarfs, and which have orbital periods on the order of days. The fourth pulsar is isolated. All four pulsars have large dispersion measures (DM > 100 pc cm-3), are distant (> 3.4 kpc), faint at 1.4 GHz (< 0.2 mJy), and are fully recycled (with spin periods P between 3.5 and 4.9 ms). The three binaries also have very small orbital eccentricities, as expected for tidally circularized, fully recycled systems with low-mass companions. These four pulsars have DM/P ratios that are among the highest values for field MSPs in the Galaxy. These discoveries bring the total number of confirmed MSPs from the PALFA survey to fifteen. The discovery of these MSPs illustrates the power of PALFA for finding weak, distant MSPs at low-Galactic latitudes. This is important for accurate estimates of the Galactic MSP population and for the number of MSPs that the Square Kilometer Array can be expected to detect.

We present two millisecond pulsar discoveries from the PALFA survey of the Galactic plane with the Arecibo telescope. PSR J1955+2527 is an isolated pulsar with a period of 4.87 ms, and PSR J1949+3106 has a period of 13.14 ms and is in a 1.9-day binary system with a massive companion. Their timing solutions, based on 4 years of timing measurements with the Arecibo, Green Bank, Nan\c{c}ay and Jodrell Bank telescopes, allow precise determination of spin and astrometric parameters, including precise determinations of their proper motions. For PSR J1949+3106, we can clearly detect the Shapiro delay. From this we measure the pulsar mass to be 1.47(+0.43/-0.31) solar masses, the companion mass to be 0.85(+0.14/-0.11) solar masses and the orbital inclination to be i = 79.9(+1.6/-1.9) degrees, where uncertainties correspond to +/- 1-\sigma\ confidence levels. With continued timing, we expect to also be able to detect the advance of periastron for the J1949+3106 system. This effect, combined with the Shapiro delay, will eventually provide very precise mass measurements for this system and a test of general relativity.

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.

We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7x10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.

We present a multimoment technique for signal classification and apply it to the detection of fast radio transients in incoherently dedispersed data. Specifically, we define a spectral modulation index in terms of the fractional variation in intensity across a spectrum. A signal whose intensity is distributed evenly across the entire band has a much lower modulation index than a spectrum with the same intensity localized in a single channel. We are interested in broadband pulses and use the modulation index to excise narrowband radio frequency interference (RFI) by applying a modulation index threshold above which candidate events are removed. The technique is tested both with simulations and using data from sources of known radio pulses (RRAT J1928+15 and giant pulses from the Crab pulsar). We find that our technique is effective at eliminating not only narrowband RFI but also spurious signals from bright, real pulses that are dedispersed at incorrect dispersion measures. The method is generalized to coherent dedispersion, image cubes, and astrophysical narrowband signals that are steady in time. We suggest that the modulation index, along with other statistics using higher-order moments, should be incorporated into signal detection pipelines to characterize and classify signals.

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.

Millisecond pulsars (MSPs) have been studied in detail since their discovery in 1982. The integrated pulse profiles of MSPs appear to be stable, which enables precision monitoring of the pulse times of arrival (TOAs). However, for individual pulses the shape and arrival phase can vary dramatically, which is known as pulse jitter. In this paper, we investigate the stability of integrated pulse profiles for 5 MSPs, and estimate the amount of jitter for PSR J0437-4715. We do not detect intrinsic profile shape variation based on integration times from ~10 to ~100 s with the provided instrumental sensitivity. For PSR J0437-4715 we calculate the jitter parameter to be f_J=0.067+-0.002, and demonstrate that the result is not significantly affected by instrumental TOA uncertainties. Jitter noise is also found to be independent of observing frequency and bandwidth around 1.4 GHz on frequency scales of <100 MHz, which supports the idea that pulses within narrow frequency scale are equally jittered. In addition, we point out that pulse jitter would limit TOA calculation for the timing observations with future telescopes like the Square Kilometre Array and the Five hundred metre Aperture Spherical Telescope. A quantitative understanding of pulse profile stability and the contribution of jitter would enable improved TOA calculations, which are essential for the ongoing endeavours in pulsar timing, such as the detection of the stochastic gravitational wave background.

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.

The relatively unexplored fast radio transient parameter space is known to be home to a variety of interesting sources, including pulsars, pulsar giant pulses and non-thermal emission from planetary magnetospheres. In addition, a variety of hypothesized but as-yet-unobserved phenomena, such as primordial black hole evaporation and prompt emission associated with coalescing massive objects have been suggested. The 2007 announcement by Lorimer et al. of the detection of a bright (30 Jy) radio pulse that was inferred to be of extragalactic origin and the subsequent consternation have demonstrated both the potential utility of bright radio pulses as probes of the interstellar medium and intergalactic medium, as well as the need for wide-field surveys characterizing the fast-transient parameter space. Here we present results from the 450 hour, 150 deg^2 Fly's Eye survey for bright dispersed radio pulses at the Allen Telescope Array (ATA). The Fly's Eye spectrometer produces 128 channel power spectra over a 209 MHz bandwidth, centered at 1430 MHz, on 44 independent signals paths originating with 30 independent ATA antennas. Data were dedispersed between 0 and 2000 pc cm^-3 and searched for pulses with dispersion measures greater than 50 pc cm^-3 between 625 us and 5 s in duration. No pulses were detected in the survey, implying a limiting rate of less than 2 sky^-1 hour^-1 for 10 millisecond duration pulses having apparent energy densities greater than 440 kJy us, or mean flux densities greater than 44 Jy. Here we present details of the instrument, experiment and observations, including a discussion of our results in light of other single pulse searches.

Here we describe our ongoing efforts to develop high-performance and sensitive instrumentation for use in the search for extra-terrestrial intelligence (SETI). These efforts include our recently deployed Search for Extraterrestrial Emissions from Nearby Developed Intelligent Populations Spectrometer (SERENDIP V.v) and two instruments currently under development; the Heterogeneous Radio SETI Spectrometer (HRSS) for SETI observations in the radio spectrum and the Optical SETI Fast Photometer (OSFP) for SETI observations in the optical band. We will discuss the basic SERENDIP V.v instrument design and initial analysis methodology, along with instrument architectures and observation strategies for OSFP and HRSS. In addition, we will demonstrate how these instruments may be built using low-cost, modular components and programmed and operated by students using common languages, e.g. ANSI C.

Timing pulses of pulsars has proved to be a most powerful technique useful to a host of research areas in astronomy and physics. Importantly, the precision of this timing is not only affected by radiometer noise, but also by intrinsic pulse shape changes, interstellar medium (ISM) evolution, instrumental distortions, etc. In this paper we review the known causes of pulse shape variations and assess their effect on the precision and accuracy of a single measurement of pulse arrival time with current instrumentation. Throughout this analysis we use the brightest and most precisely timed millisecond pulsar (MSP), PSR J0437-4715, as a case study, and develop a set of diagnostic tools to evaluate profile stability in timing observations. We conclude that most causes of distortion can be either corrected by state-of-the-art techniques or taken into account in the estimation of time-of-arrival (TOA) uncertainties. The advent of a new generation of radio telescopes (e.g. the Square Kilometre Array, SKA), and their increase in collecting area has sparked speculation about the timing precision achievable through increases in gain. Based on our analysis of current data, we predict that for normal-brightness MSPs a TOA precision of between 80 and 230 ns can be achieved at 1.4 GHz with 10-minute integrations by the SKA. The actual rms timing residuals for each pulsar will be approximately at the same level only if all the other influences on timing precision (e.g. ISM, spin noise) are either corrected, modelled, or negligible.

We assess the detectability of a nanohertz gravitational wave (GW) background with respect to additive red and white noise in the timing of millisecond pulsars. We develop detection criteria based on the cross-correlation function summed over pulsar pairs in a pulsar timing array. The distribution of correlation amplitudes is found to be non-Gaussian and highly skewed, which significantly influences detection and false-alarm probabilities. When only white noise and GWs contribute, our detection results are consistent with those found by others. Red noise, however, drastically alters the results. We discuss methods to meet the challenge of GW detection ("climbing mount significance") by distinguishing between GW-dominated and red or white-noise limited regimes. We characterize detection regimes by evaluating the number of millisecond pulsars that must be monitored in a high-cadence, 5-year timing program for a GW background spectrum $h_c(f) = A f^{-2/3}$ with $A = 10^{-15}$ yr$^{-2/3}$. Unless a sample of 20 super-stable millisecond pulsars can be found --- those with timing residuals from red-noise contributions $\sigma_r \lesssim 20$ ns --- a much larger timing program on $\gtrsim 50 - 100$ MSPs will be needed. For other values of $A$, the constraint is $\sigma_r \lesssim 20 {\rm ns} (A/10^{-15} {\rm yr}^{-2/3})$. Identification of suitable MSPs itself requires an aggressive survey campaign followed by characterization of the level of spin noise in the timing residuals of each object. The search and timing programs will likely require substantial fractions of time on new array telescopes in the southern hemisphere as well as on existing ones.

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.

We report the discovery of the 20.7 ms binary pulsar J1952+2630, made using the distributed computing project Einstein@Home in Pulsar ALFA survey observations with the Arecibo telescope. Follow-up observations with the Arecibo telescope confirm the binary nature of the system. We obtain a circular orbital solution with an orbital period of 9.4 hr, a projected orbital radius of 2.8 lt-s, and a mass function of f = 0.15 solar masses by analysis of spin period measurements. No evidence of orbital eccentricity is apparent; we set a 2-sigma upper limit e < 1.7e-3. The orbital parameters suggest a massive white dwarf companion with a minimum mass of 0.95 solar masses, assuming a pulsar mass of 1.4 solar masses. Most likely, this pulsar belongs to the rare class of intermediate mass binary pulsars. Future timing observations will aim to determine the parameters of this system further, measure relativistic effects, and elucidate the nature of the companion star.

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.

(abridged) PSR J1903+0327, a millisecond pulsar in an eccentric (e = 0.44) 95-day orbit with a (~ 1Msun) companion poses a challenge to our understanding of stellar evolution in binary and multiple-star systems. Here we describe optical and radio observations which rule out most of the scenarios proposed to explain formation of this system. Radio timing measurements of three post-Keplerian effects yield the most precise measurement of the mass of a millisecond pulsar to date: 1.667 +/- 0.021 solar masses (99.7% confidence limit) (...). Optical spectroscopy of a proposed main sequence counterpart show that its orbital motion mirrors the pulsar's 95-day orbit; being therefore its binary companion (...) The optical detection also provides a measurement of the systemic radial velocity of the binary; this and the proper motion measured from pulsar timing allow the determination of the systemic 3-D velocity in the Galaxy. We find that the system is always within 270 pc of the plane of the Galaxy, but always more than 3 kpc away from the Galactic centre. Thus an exchange interaction in a dense stellar environment (like a globular cluster or the Galactic centre) is not likely to be the origin of this system. We suggest that after the supernova that formed it, the neutron star was in a tight orbit with a main-sequence star, the present companion was a tertiary farther out. The neutron star then accreted matter from its evolving inner MS companion, forming a millisecond pulsar. The former donor star then disappears, either due to a chaotic 3-body interaction with the outer star (caused by the expansion of the inner orbit that necessarily results from mass transfer), or in the case of a very compact inner system, due to ablation/accretion by the newly formed millisecond pulsar.

We investigate rotational spin noise (referred to as timing noise) in non-accreting pulsars: millisecond pulsars, canonical pulsars, and magnetars. Particular attention is placed on quantifying the strength and non-stationarity of timing noise in millisecond pulsars because the long-term stability of these objects is required to detect nanohertz gravitational radiation. We show that a single scaling law is sufficient to characterize timing noise in millisecond and canonical pulsars while the same scaling law underestimates the levels of timing noise in magnetars. The scaling law, along with a detailed study of the millisecond pulsar B1937+21, leads us to conclude that timing noise is latent in most millisecond pulsars and will be measurable in many objects when better arrival time estimates are obtained over long data spans. The sensitivity of a pulsar timing array to gravitational radiation is strongly affected by any timing noise. We conclude that detection of proposed gravitational wave backgrounds will require the analysis of more objects than previously suggested over data spans that depend on the spectra of both the gravitational wave background and of the timing noise. It is imperative to find additional millisecond pulsars in current and future surveys in order to reduce the effects of timing noise.

This paper describes a comprehensive measurement model for the error budget of pulse arrival times with emphasis on intrinsic pulse jitterand plasma propagation effects (particularly interstellar scattering), which are stochastic in time and have diverse dependences on radio frequency. To reduce their contribution, timing measurements can be made over a range of frequencies that depends on a variety of pulsar and instrumentation-dependent factors that we identify. A salient trend for high signal-to-noise measurements of millisecond pulsars is that time-of-arrival precision is limited either by irreducible interstellar scattering or by pulse-phase jitter caused by variable emission within pulsar magnetospheres. A cap on timing errors implies that pulsars must be confined to low dispersion measures (DMs) and observed at high frequencies. Use of wider bandwidths that increase signal-to-noise ratios will degrade timing precision if nondispersive chromatic effects are not mitigated. The allowable region in the DM-frequency plane depends on how chromatic timing perturbations are addressed. Without mitigation, observations at 1.4~GHz or 5~GHz are restricted to $\DM\lesssim 30$ and $\lesssim 100~\DMu$, respectively. With aggressive mitigation of interstellar scattering and use of large telescopes to provide adequate sensitivity at high frequencies (e.g. Arecibo, FAST, phase 1 of the SKA, and the SKA), pulsars with DMs up to 500~$\DMu$ can be used in precision timing applications. We analyze methods that fit arrival times vs. frequency at a given epoch prior to multi-epoch fitting. While the terms of greatest astrophysical interest are achromatic (e.g. orbital and gravitational wave perturbations), measurements may ultimately be limited by similarly achromatic stochasticity in a pulsar's spin rate.

The SKA at mid and low frequencies will be constructed in two distinct phases, the first being a subset of the second. This document defines the main scientific goals and baseline technical concept for the SKA Phase 1 (SKA_1). The major science goals for SKA_1 will be to study the history and role of neutral Hydrogen in the Universe from the dark ages to the present-day, and to employ pulsars as probes of fundamental physics. The baseline technical concept of SKA_1 will include a sparse aperture array operating at frequencies up to 450 MHz, and an array of dishes, initially operating at frequencies up to 3 GHz but capable of 10 GHz in terms of antenna surface accuracy. An associated Advanced Instrumentation Program (AIP) allows further development of new technologies currently under investigation. Construction will take place in 2016-2019 at a total capital cost of 350M\texteuro, including an element for contingency. The cost estimates of the SKA_1 telescope are now the subject of a more detailed and thorough costing exercise led by the SKA Project Development Office (SPDO). The 350 M\texteuro total for SKA_1 is a cost-constrained cap; an additional contingency is to reduce the overall scope of the project. The design of the SKA_1 is expected to evolve as the major cost estimates are refined, in particular the infrastructure costs at the two sites. The SKA_1 facility will represent a major step forward in terms of sensitivity, survey speed, image fidelity, temporal resolution and field-of-view. It will open up new areas of discovery space and demonstrate the science and technology underpinning the SKA Phase 2 (SKA_2).

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 192 countries to "mine" large data sets. It has now found a 40.8 Hz isolated pulsar in radio survey data from the Arecibo Observatory taken in February 2007. Additional timing observations indicate that this pulsar is likely a disrupted recycled pulsar. PSR J2007+2722's pulse profile is remarkably wide with emission over almost the entire spin period; the pulsar likely has closely aligned magnetic and spin axes. The massive computing power provided by volunteers should enable many more such discoveries.

In order to study precession and interstellar magnetic field variations, we measured the polarized position angle of 81 pulsars at several-month intervals for four years. We show that the uncertainties in a single-epoch measurement of position angle is usually dominated by random pulse-to-pulse jitter of the polarized subpulses. Even with these uncertainties, we find that the position angle variations in 19 pulsars are significantly better fitted (at the 3 {\sigma} level) by a sinusoid than by a constant. Such variations could be caused by precession, which would then indicate periods of ~ (200 - 1300) d and amplitudes of ~(1 - 12) degrees. We narrow this collection to four pulsars that show the most convincing evidence of sinusoidal variation in position angle. Also, in a handful of pulsars, single discrepant position angle measurements are observed which may result from the line of sight passing across a discrete ionized, magnetized structure. We calculate the standard deviation of position angle measurements from the mean for each pulsar, and relate these to limits on precession and interstellar magnetic field variations.

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.

The discovery of cosmic radio emission by Karl Jansky in the course of searching for the source of interference to telephone communications and the instrumental advances which followed, have led to a series of new paradigm changing astronomical discoveries. These discoveries, which to a large extent define much of modern astrophysical research were the result of the right people being in the right place at the right time using powerful new instruments, which in many cases they had designed and built. They were not the result of trying to test any particular theoretical model or trying to answer previously posed questions, but they opened up whole new areas of exploration and discovery. Rather many important discoveries came from military or communications research; others while looking for something else; and yet others from just looking. Traditionally, the designers of big telescopes invariably did not predict what the telescopes would ultimately be known for. The place in history of the next generation of telescopes will not likely be found in the science case created to justify their construction, but in the unexpected new phenomena, new theories, and new ideas which will emerge from these discoveries. It is important that those who are in a position to filter research proposals and plans not dismiss as butterfly collecting, investigations which explore new areas without having predefined the result they are looking for. Progress must also allow for new discoveries, as well as for the explanation of old discoveries. New telescopes need to be designed with the flexibility to make new discoveries which will invariably raise new questions and new problems.

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) is a consortium of astronomers whose goal is the creation of a galactic scale gravitational wave observatory sensitive to gravitational waves in the nHz-microHz band. It is just one component of an international collaboration involving similar organizations of European and Australian astronomers who share the same goal. Gravitational waves, a prediction of Einstein's general theory of relativity, are a phenomenon of dynamical space-time generated by the bulk motion of matter, and the dynamics of space-time itself. They are detectable by the small disturbance they cause in the light travel time between some light source and an observer. NANOGrav exploits radio pulsars as both the light (radio) source and the clock against which the light travel time is measured. In an array of radio pulsars gravitational waves manifest themselves as correlated disturbances in the pulse arrival times. The timing precision of today's best measured pulsars is less than 100 ns. With improved instrumentation and signal-to-noise it is widely believed that the next decade could see a pulsar timing network of 100 pulsars each with better than 100 ns timing precision. Such a pulsar timing array (PTA), observed with a regular cadence of days to weeks, would be capable of observing supermassive black hole binaries following galactic mergers, relic radiation from early universe phenomena such as cosmic strings, cosmic superstrings, or inflation, and more generally providing a vantage on the universe whose revolutionary potential has not been seen in the 400 years since Galileo first turned a telescope to the heavens.

We present radio transient search algorithms, results, and statistics from the ongoing Arecibo Pulsar ALFA (PALFA) Survey of the Galactic plane. We have discovered seven objects by detecting isolated dispersed pulses and one of the new discoveries has a duty cycle of 0.01%, the smallest known. The impact of selection effects on the detectability and classification of intermittent sources is discussed, and the relative efficiencies of periodicity vs. single pulse searches are compared for various pulsar classes. We find that scintillation, off-axis detection and few rotation periods within an observation may misrepresent normal periodic pulsars as intermittent sources. Finally, we derive constraints on transient pulse rate and flux density from the PALFA survey parameters and results.

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.

The time domain of the sky has been only sparsely explored. Nevertheless, recent discoveries from limited surveys and serendipitous discoveries indicate that there is much to be found on timescales from nanoseconds to years and at wavelengths from meters to millimeters. These observations have revealed unexpected phenomena such as rotating radio transients and coherent pulses from brown dwarfs. Additionally, archival studies have found not-yet identified radio transients without optical or high-energy hosts. In addition to the known classes of radio transients, possible other classes of objects include extrapolations from known classes and exotica such as orphan gamma-ray burst afterglows, radio supernovae, tidally-disrupted stars, flare stars, magnetars, and transmissions from extraterrestrial civilizations. Over the next decade, meter- and centimeter-wave radio telescopes with improved sensitivity, wider fields of view, and flexible digital signal processing will be able to explore radio transient parameter space more comprehensively and systematically.

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.

We highlight recent theoretical and observational progress in several areas of neutron star astrophysics, and discuss the prospect for advances in the next decade.

The characteristic physical timescales near stellar-mass compact objects are measured in milliseconds. These timescales -- the free-fall time, the fastest stable orbital period, and stellar spin periods -- encode the fundamental physical properties of compact objects: mass, radius, and angular momentum. The characteristic temperature of matter in the vicinity of neutron stars is such that the principal electromagnetic window into their realms is the X-ray band. Because of these connections to the fundamental properties of neutron stars, X-ray timing studies remain today the most direct means of probing their structure and dynamics. While current X-ray observatories have revealed many relevant and fascinating phenomena, they lack the sensitivity to fully exploit them to uncover the fundamental properties of compact objects and their extreme physics. With this white paper, we summarize and highlight the science opportunities that will accompany an order-of-magnitude improvement in X-ray timing sensitivity, a goal attainable in the coming decade.

More than four decades after the discovery of pulsars, the composition of matter at their cores is still a mystery. This white paper summarizes how recent high-precision measurements of millisecond pulsar masses have introduced new experimental constraints on the properties of super-dense matter, and how continued timing of intriguing new objects, coupled with radio telescope surveys to discover more pulsars, might introduce significantly more stringent constraints.

The Square Kilometre Array (SKA) is a planned multi purpose radio telescope with a collecting area approaching 1 million square metres. One of the key science objectives of the SKA is to provide exquisite strong-field tests of gravitational physics by finding and timing pulsars in extreme binary systems such as a pulsar-black hole binary. To find out how three preliminary SKA configurations will affect a pulsar survey, we have simulated SKA pulsar surveys for each configuration. We estimate that the total number of normal pulsars the SKA will detect, using only the 1-km core and 30 minutes integration time, is around 14000 normal pulsar and 6000 millisecond pulsars. We describe a simple strategy for follow-up timing observations and find that, depending on the configuration, it would take 1-6 days to obtain a single timing point for 14000 pulsars. Obtaining a single timing point for the high-precision timing projects of the SKA, will take less than 14 hours, 2 days, or 3 days, depending on the configuration. The presence of aperture arrays will be of great benefit here. We also study the computational requirements for beam forming and data analysis for a pulsar survey. Beam forming of the full field of view of the single-pixel feed 15-m dishes using the 1-km core of the SKA requires about 2.2*10^15 operations per second. The corresponding data rate from such a pulsar survey is about 4.7*10^11 bytes per second. The required computational power for a deep real time analysis is estimated to be 1.2*10^16 operations per second. For an aperture array or dishes equipped with phased array feeds, the survey can be performed faster, but the computational requirements and data rates will go up.

We present the discovery of the Vela-like radio pulsar J1856+0245 in the Arecibo PALFA survey. PSR J1856+0245 has a spin period of 81ms, a characteristic age of 21kyr, and a spin-down luminosity Edot = 4.6 x 10^36 ergs/s. It is positionally coincident with the TeV gamma-ray source HESS J1857+026, which has no other known counterparts. Young, energetic pulsars create wind nebulae, and more than a dozen pulsar wind nebulae have been associated with very-high-energy (100GeV-100TeV) gamma-ray sources discovered with the HESS telescope. The gamma-ray emission seen from HESS J1857+026 is potentially produced by a pulsar wind nebula powered by PSR J1856+0245; faint X-ray emission detected by ASCA at the pulsar's position supports this hypothesis. The inferred gamma-ray efficiency is epsilon_gamma = L_gamma/Edot = 3.1% (1-10TeV, for a distance of 9kpc), comparable to that observed in similar associations.

Binary pulsar systems are superb probes of stellar and binary evolution and the physics of extreme environments. In a survey with the Arecibo telescope, we have found PSR J1903+0327, a radio pulsar with a rotational period of 2.15 ms in a highly eccentric (e = 0.44) 95-day orbit around a solar mass companion. Infrared observations identify a possible main-sequence companion star. Conventional binary stellar evolution models predict neither large orbital eccentricities nor main-sequence companions around millisecond pulsars. Alternative formation scenarios involve recycling a neutron star in a globular cluster then ejecting it into the Galactic disk or membership in a hierarchical triple system. A relativistic analysis of timing observations of the pulsar finds its mass to be 1.74+/-0.04 Msun, an unusually high value.

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.

In this second of a two-paper sequence, we present Monte Carlo population simulation results of brown dwarf companion data collected during the Cornell High-order Adaptive Optics Survey for brown dwarf companions (CHAOS). Making reasonable assumptions of orbital parameters (random inclination, random eccentricity and random longitude of pericentre) and age distributions, and using published mass functions, we find that the brown dwarf companion fraction around main sequence stars is 0.0%-9.3% for the 25-100 AU semi-major axis region. We find a corresponding L-dwarf companion fraction of 0.0%-3.3%. We compare our population analysis methods and results with techniques and results presented by several other groups. In this comparison we discover that systematic errors (most notably resulting from orbital projection effects) occur in the majority of previously published brown dwarf companion population estimates, leading authors to claim results not supported by the observational data.

We propose that neutral, circumpulsar debris entering the light cylinder can account for many time-dependent pulsar phenomena that are otherwise difficult to explain. Neutral material avoids propeller ejection and injects sufficient charges -- after heating, evaporation, and ionization -- to alter current flows and pair-production and thus trigger, detune, or extinguish coherent emission. Relevant phenomena, with time scales from seconds to months, include nulls, rotating radio transients (RRATs), rapid changes in pulse profile (``mode changes''), variable subpulse drift rates, quasi-periodic bursts from B1931+24, and torque variations. Over the 10 Myr lifetime of a canonical pulsar with trillion-gauss surface magnetic field, less than a millionth of an Earth mass of material is needed to modulate the Goldreich-Julian current by 100%. Circumpulsar material originates from metal-rich, supernova fallback gas that aggregates into asteroids. Debris disks can inject sufficient material on time scales of interest, yet be too tenuous to form large planets detectable in pulse timing data. Asteroid migration results from collisions and the radiation-driven Yarkovsky and Poynting-Robertson effects. For B1931+24, an asteroid in a $\sim 40$~day elliptical orbit pollutes the magnetosphere stochastically through collisions with other debris. Injection is less likely for hot, young and highly magnetized pulsars or millisecond pulsars that pre-ionize any debris material well outside their small magnetospheres. Injection effects will therefore be most prominent in long-period, cooler pulsars, consistent with the distribution of relevant objects in perid and period derivative. A pulsar's spin history and its radiation-beam orientation may influence whether it displays nulling, RRATs and other effects.

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.

The `radio sky' is relatively unexplored for transient signals, although the potential of radio-transient searches is high, as demonstrated recently by the discovery of a previously unknown type of source which varies on timescales of minutes to hours. Here we report a new large-scale search for radio sources varying on much shorter timescales. This has revealed 11 objects characterized by single, dispersed bursts having durations between 2 and 30 ms. The average time intervals between bursts range from 4 minutes to 3 hours, with radio emission typically detectable for < 1 s per day. From an analysis of the burst arrival times, we have identified periodicities in the range 0.4 - 7 s for ten of the 11 sources, suggesting a rotating neutron star origin. Despite the small number of sources presently detected, their ephemeral nature implies a total Galactic population which significantly exceeds that of the regularly pulsing radio pulsars. Five of the ten sources have periods greater than 4 s, and period derivatives have been measured for three of the sources, with one having a very high inferred magnetic field of 5e13 G, suggesting that this new population is related to other classes of isolated neutron stars observed at X-ray and gamma-ray wavelengths.

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.

We report results from the initial stage of a long-term pulsar survey of the Galactic plane using the Arecibo L-band Feed Array (ALFA), a seven-beam receiver operating at 1.4 GHz with 0.3 GHz bandwidth. The search targets Galactic latitudes |b| < 5 deg in the longitude ranges 32 deg < l < 77 deg and 168 deg < l < 77 deg. Data discussed here were collected over a 100 MHz passband centered on 1.42 GHz using a spectrometer that recorded 256 channels every 64 microsec. In a preliminary, standard period-DM analysis, we have detected 29 previously known pulsars and discovered 11 new ones. One of these, with a period of 69 ms and a low characteristic age of 82 kyr, is a plausible candidate for association with the unidentified EGRET source 3EG J1928+1733. Another is a non-recycled pulsar in a relativistic binary with orbital period of 3.98 hr. We also search the data for isolated dispersed pulses, a technique that yielded discovery of an extremely sporadic radio emitter with a spin period of 1.2 s. Simulations we have carried out indicate that about 1000 new pulsars will be found in the ALFA survey. In addition to providing a large sample for use in population analyses and for probing the magnetoionic interstellar medium, the survey maximizes the chances of finding rapidly spinning millisecond pulsars and pulsars in compact binary systems. Our search algorithms will exploit the multiple data streams from ALFA to discriminate between radio frequency interference and celestial signals, including pulsars and possibly new classes of transient radio sources.

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.

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.

We report on timing and single-pulse observations of 17 pulsars discovered at the Arecibo observatory. The highlights of our sample are the recycled pulsars J1829+2456, J1944+0907 and the drifting subpulses observed in PSR J0815+0939. For the double neutron star binary J1829+2456, in addition to improving upon our existing measurement of relativistic periastron advance, we have now measured the pulsar's spin period derivative. This new result sets an upper limit on the transverse speed of 120 km/s and a lower limit on the characteristic age of 12.4 Gyr. From our measurement of proper motion of the isolated 5.2-ms pulsar J1944+0907, we infer a transverse speed of 188 +/- 65 km/s. This is higher than that of any other isolated millisecond pulsar. An estimate of the speed, using interstellar scintillation, of 235 +/- 45 km/s indicates that the scattering medium along the line of sight is non-uniform. We discuss the drifting subpulses detected from three pulsars in the sample, in particular the remarkable drifting subpulse properties of the 645-ms pulsar J0815+0939. Drifting is observed in all four components of the pulse profile, with the sense of drift varying among the different components. This unusual `bi-drifting'' behaviour challenges standard explanations of the drifting subpulse phenomenon.

PSR J0751+1807 is a millisecond pulsar in a circular 6 hr binary system with a helium white dwarf secondary. Through high precision pulse timing measurements with the Arecibo and Effelsberg radio telescopes, we have detected the decay of its orbit due to emission of gravitational radiation. This is the first detection of the relativistic orbital decay of a low-mass, circular binary pulsar system. The measured rate of change in orbital period, corrected for acceleration biases, is dP_b/dt=(-6.4+-0.9)x10^-14. Interpreted in the context of general relativity, and combined with measurement of Shapiro delay, it implies a pulsar mass of 2.1+-0.2 solar masses, the most massive pulsar measured. This adds to the emerging trend toward relatively high neutron star masses in neutron star--white dwarf binaries. Additionally, there is some evidence for an inverse correlation between pulsar mass and orbital period in these systems. We consider alternatives to the general relativistic analysis of the data, and we use the pulsar timing data to place limits on violations of the strong equivalence principle.

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.

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.

Our theory relates the secondary spectrum, the 2D power spectrum of the radio dynamic spectrum, to the scattered pulsar image in a thin scattering screen geometry. Recently discovered parabolic arcs in secondary spectra are generic features for media that scatter radiation at angles much larger than the rms scattering angle. Each point in the secondary spectrum maps particular values of differential arrival-time delay and fringe rate (or differential Doppler frequency) between pairs of components in the scattered image. Arcs correspond to a parabolic relation between these quantities through their common dependence on the angle of arrival of scattered components. Arcs appear even without consideration of the dispersive nature of the plasma. Arcs are more prominent in media with negligible inner scale and with shallow wavenumber spectra, such as the Kolmogorov spectrum, and when the scattered image is elongated along the velocity direction. The arc phenomenon can be used, therefore, to constrain the inner scale and the anisotropy of scattering irregularities for directions to nearby pulsars. Arcs are truncated by finite source size and thus provide sub micro arc sec resolution for probing emission regions in pulsars and compact active galactic nuclei. Multiple arcs sometimes seen signify two or more discrete scattering screens along the propagation path, and small arclets oriented oppositely to the main arc persisting for long durations indicate the occurrence of long-term multiple images from the scattering screen.

The radio pulsar B1951+32 and the supernova remnant CTB 80 provide a rich laboratory for the study of neutron stars and supernova remnants. Here, we present ground-based optical and near-infrared observations of them, along with X-ray observations with Chandra and a re-analysis of archival data obtained with the Hubble Space Telescope. The X-ray observations reveal a cometary pulsar wind nebula which appears to be confined by a bow shock produced by high-velocity motion of the pulsar, making PSR B1951+32 a rare pulsar exhibiting both an H alpha bow shock and a shocked X-ray pulsar wind nebula. The distribution of H alpha and radio continuum emission is indicative of a contact discontinuity of the shocked pulsar winds and shocked ambient medium at \~0.05 pc. On the other hand, the optical synchrotron knot of PSR B1951+32 likely has a flat spectrum in the optical and near-infrared wavebands, and our astrometry is consistent with only one of the two reported optical counterpart candidates for the pulsar.

We discuss four recycled pulsars found in Arecibo drift-scan searches. PSR J1944+0907 has a spin period of 5.2 ms and is isolated. The 5.8-ms pulsar J1453+19 may have a low-mass companion. We discuss these pulsars in the context of isolated millisecond pulsar formation and the minimum spin period of neutron stars. The isolated 56-ms pulsar J0609+2130 is possibly the remnant of a disrupted double neutron star binary. The 41-ms pulsar J1829+2456 is in a relativistic orbit. Its companion is most likely another neutron star, making this the eighth known double neutron star binary system.

High precision astrometric data have enabled us to determine the trajectories through the Galactic potential for a growing number of pulsars. This has resulted in the discovery of a pulsar pair (B2020+28 and B2021+51) that has a common origin in the Cygnus Superbubble or in one of its related OB associations at an epoch which is comparable with the spin-down ages of the pulsars. Analysis of the Galactic orbits indicates that the progenitors of the pulsars had similar masses and were in a binary system, which was disrupted after the second supernova explosion. The implied pulsar birth velocities are consistent with the high velocities of neutron stars in general. The initial spin period of the pulsar that was formed in the second supernova explosion was approximately 200 ms. A further increase in astrometric accuracy will allow us to more tightly constrain the birth velocities and the kick velocities that were imparted by the two respective supernova explosions. Two additional pulsars in the sample of 24 with parallax measurements may also have originated in the Cygnus Superbubble.

We report the discovery of a new binary pulsar, PSR J1829+2456, found during a mid-latitude drift-scan survey with the Arecibo telescope. Our initial timing observations show the 41-ms pulsar to be in a 28-hr, slightly eccentric, binary orbit. The advance of periastron, omegadot = 0.28 +/- 0.01 deg/yr is derived from our timing observations spanning 200 days. Assuming that the advance of periastron is purely relativistic and a reasonable range of neutron star masses for PSR J1829+2456 we constrain the companion mass to be between 1.22 Msun and 1.38 Msun, making it likely to be another neutron star. We also place a firm upper limit on the pulsar mass of 1.38 Msun. The expected coalescence time due to gravitational-wave emission is long (~60 Gyr) and this system will not significantly impact upon calculations of merger rates that are relevant to upcoming instruments such as LIGO.

We have made observations of 98 low-Galactic-latitude pulsars to measure pulse broadening caused by multipath propagation through the interstellar medium. Data were collected with the 305-m Arecibo telescope at four radio frequencies between 430 and 2380 MHz. We used a CLEAN-based algorithm to deconvolve interstellar pulse broadening from the measured pulse shapes. We employed two distinct pulse broadening functions (PBFs): PBF$_1$ is appropriate for a thin screen of scattering material between the Earth and a pulsar, while PBF$_2$ is appropriate for scattering material uniformly distributed along the line of sight from the Earth to a pulsar. We found that some observations were better fit by PBF$_1$ and some by PBF$_2$. Pulse broadening times ($\tau_d$) are derived from fits of PBFs to the data, and are compared with the predictions of a smoothed model of the Galactic electron distribution. Several lines of sight show excess broadening, which we model as clumps of high density scattering material. A global analysis of all available data finds that the pulse broadening scales with frequency, $\nu$, as $\taud \propto\nu^{-\alpha}$ where $\alpha\sim 3.9\pm 0.2$. This is somewhat shallower than the value $\alpha=4.4$ expected from a Kolmogorov medium, but could arise if the spectrum of turbulence has an inner cutoff at $\sim $300--800 km. A few objects follow particularly shallow scaling laws (the mean scaling index $\meanalpha \sim 3.1 \pm 0.1$ and $ \sim 3.8 \pm 0.2$ respectively for the case of PBF$_1$ and PBF$_2$), which may arise from large scale refraction or from the truncation of scattering screens transverse to the Earth--pulsar line of sight.

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.

The Guitar nebula is a spectacular example of an H-alpha bow shock nebula produced by the interaction of a neutron star with its environment. The radio pulsar B2224+65 is traveling at ~800--1600 km/s (for a distance of 1--2 kpc), placing it on the high-velocity tail of the pulsar velocity distribution. Here we report time evolution in the shape of the Guitar nebula, the first such observations for a bow shock nebula, as seen in H-alpha imaging with the Hubble Space Telescope. The morphology of the nebula provides no evidence for anisotropy in the pulsar wind, nor for fluctuations in the pulsar wind luminosity. The nebula shows morphological changes over two epochs spaced by seven years that imply the existence of significant gradients and inhomogeneities in the ambient interstellar medium. These observations offer astrophysically unique, in situ probes of length scales between 5E-4 pc and 0.012 pc. Model fitting suggests that the nebula axis -- and thus the three-dimensional velocity vector -- lies within 20 degrees of the plane of the sky, and also jointly constrains the distance to the neutron star and the ambient density.

We report the discovery and initial timing observations of a 55.7-ms pulsar, J0609+2130, found during a 430-MHz drift-scan survey with the Arecibo radio telescope. With a spin-down rate of $3.1 \times 10^{-19}$ s s$^{-1}$ and an inferred surface dipole magnetic field of only $4.2 \times 10^{9}$ G, J0609+2130 has very similar spin parameters to the isolated pulsar J2235+1506 found by Camilo, Nice & Taylor (1993). While the origin of these weakly magnetized isolated neutron stars is not fully understood, one intriguing possibility is that they are the remains of high-mass X-ray binary systems which were disrupted by the supernova explosion of the secondary star.

Using a likelihood analysis (updated since McLaughlin & Cordes 2000) and EGRET detections, upper limits and diffuse background measurements, we find a best-fit luminosity law for the gamma-ray pulsar population. We find that roughly 30 of the 170 unidentified EGRET sources are likely to be pulsars. This is roughly twice the number of known radio pulsars which are plausibly associated with unidentified EGRET sources. We predict that AGILE will detect roughly 70 pulsars as point sources, including 12 which will be able to be detected in blind periodicity searches. GLAST should detect roughly 1200 pulsars (including only 200 currently known radio pulsars), 210 of which will be able to be detected in blind searches. We discuss methods of searching for pulsars in gamma-ray data and present results from our searches for gamma-ray periodicities from new radio pulsars associated with unidentified EGRET sources.

We report on new pulsars discovered in Arecibo drift-scan data. Processing of 2200 square degrees of data has resulted in the detection of 41 known and 12 new pulsars. New pulsars include two millisecond pulsars, one solitary and one binary recycled pulsar, and one pulsar with very unusual pulse profile morphology and complex drifting subpulse behavior.

We have measured Faraday Rotation Measures (RMs) at Arecibo Observatory for 36 pulsars, 17 of them new. We combine these and earlier measurements to study the galactic magnetic field and its possible temporal variations. Many RM values have changed significantly on several-year timescales, but these variations probably do not reflect interstellar magnetic field changes. By studying the distribution of pulsar RMs near the plane in conjunction with the new NE2001 electron density model, we note the following structures in the first galactic longitude quadrant: (1) The local field reversal can be traced as a null in RM in a 0.5-kpc wide strip interior to the Solar Circle, extending \~7 kpc around the Galaxy. (2) Steadily increasing RMs in a 1-kpc wide strip interior to the local field reversal, and also in the wedge bounded by 42<l<52 deg, indicate that the large-scale field is approximately steady from the local reversal in to the Sagittarius arm. (3) The RMs in the 1-kpc wide strip interior to the Sagittarius arm indicate another field reversal in this strip. (4) The RMs in a final 1-kpc wide interior strip, straddling the Scutum arm, also support a second field reversal interior to the Sun,between the Sagittarius and Scutum arms. (5) Exterior to the nearby reversal, RMs from 60<l<78 deg show evidence for two reversals, on the near and far side of the Perseus arm. (6) In general, the maxima in the large-scale fields tend to lie along the spiral arms, while the field minima tend to be found between them. We have also determined polarized profiles of 48 pulsars at 430 MHz. We present morphological pulse profile classifications of the pulsars, based on our new measurements and previously published data.

Pulsars have mean space velocities >~500 km/s. The consequent ram pressure results in tight confinement of the star's energetic wind, driving a bow shock into the surrounding medium. Pulsar bow shocks have long been regarded as a curiosity, but new optical and X-ray observations are both rapidly expanding the sample of such sources, and are offering new ways to probe the interaction between pulsars and their environments. Here we discuss some of these new results, and explain how these data can be used to probe the density and structure of neutral gas in the interstellar medium.

We discuss optimal detection of fast radio transients from astrophysical objects while taking into account the effects of propagation through intervening ionized media, including dispersion, scattering and scintillation.Our analysis applies to the giant-pulse phenomenon exhibited by some pulsars, for which we show examples, and to radio pulses from other astrophysical sources, such as prompt radio emission from gamma-ray burst sources and modulated signals from extra-terrestrial civilizations.

We discuss the giant-pulse phenomenon exhibited by pulsars and the distances to which giant pulses might be detected from extragalactic pulsars. We describe the conditions under which a single-pulse search is more sensitive than a standard periodicity search. We find that, for certain pulse-amplitude distribution power laws and time series lengths, single-pulse searches can be superior. We present the results of searches toward several extragalactic targets, including M33, the LMC (PSR B0540-69) and several other galaxies. While we have not conclusively detected giant pulses from any of these targets, these searches illustrate the methodology of, issues related to and difficulties in these types of searches.

We analyze the Crab pulsar at ten frequencies from 0.43 to 8.8 GHz using data obtained at the Arecibo Observatory. Giant pulses occur only in the main and interpulse components manifest from radio frequencies to gamma-ray energies. Individual giant pulses reach brightness temperatures of at least $10^{32}$K in our data, which do not resolve the narrowest pulses, and are known to reach $10^{37}$K in nanosecond-resolution observations (Hankins et al 2003). Giant pulses are therefore the brightest known in the observable universe and represent an important milestone for theories of the pulsar emission mechanism to explain. Their short durations allow them to serve as especially sensitive probes of the Crab Nebula and the interstellar medium. We analyze frequency structure in individual giant pulses using a scintillating, amplitude-modulated,polarized shot-noise model. The frequency structure associated with multipath propagation decorrelates on a time scale of 25 sec at 1.5 GHz, which requires that multipath propagation be strongly influenced by material within the Crab Nebula. Additional frequency structure decorrelates faster than one spin period, as would be expected from the shot-noise pattern of nanosecond duration pulses emitted by the pulsar. Taking into account the Crab pulsar's locality inside a bright supernova remnant, we conclude that the brightest pulse in a typical 1-hour observation would be most easily detectable in our lowest frequency band (0.43 GHz) to a distance of 1.6 Mpc. We also discuss the detection of such pulses using future instruments such as LOFAR and the SKA.

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.

We report new pulsars discovered in drift-scan data taken by two collaborations (Berkeley/Cornell and STScI/NAIC) during the latter stages of the Arecibo upgrade period. The data were taken with the Penn State Pulsar Machine and are being processed on the COBRA cluster at Jodrell Bank. Processing is roughly 70% complete and has resulted in the detection of 10 new and 31 known pulsars, in addition to a number of pulsar candidates. The 10 new pulsars include one pulsar with a spin-period of 55 ms and another with a spin period of 5.8 ms. At the completion of the processing, we expect to have discovered roughly 20 new pulsars. All new pulsars are being subjected to a program of followup observations at Arecibo to determine spin and astrometric parameters.

Multipath propagation in the interstellar medium distorts radio pulses, an effect predominant for distant pulsars observed at low frequencies. Typically, broadened pulses are analyzed to determine the amount of propagation-induced pulse broadening, but with little interest in determining the undistorted pulse shapes. In this paper we develop and apply a method that recovers both the intrinsic pulse shape and the pulse broadening function that describes the scattering of an impulse. The method resembles the CLEAN algorithm used in synthesis imaging applications, although we search for the best pulse broadening function, and perform a true deconvolution to recover intrinsic pulse structre. As figures of merit to optimize the deconvolution, we use the positivity and symmetry of the deconvolved result along with the mean square residual and the number of points below a given threshold. Our method makes no prior assumptions about the intrinsic pulse shape and can be used for a range of scattering functions for the interstellar medium. It can therefore be applied to a wider variety of measured pulse shapes and degrees of scattering than the previous approaches. We apply the technique to both simulated data and data from Arecibo observations.

We describe a new and efficient technique which we call sideband or phase-modulation searching that allows one to detect short period binary pulsars in observations longer than the orbital period. The orbital motion of the pulsar during long observations effectively modulates the phase of the pulsar signal causing sidebands to appear around the pulsar spin frequency and its harmonics in the Fourier transform. For the majority of binary radio pulsars or Low-Mass X-ray Binaries (LMXBs), large numbers of sidebands are present, allowing efficient searches using Fourier transforms of short portions of the original power spectrum. Analysis of the complex amplitudes and phases of the sidebands can provide enough information to solve for the Keplerian orbital parameters. This technique is particularly applicable to radio pulsar searches in globular clusters and searches for coherent X-ray pulsations from LMXBs and is complementary to more standard ``acceleration'' searches.

The interaction of high-velocity neutron stars with the interstellar medium produces bow shock nebulae, where the relativistic neutron star wind is confined by ram pressure. We present multi-wavelength observations of the Guitar Nebula, including narrow-band H-alpha imaging with HST/WFPC2, which resolves the head of the bow shock. The HST observations are used to fit for the inclination of the pulsar velocity vector to the line of sight, and to determine the combination of spindown energy loss, velocity, and ambient density that sets the scale of the bow shock. We find that the velocity vector is most likely in the plane of the sky. We use the Guitar Nebula and other observed neutron star bow shocks to test scaling laws for their size and H-alpha emission, discuss their prevalence, and present criteria for their detectability in targeted searches. The set of H-alpha bow shocks shows remarkable consistency, in spite of the expected variation in ambient densities and orientations. Together, they support the assumption that a pulsar's spindown energy losses are carried away by a relativistic wind that is indistinguishable from being isotropic. Comparison of H-alpha bow shocks with X-ray and nonthermal, radio-synchrotron bow shocks produced by neutron stars indicates that the overall shape and scaling is consistent with the same physics. It also appears that nonthermal radio emission and H-alpha emission are mutually exclusive in the known objects and perhaps in all objects.

The pulsar PSR J0631+1036 was discovered during a radio search of Einstein X-ray source error circles. A detection of a 288ms sinusoidal modulation in the ASCA lightcurve, the same period as the radio pulsar, appeared to confirm the association of the X-ray source and the pulsar. Its X-ray spectrum was said to be similar to that of middle aged gamma-ray pulsars such as Geminga. However, an XMM-Newton observation of the PSR J0631+1035 field, along with a re-analysis of VLA data confirming the timing position of the pulsar, show a 75'' discrepancy in location of the X-ray source and the pulsar, and therefore these cannot be the same object. The X-ray source appears to be the counterpart of an A0 star, detected by the XMM-Newton Optical Monitor. No 288ms period was detected from either the area around the pulsar or the bright X-ray source. The upper limit on the X-ray luminosity with relation to the empirically observed correlation between radio measured dE/dt and X-ray luminosity is discussed.

(Abridged) We infer the velocity distribution of radio pulsars by modelling their birth, evolution, and detection in large-scale 0.4 GHz pulsar surveys, and by comparing model distributions of measurable pulsar properties with survey data using a likelihood function. We test models that characterize a population's birth rate, luminosity, shutoff of radio emission, birth locations, and birth velocities. We infer that the radio beam luminosity (i) is comparable to the energy flux of relativistic particles in models for spin-driven magnetospheres, signifying that radio emission losses reach nearly 100% for the oldest pulsars; and (ii) scales approximately as sqrt(Edot) which, in magnetosphere models, is proportional to the voltage drop available for acceleration of particles. We find that a two-component velocity distribution with characteristic velocities of 90 km/s and 500 km/s is greatly preferred to any one-component distribution. We explore some consequences of the preferred birth velocity distribution: (i) roughly 50% of pulsars in the solar neighborhood will escape the Galaxy, while ~15% have velocities greater than 1000 km/s; (ii) an important low-velocity population exists that increases the fraction of neutron stars retained by globular clusters and is consistent with the number of old objects that accrete from the interstellar medium; (iii) under standard assumptions for supernova remnant expansion and pulsar spin-down, ~10% of pulsars younger than 20 kyr will appear to lie outside of their host remnants. Finally, we comment on the ramifications of our birth velocity distribution for binary survival and the population of inspiraling binary neutron stars relevant to some GRB models and potential sources for LIGO.

We discuss PSR J1740+1000, one of five pulsars recently discovered in a search of 470 square degrees at 430 MHz during the upgrade of the 305-m Arecibo telescope. The period of 154 ms and period derivative of 2.1 x 10^-14 s/s imply a spin-down age of 114 kyr that is smaller than 95% of all known pulsars. The youth and proximity of this pulsar make it a good candidate for detection at X-ray and gamma-ray energies. Its high Galactic latitude of 20.4 degrees suggests a very high velocity if the pulsar was born in the midplane of the Galaxy and if its kinematic age equals its spindown age. Interstellar scintillations, however, suggest a much lower velocity. We discuss possible explanations for this discrepancy, taking into account (a) possible birth sites away from the midplane; (b) contributions from the unmeasured radial velocity; (c) a kinematic age different from the spin-down age; and (d) biasing of the scintillation velocity by enhanced scattering from the North Polar Spur.

We discuss properties of the atmospheric water vapor above the high Andean plateau region known as the Llano de Chajnantor, in the Atacama Desert of northern Chile. A combination of radiometric and radiosonde measurements indicate that the median column of precipitable water vapor (PWV) above the plateau at an elevation of 5000 m is approximately 1.2 mm. The exponential scaleheight of the water vapor density in the median Chajnantor atmosphere is 1.13 km; the median PWV is 0.5 mm above an elevation of 5750 m. Both of these numbers appear to be lower at night. Annual, diurnal and other dependences of PWV and its scaleheight are discussed, as well as the occurrence of temperature inversion layers below the elevation of peaks surrounding the plateau. We estimate the background for infrared observations and sensitivities for broad band and high resolution spectroscopy. The results suggest that exceptional atmospheric conditions are present in the region, yielding high infrared transparency and high sensitivity for future ground-based infrared telescopes.

The region surrounding the Llano de Chajnantor, a high altitude plateau in the Atacama Desert in northern Chile, has caught the attention of the astronomical community for its potential as an observatory site. Combining high elevation and extremely low atmospheric water content, the Llano has been chosen as the future site of the Atacama Large Millimeter Array. We have initiated a campaign to investigate the astronomical potential of the region in the optical/infrared. Here, we report on an aspect of our campaign aimed at establishing a seeing benchmark to be used as a reference for future activities in the region. After a brief description of the region and its climate, we describe the results of an astronomical seeing campaign, carried out with a Differential Image Motion Monitor that operates at 0.5 micron wavelength. The seeing at the Llano level of 5000 m, measured over 7 nights in May 1998, yielded a median FWHM of 1.1". However, the seeing decreased to 0.7" at a modest 100 m gain above the plateau (Cerro Chico), as measured over 38 nights spread between July 1998 and October 2000. Neither of these represents the best seeing expected in the region; the set of measurements provides a reference base for simultaneous dual runs at Cerro Chico and at other sites of interest in the region, currently underway. A comparison between simultaneous measurements at Cerro Chico and Cerro Paranal indicates that seeing at Cerro Chico is about 12% better than at Paranal. The percentage of optically photometric nights in the Chajnantor region is about 60%, while that of nights useful for astronomical work is near 80%.

The Chandra X-ray Observatory recently discovered an X-ray point source near the center of Cassiopeia A, the youngest known Galactic supernova remnant. We have conducted a sensitive search for radio pulsations from this source with the Very Large Array, taking advantage of the high angular resolution of the array to resolve out the emission from the remnant itself. No convincing signatures of a dispersed, periodic source or of isolated dispersed pulses were found, whether for an isolated or a binary source. We derive upper limits of 30 and 1.3 mJy at 327 and 1435 MHz for the phase-averaged pulsed flux density from this source. The corresponding luminosity limits are lower than those for any pulsar with age less than 10^4 years. The sensitivities of our search to single pulses were 25 and 1.0 Jy at 327 and 1435 MHz. For comparison, the Crab pulsar emits roughly 80 pulses per minute with flux densities greater than 100 Jy at 327 MHz and 8 pulses per minute with flux densities greater than 50 Jy at 1435 MHz. These limits are consistent with the suggestion that the X-ray point source in Cas A adds to the growing number of neutron stars which are not radio pulsars.

We have made high-resolution, high-sensitivity dynamic spectra of a sample of strong pulsars at 430 MHz with the Arecibo radiotelescope. For 4 pulsars we find faint but sharply delineated features in the secondary spectra. These are examples of the previously observed ``crisscross'' or ``multiple drift slope'' phenomenon presumed to be due to multiple imaging of the pulsar by the interstellar medium. The unprecedented resolution and dynamic range of our observations allow a deeper level of analysis. Distances to the dominant scattering screen along the line of sight are determined and are shown to agree well with those inferred from other scintillation phenomena. Multiple imaging of the pulsar by the ISM is required. A compact central image surrounded by a faint scattering halo, roughly circularly symmetric, is consistent with the data. Scattering from filaments may also be consistent. The angular extent of the scattering material parallel to the direction of the pulsar velocity is roughly 5 mas, corresponding to a linear extent of about 2 AU. Further observations of these features should allow better discrimination between models and an identification of the scattering structures.

We study implications for the apparent alignment of the spin axes, proper-motions, and polarization vectors of the Crab and Vela pulsars. The spin axes are deduced from recent Chandra X-ray Observatory images that reveal jets and nebular structure having definite symmetry axes. The alignments indicate these pulsars were born either in isolation or with negligible velocity contributions from binary motions. We examine the effects of rotation and the conditions under which spin-kick alignment is produced for various models of neutron star kicks. If the kick is generated when the neutron star first forms by asymmetric mass ejection or/and neutrino emission, then the alignment requires that the protoneutron star possesses an original spin with period $P_s$ much less than the kick timescale, thus spin-averaging the kick forces. The kick timescale ranges from 100 ms to 10 s depending on whether the kick is hydrodynamically driven or neutrino-magnetic field driven. For hydrodynamical models, spin-kick alignment further requires the rotation period of an asymmetry pattern at the radius near shock breakout (>100 km) to be much less than ~100 ms; this is difficult to satisfy unless rotation plays a dynamically important role in the core collapse and explosion ($P_s\lo 1$ ms). Aligned kick and spin vectors are inherent to the slow process of asymmetric electromagnetic radiation from an off-centered magnetic dipole. We reassess the viability of this effect, correcting a factor of 4 error in Harrison and Tademaru's calculation that increases the size of the effect. To produce a kick velocity of order a few hundred km/s requires that the neutron star be born with an initial spin close to 1 ms and that spindown due to r-mode driven gravitational radiation be inefficient compared to standard magnetic braking.

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 the discovery of pulsars J0030+0451, J0711+0931, and J1313+0931 that were found in a search of 470 square degrees at 430 MHz using the 305m Arecibo telescope. The search has an estimated sensitivity for long period, low dispersion measure, low zenith angle, and high Galactic latitude pulsars of ~1 mJy, comparable to previous Arecibo searches. Spin and astrometric parameters for the three pulsars are presented along with polarimetry at 430 MHz. PSR J0030+0451, a nearby pulsar with a period of 4.8 ms, belongs to the less common category of isolated millisecond pulsars. We have measured significant polarization in PSR J0030+0451 over more than 50% of the period, and use these data for a detailed discussion of its magnetospheric geometry. Scintillation observations of PSR J0030+0451 provide an estimate of the plasma turbulence level along the line of sight through the local interstellar medium.

I use a method based on interstellar scintillations for discerning information about source sizes on scales less than one micro-arc sec. I use a comprehensive model for a pulsar signal, scintillated amplitude modulated noise, that includes source fluctuations and noise statistics. The method takes into account time-frequency averaging in the signal processing as well as effects due to source structure. The method is applied to interferometric visibility data on the Vela pulsar which show slightly less scintillation modulation than expected for a point source in the strong scattering regime. The decreased scintillation modulation is likely to be due exclusively to time-frequency averaging rather than from any source size effects. The implied upper limit on source extent, derived through Bayesian inference, is compared to predictions that take into account beaming from the relativistic plasma flow in neutron star magnetospheres. The upper limit for the transverse source size ($\lesssim 400$ km at 95% confidence for a pulsar distance of 0.5 kpc) is easily consistent with conventional models for radio emission regions in pulsar magnetospheres that place them well inside the light cylinder at only a few neutron-star radii from the star's surface.

Interstellar scintillation can be used to probe transverse sizes of radio sources on scales inaccessible to the nominal resolution of any terrestrial telescope, e.g. $\lesssim 10^{-6}$ arc sec. Methodology is presented that exploits this superresolution phenomenon for both single aperture and interferometer observations. The treatment applies to the saturated (strong-scattering) regime and holds for both thin screens and extended media. A general signal model for radio sources is presented, {\it scintillated amplitude modulated noise}, which applies to compact, incoherent synchrotron sources such as AGNs and gamma-ray burst sources and also to known, coherent sources such as masers and pulsars. The exact probability density function for measured intensities and interferometric visibilities is obtained by solving a general Fredholm problem. An approximate density function is also obtained by using the equivalent number of degrees of freedom in scintillation modulations. The scintillation modulation variance is presented, which includes the effects of source structure and time-bandwidth averaging in the signal processing. Two Bayesian methods are outlined for inferring the sizes of emission regions that use first order statistics of the intensity and visibility. Intensity cross-correlation methods for inferring source sizes are also given. Intensity interferometry in the the radio context is compared to the optical intensity interferometry of Hanbury-Brown and Twiss.

We apply a likelihood analysis to pulsar detections, pulsar upper limits, and diffuse background measurements from the OSSE and EGRET instruments on the Compton Gamma Ray Observatory to constrain the luminosity law for gamma-ray pulsars and some properties of the gamma-ray pulsar population. We find that the dependence of luminosity on spin period and dipole magnetic field is much steeper at OSSE than at EGRET energies (50-200 keV and >100 MeV, respectively), suggesting that different emission mechanisms are responsible for low- and high-energy gamma-ray emission. Incorporating a spin-down model and assuming a pulsar spatial distribution, we estimate the fraction of the Galactic gamma-ray background due to unidentified pulsars and find that pulsars may be an important component of the OSSE diffuse flux, but are most likely not important at EGRET energies. Using measurements of the diffuse background flux from these instruments, we are able to place constraints on the braking index, initial spin period, and magnetic field of the Galactic pulsar population. We are also able to constrain the pulsar birthrate to be between 1/(25 yr) and 1/(500 yr). Our results are based on a large gamma-ray beam, but they do not scale in a simple way with beam size. With our assumed beam size, the implied gamma-ray efficiency for the EGRET detections is no more than 20%. We estimate that about 20 of the 169 unidentified EGRET sources are probably gamma-ray pulsars. We use our model to predict the pulsar population that will be seen by future gamma-ray instruments and estimate that GLAST will detect roughly 750 gamma-ray pulsars as steady sources, only 120 of which are currently known radio pulsars.

We derive an upper limit of 3 mJy (95% confidence) for the flux density at 317 MHz of the Geminga pulsar (J0633+1746). Our results are based on 7 hours of fast-sampled VLA data, which we averaged synchronously with the pulse period using a period model based on CGRO/EGRET gamma-ray data. Our limit accounts for the fact that this pulsar is most likely subject to interstellar scintillations on a timescale much shorter than our observing span. Our Bayesian method is quite general and can be applied to calculate the fluxes of other scintillated sources. We also present a Bayesian technique for calculating the flux in a pulsed signal of unknown width and phase. Comparing our upper limit of 3 mJy with the quoted flux density of Geminga at 102 MHz, we calculate a lower limit to its spectral index of 2.7. We discuss some possible reasons for Geminga's weakness at radio wavelengths, and the likelihood that many of the unidentified EGRET sources are also radio-quiet or radio-weak Geminga-like pulsars.

We show that the crust-core interface in neutron stars acts as a potential barrier to the peripheral neutron vortices approaching the interface in the model in which these are coupled to the proton vortex clusters. This elementary barrier because of the interaction of vortex magnetic flux with the Meissner currents set up by the crustal magnetic field at the interface. The dominant part of the force is derived from to the cluster-interface interaction. As a result of the stopping of the continuous neutron vortex current through the interface, angular momentum is stored in the superfluid layers in the vicinity of the crust-core interface during the interglitch period. Discontinuous annihilation of proton vortices on the boundary restores the neutron vortex current and spins up the observable crust on short time-scales, leading to a glitch in the spin characteristics of a pulsar.

From considerations of spin evolution and kinematics in the galactic potential, we argue that the pulsars B1913+16, B1534+12, and B2127+11C may be younger than previously assumed, and we find that a lower bound on the formation and merger rate of close double neutron star binaries is $10^{-6.7} f_b^{-1} f_d^{-1}$ yr$^{-1}$, where $f_b$ is the beaming fraction and $f_d$ accounts for the possibility that the known NS-NS binaries are atypical of the underlying population (e.g., if most such binaries are born with shorter orbital periods). If we assume no prior knowledge of the detectable lifetimes of such systems, the rate could be as large as $\simeq 10^{-5.0} f_b^{-1} f_d^{-1}$ yr$^{-1}$. From both plausible bounds on $f_b$ and $f_d$, and a revision of the independently derived limit proposed by Bailes (1996), we argue that a firm upper bound is $10^{-4}$ yr$^{-1}$. We also present a unifying empirical overview of the spin-up of massive binary pulsars based on their distribution in spin period $P$ and spin-down rate $\dot P$, finding evidence for two distinct spin-up processes, one dominated by disk accretion, the other by wind accretion. We argue that the positions of binary pulsars in the $P$-$\dot P$ diagram can be understood if (1) there exists a Galactic population of pulsars in double neutron star systems with combinations of spin and orbital periods that have prevented their detection in surveys to date; (2) recycled pulsars in wide-orbit binaries are not born near the canonical spin-up line in the $P$-$\dot P$ diagram because they were predominantly spun up through wind accretion; and (3) there exists a disfavored evolutionary endpoint for radio pulsars with spin periods 5--30 ms and $\dot P > 10^{-19}$ s-s$^{-1}$.

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.

We consider the precession of isolated neutron stars in which superfluid is not pinned to the stellar crust perfectly. In the case of perfect pinning, Shaham (1977) showed that there are no slowly oscillatory, long-lived modes. When the assumption of perfect pinning is relaxed, new modes are found that can be long-lived, but are expected to be damped rather than oscillatory, unless the drag force on moving superfluid vortex lines has a substantial component perpendicular to the direction of relative motion. The response of a neutron star to external torques, such as the spindown torque, is also treated. We find that when computing the response of a star to perturbations, assuming perfect coupling of superfluid to normal matter from the start can miss some effects.

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.

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.

We discuss the spin evolution of pulsars in the case where a superfluid component of the star is coupled to the observable crust on long, spindown timescales. The momentum transfer from the superfluid interior results in an apparent decay of the external torque and, after a dramatic increase, to an asymptotic decrease of the generic value of the braking index, e.g. $n=3$, to values $n\sim 2.5$ if the magnetic field of the star does not decay over its lifetime. In the case where an exponential decay of the magnetic field towards a residual value occurs, the star undergoes a spin-up phase after which it could emerge in the millisecond sector of the $P$-$\dot P$ diagram.

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 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 show that the neutron star's crust-core interface acts as a potential barrier on the peripheral neutron vortices approaching this interface and thus prevents their continuous decay on it in the course of equilibrium state deceleration. The barrier arises due to the interaction of vortex magnetic flux with the Meissner currents set up by the crustal magnetic field at the interface. When the non-balanced part of the Magnus force reaches the value at which the vortices are able to annihilate at the interface, the rapid transfer of angular momentum from the superfluid spins-up the observable crust on short dynamical coupling times.

We use astrometric, distance and spindown data on pulsars to: (1) estimate three-dimensional velocity components, birth distances from the galactic plane, and ages of individual objects; (2) determine the distribution of space velocities and the scale height of pulsar progenitors; (3) test spindown laws for pulsars; (4) test for correlations between space velocities and other pulsar parameters; and (5) place empirical requirements on mechanisms than can produce high velocity neutron stars. Our approach incorporates measurement errors, uncertainties in distances, deceleration in the Galactic potential, and differential galactic rotation. We find that the scale height of the progenitors is approximately 0.13 kpc, that the 3D velocities are distributed in two components with characteristic speeds of 175(+20,-30) km/s and 700(+200,-150) km/s representing 83% and 17% of the population respectively. These results are insensitive to the explicit relation of chronological and spindown ages. We infer that the most probable chronological ages are typically smaller than conventional spindown ages by factors as large as two. We assess mechanisms for producing high-velocity neutron stars in view of the derived velocity distribution function.

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.

We have detected weak soft X-ray emission from the Pulsar Wind Nebula trailing the high velocity star PSR 2224+65 (the `Guitar Nebula'). This X-ray flux gives evidence of \gamma~10^7 eV particles in the pulsar wind and constrains the properties of the post-shock flow. The X-ray emission is most easily understood if the shocked pulsar wind is partly confined in the nebula and if magnetic fields in this zone can grow to near equipartition values.

We study the field millisecond pulsar (MSP) population to infer its intrinsic distribution in spin period and luminosity and to determine its spatial distribution within the Galaxy. Our likelihood analysis on data from extant surveys (22 pulsars with periods <20 ms) accounts for important selection effects. We infer a minimum period cutoff P(min) > 0.65 ms (99% confidence), a period distribution proportional to P^{-2.0 +- 0.33} and a pseudo-luminosity distribution proportional to L_p^{-2.0 +- 0.2} (where L_p = flux density times distance^2, for L_p >= 1.1 mJy kpc^2). We find a vertical scale height 0.65{+0.16,-0.12} kpc. We use our results to estimate the total number and birthrate of MSPs in the disk of the Galaxy. We limit the density contribution of a diffuse halo-like component to <1% of the midplane value. The MSP velocity dispersion is smaller that that of young, long-period pulsars by about a factor of 5. Our best estimate of the 1D velocity kick that is unique to MSP evolution is approximately 40 km s^-1. We discuss the evolutionary relationship of MSPs and low-mass X-ray binaries and prospects for future searches for MSPs.

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 made angular broadening measurements in the satellite line of the OH molecule at 1612 MHz of 8 OH/IR stars located near the Galactic Center with the Very Large Array. Early observations with the Very Long Baseline Array had revealed a region of pronounced scattering approximately centered on SgrA*. These new observations show that the compact line emission from these stars is scattered by an amount that is two to three times larger than had been previously seen. Additionally, we measure non-circular scattering disks for these objects, evidence that the electron density variations are distributed anisotropically. The present data suggest that the enhanced turbulence is located near the Galactic Center. However, we cannot completely rule out that the scattering is caused by an unrelated region along the line of sight.