Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.

Using commissioning data from the first year of operation of the Canadian Hydrogen Intensity Mapping Experiment's (CHIME) Pulsar backend system, we conduct a systematic analysis of the Faraday Rotation Measure (RM) of the northern hemisphere pulsars detected by CHIME. We present 55 new RMs as well as obtain improved RM uncertainties for 25 further pulsars. CHIME's low observing frequency and wide bandwidth between 400-800 MHz contribute to the precision of our measurements, whereas the high cadence observation provide extremely high signal-to-noise co-added data. Our results represent a significant increase of the pulsar RM census, particularly regarding the northern hemisphere. These new RMs are for sources that are located in the Galactic plane out to 10 kpc, as well as off the plane to a scale height of ~16 kpc. This improved knowledge of the Faraday sky will contribute to future Galactic large-scale magnetic structure and ionosphere modelling.

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.

We report on the detection of seven bursts from the periodically active, repeating fast radio burst (FRB) source FRB 180916.J0158+65 in the 300-400-MHz frequency range with the Green Bank Telescope (GBT). Emission in multiple bursts is visible down to the bottom of the GBT band, suggesting that the cutoff frequency (if it exists) for FRB emission is lower than 300 MHz. Observations were conducted during predicted periods of activity of the source, and had simultaneous coverage with the Low Frequency Array (LOFAR) and the FRB backend on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. We find that one of the GBT-detected bursts has potentially associated emission in the CHIME band (400-800 MHz) but we detect no bursts in the LOFAR band (110-190 MHz), placing a limit of $\alpha > -1.0$ on the spectral index of broadband emission from the source. We also find that emission from the source is severely band-limited with burst bandwidths as low as $\sim$40 MHz. In addition, we place the strictest constraint on observable scattering of the source, $<$ 1.7 ms, at 350 MHz, suggesting that the circumburst environment does not have strong scattering properties. Additionally, knowing that the circumburst environment is optically thin to free-free absorption at 300 MHz, we find evidence against the association of a hyper-compact HII region or a young supernova remnant (age $<$ 50 yr) with the source.

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

We report on the discovery and analysis of bursts from nine new repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 195 to 1380 pc cm$^{-3}$. We detect two bursts from three of the new sources, three bursts from four of the new sources, four bursts from one new source, and five bursts from one new source. We determine sky coordinates of all sources with uncertainties of $\sim$10$^\prime$. We detect Faraday rotation measures for two sources, with values $-20(1)$ and $-499.8(7)$ rad m$^{-2}$, that are substantially lower than the RM derived from bursts emitted by FRB 121102. We find that the DM distribution of our events, combined with the nine other repeaters discovered by CHIME/FRB, is indistinguishable from that of thus far non-repeating CHIME/FRB events. However, as previously reported, the burst widths appear statistically significantly larger than the thus far non-repeating CHIME/FRB events, further supporting the notion of inherently different emission mechanisms and/or local environments. These results are consistent with previous work, though are now derived from 18 repeating sources discovered by CHIME/FRB during its first year of operation. We identify candidate galaxies that may contain FRB 190303.J1353+48 (DM = 222.4 pc cm$^{-3}$).

We report on multi-band photometric observations of four binary millisecond pulsars with the Gran Telescopio Canarias. The observations led to detection of binary companions to PSRs J1630+3734, J1741+1351 and J2042+0246 in the Sloan g', r' and i' bands. Their magnitudes in the r' band are $\approx$24.4, 24.4 and 24.0, respectively. We also set a 3$\sigma$ upper limit on the brightness of the PSR J0557+1550 companion in the r' band of $\approx$25.6 mag. Combining the optical data with the radio timing measurements and white dwarf cooling models, we show that the detected companions are cool low-mass white dwarfs with temperatures and ages in the respective ranges of (4-7)$\times 10^3$ K and 2-5 Gyr. All the detected white dwarfs are found to likely have either pure hydrogen or mixed helium-hydrogen atmospheres.

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

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.

