Higher-order statistics of the stochastic gravitational wave background from supermassive black hole binaries
Abstract
Recent progress in gravitational wave observations has positioned pulsar timing arrays as a key tool for detecting the stochastic gravitational wave background in the nanohertz band. It is widely believed that this background is primarily attributed to the cosmic ensemble of inspiraling supermassive black hole binaries. While traditional analyses have predominantly focused on the spectral amplitude and frequency dependence of the gravitational wave background, higher-order statistics such as variance, skewness, and excess kurtosis could potentially be useful for extracting further physical information. However, these statistical moments are known to diverge when the redshift integration is extended down to z=0. In this study, we propose a strategy to resolve this issue by introducing a physically motivated lower integration limit, z_min, defined by the sensitivity for detecting individual sources. Since higher-order statistics are primarily determined by local sources, we may adopt the lowest-order approximation with respect to redshift in their computations. Under this approximation, we demonstrate that all higher-order statistics beyond the expectation value depend on the mass function only through a weighted average of the chirp mass, <\mathcal{M}^{10/3}>, irrespective of the redshift evolution model. We show that the ratio of the variance to the expectation value provides information on <\mathcal{M}^{10/3}>/<\mathcal{M}^{5/3}> independently of the total number of mergers. We also find a consistency relation between the excess kurtosis and the squared skewness, paving the way for testing the binary-origin hypothesis of the gravitational wave background. Our findings demonstrate that higher-order statistics provide a new window for interpreting the gravitational wave background, offering a methodology to break existing degeneracies and refine our understanding of the mass function.
Versions
Related
- Bridging the Population Synthesis of Supermassive Binary Black Holes and the Gravitational Wave Background
- Gravitational-wave background from supermassive black holes without merger trees: Tidally regulated mergers and direct inference from the NANOGrav 15 yr data
- Approximating the statistics of a gravitational wave background
- Correcting for Evidence Uncertainty in the Templated Background Search for the Stochastic Gravitational-Wave Background
- Polarized Anisotropic Stochastic Gravitational Wave Background Search with Ground-Based Detector Networks
BibTeX
@misc{hisamatsu2026higherorder,
title = {{Higher-order statistics of the stochastic gravitational wave background from supermassive black hole binaries}},
author = {Hinano Hisamatsu and Koutarou Kyutoku},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2605.17983},
eprint = {2605.17983},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2605.17983},
}
RIS
TY - UNPB
TI - Higher-order statistics of the stochastic gravitational wave background from supermassive black hole binaries
AU - Hinano Hisamatsu
AU - Koutarou Kyutoku
PY - 2026
DA - 2026/09/29
PB - arXiv
DO - 10.48550/arxiv.2605.17983
UR - https://arxiv.org/abs/2605.17983
AB - Recent progress in gravitational wave observations has positioned pulsar timing arrays as a key tool for detecting the stochastic gravitational wave background in the nanohertz band. It is widely believed that this background is primarily attributed to the cosmic ensemble of inspiraling supermassive black hole binaries. While traditional analyses have predominantly focused on the spectral amplitude and frequency dependence of the gravitational wave background, higher-order statistics such as variance, skewness, and excess kurtosis could potentially be useful for extracting further physical information. However, these statistical moments are known to diverge when the redshift integration is extended down to z=0. In this study, we propose a strategy to resolve this issue by introducing a physically motivated lower integration limit, z_min, defined by the sensitivity for detecting individual sources. Since higher-order statistics are primarily determined by local sources, we may adopt the lowest-order approximation with respect to redshift in their computations. Under this approximation, we demonstrate that all higher-order statistics beyond the expectation value depend on the mass function only through a weighted average of the chirp mass, <\mathcal{M}^{10/3}>, irrespective of the redshift evolution model. We show that the ratio of the variance to the expectation value provides information on <\mathcal{M}^{10/3}>/<\mathcal{M}^{5/3}> independently of the total number of mergers. We also find a consistency relation between the excess kurtosis and the squared skewness, paving the way for testing the binary-origin hypothesis of the gravitational wave background. Our findings demonstrate that higher-order statistics provide a new window for interpreting the gravitational wave background, offering a methodology to break existing degeneracies and refine our understanding of the mass function.
ER -