Work

Generalized stepping-stone sampling: efficient marginal likelihood estimation in gravitational wave analysis of pulsar timing array data

Year: 2025

Type: article

Abstract: ABSTRACT Globally, pulsar timing array (PTA) experiments have revealed tentative evidence of a stochastic gravitational wave background (GWB) signal in PTA data. Apart from acquiring more observations, the sensitivity of PTA experiments can be increased by improving the accuracy of noise modelling. In PTA data analysis, noise modelling is conducted primarily using Bayesian statistics, relying on the marginal likelihoods and Bayes factors to assess evidence. We introduce generalized stepping-stone sampling (GSS) as an efficient and accurate marginal likelihood estimation method for the PTA–Bayesian framework. This method enables low-cost estimates with high accuracy, especially when comparing expensive models such as the Hellings–Downs model or the overlap reduction function model. We demonstrate the efficiency and the accuracy of GSS for model selection and evidence calculation by reevaluating the evidence of previous analyses from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 15-yr data set and the European Pulsar Timing Array (EPTA) second data release. We find similar evidence for the GWB compared to the one reported by the NANOGrav 15-yr data set. Compared to the evidence reported for the EPTA second data release, we find a substantial increase in evidence for a GWB signal across all data sets.

Source: Monthly Notices of the Royal Astronomical Society

DOI: 10.1093/mnras/staf953

Authors: El Mehdi Zahraoui, Patricio Maturana-Russel, Willem van Straten, Renate Meyer, Sergei Gulyaev

Institutions: Auckland University of Technology, University of Auckland

Language: English

Topic: Pulsars and Gravitational Waves Research

Subfield: Astronomy and Astrophysics

Field: Physics and Astronomy

Domain: Physical Sciences

Open Access status: gold

Estimated APC paid: $3,046