For large but finite source counts, the PDF of rescaled GWB characteristic strain squared follows the universal form N^{1/3} times the reflected map-Airy distribution evaluated at N^{1/3}(y-1), fully determined by the mean strain and a new cubic shot-noise scale.
The Heavy Tailed Non-Gaussianity of the Supermassive Black Hole Gravitational Wave Background
3 Pith papers cite this work. Polarity classification is still indexing.
3
Pith papers citing it
abstract
We study the non-Gaussian features of the gravitational wave (GW) background generated by a population of inspiraling supermassive black hole (SMBH) binaries. We show that the SMBH GW amplitude distribution (GWAD) features a universal heavy power-law tail $\propto A^{-4}$, while the low-amplitude tail depends on the SMBH merger rate and the energy-loss mechanisms of the binaries. The distribution of the induced timing residuals inherits this heavy tail. As a result, the ensemble averaged statistical moments of order three and higher diverge, limiting their usefulness as measures of non-Gaussianity, and the GW background from SMBH binaries exhibits the single loud source principle, according to which the strongest signals are more likely to be caused by a small number of loud sources. We confirm that the variance-averaged Gaussian approximation accurately describes the timing residual statistics. This approximation justifies a factored likelihood structure that combines standard Gaussian-process PTA posteriors with the non-Gaussian population prior, enabling consistent incorporation of non-Gaussian effects into SMBH model inference. We provide a fast and flexible Python implementation to compute the distribution of timing residuals from a given SMBH merger rate or GWAD.
citation-role summary
background 2
citation-polarity summary
years
2026 3verdicts
UNVERDICTED 3roles
background 2polarities
background 2representative citing papers
With a physically motivated z_min cutoff, higher-order moments of the SGWB from SMBH binaries depend on the mass function solely via <M^{10/3}>, giving a variance-to-mean ratio for <M^{10/3}>/<M^{5/3}> and a kurtosis-skewness consistency relation independent of redshift evolution.
PTA statistical tests lose sensitivity to non-Gaussian GW features after decorrelation and cannot distinguish them model-agnostically.
citing papers explorer
-
A practical theorem on gravitational-wave background statistics
For large but finite source counts, the PDF of rescaled GWB characteristic strain squared follows the universal form N^{1/3} times the reflected map-Airy distribution evaluated at N^{1/3}(y-1), fully determined by the mean strain and a new cubic shot-noise scale.
-
Higher-order statistics of the stochastic gravitational wave background from supermassive black hole binaries
With a physically motivated z_min cutoff, higher-order moments of the SGWB from SMBH binaries depend on the mass function solely via <M^{10/3}>, giving a variance-to-mean ratio for <M^{10/3}>/<M^{5/3}> and a kurtosis-skewness consistency relation independent of redshift evolution.
-
Are PTA measurements sensitive to gravitational wave non-Gaussianities?
PTA statistical tests lose sensitivity to non-Gaussian GW features after decorrelation and cannot distinguish them model-agnostically.