We study the timing stability of three black widow pulsars, both in terms of their long-term spin evolution and their shorter-term orbital stability. The erratic timing behaviour and radio eclipses of the first two black widow pulsar systems discovered (PSRs B1957+20 and J2051$-$0827) was assumed to be representative for this class of pulsars. With several new black widow systems added to this population in the last decade, there are now several systems known that do not show these typical orbital variations or radio eclipses. We present timing solutions using 7$-$8 yrs of observations from four of the European Pulsar Timing Array telescopes for PSRs\,J0023+0923, J2214+3000 and J2234+0944, and confirm that two of these systems do not show any significant orbital variability over our observing time span, both in terms of secular or orbital parameters. The third pulsar PSR J0023+0923 shows orbital variability and we discuss the implications for the timing solution. Our results from the long-term timing of these pulsars provide several new or improved parameters compared to earlier works. We discuss our results regarding the stability of these pulsars, and the stability of the class of black widow pulsars in general, in the context of the binary parameters, and discuss the potential of the Roche-lobe filling factor of the companion star being an indicator for stability of these systems.

We are using the LOw-Frequency ARray (LOFAR) to perform the LOFAR Tied-Array All-Sky (LOTAAS) survey for pulsars and fast transients. Here we present the astrometric and rotational parameters of 20 pulsars discovered as part of LOTAAS. These pulsars have regularly been observed with LOFAR at 149 MHz and the Lovell telescope at 1532 MHz, supplemented by some observations with the Lovell telescope at 334 MHz and the Nancay Radio Telescope at 1484 MHz. Timing models are calculated for the 20 pulsars, some of which are among the slowest-spinning pulsars known. PSR J1236-0159 rotates with a period P ~ 3.6 s, while 5 additional pulsars show P > 2 s. Also, the spin-down rates Pdot are, on average, low, with PSR J0815+4611 showing Pdot ~ 4E-18. Some of the pulse profiles, generically single-peaked, present complex shapes evolving with frequency. Multi-frequency flux measurements show that these pulsars have generically relatively steep spectra but exceptions are present, with values ranging between ~ -4 and -1. Among the pulsar sample, a large fraction shows large single-pulse variability, with 4 pulsars being undetectable more than 15% of the time and one tentatively classified as a Rotating Radio Transient. Two single-peaked pulsars show drifting sub-pulses.

The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_\gamma = 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.

Dispersive delays due to the Solar wind introduce excess noise in high-precision pulsar timing experiments, and must be removed in order to achieve the accuracy needed to detect, e.g., low-frequency gravitational waves. In current pulsar timing experiments, this delay is usually removed by approximating the electron density distribution in the Solar wind either as spherically symmetric, or with a two-phase model that describes the contributions from both high- and low-speed phases of the Solar wind. However, no dataset has previously been available to test the performance and limitations of these models over extended timescales and with sufficient sensitivity. Here we present the results of such a test with an optimal dataset of observations of pulsar J0034-0534, taken with the German stations of LOFAR. We conclude that the spherical approximation performs systematically better than the two-phase model at almost all angular distances, with a residual root-mean-square (rms) given by the two-phase model being up to 28% larger than the result obtained with the spherical approximation. Nevertheless, the spherical approximation remains insufficiently accurate in modelling the Solar-wind delay (especially within 20 degrees of angular distance from the Sun), as it leaves timing residuals with rms values that reach the equivalent of 0.3 microseconds at 1400 MHz. This is because a spherical model ignores the large daily variations in electron density observed in the Solar wind. In the short term, broadband observations or simultaneous observations at low frequencies are the most promising way forward to correct for Solar-wind induced delay variations.

Context. High-precision pulsar-timing experiments are affected by temporal variations of the Dispersion Measure (DM), which are related to spatial variations in the interstellar electron content. Correcting for DM variations relies on the cold-plasma dispersion law which states that the dispersive delay varies with the squared inverse of the observing frequency. This may however give incorrect measurements if the probed electron content (and therefore the DM) varies with observing frequency, as is predicted theoretically. Aims. We study small-scale density variations in the ionised interstellar medium. These structures may lead to frequency-dependent DMs in pulsar signals and could inhibit the use of lower-frequency pulsar observations to correct time-variable interstellar dispersion in higher-frequency pulsar-timing data. Methods. We used high-cadence, low-frequency observations with three stations from the German LOng-Wavelength (GLOW) consortium, which are part of the LOw Frequency ARray (LOFAR). Specifically, 3.5 years of weekly observations of PSR J2219+4754 are presented. Results. We present the first detection of frequency-dependent DMs towards any interstellar object and a precise multi-year time-series of the time- and frequency-dependence of the measured DMs. The observed DM variability is significant and may be caused by extreme scattering events. Potential causes for frequency-dependent DMs are quantified and evaluated. Conclusions. We conclude that frequency-dependence of DMs has been reliably detected and is caused by small-scale (up to 10s of AUs) but steep density variations in the interstellar electron content. We find that long-term trends in DM variability equally affect DMs measured at both ends of our frequency band and hence the negative impact on long-term high-precision timing projects is expected to be limited.

