
Recent nondetection of gravitationalwave backgrounds from pulsar timing
arrays casts further uncertainty on the evolution of supermassive black hole
binaries. We study the capabilities of current gravitationalwave observatories
to detect individual binaries and demonstrate that, contrary to conventional
wisdom, some are in principle detectable throughout the Universe. In
particular, a binary with restframe mass $\gtrsim10^{10}\,M_\odot$ can be
detected by current timing arrays at arbitrarily high redshifts. The same claim
will apply for less massive binaries with more sensitive future arrays. As a
consequence, future searches for nanohertz gravitational waves could be
expanded to target evolving highredshift binaries. We calculate the maximum
distance at which binaries can be observed with pulsar timing arrays and other
detectors, properly accounting for redshift and using realistic binary
waveforms.

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.

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.

In this paper we attempt to investigate the nature of the first gravitational
wave (GW) signal to be detected by pulsar timing arrays (PTAs): will it be an
individual, resolved supermassive black hole binary (SBHB), or a stochastic
background made by the superposition of GWs produced by an ensemble of SBHBs?
To address this issue, we analyse a broad set of simulations of the
cosmological population of SBHBs, that cover the entire parameter space allowed
by current electromagnetic observations in an unbiased way. For each
simulation, we construct the expected GW signal and identify the loudest
individual sources. We then employ appropriate detection statistics to evaluate
the relative probability of detecting each type of source as a function of time
for a variety of PTAs; we consider the current International PTA, and speculate
into the era of the Square Kilometre Array. The main properties of the first
detectable individual SBHBs are also investigated. Contrary to previous work,
we cast our results in terms of the detection probability (DP), since the
commonly adopted criterion based on a signaltonoise ratio threshold is
statisticdependent and may result in misleading conclusions for the statistics
adopted here. Our results confirm quantitatively that a stochastic signal is
more likely to be detected first (with between 75 to 93 per cent probability,
depending on the array), but the DP of singlesources is not negligible. Our
framework is very flexible and can be easily extended to more realistic arrays
and to signal models including environmental coupling and SBHB eccentricity.

Pulsar timing arrays (PTAs) can be used to search for very low frequency
($10^{9}$$10^{7}$ Hz) gravitational waves (GWs). In this paper we present a
general method for the detection and localization of singlesource GWs using
PTAs. We demonstrate the effectiveness of this new method for three types of
signals: monochromatic waves as expected from individual supermassive binary
black holes in circular orbits, GWs from eccentric binaries and GW bursts. We
also test its implementation in realistic data sets that include effects such
as uneven sampling and heterogeneous data spans and measurement precision. It
is shown that our method, which works in the frequency domain, performs as well
as published timedomain methods. In particular, we find it equivalent to the
$\mathcal{F}_{e}$statistic for monochromatic waves. We also discuss the
construction of null streams  data streams that have null response to GWs,
and the prospect of using null streams as a consistency check in the case of
detected GW signals. Finally, we present sensitivities to individual
supermassive binary black holes in eccentric orbits. We find that a
monochromatic search that is designed for circular binaries can efficiently
detect eccentric binaries with both high and low eccentricities, while a
harmonic summing technique provides greater sensitivities only for binaries
with moderate eccentricities.

This paper presents a technique to search for supermassive black hole
binaries (MBHBs) in the Sloan Digital Sky Survey (SDSS). The search is based on
the peculiar properties of merging galaxies that are found in a mock galaxy
catalog from the Millennium Simulation. MBHBs are expected to be the main
gravitational wave (GW) sources for pulsar timing arrays (PTAs); however, it is
still unclear if the observed GW signal will be produced by a few single MBHBs,
or if it will have the properties of a stochastic background. The goal of this
work is to produce a map of the sky in which each galaxy is assigned a
probability of having suffered a recent merger, and of hosting a MBHB that
could be detected by PTAs. This constitutes a step forward in the understanding
of the expected PTA signal: the skymap can be used to investigate the
clustering properties of PTA sources and the spatial distribution of the
observable GW signal power; moreover, galaxies with the highest probabilities
could be used as inputs in targeted searches for individual GW sources. We also
investigate the distribution of neighboring galaxies around galaxies hosting
MBHBs, finding that the most likely detectable PTA sources are located in dense
galaxy environments. Different techniques are used in the search, including
Bayesian and Machine Learning algorithms, with consistent outputs. Our method
generates a list of galaxies classified as MBHB hosts, that can be combined
with other searches to effectively reduce the number of misclassifications. The
spectral coverage of the SDSS reaches less than a fifth of the sky, and the
catalog becomes severely incomplete at large redshifts; however, this technique
can be applied in the future to larger catalogs to obtain complete,
observationallybased information of the expected GW signal detectable by PTAs.

The background of gravitational waves produced by the ensemble of rotating
neutron stars (which includes pulsars, magnetars and gravitars) is
investigated. A formula for \Omega(f) (commonly used to quantify the
background) is derived, properly taking into account the time evolution of the
systems since their formation until the present day. Moreover, the formula
allows one to distinguish the different parts of the background: the
unresolvable (which forms a stochastic background) and the resolvable. Several
estimations of the background are obtained, for different assumptions on the
parameters that characterize neutron stars and their population. In particular,
different initial spin period distributions lead to very different results. For
one of the models, with slow initial spins, the detection of the background can
be rejected. However, other models do predict the detection of the background
by the future groundbased gravitational wave detector ET. A robust upper limit
for the background of rotating neutron stars is obtained; it does not exceed
the detection threshold of two crosscorrelated Advanced LIGO interferometers.
If gravitars exist and constitute more than a few percent of the neutron star
population, then they produce an unresolvable background that could be detected
by ET. Under the most reasonable assumptions on the parameters characterizing a
neutron star, the background is too faint. Previous papers have suggested
neutron star models in which large magnetic fields (like the ones that
characterize magnetars) induce big deformations in the star, which produce a
stronger emission of gravitational radiation. Considering the most optimistic
(in terms of the detection of gravitational waves) of these models, an upper
limit for the background produced by magnetars is obtained; it could be
detected by ET, but not by BBO or DECIGO.

Basic aspects of the background of gravitational waves and its mathematical
characterization are reviewed. The spectral energy density parameter
$\Omega(f)$, commonly used as a quantifier of the background, is derived for an
ensemble of many identical sources emitting at different times and locations.
For such an ensemble, $\Omega(f)$ is generalized to account for the duration of
the signals and of the observation, so that one can distinguish the resolvable
and unresolvable parts of the background. The unresolvable part, often called
confusion noise or stochastic background, is made by signals that cannot be
either individually identified or subtracted out of the data. To account for
the resolvability of the background, the overlap function is introduced. This
function is a generalization of the duty cycle, which has been commonly used in
the literature, in some cases leading to incorrect results. The spectra
produced by binary systems (stellar binaries and massive black hole binaries)
are presented over the frequencies of all existing and planned detectors. A
semianalytical formula for $\Omega(f)$ is derived in the case of stellar
binaries (containing white dwarfs, neutron stars or stellarmass black holes).
Besides a realistic expectation of the level of background, upper and lower
limits are given, to account for the uncertainties in some astrophysical
parameters such as binary coalescence rates. One interesting result concerns
all current and planned groundbased detectors (including the Einstein
Telescope). In their frequency range, the background of binaries is resolvable
and only sporadically present. In other words, there is no stochastic
background of binaries for groundbased detectors.