[Paper Review] Detecting the Stochastic Gravitational Wave Background from Massive Gravity with Pulsar Timing Arrays

Qiuyue Liang, Mark Trodden|arXiv (Cornell University)|Aug 11, 2021

Cosmology and Gravitation TheoriesPhysics and Astronomy28 references36 citations

TL;DR

This paper derives the analytical form of the overlap reduction function (ORF) for the stochastic gravitational wave background (SGWB) in ghost-free massive gravity using pulsar timing arrays (PTAs). It computes corrections to the Hellings-Downs curve due to the graviton mass and five polarization modes (two tensor, two vector, one scalar), providing a framework to test massive gravity against PTA data, with validation in the massless and stationary limits.

ABSTRACT

We explore the potential of Pulsar Timing Arrays (PTAs) such as NANOGrav, EPTA, and PPTA to detect the Stochastic Gravitational Wave Background (SGWB) in theories of massive gravity. In General Relativity, the function describing the dependence of the correlation between the arrival times of signals from two pulsars on the angle between them is known as the Hellings-Downs curve. We compute the analogous overlap reduction function for massive gravity, including the additional polarization states and the correction due to the mass of the graviton, and compare the result with the Hellings-Downs curve. The primary result is a complete analytical form for the analog Hellings-Downs curve, providing a starting point for future numerical studies aimed at a detailed comparison between PTA data and the predictions of massive gravity. We study both the massless limit and the stationary limit as checks on our calculation, and discuss how our formalism also allows us to study the impact of massive spin-2 dark matter candidates on data from PTAs.

Motivation & Objective

• To extend the standard Hellings-Downs curve for PTAs to the case of massive gravity, which includes additional polarization states and a massive graviton.
• To compute the overlap reduction function (ORF) for all five polarization modes (tensor, vector, scalar) in ghost-free massive gravity, accounting for the graviton mass and modified propagation.
• To validate the analytical approximation used in the ORF calculation within the relevant frequency and distance scales of current PTA experiments.
• To provide a model-independent tool for testing massive gravity and related spin-2 dark matter candidates using PTA data.
• To compare the modified ORF with the standard Hellings-Downs curve in General Relativity, highlighting observable differences due to mass and polarization structure.

Proposed method

• Derives polarization tensors for massive spin-2 fields using a formalism analogous to spin-1 polarization vectors, defining five modes: two transverse (h+ and h×), two vector (hx, hy), and one scalar (hl).
• Computes the frequency shift in pulsar signals induced by massive gravitational waves using the metric perturbation hμν, derived from the polarization tensors and the graviton dispersion relation.
• Constructs the overlap reduction function (ORF) as the correlation of signal arrival time shifts between two pulsars, integrating over the wavevector k with a directional dependence on the angle between pulsars.
• Applies the Hellings-Downs approximation by assuming a stochastic, isotropic background and integrating over the wave vector direction, using the standard form for the correlation function.
• Performs numerical checks on the validity of the approximation by comparing the full integral with the approximate form across relevant pulsar distances and PTA frequency bands (1–100 nHz).
• Validates the analytical results in two limits: the massless limit (recovering the standard Hellings-Downs curve) and the stationary limit (relevant for graviton masses near 10−24–10−23 eV).

Experimental results

Research questions

• RQ1 How does the overlap reduction function (ORF) for the stochastic gravitational wave background (SGWB) in massive gravity differ from the standard Hellings-Downs curve in General Relativity?
• RQ2 What are the contributions of the five polarization modes (two tensor, two vector, one scalar) to the ORF in ghost-free massive gravity?
• RQ3 How do the graviton mass and modified dispersion relation affect the shape and amplitude of the ORF in PTA observations?
• RQ4 Is the standard Hellings-Downs approximation valid for massive gravity in the frequency and distance regimes relevant to current PTAs?
• RQ5 Can the formalism be extended to probe massive spin-2 dark matter candidates that may contribute to PTA signals?

Key findings

• The paper derives a complete analytical expression for the overlap reduction function (ORF) in massive gravity, valid for all five polarization modes, which generalizes the Hellings-Downs curve.
• In the massless limit, the tensor-mode ORF reduces to the standard Hellings-Downs curve, confirming consistency with General Relativity.
• The vector and scalar polarization modes contribute non-trivially to the ORF, introducing angular and frequency-dependent deviations from the GR curve.
• The stationary limit of the ORF is analytically computed and shown to be relevant for graviton masses in the range 10−24–10−23 eV, matching the sensitivity window of current PTAs.
• Numerical validation confirms that the key approximation used in the derivation (e.g., the stationary phase or long-wavelength limit) is valid for typical pulsar distances and PTA frequencies (1–100 nHz).
• The formalism provides a direct tool for comparing PTA data (e.g., from NANOGrav, EPTA, PPTA) with predictions from massive gravity and related spin-2 dark matter models.