pith. machine review for the scientific record. sign in

arxiv: 2604.19701 · v1 · submitted 2026-04-21 · 🌌 astro-ph.CO · gr-qc

Recognition: unknown

A practical theorem on gravitational-wave background statistics

Yacine Ali-Ha\"imoud (NYU)

Authors on Pith no claims yet

Pith reviewed 2026-05-10 01:36 UTC · model grok-4.3

classification 🌌 astro-ph.CO gr-qc
keywords gravitational wave backgroundpulsar timing arrayssupermassive black hole binariesprobability distributionshot noisefinite source countmap-Airy distribution
0
0 comments X

The pith

The probability distribution of the gravitational wave background strain squared takes a universal self-similar form for large but finite source numbers.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper derives a simple analytic expression for the PDF of the characteristic strain squared of the nanohertz gravitational wave background produced by supermassive black hole binaries. In the regime of large but finite effective source count N, the rescaled variable y equal to the strain squared divided by its mean has a PDF given by N to the one third times a reflected map-Airy distribution evaluated at N to the one third times (y minus one). The effective N is fixed entirely by the mean strain squared and a new cubic shot-noise strain scale that depends only on local population properties. This matters for pulsar timing array observations because those arrays are most sensitive in the lowest frequency bands where the finite number of sources produces measurable deviations from the infinite-source limit.

Core claim

For N much greater than one, the PDF of y defined as h_c squared over its mean takes the universal self-similar form approximately N to the one third times P of N to the one third times (y minus one), where P is the reflected map-Airy distribution. The effective number of in-band sources N is fully specified by the mean h_c squared and the cubic shot-noise strain scale h_0 cubed. This result holds for any population of supermassive black hole binaries, whether circular or eccentric and independent of the dominant orbital hardening mechanism.

What carries the argument

The rescaled variable y equals h_c squared divided by its mean, combined with the effective source number N fixed by the cubic strain scale, which collapses the full PDF onto the reflected map-Airy distribution.

If this is right

  • Pulsar timing array analyses can insert this PDF directly into likelihoods for the lowest frequency bins instead of assuming infinitely many sources.
  • The cubic strain scale becomes a new summary statistic for the background that depends only on local properties and can be measured independently of the full distribution.
  • The large-N approximation can be tested for accuracy against any specific supermassive black hole binary population model.
  • The same self-similar form applies regardless of whether binaries are circular or eccentric and independent of the orbital evolution mechanism.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Incorporating this non-Gaussian PDF into PTA pipelines could tighten constraints on binary populations by properly accounting for shot-noise fluctuations.
  • Analogous universal distributions may appear in other stochastic backgrounds with finite Poisson sources whenever the cubic moment dominates the scaling.
  • The independence assumption could be relaxed in future work to include source correlations and test how they modify the reflected map-Airy form.

Load-bearing premise

The background is produced by a Poisson process of independent sources whose statistics are fully captured by local population properties through the mean and cubic strain scales.

What would settle it

Pulsar timing array measurements showing that the spread of h_c squared values around the mean fails to follow the predicted N to the one third scaling when frequency band or population parameters are varied would falsify the universal form.

Figures

Figures reproduced from arXiv: 2604.19701 by Yacine Ali-Ha\"imoud (NYU).

