Approximating the statistics of a gravitational wave background

Mikel Falxa

REVIEW

Approximating the statistics of a gravitational wave background

Mikel Falxa

Instrumentation and Methods for Astrophysics · Submitted 7 Sep 2026 · 48 references

Pith's verdict: Useful tool for non-Gaussian GWB statistics, but validation is too thin to support its central accuracy claim.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 9, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 8 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv ID2609.07686 (version 1)

Pith IDJSOVQTKK, Pith's permanent record of this paper and its review

CitationMikel Falxa. "Approximating the statistics of a gravitational wave background." arXiv:2609.07686 (2026). https://pith.science/paper/JSOVQTKK

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

The paper delivers what it says: a publicly available code (zamari) that computes the saddlepoint approximation for the distribution of characteristic strain from an arbitrary SMBHB population model. The math is correct, the connection to normal–variance–mean mixtures is clearly explained, and the code is a practical resource for anyone working on hierarchical PTA inference. The demonstration that a moment-matched log-normal distribution can bias parameter recovery is a useful caution.

The honest issue is the strength of the validation. The Hellinger distance between the saddlepoint approximation and the true compound-Poisson distribution is shown for one frequency (f = 1/yr) and can reach 0.3. That is not small, and it tells us nothing about performance across the full frequency grid used in inference (30 bins). The paper does not test whether posteriors derived from the approximation are consistent with those from the true distribution. Without that check, the central claim that the method is “accurate” for any population is unsupported. The Gaussian ensemble approximation (Appendix B) adds a further uncertainty at the 8% level for high-mass populations, which the paper flags but does not propagate into the inference.

For whom is this paper? It is a methods note aimed at practitioners who want to include non-Gaussian GWB statistics in PTA analysis. The code will be used. But as a scientific argument for the accuracy of the saddlepoint approximation, the evidence is currently insufficient. The paper would need a multi-frequency validation, a demonstration that inference biases are tolerable, and ideally a test on more realistic simulated data (not just white noise and 60 pulsars) before I would rely on its results.

Recommendation: send to peer review. The referee can request the missing validation, and the community will benefit from having this tool on arXiv in the meantime. But the paper in its current form should not be used to draw astrophysical conclusions.

Referee report

Major objections (3)

[Appendix C] The Hellinger-distance validation is performed only at a single frequency ($f=1/\mathrm{yr}$). The Bayesian inference in Sect. VI B uses 30 frequency bins (from $1/20~\mathrm{yr}$ to $30/20~\mathrm{yr}$), and the accuracy of the saddlepoint approximation at higher frequencies (where the expected number of sources drops) is not quantified. Given that Hellinger distances up to $0.3$ are already seen at $f=1/\mathrm{yr}$, it is unclear whether the approximation is sufficiently accurate across the full band, and whether errors accumulate or cancel in the joint likelihood. The central claim of accuracy for 'any population model' is not supported without multi-frequency validation.
[Section VI B] The inference results (Figs. 4 and 5) compare the saddlepoint-based posterior only to a log-normal posterior, not to the true posterior obtained by using the exact compound Poisson distribution (or a high-resolution Monte Carlo estimate). Without this baseline, the reader cannot assess whether the saddlepoint method yields unbiased or consistent parameter recovery. The paper should either perform a coverage test (e.g., check that the true injected values fall within credible intervals at the expected rate) or directly sample from the compound Poisson distribution in the likelihood to establish that the approximation does not introduce systematic bias.
[Appendix A, Section II] The paper adopts the long-detector-arm limit for the overlap reduction function (Eq. A2) and notes that this 'can have consequences on the inference' citing Refs. [26,27]. The impact of this approximation on the saddlepoint prediction and on the subsequent Bayesian inference is not investigated. A robustness test (e.g., using a full frequency-dependent response for a subset of simulations) would be needed to ensure that the method is reliable in the realistic scenario where the approximation may break down.

Minor objections (8)

[Abstract] The phrase 'heavy-tailed behavior distributions that is expected' should be reworded; e.g., 'heavy-tailed distributions that are expected' or 'heavy-tailed behavior that is expected'.
[Introduction] 'developped' should be 'developed'.
[Eq. (5)] The explicit dependence of $w^\pm_j(f,f_s)$ on the source frequency $f_s$ is introduced but later the approximation $w^\pm_j(f,f_s)\approx\delta(f\pm f_s)$ is used without comment on the validity of this replacement for the discrete sum. Clarify that this is a standard simplifying assumption.
[Fig. 2] The log-scale panel shows the probability density on the y-axis but the units are labelled 'Log scale' without indicating that the density itself is not log-transformed. Consider using a log-log plot or explicitly stating that the axes are logarithmic but the pdf is linear.
[Appendix C] The Hellinger distance is computed using a kernel density estimate from 5000 Monte Carlo samples. The uncertainty in the KDE itself (due to finite sample size and bandwidth choice) is not discussed. Adding error bars or a note on stability would strengthen the validation.
[Eq. (B6) and Fig. 6] The definition of $\Delta p$ involves an integral over $s$ but the numerical integration limits are not specified. Also, the caption of Fig. 6 refers to 'massive and light' populations while the text uses 'high mass' and 'low mass'. Unify terminology.
[Table II] The prior range for $\log_{10}\dot{n}_0$ is $[-5,0]$, but the injected values in Table I are $-2.9$ and $-1.74$, which are both within the range. Consider whether a wider prior would be more realistic for a general astrophysical model.
[Section III] The list of example distributions (log-normal, inverse Gaussian, gamma) is helpful, but the paper does not explain why the log-normal is chosen as the comparison distribution in Sect. VI B rather than, say, the inverse Gaussian which is analytically more tractable for NVM. A brief justification would be useful.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Saddlepoint approximation accurately estimates the non-Gaussian statistics of the gravitational wave background from supermassive black hole binaries.

