State-space algorithm for detecting the nanohertz gravitational wave background

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Kimpson, Tom, Melatos, Andrew, O'Leary, Joseph, Carlin, Julian B., Evans, Robin J., Moran, William, Cheunchitra, Tong, Dong, Wenhao, Dunn, Liam, Greentree, Julian, O'Neill, Nicholas J., Suvorova, Sofia, Thong, Kok Hong, Vargas, Andrés F.
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866915130661404672
author Kimpson, Tom
Melatos, Andrew
O'Leary, Joseph
Carlin, Julian B.
Evans, Robin J.
Moran, William
Cheunchitra, Tong
Dong, Wenhao
Dunn, Liam
Greentree, Julian
O'Neill, Nicholas J.
Suvorova, Sofia
Thong, Kok Hong
Vargas, Andrés F.
author_facet Kimpson, Tom
Melatos, Andrew
O'Leary, Joseph
Carlin, Julian B.
Evans, Robin J.
Moran, William
Cheunchitra, Tong
Dong, Wenhao
Dunn, Liam
Greentree, Julian
O'Neill, Nicholas J.
Suvorova, Sofia
Thong, Kok Hong
Vargas, Andrés F.
contents The stochastic gravitational wave background (SGWB) can be observed in the nanohertz band using a pulsar timing array (PTA). Here a computationally efficient state-space framework is developed for analysing SGWB data, in which the stochastic gravitational wave strain at Earth is tracked with a non-linear Kalman filter and separated simultaneously from intrinsic, achromatic pulsar spin wandering. The filter is combined with a nested sampler to estimate the parameters of the model, and to calculate a Bayes factor for selecting between models with and without a SGWB. The procedure extends previous state-space formulations of PTA data analysis applied to individually resolvable binary black hole sources. The performance of the new algorithm is tested on synthetic data from the first International PTA Mock Data Challenge. It is shown that the algorithm distinguishes a SGWB from pure noise for $A_{\rm gw} \geq 3 \times 10^{-14}$, where $A_{\rm gw}$ denotes the standard normalization factor for a power spectral density with power-law exponent $-13/3$. Additional, systematic validation tests are also performed with synthetic data generated independently by adjusting the injected parameters to cover astrophysically plausible ranges. Full posterior distributions are recovered and tested for accuracy. The state-space procedure is memory-light and evaluates the likelihood for a standard-sized PTA dataset in $\lesssim 10^{-1}$ s without optimization on a standard central processing unit.
format Preprint
id arxiv_https___arxiv_org_abs_2501_06990
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle State-space algorithm for detecting the nanohertz gravitational wave background
Kimpson, Tom
Melatos, Andrew
O'Leary, Joseph
Carlin, Julian B.
Evans, Robin J.
Moran, William
Cheunchitra, Tong
Dong, Wenhao
Dunn, Liam
Greentree, Julian
O'Neill, Nicholas J.
Suvorova, Sofia
Thong, Kok Hong
Vargas, Andrés F.
Instrumentation and Methods for Astrophysics
High Energy Astrophysical Phenomena
General Relativity and Quantum Cosmology
The stochastic gravitational wave background (SGWB) can be observed in the nanohertz band using a pulsar timing array (PTA). Here a computationally efficient state-space framework is developed for analysing SGWB data, in which the stochastic gravitational wave strain at Earth is tracked with a non-linear Kalman filter and separated simultaneously from intrinsic, achromatic pulsar spin wandering. The filter is combined with a nested sampler to estimate the parameters of the model, and to calculate a Bayes factor for selecting between models with and without a SGWB. The procedure extends previous state-space formulations of PTA data analysis applied to individually resolvable binary black hole sources. The performance of the new algorithm is tested on synthetic data from the first International PTA Mock Data Challenge. It is shown that the algorithm distinguishes a SGWB from pure noise for $A_{\rm gw} \geq 3 \times 10^{-14}$, where $A_{\rm gw}$ denotes the standard normalization factor for a power spectral density with power-law exponent $-13/3$. Additional, systematic validation tests are also performed with synthetic data generated independently by adjusting the injected parameters to cover astrophysically plausible ranges. Full posterior distributions are recovered and tested for accuracy. The state-space procedure is memory-light and evaluates the likelihood for a standard-sized PTA dataset in $\lesssim 10^{-1}$ s without optimization on a standard central processing unit.
title State-space algorithm for detecting the nanohertz gravitational wave background
topic Instrumentation and Methods for Astrophysics
High Energy Astrophysical Phenomena
General Relativity and Quantum Cosmology
url https://arxiv.org/abs/2501.06990