Spatial correlation between pulsars from interfering gravitational-wave sources in massive gravity

In the nanohertz band, the spatial correlations in pulsar timing arrays (PTAs) produced by interfering gravitational waves (GWs) from multiple sources likely deviate from the traditional ones without interference under the assumption of an isotropic Gaussian ensemble. This work investigates the impact of such interference within the framework of massive gravity. Through simulations, we show that while the resulting correlation patterns can be described by Legendre expansions with coefficients that depend on the interference configuration, they remain predominantly quadrupolar (l = 2), with this feature becoming more pronounced as the graviton mass increases--reflecting both the tensorial polarizations and the modified GW dispersion. However, the interference introduces significant variability in the angular correlation, making it difficult to distinguish massive gravity from general relativity based on a single realization of the Universe. We conclude that beyond a fundamental constraint set by the PTA observation time, achieving a substantially tighter bound on the graviton mass is statistically challenging and observationally limited under realistic conditions.