The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentially explained by a progenitor system that is asymmetric in terms of the stellar component masses, i.e. with a mass ratio $q$ of 0.7-0.8, the known Galactic population of merging double neutron star (DNS) systems (i.e. those that will coalesce within billions of years or less) has, until now, only consisted of nearly equal-mass ($q > 0.9$) binaries. PSR J1913+1102 is a DNS system in a 5-hour, low-eccentricity ($e = 0.09$) orbit, implying an orbital separation of 1.8 solar radii, with the two neutron stars predicted to coalesce in 470 million years due to gravitational-wave emission. Here we report that the masses of the two neutron stars, as measured by a dedicated pulsar timing campaign, are $1.62 \pm 0.03$ and $1.27 \pm 0.03$ solar masses for the pulsar and companion neutron star, respectively; with a measured mass ratio $q = 0.78 \pm 0.03$, it is the most asymmetric DNS among known merging systems. Based on this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30% (90% confidence) of the total population of merging DNS binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

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

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

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 on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo PALFA survey. Phase-coherent timing solutions derived from 2.5 to 5 years of observations carried out at Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65-ms pulsar in a 39.6-day orbit for which we detect $\gamma$-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mass-function system with an orbital period of 3.3 hours that belongs to the non-eclipsing black widow population. We also present PSR J1932+1756, the longest-orbital-period (41.5 days) intermediate-mass binary pulsar known to date. In light of the numerous discoveries of binary MSPs over the past years, we characterize the Galactic distribution of known MSP binaries in terms of binary class. Our results support and strengthen previous claims that the scatter in the Galactic scale height distribution correlates inversely with the binary mass function. We provide evidence of observational biases against detecting the most recycled pulsars near the Galactic plane, which overestimates the scale height of lighter systems. A possible bimodality in the mass function of MSPs with massive white dwarfs is also reported.

In this paper, we present the results of timing observations of PSRs J1949+3106 and J1950+2414, two binary millisecond pulsars discovered in data from the Arecibo ALFA pulsar survey (PALFA). The timing parameters include precise measurements of the proper motions of both pulsars, which show that PSR J1949+3106 has a transversal motion very similar to that of an object in the local standard of rest. The timing also includes measurements of the Shapiro delay and the rate of advance of periastron for both systems. Assuming general relativity, these allow estimates of the masses of the components of the two systems; for PSR J1949+3106, the pulsar mass is $M_p \, = \, 1.34^{+0.17}_{-0.15} \, M_{\odot}$ and the companion mass $M_c \, = \, 0.81^{+0.06}_{-0.05}\, M_{\odot}$; for PSR J1950+2414 $M_p \, = \, 1.496 \, \pm \, 0.023\, M_{\odot}$ and $M_c \, = \, 0.280^{+0.005}_{-0.004}\, M_{\odot}$ (all values 68.3 % confidence limits). We use these masses and proper motions to investigate the evolutionary history of both systems: PSR J1949+3106 is likely the product of a low-kick supernova; PSR J1950+2414 is a member of a new class of eccentric millisecond pulsar binaries with an unknown formation mechanism. We discuss the proposed hypotheses for the formations of these systems in light of our new mass measurements.

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.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

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.

We present a newly implemented single-pulse pipeline for the PALFA survey to efficiently identify single radio pulses from pulsars, Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs). We have conducted a sensitivity analysis of this new pipeline in which multiple single pulses with a wide range of parameters were injected into PALFA data sets and run through the pipeline. Based on the recovered pulses, we find that for pulse widths $\rm < 5\ ms$ the sensitivity of the PALFA pipeline is at most a factor of $\rm \sim 2$ less sensitive to single pulses than our theoretical predictions. For pulse widths $\rm > 10\ ms$, as the $\rm DM$ decreases, the degradation in sensitivity gets worse and can increase up to a factor of $\rm \sim 4.5$. Using this pipeline, we have thus far discovered 7 pulsars and 2 RRATs and identified 3 candidate RRATs and 1 candidate FRB. The confirmed pulsars and RRATs have DMs ranging from 133 to 386 pc cm$^{-3}$ and flux densities ranging from 20 to 160 mJy. The pulsar periods range from 0.4 to 2.1 s. We report on candidate FRB 141113, which we argue is likely astrophysical and extragalactic, having $\rm DM \simeq 400\ pc~cm^{-3}$, which represents an excess over the Galactic maximum along this line of sight of $\rm \sim$ 100 - 200 pc cm$^{-3}$. We consider implications for the FRB population and show via simulations that if FRB 141113 is real and extragalactic, the slope $\alpha$ of the distribution of integral source counts as a function of flux density ($N (>S) \propto S^{-\alpha}$) is $1.4 \pm 0.5$ (95% confidence range). However this conclusion is dependent on several assumptions that require verification.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

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

