
Symmetries play an important role in modern theories of gravity. The strong
equivalence principle (SEP) constitutes a collection of gravitational
symmetries which are all implemented by general relativity. Alternative
theories, however, are generally expected to violate some aspects of SEP. We
test three aspects of SEP using observed change rates in the orbital period and
eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's
constancy as part of locational invariance of gravitation; 2. the
postNewtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance;
3. the universality of free fall (UFF) for strongly selfgravitating bodies.
Based on the pulsar timing result of the combined dataset from the North
American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar
Timing Array (EPTA), we find $\dot{G}/G = (0.1 \pm 0.9) \times 10^{12}\,{\rm
yr}^{1}$, which is weaker than Solar system limits, but applies for strongly
selfgravitating objects. Furthermore, we obtain the constraints $\Delta<
0.002$ for the UFF test and $3\times10^{20} < \hat{\alpha}_3 <
4\times10^{20}$ at 95% confidence. These are the first direct UFF and
$\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various
limitations of previous tests.

We search for an isotropic stochastic gravitationalwave background (GWB) in
the newly released $11$year dataset from the North American Nanohertz
Observatory for Gravitational Waves (NANOGrav). While we find no significant
evidence for a GWB, we place constraints on a GWB from a population of
supermassive blackhole binaries, from cosmic strings, and from a primordial
GWB. For the first time, we find that the GWB upper limits and detection
statistics are sensitive to the Solar System ephemeris (SSE) model used, and
that SSE errors can mimic a GWB signal. To mitigate this effect, we developed
and implemented a novel approach that bridges systematic SSE differences,
producing the first PTA constraints that are robust against SSE uncertainties.
We place a $95\%$ upper limit on the GW strain amplitude of
$A_\mathrm{GWB}<1.45\times 10^{15}$ at a frequency of $f=1$ yr$^{1}$ for a
fiducial $f^{2/3}$ powerlaw spectrum, and with interpulsar correlations
modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$year
limit, calculated using the same procedure. Previous PTA upper limits on the
GWB will need revision in light of SSE systematic uncertainties. We also
characterize the combined influence of the massdensity of stars in galactic
cores, the eccentricity of binaries at formation, and the relation between the
mass of the central supermassive black hole and the galactic bulge (the
$M_\mathrm{BH}M_\mathrm{bulge}$ relation). We constrain cosmicstring tension
on the basis of recent simulations, yielding an SSEmarginalized 95\% upper
limit on the cosmic string tension of $G\mu < 5.3\times 10^{11}$a factor of
$\sim 2$ better than the NANOGrav $9$year constraints. We then use our new
Bayesian SSE model to limit the energy density of primordial GWBs,
corresponding to $\Omega_\mathrm{GWB}(f)h^2<3.4 \times 10^{10}$ for a
radiationdominated inflationary era. [ABRIDGED]

We analyze dispersion measure (DM) variations of 37 millisecond pulsars in
the 9year NANOGrav data release and constrain the sources of these variations.
Variations are significant for nearly all pulsars, with characteristic
timescales comparable to or even shorter than the average spacing between
observations. Five pulsars have periodic annual variations, 14 pulsars have
monotonically increasing or decreasing trends, and 13 pulsars show both
effects. Several pulsars show correlations between DM excesses and lines of
sight that pass close to the Sun. Mapping of the DM variations as a function of
the pulsar trajectory can identify localized ISM features and, in one case, an
upper limit to the size of the dispersing region of 13.2 AU. Finally, five
pulsars show very nearly quadratic structure functions, which could be
indicative of an underlying Kolmogorov medium. Four pulsars show roughly
Kolmogorov structure functions and another four show structure functions less
steep than Kolmogorov. One pulsar has too large an uncertainty to allow
comparisons. We discuss explanations for apparent departures from a
Kolmogorovlike spectrum, and show that the presence of other trends in the
data is the most likely cause.

Gravitational wave astronomy using a pulsar timing array requires
highquality millisecond pulsars, correctable interstellar propagation delays,
and highprecision measurements of pulse times of arrival. Here we identify
noise in timing residuals that exceeds that predicted for arrival time
estimation for millisecond pulsars observed by the North American Nanohertz
Observatory for Gravitational Waves. We characterize the excess noise using
variance and structure function analyses. We find that 26 out of 37 pulsars
show inconsistencies with a whitenoiseonly model based on the short timescale
analysis of each pulsar and we demonstrate that the excess noise has a red
power spectrum for 15 pulsars. We also decompose the excess noise into
chromatic (radiofrequencydependent) and achromatic components. Associating
the achromatic rednoise component with spin noise and including additional
powerspectrumbased estimates from the literature, we estimate a scaling law
in terms of spin parameters (frequency and frequency derivative) and dataspan
length and compare it to the scaling law of Shannon \& Cordes (2010). We
briefly discuss our results in terms of detection of gravitational waves at
nanohertz frequencies.

The use of pulsars as astrophysical clocks for gravitational wave experiments
demands the highest possible timing precision. Pulse times of arrival (TOAs)
are limited by stochastic processes that occur in the pulsar itself, along the
line of sight through the interstellar medium, and in the measurement process.
On timescales of seconds to hours, the TOA variance exceeds that from
templatefitting errors due to additive noise. We assess contributions to the
total variance from two additional effects: amplitude and phase jitter
intrinsic to single pulses and changes in the interstellar impulse response
from scattering. The three effects have different dependencies on time,
frequency, and pulse signaltonoise ratio. We use data on 37 pulsars from the
North American Nanohertz Observatory for Gravitational Waves to assess the
individual contributions to the overall intraday noise budget for each pulsar.
We detect jitter in 22 pulsars and estimate the average value of RMS jitter in
our pulsars to be $\sim 1\%$ of pulse phase. We examine how jitter evolves as a
function of frequency and find evidence for evolution. Finally, we compare our
measurements with previous noise parameter estimates and discuss methods to
improve gravitational wave detection pipelines.

