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$.

10 pages, 5 figures

References in corpus (20)

• Array Programming with NumPy
• The NANOGrav 12.5-year Data Set: Search For An Isotropic 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
• 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
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Identifying and mitigating noise sources in precision pulsar timing data sets
• Constraining the Solution to the Last Parsec Problem with Pulsar Timing
• The impact of Solar wind variability on pulsar timing
• 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
• 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
• Bayesian Solar Wind Modeling with Pulsar Timing Arrays
• Ephemeris Errors and the Gravitational Wave Signal: Harmonic Mode Coupling in Pulsar Timing Array Searches
• Efficient Gravitational Wave Searches with Pulsar Timing Arrays using Hamiltonian Monte Carlo
• 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