The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fourier analyses (typically $\sim$16), or due to the strong effect of red noise at low modulation frequencies. We address this reduction in sensitivity by using a time-domain search technique: the Fast-Folding Algorithm (FFA). We designed a program that implements a FFA-based search in the PALFA processing pipeline, and tested the efficiency of the algorithm by performing tests under both ideal, white noise conditions, as well as with real PALFA observational data. In the two scenarios, we show that the time-domain algorithm has the ability to outperform the FFT-based periodicity search implemented in the survey. We perform simulations to compare the previously reported PALFA sensitivity with that obtained using our new FFA implementation. These simulations show that for a pulsar having a pulse duty cycle of roughly 3%, the performance of our FFA pipeline exceeds that of our FFT pipeline for pulses with DM $\lesssim$ 40 pc cm$^{-3}$ and for periods as short as $\sim$500 ms, and that the survey sensitivity is improved by at least a factor of two for periods $\gtrsim$ 6 sec. Discoveries from the implementation of the algorithm in PALFA are also presented in this paper.

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.

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 report the discovery and initial follow-up of a double neutron star (DNS) system, PSR J1946$+$2052, with the Arecibo L-Band Feed Array pulsar (PALFA) survey. PSR J1946$+$2052 is a 17-ms pulsar in a 1.88-hour, eccentric ($e \, =\, 0.06$) orbit with a $\gtrsim 1.2 \, M_\odot$ companion. We have used the Jansky Very Large Array to localize PSR J1946$+$2052 to a precision of 0.09 arcseconds using a new phase binning mode. We have searched multiwavelength catalogs for coincident sources but did not find any counterparts. The improved position enabled a measurement of the spin period derivative of the pulsar ($\dot{P} \, = \, 9\,\pm \, 2 \,\times 10^{-19}$); the small inferred magnetic field strength at the surface ($B_S \, = \, 4 \, \times \, 10^9 \, \rm G$) indicates that this pulsar has been recycled. This and the orbital eccentricity lead to the conclusion that PSR J1946$+$2052 is in a DNS system. Among all known radio pulsars in DNS systems, PSR J1946$+$2052 has the shortest orbital period and the shortest estimated merger timescale, 46 Myr; at that time it will display the largest spin effects on gravitational wave waveforms of any such system discovered to date. We have measured the advance of periastron passage for this system, $\dot{\omega} \, = \, 25.6 \, \pm \, 0.3\, \deg \rm yr^{-1}$, implying a total system mass of only 2.50 $\pm$ 0.04 $M_\odot$, so it is among the lowest mass DNS systems. This total mass measurement combined with the minimum companion mass constrains the pulsar mass to $\lesssim 1.3 \, M_\odot$.

We present a timing and glitch analysis of the young X-ray pulsar PSR J0537$-$6910, located within the Large Magellanic Cloud, using 13 years of data from the now decommissioned Rossi X-ray Timing Explorer. Rotating with a spin period of 16 ms, PSR J0537$-$6910 is the fastest spinning and most energetic young pulsar known. It also displays the highest glitch activity of any known pulsar. We have found 42 glitches over the data span, corresponding to a glitch rate of 3.2 yr$^{-1}$, with an overall glitch activity rate of $8.8\times 10^{-7}\,$yr$^{-1}$. The high glitch frequency has allowed us to study the glitch behavior in ways that are inaccessible in other pulsars. We observe a strong linear correlation between spin frequency glitch magnitude and wait time to the following glitch. We also find that the post-glitch spin-down recovery is well described by a single two-component model fit to all glitches for which we have adequate input data. This consists of an exponential amplitude $A = (7.6 \pm 1.0)\times 10^{-14}\,$s$^{-2}$ and decay timescale $\tau = 27_{-6}^{+7}\,$d, and linear slope $m = (4.1\pm 0.4)\times 10^{-16}\,$s$^{-2}\,$d$^{-1}$. The latter slope corresponds to a second frequency derivative $\ddot{\nu} = (4.7\pm 0.5) \times 10^{-22}\,$s$^{-3}$, from which we find an implied braking index $n=7.4 \pm 0.8$. We also present a maximum-likelihood technique for searching for periods in event-time data, which we used to both confirm previously published values and determine rotation frequencies in later observations. We discuss the implied constraints on glitch models from the observed behavior of this system, which we argue cannot be fully explained in the context of existing theories.

We report on the discovery and timing observations of 29 distant long-period pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, confirmation and timing observations of these pulsars over several years at Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation and radiation properties. We have used multi-frequency data to measure the interstellar scattering properties of some of these pulsars. Most of the pulsars have properties that mirror those of the previously known pulsar population, although four show some notable characteristics. PSRs J1907+0631 and J1925+1720 are young and are associated with supernova remnants or plerionic nebulae: J1907+0631 lies close to the center of SNR G40.5-0.5, while J1925+1720 is coincident with a high-energy Fermi gamma-ray source. One pulsar, J1932+1500, is in a surprisingly eccentric, 199-day binary orbit with a companion having a minimum mass of 0.33 solar masses. Several of the sources exhibit timing noise, and two, PSRs J0611+1436 and J1907+0631, have both suffered large glitches, but with very different post-glitch rotation properties. In particular, the rotational period of PSR J0611+1436 will not recover to its pre-glitch value for about 12 years, a far greater recovery timescale than seen following any other large glitches.

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.

