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.

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.

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

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

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

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.

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.

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

We report on Westerbork Synthesis Radio Telescope observations of the Crab pulsar at 350 MHz from 2012 November 24 until 2015 June 21. During this period we consistently observe variations in the pulse profile of the Crab. Both variations in the scattering width of the pulse profile as well as delayed copies, also known as echoes, are seen regularly. These observations support the classification of two types of echoes: those that follow the truncated exponential shape expected for the thin-screen scattering approximation, and echoes that show a smoother, more Gaussian shape. During a sequence of high-cadence observations in 2015, we find that these non-exponential echoes evolve in time by approaching the main pulse and interpulse in phase, overlapping the main pulse and interpulse, and later receding. We find a pulse scatter-broadening time scale, $\tau$, scaling with frequency as $\nu^{\alpha}$, with $\alpha=-3.9\pm0.5$, which is consistent with expected values for thin-screen scattering models

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

Scintillation arcs have become a powerful tool for exploring scattering in the ionized interstellar medium. There is accumulating evidence that the scattering from many pulsars is extremely anisotropic resulting in highly elongated, linear brightness functions. We present a three-frequency (327~MHz, 432~MHz, 1450~MHz) Arecibo study of scintillation arcs from one nearby, bright, high-velocity pulsar, PSR~B1133+16. We show that a one-dimensional (1D), linear brightness function is in good agreement with the data at all three observing frequencies. We use two methods to explore the broadening of the 1D brightness function $B(\theta)$ as a function of frequency: 1) crosscuts of the forward arc at constant delay and 2) a 1D modeling of $B(\theta)$ using a comparison between model and observed secondary spectrum as a goodness-of-fit metric. Both methods show that the half-power width of $B(\theta)$ deviates from the expected dependence $\propto \nu^{-a}$, where $\nu$ is the observing frequency . Our estimates of $a$ have moderately large uncertainties but imply $a \lesssim1.8$, and so are inconsistent with the expected $a = 2.0$ for plasma refraction or $a = 2.2$ for Kolmogorov turbulence. In addition the shape of $B(\theta)$ cuts off more steeply than predicted for Kolmogorov turbulence. Ultimately, we conclude that the underlying physics of the broadening mechanism remains unexplained. Our results place the scattering screen at a distance that is broadly consistent with an origin at the boundary of the Local Bubble.

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

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

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

We present high-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 have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar medium (IISM) are the most likely cause. The putative structures in the IISM causing the profile variation are roughly half-way between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We consider different models for the structures. Under the assumption of spherical symmetry, their implied average electron density is $\overline{n}_e \sim 100$ cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than the average pulsar population visible to LOFAR, similar profile variations would not have been identified in most pulsars, suggesting that subtle profile variations in low-frequency profiles might be more common than we have observed to date. Systematic studies of these variations at low frequencies can provide a new tool to investigate the proprieties of the IISM and the limits to the precision of pulsar timing.

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.

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.

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.

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.

Light travel time changes due to gravitational waves may be detected within the next decade through precision timing of millisecond pulsars. Removal of frequency-dependent interstellar medium (ISM) delays due to dispersion and scattering is a key issue in the detection process. Current timing algorithms routinely correct pulse times of arrival (TOAs) for time-variable delays due to cold plasma dispersion. However, none of the major pulsar timing groups correct for delays due to scattering from multi-path propagation in the ISM. Scattering introduces a frequency-dependent phase change in the signal that results in pulse broadening and arrival time delays. Any method to correct the TOA for interstellar propagation effects must be based on multi-frequency measurements that can effectively separate dispersion and scattering delay terms from frequency-independent perturbations such as those due to a gravitational wave. Cyclic spectroscopy, first described in an astronomical context by Demorest (2011), is a potentially powerful tool to assist in this multi-frequency decomposition. As a step toward a more comprehensive ISM propagation delay correction, we demonstrate through a simulation that we can accurately recover impulse response functions (IRFs), such as those that would be introduced by multi-path scattering, with a realistic signal-to-noise ratio. We demonstrate that timing precision is improved when scatter-corrected TOAs are used, under the assumptions of a high signal-to-noise and highly scattered signal. We also show that the effect of pulse-to-pulse "jitter" is not a serious problem for IRF reconstruction, at least for jitter levels comparable to those observed in several bright pulsars.

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