We determined Faraday rotation measures (RMs) towards 137 pulsars in the northern sky, using Low-Frequency Array (LOFAR) observations at 110-190 MHz. This low-frequency RM catalogue, the largest to date, improves the precision of existing RM measurements on average by a factor of 20 - due to the low frequency and wide bandwidth of the data, aided by the RM synthesis method. We report RMs towards 25 pulsars for the first time. The RMs were corrected for ionospheric Faraday rotation to increase the accuracy of our catalogue to approximately 0.1 rad m$^{\rm -2}$. The ionospheric RM correction is currently the largest contributor to the measurement uncertainty. In addition, we find that the Faraday dispersion functions towards pulsars are extremely Faraday thin - mostly less than 0.001 rad m$^{\rm -2}$. We use these new precise RM measurements (in combination with existing RMs, dispersion measures, and distance estimates) to estimate the scale height of the Galactic halo magnetic field: 2.0$\pm$0.3 kpc for Galactic quadrants I and II above and below the Galactic plane (we also evaluate the scale height for these regions individually). Overall, our initial low-frequency catalogue provides valuable information about the 3-D structure of the Galactic magnetic field.

Faraday rotation of polarized emission from pulsars measured at radio frequencies provides a powerful tool to investigate the interstellar and interplanetary magnetic fields. However, besides being sensitive to the astrophysical media, pulsar observations in radio are affected by the highly time-variable ionosphere. In this article, the amount of ionospheric Faraday rotation has been computed by assuming a thin layer model. For this aim, ionospheric maps of the free electron density (based on Global Positioning System data) and semi-empirical geomagnetic models are needed. Through the data of five highly polarized pulsars observed with the individual German LOw-Frequency ARray stations, we investigate the performances of the ionospheric modelling. In addition, we estimate the parameters of the systematics and the correlated noise generated by the residual unmodelled ionospheric effects, and show the comparison of the different free-electron density maps. For the best ionospheric maps, we have found that the rotation measure corrections on one-year timescales after subtraction of diurnal periodicity are accurate to $\sim$ 0.06--0.07 rad m$^{-2}$.

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.

We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $\mu$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz.

We operate the six German stations of the LOw Frequency ARray as standalone telescopes to observe more than 100 pulsars every week. To date, we have collected almost four years of high-quality data at an unprecedented weekly cadence. This allows us to perform a wide variety of analyses aimed at characterising the magnetoionic plasma crossed by pulsar radiation. In particular, our studies are focused on electron density variations in the interstellar and interplanetary media, the Galactic and interplanetary magnetic field, scintillation, and extreme scattering events. Here we report the first results from our Solar wind monitoring campaign.

PSR J0922+0638 (B0919+06) shows unexplained anomalous variations in the on-pulse phase, where the pulse appears to episodically move to an earlier longitude for a few tens of rotations before reverting to the usual phase for approximately several hundred to more than a thousand rotations. These events, where the pulse moves in phase by up to 5$^{\circ}$, have been previously detected in observations from $\sim$300 to 2000 MHz. We present simultaneous observations from the Effelsberg 100-m radio telescope at 1350 MHz and the Bornim (Potsdam) station of the LOw Frequency ARray at 150 MHz. Our observations present the first evidence for an absence of the anomalous phase-shifting behaviour at 150 MHz. Instead, the observed intensity at the usual pulse-phase typically decreases, often showing a pseudo-nulling feature corresponding to the times when phase shifts are observed at 1350 MHz. The presence of weak emission at the usual pulse-phase supports the theory that these shifts may result from processes similar to the 'profile-absorption' expected to operate for PSR J0814+7429 (B0809+74). A possible mechanism for this could be intrinsic variations of the emission within the pulsar's beam combined with absorption by expanding shells of electrons in the line of sight.

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 have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar medium (IISM) are the most likely cause. The putative structures in the IISM causing the profile variation are roughly half-way between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We consider different models for the structures. Under the assumption of spherical symmetry, their implied average electron density is $\overline{n}_e \sim 100$ cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than the average pulsar population visible to LOFAR, similar profile variations would not have been identified in most pulsars, suggesting that subtle profile variations in low-frequency profiles might be more common than we have observed to date. Systematic studies of these variations at low frequencies can provide a new tool to investigate the proprieties of the IISM and the limits to the precision of pulsar timing.

We present a high-precision timing analysis of PSR J1933$-$6211, a millisecond pulsar (MSP) with a 3.5-ms spin period and a white dwarf (WD) companion, using data from the Parkes radio telescope. Since we have accurately measured the polarization properties of this pulsar we have applied the matrix template matching approach in which the times of arrival are measured using full polarimetric information. We achieved a weighted root-mean-square timing residuals (rms) of the timing residuals of 1.23 $\rm \mu s$, 15.5$\%$ improvement compared to the total intensity timing analysis. After studying the scintillation properties of this pulsar we put constraints on the inclination angle of the system. Based on these measurements and on $\chi^2$ mapping we put a 2-$\sigma$ upper limit on the companion mass (0.44 M$_\odot$). Since this mass limit cannot reveal the nature of the companion we further investigate the possibility of the companion to be a He WD. Applying the orbital period-mass relation for such WDs, we conclude that the mass of a He WD companion would be about 0.26$\pm$0.01 M$_\odot$ which, combined with the measured mass function and orbital inclination limits, would lead to a light pulsar mass $\leqslant$ 1.0 M$_\odot$. This result seems unlikely based on current neutron star formation models and we therefore conclude that PSR J1933$-$6211 most likely has a CO WD companion, which allows for a solution with a more massive pulsar.

Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.