We report the discovery of two long-term intermittent radio pulsars in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, extended observations of these pulsars over several years at Jodrell Bank Observatory have revealed the details of their rotation and radiation properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bi-modal intermittency, switching between active (ON) and inactive (OFF) emission states and indicating the presence of a large, hitherto unrecognised, underlying population of such objects. For PSR J1929+1357, the initial duty cycle was fON=0.008, but two years later this changed, quite abruptly, to fON=0.16. This is the first time that a significant evolution in the activity of an intermittent pulsar has been seen and we show that the spin-down rate of the pulsar is proportional to the activity. The spin-down rate of PSR J1929+1357 is increased by a factor of 1.8 when it is in active mode, similar to the increase seen in the other three known long-term intermittent pulsars.

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.

We report here the Einstein@Home discovery of PSR J1913+1102, a 27.3-ms pulsar found in data from the ongoing Arecibo PALFA pulsar survey. The pulsar is in a 4.95-hr double neutron star (DNS) system with an eccentricity of 0.089. From radio timing with the Arecibo 305-m telescope, we measure the rate of advance of periastron to be 5.632(18) deg/yr. Assuming general relativity accurately models the orbital motion, this corresponds to a total system mass of 2.875(14) solar masses, similar to the mass of the most massive DNS known to date, B1913+16, but with a much smaller eccentricity. The small eccentricity indicates that the second-formed neutron star (the companion of PSR J1913+1102) was born in a supernova with a very small associated kick and mass loss. In that case this companion is likely, by analogy with other systems, to be a light (1.2 solar mass) neutron star; the system would then be highly asymmetric. A search for radio pulsations from the companion yielded no plausible detections, so we can't yet confirm this mass asymmetry. By the end of 2016, timing observations should permit the detection of two additional post-Keplerian parameters: the Einstein delay, which will enable precise mass measurements and a verification of the possible mass asymmetry of the system, and the orbital decay due to the emission of gravitational waves, which will allow another test of the radiative properties of gravity. The latter effect will cause the system to coalesce in ~0.5 Gyr.

We report the discovery and timing results for five millisecond pulsars (MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726, J1938+2516, and J1957+2516. Timing observations of the 5 pulsars were conducted with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3 yr. All of the MSPs except one (PSR J1914+0659) are in binary systems with low eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ~0.1 $M_\odot$ companion and possible eclipses that last ~10% of the orbit. The position of PSR J1957+2516 is also coincident with a NIR source. All 5 MSPs are distant (>3.1 kpc) as determined from their dispersion measures, and none of them show evidence of $\gamma$-ray pulsations in a search of Fermi Gamma-Ray Space Telescope data. These 5 MSPs bring the total number of MSPs discovered by the PALFA survey to 26 and further demonstrate the power of this survey in finding distant, highly dispersed MSPs deep in the Galactic plane.

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.

Fast Radio Bursts are millisecond-duration astronomical radio pulses of unknown physical origin that appear to come from extragalactic distances. Previous follow-up observations have failed to find additional bursts at the same dispersion measures (i.e. integrated column density of free electrons between source and telescope) and sky position as the original detections. The apparent non-repeating nature of the fast radio bursts has led several authors to hypothesise that they originate in cataclysmic astrophysical events. Here we report the detection of ten additional bursts from the direction of FRB121102, using the 305-m Arecibo telescope. These new bursts have dispersion measures and sky positions consistent with the original burst. This unambiguously identifies FRB121102 as repeating and demonstrates that its source survives the energetic events that cause the bursts. Additionally, the bursts from FRB121102 show a wide range of spectral shapes that appear to be predominantly intrinsic to the source and which vary on timescales of minutes or shorter. While there may be multiple physical origins for the population of fast radio bursts, the repeat bursts with high dispersion measure and variable spectra specifically seen from FRB121102 support models that propose an origin in a young, highly magnetised, extragalactic neutron star.

We present a phase-coherent timing solution for PSR J1640-4631, a young 206 ms pulsar using X-ray timing observations taken with NuSTAR. Over this timing campaign, we have measured the braking index of PSR J1640-4631 to be n = 3.15+/-0.03. Using a series of simulations, we argue that this unusually high braking index is not due to timing noise, but is intrinsic to the pulsar's spin-down. We cannot, however, rule out contamination due to an unseen glitch recovery, although the recovery timescale would have to be longer than most yet observed. If this braking index is eventually proven to be stable, it demonstrates that pulsar braking indices greater than 3 are allowed in nature, hence other physical mechanisms such as mass or magnetic quadrupoles are important in pulsar spin-down. We also present a 3-sigma upper limit on the pulsed flux at 1.4 GHz of 0.018 mJy.

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 present timing solutions and spin properties of the young pulsar PSR B0540-69 from analysis of 15.8 yr of data from the Rossi X-Ray Timing Explorer. We perform a partially phase-coherent timing analysis in order to mitigate the pronounced effects of timing noise in this pulsar. We also perform fully coherent timing over large subsets of the data set in order to arrive at a more precise solution. In addition to the previously reported first glitch undergone by this pulsar, we find a second glitch, which occurred at MJD 52927 $\pm$ 4, with fractional changes in spin frequency $\Delta\nu/\nu = (1.64 \pm 0.05) \times 10^{-9}$ and spin-down rate $\Delta\dot{\nu}/\dot{\nu} = (0.930 \pm 0.011) \times 10^{-4}$ (taken from our fully coherent analysis). We measure a braking index that is consistent over the entire data span, with a mean value $n = 2.129 \pm 0.012$, from our partially coherent timing analysis. We also investigated the emission behavior of this pulsar, and have found no evidence for significant flux changes, flares, burst-type activity, or pulse profile shape variations. While there is strong evidence for the much-touted similarity of PSR B0540-69 to the Crab pulsar, they nevertheless differ in several aspects, including glitch activity, where PSR B0540-69 can be said to resemble certain other very young pulsars. It seems clear that the specific processes governing the formation, evolution, and interiors of this population of recently born neutron stars can vary significantly, as reflected in their observed properties.

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.

We report the discovery of the millisecond pulsar (MSP) PSR J1950+2414 ($P=4.3$ ms) in a binary system with an eccentric ($e=0.08$) 22-day orbit in Pulsar ALFA survey observations with the Arecibo telescope. Its companion star has a median mass of 0.3 $M_\odot$ and is most likely a white dwarf. Fully recycled MSPs like this one are thought to be old neutron stars spun-up by mass transfer from a companion star. This process should circularize the orbit, as is observed for the vast majority of binary MSPs, which predominantly have orbital eccentricities $e < 0.001$. However, four recently discovered binary MSPs have orbits with $0.027 < e < 0.44$; PSR J1950+2414 is the fifth such system to be discovered. The upper limits for its intrinsic spin period derivative and inferred surface magnetic field strength are comparable to those of the general MSP population. The large eccentricities are incompatible with the predictions of the standard recycling scenario: something unusual happened during their evolution. Proposed scenarios are a) initial evolution of the pulsar in a triple system which became dynamically unstable, b) origin in an exchange encounter in an environment with high stellar density, c) rotationally delayed accretion-induced collapse of a super-Chandrasekhar white dwarf, and d) dynamical interaction of the binary with a circumbinary disk. We compare the properties of all five known eccentric MSPs with the predictions of these formation channels. Future measurements of the masses and proper motion might allow us to firmly exclude some of the proposed formation scenarios.