Dense, continuous pulsar timing observations over a 24hr period provide a
method for probing intermediate gravitational wave (GW) frequencies from 10
microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North
American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes
Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array
(IPTA) all use millisecond pulsar observations to detect or constrain GWs
typically at nanohertz frequencies. In the case of the IPTA's ninetelescope
24Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the
intermediate frequency regime can be produced. The negligible change in
dispersion measure during the observation minimizes red noise in the timing
residuals, constraining any contributions from GWs due to individual sources.
At 10$^{5}$Hz, the 95% upper limit on strain is 10$^{11}$ for GW sources in
the pulsar's direction.

PSR J1024$$0719 is a millisecond pulsar that was long thought to be
isolated. However, puzzling results concerning its velocity, distance, and low
rotational period derivative have led to reexamination of its properties. We
present updated radio timing observations along with new and archival optical
data that show PSR J1024$$0719 is most likely in a long period (2$$20 kyr)
binary system with a lowmass ($\approx 0.4\,M_\odot$) lowmetallicity ($Z
\approx 0.9\,$ dex) main sequence star. Such a system can explain most of the
anomalous properties of this pulsar. We suggest that this system formed through
a dynamical exchange in a globular cluster that ejected it into a halo orbit,
consistent with the low observed metallicity for the stellar companion. Further
astrometric and radio timing observations such as measurement of the third
period derivative could strongly constrain the range of orbital parameters.

We present highprecision timing observations spanning up to nine years for
37 millisecond pulsars monitored with the Green Bank and Arecibo radio
telescopes as part of the North American Nanohertz Observatory for
Gravitational Waves (NANOGrav) project. We describe the observational and
instrumental setups used to collect the data, and methodology applied for
calculating pulse times of arrival; these include novel methods for measuring
instrumental offsets and characterizing low signaltonoise ratio timing
results. The time of arrival data are fit to a physical timing model for each
source, including terms that characterize timevariable dispersion measure and
frequencydependent pulse shape evolution. In conjunction with the timing model
fit, we have performed a Bayesian analysis of a parameterized timing noise
model for each source, and detect evidence for excess lowfrequency, or "red,"
timing noise in 10 of the pulsars. For 5 of these cases this is likely due to
interstellar medium propagation effects rather than intrisic spin variations.
Subsequent papers in this series will present further analysis of this data set
aimed at detecting or limiting the presence of nanohertzfrequency
gravitational wave signals.

We report 21yr timing of one of the most precise pulsars: PSR J1713+0747.
Its pulse times of arrival are well modeled by a comprehensive pulsar binary
model including its threedimensional orbit and a noise model that incorporates
correlated noise such as jitter and red noise. Its timing residuals have
weighted root mean square $\sim 92$ ns. The new dataset allows us to update and
improve previous measurements of the system properties, including the masses of
the neutron star ($1.31\pm0.11$ $M_{\odot}$) and the companion white dwarf
($0.286\pm0.012$ $M_{\odot}$) and the parallax distance $1.15\pm0.03$ kpc. We
measured the intrinsic change in orbital period, $\dot{P}^{\rm Int}_{\rm b}$,
is $0.20\pm0.17$ ps s$^{1}$, which is not distinguishable from zero. This
result, combined with the measured $\dot{P}^{\rm Int}_{\rm b}$ of other
pulsars, can place a generic limit on potential changes in the gravitational
constant $G$. We found that $\dot{G}/G$ is consistent with zero
[$(0.6\pm1.1)\times10^{12}$ yr$^{1}$, 95\% confidence] and changes at least
a factor of $31$ (99.7\% confidence) more slowly than the average expansion
rate of the Universe. This is the best $\dot{G}/G$ limit from pulsar binary
systems. The $\dot{P}^{\rm Int}_{\rm b}$ of pulsar binaries can also place
limits on the putative coupling constant for dipole gravitational radiation
$\kappa_D=(0.9\pm3.3)\times10^{4}$ (95\% confidence). Finally, the nearly
circular orbit of this pulsar binary allows us to constrain statistically the
strongfield postNewtonian parameters $\Delta$, which describes the violation
of strong equivalence principle, and $\hat{\alpha}_3$, which describes a
breaking of both Lorentz invariance in gravitation and conservation of
momentum. We found, at 95\% confidence, $\Delta<0.01$ and
$\hat{\alpha}_3<2\times10^{20}$ based on PSR J1713+0747.

The radio millisecond pulsar J1713+0747 is regarded as one of the
highestprecision clocks in the sky, and is regularly timed for the purpose of
detecting gravitational waves. The International Pulsar Timing Array
collaboration undertook a 24hour global observation of PSR J1713+0747 in an
effort to better quantify sources of timing noise in this pulsar, particularly
on intermediate (1  24 hr) timescales. We observed the pulsar continuously
over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell,
Nancay, Parkes, and WSRT radio telescopes. The combined pulse timesofarrival
presented here provide an estimate of what sources of timing noise, excluding
DM variations, would be present as compared to an idealized rootN improvement
in timing precision, where N is the number of pulses analyzed. In the case of
this particular pulsar, we find that intrinsic pulse phase jitter dominates
arrival time precision when the S/N of single pulses exceeds unity, as measured
using the eight telescopes that observed at Lband/1.4 GHz. We present first
results of specific phenomena probed on the unusually long timescale (for a
single continuous observing session) of tens of hours, in particular
interstellar scintillation, and discuss the degree to which scintillation and
profile evolution affect precision timing. This paper presents the data set as
a basis for future, deeper studies.