Robust Bayesian non-linear pulsar timing and noise analysis using a sequential data-tempered ensemble sampler
Year: 2026
Type: preprint
Abstract: We develop a sequential ensemble MCMC scheme for Bayesian pulsar timing and noise analysis that uses data size as the tempering parameter. The method enables gradual and reliable convergence from potentially uninformed initial points by iteratively increasing the data size, thereby increasing sensitivity and progressively narrowing the target distribution until the posterior of the full dataset is reached. The initial iterations use sub-datasets that can be orders of magnitude smaller than the full dataset, substantially reducing computational cost, which is further lowered by employing approximate lower-rank models justified by the reduced sensitivity of these sub-datasets. Bayesian inference is performed using the final MCMC chain sampled from the true target distribution after applying appropriate burn-in. We demonstrate the method using the NANOGrav 12.5-year data of PSR B1855+09, showing progressive convergence to the target posterior and that most of the computational time is spent sampling the true target distribution rather than in the preceding tempered iterations.
Source: arXiv (Cornell University)
Author: Susobhanan, Abhimanyu
Language: English
Subfield: Astronomy and Astrophysics
Field: Physics and Astronomy
Domain: Physical Sciences
Open Access status: green
OpenAlex ID: W7220467610