The on-going PALFA survey at the Arecibo Observatory began in 2004 and is searching for radio pulsars in the Galactic plane at 1.4 GHz. Observations since 2009 have been made with new wider-bandwidth spectrometers than were previously employed in this survey. A new data reduction pipeline has been in place since mid-2011 which consists of standard methods using dedispersion, searches for accelerated periodic sources, and search for single pulses, as well as new interference-excision strategies and candidate selection heuristics. This pipeline has been used to discover 41 pulsars, including 8 millisecond pulsars (MSPs; P < 10 ms), bringing the PALFA survey's discovery totals to 145 pulsars, including 17 MSPs, and one Fast Radio Burst (FRB). The pipeline presented here has also re-detected 188 previously known pulsars including 60 found in PALFA data by re-analyzing observations previously searched by other pipelines. A comprehensive description of the survey sensitivity, including the effect of interference and red noise, has been determined using synthetic pulsar signals with various parameters and amplitudes injected into real survey observations and subsequently recovered with the data reduction pipeline. We have confirmed that the PALFA survey achieves the sensitivity to MSPs predicted by theoretical models. However, we also find that compared to theoretical survey sensitivity models commonly used there is a degradation in sensitivity to pulsars with periods P >= 100 ms that gradually becomes up to a factor of ~10 worse for P > 4 s at DM < 150 pc/cc. This degradation of sensitivity at long periods is largely due to red noise. We find that 35 +- 3% of pulsars are missed despite being bright enough to be detected in the absence of red noise. This reduced sensitivity could have implications on estimates of the number of long-period pulsars in the Galaxy.

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 the discovery of five millisecond pulsars (MSPs) from the PALFA Galactic plane survey using Arecibo. Four of these (PSRs J0557+1551, J1850+0244, J1902+0300, and J1943+2210) are binary pulsars whose companions are likely white dwarfs, and one (PSR J1905+0453) is isolated. Phase-coherent timing solutions, ranging from $\sim$1 to $\sim$3 years in length, and based on observations from the Jodrell Bank and Arecibo telescopes, provide precise determinations of spin, orbital, and astrometric parameters. All five pulsars have large dispersion measures ($>100$ pc cm$^{-3}$, within the top 20% of all known Galactic field MSPs) and are faint (1.4 GHz flux density < 0.1 mJy, within the faintest 5% of all known Galactic field MSPs), illustrating PALFA's ability to find increasingly faint, distant MSPs in the Galactic plane. In particular, PSR J1850+0244 has a dispersion measure of 540 pc cm$^{-3}$, the highest of all known MSPs. Such distant, faint MSPs are important input for accurately modeling the total Galactic MSP population.

