Testing strengths, limitations, and biases of current pulsar timing arrays'detection analyses on realistic data

Research output: Contribution to journalArticleResearchpeer review

Authors

  • Serena Valtolina
  • Golam Shaifullah
  • Anuradha Samajdar
  • Alberto Sesana

Research Organisations

External Research Organisations

  • University of Milan - Bicocca
  • Istituto Nazionale di Fisica Nucleare (INFN)
  • Istituto Nazionale di Astrofisica (INAF)
  • Utrecht University
  • University of Potsdam
  • Osservatorio Astronomico di Brera
  • Max Planck Institute for Gravitational Physics (Albert Einstein Institute)
View graph of relations

Details

Original languageEnglish
Article numberA201
Number of pages14
JournalAstronomy and astrophysics
Volume683
Early online date20 Mar 2024
Publication statusPublished - Mar 2024

Abstract

State-of-the-art searches for gravitational waves (GWs) in pulsar timing array (PTA) datasets model the signal as an isotropic, Gaussian, and stationary process described by a power law. In practice, none of these properties are expected to hold for an incoherent superposition of GWs generated by a cosmic ensemble of supermassive black hole binaries (SMBHBs). This stochastic signal is usually referred to as the GW background (GWB) and is expected to be the primary signal in the PTA band. We performed a systematic investigation of the performance of current search algorithms, using a simple power-law model to characterise GW signals in realistic datasets. We used, as the baseline dataset, synthetic realisations of timing residuals mimicking the European PTA (EPTA) second data release (DR2). Thus, we included in the dataset uneven time stamps, achromatic and chromatic red noise, and multi-frequency observations. We then injected timing residuals from an ideal isotropic, Gaussian, single power-law stochastic process and from a realistic population of SMBHBs, performing a methodical investigation of the recovered signal. We found that current search models are efficient at recovering the GW signal, but several biases can be identified due to the signal-template mismatch, which we identified via probability-probability (P- P) plots and quantified using Kolmogorov-Smirnov (KS) statistics. We discuss our findings in light of the signal observed in the EPTA DR2 and corroborate its consistency with a SMBHB origin.

Keywords

    Black hole physics, Gravitation, Gravitational waves, Methods: data analysis, Pulsars: general

ASJC Scopus subject areas

Cite this

Testing strengths, limitations, and biases of current pulsar timing arrays'detection analyses on realistic data. / Valtolina, Serena; Shaifullah, Golam; Samajdar, Anuradha et al.
In: Astronomy and astrophysics, Vol. 683, A201, 03.2024.

Research output: Contribution to journalArticleResearchpeer review

Valtolina S, Shaifullah G, Samajdar A, Sesana A. Testing strengths, limitations, and biases of current pulsar timing arrays'detection analyses on realistic data. Astronomy and astrophysics. 2024 Mar;683:A201. Epub 2024 Mar 20. doi: 10.48550/arXiv.2309.13117, 10.1051/0004-6361/202348084
Valtolina, Serena ; Shaifullah, Golam ; Samajdar, Anuradha et al. / Testing strengths, limitations, and biases of current pulsar timing arrays'detection analyses on realistic data. In: Astronomy and astrophysics. 2024 ; Vol. 683.
Download
@article{ac4da8b278404f4f8ff10f54b347c9a3,
title = "Testing strengths, limitations, and biases of current pulsar timing arrays'detection analyses on realistic data",
abstract = "State-of-the-art searches for gravitational waves (GWs) in pulsar timing array (PTA) datasets model the signal as an isotropic, Gaussian, and stationary process described by a power law. In practice, none of these properties are expected to hold for an incoherent superposition of GWs generated by a cosmic ensemble of supermassive black hole binaries (SMBHBs). This stochastic signal is usually referred to as the GW background (GWB) and is expected to be the primary signal in the PTA band. We performed a systematic investigation of the performance of current search algorithms, using a simple power-law model to characterise GW signals in realistic datasets. We used, as the baseline dataset, synthetic realisations of timing residuals mimicking the European PTA (EPTA) second data release (DR2). Thus, we included in the dataset uneven time stamps, achromatic and chromatic red noise, and multi-frequency observations. We then injected timing residuals from an ideal isotropic, Gaussian, single power-law stochastic process and from a realistic population of SMBHBs, performing a methodical investigation of the recovered signal. We found that current search models are efficient at recovering the GW signal, but several biases can be identified due to the signal-template mismatch, which we identified via probability-probability (P- P) plots and quantified using Kolmogorov-Smirnov (KS) statistics. We discuss our findings in light of the signal observed in the EPTA DR2 and corroborate its consistency with a SMBHB origin.",
keywords = "Black hole physics, Gravitation, Gravitational waves, Methods: data analysis, Pulsars: general",
author = "Serena Valtolina and Golam Shaifullah and Anuradha Samajdar and Alberto Sesana",
note = "Funding Information: Open Access funding provided by Max Planck Society.",
year = "2024",
month = mar,
doi = "10.48550/arXiv.2309.13117",
language = "English",
volume = "683",
journal = "Astronomy and astrophysics",
issn = "0004-6361",
publisher = "EDP Sciences",

}

