Accelerating pulsar timing data analysis

Abstract

The analysis of pulsar timing data, especially in pulsar timing array (PTA) projects, has encountered practical difficulties: evaluating the likelihood and/or correlation-based statistics can become prohibitively computationally expensive for large datasets. In situations where a stochastic signal of interest has a power spectral density that dominates the noise in a limited bandwidth of the total frequency domain (e.g. the isotropic background of gravitational waves), a linear transformation exists that transforms the timing residuals to a basis in which virtually all the information about the stochastic signal of interest is contained in a small fraction of basis vectors. By only considering such a small subset of these "generalised residuals", the dimensionality of the data analysis problem is greatly reduced, which can cause a large speedup in the evaluation of the likelihood: the ABC-method (Acceleration By Compression). The compression fidelity, calculable with crude estimates of the signal and noise, can be used to determine how far a dataset can be compressed without significant loss of information. Both direct tests on the likelihood, and Bayesian analysis of mock data, show that the signal can be recovered as well as with an analysis of uncompressed data. In the analysis of IPTA Mock Data Challenge datasets, speedups of a factor of three orders of magnitude are demonstrated. For realistic PTA datasets the acceleration may become greater than six orders of magnitude due to the low signal to noise ratio.

9 pages, 5 figures, accepted by MNRAS

References in corpus (16)

• TEMPO2, a new pulsar timing package. I: Overview
• Tests of general relativity from timing the double pulsar
• Tempo2, a new pulsar timing package. II: The timing model and precision estimates
• Use and Abuse of the Fisher Information Matrix in the Assessment of Gravitational-Wave Parameter-Estimation Prospects
• Gravitational radiation from cosmic (super)strings: bursts, stochastic background, and observational windows
• The stochastic gravitational-wave background from massive black hole binary systems: implications for observations with Pulsar Timing Arrays
• Self-Regulated Growth of Supermassive Black Holes in Galaxies as the Origin of the Optical and X-ray Luminosity Functions of Quasars
• Gravitational Waves Probe the Coalescence Rate of Massive Black Hole Binaries
• Gravitational-Wave Stochastic Background from Kinks and Cusps on Cosmic Strings
• Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data
• Assessing the Role of Spin Noise in the Precision Timing of Millisecond Pulsars
• On measuring the gravitational-wave background using Pulsar Timing Arrays
• Constraints on cosmic string tension imposed by the limit on the stochastic gravitational wave background from the European Pulsar Timing Array
• Understanding and analysing time-correlated stochastic signals in pulsar timing
• Relic Gravitational Waves and Cosmology
• The optimal schedule for pulsar timing array observations

Cited by in corpus (13)

• The stochastic background: scaling laws and time to detection for pulsar timing arrays
• PINT: A Modern Software Package for Pulsar Timing
• Characterising gravitational wave stochastic background anisotropy with Pulsar Timing Arrays
• NANOGrav Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries in Circular Orbits
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• Searching For Anisotropic Gravitational-wave Backgrounds Using Pulsar Timing Arrays
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Dust Reddened Quasars in FIRST and UKIDSS: Beyond the Tip of the Iceberg
• An Efficient Approximation to the Likelihood for Gravitational Wave Stochastic Background Detection Using Pulsar Timing Data
• Use Model Averaging instead of Model Selection in Pulsar Timing
• A Bayesian method for pulsar template generation
• Analysis of the first IPTA Mock Data Challenge by the EPTA timing data analysis working group
• Rapid inference for individual binaries and a stochastic background with pulsar timing array data