Total: 1
Using public data from five pulsar timing arrays (PTAs), we search for a stationary, isotropic, and unpolarized nHz stochastic gravitational-wave background (SGWB). We use pulse time-of-arrival data from 121 pulsars, which is more informative than previous searches, carried out separately by the individual PTA collaborations using only their own data. For pulsars observed by more than one PTA, we employ a new "direct combination" method to merge their astrophysical models and data. This avoids the challenge of reconciling the different PTA timing models to obtain a single "best" model. A central result of our analysis is posterior probability distributions for the amplitude $A_{gw}$ and exponent $γ_{gw}$ of the putative SGWB energy-density spectrum, modeled as a power law in frequency. While these results are consistent with a nonzero SGWB amplitude $A_{gw}$, the statistical significance--assessed via a Bayesian odds ratio and through noise-marginalized false-alarm probabilities ($p$-values) for three different detection statistics--remains below the conventional $5σ$ threshold for a confident detection. We also reconstruct the inter-pulsar timing-residual correlation as a function of the angle $θ$ between the pulsar lines of sight. This is consistent with the curve predicted by Hellings and Downs (HD).