TEMPO2, a new pulsar timing package. III: Gravitational wave simulation

Analysis of pulsar timing data-sets may provide the first direct detection of gravitational waves. This paper, the third in a series describing the mathematical framework implemented into the tempo2 pulsar timing package, reports on using tempo2 to simulate the timing residuals induced by gravitational waves. The tempo2 simulations can be used to provide upper bounds on the amplitude of an isotropic, stochastic, gravitational wave background in our Galaxy and to determine the sensitivity of a given pulsar timing experiment to individual, supermassive, binary black hole systems.

💡 Research Summary

The paper presents a comprehensive implementation of gravitational‑wave (GW) simulation capabilities within the TEMPO2 pulsar‑timing software, extending the mathematical framework described in the earlier two papers of the series. The authors begin by reformulating the standard pulsar timing model to incorporate the effect of a passing GW on the propagation time of radio pulses. By linearising the metric perturbation h_{ij}(t, x) and projecting it onto the line‑of‑sight vector between the pulsar and the Earth, they derive an explicit expression for the induced timing residual as an integral over the GW strain along the photon trajectory. This formulation naturally separates the two polarisation modes (plus and cross) and introduces the GW propagation direction as a set of free parameters.

Two principal simulation scenarios are explored. The first concerns an isotropic, stochastic GW background, typically modeled as the superposition of many super‑massive black‑hole binaries, with a characteristic strain spectrum h_c(f)∝f^{-2/3}. Using TEMPO2, the same synthetic background is injected into the timing models of a network of pulsars. The resulting residuals are cross‑correlated and compared with the Hellings‑Downs angular correlation curve, allowing the authors to place upper limits on the dimensionless energy density Ω_gw. Their injection‑recovery tests show that current PTA datasets (e.g., the Parkes Pulsar Timing Array) can constrain Ω_gw to ≲10^{-9} in the nanohertz band.

The second scenario targets an individual continuous‑wave source, such as a binary super‑massive black hole with masses of order 10^9 M_⊙, orbital period of a few years, and located within a few hundred megaparsecs. The source is fully parameterised by its chirp mass, distance, sky location, orbital phase, inclination, and polarisation angle. TEMPO2 computes the corresponding deterministic residual waveform for each pulsar in the array, enabling the construction of sensitivity curves that map detectable strain amplitudes as a function of source sky position and frequency. The authors demonstrate that, for optimistic source parameters, a 5‑σ detection is achievable with existing timing precision (tens of nanoseconds) and a decade‑long data span.

From a software engineering perspective, the GW module is delivered as a plug‑in that reads a simple text‑based “.gw” configuration file. Users can specify either a stochastic spectrum (via amplitude and spectral index) or a deterministic binary source (via orbital elements). The command “add_gw” automatically augments the standard TEMPO2 timing model, and the resulting residuals are written in the familiar .tim format, ensuring seamless integration with existing analysis pipelines. Parallelisation options allow simultaneous processing of dozens of pulsars, making large‑scale injection campaigns computationally tractable.

Statistical analysis of the simulated data is performed using both frequentist (maximum‑likelihood) and Bayesian (Markov‑Chain Monte Carlo, evidence‑ratio) techniques. The authors illustrate how posterior distributions for Ω_gw or binary parameters can be recovered from the injected data, and they quantify detection efficiencies as a function of timing noise, cadence, and data span. The study highlights that the dominant limitation for current PTAs is the long‑term stability of pulsar clocks (≈10–30 ns over ten years) rather than instrumental noise.

In conclusion, the paper establishes TEMPO2 as a unified platform for both real‑world pulsar timing analysis and realistic GW signal simulation. This capability enables researchers to design optimal PTA experiments, forecast sensitivity improvements for upcoming facilities such as the Square Kilometre Array, and ultimately to extract or constrain nanohertz gravitational‑wave signals from the stochastic background and from individual super‑massive black‑hole binaries.