The reading

The paper presents the saddlepoint approximation as a fast, flexible method to compute the full probability distribution of the gravitational wave background (GWB) strain from a population of supermassive black hole binaries. Unlike the common Gaussian ensemble approximation, this method captures the heavy-tailed non-Gaussianity caused by the discreteness of the source population. The resulting distribution can be used as a hierarchical prior in pulsar timing array (PTA) likelihoods, effectively transforming a Gaussian likelihood into a non-Gaussian one. The author argues that including these statistics is essential for unbiased inference of the binary mass function and other astrophysical parameters, especially when a few bright sources dominate the signal.

What carries the argument

The cumulant generating function (CGF) of a compound Poisson process that sums individual gravitational wave contributions over a grid of chirp mass and redshift bins. The saddlepoint approximation inverts the CGF via Laplace's method around the stationary point defined by K'(s) = $h_c^{2}$, giving a closed-form estimate of the probability density of the total characteristic strain. This machinery replaces costly Monte Carlo simulations and requires only the source distribution model and its moments.

Extended reading notes

Core claim

The central claim is that the saddlepoint approximation, derived from the cumulant generating function of a compound Poisson process, provides an accurate and computationally efficient estimate of the distribution of characteristic strain h_c for any binary population model. The paper shows through simulations that this approximation closely matches the true distribution obtained from Monte Carlo sampling, with Hellinger distances comparable to those of machine-learning density estimators (normalizing flows). By using this approximation as a prior on the free spectrum coefficients in a Gaussian PTA likelihood, one can perform Bayesian inference on population hyperparameters without constructing an explicit non-Gaussian likelihood.

Load-bearing premise

The long detector arm limit (Appendix A, Eq. A2) neglects the frequency-dependent response from pulsar terms in the cross-correlation; if this approximation fails at PTA observation frequencies, the inferred distribution of characteristic strain could be systematically biased.

Editorial extensions

If this is right

PTA analyses can incorporate the full non-Gaussian statistics of the GWB, enabling more accurate inference of the SMBHB mass function and merger rate from the observed spectrum.
The method can be used to produce realistic visualizations of the GWB cosmic variance, going beyond the traditional mean-strain curves.
The flexible code implementation allows easy coupling with semi-analytical galaxy formation models to predict the expected GWB statistics.
For populations with high characteristic masses (few bright binaries), the saddlepoint approximation significantly outperforms a simple moment-matched log-normal distribution, reducing bias in recovered parameters.
The approach can be extended to include environmental effects (e.g., eccentricity, stellar hardening) by modifying the frequency-evolution kernel in the CGF.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The same saddlepoint framework could be applied to other stochastic gravitational-wave backgrounds, such as those from cosmic strings or first-order phase transitions, where source discreteness matters.
If the long-arm-limit approximation (Eq. A2) breaks down at the lowest PTA frequencies, the inferred distribution of h_c could be systematically shifted; a frequency-dependent response function would need to be folded into the CGF.
The saddlepoint approximation's performance relative to machine learning methods degrades for strongly eccentric binaries, but it offers the advantage of interpretability and trivial portability to new population models.
A direct observational test would involve comparing posterior distributions of population parameters inferred with and without the saddlepoint prior on real PTA data, checking for significant shifts in the mass function parameters.

Original abstract

The astrophysical origin of the gravitational wave background (GWB) reported by pulsar timing array (PTA) collaborations has yet to be confirmed. A GWB signal made of the sum of individual gravitational wave (GW) from the population of supermassive black hole binaries (SMBHB) would show the imprint of a discrete Poissonian statistics in its spectral properties. In this work, we propose a tool based on the saddlepoint approximation method to estimate the distribution of characteristic strain for any given population model. This tool can be used for Bayesian inference or for a quick visualization of the statistics of the GWB from the output of more realistic semi-analytical models. We introduce the general family of variance mixture Gaussian distributions that models heavy-tailed behavior distributions that is expected for a non-Gaussian GWB signal. We show that setting the correct hierarchical priors to a Gaussian free spectrum PTA likelihood is effectively equivalent to a non-Gaussian likelihood. Using ideal simulations, we compare the performance of the saddlepoint approximation with a log-Normal distribution to infer the parameters of the astrophysical model from the statistics of the GWB and show that correctly modeling higher order moments is essential. Future PTA analyses should include the statistics of the GWB in their pipelines.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,308 input tokens.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.29 · human participants 0.01 · controlled comparison 0.29 · randomized assignment 0.01 · states quantitative result 0.09 · compares with prior results 0.87 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

