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.

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.

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.

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.

We present a detailed analysis of the International Pulsar Timing Array (IPTA) Second Mock Data Challenge. We tested our analysis methods using the open datasets, and then analyzed the closed datasets. In both the open and the closed datasets, we were able to detect some, but not all, of the injected gravitational wave signals. This work presents two search cases that are not well explored in the pulsar timing array (PTA) literature: a simultaneous search for a stochastic GW background and an individual loud super-massive black hole binary (SMBHB) and a simultaneous search for two SMBHB sources. While we have constructed a cohesive framework for performing these GW searches, our analyses required fine-tuning of the sampling method used in order to appropriately converge. Given the nature of real PTA data in which multiple sources will be present in data, improved techniques will be required in the future to accurately detect and characterize these GW signals.

Rotating radio transients (RRATs) are a sub-class of pulsars characterized by sporadic emission and thus can generally only be studied by analysis of their single-pulses. Here we present a single-pulse analysis using 11 years of timing data at 1400~MHz of three RRATs, PSRs~J1819$-$1458, J1317$-$5759, and J1913$+$1330. We perform a spectral analysis on the single-pulses of these RRATs for the first time, finding their mean spectral indices to be $-1.1 \pm 0.1$, $-0.6 \pm 0.1$, and $-1.2 \pm 0.2$ respectively, within the known range of pulsar spectral indices. We find no evidence for narrowband features as seen for FRB~121102. However, we find the spread of single-pulse spectral indices for these RRATs (ranging from $-7$ to $+4$) to be larger than has been seen in other pulsars, with the exception of the Crab pulsar. We also analyze the time between detected pulses, or wait-time, and find that the pulses are not random and cluster around wait-times of a few pulse periods as well as $\sim 25$ pulse periods for PSRs~J1819$-$1458 and J1317$-$5759. Additionally we find that there is no correlation between the wait-time and pulse flux density. Finally we find that the distribution of the pulse energy for PSRs~J1317$-$5759 and J1913$+$1330 are log-normal, while that of PSR~J1819$-$1458 is log-normal with possible evidence of an additional power-law component.