Graph-based Summary Statistics for Revealing the Stochastic Gravitational Wave Background in Pulsar Timing Arrays

Abstract

In this work, we propose a graph-based method implemented on the pulsar timing residuals (PTRs) for stochastic gravitational wave background (SGWB) detection within the nano-Hertz frequency regime and examining uncertainties of its parameters. We construct a correlation graph with pulsars as its nodes, and analyze the graph-based summary statistics, including structural characteristics of complex network, for identifying SGWB in the real and synthetic datasets. The effect of the number of pulsars, the observation time span, and the strength of the SGWB on the graph-based feature vector is evaluated. Our results demonstrate that the Discriminative Summary Statistics for common signal detection consists of the average clustering coefficient and the edge weight fluctuation. The SGWB detection conducted after the observation of a common signal and then exclusion of non-Hellings \& Downs templates is performed by the second cumulant of edge weight for angular separation thresholds $\bar{\hat{I}}\P \gtrsim 40^{\circ }$. The lowest detectable value of SGWB strain amplitude utilizing our graph-based measures at the current PTAs sensitivity is $A_{\mathrm{SGWB}}\gtrsim 1.2\times 10^{-15}$. Fisher forecasts confirmed that the uncertainty levels of ${\log }_{10}{A}_{\mathrm{SGWB}}$ and spectral index reach $1.5\%$ and $19.5\%$, respectively, at $2\ddot{I}ƒ$ confidence interval. A weak evidence for an SGWB at $\sim 2.3\ddot{I}ƒ$ level is obtained by applying our graph-based method to the NANOGrav 15-year dataset.

29 pages, 15 figures, 1 table. Matched with the published version. Including the revision in a part of method