Download

TY - JOUR

T1 - Testing strengths, limitations, and biases of current pulsar timing arrays'detection analyses on realistic data

AU - Valtolina, Serena

AU - Shaifullah, Golam

AU - Samajdar, Anuradha

AU - Sesana, Alberto

N1 - Funding Information: Open Access funding provided by Max Planck Society.

PY - 2024/3

Y1 - 2024/3

N2 - State-of-the-art searches for gravitational waves (GWs) in pulsar timing array (PTA) datasets model the signal as an isotropic, Gaussian, and stationary process described by a power law. In practice, none of these properties are expected to hold for an incoherent superposition of GWs generated by a cosmic ensemble of supermassive black hole binaries (SMBHBs). This stochastic signal is usually referred to as the GW background (GWB) and is expected to be the primary signal in the PTA band. We performed a systematic investigation of the performance of current search algorithms, using a simple power-law model to characterise GW signals in realistic datasets. We used, as the baseline dataset, synthetic realisations of timing residuals mimicking the European PTA (EPTA) second data release (DR2). Thus, we included in the dataset uneven time stamps, achromatic and chromatic red noise, and multi-frequency observations. We then injected timing residuals from an ideal isotropic, Gaussian, single power-law stochastic process and from a realistic population of SMBHBs, performing a methodical investigation of the recovered signal. We found that current search models are efficient at recovering the GW signal, but several biases can be identified due to the signal-template mismatch, which we identified via probability-probability (P- P) plots and quantified using Kolmogorov-Smirnov (KS) statistics. We discuss our findings in light of the signal observed in the EPTA DR2 and corroborate its consistency with a SMBHB origin.

AB - State-of-the-art searches for gravitational waves (GWs) in pulsar timing array (PTA) datasets model the signal as an isotropic, Gaussian, and stationary process described by a power law. In practice, none of these properties are expected to hold for an incoherent superposition of GWs generated by a cosmic ensemble of supermassive black hole binaries (SMBHBs). This stochastic signal is usually referred to as the GW background (GWB) and is expected to be the primary signal in the PTA band. We performed a systematic investigation of the performance of current search algorithms, using a simple power-law model to characterise GW signals in realistic datasets. We used, as the baseline dataset, synthetic realisations of timing residuals mimicking the European PTA (EPTA) second data release (DR2). Thus, we included in the dataset uneven time stamps, achromatic and chromatic red noise, and multi-frequency observations. We then injected timing residuals from an ideal isotropic, Gaussian, single power-law stochastic process and from a realistic population of SMBHBs, performing a methodical investigation of the recovered signal. We found that current search models are efficient at recovering the GW signal, but several biases can be identified due to the signal-template mismatch, which we identified via probability-probability (P- P) plots and quantified using Kolmogorov-Smirnov (KS) statistics. We discuss our findings in light of the signal observed in the EPTA DR2 and corroborate its consistency with a SMBHB origin.

KW - Black hole physics

KW - Gravitation

KW - Gravitational waves

KW - Methods: data analysis

KW - Pulsars: general

UR - http://www.scopus.com/inward/record.url?scp=85188431655&partnerID=8YFLogxK

U2 - 10.48550/arXiv.2309.13117

DO - 10.48550/arXiv.2309.13117

M3 - Article

AN - SCOPUS:85188431655

VL - 683

JO - Astronomy and astrophysics

JF - Astronomy and astrophysics

SN - 0004-6361

M1 - A201

ER -