Pulsar timing arrays are sensitive to gravitational wave perturbations produced by individual supermassive black hole binaries during their early inspiral phase. Modified gravity theories allow for the emission of gravitational dipole radiation, which is enhanced relative to the quadrupole contribution for low orbital velocities, making the early inspiral an ideal regime to test for the presence of modified gravity effects. Using a theory-agnostic description of modified gravity theories based on the parametrized post-Einsteinian framework, we explore the possibility of detecting deviations from General Relativity using simulated pulsar timing array data, and provide forecasts for the constraints that can be achieved. We generalize the {\tt enterprise} pulsar timing software to account for possible additional polarization states and modifications to the phase evolution, and study how accurately the parameters of simulated signals can be recovered. We find that while a pure dipole model can partially recover a pure quadrupole signal, there is little possibility for confusion when the full model with all polarization states is used. With no signal present, and using noise levels comparable to those seen in contemporary arrays, we produce forecasts for the upper limits that can be placed on the amplitudes of alternative polarization modes as a function of the sky location of the source.

Gravitational waves perturb the paths of photons, impacting both the time-of-flight and the arrival direction of light from stars. Pulsar timing arrays can detect gravitational waves by measuring the variations in the time of flight of radio pulses, while astrometry missions such as Gaia can detect gravitational waves from the time-varying changes in the apparent position of a field of stars. Just as gravitational waves impart a characteristic correlation pattern in the arrival times of pulses from pulsars at different sky locations, the deflection of starlight is similarly correlated across the sky. Here we compute the astrometric correlation patterns for the full range of polarization states found in alternative theories of gravity, and decompose the sky-averaged correlation patterns into vector spherical harmonics. We find that the tensor and vector polarization states produce equal power in the electric- and magnetic-type vector spherical harmonics, while the scalar modes produce only electric-type correlations. Any difference in the measured electric and magnetic-type correlations would represent a clear violation of Einstein gravity. The angular correlations functions for the vector and scalar longitudinal modes show the same enhanced response at small angular separations that is familiar from pulsar timing.