The pulsar PSR J1756$-$2251 resides in a relativistic double neutron star (DNS) binary system with a 7.67-hr orbit. We have conducted long-term precision timing on more than 9 years of data acquired from five telescopes, measuring five post-Keplerian parameters. This has led to several independent tests of general relativity (GR), the most constraining of which shows agreement with the prediction of GR at the 4% level. Our measurement of the orbital decay rate disagrees with that predicted by GR, likely due to systematic observational biases. We have derived the pulsar distance from parallax and orbital decay measurements to be 0.73$_{-0.24}^{+0.60}$ kpc (68%) and < 1.2 kpc (95% upper limit), respectively; these are significantly discrepant from the distance estimated using Galactic electron density models. We have found the pulsar mass to be 1.341$\pm$0.007 M$_\odot$, and a low neutron star (NS) companion mass of 1.230$\pm$0.007 M$_\odot$. We also determined an upper limit to the spin-orbit misalignment angle of 34{\deg} (95%) based on a system geometry fit to long-term profile width measurements. These and other observed properties have led us to hypothesize an evolution involving a low mass loss, symmetric supernova progenitor to the second-formed NS companion, as is thought to be the case for the double pulsar system PSR J0737$-$3039A/B. This would make PSR J1756$-$2251 the second compact binary system providing concrete evidence for this type of NS formation channel.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

The Double Pulsar, PSR J$0737$$-$$3039$A/B, is a unique system in which both neutron stars have been detected as radio pulsars. As shown in Ferdman et al., there is no evidence for pulse profile evolution of the A pulsar, and the geometry of the pulsar was fit well with a double-pole circular radio beam model. Assuming a more realistic polar cap model with a vacuum retarded dipole magnetosphere configuration including special relativistic effects, we create synthesized pulse profiles for A given the best-fit geometry from the simple circular beam model. By fitting synthesized pulse profiles to those observed from pulsar A, we constrain the geometry of the radio beam, namely the half-opening angle and the emission altitude, to be $30^\circ$ and $10$ neutron star radii, respectively. Combining the observational constraints of PSR J$0737$$-$$3039$A/B, we are able to construct the full three-dimensional orbital geometry of the Double Pulsar. The relative angle between the spin axes of the two pulsars ($\Delta_S$) is estimated to be ($138^\circ \pm 5^\circ$) at the current epoch and will likely remain constant until tidal interactions become important in $\sim$$85$ Myr, at merger.

We report here a series of observations of the interstellar scintillation (ISS) of the double pulsar J0737$-$3039 over the course of 18 months. As in earlier work (Coles et al., 2005) the basic phenomenon is the variation in the ISS caused by the changing transverse velocities of each pulsar, the ionized interstellar medium (IISM), and the Earth. The transverse velocity of the binary system can be determined both by VLBI and timing observations. The orbital velocity and inclination is almost completely determined from timing observations, but the direction of the orbital angular momentum is not known. Since the Earth's velocity is known, and can be compared with the orbital velocity by its effect on the timescale of the ISS, we can determine the orientation $\Omega$ of the pulsar orbit with respect to equatorial coordinates ($\Omega = 65\pm2$ deg). We also resolve the ambiguity ($i= 88.7$ or $91.3$ deg) in the inclination of the orbit deduced from the measured Shapiro delay by our estimate $i=88.1\pm0.5$ deg. This relies on analysis of the ISS over both frequency and time and provides a model for the location, anisotropy, turbulence level and transverse phase gradient of the IISM. We find that the IISM can be well-modeled during each observation, typically of a few orbital periods, but its turbulence level and mean velocity vary significantly over the 18 months.

In this paper, we present a novel artificial intelligence (AI) program that identifies pulsars from recent surveys using image pattern recognition with deep neural nets---the PICS (Pulsar Image-based Classification System) AI. The AI mimics human experts and distinguishes pulsars from noise and interferences by looking for patterns from candidate. The information from each pulsar candidate is synthesized in four diagnostic plots, which consist of up to thousands pixel of image data. The AI takes these data from each candidate as its input and uses thousands of such candidates to train its ~9000 neurons. Different from other pulsar selection programs which use pre-designed patterns, the PICS AI teaches itself the salient features of different pulsars from a set of human-labeled candidates through machine learning. The deep neural networks in this AI system grant it superior ability in recognizing various types of pulsars as well as their harmonic signals. The trained AI's performance has been validated with a large set of candidates different from the training set. In this completely independent test, PICS ranked 264 out of 277 pulsar-related candidates, including all 56 previously known pulsars, to the top 961 (1%) of 90008 test candidates, missing only 13 harmonics. The first non-pulsar candidate appears at rank 187, following 45 pulsars and 141 harmonics. In other words, 100% of the pulsars were ranked in the top 1% of all candidates, while 80% were ranked higher than any noise or interference. The performance of this system can be improved over time as more training data are accumulated. This AI system has been integrated into the PALFA survey pipeline and has discovered six new pulsars to date.