Timing results for the black-widow pulsar J2051-0827 are presented, using a 21-year dataset from four European Pulsar Timing Array telescopes and the Parkes radio telescope. This dataset, which is the longest published to date for a black-widow system, allows for an improved analysis that addresses previously unknown biases. While secular variations, as identified in previous analyses, are recovered, short-term variations are detected for the first time. Concurrently, a significant decrease of approx. 2.5x10-3 cm-3 pc in the dispersion measure associated with PSR J2051-0827 is measured for the first time and improvements are also made to estimates of the proper motion. Finally, PSR J2051-0827 is shown to have entered a relatively stable state suggesting the possibility of its eventual inclusion in pulsar timing arrays.

We present evidence for a small glitch in the spin evolution of the millisecond pulsar J0613$-$0200, using the EPTA Data Release 1.0, combined with Jodrell Bank analogue filterbank TOAs recorded with the Lovell telescope and Effelsberg Pulsar Observing System TOAs. A spin frequency step of 0.82(3) nHz and frequency derivative step of ${-1.6(39) \times 10^{-19}\,\text{Hz} \ \text{s}^{-1}}$ are measured at the epoch of MJD 50888(30). After PSR B1821$-$24A, this is only the second glitch ever observed in a millisecond pulsar, with a fractional size in frequency of ${\Delta \nu/\nu=2.5(1) \times 10^{-12}}$, which is several times smaller than the previous smallest glitch. PSR J0613$-$0200 is used in gravitational wave searches with pulsar timing arrays, and is to date only the second such pulsar to have experienced a glitch in a combined 886 pulsar-years of observations. We find that accurately modelling the glitch does not impact the timing precision for pulsar timing array applications. We estimate that for the current set of millisecond pulsars included in the International Pulsar Timing Array, there is a probability of $\sim 50$% that another glitch will be observed in a timing array pulsar within 10 years.

We present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations ${\rm Dec}> 8^\circ$ and Galactic latitudes $|{\rm Gb}| > 3^\circ$, regardless of their expected flux densities and scattering times. Each pulsar was observed for $\geq 20$ minutes in the contiguous frequency range of 110--188 MHz. Full-Stokes data were recorded. We present the dispersion measures, flux densities, and calibrated total intensity profiles for the 158 pulsars detected in the sample. The median uncertainty in census dispersion measures ($1.5 \times 10^{-3}$ pc cm$^{-3}$) is ten times smaller, on average, than in the ATNF pulsar catalogue. We combined census flux densities with those in the literature and fitted the resulting broadband spectra with single or broken power-law functions. For 48 census pulsars such fits are being published for the first time. Typically, the choice between single and broken power-laws, as well as the location of the spectral break, were highly influenced by the spectral coverage of the available flux density measurements. In particular, the inclusion of measurements below 100 MHz appears essential for investigating the low-frequency turnover in the spectra for most of the census pulsars. For several pulsars, we compared the spectral indices from different works and found the typical spread of values to be within 0.5--1.5, suggesting a prevailing underestimation of spectral index errors in the literature. The census observations yielded some unexpected individual source results, as we describe in the paper. Lastly, we will provide this unique sample of wide-band, low-frequency pulse profiles via the European Pulsar Network Database.

We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.

We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.

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.

The PSRIX backend is the primary pulsar timing instrument of the Effelsberg 100-m radio telescope since early 2011. This new ROACH-based system enables bandwidths up to 500 MHz to be recorded, significantly more than what was possible with its predecessor, the Effelsberg-Berkeley Pulsar Processor (EBPP). We review the first four years of PSRIX timing data for 33 pulsars collected as part of the monthly European Pulsar Timing Array (EPTA) observations. We describe the automated data analysis pipeline, CoastGuard, that we developed to reduce these observations. We also introduce TOASTER, the EPTA timing database used to store timing results, processing information and observation metadata. Using these new tools, we measure the phase-averaged flux densities at 1.4 GHz of all 33 pulsars. For 7 of these pulsars, our flux density measurements are the first values ever reported. For the other 26 pulsars, we compare our flux density measurements with previously published values. By comparing PSRIX data with EBPP data, we find an improvement of ~2-5 times in signal-to-noise ratio achievable, which translates to an increase of ~2-5 times in pulse time-of-arrival (TOA) precision. We show that such an improvement in TOA precision will improve the sensitivity to the stochastic gravitational wave background. Finally, we showcase the flexibility of the new PSRIX backend by observing several millisecond-period pulsars (MSPs) at 5 and 9 GHz. Motivated by our detections, we discuss the potential for complementing existing pulsar timing array data sets with MSP monitoring campaigns at these frequencies.

The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.

We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of LOFAR profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves.

