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

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 p = 1.18 +0.10 −0.09 M for PSR J1918−0642 and as high as m p = 1.928 +0.017 −0.017 M 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 p = 1.41 +0.21 −0.18 M (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 p = 1.65 +0.02 −0.02 M (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.

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr −1 ) for a GWB from supermassive black hole binaries of <´-A 1.5 10 gw 15 . We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. 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. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, W <´-f h 4.2 10 gw 2 1 0 ( ) . Our limit on the cosmic string GWB, W <´-f h 2.2 10 gw 2 1 0 ( ) , translates to a conservative limit on the cosmic string tension with m <´-G 3.3 10 8 , a factor of four better than the joint Planck and high-lcosmic microwave background data from other experiments.