We analyze timing noise from five years of Arecibo and Green Bank observations of the seventeen millisecond pulsars of the North-American Nanohertz Observatory for Gravitational Waves (NANOGrav) pulsar timing array. The weighted autocovariance of the timing residuals was computed for each pulsar and compared against two possible models for the underlying noise process. The first model includes red noise and predicts the autocovariance to be a decaying exponential as a function of time lag. The second model is Gaussian white noise whose autocovariance would be a delta function. We also perform a ``nearest-neighbor" correlation analysis. We find that the exponential process does not accurately describe the data. Two pulsars, J1643-1224 and J1910+1256, exhibit weak red noise, but the rest are well described as white noise. The overall lack of evidence for red noise implies that sensitivity to a (red) gravitational wave background signal is limited by statistical rather than systematic uncertainty. In all pulsars, the ratio of non-white noise to white noise is low, so that we can increase the cadence or integration times of our observations and still expect the root-mean-square of timing residual averages to decrease by the square-root of observation time, which is key to improving the sensitivity of the pulsar timing array.

Modern radio pulsar surveys produce a large volume of prospective candidates, the majority of which are polluted by human-created radio frequency interference or other forms of noise. Typically, large numbers of candidates need to be visually inspected in order to determine if they are real pulsars. This process can be labor intensive. In this paper, we introduce an algorithm called PEACE (Pulsar Evaluation Algorithm for Candidate Extraction) which improves the efficiency of identifying pulsar signals. The algorithm ranks the candidates based on a score function. Unlike popular machine-learning based algorithms, no prior training data sets are required. This algorithm has been applied to data from several large-scale radio pulsar surveys. Using the human-based ranking results generated by students in the Arecibo Remote Command enter programme, the statistical performance of PEACE was evaluated. It was found that PEACE ranked 68% of the student-identified pulsars within the top 0.17% of sorted candidates, 95% within the top 0.34%, and 100% within the top 3.7%. This clearly demonstrates that PEACE significantly increases the pulsar identification rate by a factor of about 50 to 1000. To date, PEACE has been directly responsible for the discovery of 47 new pulsars, 5 of which are millisecond pulsars that may be useful for pulsar timing based gravitational-wave detection projects.

We report the Fermi Large Area Telescope discovery of gamma-ray pulsations from the 22.7 ms pulsar A in the double pulsar system J0737-3039A/B. This is the first mildly recycled millisecond pulsar (MSP) detected in the GeV domain. The 2.7 s companion object PSR J0737-3039B is not detected in gamma rays. PSR J0737-3039A is a faint gamma-ray emitter, so that its spectral properties are only weakly constrained; however, its measured efficiency is typical of other MSPs. The two peaks of the gamma-ray light curve are separated by roughly half a rotation and are well offset from the radio and X-ray emission, suggesting that the GeV radiation originates in a distinct part of the magnetosphere from the other types of emission. From the modeling of the radio and the gamma-ray emission profiles and the analysis of radio polarization data, we constrain the magnetic inclination $\alpha$ and the viewing angle $\zeta$ to be close to 90$^\circ$, which is consistent with independent studies of the radio emission from PSR J0737-3039A. A small misalignment angle between the pulsar's spin axis and the system's orbital axis is therefore favored, supporting the hypothesis that pulsar B was formed in a nearly symmetric supernova explosion as has been discussed in the literature already.

The double pulsar system PSR J0737-3039A/B is a double neutron star binary, with a 2.4-hour orbital period, which has allowed measurement of relativistic orbital perturbations to high precision. The low mass of the second-formed neutron star, as well as the low system eccentricity and proper motion, point to a different evolutionary scenario compared to most other known double neutron star systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to be a near-orthogonal rotator, with an average separation between its spin and magnetic axes of 90 +/- 11 +/- 5 deg. Furthermore, we find a mean 95% upper limit on the misalignment between its spin and orbital angular momentum axes of 3.2 deg, assuming that the observed emission comes from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7x10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.

The double pulsar PSR J0737-3039A/B displays short, 30 s eclipses that arise around conjunction when the radio waves emitted by pulsar A are absorbed as they propagate through the magnetosphere of its companion pulsar B. These eclipses offer a unique opportunity to probe directly the magnetospheric structure and the plasma properties of pulsar B. We have performed a comprehensive analysis of the eclipse phenomenology using multi-frequency radio observations obtained with the Green Bank Telescope. We have characterized the periodic flux modulations previously discovered at 820 MHz by McLaughlin et al., and investigated the radio frequency dependence of the duration and depth of the eclipses. Based on their weak radio frequency evolution, we conclude that the plasma in pulsar B's magnetosphere requires a large multiplicity factor (~ 10^5). We also found that, as expected, flux modulations are present at all radio frequencies in which eclipses can be detected. Their complex behavior is consistent with the confinement of the absorbing plasma in the dipolar magnetic field of pulsar B as suggested by Lyutikov & Thompson and such a geometric connection explains that the observed periodicity is harmonically related to pulsar B's spin frequency. We observe that the eclipses require a sharp transition region beyond which the plasma density drops off abruptly. Such a region defines a plasmasphere which would be well inside the magnetospheric boundary of an undisturbed pulsar. It is also two times smaller than the expected standoff radius calculated using the balance of the wind pressure from pulsar A and the nominally estimated magnetic pressure of pulsar B.

