
Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion
$M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large
and poorly constrained concentrations of gas and stars. There are currently no
observations of merging SMBHBs it is in fact possible that they stall at
their final parsec of separation and never merge. While LIGO has detected high
frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is
inaccessible to groundbased interferometers, but possible with Pulsar Timing
Arrays (PTAs). Using data from local galaxies in the 2 Micron AllSky Survey,
together with galaxy merger rates from Illustris, we find that there are on
average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries
which will never merge. Local unresolved SMBHBs can contribute to GW background
anisotropy at a level of $\sim20\%$, and if the GW background can be
successfully isolated, GWs from at least one local SMBHB can be detected in 10
years.

Quantum fluctuations of the gravitational field in the early Universe,
amplified by inflation, produce a primordial gravitationalwave background
across a broad frequency band. We derive constraints on the spectrum of this
gravitational radiation, and hence on theories of the early Universe, by
combining experiments that cover 29 orders of magnitude in frequency. These
include Planck observations of cosmic microwave background temperature and
polarization power spectra and lensing, together with baryon acoustic
oscillations and big bang nucleosynthesis measurements, as well as new pulsar
timing array and groundbased interferometer limits. While individual
experiments constrain the gravitationalwave energy density in specific
frequency bands, the combination of experiments allows us to constrain
cosmological parameters, including the inflationary spectral index, $n_t$, and
the tensortoscalar ratio, $r$. Results from individual experiments include
the most stringent nanohertz limit of the primordial background to date from
the Parkes Pulsar Timing Array, $\Omega_{\rm gw}(f)<2.3\times10^{10}$.
Observations of the cosmic microwave background alone limit the
gravitationalwave spectral index at 95\% confidence to $n_t\lesssim5$ for a
tensortoscalar ratio of $r = 0.11$. However, the combination of all the above
experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected
to further constrain $n_t<0.34$ by 2020. When cosmic microwave background
experiments detect a nonzero $r$, our results will imply even more stringent
constraints on $n_t$ and hence theories of the early Universe.

Most black holes (BHs) will absorb a neutron star (NS) companion fully
intact, without tidal disruption, suggesting the pair will remain dark to
telescopes. Even without tidal disruption, electromagnetic luminosity is
generated from the battery phase of the binary when the BH interacts with the
NS magnetic field. Originally the luminosity was expected in highenergy Xrays
or gammarays, however we conjecture that some of the battery power is emitted
in the radio bandwidth. While the luminosity and timescale are suggestive of
fast radio bursts (FRBs; millisecondscale radio transients) NSBH coalescence
rates are too low to make these a primary FRB source. Instead, we propose the
transients form a FRB subpopulation, distinguishable by a double peak with a
precursor. The rapid rampup in luminosity manifests as a precursor to the
burst which is $20\%80\%$ as luminous, given 0.5~ms timing resolution. The
main burst is from the peak luminosity before merger. The postmerger burst
follows from the NS magnetic field migration to the BH, causing a shock. NSBH
pairs are especially desirable for groundbased gravitational wave (GW)
observatories since the pair might not otherwise be detected, with
electromagnetic counterparts greatly augmenting the scientific leverage beyond
the GW signal. Valuably, the electromagnetic signal can break degeneracies in
the parameters encoded in the GW as well as probe the NS magnetic field
strength, yielding insights into open problems in NS magnetic field decay.

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 ultraprecise timing data on 41
millisecond pulsars. We develop frequentist and Bayesian detection algorithms
to search both for monochromatic and frequencyevolving 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
skyaveraged 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
subcentiparsec 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 stateoftheart SMBHB
population models predict $<1\%$ probability of detecting a CGW with the
current EPTA dataset, consistent with the reported nondetection. 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 gravitationalwave
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 supermassive 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
NambuGoto 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.

We extend the formalisms developed in Gair et al. and Cornish and van
Haasteren to create maps of gravitationalwave backgrounds using a network of
groundbased laser interferometers. We show that in contrast to pulsar timing
arrays, which are insensitive to half of the gravitationalwave sky (the curl
modes), a network of groundbased interferometers is sensitive to both the
gradient and curl components of the background. The spatial separation of a
network of interferometers, or of a single interferometer at different times
during its rotational and orbital motion around the Sun, allows for recovery of
both components. We derive expressions for the response functions of a laser
interferometer in the smallantenna limit, and use these expressions to
calculate the overlap reduction function for a pair of interferometers. We also
construct maximumlikelihood estimates of the + and xpolarization modes of the
gravitationalwave sky in terms of the response matrix for a network of
groundbased interferometers, evaluated at discrete times during Earth's
rotational and orbital motion around the Sun. We demonstrate the feasibility of
this approach for some simple simulated backgrounds (a single point source and
spatiallyextended distributions having only grad or curl components),
calculating maximumlikelihood sky maps and uncertainty maps based on the
(pseudo)inverse of the response matrix. The distinction between this approach
and standard methods for mapping gravitationalwave power is also discussed.

