
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.

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

We investigate the constraints that can be placed on the cosmic string
tension by using the current Pulsar Timing Array limits on the stochastic
gravitational wave background (SGWB). We have developed a code to compute the
spectrum of gravitational waves (GWs) based on the widely accepted onescale
model. In its simplest form the onescale model allows one to vary: (i) the
string tension, G\mu/c^2; (ii) the size of cosmic string loops relative to the
horizon at birth, \alpha; (iii) the spectral index of the emission spectrum, q;
(iv) the cutoff in the emission spectrum, n_*; and (v) the intercommutation
probability, p. The amplitude and slope of the spectrum in the nHz frequency
range is very sensitive to these unknown parameters. We have also investigated
the impact of more complicated scenarios with multiple initial loop sizes, in
particular the 2\alpha models proposed in the literature and a lognormal
distribution for \alpha. We have computed the constraint on G\mu/c^2 due to the
limit on a SGWB imposed by data from the European Pulsar Timing Array. Taking
into account all the possible uncertainties in the parameters we find a
conservative upper limit of G\mu/c^2<5.3x 10^{7} which typically occurs when
the loop production scale is close to the gravitational backreaction scale,
\alpha\approx\Gamma G\mu/c^2. Stronger limits are possible for specific values
of the parameters which typically correspond to the extremal cases \alpha\ll
\Gamma G\mu/c^2 and \alpha\gg \Gamma G\mu/c^2. This limit is less stringent
than the previously published limits which are based on cusp emission, an
approach which does not necessarily model all the possible uncertainties. We
discuss the prospects for lowering this limit by two orders of magnitude, or
even a detection of the SGWB, in the very near future in the context of the
Large European Array for Pulsars and the Square Kilometre Array.