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.

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3).

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.

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing - the process of accounting for every rotation of a pulsar over long time periods - can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system (Lorimer & Kramer 2005). One of these, called relativistic Shapiro delay (Shapiro 1964), can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-year data set with recent orbital-phase-specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be $\mathbf{2.14^{+0.10}_{-0.09}}$ solar masses (68.3% credibility interval; 95.4% credibility interval is $\mathbf{2.14^{+0.20}_{-0.18}}$ solar masses). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

Recycled pulsars are old ($\gtrsim10^{8}$ yr) neutron stars that are descendants from close, interacting stellar systems. In order to understand their evolution and population, we must find and study the largest number possible of recycled pulsars in a way that is as unbiased as possible. In this work, we present the discovery and timing solutions of five recycled pulsars in binary systems (PSRs J0509$+$0856, J0709$+$0458, J0732$+$2314, J0824$+$0028, J2204$+$2700) and one isolated millisecond pulsar (PSR J0154$+$1833). These were found in data from the Arecibo 327-MHz Drift-Scan Pulsar Survey (AO327). All these pulsars have a low dispersion measure (DM) ($\lesssim 45 \, \rm{pc}\, cm^{-3}$), and have a DM-determined distance of $\lesssim$ 3 kpc. Their timing solutions, have data spans ranging from 1 to $\sim$ 7 years, include precise estimates of their spin and astrometric parameters, and for the binaries, precise estimates of their Keplerian binary parameters. Their orbital periods range from about 4 to 815 days and the minimum companion masses (assuming a pulsar mass of 1.4 $\rm{M_{\odot}}$) range from $\sim$ 0.06--1.11 $\rm{M_{\odot}}$. For two of the binaries we detect post-Keplerian parameters; in the case of PSR~J0709$+$0458 we measure the component masses but with a low precision, in the not too distant future the measurement of the rate of advance of periastron and the Shapiro delay will allow very precise mass measurements for this system. Like several other systems found in the AO327 data, PSRs J0509$+$0854, J0709$+$0458 and J0732$+$2314 are now part of the NANOGrav timing array for gravitational wave detection.

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.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.

We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.

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.

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.

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the most stable millisecond pulsars. Several known sources of noise on short timescales in single radio-pulsar observations are well described by a simple model of three components: template-fitting from a finite signal-to-noise ratio, pulse phase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellar scattering variations. Currently template-fitting errors dominate, but as radio telescopes push towards higher signal-to-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics of jitter becomes crucial for properly characterizing arrival-time uncertainties. We characterize the radio-frequency dependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for its frequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations of rms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter is typically ~1% of pulse phase. The average frequency dependence for all pulsars using a power-law model has index -0.42. We investigate the jitter variations for the interpulse of PSR~B1937+21 and find no significant deviations from the main pulse rms jitter. We also test the time-variation of jitter in two pulsars and find that systematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly model jitter as a significant component of the noise within the detector.

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.

Current single dish, low-frequency radio pulsar surveys provide efficient sky coverage, but poor localization of new discoveries. Here, we describe a practical technique for rapidly localizing pulsars discovered in these surveys with on-the-fly mapping and provide code to facilitate and formalize its implementation. As a proof of concept, we alter the positions of four test sources and use the Green Bank Telescope (GBT) 350 MHz receiver to recover source positions within $\approx$1$-$3' of their true values, compared to an 18' error radius for new discoveries. Achieving similar precision with a traditional gridding strategy using the GBT requires 2$-$3 times as much telescope time (including overhead), multiple receivers and relies on assumptions about the pulsars' spectral indices. For one of our test sources (PSR J1400$-$1431), this method revealed a discrepancy with the initial, published position, prompting additional follow-up and an improved timing solution. Rapid localization is important for improving data quality and providing flexibility in choice of center frequency for future timing observations - both of which facilitate evaluating new millisecond pulsars for potential inclusion in pulsar timing arrays.

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 present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.

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.

