Testing General Relativity with Individual Supermassive Black Hole Binaries

TL;DR

This paper develops a comprehensive framework to test deviations from General Relativity using continuous gravitational waves from supermassive black hole binaries, enhancing the potential for precision gravity tests with pulsar timing arrays.

Contribution

It introduces a unified method to analyze various beyond-GR effects on gravitational wave signals, including polarization, dispersion, and birefringence, validated through simulations.

Findings

PTAs can detect non-tensorial polarizations with linear scaling.

Modified dispersion relations are effectively constrained at nanohertz frequencies.

Birefringence effects are generally suppressed at the studied frequencies.

Abstract

We develop a unified framework for testing gravity beyond General Relativity (GR) with continuous gravitational waves (CWs) from individual supermassive black hole binaries (SMBHBs). These long-lived, nearly monochromatic nanohertz signals offer unique strengths for precision tests of gravity, since their coherent phase evolution and inter-pulsar correlations in pulsar timing arrays (PTAs) retain detailed information about departures from GR over cosmological propagation distances. We consider three representative classes of deviations from GR: additional polarization states, modified dispersion relations, and parity-violating birefringence. For each, we derive the inter-pulsar cross correlation, the modified antenna response, and the propagation-induced pulsar-term phase delay. For non-tensorial polarizations, the CW cross correlation scales linearly in the alternative-polarization amplitude, compared to the quadratic scaling of the gravitational-wave background (GWB), provided the beyond-GR modes are sub-dominant. PTAs are also competitive for modified dispersion relations, where low frequencies enhance both the antenna-pattern modification and the pulsar-term phase delay. Birefringence, by contrast, is suppressed at nanohertz frequencies for most parity-violating theories. We validate the framework with injection-and-recovery simulations for breathing-mode and massive-graviton signals at current observational limits, recovering the injected beyond-GR parameters and distinguishing the CW signal from both correlated and uncorrelated background models. We further show that a pure-GR CW template recovers source parameters without significant bias when beyond-GR physics is present in the data, supporting a two-stage analysis strategy: identify candidates under GR, then test for deviations.