The astrophysical origin of the gravitational wave background (GWB) reported by pulsar timing array (PTA) collaborations has yet to be confirmed. A GWB signal made of the sum of individual gravitational wave (GW) from the population of supermassive black hole binaries (SMBHB) would show the imprint of a discrete Poissonian statistics in its spectral properties. In this work, we propose a tool based on the saddlepoint approximation method to estimate the distribution of characteristic strain for any given population model. This tool can be used for Bayesian inference or for a quick visualization of the statistics of the GWB from the output of more realistic semi-analytical models. We introduce the general family of variance mixture Gaussian distributions that models heavy-tailed behavior distributions that is expected for a non-Gaussian GWB signal. We show that setting the correct hierarchical priors to a Gaussian free spectrum PTA likelihood is effectively equivalent to a non-Gaussian likelihood. Using ideal simulations, we compare the performance of the saddlepoint approximation with a log-Normal distribution to infer the parameters of the astrophysical model from the statistics of the GWB and show that correctly modeling higher order moments is essential. Future PTA analyses should include the statistics of the GWB in their pipelines.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	The astrophysical origin of the gravitational wave background (GWB) reported by pulsar timing array (PTA) collaborations has yet to be confirmed.	background 0.80	0.05	0.03	asserted 0.92
s1	A GWB signal made of the sum of individual gravitational wave (GW) from the population of supermassive black hole binaries (SMBHB) would show the imprint of a discrete Poissonian statistics in its spectral properties.	background 0.99	0.13	0.09	asserted 0.66
s2	In this work, we propose a tool based on the saddlepoint approximation method to estimate the distribution of characteristic strain for any given population model.	contribution 0.67	0.92	0.12	asserted 0.92
s3	This tool can be used for Bayesian inference or for a quick visualization of the statistics of the GWB from the output of more realistic semi-analytical models.	contribution 0.71	0.45	0.09	hedged possibility 0.52
s4	We introduce the general family of variance mixture Gaussian distributions that models heavy-tailed behavior distributions that is expected for a non-Gaussian GWB signal.	contribution 0.78	0.85	0.08	asserted 0.94
s5	We show that setting the correct hierarchical priors to a Gaussian free spectrum PTA likelihood is effectively equivalent to a non-Gaussian likelihood.	result 0.83	0.83	0.21	asserted 1.00
s6	Using ideal simulations, we compare the performance of the saddlepoint approximation with a log-Normal distribution to infer the parameters of the astrophysical model from the statistics of the GWB and show that correctly modeling higher order moments is essential.	result 0.86	0.96	0.31	asserted 1.00
s7	Future PTA analyses should include the statistics of the GWB in their pipelines.	future work 1.00	0.05	0.33	asserted 0.56

What kind of paper the classifier reports this is: method or tool 1.00 · evidence from simulation or computation 0.89.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 1.** FIG. 1. Saddlepoint approximation for the characteristic

**Figure 2.** FIG. 2. Compound poisson Monte Carlo versus saddlepoint

**Figure 3.** FIG. 3. Posterior distribution for the simplistic GWB model

**Figure 5.** FIG. 5. Posterior distributions of population hyperparame

**Figure 6.** FIG. 6. Relative difference between the arcsine and Gaussian

**Figure 7.** FIG. 7. Hellinger distance between a KDE of compound pois

Showing 6 of the paper's 7 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.07686 by the authors.

References

Works cited (48)

[1]