Black widows and redbacks are binary systems consisting of a millisecond pulsar in a close binary with a companion having matter driven off of its surface by the pulsar wind. X-rays due to an intra-binary shock have been observed from many of these systems, as well as orbital variations in the optical emission from the companion due to heating and tidal distortion. We have been systematically studying these systems in radio, optical and X-rays. Here we will present an overview of X-ray and optical studies of these systems, including new XMM-Newton and NuStar data obtained from several of them, along with new optical photometry.

In 2012, five high school students involved in the Pulsar Search Collaboratory discovered the millisecond pulsar PSR J1400$-$1431 and initial timing parameters were published in Rosen et al. (2013) a year later. Since then, we have obtained a phase-connected timing solution spanning five years, resolving a significant position discrepancy and measuring $\dot{P}$, proper motion, parallax, and a monotonic slope in dispersion measure over time. Due to PSR J1400$-$1431's proximity and significant proper motion, we use the Shklovskii effect and other priors to determine a 95% confidence interval for PSR J1400$-$1431's distance, $d=270^{+130}_{-80}$ pc. With an improved timing position, we present the first detection of the pulsar's low-mass white dwarf (WD) companion using the Goodman Spectrograph on the 4.1-m SOAR telescope. Deeper imaging suggests that it is a cool DA-type WD with $T_{\rm eff}=3000\pm100$ K and $R/R_\odot=(2.19\pm0.03)\times10^{-2}\,(d/270\,{\rm pc})$. We show a convincing association between PSR J1400$-$1431 and a $\gamma$-ray point source, 3FGL J1400.5$-$1437, but only weak (3.3-$\sigma$) evidence of pulsations after folding $\gamma$-ray photons using our radio timing model. We detect an X-ray counterpart with XMM-Newton but the measured X-ray luminosity ($1\times10^{29}$ ergs s$^{-1}$) makes PSR J1400$-$1431 the least X-ray luminous rotation-powered millisecond pulsar (MSP) detected to date. Together, our findings present a consistent picture of a nearby ($d\approx230$ pc) MSP in a 9.5-day orbit around a cool, $\sim$0.3 M$_\odot$ WD companion, with orbital inclination, $i\gtrsim60^\circ$.

PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.

We present the first X-ray observations of three recently discovered millisecond pulsars (MSPs) with interesting characteristics: PSR J0337+1715, PSR J0636+5129, and PSR J0645+5158. PSR J0337+1715 is a fast-spinning, bright, and so-far unique MSP in a hierarchical triple system with two white dwarf (WD) companions. PSR J0636+5129 is a MSP in a very tight 96-min orbit with a low-mass, 8 $M_J$ companion. PSR J0645+5158 is a nearby, isolated MSP with a very small duty cycle (1-2%), which has led to its inclusion in high-precision pulsar timing programs. Using data from XMM-Newton, we have analyzed X-ray spectroscopy for these three objects, as well as optical/ultraviolet photometry for PSR J0337+1715. The X-ray data for each are largely consistent with expectations for most MSPs with regards to the ratios of thermal and non-thermal emission. We discuss the implications of these data on the pulsar population, and prospects for future observations of these pulsars.

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 consider the X-ray properties of the redback class of eclipsing millisecond pulsars. These are transitional systems between accreting low-mass X-ray binaries and binary millisecond pulsars orbiting white dwarfs, and hence their companions are non-degenerate and nearly Roche-lobe filling. The X-ray luminosity seems to scale with the fraction of the pulsar sky subtended by the companion, suggesting the shock region is not much larger than the companion, which is supported by modeling of the orbital light curves. The typical X-ray photon spectral index is $\sim 1$ and the typical 0.3-8 keV X-ray efficiency, assuming a shock size on the order of the companion's Roche lobe cross-section, is on the order of 10%. We present an overview of previous investigations, and present new observations of two redbacks, a Chandra observation of PSR J1628$-$3205 and a XMM-Newton observation of PSR J2129$-$0429. The latter shows a clearly double peaked orbital light curve with variation of the non-thermal flux by a factor of $\sim 11$, with peaks around orbital phases 0.6 and 0.9. We suggest the magnetic field of the companion plays a significant role in the X-ray emission from intrabinary shocks in redbacks.