High-precision timing of millisecond pulsars (MSPs) over years to decades is a promising technique for direct detection of gravitational waves at nanohertz frequencies. Time-variable, multi-path scattering in the interstellar medium is a significant source of noise for this detector, particularly as timing precision approaches 10 ns or better for MSPs in the pulsar timing array. For many MSPs the scattering delay above 1 GHz is at the limit of detectability; therefore, we study it at lower frequencies. Using the LOFAR (LOw-Frequency ARray) radio telescope we have analyzed short (5-20 min) observations of three MSPs in order to estimate the scattering delay at 110-190 MHz, where the number of scintles is large and, hence, the statistical uncertainty in the scattering delay is small. We used cyclic spectroscopy, still relatively novel in radio astronomy, on baseband-sampled data to achieve unprecedented frequency resolution while retaining adequate pulse phase resolution. We detected scintillation structure in the spectra of the MSPs PSR B1257+12, PSR J1810+1744, and PSR J2317+1439 with diffractive bandwidths of $6\pm 3$, $2.0\pm 0.3$, and $\sim 7$ kHz, respectively, where the estimate for PSR J2317+1439 is reliable to about a factor of 2. For the brightest of the three pulsars, PSR J1810+1744, we found that the diffractive bandwidth has a power-law behavior $\Delta\nu_d \propto \nu^{\alpha}$, where $\nu$ is the observing frequency and $\alpha = 4.5\pm 0.5$, consistent with a Kolmogorov inhomogeneity spectrum. We conclude that this technique holds promise for monitoring the scattering delay of MSPs with LOFAR and other high-sensitivity, low-frequency arrays like SKA-Low.

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.

Time variable delays due to radio wave propagation in the ionized interstellar medium are a substantial source of error in pulsar timing array efforts. We describe the physical origin of these effects, discussing dispersive and scattering effects separately. Where possible, we give estimates of the magnitude of timing errors produced by these effects and their scaling with radio frequency. Although there is general understanding of the interstellar medium propagation errors to be expected with pulsar timing array observations, detailed comparison between theory and practice is still in its infancy, particularly with regard to scattering effects.

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 compare signal and noise for observations of the scintillating pulsar B0834+06, using very-long baseline interferometry and a single-dish spectrometer. Comparisons between instruments and with models suggest that amplitude variations of the pulsar strongly affect the amount and distribution of self-noise. We show that noise follows a quadratic polynomial with flux density, in spectral observations. Constant coefficients, indicative of background noise, agree well with expectation; whereas second-order coefficients, indicative of self-noise, are about 3 times values expected for a pulsar with constant on-pulse flux density. We show that variations in flux density during the 10-sec integration account for the discrepancy. In the secondary spectrum, about 97% of spectral power lies within the pulsar's typical scintillation bandwidth and timescale; an extended scintillation arc contains about 3%. For a pulsar with constant on-pulse flux density, noise in the dynamic spectrum will appear as a uniformly-distributed background in the secondary spectrum. We find that this uniform noise background contains 95% of noise in the dynamic spectrum for interferometric observations; but only 35% of noise in the dynamic spectrum for single-dish observations. Receiver and sky dominate noise for our interferometric observations, whereas self-noise dominates for single-dish. We suggest that intermittent emission by the pulsar, on timescales < 300 microseconds, concentrates self-noise near the origin in the secondary spectrum, by correlating noise over the dynamic spectrum. We suggest that intermittency sets fundamental limits on pulsar astrometry or timing. Accounting of noise may provide means for detection of intermittent sources, when effects of propagation are unknown or impractical to invert.

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.

The dynamic spectrum of a radio pulsar is an in-line digital hologram of the ionised interstellar medium. It has previously been demonstrated that such holograms permit image reconstruction, in the sense that one can determine an approximation to the complex electric field values as a function of Doppler-shift and delay, but to date the quality of the reconstructions has been poor. Here we report a substantial improvement in the method which we have achieved by simultaneous optimisation of the thousands of coefficients that describe the electric field. For our test spectrum of PSR B0834+06 we find that the model provides an accurate representation of the data over the full 63 dB dynamic range of the observations: residual differences between model and data are noise-like. The advent of interstellar holography enables detailed quantitative investigation of the interstellar radio-wave propagation paths for a given pulsar at each epoch of observation; we illustrate this using our test data which show the scattering material to be structured and highly anisotropic. The temporal response of the medium exhibits a scattering tail out to beyond 100 microsec and a pulse arrival time measurement at this frequency and this epoch of observation would be affected by a mean delay of 15 microsec due to multipath propagation in the interstellar medium.

We present a highly personal and biased review of the observational portion of this meeting. We cover Tiny-Scale Atomic, Molecular, and Ionized Structure (TSAS, TSMS, TSIS), emphasizing the physical aspects involved when the structures are considered discrete. TSAS includes time- and angular-variable 21-cm and optical absorption lines. TSMS includes the Heithausen revolutionary component called Small Area Molecular Structures. TSIS includes classical scintillation and discrete structures that produce Arclets, Intra-Day and Intra-Hour Variables, and Extreme Scattering Events. We conclude with reflections on the relationship of Tiny-Scale Structure to the Local Bubble, TSAS to TSIS, and the use of globular clusters as an illuminating backdrop.

