Abstract:Pulsar timing arrays (PTAs) have recently found evidence for a nanohertz-frequency stochastic gravitational-wave background (SGWB). Constraining its spectral characteristics will reveal its origins. To achieve this, we must understand how data and modeling conditions in each pulsar influence the precision and accuracy of SGWB spectral recovery, typically requiring many Bayesian analyses on real data sets and large-scale simulations that are slow and computationally taxing. To combat this, we have developed several new rapid approaches that operate on intermediate SGWB analysis products. These techniques refit SGWB spectral models to previously computed Bayesian posterior estimates of the timing power spectra. We test our new techniques on simulated PTA data sets and the NANOGrav 12.5-year data set, where in the latter our refit posterior achieves a Hellinger distance -- bounded between 0 for identical distributions and 1 for zero overlap -- from the current full production-level pipeline that is < 0.1. Our techniques are ~ $10^2$--$10^4$ times faster than the production-level likelihood and scale much more favorably (sub-linearly) as a PTA is expanded with new pulsars or observations. Our techniques also allow us to demonstrate conclusively that SGWB spectral characterization in PTA data sets is driven by the longest-timed pulsars and the best-measured power spectral densities. Indeed, the common-process spectral properties found in the NANOGrav 12.5-year data set are given by analyzing only the ~14 longest-timed pulsars out of the full 45 pulsar array, and we find that the 'shallowing' of the common-process power-law model occurs when gravitational-wave frequencies higher than ~50 nanohertz are included. The implementation of our techniques is openly available as a software suite to allow fast and flexible PTA SGWB spectral characterization and model selection.

Efficient prescription to search for linear gravitational wave memory from hyperbolic black hole encounters and its application to the NANOGrav 12.5-year dataset

Efficient prescription to search for linear gravitational wave memory from hyperbolic black hole encounters and its application to the NANOGrav 12.5-year dataset

The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory

The NANOGrav 12.5 yr Data Set: Search for Gravitational Wave Memory

Search for the gravitational wave memory effect with the Parkes Pulsar Timing Array

Searching for gravitational wave burst in PTA data with piecewise linear functions

The Need For Speed: Rapid Refitting Techniques for Bayesian Spectral Characterization of the Gravitational Wave Background Using PTAs

The Gravitational Wave Background from Massive Black Hole Binaries in Illustris: spectral features and time to detection with pulsar timing arrays

NANOGrav 15-year gravitational-wave background methods

Impact of memory-burdened black holes on primordial gravitational waves in light of Pulsar Timing Array

The NANOGrav 15 yr data set: Posterior predictive checks for gravitational-wave detection with pulsar timing arrays

Modeling nonstationary noise in pulsar timing array data analysis

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array

Characterizing the nanohertz gravitational wave background using a $t$-process power spectral density

Modeling non stationary noise in pulsar timing array data analysis

The NANOGrav 12.5 yr Data Set: Search for an Isotropic Stochastic Gravitational-wave Background

Efficient Bayesian inference and model selection for continuous gravitational waves in pulsar timing array data