Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
read the original abstract
A new model independent method is presented for the analysis of pulsar timing data and the estimation of the spectral properties of an isotropic gravitational wave background (GWB). We show that by rephrasing the likelihood we are able to eliminate the most costly aspects of computation normally associated with this type of data analysis. When applied to the International Pulsar Timing Array Mock Data Challenge data sets this results in speedups of approximately 2 to 3 orders of magnitude compared to established methods. We present three applications of the new likelihood. In the low signal to noise regime we sample directly from the power spectrum coefficients of the GWB signal realization. In the high signal to noise regime, where the data can support a large number of coefficients, we sample from the joint probability density of the power spectrum coefficients for the individual pulsars and the GWB signal realization. Critically in both these cases we need make no assumptions about the form of the power spectrum of the GWB, or the individual pulsars. Finally we present a method for characterizing the spatial correlation between pulsars on the sky, making no assumptions about the form of that correlation, and therefore providing the only truly general Bayesian method of confirming a GWB detection from pulsar timing data.
This paper has not been read by Pith yet.
Forward citations
Cited by 7 Pith papers
-
A practical theorem on gravitational-wave background statistics
For large but finite source counts, the PDF of rescaled GWB characteristic strain squared follows the universal form N^{1/3} times the reflected map-Airy distribution evaluated at N^{1/3}(y-1), fully determined by the...
-
Temporal Correlation Statistic for Intrinsic Phase Fluctuation in Double White Dwarf Gravitational-Wave Signals
Introduces a minimal quadratic statistic to isolate temporal correlations in intrinsic phase fluctuations of double white dwarf binaries and derives an analytic SNR scaling with total observation time and intrinsic co...
-
Finding Supermassive Black Hole Binary Mergers in Pulsar Timing Array Data
A complete SMBHB waveform model enables unified PTA searches for mergers and memory signals, with parameter recovery shown on simulated data for 10^8-10^10 solar mass systems.
-
The NANOGrav 15 yr Data Set: Customized Chromatic Noise Models
Customized chromatic noise models for 67 pulsars detect non-dispersive delays in 21 cases, alter achromatic noise inferences in 19, and enable solar wind density estimates over 1.5 cycles.
-
Expectations for the first supermassive black-hole binary resolved by PTAs I: Model efficacy
Simulations of PTA data show that a full gravitational-wave signal template achieves the highest Bayes factors and most robust parameter estimation for individual supermassive black hole binaries compared to an Earth-...
-
Evaluating the Fourier Approximation in Pulsar Timing Array Analysis
The Fourier approximation in PTA analysis for power-law PSDs produces marginal likelihoods roughly twice as large as more accurate calculations, but the key comparison between uncorrelated and Hellings-Downs correlate...
-
The NANOGrav 15 yr Data Set: Impacts of Customized Chromatic Noise Models on Gravitational Wave Analyses
Customized chromatic noise models applied to NANOGrav 15 yr data raise the Bayes factor for Hellings-Downs GWB correlations by a factor of ~8, lower the amplitude to 2.1e-15, and increase the spectral index to 3.5.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.