Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, A c,yr , to be < 1.0×10-15 with 95% confidence. This limit excludes predicted ranges for A c,yr from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.Studies of the dynamics of stars and gas in nearby galaxies provide strong evidence for the ubiquity of supermassive (> 10 6 solar mass) black holes (SMBHs) (1). Observations of luminous quasars indicate that SMBHs are hosted by galaxies throughout the history of the universe (2) and affect global properties of the host galaxies (3). The prevailing dark energycold dark matter cosmological paradigm predicts that large galaxies are assembled through the hierarchical merging of smaller galaxies. The remnants of mergers can host gravitationally bound binary SMBHs with orbits decaying through the emission of gravitational waves (GWs) (4).Gravitational waves from binary SMBHs, with periods between ~ 0.1 and 30 yr (5), can be detected or constrained by monitoring, for years to decades, a set of rapidly rotating millisecond pulsars (MSPs) distributed throughout our galaxy. Radio emission beams from MSPs are observed as pulses that can be time-tagged with as small as 20 ns precision (6). When traveling across the pulsar-Earth line of sight, GWs induce variations in the arrival times of the pulses (7).The superposition of GWs from the binary SMBH population is a stochastic background (GWB), which is typically characterized by the strain-amplitude spectrum h c (f)=A c,yr [f/(1 yr -1 )] -2/3 , where f is the GW frequency, A c,yr is the characteristic amplitude of the GWB measured at f = 1 yr -1 , predicted to be A c,yr > 10 -15 (5,(8)(9)(10)(11)(12), and -2/3 is the predicted spectral index (5,(8)(9)(10)(11)(12). The GWB will add low-frequency perturbations to pulse arrival times. While the detection of the GWB would confirm the presence of a cosmological population of binary SMBHs, limits on its amplitude constrain models of galaxy and SMBH evolution (8).As part of the Parkes Pulsar Timing Array project to detect GWs (6), we have been monitoring 24 pulsars with the 64-m Parkes radio telescope. We have produced a new data set, using observations taken at a central wavelength of 10 cm and previously reported methods (6,8), that spans 11 yr, which is 3 yr longer than previous data sets analyzed at this wavelength. In addition to having greater sensitivity to the GWB because of the longer duration, the data set was improved by identifying and correc...

The analysis of high-precision timing observations of an array of ∼20 millisecond pulsars (a so-called 'timing array') may ultimately result in the detection of a stochastic gravitationalwave background. The feasibility of such a detection and the required duration of this type of experiment are determined by the achievable rms of the timing residuals and the timing stability of the pulsars involved. We present results of the first long-term, high-precision timing campaign on a large sample of millisecond pulsars used in gravitational-wave detection projects. We show that the timing residuals of most pulsars in our sample do not contain significant low-frequency noise that could limit the use of these pulsars for decade-long gravitational-wave detection efforts. For our most precisely timed pulsars, intrinsic instabilities of the pulsars or the observing system are shown to contribute to timing irregularities on a 5-year time-scale below the 100 ns level. Based on those results, realistic sensitivity curves for planned and ongoing timing array efforts are determined. We conclude that prospects for detection of a gravitational-wave background through pulsar timing array efforts within 5 years to a decade are good.

We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024−0719, J1045−4509, J1600−3053, J1603−7202, and J1730−2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437−4715 and J1909−3744 with M p = 1.44 ± 0.07 M and M p = 1.47 ± 0.03 M respectively. The improved orbital period-derivative measurement for PSR J0437−4715 results in a derived distance measurement at the 0.16% level of precision, D = 156.79 ± 0.25 pc, one of the most fractionally precise distance measurements of any star to date.

Interstellar scintillation (ISS), observed as time variation in the intensity of a compact radio source, is caused by small-scale structure in the electron density of the interstellar plasma. Dynamic spectra of ISS show modulation in radio frequency and time. Here we relate the (two-dimensional) power spectrum of the dynamic spectrum-the secondary spectrum-to the scattered image of the source. Recent work has identified remarkable parabolic arcs in secondary spectra. Each point in a secondary spectrum corresponds to interference between points in the scattered image with a certain Doppler shift and a certain delay. The parabolic arc corresponds to the quadratic relation between differential Doppler shift and delay through their common dependence on scattering angle. We show that arcs will occur in all media that scatter significant power at angles larger than the rms angle. Thus, effects such as source diameter, steep spectra, and dissipation scales, which truncate high angle scattering, also truncate arcs. Arcs are equally visible in simulations of nondispersive scattering. They are enhanced by anisotropic scattering when the spatial structure is elongated perpendicular to the velocity. In weak scattering the secondary spectrum is directly mapped from the scattered image, and this mapping can be inverted. We discuss additional observed phenomena including multiple arcs and reverse arclets oriented oppositely to the main arc. These phenomena persist for many refractive scattering times, suggesting that they are due to large-scale density structures, rather than low-frequency components of Kolmogorov turbulence.

Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparallelled measurement of the turbulent interstellar plasma at au scales.In this paper we demonstrate that without consideration of wavelengthindependent red-noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimise this additional noise, and factors, such as scintillation, that can exacerbate the problem.Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603−7202 of ∼ 2×10 −3 cm −3 pc for about 250 days. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously.