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

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.

arXiv comment22 pages, 13 figures. Accepted to Physical Review D

References in corpus (47)

• 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 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background
• The stochastic gravitational-wave background from massive black hole binary systems: implications for observations with Pulsar Timing Arrays
• The NANOGrav 11-year Data Set: Pulsar-timing Constraints On The Stochastic Gravitational-wave Background
• The CHIME Fast Radio Burst Project: System Overview
• The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
• Science with the Square Kilometer Array: Motivation, Key Science Projects, Standards and Assumptions
• Gravitational Waves From a Dark (Twin) Phase Transition
• The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
• On the evidence for a common-spectrum process in the search for the nanohertz gravitational-wave background with the Parkes Pulsar Timing Array
• The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background
• The Astrophysics of Nanohertz Gravitational Waves
• Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search
• Expected properties of the first gravitational wave signal detected with pulsar timing arrays
• Measurement and correction of variations in interstellar dispersion in high-precision pulsar timing
• Gravitational-wave cosmology across 29 decades in frequency
• Pulsar timing analysis in the presence of correlated noise
• The Gravitational Wave Background from Massive Black Hole Binaries in Illustris: spectral features and time to detection with pulsar timing arrays
• The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Astrophysics Milestones For Pulsar Timing Array Gravitational Wave Detection
• Insights on the astrophysics of supermassive black hole binaries from pulsar timing observations
• Nested Sampling Methods
• Realistic sensitivity curves for pulsar timing arrays
• Whispers from the dark side: Confronting light new physics with NANOGrav data
• Massive black hole binary systems and the NANOGrav 12.5 year results
• The CHIME Pulsar Project: System Overview
• Bilby-MCMC: An MCMC sampler for gravitational-wave inference
• Constraining the Solution to the Last Parsec Problem with Pulsar Timing
• Towards Robust Gravitational Wave Detection with Pulsar Timing Arrays
• Constraints On The Dynamical Environments Of Supermassive Black-hole Binaries Using Pulsar-timing Arrays
• Estimating the sensitivity of pulsar timing arrays
• Constraining alternative theories of gravity using pulsar timing arrays
• Low-rank approximations for large stationary covariance matrices, as used in the Bayesian and generalized-least-squares analysis of pulsar-timing data
• Common-spectrum process versus cross-correlation for gravitational-wave searches using pulsar timing arrays
• Accurate characterization of the stochastic gravitational-wave background with pulsar timing arrays by likelihood reweighting
• Cosmic variance in the nanohertz gravitational wave background
• Weighing The Evidence For A Gravitational-Wave Background In The First International Pulsar Timing Array Data Challenge
• Disentangling Multiple Stochastic Gravitational Wave Background Sources in PTA Datasets
• Exploring the Capabilities of Gibbs Sampling in Pulsar Timing Arrays
• Implementation of an efficient Bayesian search for gravitational wave bursts with memory in pulsar timing array data

Cited by in corpus (39)

• The Detected Stochastic Gravitational Waves and Subsolar-Mass Primordial Black Holes
• Induced Gravitational Wave interpretation of PTA data: a complete study for general equation of state
• Broken blue-tilted inflationary gravitational waves: a joint analysis of NANOGrav 15-year and BICEP/Keck 2018 data
• Probing the speed of scalar-induced gravitational waves with pulsar timing arrays
• Can the gravitational wave background feel wiggles in spacetime?
• Towards distinguishing Dirac from Majorana neutrino mass with gravitational waves
• NANOGrav hints for first-order confinement-deconfinement phase transition in different QCD-matter scenarios
• Practical approaches to analyzing PTA data: Cosmic strings with six pulsars
• Null energy condition violation during inflation and pulsar timing array observations
• Unraveling the early universe's equation of state and primordial black hole production with PTA, BBN, and CMB observations
• Spectral Variance in a Stochastic Gravitational-Wave Background From a Binary Population
• Status Report on Global Pulsar-Timing-Array Efforts to Detect Gravitational Waves
• The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background
• Impact of memory-burdened black holes on primordial gravitational waves in light of Pulsar Timing Array
• Pulsar Timing Arrays require hierarchical models
• How to Detect an Astrophysical Nanohertz Gravitational-Wave Background
• Spatial and Spectral Characterization of the Gravitational-wave Background with the PTA Optimal Statistic
• Use Model Averaging instead of Model Selection in Pulsar Timing
• Sub-GeV dark matter and nano-Hertz gravitational waves from a classically conformal dark sector
• Kinetic Mixing, Proton Decay and Gravitational Waves in SO(10)
• Signals of merging supermassive primordial black holes in pulsar timing arrays
• Solving the PTA Data Analysis Problem with a Global Gibbs Scheme
• Gravitational waves from cosmic superstrings and gauge strings
• Constraining string cosmology with the gravitational-wave background using the NANOGrav 15-year data set
• Gravitational Waves from Metastable Cosmic Strings in Supersymmetric New Inflation Model
• Regularizing the Pulsar Timing Array likelihood: A path towards Fourier Space
• Supercooled phase transitions in conformal dark sectors explain NANOGrav data
• Bayesian reconstruction of primordial perturbations from induced gravitational waves
• Can the NANOGrav observations constrain the geometry of the universe?
• Inferring Mbh-Mbulge Evolution from the Gravitational Wave Background
• Constraints on ultra-slow-roll inflation with the NANOGrav 15-Year Dataset
• Implications of pulsar timing arrays for Gauss-Bonnet Inflation
• Is PSR J0514$-$4002E in a PBH-NS binary?
• Finite Populations & Finite Time: The Non-Gaussianity of a Gravitational Wave Background
• The NANOGrav 15 yr Data Set: Customized Chromatic Noise Models
• Testing No slip model with pulsar timing arrays: NANOGrav and IPTA
• Purely quadratic non-Gaussianity from tachyonic instability: Primordial black holes and scalar-induced gravitational waves
• Bayesian leave-one-out cross-validation for astrophysical model comparison using gravitational-wave background data
• Full analytic expressions of overlap reduction functions for anisotropies of the stochastic gravitational-wave background with pulsar timing arrays