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 wavelength- independent 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 ~2x10^{-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.

A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.

We investigate the effects of gravitational waves (GWs) from a simulated population of binary super-massive black holes (SMBHs) on pulsar timing array datasets. We construct a distribution describing the binary SMBH population from an existing semi-analytic galaxy formation model. Using realizations of the binary SMBH population generated from this distribution, we simulate pulsar timing datasets with GW-induced variations. We find that the statistics of these variations do not correspond to an isotropic, stochastic GW background. The "Hellings & Downs" correlations between simulated datasets for different pulsars are recovered on average, though the scatter of the correlation estimates is greater than expected for an isotropic, stochastic GW background. These results are attributable to the fact that just a few GW sources dominate the GW-induced variations in every Fourier frequency bin of a 5-year dataset. Current constraints on the amplitude of the GW signal from binary SMBHs will be biased. Individual binary systems are likely to be detectable in 5-year pulsar timing array datasets where the noise is dominated by GW-induced variations. Searches for GWs in pulsar timing array data therefore need to account for the effects of individual sources of GWs.

Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.

We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the International Reference Ionosphere gave the best results. Little or no significant long-term variation in interstellar RM was found with limits typically about 0.1 rad m$^{-2}$ yr$^{-1}$ in absolute value. In a few cases, apparently significant RM variations over timescales of a few 100 days or more were seen. These are unlikely to be due to localised magnetised regions crossing the line of sight since the implied magnetic fields are too high. Most probably they are statistical fluctuations due to random spatial and temporal variations in the interstellar electron density and magnetic field along the line of sight.

Pulsar timing observations are used to place constraints on the rate of coalescence of supermassive black-hole (SMBH) binaries as a function of mass and redshift. In contrast to the indirect constraints obtained from other techniques, pulsar timing observations provide a direct constraint on the black-hole merger rate. This is possible since pulsar timing is sensitive to the gravitational waves (GWs) emitted by these sources in the final stages of their evolution. We find that upper bounds calculated from the recently published Parkes Pulsar Timing Array data are just above theoretical predictions for redshifts below 10. In the future, with improved timing precision and longer data spans, we show that a non-detection of GWs will rule out some of the available parameter space in a particular class of SMBH binary merger models. We also show that if we can time a set of pulsars to 10ns timing accuracy, for example, using the proposed Square Kilometre Array, it should be possible to detect one or more individual SMBH binary systems.

Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/noise ratios and show many previously undetected profile features. Thirteen of the 20 pulsars show emission over more than half of the pulse period. Polarization variations across the profiles are complex and the observed position angle variations are generally not in accord with the rotating-vector model for pulsar polarization. Never-the-less, the polarization properties are broadly similar to those of normal (non-millisecond) pulsars, suggesting that the basic radio emission mechanism is the same in both classes of pulsar. The results support the idea that radio emission from millisecond pulsars originates high in the pulsar magnetosphere, probably close to the emission regions for high-energy X-ray and gamma-ray emission. Rotation measures were obtained for all 20 pulsars, eight of which had no previously published measurements.

We search for the signature of an isotropic stochastic gravitational-wave background in pulsar timing observations using a frequency-domain correlation technique. These observations, which span roughly 12 yr, were obtained with the 64-m Parkes radio telescope augmented by public domain observations from the Arecibo Observatory. A wide range of signal processing issues unique to pulsar timing and not previously presented in the literature are discussed. These include the effects of quadratic removal, irregular sampling, and variable errors which exacerbate the spectral leakage inherent in estimating the steep red spectrum of the gravitational-wave background. These observations are found to be consistent with the null hypothesis, that no gravitational-wave background is present, with 76 percent confidence. We show that the detection statistic is dominated by the contributions of only a few pulsars because of the inhomogeneity of this data set. The issues of detecting the signature of a gravitational-wave background with future observations are discussed.

We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed.

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.

The Parkes Pulsar Timing Array project aims to make a direct detection of a gravitational-wave background through timing of millisecond pulsars. In this article, the main requirements for that endeavour are described and recent and ongoing progress is outlined. We demonstrate that the timing properties of millisecond pulsars are adequate and that technological progress is timely to expect a successful detection of gravitational waves within a decade, or alternatively to rule out all current predictions for gravitational wave backgrounds formed by supermassive black-hole mergers.

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.

Analysis of high-precision timing observations of an array of approx. 20 millisecond pulsars (a so-called "timing array") may ultimately result in the detection of a stochastic gravitational-wave 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 five-year timescale 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 five years to a decade are good.

The PULSE@Parkes project has been designed to monitor the rotation of radio pulsars over time spans of days to years. The observations are obtained using the Parkes 64-m and 12-m radio telescopes by Australian and international high school students. These students learn the basis of radio astronomy and undertake small projects with their observations. The data are fully calibrated and obtained with the state-of-the-art pulsar hardware available at Parkes. The final data sets are archived and are currently being used to carry out studies of 1) pulsar glitches, 2) timing noise, 3) pulse profile stability over long time scales and 4) the extreme nulling phenomenon. The data are also included in other projects such as gamma-ray observatory support and for the Parkes Pulsar Timing Array project. In this paper we describe the current status of the project and present the first scientific results from the Parkes 12-m radio telescope. We emphasise that this project offers a straightforward means to enthuse high school students and the general public about radio astronomy while obtaining scientifically valuable data sets.

Analysis of pulsar timing data-sets may provide the first direct detection of gravitational waves. This paper, the third in a series describing the mathematical framework implemented into the tempo2 pulsar timing package, reports on using tempo2 to simulate the timing residuals induced by gravitational waves. The tempo2 simulations can be used to provide upper bounds on the amplitude of an isotropic, stochastic, gravitational wave background in our Galaxy and to determine the sensitivity of a given pulsar timing experiment to individual, supermassive, binary black hole systems.

The first direct detection of gravitational waves may be made through observations of pulsars. The principal aim of pulsar timing array projects being carried out worldwide is to detect ultra-low frequency gravitational waves (f ~ 10^-9 to 10^-8 Hz). Such waves are expected to be caused by coalescing supermassive binary black holes in the cores of merged galaxies. It is also possible that a detectable signal could have been produced in the inflationary era or by cosmic strings. In this paper we review the current status of the Parkes Pulsar Timing Array project (the only such project in the Southern hemisphere) and compare the pulsar timing technique with other forms of gravitational-wave detection such as ground- and space-based interferometer systems.