We report on the discovery of eight repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 103.5 to 1281 pc cm$^{-3}$. They display varying degrees of activity: six sources were detected twice, another three times, and one ten times. These eight repeating FRBs likely represent the bright and/or high-rate end of a distribution of infrequently repeating sources. For all sources, we determine sky coordinates with uncertainties of $\sim$10$^\prime$. FRB 180916.J0158+65 has a burst-averaged DM = $349.2 \pm 0.3$ pc cm$^{-3}$ and a low DM excess over the modelled Galactic maximum (as low as $\sim$20 pc cm$^{-3}$); this source also has a Faraday rotation measure (RM) of $-114.6 \pm 0.6$ rad m$^{-2}$, much lower than the RM measured for FRB 121102. FRB 181030.J1054+73 has the lowest DM for a repeater, $103.5 \pm 0.3$ pc cm$^{-3}$, with a DM excess of $\sim$ 70 pc cm$^{-3}$. Both sources are interesting targets for multi-wavelength follow-up due to their apparent proximity. The DM distribution of our repeater sample is statistically indistinguishable from that of the first 12 CHIME/FRB sources that have not repeated. We find, with 4$\sigma$ significance, that repeater bursts are generally wider than those of CHIME/FRB bursts that have not repeated, suggesting different emission mechanisms. Our repeater events show complex morphologies that are reminiscent of the first two discovered repeating FRBs. The repetitive behavior of these sources will enable interferometric localizations and subsequent host galaxy identifications.

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.

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

We report the detection of a single burst from the first-discovered repeating Fast Radio Burst source, FRB 121102, with CHIME/FRB, which operates in the frequency band 400-800 MHz. The detected burst occurred on 2018 November 19 and its emission extends down to at least 600 MHz, the lowest frequency detection of this source yet. The burst, detected with a significance of 23.7$\sigma$, has fluence 12$\pm$3 Jy ms and shows complex time and frequency morphology. The 34 ms width of the burst is the largest seen for this object at any frequency. We find evidence of sub-burst structure that drifts downward in frequency at a rate of -3.9$\pm$0.2 MHz ms$^{-1}$. Our best fit tentatively suggests a dispersion measure of 563.6$\pm$0.5 pc cm$^{-3}$, which is ${\approx}$1% higher than previously measured values. We set an upper limit on the scattering time at 500 MHz of 9.6 ms, which is consistent with expectations from the extrapolation from higher frequency data. We have exposure to the position of FRB 121102 for a total of 11.3 hrs within the FWHM of the synthesized beams at 600 MHz from 2018 July 25 to 2019 February 25. We estimate on the basis of this single event an average burst rate for FRB 121102 of 0.1-10 per day in the 400-800 MHz band for a median fluence threshold of 7 Jy ms in the stated time interval.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

We report detection of the binary companion to the millisecond pulsar J1741+1351 with the Gran Telescopio Canarias. The optical source position coincides with the pulsar coordinates and its magnitudes are g' = 24.84(5), r' = 24.38(4) and i' = 24.17(4). Comparison of the data with the white dwarf evolutionary models shows that the source can be a He-core white dwarf with a temperature of $\approx 6000$ K and a mass of $\approx 0.2$ M$_{\odot}$. The latter is in excellent agreement with the companion mass obtained from the radio timing solution for PSR J1741+1351.

We summarize the state of the art and future directions in using millisecond radio pulsars to test gravitation and measure intrinsic, fundamental parameters of the pulsar systems. As discussed below, such measurements continue to yield high-impact constraints on viable nuclear processes that govern the theoretically-elusive interior of neutron stars, and place the most stringent limits on the validity of general relativity in extreme environments. Ongoing and planned pulsar-timing measurements provide the greatest opportunities for measurements of compact-object masses and new general-relativistic variations in orbits.

Our understanding of the neutron star population is informed to a great degree by large-scale surveys that have been carried out by radio facilities during the past fifty years. We summarize some of the recent breakthroughs in our understanding of the radio pulsar population and look ahead to future yields from upcoming experiments planned for the next decade. By the end of the 2020s, we anticipate a much more complete census of the Galactic population and being able to probe the populations of radio-emitting neutron stars more effectively in external galaxies. Among the anticipated discoveries are pulsar--black hole binary systems that will provide further tests of strong-field gravity, as well as large numbers of millisecond pulsars that are crucial to enhancing the sensitivity of timing arrays for low-frequency gravitational waves.

We present timing solutions for four pulsars discovered in the Green Bank Northern Celestial Cap (GBNCC) survey. All four pulsars are isolated with spin periods between 0.26$\,$s and 1.84$\,$s. PSR J0038$-$2501 has a 0.26$\,$s period and a period derivative of ${7.6} \times {10}^{-19}\,{\rm s\,s}^{-1}$, which is unusually low for isolated pulsars with similar periods. This low period derivative may be simply an extreme value for an isolated pulsar or it could indicate an unusual evolution path for PSR J0038$-$2501, such as a disrupted recycled pulsar (DRP) from a binary system or an orphaned central compact object (CCO). Correcting the observed spin-down rate for the Shklovskii effect suggests that this pulsar may have an unusually low space velocity, which is consistent with expectations for DRPs. There is no X-ray emission detected from PSR J0038$-$2501 in an archival swift observation, which suggests that it is not a young orphaned CCO. The high dispersion measure of PSR J1949+3426 suggests a distance of 12.3$\,$kpc. This distance indicates that PSR J1949+3426 is among the most distant 7% of Galactic field pulsars, and is one of the most luminous pulsars.

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.

Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio flashes likely arriving from far outside the Milky Way galaxy. This phenomenon was discovered at radio frequencies near 1.4 GHz and to date has been observed in one case at as high as 8 GHz, but not below 700 MHz in spite of significant searches at low frequencies. Here we report detections of FRBs at radio frequencies as low as 400 MHz, on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) using the CHIME/FRB instrument. We present 13 FRBs detected during a telescope pre-commissioning phase, when our sensitivity and field-of-view were not yet at design specifications. Emission in multiple events is seen down to 400 MHz, the lowest radio frequency to which we are sensitive. The FRBs show a variety of temporal scattering behaviours, with the majority significantly scattered, and some apparently unscattered to within measurement uncertainty even at our lowest frequencies. Of the 13 reported here, one event has the lowest dispersion measure yet reported, implying it is among the closest yet known, and another has shown multiple repeat bursts, as described in a companion paper. Our low-scattering events suggest that efforts to detect FRBs at radio frequencies below 400 MHz will eventually be successful. The overall scattering properties of our sample suggest that FRBs as a class are preferentially located in environments that scatter radio waves more strongly than the diffuse interstellar medium (ISM) in the Milky Way.

The discovery of a repeating Fast Radio Burst (FRB) source, FRB 121102, eliminated models involving cataclysmic events for this source. No other repeating FRB has yet been detected in spite of many recent FRB discoveries and follow-ups, suggesting repeaters may be rare in the FRB population. Here we report the detection of six repeat bursts from FRB 180814.J0422+73, one of the 13 FRBs detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project during its pre-commissioning phase in July and August 2018. These repeat bursts are consistent with originating from a single position on the sky, with the same dispersion measure (DM), ~189 pc cm-3. This DM is approximately twice the expected Milky Way column density, and implies an upper limit on the source redshift of 0.1, at least a factor of ~2 closer than FRB 121102. In some of the repeat bursts, we observe sub-pulse frequency structure, drifting, and spectral variation reminiscent of that seen in FRB 121102, suggesting similar emission mechanisms and/or propagation effects. This second repeater, found among the first few CHIME/FRB discoveries, suggests that there exists -- and that CHIME/FRB and other wide-field, sensitive radio telescopes will find -- a substantial population of repeating FRBs.

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.

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 provide timing solutions for 45 radio pulsars discovered by the Robert C. Byrd Green Bank Telescope. These pulsars were found in the Green Bank North Celestial Cap pulsar survey, an all-GBT-sky survey being carried out at a frequency of 350 MHz. We include pulsar timing data from the Green Bank Telescope and Low Frequency Array. Our sample includes five fully recycled millisecond pulsars (MSPs, three of which are in a binary system), a new relativistic double neutron star system, an intermediate mass binary pulsar, a mode-changing pulsar, a 138-ms pulsar with a very low magnetic field, and several nulling pulsars. We have measured two post-Keplerian parameters and thus the masses of both objects in the double neutron star system. We also report a tentative companion mass measurement via Shapiro delay in a binary MSP. Two of the MSPs can be timed with high precision and have been included in pulsar timing arrays being used to search for low-frequency gravitational waves, while a third MSP is a member of the black widow class of binaries. Proper motion is measurable in five pulsars and we provide an estimate of their space velocity. We report on an optical counterpart to a new black widow system and provide constraints on the optical counterparts to other binary MSPs. We also present a preliminary analysis of nulling pulsars in our sample. These results demonstrate the scientific return of long timing campaigns on pulsars of all types.

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.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400-800-MHz band. CHIME is comprised of four 20-m x 100-m semi-cylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a >200 square degree field-of-view. This, combined with wide bandwidth, high sensitivity, and a powerful correlator makes CHIME an excellent instrument for the detection of Fast Radio Bursts (FRBs). The CHIME Fast Radio Burst Project (CHIME/FRB) will search beam-formed, high time-and frequency-resolution data in real time for FRBs in the CHIME field-of-view. Here we describe the CHIME/FRB backend, including the real-time FRB search and detection software pipeline as well as the planned offline analyses. We estimate a CHIME/FRB detection rate of 2-42 FRBs/sky/day normalizing to the rate estimated at 1.4-GHz by Vander Wiel et al. (2016). Likely science outcomes of CHIME/FRB are also discussed. CHIME/FRB is currently operational in a commissioning phase, with science operations expected to commence in the latter half of 2018.

