
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 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.

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 suboptimal. We present the first
constraints on the angular distribution of a nanohertz stochastic GWB from
circular, inspiraldriven 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.