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.

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 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 present the results of simultaneous radio and X-ray observations of PSR J1819-1458. Our 94-ks XMM-Newton observation of the high magnetic field 5*10^13 G pulsar reveals a blackbody spectrum (kT~130 eV) with a broad absorption feature, possibly composed of two lines at ~1.0 and ~1.3 keV. We performed a correlation analysis of the X-ray photons with radio pulses detected in 16.2 hours of simultaneous observations at 1-2 GHz with the Green Bank, Effelsberg, and Parkes telescopes, respectively. Both the detected X-ray photons and radio pulses appear to be randomly distributed in time. We find tentative evidence for a correlation between the detected radio pulses and X-ray photons on timescales of less than 10 pulsar spin periods, with the probability of this occurring by chance being 0.46%. This suggests that the physical process producing the radio pulses may also heat the polar-cap.

Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.

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.

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.

We describe our studies of the radio and high-energy properties of Rotating Radio Transients (RRATs). We find that the radio pulse intensity distributions are log-normal, with power-law tails evident in two cases. For the three RRATs with coverage over a wide range of frequency, the mean spectral index is -1.7\pm0.1, roughly in the range of normal pulsars. We do not observe anomalous magnetar-like spectra for any RRATs. Our 94-ks XMM-Newton observation of the high magnetic field RRAT J1819-1458 reveals a blackbody spectrum (kT ~130 eV) with an unusual absorption feature at ~1 keV. We find no evidence for X-ray bursts or other X-ray variability. We performed a correlation analysis of the X-ray photons with radio pulses detected in concurrent observations with the Green Bank, Effelsberg, and Parkes telescopes. We find no evidence for any correlations between radio pulse emission and X-ray photons, perhaps suggesting that sporadicity is not due to variations in magnetospheric particle density but to changes in beaming or coherence.

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

The European Pulsar Timing Array (EPTA) network is a collaboration between the five largest radio telescopes in Europe aiming to study the astrophysics of millisecond pulsars and to detect cosmological gravitational waves in the nano-Hertz regime. The advantages and techniques of handling the multi-telescope datasets of a number of sources will be presented. In addition, the results of the EPTA timing analysis of the pulsar-white dwarf binary PSR J1012+5307 will be reported. Specifically, the measurements for the first time for this system, of the parallax, the variation of the projected semi-major axis and of the orbital period. Finally, the derived stringent, theory independent limits on alternative theories of gravity, with the use of this ideal laboratory for strong- field gravity tests, will be presented.

We report on the high precision timing analysis of the pulsar-white dwarf binary PSR J1012+5307. Using 15 years of multi-telescope data from the European Pulsar Timing Array (EPTA) network, a significant measurement of the variation of the orbital period is obtained. Using this ideal strong-field gravity laboratory we derive theory independent limits for both the dipole radiation and the variation of the gravitational constant.

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.

We present results from the high precision timing analysis of the pulsar-white dwarf (WD) binary PSR J1012+5307 using 15 years of multi-telescope data. Observations were performed regularly by the European Pulsar Timing Array (EPTA) network, consisting of Effelsberg, Jodrell Bank, Westerbork and Nan\c{c}ay. All the timing parameters have been improved from the previously published values, most by an order of magnitude. In addition, a parallax measurement of $\pi = 1.2(3)$ mas is obtained for the first time for PSR J1012+5307, being consistent with the optical estimation from the WD companion. Combining improved 3D velocity information and models for the Galactic potential the complete evolutionary Galactic path of the system is obtained. A new intrinsic eccentricity upper limit of $e<8.4\times 10^{-7}$ is acquired, one of the smallest calculated for a binary system and a measurement of the variation of the projected semi-major axis also constrains the system's orbital orientation for the first time. It is shown that PSR J1012+5307 is an ideal laboratory for testing alternative theories of gravity. The measurement of the change of the orbital period of the system of $\dot{P}_{b} = 5(1)\times 10^{-14}$ is used to set an upper limit on the dipole gravitational wave emission that is valid for a wide class of alternative theories of gravity. Moreover, it is shown that in combination with other binary pulsars PSR J1012+5307 is an ideal system to provide self-consistent, generic limits, based only on millisecond pulsar data, for the dipole radiation and the variation of the gravitational constant $\dot{G}$.

We have used the 76-m Lovell, 94-m equivalent WSRT and 100-m Effelsberg radio telescopes to investigate the simultaneous single-pulse properties of the radio emitting magnetar AXP XTE J1810-197 at frequencies of 1.4, 4.8 and 8.35 GHz during May and July 2006. We study the magnetar's pulse-energy distributions which are found to be very peculiar as they are changing on time-scales of days and cannot be fit by a single statistical model. The magnetar exhibits strong spiky single giant-pulse-like subpulses, but they do not fit the definition of the giant pulse or giant micropulse phenomena. Measurements of the longitude-resolved modulation index reveal a high degree of intensity fluctuations on day-to-day time-scales and dramatic changes across pulse phase. We find the frequency evolution of the modulation index values differs significantly from what is observed in normal radio pulsars. We find that no regular drifting subpulse phenomenon is present at any of the observed frequencies at any observing epoch. However, we find a quasi-periodicity of the subpulses present in the majority of the observing sessions. A correlation analysis indicates a relationship between components from different frequencies. We discuss the results of our analysis in light of the emission properties of normal radio pulsars and a recently proposed model which takes radio emission from magnetars into consideration.

As part of a European Pulsar Network (EPN) multi-telescope observing campaign, we performed simultaneous multi-frequency observations at 1.4, 4.9 and 8.4 GHz during July 2006 and quasi-simultaneous multi-frequency observations from Decem- ber 2006 until July 2007 at 2.7, 4.9, 8.4, 14.6 and 32 GHz, in order to obtain flux density measurements and spectral features of the 5.5-sec radio-emitting magnetar AXP J1810-197. We monitored the spectral evolution of its pulse shape which consists of a main pulse (MP) and an interpulse (IP). We present the flux density spectrum of the average profile and of the separate pulse components of this first-known radio-emitting transient anomalous X-ray pulsar. We observe a decrease of the flux density by a factor of 10 within 8 months and follow the disappearance of one of the two main components. Although the spectrum is generally flat, we observe large fluctuations of the spectral index with time. For that reason we have made some measurements of modulation indices for individual pulses in order to also investigate the origin of these fluctuations.