Eclipsing millisecond pulsars in close ($P_b < 1$~day) binary systems provide a different view of pulsar winds and shocks than do isolated pulsars. Since 2009, the numbers of these systems known in the Galactic field has increased enormously. We have been systematically studying many of these newly discovered systems at multiple wavelengths. Typically, the companion is nearly Roche-lobe filling and heated by the pulsar which drives mass loss from the companion. The pulsar wind shocks with this material just above the surface of the companion. We discuss various observational properties of this shock, including radio eclipses, orbitally modulated X-ray emission, and the potential for $\gamma$-ray emission. Redbacks, whose companions are likely non-degenerate and significantly more massive, generally have more luminous shocks than black widows which have very low mass companions. This is expected since the more massive redback companions intercept a greater fraction of the pulsar wind. We also compare these systems to accreting millisecond pulsars, which may be progenitors of black widows and in some cases can pass back and forth between redback and accretion phases.

We describe the first X-ray observations of five short orbital period ($P_B < 1$ day), $\gamma$-ray emitting, binary millisecond pulsars. Four of these, PSRs J0023+0923, J1124$-$3653, J1810+1744, and J2256$-$1024 are `black-widow' pulsars, with degenerate companions of mass $\ll0.1 M_{\odot}$, three of which exhibit radio eclipses. The fifth source, PSR J2215+5135, is an eclipsing `redback' with a near Roche-lobe filling $\sim$0.2 solar mass non-degenerate companion. Data were taken using the \textit{Chandra X-Ray Observatory} and covered a full binary orbit for each pulsar. Two pulsars, PSRs J2215+5135 and J2256$-$1024, show significant orbital variability while PSR J1124$-$3653 shows marginal orbital variability. The lightcurves for these three pulsars have X-ray flux minima coinciding with the phases of the radio eclipses. This phenomenon is consistent with an intrabinary shock emission interpretation for the X-rays. The other two pulsars, PSRs J0023+0923 and J1810+1744, are fainter and do not demonstrate variability at a level we can detect in these data. All five spectra are fit with three separate models: a power-law model, a blackbody model, and a combined model with both power-law and blackbody components. The preferred spectral fits yield power-law indices that range from 1.3 to 3.2 and blackbody temperatures in the hundreds of eV. The spectrum for PSR J2215+5135 shows a significant hard X-ray component, with a large number of counts above 2 keV, which is additional evidence for the presence of intrabinary shock emission and is similar to what has been detected in the low-mass X-ray binary to millisecond pulsar transition object PSR J1023+0038.

We describe the first X-ray observations of binary millisecond pulsars PSRs J0023+0923, J1810+1744, J2215+5135, and J2256-1024. All four are Fermi gamma-ray sources and three are 'black-widow' pulsars, with companions of mass < 0.1 solar masses. Data were taken using the Chandra X-Ray Observatory and covered a full binary orbit for each pulsar. Two pulsars, PSRs J2215+5135 and J2256-1024, show significant orbital variability and X-ray flux minima at the times of eclipses observed at radio wavelengths. This phenomenon is consistent with intrabinary shock emission characteristic of black-widow pulsars. The other two pulsars, PSRs J0023+0923 and J1810+1744, do not demonstrate significant variability, but are fainter than the other two sources. Spectral fits yield power-law indices that range from 1.4 to 2.3 and blackbody temperatures in the hundreds of eV. The spectrum for PSR J2215+5135 shows a significant hard X-ray component (41% of counts are above 2 keV), which is additional evidence for the presence of intrabinary shock emission.