We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $\Omega_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $G\mu/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $G\mu/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $\Omega^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

Measurements of pulsar flux densities are of great importance for understanding the pulsar emission mechanism and for predictions of pulsar survey yields and the pulsar population at large. Typically these flux densities are determined from phase-averaged "pulse profiles", but this method has limited applicability at low frequencies because the observed pulses can easily be spread out by interstellar effects like scattering or dispersion, leading to a non-pulsed continuum component that is necessarily ignored in this type of analysis. In particular for the class of the millisecond pulsars (MSPs) at frequencies below 200MHz, such interstellar effects can seriously compromise de- tectability and measured flux densities. In this paper we investigate MSP spectra based on a complementary approach, namely through investigation of archival con- tinuum imaging data. Even though these images lose sensitivity to pulsars since the on-pulse emission is averaged with off-pulse noise, they are insensitive to effects from scattering and provide a reliable way to determine the flux density and spectral indices of MSPs based on both pulsed and unpulsed components. Using the 74MHz VLSSr as well as the 325MHz WENSS and 1.4GHz NVSS catalogues, we investigate the imaging flux densities of MSPs and evaluate the likelihood of spectral turn-overs in this population. We determine three new MSP spectral indices and identify six new MSPs with likely spectral turn-overs.

The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.

We present the highest-quality polarisation profiles to date of 16 non-recycled pulsars and four millisecond pulsars, observed below 200 MHz with the LOFAR high-band antennas. Based on the observed profiles, we perform an initial investigation of expected observational effects resulting from the propagation of polarised emission in the pulsar magnetosphere and the interstellar medium. The predictions of magnetospheric birefringence in pulsars have been tested using spectra of the pulse width and fractional polarisation from multifrequency data. The derived spectra offer only partial support for the expected effects of birefringence on the polarisation properties, with only about half of our sample being consistent with the model's predictions. It is noted that for some pulsars these measurements are contaminated by the effects of interstellar scattering. For a number of pulsars in our sample, we have observed significant variations in the amount of Faraday rotation as a function of pulse phase, which is possibly an artefact of scattering. These variations are typically two orders of magnitude smaller than that observed at 1400 MHz by Noutsos et al. (2009), for a different sample of southern pulsars. In this paper we present a possible explanation for the difference in magnitude of this effect between the two frequencies, based on scattering. Finally, we have estimated the magnetospheric emission heights of low-frequency radiation from four pulsars, based on the phase lags between the flux-density and the PA profiles, and the theoretical framework of Blaskiewicz, Cordes & Wasserman (1991). These estimates yielded heights of a few hundred km; at least for PSR B1133+16, this is consistent with emission heights derived based on radius-to-frequency mapping, but is up to a few times larger than the recent upper limit based on pulsar timing.

On a time scale of years to decades, gravitational wave (GW) astronomy will become a reality. Low frequency (nanoHz) GWs are detectable through long-term timing observations of the most stable pulsars. Radio observatories worldwide are currently carrying out observing programmes to detect GWs, with data sets being shared through the International Pulsar Timing Array project. One of the most likely sources of low frequency GWs are supermassive black hole binaries (SMBHBs), detectable as a background due to a large number of binaries, or as continuous or burst emission from individual sources. No GW signal has yet been detected, but stringent constraints are already being placed on galaxy evolution models. The SKA will bring this research to fruition. In this chapter, we describe how timing observations using SKA1 will contribute to detecting GWs, or can confirm a detection if a first signal already has been identified when SKA1 commences observations. We describe how SKA observations will identify the source(s) of a GW signal, search for anisotropies in the background, improve models of galaxy evolution, test theories of gravity, and characterise the early inspiral phase of a SMBHB system. We describe the impact of the large number of millisecond pulsars to be discovered by the SKA; and the observing cadence, observation durations, and instrumentation required to reach the necessary sensitivity. We describe the noise processes that will influence the achievable precision with the SKA. We assume a long-term timing programme using the SKA1-MID array and consider the implications of modifications to the current design. We describe the possible benefits from observations using SKA1-LOW. Finally, we describe GW detection prospects with SKA1 and SKA2, and end with a description of the expectations of GW astronomy.

The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

Recent VLBI measurements of the astrometric parameters of the millisecond pulsar J0218+4232 by Du et al. have suggested this pulsar is as distant as 6.3 kpc. At such a large distance, the large {\gamma}-ray flux observed from this pulsar would make it the most luminous {\gamma}-ray pulsar known. This luminosity would exceed what can be explained by the outer gap and slot-gap pulsar emission models, potentially placing important and otherwise elusive constraints on the pulsar emission mechanism. We show that the VLBI parallax measurement is dominated by the Lutz-Kelker bias. When this bias is corrected for, the most likely distance for this pulsar is 3.15(+0.85/-0.60) kpc. This revised distance places the luminosity of PSR J0218+4232 into a range where it does not challenge any of the standard theories of the pulsar emission mechanism.