We present high-precision timing of five millisecond pulsars (MSPs) carried out for more than seven years; four pulsars are in binary systems and one is isolated. We are able to measure the pulsars' proper motions and derive an estimate for their space velocities. The measured two-dimensional velocities are in the range 70-210 km/s, consistent with those measured for other MSPs. We also use all the available proper motion information for isolated and binary MSPs to update the known velocity distribution for these populations. As found by earlier works, we find that the velocity distribution of binary and isolated MSPs are indistinguishable with the current data. Four of the pulsars in our observing program are highly recycled with low-mass white dwarf companions and we are able to derive accurate binary parameters for these systems. For three of these binary systems we are able to place initial constraints on the pulsar masses with best-fit values in the range 1.0-1.6 M_sun. The implications of the results presented here to our understanding of binary pulsar evolution are discussed. The updated parameters for the binary systems studied here, together with recently discovered similar systems, allowed us to update previous limits on the the violation of the strong equivalence principle through the parameter |Delta| to 4.6x10^-3 (95% confidence) and the violation of Lorentz-invariance/momentum-conservation through the parameter |hat{alpha}_3| to 5.5x10^-20 (95% confidence).

Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^{\alpha}$, where $f$ is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude $A$ as a function of the unknown spectral slope $\alpha$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $\alpha=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.

High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.

We present the evolution of the radio emission from the 2.8-s pulsar of the double pulsar system PSR J0737-3039A/B. We provide an update on the Burgay et al. (2005) analysis by describing the changes in the pulse profile and flux density over five years of observations, culminating in the B pulsar's radio disappearance in 2008 March. Over this time, the flux density decreases by 0.177 mJy/yr at the brightest orbital phases and the pulse profile evolves from a single to a double peak, with a separation rate of 2.6 deg/yr. The pulse profile changes are most likely caused by relativistic spin precession, but can not be easily explained with a circular hollow-cone beam as in the model of Clifton & Weisberg (2008). Relativistic spin precession, coupled with an elliptical beam, can model the pulse profile evolution well. This particular beam shape predicts geometrical parameters for the two bright orbital phases which are consistent and similar to those derived by Breton et al. (2008). However, the observed decrease in flux over time and B's eventual disappearance cannot be easily explained by the model and may be due to the changing influence of A on B.

The European Pulsar Timing Array (EPTA) is a multi-institutional, multi-telescope collaboration, with the goal of using high-precision pulsar timing to directly detect gravitational waves. In this article we discuss the EPTA member telescopes, current achieved timing precision, and near-future goals. We report a preliminary upper limit to the amplitude of a gravitational wave background. We also discuss the Large European Array for Pulsars, in which the five major European telescopes involved in pulsar timing will be combined to provide a coherent array that will give similar sensitivity to the Arecibo radio telescope, and larger sky coverage.

PSR J1802-2124 is a 12.6-ms pulsar in a 16.8-hour binary orbit with a relatively massive white dwarf (WD) companion. These properties make it a member of the intermediate-mass class of binary pulsar (IMBP) systems. We have been timing this pulsar since its discovery in 2002. Concentrated observations at the Green Bank Telescope, augmented with data from the Parkes and Nancay observatories, have allowed us to determine the general relativistic Shapiro delay. This has yielded pulsar and white dwarf mass measurements of 1.24(11) and 0.78(4) solar masses (68% confidence), respectively. The low mass of the pulsar, the high mass of the WD companion, the short orbital period, and the pulsar spin period may be explained by the system having gone through a common-envelope phase in its evolution. We argue that selection effects may contribute to the relatively small number of known IMBPs.

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

Nancay radiotelescope is involved in high precision timing since 20 years. Since 2004, a coherent dedispersion instrumentation enables numerous routine observations on more than 200 pulsars using half of the time if this 100-meters class radiotelescope. Two main programs are currently conducted. A large set of young and old pulsars is timed for a multi-wavelength approach, complementary to the very successful high energy observations of pulsars done by FERMI. A set of highly stable millisecond pulsars is monitored as our contribution to the European Pulsar Timing Array in order to probe the cosmological Gravitational Wave Background.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) is a consortium of astronomers whose goal is the creation of a galactic scale gravitational wave observatory sensitive to gravitational waves in the nHz-microHz band. It is just one component of an international collaboration involving similar organizations of European and Australian astronomers who share the same goal. Gravitational waves, a prediction of Einstein's general theory of relativity, are a phenomenon of dynamical space-time generated by the bulk motion of matter, and the dynamics of space-time itself. They are detectable by the small disturbance they cause in the light travel time between some light source and an observer. NANOGrav exploits radio pulsars as both the light (radio) source and the clock against which the light travel time is measured. In an array of radio pulsars gravitational waves manifest themselves as correlated disturbances in the pulse arrival times. The timing precision of today's best measured pulsars is less than 100 ns. With improved instrumentation and signal-to-noise it is widely believed that the next decade could see a pulsar timing network of 100 pulsars each with better than 100 ns timing precision. Such a pulsar timing array (PTA), observed with a regular cadence of days to weeks, would be capable of observing supermassive black hole binaries following galactic mergers, relic radiation from early universe phenomena such as cosmic strings, cosmic superstrings, or inflation, and more generally providing a vantage on the universe whose revolutionary potential has not been seen in the 400 years since Galileo first turned a telescope to the heavens.

