Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: II. Posterior predictive distributions and pseudo Bayes factors

Abstract

The detection of nanoHertz gravitational waves through pulsar timing arrays hinges on identifying a common stochastic process affecting all pulsars in a correlated way across the sky. In the presence of other deterministic and stochastic processes affecting the time-of-arrival of pulses, a detection claim must be accompanied by a detailed assessment of the various physical or phenomenological models used to describe the data. In this study, we propose posterior predictive checks as a model-checking tool that relies on the predictive performance of the models with regards to new data. We derive and study predictive checks based on different components of the models, namely the Fourier coefficients of the stochastic process, the correlation pattern, and the timing residuals. We assess the ability of our checks to identify model misspecification in simulated datasets. We find that they can accurately flag a stochastic process spectral shape that deviates from the common power-law model as well as a stochastic process that does not display the expected angular correlation pattern. Posterior predictive likelihoods derived under different assumptions about the correlation pattern can further be used to establish detection significance. In the era of nanoHertz gravitational wave detection from different pulsar-timing datasets, such tests represent an essential tool in assessing data consistency and supporting astrophysical inference.

arXiv comment18 pages, 9 figures

References in corpus (61)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• 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
• Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo
• 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
• Tempo2, a new pulsar timing package. II: The timing model and precision estimates
• Extracting distribution parameters from multiple uncertain observations with selection biases
• The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
• The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
• The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background
• Assessing the Role of Spin Noise in the Precision Timing of Millisecond Pulsars
• Optimal strategies for gravitational wave stochastic background searches in pulsar timing data
• The NANOGrav Nine-year Data Set: Observations, Arrival Time Measurements, and Analysis of 37 Millisecond Pulsars
• On measuring the gravitational-wave background using Pulsar Timing Arrays
• PINT: A Modern Software Package for Pulsar Timing
• The NANOGrav 12.5 yr Data Set: Observations and Narrowband Timing of 47 Millisecond Pulsars
• The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• Understanding and analysing time-correlated stochastic signals in pulsar timing
• Mapping gravitational-wave backgrounds using methods from CMB analysis: Application to pulsar timing arrays
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• A study of spatial correlations in pulsar timing array data
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Astrophysics Milestones For Pulsar Timing Array Gravitational Wave Detection
• The NANOGrav Nine-Year Data Set: Measurement and Interpretation of Variations in Dispersion Measures
• Accounting for Source Uncertainties in Analyses of Astronomical Survey Data
• Identifying and mitigating noise sources in precision pulsar timing data sets
• Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data
• The Most Massive Binary Black Hole Detections and the Identification of Population Outliers
• Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays
• Massive black hole binary systems and the NANOGrav 12.5 year results
• All correlations must die: Assessing the significance of a stochastic gravitational-wave background in pulsar-timing arrays
• Constraining the Solution to the Last Parsec Problem with Pulsar Timing
• Variance of the Hellings-Downs Correlation
• Towards Robust Gravitational Wave Detection with Pulsar Timing Arrays
• A Second Chromatic Timing Event of Interstellar Origin toward PSR J1713+0747
• Noise-marginalized optimal statistic: A robust hybrid frequentist-Bayesian statistic for the stochastic gravitational-wave background in pulsar timing arrays
• Hellings and Downs correlation of an arbitrary set of pulsars
• Constraints On The Dynamical Environments Of Supermassive Black-hole Binaries Using Pulsar-timing Arrays
• The NANOGrav 12.5-year data set: Search for Non-Einsteinian Polarization Modes in theGravitational-Wave Background
• Studying the solar system with the International Pulsar Timing Array
• Assigning confidence to inspiral gravitational wave candidates with Bayesian model selection
• Testing general relativity with gravitational waves: a reality check
• Enabling high confidence detections of gravitational-wave bursts
• Pulsar and cosmic variances of pulsar timing-array correlation measurements of the stochastic gravitational wave background
• 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
• Evaluating the prevalence of spurious correlations in pulsar timing array datasets
• Consistency of the Parkes Pulsar Timing Array Signal with a Nanohertz Gravitational-wave Background
• Model Dependence of Bayesian Gravitational-Wave Background Statistics for Pulsar Timing Arrays
• Hunting the stochastic gravitational wave background in pulsar timing array cross correlations through theoretical uncertainty
• Enhancing confidence in the detection of gravitational waves from compact binaries using signal coherence
• Bayesian Solar Wind Modeling with Pulsar Timing Arrays
• Ephemeris Errors and the Gravitational Wave Signal: Harmonic Mode Coupling in Pulsar Timing Array Searches
• Weighing The Evidence For A Gravitational-Wave Background In The First International Pulsar Timing Array Data Challenge
• Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: I. The optimal statistic
• Bayesian cross validation for gravitational-wave searches in pulsar-timing array data
• Gravitational Wave Searches with Pulsar Timing Arrays. I: Cancellation of Clock and Ephemeris Noises

Cited by in corpus (6)

• Spatial and Spectral Characterization of the Gravitational-wave Background with the PTA Optimal Statistic
• Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: I. The optimal statistic
• Choosing suitable noise models for nanohertz gravitational-wave astrophysics
• The NANOGrav 15 yr data set: Posterior predictive checks for gravitational-wave detection with pulsar timing arrays
• Rapid Construction of Joint Pulsar Timing Array Datasets: The Lite Method
• The NANOGrav 15 yr Data Set: Customized Chromatic Noise Models