
We present highprecision timing data over time spans of up to 11 years for
45 millisecond pulsars observed as part of the North American Nanohertz
Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and
characterizing lowfrequency gravitational waves. The pulsars were observed
with the Arecibo Observatory and/or the Green Bank Telescope at frequencies
ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately
monthly cadence, with six hightimingprecision pulsars observed weekly, and
all were observed at widely separated frequencies at each observing epoch in
order to fit for timevariable dispersion delays. We describe our methods for
data processing, timeofarrival (TOA) calculation, and the implementation of a
new, automated method for removing outlier TOAs. We fit a timing model for each
pulsar that includes spin, astrometric, and, if necessary, binary parameters,
in addition to timevariable dispersion delays and parameters that quantify
pulseprofile evolution with frequency. The new timing solutions provide three
new parallax measurements, two new Shapiro delay measurements, and two new
measurements of large orbitalperiod variations. We fit models that
characterize sources of noise for each pulsar. We find that 11 pulsars show
significant red noise, with generally smaller spectral indices than typically
measured for nonrecycled pulsars, possibly suggesting a different origin.
Future papers will use these data to constrain or detect the signatures of
gravitationalwave signals.

We compute upper limits on the nanohertzfrequency isotropic stochastic
gravitational wave background (GWB) using the 9year 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 powerlaw analysis we construct posterior
distributions on environmental factors that drive the binary to the GWdriven
regime including the stellar mass density for stellarscattering, 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.

The 2012 International Pulsar Timing Array (IPTA) Mock Data Challenge (MDC)
is designed to test current Gravitational Wave (GW) detection algorithms. Here
we will briefly outline two detection algorithms for a stochastic background of
gravitational waves, namely, a firstorder likelihood method and an optimal
statistic method and present our results from the closed MDC data sets.