We describe an alternative approach to the analysis of gravitationalwave
backgrounds, based on the formalism used to characterise the polarisation of
the cosmic microwave background. In contrast to standard analyses, this
approach makes no assumptions about the nature of the background and so has the
potential to reveal much more about the physical processes that generated it.
An arbitrary background can be decomposed into modes whose angular dependence
on the sky is given by gradients and curls of spherical harmonics. We derive
the pulsar timing overlap reduction functions for the individual modes, which
are given by simple combinations of spherical harmonics evaluated at the pulsar
locations. We show how these can be used to recover the components of an
arbitrary background, giving explicit results for both isotropic and
anisotropic uncorrelated backgrounds. We also find that the response of a
pulsar timing array to curl modes is identically zero, so half of the
gravitationalwave sky will never be observed using pulsar timing, no matter
how many pulsars are included in the array. An isotropic, unpolarised and
uncorrelated background can be accurately represented using only three modes,
and so a search of this type will be only slightly more complicated than the
standard crosscorrelation search using the Hellings and Downs overlap
reduction function. However, by measuring the components of individual modes of
the background and checking for consistency with isotropy, this approach has
the potential to reveal much more information. Each individual mode on its own
describes a background that is correlated between different points on the sky.
A measurement of the components that indicates the presence of correlations in
the background on large angular scales would suggest startling new physics.

One of the primary objectives for Pulsar Timing Arrays (PTAs) is to detect a
stochastic background generated by the incoherent superposition of
gravitational waves (GWs), in particular from the cosmic population of
supermassive black hole binaries. Current stochastic background searches assume
that pulsars in a PTA are separated from each other and the Earth by many GW
wavelengths. As more millisecond pulsars are discovered and added to PTAs, some
may be separated by only a few radiation wavelengths or less, resulting in
correlated GW phase changes between close pulsars in the array. Here we
investigate how PTA overlap reduction functions (ORFs), up to quadrupole order,
are affected by these additional correlated phase changes, and how they are in
turn affected by relaxing the assumption that all pulsars are equidistant from
the solar system barycenter. We find that in the low frequency GW background
limit of $f\sim10^{9}$~Hz, and for pulsars at varying distances from the
Earth, that these additional correlations only affect the ORFs by a few percent
for pulsar pairs at large angular separations, as expected. However when nearby
(order 100 pc) pulsars are separated by less than a few degrees, the correlated
phase changes can introduce variations of a few tens of percent in the
magnitude of the isotropic ORF, and much larger fractional differences in the
anisotropic ORFs up to 188 in the $m=0$, $l=2$ ORF for equidistant pulsars
separated by 3 degrees. In fact, the magnitude of most of the anisotropic ORFs
is largest at small, but nonzero, pulsar separations. Finally, we write down a
small angle approximation for the correlated phase changes which can easily be
implemented in search pipelines, and for completeness, examine the behavior of
the ORFs for pulsars which lie at a radiation wavelength from the Earth.

Detecting a stochastic gravitational wave background, particularly radiation
from individually unresolvable supermassive black hole binary systems, is one
of the primary targets for Pulsar Timing Arrays. Increasingly more stringent
upper limits are being set on these signals under the assumption that the
background radiation is isotropic. However, some level of anisotropy may be
present and the characterisation of the power at different angular scales
carries important information. We show that the standard analysis for isotropic
backgrounds can be generalised in a conceptually straightforward way to the
case of generic anisotropic background radiation by decomposing the angular
distribution of the gravitational wave power on the sky into multipole moments.
We introduce the concept of generalised overlap reduction functions which
characterise the effect of the anisotropy multipoles on the correlation of the
timing residuals from the pulsars timed by a Pulsar Timing Array. In a search
for a signal characterised by a generic anisotropy, the generalised overlap
reduction functions play the role of the socalled Hellings and Downs curve
used for isotropic radiation. We compute the generalised overlap reduction
functions for a generic level of anisotropy and Pulsar Timing Array
configuration. We also provide an order of magnitude estimate of the level of
anisotropy that can be expected in the background generated by supermassive
black hole binary systems.

This is a summary of the methods we used to analyse the first IPTA Mock Data
Challenge (MDC), and the obtained results. We have used a Bayesian analysis in
the time domain, accelerated using the recently developed ABCmethod which
consists of a form of lossy linear data compression. The TOAs were first
processed with Tempo2, where the design matrix was extracted for use in a
subsequent Bayesian analysis. We used different noise models to analyse the
datasets: no red noise, red noise the same for all pulsars, and individual red
noise per pulsar. We sampled from the likelihood with four different samplers:
"emcee", "twalk", "MetropolisHastings", and "pyMultiNest". All but emcee
agreed on the final result, with emcee failing due to artefacts of the
highdimensionality of the problem. An interesting issue we ran into was that
the prior of all the 36 (red) noise amplitudes strongly affects the results. A
flat prior in the noise amplitude biases the inferred GWB amplitude, whereas a
flat prior in logamplitude seems to work well. This issue is only apparent
when using a noise model with individually modelled red noise for all pulsars.
Our results for the blind challenges are in good agreement with the injected
values. For the GWB amplitudes we found h_c = 1.03 +/ 0.11 [10^{14}], h_c =
5.70 +/ 0.35 [10^{14}], and h_c = 6.91 +/ 1.72 [10^{15}], and for the GWB
spectral index we found gamma = 4.28 +/ 0.20, gamma = 4.35 +/ 0.09, and gamma
= 3.75 +/ 0.40. We note that for closed challenge 3 there was quite some
covariance between the signal and the red noise: if we constrain the GWB
spectral index to the usual choice of gamma = 13/3, we obtain the estimates:
h_c = 10.0 +/ 0.64 [10^{15}], h_c = 56.3 +/ 2.42 [10^{15}], and h_c = 4.83
+/ 0.50 [10^{15}], with onesided 2 sigma upperlimits of: h_c <= 10.98
[10^{15}], h_c <= 60.29 [10^{15}], and h_c <= 5.65 [10^{15}].