In recent years, instrumentation enabling pulsar observations with unprecedentedly high fractional bandwidth has been under development which can be used to substantially improve the precision of pulsar timing experiments. The traditional template-matching method used to calculate pulse times-of-arrival (ToAs), may not function effectively on these broadband data due to a variety of effects such as diffractive scintillation in the interstellar medium, profile variation as a function of frequency, dispersion measure (DM) evolution and so forth. In this paper, we describe the channelised Discrete Fourier Transform method that can greatly mitigate the influence of the aforementioned effects when measuring ToAs from broadband timing data. The method is tested on simulated data, and its potential in improving timing precision is shown. We further apply the method to PSR J1909$-$3744 data collected at the Nan\c{c}ay Radio Telescope with the Nan\c{c}ay Ultimate Pulsar Processing Instrument. We demonstrate a removal of systematics due to the scintillation effect as well as improvement on ToA measurement uncertainties. Our method also determines temporal variations in dispersion measure, which are consistent with multi-channel timing approaches used earlier.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be ${\Omega}_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

We report on the setup and initial discoveries of the Northern High Time Resolution Universe survey for pulsars and fast transients, the first major pulsar survey conducted with the 100-m Effelsberg radio telescope and the first in 20 years to observe the whole northern sky at high radio frequencies. Using a newly developed 7-beam receiver system combined with a state-of-the-art polyphase filterbank, we record an effective bandwidth of 240 MHz in 410 channels centred on 1.36 GHz with a time resolution of 54 $\mu$s. Such fine time and frequency resolution increases our sensitivity to millisecond pulsars and fast transients, especially deep inside the Galaxy, where previous surveys have been limited due to intra-channel dispersive smearing. To optimise observing time, the survey is split into three integration regimes dependent on Galactic latitude, with 1500-s, 180-s and 90-s integrations for latitude ranges $|b|<3.5^{\circ}$, $|b|<15^{\circ}$ and $|b|>15^{\circ}$, respectively. The survey has so far resulted in the discovery of 15 radio pulsars, including a pulsar with a characteristic age of $\sim18$ kyr, {PSR J2004+3429}, and a highly eccentric, binary millisecond pulsar, {PSR J1946+3417}. All newly discovered pulsars are timed using the 76-m Lovell radio telescope at the Jodrell Bank Observatory and the Effelsberg radio telescope. We present timing solutions for all newly discovered pulsars and discuss potential supernova remnant associations for {PSR J2004+3429}.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

Many physically motivated extensions to general relativity (GR) predict significant deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 +/- 0.04 solar mass pulsar in a 2.46-hr orbit with a 0.172 +/- 0.003 solar mass white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.

The flux density spectra of radio pulsars are known to be steep and, to first order, described by a power-law relationship of the form S_{\nu} \propto \nu^{\alpha}, where S_{\nu} is the flux density at some frequency \nu and \alpha is the spectral index. Although measurements of \alpha have been made over the years for several hundred pulsars, a study of the intrinsic distribution of pulsar spectra has not been carried out. From the result of pulsar surveys carried out at three different radio frequencies, we use population synthesis techniques and a likelihood analysis to deduce what underlying spectral index distribution is required to replicate the results of these surveys. We find that in general the results of the surveys can be modelled by a Gaussian distribution of spectral indices with a mean of -1.4 and unit standard deviation. We also consider the impact of the so-called "Gigahertz-peaked spectrum" pulsars. The fraction of peaked spectrum sources in the population with significant turn-over at low frequencies appears to be at most 10%. We demonstrate that high-frequency (>2 GHz) surveys preferentially select flatter-spectrum pulsars and the converse is true for lower-frequency (<1 GHz) surveys. This implies that any correlations between \alpha and other pulsar parameters (for example age or magnetic field) need to carefully account for selection biases in pulsar surveys. We also expect that many known pulsars which have been detected at high frequencies will have shallow, or positive, spectral indices. The majority of pulsars do not have recorded flux density measurements over a wide frequency range, making it impossible to constrain their spectral shapes. We also suggest that such measurements would allow an improved description of any populations of pulsars with 'non-standard' spectra.