The LOw Frequency ARray, LOFAR, will have the sensitivity, bandwidth, frequency range and processing power to revolutionise low-frequency pulsar studies. We present results of simulations that indicate that a LOFAR survey will find approximately 1500 new pulsars. These new pulsars will give us a much better understanding of the low end of the luminosity function and thus allow for a much more precise estimate of the true total local pulsar population. The survey will also be very sensitive to the ultra-steep spectrum pulsars, RRATs, and the pulsed radio emission from objects like Geminga and AXPs. We will also show that by enabling us to observe single pulses from hundreds of pulsars, including many millisecond pulsars, LOFAR opens up new possibilities for the study of emission physics.

Topics of the Colloquium: a) Interplanetary scintillation b) Interstellar scintillation c) Modeling and physical origin of the interplanetary and the interstellar plasma turbulence d) Scintillation as a tool for investigation of radio sources e) Seeing through interplanetary and interstellar turbulent media Ppt-presentations are available on the Web-site: http://www.prao.ru/conf/Colloquium/main.html

Pulsar dynamic spectra exhibit high visibility fringes arising from interference between scattered radio waves. These fringes may be random or highly ordered patterns, depending on the nature of the scattering or refraction. Here we consider the possibility of decomposing pulsar dynamic spectra -- which are intensity measurements -- into their constituent scattered waves, i.e. electric field components. We describe an iterative method of achieving this decomposition and show how the algorithm performs on data from the pulsar B0834+06. The match between model and observations is good, although not formally acceptable as a representation of the data. Scattered wave components derived in this way are immediately useful for qualitative insights into the scattering geometry. With some further development this approach can be put to a variety of uses, including: imaging the scattering and refracting structures in the interstellar medium; interstellar interferometric imaging of pulsars at very high angular resolution; and mitigating pulse arrival time fluctuations due to interstellar scattering.

We have detected a strong deflection of radio waves from the pulsar PSR B0834+06 in scintillation observations. Interference between the undeflected pulsar image and deflected subimages allows single dish interferometry of the interstellar medium with sub-milliarcsecond resolution. We infer the presence of scattering structure(s) similar to those that are thought to cause Extreme Scattering Events in quasar flux monitoring programs: size ~ 0.2 AU (an angular size of 0.1 mas) with an electron density of > 10^3 compared to the warm ionized medium. The deflectors are nearly stationary in a scattering screen that is thin (< 5% of the pulsar-observer distance in extent), is located 70% of the way from the Earth to the pulsar, and has been seen consistently in observations dating back 20 years. The pulsar scans the scattering screen at a velocity of 110 km/s with a detection radius of 15 mas. Pulsar observations such as these -- particularly with a new generation of low-frequency radio telescopes with large collecting areas -- hold promise for improving constraints on the poorly understood physical characteristics and space density of the deflecting structures. Such observations may also prove useful in correcting deviations the deflectors produce in high-precision timing of millarcsecond pulsars.

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.

Pulsar dynamic spectra sometimes show organised interference patterns; these patterns have been shown to have power spectra which often take the form of parabolic arcs, or sequences of inverted parabolic arclets whose apexes themselves follow a parabolic locus. Here we consider the interpretation of these arc and arclet features. We give a statistical formulation for the appearance of the power spectra, based on the stationary phase approximation to the Fresnel-Kirchoff integral. We present a simple analytic result for the power-spectrum expected in the case of highly elongated images, and a single-integral analytic formulation appropriate to the case of axisymmetric images. Our results are illustrated in both the ensemble-average and snapshot regimes. Highly anisotropic scattering appears to be an important ingredient in the formation of the observed arclets.

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 have measured the properties of the diffractive scintillation toward the Vela pulsar under the extremely strong scattering conditions encountered at 660 MHz. We obtain a decorrelation bandwidth of $\nu_d = 244 \pm 4$ Hz and diffractive decorrelation timescale of $t_{\rm diff} = 3.3\pm 0.3$ s. Our measurement of the modulation indices $m=0.87\pm 0.003\pm 0.05$ and $m=0.93\pm 0.03 \pm 0.05$ (one for each polarization stream), are at variance with the modulation index of the Vela pulsar obtained at 2.3 GHz by Gwinn et al. (1997) {\it if} the deviation from a modulation index of unity is ascribed to a source size effect.

PSR B1259-63 is in a highly eccentric 3.4 yr orbit around the Be-star SS2883. The system is located in the direction of the Sagittarius-Carina spiral arm at a distance of $\sim$1.5 kpc. We have made scintillation observations of the pulsar far from periastron at 4.8 GHz and 8.4 GHz, determining the diffractive bandwidth, $\Delta \nu_{\rm d}$, and timescale, $\Delta t_{\rm d}$, at both frequencies. We find no dependency on orbital phase until within 30 days of periastron. The data indicate that the scintillation is caused, not by the circumstellar environment, but by an HII region within the Sagittarius-Carina spiral arm located at least three-quarters of the way to the pulsar. Close to periastron, when the line-of-sight to the pulsar intersects the disk of the Be star, the electron densities within the disk are sufficient to overcome the `lever-arm effect' and produce a reduction in the scintillation bandwidth by six orders of magnitude.