We discuss the theory of pulsar-timing and astrometry probes of a stochastic gravitational-wave background with a recently developed total-angular-momentum (TAM) formalism for cosmological perturbations. We review the formalism, emphasizing in particular the features relevant for this work and describe the observables we consider (i.e. the pulsar redshift and stellar angular displacement). Using the TAM approach, we calculate the angular power spectra for the observables and from them derive angular auto- and cross-correlation functions. We provide the full set of power spectra and correlation functions not only for the standard transverse-traceless propagating degrees of freedom in general relativity, but also for the four additional non-Einsteinian polarizations that may arise in alternative-gravity theories. We discuss how pulsar-timing and astrometry surveys can complement and serve as cross checks to one another and comment on the importance of testing the chirality of the gravitational-wave background as a tool to understand the nature of its sources. A simple rederivation of the power spectra from the plane-wave formalism is provided in an Appendix.
The detection of a stochastic background of low-frequency gravitational waves by pulsar-timing and astrometric surveys will enable tests of gravitational theories beyond general relativity. These theories generally permit gravitational waves with non-Einsteinian polarization modes, which may propagate slower than the speed of light. We use the total-angular-momentum wave formalism to derive the angular correlation patterns of observables relevant for pulsar timing arrays and astrometry that arise from a background of subluminal gravitational waves with scalar, vector, or tensor polarizations. We find that the pulsar timing observables for the scalar longitudinal mode, which diverge with source distance in the luminal limit, are finite in the subluminal case. Furthermore, we apply our results to $f(R)$ gravity, which contains a massive scalar degree of freedom in addition to the standard transverse-traceless modes. The scalar mode in this $f(R)$ theory is a linear combination of the scalar-longitudinal and scalar-transverse modes, exciting only the monopole and dipole for pulsar timing arrays and only the dipole for astrometric surveys.
We have begun an exciting era for gravitational wave detection, as several world-leading experiments are breaching the threshold of anticipated signal strengths. Pulsar timing arrays (PTAs) are pan-Galactic gravitational wave detectors that are already cutting into the expected strength of gravitational waves from cosmic strings and binary supermassive black holes in the nHz-$mu$Hz gravitational wave band. These limits are leading to constraints on the evolutionary state of the Universe. Here, we provide a broad review of this field, from how pulsars are used as tools for detection, to astrophysical sources of uncertainty in the signals PTAs aim to see, to the primary current challenge areas for PTA work. This review aims to provide an up-to-date reference point for new parties interested in the field of gravitational wave detection via pulsar timing.
In this paper, we focus on testing gravity theories in the radiative regime using pulsar timing array observations. After reviewing current techniques to measure the dispersion and alternative polarization of gravitational waves, we extend the framework to the most general situations, where the combinations of a massive graviton and alternative polarization modes are considered. The atlas of the Hellings-Downs functions is completed by the new calculations for these dispersive alternative polarization modes. We find that each mode and corresponding graviton mass introduce characteristic features in the Hellings-Downs function. Thus, in principal, we can not only detect each polarization mode, measure the corresponding graviton mass, but also discriminate the different scenarios. In this way, we can test gravity theories in the radiative regime in a generalized fashion, and such method is a direct experiment, where one can address the gauge symmetry of the gravity theories in their linearised limits. Although current pulsar timing still lacks enough stable pulsars and sensitivity for such practices, we expect that future telescopes with larger collecting area could make such experiments be feasible.
Within the next several years, pulsar-timing array programs will likely usher in the next era of gravitational-wave astronomy through the detection of a stochastic background of nanohertz-frequency gravitational waves, originating from a cosmological population of inspiraling supermassive binary black holes. While the source positions will likely be isotropic to a good approximation, the gravitational-wave angular power distribution will be anisotropic, with the most massive and/or nearby binaries producing signals that may resound above the background. We study such a realistic angular power distribution, developing fast and accurate sky-mapping strategies to localize pixels and extended regions of excess power while simultaneously modeling the background signal from the less massive and more distant ensemble. We find that power anisotropy will be challenging to discriminate from isotropy for realistic gravitational-wave skies, requiring SNR $>10$ in order to favor anisotropy with $10:1$ posterior odds in our case study. Amongst our techniques, modeling the population signal with multiple point sources in addition to an isotropic background provides the most physically-motivated and easily interpreted maps, while spherical-harmonic modeling of the square-root power distribution, $P(hatOmega)^{1/2}$, performs best in discriminating from overall isotropy. Our techniques are modular and easily incorporated into existing pulsar-timing array analysis pipelines.
The maximum frequency of gravitational waves (GWs) detectable with traditional pulsar timing methods is set by the Nyquist frequency ($f_{rm{Ny}}$) of the observation. Beyond this frequency, GWs leave no temporal-correlated signals; instead, they appear as white noise in the timing residuals. The variance of the GW-induced white noise is a function of the position of the pulsars relative to the GW source. By observing this unique functional form in the timing data, we propose that we can detect GWs of frequency $>$ $f_{rm{Ny}}$ (super-Nyquist frequency GWs; SNFGWs). We demonstrate the feasibility of the proposed method with simulated timing data. Using a selected dataset from the Parkes Pulsar Timing Array data release 1 and the North American Nanohertz Observatory for Gravitational Waves publicly available datasets, we try to detect the signals from single SNFGW sources. The result is consistent with no GW detection with 65.5% probability. An all-sky map of the sensitivity of the selected pulsar timing array to single SNFGW sources is generated, and the position of the GW source where the selected pulsar timing array is most sensitive to is $lambda_{rm{s}}=-0.82$, $beta_{rm{s}}=-1.03$ (rad); the corresponding minimum GW strain is $h=6.31times10^{-11}$ at $f=1times10^{-5}$ Hz.