Figure 1
Figure 1. Figure 1: FIG. 1. Universal shape of the large-source-count probability [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Probability distribution of the band-averaged GWB characteristic strain squared in the large-source-count limit, given [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Numerical PDF of the band-averaged characteristic strain squared for the simple SMBHB model described in Section [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of the exact PDF [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
read the original abstract

Inspiralling supermassive black-hole binaries (SMBHBs) are expected to be the main source of the nanohertz gravitational-wave background (GWB) targeted by pulsar timing arrays (PTAs). We provide a simple and general analytic expression for the probability distribution function (PDF) of the GWB characteristic strain squared $h_c^2$ in the limit of a large but finite effective number of sources, $N$, relevant for the lowest-frequency bands where PTAs are most sensitive. Explicitly, we show that for $N \gg 1$, the PDF of the rescaled variable $y \equiv h_c^2/\overline{h_c^2}$ takes the universal self-similar form $P(y) \simeq N^{1/3} \mathcal{P}(N^{1/3} (y -1))$, where $\mathcal{P}$ is the reflected map-Airy distribution. The effective number of in-band sources $N$ is fully specified by the mean $\overline{h_c^2}$ and the cubic shot-noise strain scale $\overline{h_0^3}$, a new summary statistic of the GWB that depends only on the local properties of the SMBHB population. This result is universal: it applies to any population of SMBHBs, regardless of whether they are circular or eccentric, and of the mechanism dominating orbital hardening. We explicitly quantify the accuracy of the large-source-count PDF for a simple but physically realistic SMBHB model, and outline its practical application to PTA data analysis.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript derives an analytic expression for the PDF of the GWB characteristic strain squared h_c² from a compound Poisson process of independent SMBHB sources. For large but finite effective source count N, the rescaled variable y ≡ h_c² / mean(h_c²) has the universal self-similar form P(y) ≃ N^{1/3} P(N^{1/3}(y-1)), where P is the reflected map-Airy distribution. The effective N is fully determined by the mean strain and the new cubic shot-noise strain scale, which depends only on local population moments. The result is claimed to hold for any SMBHB population (circular or eccentric, any hardening mechanism), with accuracy quantified via direct numerical comparison to a realistic population model.

Significance. If the central derivation holds, the result supplies a practical, parameter-light analytic tool for non-Gaussian GWB statistics in PTA analyses, particularly at the lowest frequencies where the Gaussian limit is least accurate. Strengths include the explicit derivation from the compound Poisson process, the independence of the cubic scale from global PDF details, the universality statement, and the numerical validation showing convergence at relevant N. This advances beyond the central-limit Gaussian approximation while remaining computationally simple for data analysis.

minor comments (2)
  1. The definition and explicit formula for the cubic shot-noise strain scale should appear in the main text before its use in the effective-N expression.
  2. The numerical validation section would benefit from a short table or inset quantifying the fractional error of the map-Airy approximation as a function of N for the eccentric and circular cases.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive evaluation of the manuscript, their recognition of its potential utility for PTA analyses, and their recommendation to accept. No major comments were raised in the report.

Circularity Check

0 steps flagged

Derivation is self-contained from Poisson process assumptions

full rationale

The paper derives the universal self-similar PDF form P(y) ≃ N^{1/3} P(N^{1/3}(y-1)) directly from the compound Poisson statistics of independent sources in the large-N limit, with the effective N and cubic scale h0^3 defined explicitly from local population moments (mean strain and third moment) without reference to the global PDF shape. The reflected map-Airy distribution is introduced via its own integral or characteristic-function definition, and the result is validated by direct numerical sampling for concrete SMBHB populations rather than by fitting. No equation reduces to a tautology, no parameter is fitted then relabeled as a prediction, and no load-bearing step relies on a self-citation whose content is itself unverified; the chain is therefore independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the large-N limit for a Poisson sum of independent sources and on the locality of the population statistics that define the cubic scale; these are domain assumptions rather than new postulates.

axioms (1)
  • domain assumption The effective number of in-band sources N is large but finite and the sources contribute as a Poisson process of independent emitters
    Invoked to obtain the self-similar Airy form; stated in the abstract as the relevant limit for lowest-frequency PTA bands.

pith-pipeline@v0.9.0 · 5575 in / 1437 out tokens · 126305 ms · 2026-05-10T01:36:06.345604+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Are PTA measurements sensitive to gravitational wave non-Gaussianities?

    astro-ph.CO 2026-05 unverdicted novelty 6.0

    PTA statistical tests lose sensitivity to non-Gaussian GW features after decorrelation and cannot distinguish them model-agnostically.

Reference graph

Works this paper leans on

39 extracted references · 28 canonical work pages · cited by 1 Pith paper · 4 internal anchors

  1. [1]

    typical”, distant sources. On the other hand, the small-h 2 c,k limit corresponds to a collective downward fluctuation of a large number of “typical

    = 1, we may replaceJ(v/N k) with its limit J(0) = 1 4√π Z ∞ 0 dt eit −1−it t5/2 =− 1 +i 3 √ 2 ,(32) so that L(v≪N k) =− 1 +i 3 √ 2 v3/2 √Nk .(33) Inserting this expression into Eq. (28), and using L(−v) =L(v) ∗, we find that the PDF ofy=h 2 c,k/h2 c,k is approximately P(y)≃ Z ∞ 0 dv π cos (y−1)v+ v3/2 3√2Nk ×exp − v3/2 3√2Nk .(34) Finally, changing variab...

    2018

  2. [2]

    (37) is exactly equivalent to the reflected map-Airy function, Eq

    This appendix is thus dedicated to explicitly proving that Eq. (37) is exactly equivalent to the reflected map-Airy function, Eq. (38). We recall that the Airy function can be defined by its Fourier transform [37] Ai(x) = Z ∞ −∞ dz 2π eixz+iz 3/3.(73) This implies that the reflected map-Airy distribution given in Eq. (38), can be rewritten as PEq(38)(x) =...

  3. [3]

    The Nanograv Collaboration, Astrophys. J. Lett.952, L37 (2023), arXiv:2306.16220

    work page arXiv 2023

  4. [4]

    Antoniadis et al

    EPTA Collaboration, Astron. Astrophys.685, A94 (2024), arXiv:2306.16227

    work page arXiv 2024

  5. [5]

    Searching for the nano-Hertz stochastic gravitational wave background with the Chinese Pulsar Timing Array Data Release I

    H. Xuet al., Research in Astronomy and Astrophysics 23, 075024 (2023), arXiv:2306.16216 [astro-ph.HE]

    work page internal anchor Pith review arXiv 2023

  6. [6]

    D. J. Reardonet al., Astrophys. J. Lett.951, L6 (2023), arXiv:2306.16215

    work page internal anchor Pith review arXiv 2023

  7. [7]

    R. W. Hellings and G. S. Downs, Astrophys. J. Lett.265, L39 (1983)

    1983

  8. [8]

    E. S. Phinney, arXiv e-prints (2001), arXiv:astro- ph/0108028

    work page arXiv 2001

  9. [9]

    Sesana, A

    A. Sesana, A. Vecchio, and C. N. Colacino, Mon. Not. R. Astron. Soc.390, 192 (2008), arXiv:0804.4476 [astro-ph]

    work page arXiv 2008

  10. [10]

    C. M. F. Mingarelli, T. Sidery, I. Mandel, and A. Vec- chio, Phys. Rev. D88, 062005 (2013), arXiv:1306.5394

    work page arXiv 2013

  11. [11]

    The NANOGrav Collaboration, The Astrophysical Jour- nal978, 31 (2024)

    2024

  12. [12]

    The NANOGrav Collaboration, The Astrophysical Jour- nal Letters951, L11 (2023)

    2023

  13. [13]

    V. Ravi, J. S. B. Wyithe, G. Hobbs, R. M. Shannon, R. N. Manchester, D. R. B. Yardley, and M. J. Keith, Astrophys. J.761, 84 (2012), arXiv:1210.3854

    work page arXiv 2012

  14. [14]

    W. G. Lamb, J. M. Wachter, A. Mitridate, S. C. Sarde- sai, B. B´ ecsy, E. L. Hagen, S. R. Taylor, and L. Z. Kel- ley, arXiv e-prints (2025), 10.48550/arXiv.2511.09659, arXiv:2511.09659

    work page doi:10.48550/arxiv.2511.09659 2025

  15. [15]

    X. Xue, Z. Pan, and L. Dai, Phys. Rev. D111, 043022 (2025), arXiv:2409.19516

    work page arXiv 2025

  16. [16]

    Exploring the spectrum of stochastic gravitational-wave anisotropies with pulsar timing arrays,

    G. Sato-Polito and M. Kamionkowski, Phys. Rev. D109, 123544 (2024), arXiv:2305.05690

    work page arXiv 2024

  17. [17]

    Raidal, J

    J. Raidal, J. Urrutia, V. Vaskonen, and H. Veerm¨ ae, Astron. Astrophys.706, A159 (2026), arXiv:2411.19692

    work page arXiv 2026

  18. [18]

    M.-X. Lin, A. Lidz, and C.-P. Ma, arXiv e-prints (2026), 10.48550/arXiv.2602.16808, arXiv:2602.16808

    work page doi:10.48550/arxiv.2602.16808 2026

  19. [19]

    The Heavy Tailed Non-Gaussianity of the Supermassive Black Hole Gravitational Wave Background

    J. Raidal, J. Urrutia, V. Vaskonen, and H. Veerm¨ ae, (2026), 10.48550/arXiv.2604.08506, arXiv:2604.08506

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.08506 2026

  20. [20]

    N. Laal, S. R. Taylor, L. Z. Kelley, J. Simon, K. G¨ ultekin, D. Wright, B. B´ ecsy, J. A. Casey-Clyde, S. Chen, A. Cin- goranelli, D. J. D’Orazio, E. C. Gardiner, W. G. Lamb, C. Matt, M. S. Siwek, and J. M. Wachter, Astrophys. J. 982, 55 (2025), arXiv:2411.10519. 15

    work page arXiv 2025

  21. [21]

    Ellis, M

    J. Ellis, M. Fairbairn, G. H¨ utsi, M. Raidal, J. Urrutia, V. Vaskonen, and H. Veerm¨ ae, Astron. Astrophys.676, A38 (2023), arXiv:2301.13854

    work page arXiv 2023

  22. [22]

    Sato-Polito and M

    G. Sato-Polito and M. Zaldarriaga, Phys. Rev. D111, 023043 (2025), arXiv:2406.17010

    work page arXiv 2025

  23. [23]

    Banderier, P

    C. Banderier, P. Flajolet, G. Schaeffer, and M. Soria, Random Structures & Algorithms19, 194 (2001)

    2001

  24. [24]

    Ali-Ha¨ ımoud, T

    Y. Ali-Ha¨ ımoud, T. L. Smith, and C. M. F. Mingarelli, Phys. Rev. D102, 122005 (2020), arXiv:2006.14570

    work page arXiv 2020

  25. [25]

    P. C. Peters and J. Mathews, Phys. Rev.131, 435 (1963)

    1963

  26. [26]

    Planck 2018 results. VI. Cosmological parameters

    Planck Collaboration, Astron. Astrophys.641, A6 (2020), arXiv:1807.06209

    work page internal anchor Pith review arXiv 2020

  27. [27]

    Map-airy distribution,

    E. W. Weisstein, “Map-airy distribution,”https:// mathworld.wolfram.com/Map-AiryDistribution.html (2026), from MathWorld—A Wolfram Web Resource

    2026

  28. [28]

    Ellis, M

    J. Ellis, M. Fairbairn, G. H¨ utsi, J. Raidal, J. Urru- tia, V. Vaskonen, and H. Veerm¨ ae, Phys. Rev. D109, L021302 (2024), arXiv:2306.17021

    work page arXiv 2024

  29. [29]

    L. Z. Kelley, L. Blecha, L. Hernquist, A. Sesana, and S. R. Taylor, Mon. Not. R. Astron. Soc.471, 4508 (2017), arXiv:1702.02180

    work page arXiv 2017

  30. [30]

    Raidal, J

    J. Raidal, J. Urrutia, V. Vaskonen, and H. Veerm¨ ae, Astron. Astrophys.691, A212 (2024), arXiv:2406.05125

    work page arXiv 2024

  31. [31]

    Middleton, W

    H. Middleton, W. Del Pozzo, W. M. Farr, A. Sesana, and A. Vecchio, Mon. Not. R. Astron. Soc.455, L72 (2016), arXiv:1507.00992

    work page arXiv 2016

  32. [32]

    E. R. Liepold and C.-P. Ma, Astrophys. J. Lett.971, L29 (2024), arXiv:2407.14595

    work page arXiv 2024

  33. [33]

    Sato-Polito, M

    G. Sato-Polito, M. Zaldarriaga, and E. Quataert, Phys. Rev. D110, 063020 (2024)

    2024

  34. [34]

    Sato-Polito, M

    G. Sato-Polito, M. Zaldarriaga, and E. Quataert, Phys. Rev. D112, 123018 (2025), arXiv:2501.09786

    work page arXiv 2025

  35. [35]

    NANOGrav Collaboration, ApJL951, L8 (2023)

    2023

  36. [36]

    Lentati, P

    L. Lentati, P. Alexander, M. P. Hobson, S. Taylor, J. Gair, S. T. Balan, and R. van Haasteren, Phys. Rev. D87, 104021 (2013), arXiv:1210.3578

    work page arXiv 2013

  37. [37]

    Lentati, P

    L. Lentati, P. Alexander, M. P. Hobson, F. Feroz, R. van Haasteren, K. J. Lee, and R. M. Shannon, Mon. Not. R. Astron. Soc.437, 3004 (2014), arXiv:1310.2120

    work page arXiv 2014

  38. [38]

    Addario-Berry, ALEA, Lat

    L. Addario-Berry, ALEA, Lat. Am. J. Probab. Math. Stat.16, 1 (2019), arXiv:1503.08159

    work page arXiv 2019

  39. [39]

    Airy function,

    E. W. Weisstein, “Airy function,”https: //mathworld.wolfram.com/AiryFunctions.html (2026), from MathWorld—A Wolfram Web Resource

    2026