We present timing solutions for ten pulsars discovered in 350 MHz searches with the Green Bank Telescope. Nine of these were discovered in the Green Bank Northern Celestial Cap survey and one was discovered by students in the Pulsar Search Collaboratory program in analysis of drift-scan data. Following discovery and confirmation with the Green Bank Telescope, timing has yielded phase-connected solutions with high precision measurements of rotational and astrometric parameters. Eight of the pulsars are slow and isolated, including PSR J0930$-$2301, a pulsar with nulling fraction lower limit of $\sim$30\% and nulling timescale of seconds to minutes. This pulsar also shows evidence of mode changing. The remaining two pulsars have undergone recycling, accreting material from binary companions, resulting in higher spin frequencies. PSR J0557$-$2948 is an isolated, 44 \rm{ms} pulsar that has been partially recycled and is likely a former member of a binary system which was disrupted by a second supernova. The paucity of such so-called `disrupted binary pulsars' (DRPs) compared to double neutron star (DNS) binaries can be used to test current evolutionary scenarios, especially the kicks imparted on the neutron stars in the second supernova. There is some evidence that DRPs have larger space velocities, which could explain their small numbers. PSR J1806+2819 is a 15 \rm{ms} pulsar in a 44 day orbit with a low mass white dwarf companion. We did not detect the companion in archival optical data, indicating that it must be older than 1200 Myr.

Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.

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.

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

We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make fourteen significant measurements of Shapiro delay, including new detections in four pulsar-binary systems (PSRs J0613$-$0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-binary systems. We find a wide range of binary pulsar masses, with values as low as $m_{\rm p} = 1.18^{+0.10}_{-0.09}\text{M}_{\odot}$ for PSR J1918$-$0642 and as high as $m_{\rm p} = 1.928^{+0.017}_{-0.017}\text{M}_{\odot}$ for PSR J1614$-$2230 (both 68.3\% credibility). We make an improved measurement of the Shapiro timing delay in the PSR J1918$-$0642 and J2043+1711 systems, measuring the pulsar mass in the latter system to be $m_{\rm p} = 1.41^{+0.21}_{-0.18}\text{M}_{\odot}$ (68.3\% credibility) for the first time. We measure secular variations of one or more orbital elements in many systems, and use these measurements to further constrain our estimates of the pulsar and companion masses whenever possible. In particular, we used the observed Shapiro delay and periastron advance due to relativistic gravity in the PSR J1903+0327 system to derive a pulsar mass of $m_{\rm p} = 1.65^{+0.02}_{-0.02}\text{M}_{\odot}$ (68.3\% credibility). We discuss the implications that our mass measurements have on the overall neutron-star mass distribution, and on the "mass/orbital-period" correlation due to extended mass transfer.

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

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

Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance. We compare these measurements to distances predicted by the NE2001 interstellar electron density model and find them to be in general agreement. We use measured orbital-decay rates and spin-down rates to confirm two of the parallax distances and to place distance upper limits on other sources; these distance limits agree with the parallax distances with one exception, PSR J1024-0719, which we discuss at length. Using the proper motions of the 37 NANOGrav pulsars in combination with other published measurements, we calculate the velocity dispersion of the millisecond pulsar population in Galactocentric coordinates. We find the radial, azimuthal, and perpendicular dispersions to be 46, 40, and 24 km s-1, respectively, in a model that allows for high-velocity outliers; or 81, 58, and 62 km s-1 for the full population. These velocity dispersions are far smaller than those of the canonical pulsar population, and are similar to older Galactic disk populations. This suggests that millisecond pulsar velocities are largely attributable to their being an old population rather than being artifacts of their birth and evolution as neutron star binary systems. The components of these velocity dispersions follow similar proportions to other Galactic populations, suggesting that our results are not biased by selection effects.

We report 21-yr timing of one of the most precise pulsars: PSR J1713+0747. Its pulse times of arrival are well modeled by a comprehensive pulsar binary model including its three-dimensional orbit and a noise model that incorporates correlated noise such as jitter and red noise. Its timing residuals have weighted root mean square $\sim 92$ ns. The new dataset allows us to update and improve previous measurements of the system properties, including the masses of the neutron star ($1.31\pm0.11$ $M_{\odot}$) and the companion white dwarf ($0.286\pm0.012$ $M_{\odot}$) and the parallax distance $1.15\pm0.03$ kpc. We measured the intrinsic change in orbital period, $\dot{P}^{\rm Int}_{\rm b}$, is $-0.20\pm0.17$ ps s$^{-1}$, which is not distinguishable from zero. This result, combined with the measured $\dot{P}^{\rm Int}_{\rm b}$ of other pulsars, can place a generic limit on potential changes in the gravitational constant $G$. We found that $\dot{G}/G$ is consistent with zero [$(-0.6\pm1.1)\times10^{-12}$ yr$^{-1}$, 95\% confidence] and changes at least a factor of $31$ (99.7\% confidence) more slowly than the average expansion rate of the Universe. This is the best $\dot{G}/G$ limit from pulsar binary systems. The $\dot{P}^{\rm Int}_{\rm b}$ of pulsar binaries can also place limits on the putative coupling constant for dipole gravitational radiation $\kappa_D=(-0.9\pm3.3)\times10^{-4}$ (95\% confidence). Finally, the nearly circular orbit of this pulsar binary allows us to constrain statistically the strong-field post-Newtonian parameters $\Delta$, which describes the violation of strong equivalence principle, and $\hat{\alpha}_3$, which describes a breaking of both Lorentz invariance in gravitation and conservation of momentum. We found, at 95\% confidence, $\Delta<0.01$ and $\hat{\alpha}_3<2\times10^{-20}$ based on PSR J1713+0747.

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.

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

We present updated analyses of pulse profiles and their arrival-times from PSR B1534+12, a 37.9-ms radio pulsar in orbit with another neutron star. A high-precision timing model is derived from twenty-two years of timing data, and accounts for all astrophysical processes that systematically affect pulse arrival-times. Five "post-Keplerian" parameters are measured that represent relativistic corrections to the standard Keplerian quantities of the pulsar's binary orbit. These relativistic parameters are then used to test general relativity by comparing the measurements with their predicted values. We conclude that relativity theory is confirmed to within 0.17% of its predictions. Furthermore, we derive the following astrophysical results from our timing analysis: a distance of $d_{\rm GR} = 1.051 \pm 0.005$ kpc to the pulsar-binary system, by relating the "excess" orbital decay to Galactic parameters; evidence for pulse "jitter" in PSR B1534+12 due to short-term magnetospheric activity; and evolution in pulse-dispersion properties. As a secondary study, we also present several analyses on pulse-structure evolution and its connection to relativistic precession of the pulsar's spin axis. The precession-rate measurement yields a value of $\Omega_1^{\rm spin}$ = $0.59^{+0.12}_{-0.08}$ $^{\circ}$/year (68% confidence) that is consistent with expectations, and represents an additional test of relativistic gravity.

GRB 090417B was an unusually long burst with a T_90 duration of at least 2130 s and a multi-peaked light curve at energies of 15-150 keV. It was optically dark and has been associated with a bright star-forming galaxy at a redshift of 0.345 that is broadly similar to the Milky Way. This is one of the few cases where a host galaxy has been clearly identified for a dark gamma-ray burst and thus an ideal candidate for studying the origin of dark bursts. We find that the dark nature of GRB 090417B cannot be explained by high redshift, incomplete observations, or unusual physics in the production of the afterglow. Assuming the standard relativistic fireball model for the afterglow we find that the optical flux is at least 2.5 mag fainter than predicted by the X-ray flux. The Swift/XRT X -ray data are consistent with the afterglow being obscured by a dense, localized sheet of dust approximately 30-80 pc from the burst along the line of sight. Our results suggest that this dust sheet imparts an extinction of A_V >~ 12 mag, which is sufficient to explain the missing optical flux. GRB 090417B is an example of a gamma-ray burst that is dark due to the localized dust structure in its host galaxy.