High mass population The high mass population has a characteristic mass of log 10M0 = 9.3. Then, some individual binaries are 8 High mass Low mass log10 ˙n0 -2.9 -1.74 α 0.5 0.5 log10M0 9.3 8.3 β 0.5 0.5 z0 1 1 TABLE I. Parameters used for the high mass and low mass populations. They both give a GWB with an average charac- teristic strainh c = 3×10 −15(f/...

View on Pith
[2]

Low mass population The low mass population has a characteristic mass of log10M0 = 8.3. Then, the individual binaries are less bright and many binaries are required to produce a sig- nal with amplitudeh c(f) = 3×10 −15(f/fyr)−2/3. This tends to Gaussianize the signal and reduce the spectral variance. 0.6 1.2 1.8 2.4 7.2 7.8 8.4 9.0 9.6 log10 0 0.6 1.2 1.8...

View on Pith 2026
[3]

E. S. Phinney, A practical theorem on gravitational wave backgrounds (2001), arXiv:astro-ph/0108028 [astro-ph]

View on Pith Pith review arXiv 2001
[4]

G. Agazieet al., The NANOGrav 15 yr Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational-wave Background, The Astrophysi- cal Journal Letters952, L37 (2023), arXiv:2306.16220 12 [astro-ph.HE]

View on Pith Pith review arXiv 2023
[5]

M. V. Sazhin, Opportunities for detecting ultralong grav- itational waves, Soviet Astronomy22, 36 (1978)

View on Pith 1978
[6]

S. Detweiler, Pulsar timing measurements and the search for gravitational waves, The Astrophysical Journal234, 1100 (1979)

View on Pith 1979
[7]

A. Sesanaet al., Gravitational waves from resolvable massive black hole binary systems and observations with Pulsar Timing Arrays, Monthly Notices of the Royal As- tronomical Society394, 2255 (2009), arXiv:0809.3412 [astro-ph]

View on Pith Pith review arXiv 2009
[8]

A. Sesanaet al., The stochastic gravitational-wave back- ground from massive black hole binary systems: implica- tions for observations with pulsar timing arrays, Monthly Notices of the Royal Astronomical Society390, 192–209 (2008)

View on Pith 2008
[9]

J. Antoniadiset al., The second data release from the european pulsar timing array: Iii. search for gravita- tional wave signals, Astronomy &; Astrophysics678, A50 (2023)

View on Pith 2023
[10]

G. Agazieet al., The nanograv 15 yr data set: Evidence for a gravitational-wave background, The Astrophysical Journal Letters951, L8 (2023)

View on Pith 2023
[11]

H. Xuet al., Searching for the nano-hertz stochastic grav- itational wave background with the chinese pulsar timing array data release i, Research in Astronomy and Astro- physics23, 075024 (2023)

View on Pith 2023
[12]

D. J. Reardonet al., Search for an isotropic gravitational- wave background with the parkes pulsar timing array, The Astrophysical Journal Letters951, L6 (2023)

View on Pith 2023
[13]

M. T. Mileset al., The meerkat pulsar timing array: the first search for gravitational waves with the meerkat ra- dio telescope, Monthly Notices of the Royal Astronomical Society536, 1489–1500 (2024)

View on Pith 2024
[14]

S. Bonoli, D. Izquierdo-Villalba, D. Spinoso, M. Colpi, A. Sesana, M. Polkas, and V. Springel, Constraints on the early growth of massive black holes from pta and jwst with l-galaxiesbh (2025), arXiv:2509.12325 [astro- ph.GA]

View on Pith arXiv 2025
[15]

B. Goncharov, G. Sato-Polito, X. Bi, and M. Zaldar- riaga, A joint optimal search for gravitational waves from resolved and unresolved supermassive binary black holes with pulsar timing arrays (2026), arXiv:2606.18241 [astro-ph.HE]

View on Pith Pith review arXiv 2026
[16]

G. Sato-Polito, M. Zaldarriaga, and E. Quataert, Where are the supermassive black holes measured by ptas?, Phys. Rev. D110, 063020 (2024)

View on Pith 2024
[17]

X. Xue, Z. Pan, and L. Dai, Non-gaussian statis- tics of nanohertz stochastic gravitational waves (2025), arXiv:2409.19516 [astro-ph.CO]

View on Pith Pith review arXiv 2025
[18]

G. Sato-Polito and M. Zaldarriaga, Distribution of the gravitational-wave background from supermassive black holes, Phys. Rev. D111, 023043 (2025)

View on Pith 2025
[19]

Y. Ali-Ha¨ ımoud, A practical theorem on gravitational- wave background statistics (2026), arXiv:2604.19701 [astro-ph.CO]

View on Pith Pith review arXiv 2026
[20]

J. Raidal, J. Urrutia, V. Vaskonen, and H. Veerm¨ ae, The heavy tailed non-gaussianity of the supermas- sive black hole gravitational wave background (2026), arXiv:2604.08506 [astro-ph.CO]

View on Pith Pith review arXiv 2026
[21]

C. Cecchini, J. E. Gammal, G. Franciolini, and M. Pieroni, Are pta measurements sensitive to gravita- tional wave non-gaussianities? (2026), arXiv:2605.05157 [astro-ph.CO]

View on Pith Pith review arXiv 2026
[22]

A. Kuntz, C. Smarra, and M. Vaglio, Looking for non- gaussianity in pulsar timing arrays through the four point correlator (2026), arXiv:2603.12311 [gr-qc]

View on Pith Pith review arXiv 2026
[23]

M. Falxa and A. Sesana, Modeling non-gaussianities in pulsar timing array data analysis using gaussian mixture models, Phys. Rev. D113, 043047 (2026)

View on Pith 2026
[24]

A. Gundersen, R. van Haasteren, N. J. Cornish, M. Val- lisneri, and P. M. Meyers, A new framework for lightning- fast gravitational wave analysis of pulsar timing data (2026), arXiv:2607.06834 [gr-qc]

View on Pith Pith review arXiv 2026
[25]

N. Laal, S. R. Taylor, L. Z. Kelley, J. Simon, K. Gul- tekin, D. Wright, B. Becsy, J. A. Casey-Clyde, S. Chen, A. Cingoranelli, D. J. D’Orazio, E. C. Gardiner, W. G. Lamb, C. Matt, M. S. Siwek, and J. M. Wachter, Deep neural emulation of the supermassive black-hole binary population (2024), arXiv:2411.10519 [astro-ph.IM]

View on Pith Pith review arXiv 2024
[26]

M. Bonetti, A. Franchini, B. G. Galuzzi, and A. Sesana, Neural networks unveiling the properties of gravitational wave background from supermassive black hole binaries, Astronomy &; Astrophysics687, A42 (2024)

View on Pith 2024
[27]

W. G. Lamb, J. M. Wachter, A. Mitridate, S. C. Sardesai, B. B´ ecsy, E. L. Hagen, S. R. Taylor, and L. Z. Kelley, Finite populations and finite time: The non-gaussianity of a gravitational wave background, Physical Review D 113, 10.1103/96zk-qtck (2026)

View on Pith doi:10.1103/96zk-qtck 2026
[28]

M. Crisostomiet al., Beyond diagonal approximations: improved covariance modeling for pulsar timing array data analysis (2025), arXiv:2506.13866 [astro-ph.IM]

View on Pith Pith review arXiv 2025
[29]

H. Quelquejay Leclere, K. Li, M. Volonteri, S. Babak, R. S. Beckmann, Y. Dubois, C. Laigle, and N. A. Webb, The multi-messenger view of pulsar timing array black holes with the horizon-agn simulation, Astronomy &; As- trophysics705, A56 (2026)

View on Pith 2026
[30]

S. R. Taylor, R. van Haasteren, and A. Sesana, From bright binaries to bumpy backgrounds: Mapping realistic gravitational wave skies with pulsar-timing arrays, Phys- ical Review D102, 10.1103/physrevd.102.084039 (2020)

View on Pith doi:10.1103/physrevd.102.084039 2020
[31]

R. W. Hellings and G. S. Downs, Upper limits on the isotropic gravitational radiation background from pulsar timing analysis., The Astrophysical Journal Letters265, L39 (1983)

View on Pith 1983
[32]

O. Barndorff-Nielsen, J. Kent, and M. Sørensen, Nor- mal variance-mean mixtures and z-distributions, Inter- national Statistical Review50, 145 (1982)

View on Pith 1982
[33]

Y. Yu, On normal variance-mean mixtures (2011), arXiv:1106.2333 [math.ST]

View on Pith Pith review arXiv 2011
[34]

N. Karnesis, A. Sasli, R. Buscicchio, and N. Stergioulas, Characterization of non-gaussian stochastic signals with heavier-tailed likelihoods, Phys. Rev. D111, 022005 (2025)

View on Pith 2025
[35]

H. E. Daniels, Saddlepoint Approximations in Statistics, The Annals of Mathematical Statistics25, 631 (1954)

View on Pith 1954
[36]

G. Sato-Polito and M. Kamionkowski, Exploring the spectrum of stochastic gravitational-wave anisotropies with pulsar timing arrays (2025), arXiv:2305.05690 [astro-ph.CO]

View on Pith Pith review arXiv 2025
[37]

H. Middleton, W. Del Pozzo, W. M. Farr, A. Sesana, and A. Vecchio, Astrophysical constraints on massive black hole binary evolution from pulsar timing arrays, Monthly Notices of the Royal Astronomical Society: Letters455, L72–L76 (2015)

View on Pith 2015
[38]

H. Quelquejay Leclere, Probing past mergers of super- massive black holes with pulsar timing arrays: The role of pulsar terms, Physical Review D114, 10.1103/9wq2- 13 gqkh (2026)

View on Pith doi:10.1103/9wq2- 2026
[39]

R. van Haasteren and M. Vallisneri, New advances in the gaussian-process approach to pulsar-timing data analy- sis, Physical Review D90, 10.1103/physrevd.90.104012 (2014)

View on Pith Pith review doi:10.1103/physrevd.90.104012 2014
[40]

R. van Haasteren and Y. Levin, Understanding and analysing time-correlated stochastic signals in pulsar tim- ing, Monthly Notices of the Royal Astronomical Society 428, 1147–1159 (2012)

View on Pith 2012
[41]

J. Antoniadiset al., The second data release from the european pulsar timing array: I. the dataset and timing analysis, Astronomy &; Astrophysics678, A48 (2023)

View on Pith 2023
[42]

R. van Haasteren and M. Vallisneri, Low-rank approxi- mations for large stationary covariance matrices, as used in the bayesian and generalized-least-squares analysis of pulsar-timing data, Monthly Notices of the Royal Astro- nomical Society446, 1170–1174 (2014)

View on Pith 2014
[43]

J. D. Romano and N. J. Cornish, Detection meth- ods for stochastic gravitational-wave backgrounds: a unified treatment, Living Reviews in Relativity20, 10.1007/s41114-017-0004-1 (2017)

View on Pith doi:10.1007/s41114-017-0004-1 2017
[44]

H. Quelquejay Leclereet al., Practical approaches to ana- lyzing pta data: Cosmic strings with six pulsars, Physical Review D108, 10.1103/physrevd.108.123527 (2023)

View on Pith doi:10.1103/physrevd.108.123527 2023
[45]

W. G. Lambet al., Rapid refitting techniques for bayesian spectral characterization of the gravitational wave background using pulsar timing arrays, Physical Review D108, 10.1103/physrevd.108.103019 (2023)

View on Pith doi:10.1103/physrevd.108.103019 2023
[46]

D. Phan, N. Pradhan, and M. Jankowiak, Compos- able effects for flexible and accelerated probabilistic pro- gramming in numpyro, arXiv preprint arXiv:1912.11554 (2019)

View on Pith Pith review arXiv 2019
[47]

M. Vallisneri, P. M. Meyers, D. Wright, A. D. Johnson, J. G. Baier, and R. van Haasteren, nanograv/discovery: 0.5.1 (2025)

View on Pith 2025
[48]

M. Anholm, S. Ballmer, J. D. E. Creighton, L. R. Price, and X. Siemens, Optimal strategies for gravitational wave stochastic background searches in pulsar timing data, Physical Review D79, 10.1103/physrevd.79.084030 (2009)

View on Pith doi:10.1103/physrevd.79.084030 2009

Forward citations

Cited by 2 Pith papers

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

Gravitational-wave background from supermassive black holes without merger trees: Tidally regulated mergers and direct inference from the NANOGrav 15 yr data
astro-ph.CO 2026-10 conditional novelty 7.0 of 10

A merger-tree-free population model coupling tidal stripping to dynamical friction directly fits NANOGrav 15 yr data and predicts a stalled, wandering black hole population below a critical mass ratio.
Joint Angular-Power Distributions of Interfering Nanohertz Gravitational Waves
astro-ph.CO 2026-10 conditional novelty 5.0 of 10

Closed-form statistics for the realized angular correlation of interfering nanohertz gravitational waves: Gamma for many weak sources, non-central chi-squared with a bright source, and a compact multivariate lognormal...

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.

[1] [1]

High mass population The high mass population has a characteristic mass of log 10M0 = 9.3. Then, some individual binaries are 8 High mass Low mass log10 ˙n0 -2.9 -1.74 α 0.5 0.5 log10M0 9.3 8.3 β 0.5 0.5 z0 1 1 TABLE I. Parameters used for the high mass and low mass populations. They both give a GWB with an average charac- teristic strainh c = 3×10 −15(f/...

View on Pith

[2] [2]

Low mass population The low mass population has a characteristic mass of log10M0 = 8.3. Then, the individual binaries are less bright and many binaries are required to produce a sig- nal with amplitudeh c(f) = 3×10 −15(f/fyr)−2/3. This tends to Gaussianize the signal and reduce the spectral variance. 0.6 1.2 1.8 2.4 7.2 7.8 8.4 9.0 9.6 log10 0 0.6 1.2 1.8...

View on Pith 2026

[3] [3]

E. S. Phinney, A practical theorem on gravitational wave backgrounds (2001), arXiv:astro-ph/0108028 [astro-ph]

View on Pith Pith review arXiv 2001

[4] [4]

G. Agazieet al., The NANOGrav 15 yr Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational-wave Background, The Astrophysi- cal Journal Letters952, L37 (2023), arXiv:2306.16220 12 [astro-ph.HE]

View on Pith Pith review arXiv 2023

[5] [5]

M. V. Sazhin, Opportunities for detecting ultralong grav- itational waves, Soviet Astronomy22, 36 (1978)

View on Pith 1978

[6] [6]

S. Detweiler, Pulsar timing measurements and the search for gravitational waves, The Astrophysical Journal234, 1100 (1979)

View on Pith 1979

[7] [7]

A. Sesanaet al., Gravitational waves from resolvable massive black hole binary systems and observations with Pulsar Timing Arrays, Monthly Notices of the Royal As- tronomical Society394, 2255 (2009), arXiv:0809.3412 [astro-ph]

View on Pith Pith review arXiv 2009

[8] [8]

A. Sesanaet al., The stochastic gravitational-wave back- ground from massive black hole binary systems: implica- tions for observations with pulsar timing arrays, Monthly Notices of the Royal Astronomical Society390, 192–209 (2008)

View on Pith 2008

[9] [9]

J. Antoniadiset al., The second data release from the european pulsar timing array: Iii. search for gravita- tional wave signals, Astronomy &; Astrophysics678, A50 (2023)

View on Pith 2023

[10] [10]

G. Agazieet al., The nanograv 15 yr data set: Evidence for a gravitational-wave background, The Astrophysical Journal Letters951, L8 (2023)

View on Pith 2023

[11] [11]

H. Xuet al., Searching for the nano-hertz stochastic grav- itational wave background with the chinese pulsar timing array data release i, Research in Astronomy and Astro- physics23, 075024 (2023)

View on Pith 2023

[12] [12]

D. J. Reardonet al., Search for an isotropic gravitational- wave background with the parkes pulsar timing array, The Astrophysical Journal Letters951, L6 (2023)

View on Pith 2023

[13] [13]

M. T. Mileset al., The meerkat pulsar timing array: the first search for gravitational waves with the meerkat ra- dio telescope, Monthly Notices of the Royal Astronomical Society536, 1489–1500 (2024)

View on Pith 2024

[14] [14]

S. Bonoli, D. Izquierdo-Villalba, D. Spinoso, M. Colpi, A. Sesana, M. Polkas, and V. Springel, Constraints on the early growth of massive black holes from pta and jwst with l-galaxiesbh (2025), arXiv:2509.12325 [astro- ph.GA]

View on Pith arXiv 2025

[15] [15]

B. Goncharov, G. Sato-Polito, X. Bi, and M. Zaldar- riaga, A joint optimal search for gravitational waves from resolved and unresolved supermassive binary black holes with pulsar timing arrays (2026), arXiv:2606.18241 [astro-ph.HE]

View on Pith Pith review arXiv 2026

[16] [16]

G. Sato-Polito, M. Zaldarriaga, and E. Quataert, Where are the supermassive black holes measured by ptas?, Phys. Rev. D110, 063020 (2024)

View on Pith 2024

[17] [17]

X. Xue, Z. Pan, and L. Dai, Non-gaussian statis- tics of nanohertz stochastic gravitational waves (2025), arXiv:2409.19516 [astro-ph.CO]

View on Pith Pith review arXiv 2025

[18] [18]

G. Sato-Polito and M. Zaldarriaga, Distribution of the gravitational-wave background from supermassive black holes, Phys. Rev. D111, 023043 (2025)

View on Pith 2025

[19] [19]

Y. Ali-Ha¨ ımoud, A practical theorem on gravitational- wave background statistics (2026), arXiv:2604.19701 [astro-ph.CO]

View on Pith Pith review arXiv 2026

[20] [20]

J. Raidal, J. Urrutia, V. Vaskonen, and H. Veerm¨ ae, The heavy tailed non-gaussianity of the supermas- sive black hole gravitational wave background (2026), arXiv:2604.08506 [astro-ph.CO]

View on Pith Pith review arXiv 2026

[21] [21]

C. Cecchini, J. E. Gammal, G. Franciolini, and M. Pieroni, Are pta measurements sensitive to gravita- tional wave non-gaussianities? (2026), arXiv:2605.05157 [astro-ph.CO]

View on Pith Pith review arXiv 2026

[22] [22]

A. Kuntz, C. Smarra, and M. Vaglio, Looking for non- gaussianity in pulsar timing arrays through the four point correlator (2026), arXiv:2603.12311 [gr-qc]

View on Pith Pith review arXiv 2026

[23] [23]

M. Falxa and A. Sesana, Modeling non-gaussianities in pulsar timing array data analysis using gaussian mixture models, Phys. Rev. D113, 043047 (2026)

View on Pith 2026

[24] [24]

A. Gundersen, R. van Haasteren, N. J. Cornish, M. Val- lisneri, and P. M. Meyers, A new framework for lightning- fast gravitational wave analysis of pulsar timing data (2026), arXiv:2607.06834 [gr-qc]

View on Pith Pith review arXiv 2026

[25] [25]

N. Laal, S. R. Taylor, L. Z. Kelley, J. Simon, K. Gul- tekin, D. Wright, B. Becsy, J. A. Casey-Clyde, S. Chen, A. Cingoranelli, D. J. D’Orazio, E. C. Gardiner, W. G. Lamb, C. Matt, M. S. Siwek, and J. M. Wachter, Deep neural emulation of the supermassive black-hole binary population (2024), arXiv:2411.10519 [astro-ph.IM]

View on Pith Pith review arXiv 2024

[26] [26]

M. Bonetti, A. Franchini, B. G. Galuzzi, and A. Sesana, Neural networks unveiling the properties of gravitational wave background from supermassive black hole binaries, Astronomy &; Astrophysics687, A42 (2024)

View on Pith 2024

[27] [27]

W. G. Lamb, J. M. Wachter, A. Mitridate, S. C. Sardesai, B. B´ ecsy, E. L. Hagen, S. R. Taylor, and L. Z. Kelley, Finite populations and finite time: The non-gaussianity of a gravitational wave background, Physical Review D 113, 10.1103/96zk-qtck (2026)

View on Pith doi:10.1103/96zk-qtck 2026

[28] [28]

M. Crisostomiet al., Beyond diagonal approximations: improved covariance modeling for pulsar timing array data analysis (2025), arXiv:2506.13866 [astro-ph.IM]

View on Pith Pith review arXiv 2025

[29] [29]

H. Quelquejay Leclere, K. Li, M. Volonteri, S. Babak, R. S. Beckmann, Y. Dubois, C. Laigle, and N. A. Webb, The multi-messenger view of pulsar timing array black holes with the horizon-agn simulation, Astronomy &; As- trophysics705, A56 (2026)

View on Pith 2026

[30] [30]

S. R. Taylor, R. van Haasteren, and A. Sesana, From bright binaries to bumpy backgrounds: Mapping realistic gravitational wave skies with pulsar-timing arrays, Phys- ical Review D102, 10.1103/physrevd.102.084039 (2020)

View on Pith doi:10.1103/physrevd.102.084039 2020

[31] [31]

R. W. Hellings and G. S. Downs, Upper limits on the isotropic gravitational radiation background from pulsar timing analysis., The Astrophysical Journal Letters265, L39 (1983)

View on Pith 1983

[32] [32]

O. Barndorff-Nielsen, J. Kent, and M. Sørensen, Nor- mal variance-mean mixtures and z-distributions, Inter- national Statistical Review50, 145 (1982)

View on Pith 1982

[33] [33]

Y. Yu, On normal variance-mean mixtures (2011), arXiv:1106.2333 [math.ST]

View on Pith Pith review arXiv 2011

[34] [34]

N. Karnesis, A. Sasli, R. Buscicchio, and N. Stergioulas, Characterization of non-gaussian stochastic signals with heavier-tailed likelihoods, Phys. Rev. D111, 022005 (2025)

View on Pith 2025

[35] [35]

H. E. Daniels, Saddlepoint Approximations in Statistics, The Annals of Mathematical Statistics25, 631 (1954)

View on Pith 1954

[36] [36]

G. Sato-Polito and M. Kamionkowski, Exploring the spectrum of stochastic gravitational-wave anisotropies with pulsar timing arrays (2025), arXiv:2305.05690 [astro-ph.CO]

View on Pith Pith review arXiv 2025

[37] [37]

H. Middleton, W. Del Pozzo, W. M. Farr, A. Sesana, and A. Vecchio, Astrophysical constraints on massive black hole binary evolution from pulsar timing arrays, Monthly Notices of the Royal Astronomical Society: Letters455, L72–L76 (2015)

View on Pith 2015

[38] [38]

H. Quelquejay Leclere, Probing past mergers of super- massive black holes with pulsar timing arrays: The role of pulsar terms, Physical Review D114, 10.1103/9wq2- 13 gqkh (2026)

View on Pith doi:10.1103/9wq2- 2026

[39] [39]

R. van Haasteren and M. Vallisneri, New advances in the gaussian-process approach to pulsar-timing data analy- sis, Physical Review D90, 10.1103/physrevd.90.104012 (2014)

View on Pith Pith review doi:10.1103/physrevd.90.104012 2014

[40] [40]

R. van Haasteren and Y. Levin, Understanding and analysing time-correlated stochastic signals in pulsar tim- ing, Monthly Notices of the Royal Astronomical Society 428, 1147–1159 (2012)

View on Pith 2012

[41] [41]

J. Antoniadiset al., The second data release from the european pulsar timing array: I. the dataset and timing analysis, Astronomy &; Astrophysics678, A48 (2023)

View on Pith 2023

[42] [42]

R. van Haasteren and M. Vallisneri, Low-rank approxi- mations for large stationary covariance matrices, as used in the bayesian and generalized-least-squares analysis of pulsar-timing data, Monthly Notices of the Royal Astro- nomical Society446, 1170–1174 (2014)

View on Pith 2014

[43] [43]

J. D. Romano and N. J. Cornish, Detection meth- ods for stochastic gravitational-wave backgrounds: a unified treatment, Living Reviews in Relativity20, 10.1007/s41114-017-0004-1 (2017)

View on Pith doi:10.1007/s41114-017-0004-1 2017

[44] [44]

H. Quelquejay Leclereet al., Practical approaches to ana- lyzing pta data: Cosmic strings with six pulsars, Physical Review D108, 10.1103/physrevd.108.123527 (2023)

View on Pith doi:10.1103/physrevd.108.123527 2023

[45] [45]

W. G. Lambet al., Rapid refitting techniques for bayesian spectral characterization of the gravitational wave background using pulsar timing arrays, Physical Review D108, 10.1103/physrevd.108.103019 (2023)

View on Pith doi:10.1103/physrevd.108.103019 2023

[46] [46]

D. Phan, N. Pradhan, and M. Jankowiak, Compos- able effects for flexible and accelerated probabilistic pro- gramming in numpyro, arXiv preprint arXiv:1912.11554 (2019)

View on Pith Pith review arXiv 2019

[47] [47]

M. Vallisneri, P. M. Meyers, D. Wright, A. D. Johnson, J. G. Baier, and R. van Haasteren, nanograv/discovery: 0.5.1 (2025)

View on Pith 2025

[48] [48]

M. Anholm, S. Ballmer, J. D. E. Creighton, L. R. Price, and X. Siemens, Optimal strategies for gravitational wave stochastic background searches in pulsar timing data, Physical Review D79, 10.1103/physrevd.79.084030 (2009)

View on Pith doi:10.1103/physrevd.79.084030 2009