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.

Recent work has shown that density profiles in the outskirts of dark matter halos can become extremely steep over a narrow range of radius. This behavior is produced by splashback material on its first apocentric passage after accretion. We show that the location of this splashback feature may be understood quite simply, from first principles. We present a simple model, based on spherical collapse, that accurately predicts the location of splashback without any free parameters. The important quantities that determine the splashback radius are accretion rate and redshift.

Satellites within simulated massive clusters are significantly spatially correlated with each other, even when those satellites are not gravitationally bound to each other. This correlation is produced by satellites that entered their hosts relatively recently, and is undetectable for satellites that have resided in their hosts for multiple dynamical timescales. Therefore, a measurement of clustering statistics of cluster satellites may be used to determine the typical accretion redshifts of those satellites into their observed hosts. We argue that such measurements may be used to determine the fraction of satellite galaxies that were quenched by their current hosts, thereby discriminating among models for quenching of star formation in satellite galaxies.

We investigate further a model of the accreting millisecond X-ray pulsars we proposed earlier. In this model, the X-ray-emitting regions of these pulsars are near their spin axes but move. This is to be expected if the magnetic poles of these stars are close to their spin axes, so that accreting gas is channeled there. As the accretion rate and the structure of the inner disk vary, gas is channeled along different field lines to different locations on the stellar surface, causing the X-ray-emitting areas to move. We show that this "nearly aligned moving spot model" can explain many properties of the accreting millisecond X-ray pulsars, including their generally low oscillation amplitudes and nearly sinusoidal waveforms; the variability of their pulse amplitudes, shapes, and phases; the correlations in this variability; and the similarity of the accretion- and nuclear-powered pulse shapes and phases in some. It may also explain why accretion-powered millisecond pulsars are difficult to detect, why some are intermittent, and why all detected so far are transients. This model can be tested by comparing with observations the waveform changes it predicts, including the changes with accretion rate.

We have shown previously that many of the properties of persistent accretion-powered millisecond pulsars can be understood if their X-ray emitting areas are near their spin axes and move as the accretion rate and structure of the inner disk vary. Here we show that this "nearly aligned moving spot model" may also explain the intermittent accretion-powered pulsations that have been detected in three weakly magnetic accreting neutron stars. We show that movement of the emitting area from very close to the spin axis to about 10 degrees away can increase the fractional rms amplitude from less than about 0.5 percent, which is usually undetectable with current instruments, to a few percent, which is easily detectable. The second harmonic of the spin frequency usually would not be detected, in agreement with observations. The model produces intermittently detectable oscillations for a range of emitting area sizes and beaming patterns, stellar masses and radii, and viewing directions. Intermittent oscillations are more likely in stars that are more compact. In addition to explaining the sudden appearance of accretion-powered millisecond oscillations in some neutron stars with millisecond spin periods, the model explains why accretion-powered millisecond oscillations are relatively rare and predicts that the persistent accretion-powered millisecond oscillations of other stars may become undetectable for brief intervals. It suggests why millisecond oscillations are frequently detected during the X-ray bursts of some neutron stars but not others and suggests mechanisms that could explain the occasional temporal association of intermittent accretion-powered oscillations with thermonuclear X-ray bursts.