Using the 100-m Effelsberg radio telescope operating at 1.36 GHz, we have performed a targeted radio pulsar survey of 289 unassociated gamma-ray sources discovered by the Large Area Telescope (LAT) aboard the Fermi satellite and published in the 1FGL catalogue (Abdo et al., 2010). This survey resulted in the discovery of millisecond pulsar J1745+1017, which resides in a short-period binary system with a low-mass companion, Mmin ~ 0.0137 Msun, indicative of `Black Widow' type systems. A two-year timing campaign has produced a refined radio ephemeris, accurate enough to allow for phase-folding of the LAT photons, resulting in the detection of a dual-peaked gamma-ray light-curve, proving that PSR J1745+1017 is the source responsible for the gamma-ray emission seen in 1FGL J1745.5 + 1018 (2FGL J1745.6+1015; Nolan et al., 2012). We find the gamma-ray spectrum of PSR J1745+1017 to be well modelled by an exponentially-cut-off power law with cut-off energy 3.2 GeV and photon index 1.6. The observed sources are known to contain a further 10 newly discovered pulsars which were undetected in this survey. Our radio observations of these sources are discussed and in all cases limiting flux densities are calculated. The reasons behind the seemingly low yield of discoveries are also discussed.

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.

Pulsars in close binary systems have provided some of the most stringent tests of strong-field gravity to date. The pulsar--white-dwarf binary system J1141-6545 is specifically interesting due to its gravitational asymmetry which makes it one of the most powerful probes of tensor-scalar theories of gravity. We give an overview of current gravitational tests provided by the J1141-6545 binary system and comment on how anomalous accelerations, geodetic precession and timing instabilities may be prevented from limiting future tests of gravity to come from this system.

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.

Accurate distances to pulsars can be used for a variety of studies of the Galaxy and its electron content. However, most distance measures to pulsars have been derived from the absorption (or lack thereof) of pulsar emission by Galactic HI gas, which typically implies that only upper or lower limits on the pulsar distance are available. We present a critical analysis of all measured HI distance limits to pulsars and other neutron stars, and translate these limits into actual distance estimates through a likelihood analysis that simultaneously corrects for statistical biases. We also apply this analysis to parallax measurements of pulsars in order to obtain accurate distance estimates and find that the parallax and HI distance measurements are biased in different ways, because of differences in the sampled populations. Parallax measurements typically underestimate a pulsar's distance because of the limited distance to which this technique works and the consequential strong effect of the Galactic pulsar distribution (i.e. the original Lutz-Kelker bias), in HI distance limits, however, the luminosity bias dominates the Lutz-Kelker effect, leading to overestimated distances because the bright pulsars on which this technique is applicable are more likely to be nearby given their brightness.

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

(abridged) We report the results of a 10-year timing campaign on PSR J1738+0333, a 5.85-ms pulsar in a low-eccentricity 8.5-hour orbit with a low-mass white dwarf companion (...) The measurements of proper motion and parallax allow for a precise subtraction of the kinematic contribution to the observed orbital decay; this results in a significant measurement of the intrinsic orbital decay: (-25.9 +/- 3.2) \times 10^{-15} s/s. This is consistent with the orbital decay from the emission of gravitational waves predicted by general relativity, (-27.7 +1.5/-1.9) \times 10^{-15} s/s (...). This agreement introduces a tight upper limit on dipolar gravitational wave emission, a prediction of most alternative theories of gravity for asymmetric binary systems such as this. We use this limit to derive the most stringent constraints ever on a wide class of gravity theories, where gravity involves a scalar field contribution. When considering general scalar-tensor theories of gravity, our new bounds are more stringent than the best current solar-system limits over most of the parameter space, and constrain the matter-scalar coupling constant {\alpha}_0^2 to be below the 10^{-5} level. For the special case of the Jordan-Fierz-Brans-Dicke, we obtain the one-sigma bound {\alpha}_0^2 < 2 \times 10^{-5}, which is within a factor two of the Cassini limit. We also use our limit on dipolar gravitational wave emission to constrain a wide class of theories of gravity which are based on a generalization of Bekenstein's Tensor-Vector-Scalar gravity (TeVeS), a relativistic formulation of Modified Newtonian Dynamics (MOND).

Pulsar timing observations are usually analysed with least-square-fitting procedures under the assumption that the timing residuals are uncorrelated (statistically "white"). Pulsar observers are well aware that this assumption often breaks down and causes severe errors in estimating the parameters of the timing model and their uncertainties. Ad hoc methods for minimizing these errors have been developed, but we show that they are far from optimal. Compensation for temporal correlation can be done optimally if the covariance matrix of the residuals is known using a linear transformation that whitens both the residuals and the timing model. We adopt a transformation based on the Cholesky decomposition of the covariance matrix, but the transformation is not unique. We show how to estimate the covariance matrix with sufficient accuracy to optimize the pulsar timing analysis. We also show how to apply this procedure to estimate the spectrum of any time series with a steep red power-law spectrum, including those with irregular sampling and variable error bars, which are otherwise very difficult to analyse.

Gravitational waves are predicted by Einstein's theory of general relativity as well as other theories of gravity. The rotational stability of the fastest pulsars means that timing of an array of these objects can be used to detect and investigate gravitational waves. Simultaneously, however, pulsar timing is used to estimate spin period, period derivative, astrometric, and binary parameters. Here we calculate the effects that a stochastic background of gravitational waves has on pulsar timing parameters through the use of simulations and data from the millisecond pulsars PSR J0437--4715 and PSR J1713+0747. We show that the reported timing uncertainties become underestimated with increasing background amplitude by up to a factor of $\sim10$ for a stochastic gravitational-wave background amplitude of $A=5\times 10^{-15}$, where $A$ is the amplitude of the characteristic strain spectrum at one-year gravitational wave periods. We find evidence for prominent low-frequency spectral leakage in simulated data sets including a stochastic gravitational-wave background. We use these simulations along with independent Very Long Baseline Interferometry (VLBI) measurements of parallax to set a 2--sigma upper limit of $A\le9.1\times 10^{-14}$. We find that different supermassive black hole assembly scenarios do not have a significant effect on the calculated upper limits. We also test the effects that ultralow--frequency (10$^{-12}$--10$^{-9}$ Hz) gravitational waves have on binary pulsar parameter measurements and find that the corruption of these parameters is less than those due to $10^{-9}$--$10^{-7}$ Hz gravitational waves.

Timing pulses of pulsars has proved to be a most powerful technique useful to a host of research areas in astronomy and physics. Importantly, the precision of this timing is not only affected by radiometer noise, but also by intrinsic pulse shape changes, interstellar medium (ISM) evolution, instrumental distortions, etc. In this paper we review the known causes of pulse shape variations and assess their effect on the precision and accuracy of a single measurement of pulse arrival time with current instrumentation. Throughout this analysis we use the brightest and most precisely timed millisecond pulsar (MSP), PSR J0437-4715, as a case study, and develop a set of diagnostic tools to evaluate profile stability in timing observations. We conclude that most causes of distortion can be either corrected by state-of-the-art techniques or taken into account in the estimation of time-of-arrival (TOA) uncertainties. The advent of a new generation of radio telescopes (e.g. the Square Kilometre Array, SKA), and their increase in collecting area has sparked speculation about the timing precision achievable through increases in gain. Based on our analysis of current data, we predict that for normal-brightness MSPs a TOA precision of between 80 and 230 ns can be achieved at 1.4 GHz with 10-minute integrations by the SKA. The actual rms timing residuals for each pulsar will be approximately at the same level only if all the other influences on timing precision (e.g. ISM, spin noise) are either corrected, modelled, or negligible.

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.

We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, {\mu}_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.

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.

(abridged) PSR J1903+0327, a millisecond pulsar in an eccentric (e = 0.44) 95-day orbit with a (~ 1Msun) companion poses a challenge to our understanding of stellar evolution in binary and multiple-star systems. Here we describe optical and radio observations which rule out most of the scenarios proposed to explain formation of this system. Radio timing measurements of three post-Keplerian effects yield the most precise measurement of the mass of a millisecond pulsar to date: 1.667 +/- 0.021 solar masses (99.7% confidence limit) (...). Optical spectroscopy of a proposed main sequence counterpart show that its orbital motion mirrors the pulsar's 95-day orbit; being therefore its binary companion (...) The optical detection also provides a measurement of the systemic radial velocity of the binary; this and the proper motion measured from pulsar timing allow the determination of the systemic 3-D velocity in the Galaxy. We find that the system is always within 270 pc of the plane of the Galaxy, but always more than 3 kpc away from the Galactic centre. Thus an exchange interaction in a dense stellar environment (like a globular cluster or the Galactic centre) is not likely to be the origin of this system. We suggest that after the supernova that formed it, the neutron star was in a tight orbit with a main-sequence star, the present companion was a tertiary farther out. The neutron star then accreted matter from its evolving inner MS companion, forming a millisecond pulsar. The former donor star then disappears, either due to a chaotic 3-body interaction with the outer star (caused by the expansion of the inner orbit that necessarily results from mass transfer), or in the case of a very compact inner system, due to ablation/accretion by the newly formed millisecond pulsar.

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.

We present simulations of scattering phenomena which are important in pulsar observations, but which are analytically intractable. The simulation code, which has also been used for solar wind and atmospheric scattering problems, is available from the authors. These simulations reveal an unexpectedly important role of dispersion in combination with refraction. We demonstrate the effect of analyzing observations which are shorter than the refractive scale. We examine time-of-arrival fluctuations in detail: showing their correlation with intensity and dispersion measure; providing a heuristic model from which one can estimate their contribution to pulsar timing observations; and showing that much of the effect can be corrected making use of measured intensity and dispersion. Finally, we analyze observations of the millisecond pulsar J0437$-$4715, made with the Parkes radio telescope, that show timing fluctuations which are correlated with intensity. We demonstrate that these timing fluctuations can be corrected, but we find that they are much larger than would be expected from scattering in a homogeneous turbulent plasma with isotropic density fluctuations. We do not have an explanation for these timing fluctuations.

Lutz & Kelker showed that parallax measurements are systematically overestimated because they do not properly account for the larger volume of space that is sampled at smaller parallax values. We apply their analysis to neutron stars, incorporating the bias introduced by the intrinsic radio luminosity function and a realistic Galactic population model for neutron stars. We estimate the bias for all published neutron star parallax measurements and find that measurements with less than ~95% certainty, are likely to be significantly biased. Through inspection of historic parallax measurements, we confirm the described effects in optical and radio measurements, as well as in distance estimates based on interstellar dispersion measures. The potential impact on future tests of relativistic gravity through pulsar timing and on X-ray--based estimates of neutron star radii is briefly discussed.

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.

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.

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.

(Abridged) This thesis presents long-term timing results on 20 millisecond pulsars (MSPs). It has been predicted that such timing may detect gravitational waves (GWs) - a major (but untested) prediction of general relativity. Our results demonstrate that most of the investigated MSPs have sufficient stability to enable such detection experiments. Furthermore, the timing data of the brightest few sources shows that timing at sub-100 ns precision may be achievable and that, therefore, GW detection within the next decade is likely. Finally, we use our data to place a strong limit on the strength of a predicted background of GWs.

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.

Using the recently upgraded Long Baseline Array, we have measured the trigonometric parallax of PSR J0437-4715 to better than 1% precision, the most precise pulsar distance determination made to date. Comparing this VLBI distance measurement to the kinematic distance obtained from pulsar timing, which is calculated from the pulsar's proper motion and apparent rate of change of orbital period, gives a precise limit on the unmodeled relative acceleration between the Solar System and PSR J0437-4715, which can be used in a variety of applications. Firstly, it shows that Newton's gravitational constant G is stable with time (\dot{G}/G = (-5 +- 26) x 10^{-13} yr^{-1}, 95% confidence). Secondly, if a stochastic gravitational wave background existed at the currently quoted limit, this null result would fail ~50% of the time. Thirdly, it excludes Jupiter-mass planets within 226 AU of the Sun in 50% of the sky (95% confidence). Finally, the ~1% agreement of the parallax and orbital period derivative distances provides a fundamental confirmation of the parallax distance method upon which all astronomical distances are based.

Pulsars in close binary orbit around another neutron star or a massive white dwarf make ideal laboratories for testing the predictions of gravitational radiation and self-gravitational effects. We report new timing measurements of the pulsar-white-dwarf binary PSR J1141-6545, providing strong evidence that such asymmetric systems have gravitational wave losses that are consistent with general relativity. The orbit is found to be decaying at a rate of $1.04\pm0.06$ times the general relativistic prediction and the Shapiro delay is consistent with the orbital inclination angle derived from scintillation measurements. The system provides a unique test-bed for tensor-scalar theories of gravity; our current measurements place stringent constraints in the theory space, with a limit of $\alpha_0^2 < 2.1 \times 10^{-5}$ for weakly non-linear coupling and an asymptotic limit of $\alpha_0^2 < 3.4 \times 10^{-6}$ for strongly non-linear coupling, where $\alpha_0$ is the linear coupling strength of matter to an underlying scalar field. This asymptotic limit is nearly three times smaller than the Cassini bound ($\alpha_0^2 \approx 10^{-5}$).

Analysis of ten years of high-precision timing data on the millisecond pulsar PSR J0437-4715 has resulted in a model-independent kinematic distance based on an apparent orbital period derivative, Pbdot, determined at the 1.5% level of precision (Dk = 157.0 +/- 2.4 pc), making it one of the most accurate stellar distance estimates published to date. The discrepancy between this measurement and a previously published parallax distance estimate is attributed to errors in the DE200 Solar System ephemerides. The precise measurement of Pbdot allows a limit on the variation of Newton's gravitational constant, |Gdot/G| < 23 x 10^{-12} 1/yr. We also constrain any anomalous acceleration along the line of sight to the pulsar to |a(Sun)/c| < 1.5 x 10^{-18} 1/s at 95% confidence, and derive a pulsar mass, m(psr) = 1.76 +/- 0.20 M, one of the highest estimates so far obtained.

We present an analysis of the variations seen in the dispersion measures (DMs) of 20 millisecond pulsars observed as part of the Parkes Pulsar Timing Array project. We carry out a statistically rigorous structure function analysis for each pulsar and show that the variations seen for most pulsars are consistent with those expected for an interstellar medium characterised by a Kolmogorov turbulence spectrum. The structure functions for PSRs J1045-4509 and J1909-3744 provide the first clear evidence for a large inner scale, possibly due to ion-neutral damping. We also show the effect of the solar wind on the DMs and show that the simple models presently implemented into pulsar timing packages cannot reliably correct for this effect. For the first time we clearly show how DM variations affect pulsar timing residuals and how they can be corrected in order to obtain the highest possible timing precision. Even with our presently limited data span, the residuals (and all parameters derived from the timing) for six of our pulsars have been significantly improved by correcting for the DM variations.

Using a statistically rigorous analysis method, we place limits on the existence of an isotropic stochastic gravitational wave background using pulsar timing observations. We consider backgrounds whose characteristic strain spectra may be described as a power-law dependence with frequency. Such backgrounds include an astrophysical background produced by coalescing supermassive black-hole binary systems and cosmological backgrounds due to relic gravitational waves and cosmic strings. Using the best available data, we obtain an upper limit on the energy density per unit logarithmic frequency interval of \Omega^{\rm SMBH}_g(1/8yr) h^2 <= 1.9 x 10^{-8} for an astrophysical background which is five times more stringent than the earlier Kaspi et al. (1994) limit of 1.1 x 10^{-7}. We also provide limits on a background due to relic gravitational waves and cosmic strings of \Omega^{\rm relic}_g(1/8yr) h^2 <= 2.0 x 10^{-8} and \Omega^{\rm cs}_g(1/8yr) h^2 <= 1.9 x 10^{-8} respectively. All of the quoted upper limits correspond to a 0.1% false alarm rate together with a 95% detection rate. We discuss the physical implications of these results and highlight the future possibilities of the Parkes Pulsar Timing Array project. We find that our current results can 1) constrain the merger rate of supermassive binary black hole systems at high red shift, 2) rule out some relationships between the black hole mass and the galactic halo mass, 3) constrain the rate of expansion in the inflationary era and 4) provide an upper bound on the dimensionless tension of a cosmic string background.