
We search for continuous gravitational waves (CGWs) produced by individual
supermassive blackhole binaries (SMBHBs) in circular orbits using
highcadence 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 highcadence observations are used to
improve the pulsar timing sensitivity across the GW frequency range of
$0.0085$ $\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 skyaveraged 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 skyaveraged strain amplitude of
lowfrequency CGWs to be $\lesssim1.4\times10^{14}$ at a reference frequency
of 20~nHz.

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
postNewtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance;
3. the universality of free fall (UFF) for strongly selfgravitating 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
selfgravitating 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 present a highprecision timing analysis of PSR J1933$$6211, a
millisecond pulsar (MSP) with a 3.5ms 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 rootmeansquare 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 periodmass
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.

Timing results for the blackwidow pulsar J20510827 are presented, using a
21year dataset from four European Pulsar Timing Array telescopes and the
Parkes radio telescope. This dataset, which is the longest published to date
for a blackwidow system, allows for an improved analysis that addresses
previously unknown biases. While secular variations, as identified in previous
analyses, are recovered, shortterm variations are detected for the first time.
Concurrently, a significant decrease of approx. 2.5x103 cm3 pc in the
dispersion measure associated with PSR J20510827 is measured for the first
time and improvements are also made to estimates of the proper motion. Finally,
PSR J20510827 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 pulsaryears 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 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 lowmetallicity mainsequence 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 mainsequence star and the pulsar accounts for the spin frequency
derivatives, which in turn resolves the distance discrepancy. Our observations
suggest that the pulsar and mainsequence 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 mainsequence companion star cannot have recycled the pulsar to
millisecond spin periods, an exotic formation scenario is required. We
demonstrate that this extremely wideorbit binary could have evolved from a
triple system that underwent an asymmetric supernova explosion, though find
that significant finetuning 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 highprecision timing of 42 radio millisecond pulsars (MSPs)
observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0
extends up to mid2014 and baselines range from 718 years. It forms the basis
for the stochastic gravitationalwave background, anisotropic background, and
continuouswave 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 semimajor axis. We find the NE2001
Galactic electron density model to be a better match to our parallax distances
(after correction from the LutzKelker 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 annualorbital 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 spinnoise and
dispersionmeasure (DM) variations, these models can include timing noise
unique to a single observing system, or frequency band. We show the improved
radiofrequency coverage and presence of overlapping data from different
observing systems in the IPTA data set enables us to separate both system and
banddependent 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 banddependent noise that is coherent between PTAs
which we interpret as coming from timevariable 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 longterm sensitivity of the
pulsar to a stochastic gravitationalwave (GW) background. By using a more
appropriate model, however, we can greatly improve a pulsar's sensitivity to
GWs. For example, including system and banddependent 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 spinnoise
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 longperiod 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 gravitationalwave
limits placed by individual PTAs by a factor of ~2 and provides a 2sigma 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, highquality
data they used.

The PSRIX backend is the primary pulsar timing instrument of the Effelsberg
100m radio telescope since early 2011. This new ROACHbased system enables
bandwidths up to 500 MHz to be recorded, significantly more than what was
possible with its predecessor, the EffelsbergBerkeley 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 phaseaveraged 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 ~25 times in signaltonoise ratio
achievable, which translates to an increase of ~25 times in pulse
timeofarrival (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 millisecondperiod 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 lowfrequency, 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 blackhole binaries with circular orbits.