Spatial and Spectral Characterization of the Gravitational-wave Background with the PTA Optimal Statistic

Abstract

Pulsar timing arrays (PTAs) have made tremendous progress and are now showing strong evidence for the gravitational-wave background (GWB). Further probing the origin and characteristics of the GWB will require more generalized analysis techniques. Bayesian methods are most often used but can be computationally expensive. On the other hand, frequentist methods, like the PTA Optimal Statistic (OS), are more computationally efficient and can produce results that are complementary to Bayesian methods, allowing for stronger statistical cases to be built from a confluence of different approaches. In this work we expand the capabilities of the OS through a technique we call the Per-Frequency Optimal Statistic (PFOS). The PFOS removes the underlying power-law assumption inherent in previous implementations of the OS, and allows one to estimate the GWB spectrum in a frequency-by-frequency manner. We have also adapted a recent generalization from the OS pipeline into the PFOS, making it capable of accurately characterizing the spectrum in the intermediate and strong GW signal regimes using only a small fraction of the necessary computational resources when compared with fully-correlated Bayesian methods, while also empowering many new types of analyses not possible before. We find that even in the strong GW signal regime, where the GWB dominates over noise in all frequencies, the injected value of the signal lies within the 50th-percentile of the PFOS uncertainty distribution in 41-45% of simulations, remaining 3$\sigma$-consistent with unbiased estimation.

23 pages, 9 figures, Accepted by Physical Review D

References in corpus (23)

• The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background
• Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array
• Searching for the nano-Hertz stochastic gravitational wave background with the Chinese Pulsar Timing Array Data Release I
• The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals
• The NANOGrav 15-year Data Set: Search for Signals from New Physics
• The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
• The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars
• Expected properties of the first gravitational wave signal detected with pulsar timing arrays
• The Parkes Pulsar Timing Array Third Data Release
• Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data
• The Need For Speed: Rapid Refitting Techniques for Bayesian Spectral Characterization of the Gravitational Wave Background Using PTAs
• Constraining the Solution to the Last Parsec Problem with Pulsar Timing
• Variance of the Hellings-Downs Correlation
• Hellings and Downs correlation of an arbitrary set of pulsars
• Cosmological Interpretation for the Stochastic Signal in Pulsar Timing Arrays
• Forecasting pulsar timing array sensitivity to anisotropy in the stochastic gravitational wave background
• Common-spectrum process versus cross-correlation for gravitational-wave searches using pulsar timing arrays
• The NANOGrav 15-year Gravitational-Wave Background Methods
• Accurate characterization of the stochastic gravitational-wave background with pulsar timing arrays by likelihood reweighting
• Generalized optimal statistic for characterizing multiple correlated signals in pulsar timing arrays
• Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: I. The optimal statistic
• Efficient Gravitational Wave Searches with Pulsar Timing Arrays using Hamiltonian Monte Carlo
• Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: II. Posterior predictive distributions and pseudo Bayes factors

Cited by in corpus (6)

• Solving the PTA Data Analysis Problem with a Global Gibbs Scheme
• Searching for Anisotropy in the Gravitational Wave Background Using the Parkes Pulsar Timing Array
• Dissecting the nanoHz gravitational wave sky: frequency-correlated anisotropy induced by eccentric supermassive black hole binaries
• Rapid Construction of Joint Pulsar Timing Array Datasets: The Lite Method
• On the calculation of p-values for quadratic statistics in Pulsar Timing Arrays
• Pulsar timing array analysis in a Legendre polynomial basis