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 ultra-precise timing data on 41
millisecond pulsars. We develop frequentist and Bayesian detection algorithms
to search both for monochromatic and frequency-evolving 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
sky-averaged 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
sub-centiparsec 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 state-of-the-art SMBHB
population models predict $<1\%$ probability of detecting a CGW with the
current EPTA dataset, consistent with the reported non-detection. 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 compute upper limits on the nanohertz-frequency isotropic stochastic
gravitational wave background (GWB) using the 9-year data release from the
North American Nanohertz Observatory for Gravitational Waves (NANOGrav)
collaboration. We set upper limits for a GWB from supermassive black hole
binaries under power law, broken power law, and free spectral coefficient GW
spectrum models. We place a 95\% upper limit on the strain amplitude (at a
frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times
10^{-15}$. For a broken power law model, we place priors on the strain
amplitude derived from simulations of Sesana (2013) and McWilliams et al.
(2014). We find that the data favor a broken power law to a pure power law with
odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models,
respectively. The McWilliams model is essentially ruled out by the data, and
the Sesana model is in tension with the data under the assumption of a pure
power law. Using the broken power-law analysis we construct posterior
distributions on environmental factors that drive the binary to the GW-driven
regime including the stellar mass density for stellar-scattering, mass
accretion rate for circumbinary disk interaction, and orbital eccentricity for
eccentric binaries, marking the first time that the shape of the GWB spectrum
has been used to make astrophysical inferences. We then place the most
stringent limits so far on the energy density of relic GWs,
$\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the
Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where
$m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB,
$\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a
conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the
joint Planck and high-$l$ CMB data from other experiments.