Accurate characterization of the stochastic gravitational-wave background with pulsar timing arrays by likelihood reweighting

Abstract

An isotropic stochastic background of nanohertz gravitational waves creates excess residual power in pulsar-timing-array datasets, with characteristic inter-pulsar correlations described by the Hellings-Downs function. These correlations appear as nondiagonal terms in the noise covariance matrix, which must be inverted to obtain the pulsar-timing-array likelihood. Searches for the stochastic background, which require many likelihood evaluations, are therefore quite computationally expensive. We propose a more efficient method: we first compute approximate posteriors by ignoring cross correlations, and then reweight them to exact posteriors via importance sampling. We show that this technique results in accurate posteriors and marginal likelihood ratios, because the approximate and exact posteriors are similar, which makes reweighting especially accurate. The Bayes ratio between the marginal likelihoods of the exact and approximate models, commonly used as a detection statistic, is also estimated reliably by our method, up to ratios of at least $10^6$.

arXiv comment10 pages, 5 figures

References in corpus (39)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background
• The NANOGrav 11-year Data Set: Pulsar-timing Constraints On The Stochastic Gravitational-wave Background
• 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
• GW190521: orbital eccentricity and signatures of dynamical formation in a binary black hole merger signal
• Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search
• On measuring the gravitational-wave background using Pulsar Timing Arrays
• Understanding and analysing time-correlated stochastic signals in pulsar timing
• 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
• Identifying and mitigating noise sources in precision pulsar timing data sets
• From Bright Binaries To Bumpy Backgrounds: Mapping Realistic Gravitational Wave Skies With Pulsar-Timing Arrays
• All correlations must die: Assessing the significance of a stochastic gravitational-wave background in pulsar-timing arrays
• Higher order gravitational-wave modes with likelihood reweighting
• Constraining the Solution to the Last Parsec Problem with Pulsar Timing
• Bayesian model selection without evidences: application to the dark energy equation-of-state
• On the origin of GW190425
• Towards Robust Gravitational Wave Detection with Pulsar Timing Arrays
• Constraints On The Dynamical Environments Of Supermassive Black-hole Binaries Using Pulsar-timing Arrays
• The impact of Solar wind variability on pulsar timing
• Forecasting pulsar timing array sensitivity to anisotropy in the stochastic gravitational wave background
• Studying the solar system with the International Pulsar Timing Array
• The optimal schedule for pulsar timing array observations
• Common-spectrum process versus cross-correlation for gravitational-wave searches using pulsar timing arrays
• On the usefulness of existing Solar-wind models for pulsar timing corrections
• 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
• Optimal Frequency Ranges for Sub-Microsecond Precision Pulsar Timing
• Bayesian Solar Wind Modeling with Pulsar Timing Arrays
• Ephemeris Errors and the Gravitational Wave Signal: Harmonic Mode Coupling in Pulsar Timing Array Searches
• Optimizing Pulsar Timing Arrays Observations for Detection and Characterization of Low-Frequency Gravitational Wave Sources
• Quantifying Uncertainty in Transdimensional Markov Chain Monte Carlo Using Discrete Markov Models
• Efficient Gravitational Wave Searches with Pulsar Timing Arrays using Hamiltonian Monte Carlo
• Gravitational Wave Searches with Pulsar Timing Arrays. I: Cancellation of Clock and Ephemeris Noises
• A Parallelized Bayesian Approach To Accelerated Gravitational-Wave Background Characterization

Cited by in corpus (25)

• 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
• The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals
• The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• The Need For Speed: Rapid Refitting Techniques for Bayesian Spectral Characterization of the Gravitational Wave Background Using PTAs
• The NANOGrav 15-year Gravitational-Wave Background Methods
• Practical approaches to analyzing PTA data: Cosmic strings with six pulsars
• Identifying Host Galaxies of Supermassive Black Hole Binaries Found by PTAs
• How to Detect an Astrophysical Nanohertz Gravitational-Wave Background
• The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B
• 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
• Testing strengths, limitations and biases of current Pulsar Timing Arrays detection analyses on realistic data
• Exploring the Capabilities of Gibbs Sampling in Pulsar Timing Arrays
• Post-Newtonian-accurate pulsar timing array signals induced by inspiralling eccentric binaries: accuracy, computational cost, and single-pulsar search
• Separating deterministic and stochastic gravitational wave signals in realistic pulsar timing array datasets
• Projections of the uncertainty on the compact binary population background using popstock
• Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: II. Posterior predictive distributions and pseudo Bayes factors
• Choosing suitable noise models for nanohertz gravitational-wave astrophysics
• Rapid Construction of Joint Pulsar Timing Array Datasets: The Lite Method
• Searching for continuous gravitational waves in the Parkes Pulsar Timing Array Data Release 3
• Glitches far from transient gravitational-wave events do not bias inference
• 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 15 yr Data Set: Customized Chromatic Noise Models
• Searching for Exotrojans in Pulsar Systems