A Delayed Rejection Reversible Jump Markov Chain Monte Carlo Method for Multi-Resolution Maps of the Stochastic Gravitational Wave Background with Pulsar Timing Array Data

Abstract

Pulsar timing array anisotropy analyses often use a naive counting argument to set resolutions of inferred maps of the stochastic gravitational wave background (GWB). We present a data-driven method in the form of a delayed rejection reversible jump sampler tailored to multi-resolution pixel decompositions of the angular power density of the GWB. We also consider a rapid, frequentist alternative based on information criteria statistics. We verify our methods with a series of injection-and-recove…