The double pulsar PSR J0737-3039A/B consists of two neutron stars in a highly relativistic orbit that displays a roughly 30-second eclipse when pulsar A passes behind pulsar B. Describing this eclipse of pulsar A as due to absorption occurring in the magnetosphere of pulsar B, we successfully use a simple geometric model to characterize the observed changing eclipse morphology and to measure the relativistic precession of pulsar B's spin axis around the total orbital angular momentum. This provides a test of general relativity and alternative theories of gravity in the strong-field regime. Our measured relativistic spin precession rate of 4.77 (+0.66,-0.65) degrees per year (68% confidence level) is consistent with that predicted by general relativity within an uncertainty of 13%.

The double pulsar system PSR J0737-3039A/B is a highly relativistic double neutron star (DNS) binary, with a 2.4-hour orbital period. The low mass of the second-formed NS, as well the low system eccentricity and proper motion, point to a different evolutionary scenario compared to other known DNS systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to have a low misalignment between its spin and orbital angular momentum axes, with a 68.3% upper limit of 6.1 degrees, assuming emission from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

Low-mass white dwarfs can either be produced in low-mass X-ray binaries by stable mass transfer to a neutron star, or in a common-envelope phase with a heavier white dwarf companion. We have searched 8 low-mass white dwarf candidates recently identified in the Sloan Digital Sky Survey for radio pulsations from pulsar companions, using the Green Bank Telescope at 340MHz. We have found no pulsations down to flux densities of 0.6-0.8 mJy/kpc^2 and conclude that a given low-mass helium-core white dwarf has a probability of < 0.18+-0.05 of being in a binary with a radio pulsar.

The double pulsar system, PSR J0737-3039A/B, is unique in that both neutron stars are detectable as radio pulsars. This, combined with significantly higher mean orbital velocities and accelerations when compared to other binary pulsars, suggested that the system would become the best available testbed for general relativity and alternative theories of gravity in the strong-field regime. Here we report on precision timing observations taken over the 2.5 years since its discovery and present four independent strong-field tests of general relativity. Use of the theory-independent mass ratio of the two stars makes these tests uniquely different from earlier studies. By measuring relativistic corrections to the Keplerian description of the orbital motion, we find that the ``post-Keplerian'' parameter s agrees with the value predicted by Einstein's theory of general relativity within an uncertainty of 0.05%, the most precise test yet obtained. We also show that the transverse velocity of the system's center of mass is extremely small. Combined with the system's location near the Sun, this result suggests that future tests of gravitational theories with the double pulsar will supersede the best current Solar-system tests. It also implies that the second-born pulsar may have formed differently to the usually assumed core-collapse of a helium star.

This first ever double pulsar system consists of two pulsars orbiting the common center of mass in a slightly eccentric orbit of only 2.4-hr duration. The pair of pulsars with pulse periods of 22 ms and 2.8 sec, respectively, confirms the long-proposed recycling theory for millisecond pulsars and provides an exciting opportunity to study the works of pulsar magnetospheres by a very fortunate geometrical alignment of the orbit relative to our line-of-sight. In particular, this binary system represents a truly unique laboratory for relativistic gravitational physics. This contribution serves as an update on the currently obtained results and their consequences for the test of general relativity in the strong-field regime. A complete and more up-to-date report of the timing results will be presented elsewhere shortly.

Time-series data taken with the Hubble Space Telescope of a field six core radii (~5') from the center of the globular cluster Messier 4, and covering a period of about 10 weeks in early 2001, have been analyzed in search of variable objects. Various criteria were employed to select candidate variable stars. Period searches were performed on the selected candidates using phase dispersion minimization (PDM). The reliability of the PDM search results was tested using synthetic light curves of eclipsing binary stars and sinusoidal light curves of different periods. Results from this analysis showed that there are probably no eclipsing binary stars or periodic variables in our field with periods on the order of a few hours to a few days, down to limiting magnitudes of V~25 and I~24, which is consistent with the absence of contact binaries such as W Ursae Majoris systems. However, one candidate variable star does show an increase in brightness of ~0.1 magnitudes in both bandpasses, which seems to last for a few days. Possible explanations concerning the nature of this object include a binary system with a white dwarf primary and a low-mass main sequence secondary, or a BY Draconis variable star. We are able to set an upper limit to the observed fraction of photometric variability in this dataset of 0.05%.