Characterizing Continuous Gravitational Waves from Supermassive Black Hole Binaries in Realistic Pulsar Timing Array Data (2502.16016v1)

Abstract: Pulsar timing arrays recently found evidence for a gravitational wave background (GWB), likely the stochastic overlap of GWs from many supermassive black hole binaries. Anticipating a continuous gravitational wave (CW) detection from a single binary soon to follow, we examine how well current Bayesian methods can detect CWs and characterize their binary properties by modeling the response of the NANOGrav 15-year pulsar timing array to simulated binary populations. We run Markov Chain Monte Carlo searches for CWs in these datasets and compare them to quicker detection statistics including the optimal signal-to-noise ratio, matched filter detection statistic, and reduced log-likelihood ratio between the signal and noise models calculated at the injected parameters. The latter is the best proxy for Bayesian detection fractions, corresponding to a 50% detection fraction (by Bayes factors >10 favoring a CW detection over noise-only model) at a signal-to-noise ratio of 4.6. Source confusion between the GWB and a CW, or between multiple CWs, can cause false detections and unexpected dismissals. 53% of realistic binary populations consistent with the recently observed GWB have successful CW detections. 82% of these CWs are in the 4th or 5th frequency bin of the 16.03 yr dataset (6.9 nHz and 10.8 nHz), with 95 percentile regions spanning 4nHz-12nHz frequencies, $7-20\times10^9 M_\odot$ chirp masses, 60Mpc-8Gpc luminosity distances, and 18-13,000 sq. deg 68% confidence localization areas. These successful detections often poorly recover the chirp mass, with only 29% identifying the chirp mass accurately to within 1 dex with a 68% posterior width also narrower than 1 dex.