Pith. sign in

REVIEW 2 major objections 5 minor 29 references

LISA's expected SNR for a stochastic gravitational-wave background is always bounded by the square root of observing time times bandwidth, never exceeding roughly 10,000 for typical mission parameters.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-10 07:20 UTC pith:AA5NGBYN

load-bearing objection Clean, usable fix to the standard LISA auto-correlation SNR that people have been over-applying to loud signals, plus a hard ceiling and ready contour data. the 2 major comments →

arxiv 2607.08445 v1 pith:AA5NGBYN submitted 2026-07-09 gr-qc astro-ph.IMhep-ph

Signal-to-Noise Ratio Contours for LISA

Kai Schmitz , Joseph D. Romano This is my paper
classification gr-qc astro-ph.IMhep-ph
keywords LISAstochastic gravitational-wave backgroundsignal-to-noise ratiopower-law-integrated sensitivitygravitational-wave self-noiseauto-correlation statisticSNR contours
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Most forecasts of LISA's reach for a stochastic gravitational-wave background use a signal-to-noise formula that is valid only when the signal is much weaker than the detector noise. This paper derives the correct optimal SNR for an auto-correlation measurement at any signal strength: the total noise is the sum of instrumental noise and the signal itself (gravitational-wave self-noise). The resulting expression implies a hard upper bound SNR_max = sqrt(T_obs (f_max - f_min)), which is only a few thousand for a year of data. The authors turn the formula into families of power-law-integrated sensitivity curves at fixed SNR, producing LISA SNR contours that bend, cut off at low frequency, and are bounded above by roughly e times the ordinary strain-noise curve. For any signal strong enough to approach that ceiling, the old weak-signal formula is inaccurate and must be replaced by the new one.

Core claim

The expected optimal SNR of a LISA auto-correlation measurement at arbitrary signal strength is SNR = [T_obs ∫ df (h^{2}Ω_signal / (h^{2}Ω_noise + h^{2}Ω_signal))^{2} ]^{1/2} between f_min and f_max. Consequently the SNR can never exceed sqrt(T_obs (f_max - f_min)) ≲ 10^4 for typical mission parameters, and the associated power-law-integrated sensitivity curves form a set of SNR contours that terminate at a low-frequency cutoff and lie approximately a factor of e above the ordinary LISA strain-noise curve.

What carries the argument

The matched-filter auto-correlation statistic whose variance under the signal hypothesis includes both detector noise and gravitational-wave self-noise; maximizing the ratio of mean to standard deviation yields the general SNR formula that reduces to the familiar weak-signal expression only when the signal contribution can be dropped.

Load-bearing premise

The analysis assumes perfect prior knowledge of the instrumental noise spectrum (taken from the T channel) so that the pure-noise expectation can be subtracted exactly when the mean signal is formed.

What would settle it

Compute the exact SNR for any strong, scale-invariant or power-law spectrum both with the old weak-signal formula and with the new expression that retains self-noise; if the numerical values agree beyond the claimed transition near h^{2}Ω ~ 10^{-11}, or if an actual LISA auto-correlation measurement exceeds the analytic bound sqrt(T_obs (f_max - f_min)), the central claim fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Any published LISA SNR larger than a few thousand for a stochastic background is necessarily an artifact of the weak-signal approximation and must be recomputed.
  • Power-law-integrated sensitivity curves at high SNR levels are no longer simple vertical shifts of the SNR = 1 curve; each contour carries independent information and terminates at a low-frequency cutoff.
  • For signals that dominate the noise below some frequency f0, the achievable SNR is bounded by sqrt(T_obs (f0 - f_min)), independent of how large the amplitude becomes.
  • The same self-noise-corrected SNR formula and contour construction can be repeated for other detectors or networks once their noise PSDs are known.

Where Pith is reading between the lines

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

  • Cosmological models that predict millihertz backgrounds near or above the LISA noise floor (certain cosmic-string or phase-transition scenarios) will have their claimed detection significance systematically overstated if the literature continues to quote the weak-signal formula.
  • The hard bandwidth-time bound suggests that multi-year LISA runs improve SNR only as the square root of time even for arbitrarily loud signals, so mission extensions cannot push the auto-correlation SNR into the 10^5 regime.
  • Because residual uncertainty in the noise PSD would inflate the variance, any practical noise-estimation residual will lower the realized SNR below the ideal curves presented here, making the contours an optimistic upper envelope.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper derives the expected optimal SNR for a LISA auto-correlation search for a stochastic gravitational-wave background that remains valid at arbitrary signal strength, SNR = [T_obs ∫_{f_min}^{f_max} df (h²Ω_signal/(h²Ω_noise + h²Ω_signal))²]^{1/2} (Eqs. 1.4, 3.32). The derivation retains the GWB contribution (self-noise) in the variance of the matched-filter auto-correlation statistic under Gaussian assumptions, recovers the standard weak-signal limit, and implies a hard upper bound SNR ≤ SNR_max = √[T_obs(f_max − f_min)] ≲ 10^4 for typical mission parameters. The authors recast this result as generalized power-law-integrated SNR contours, identify a low-frequency cutoff contour that is approximately e times the LISA strain noise curve, release numerical contour data on Zenodo, and illustrate the formalism with Nambu–Goto cosmic-string spectra.

Significance. If correct, the result is a useful and overdue correction to a formula that is widely used beyond its domain of validity in LISA GWB forecasts. The derivation is standard (Isserlis/Wick four-point expansion plus Cauchy–Schwarz matched filtering) once the total PSD is kept in the variance; the absolute bound SNR_max and the approximate e-factor relation for the cutoff contour are clean, falsifiable consequences. Releasing the contour data on Zenodo and demonstrating graphical SNR evaluation for non-power-law cosmic-string spectra strengthens the practical value for the LISA cosmology community. The work cleanly fills the auto-correlation gap left by the network cross-correlation treatment of Allen & Romano (1999).

major comments (2)
  1. Abstract and §1: the claim that a large number of literature SNR values are “inaccurate” needs a sharper detection-vs-estimation distinction. Eqs. (1.2)–(1.3) correctly separate null-hypothesis variance (detection significance) from signal-hypothesis variance (estimation). The weak-signal formula remains the right figure of merit for false-alarm rates near threshold; the new formula is required when the same expression is used as a figure of merit for strong signals or for parameter-estimation precision. The abstract’s wording risks being read as invalidating detection SNRs; a short clarifying sentence in the abstract and conclusions would prevent misapplication.
  2. §3.2, Eqs. (3.9)–(3.10): the optimality of Eq. (3.32) rests on exact subtraction of the noise-only expectation using perfect knowledge of the instrumental PSD (assumed available from the T channel). This is conventional in LISA TDI forecasts and is stated, but residual noise uncertainty is the weakest modeling assumption. A brief quantitative remark—e.g., that fractional noise misestimation δN/N degrades the effective SNR at leading order without changing the functional form of the self-noise correction or the existence of SNR_max—would make the domain of the central claim more transparent without altering the derivation.
minor comments (5)
  1. §2.3: “unresolved GNC background” is a typo for GCN (galactic confusion noise).
  2. §4.2: “in turns of a power law” should read “in terms of a power law”.
  3. Fig. 2 and Eqs. (3.36), (3.40): the ISR analytic estimate (~3500) and numerical fit (~3200) differ by O(10 %); a one-sentence note that the Heaviside/f^5 approximation is only schematic would help readers who compare the curves.
  4. §2.2–2.3 and Table 1: T_obs is fixed to 1 yr for the GCN fit while SNR_max is quoted for “typical mission parameters.” A short remark that contours for other T_obs follow by rescaling (or a pointer to the Zenodo files) would improve usability.
  5. References: the NANOGrav self-noise discussion [18] and the partial implementation in [19] are appropriately cited; a parenthetical note that the present auto-correlation variance (Eq. 3.21) had not previously been written down would help priority claims without overstatement.

Circularity Check

0 steps flagged

No significant circularity: the generalized SNR follows from the Gaussian four-point expansion and matched-filter Cauchy–Schwarz once total PSD is retained in the variance.

full rationale

The paper’s central result (Eqs. 1.4 / 3.32) is obtained by writing the auto-correlation statistic S, evaluating its mean µ under the signal hypothesis after exact subtraction of the known noise floor, expanding the four-point correlator with Isserlis’s theorem, retaining both noise and signal contributions in the variance, introducing the positive-definite inner product weighted by P², and maximizing via Cauchy–Schwarz. Every algebraic step is self-contained; no parameter is fitted to data and then re-used as a prediction, and no uniqueness theorem or ansatz is imported from the authors’ prior work to force the functional form. The absolute bound SNR_max = sqrt(T_obs(f_max − f_min)) is an immediate consequence of the integrand being at most unity. The approximate relation h²Ω_cut ∼ e h²Ω_noise is derived as an asymptotic expansion of a hypergeometric integral for steep power-law signals and is presented as a numerical observation, not as an input. Self-citations (Thrane & Romano 2013, Allen & Romano 1999, Schmitz 2020/2021) supply standard weak-signal formulas and noise models that are then generalized; they are not load-bearing for the new arbitrary-strength expression. The derivation is therefore independent of its own conclusions.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The central claim rests on standard Gaussian statistics plus a handful of domain assumptions about LISA noise and the availability of a pure-noise reference channel. No free parameters are fitted to external data; the numerical contours are pure evaluations of the derived integral. No new physical entities are postulated.

free parameters (3)
  • T_obs = 1 yr
    Observing time fixed by hand to 1 yr for all numerical examples and contour plots; scales the overall SNR as sqrt(T_obs).
  • f_min, f_max = 10^{-5} Hz, 1 Hz
    Integration limits fixed by hand to 10^{-5} Hz and 1 Hz; enter SNR_max directly.
  • galactic-confusion-noise fit coefficients = table values for 1 yr
    α, β, κ, γ, f_knee taken from Cornish & Robson 2017 table for T_obs = 1 yr; affect the detailed shape of h²Ω_noise but not the existence of the SNR bound.
axioms (4)
  • domain assumption Both the stochastic GW signal and the instrumental noise are stationary Gaussian random processes, so that all higher-order correlators factor via Isserlis’ theorem.
    Invoked throughout Sec. 3.2 to reduce the four-point function to products of two-point functions.
  • domain assumption The instrumental noise PSD is known perfectly a priori (e.g., from the TDI T channel) and can be subtracted exactly when forming the mean µ.
    Stated explicitly after Eq. 3.9; without it the optimal filter and the SNR formula change.
  • standard math The filter Q is real and even, and the large-T_obs approximation δ_Tobs² ≈ δ_Tobs δ holds.
    Used to obtain the compact variance expression Eq. 3.21.
  • domain assumption LISA can be idealized as two equal-arm Michelson interferometers (A and E channels) with the analytic response function of Zhang et al. 2020.
    Sec. 2.2; unequal-arm and flexing effects are acknowledged but neglected.

pith-pipeline@v1.1.0-grok45 · 25266 in / 2605 out tokens · 27929 ms · 2026-07-10T07:20:21.316169+00:00 · methodology

0 comments
read the original abstract

The Laser Interferometer Space Antenna (LISA) will search for a stochastic gravitational-wave (GW) background at millihertz frequencies, from both astrophysical and cosmological sources, and thereby open a new chapter in GW astronomy. In the literature, LISA's sensitivity to prospective GW background (GWB) signals is often quantified in terms of an expected signal-to-noise ratio (SNR) assuming perfect knowledge of the detector noise. The commonly employed expression for the SNR is, however, valid only in the limit of a weak GWB signal, which renders a large number of SNR values reported in the literature inaccurate. In this paper, we address this issue by deriving for the first time an expression for the expected optimal SNR of a LISA auto-correlation measurement that is valid at arbitrary signal strength. Based on our generalized expression, we conclude that LISA data worth an observing time of T_obs across the frequency band from f_min to f_max will never yield an SNR in excess of SNR_max = sqrt(T_obs(f_max-f_min)), which evaluates to SNR_max <~ 10^4 for typical mission parameters. We illustrate our findings in terms of generalized power-law-integrated (PLI) sensitivity curves at different SNR levels, i.e., LISA SNR contour lines in plots of the GW energy-density power spectrum. In contrast to earlier work on PLI sensitivity curves, we notably find that the LISA SNR contours are bounded from above, approximately by the LISA strain noise curve multiplied by a factor of Euler's number e. For GWB signals not much weaker than this range, the expected SNR for a LISA auto-correlation measurement needs to be evaluated based on our new expression. Our numerical results for the LISA SNR contours are available on Zenodo [https://doi.org/10.5281/zenodo.21275527].

discussion (0)

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

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages · 27 internal anchors

  1. [1]

    Laser Interferometer Space Antenna

    K. Schmitz and J. D. Romano,LISA SNR Contours, Zenodo dataset (2026), https://doi.org/10.5281/zenodo.21275527. [2]LISA, P. Amaro-Seoane et al.,Laser Interferometer Space Antenna,arXiv:1702.00786. – 22 –

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.5281/zenodo.21275527 2026

  2. [2]

    The Laser Interferometer Space Antenna: Unveiling the Millihertz Gravitational Wave Sky

    J. Baker et al.,The Laser Interferometer Space Antenna: Unveiling the Millihertz Gravitational Wave Sky,arXiv:1907.06482

    work page internal anchor Pith review Pith/arXiv arXiv 1907

  3. [3]

    Astrophysics with the Laser Interferometer Space Antenna

    M. Bailes et al.,Gravitational-wave physics and astronomy in the 2020s and 2030s,Nature Rev. Phys.3(2021), no. 5 344–366. [5]LISA, P. A. Seoane et al.,Astrophysics with the Laser Interferometer Space Antenna,Living Rev. Rel.26(2023), no. 1 2, [arXiv:2203.06016]. [6]LISA Cosmology Working Group, P. Auclair et al.,Cosmology with the Laser Interferometer Spa...

    work page internal anchor Pith review Pith/arXiv arXiv 2021

  4. [4]

    Z.-C. Chen, F. Huang, and Q.-G. Huang,Stochastic Gravitational-wave Background from Binary Black Holes and Binary Neutron Stars and Implications for LISA,Astrophys. J.871 (2019), no. 1 97, [arXiv:1809.10360]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  5. [5]

    Stochastic gravitational wave background from stellar origin binary black holes in LISA

    S. Babak, C. Caprini, D. G. Figueroa, N. Karnesis, P. Marcoccia, G. Nardini, M. Pieroni, A. Ricciardone, A. Sesana, and J. Torrado,Stochastic gravitational wave background from stellar origin binary black holes in LISA,JCAP08(2023) 034, [arXiv:2304.06368]. [9]LISA Cosmology Working Group, C. Caprini, R. Jinno, M. Lewicki, E. Madge, M. Merchand, G. Nardini...

    work page internal anchor Pith review Pith/arXiv arXiv 2023

  6. [6]

    The astrophysical gravitational wave stochastic background

    T. Regimbau,The astrophysical gravitational wave stochastic background,Res. Astron. Astrophys.11(2011) 369–390, [arXiv:1101.2762]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  7. [7]

    Cosmological Backgrounds of Gravitational Waves

    C. Caprini and D. G. Figueroa,Cosmological backgrounds of gravitational waves,Class. Quant. Grav.35(2018), no. 16 163001, [arXiv:1801.04268]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  8. [8]

    Detecting a stochastic background of gravitational radiation: Signal processing strategies and sensitivities

    B. Allen and J. D. Romano,Detecting a stochastic background of gravitational radiation: Signal processing strategies and sensitivities,Phys. Rev. D59(1999) 102001, [gr-qc/9710117]

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  9. [9]

    Gravitational Wave Experiments and Early Universe Cosmology

    M. Maggiore,Gravitational wave experiments and early universe cosmology,Phys. Rept.331 (2000) 283–367, [gr-qc/9909001]

    work page internal anchor Pith review Pith/arXiv arXiv 2000

  10. [10]

    Sensitivity curves for searches for gravitational-wave backgrounds

    E. Thrane and J. D. Romano,Sensitivity curves for searches for gravitational-wave backgrounds,Phys. Rev. D88(2013), no. 12 124032, [arXiv:1310.5300]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  11. [11]

    New Sensitivity Curves for Gravitational-Wave Signals from Cosmological Phase Transitions

    K. Schmitz,New Sensitivity Curves for Gravitational-Wave Signals from Cosmological Phase Transitions,JHEP01(2021) 097, [arXiv:2002.04615]. [18]NANOGrav, G. Agazie et al.,The NANOGrav 15 yr Data Set: Detector Characterization and Noise Budget,Astrophys. J. Lett.951(2023), no. 1 L10, [arXiv:2306.16218]

    work page internal anchor Pith review Pith/arXiv arXiv 2021

  12. [12]

    Primordial gravitational waves, precisely: The role of thermodynamics in the Standard Model

    K. Saikawa and S. Shirai,Primordial gravitational waves, precisely: The role of thermodynamics in the Standard Model,JCAP05(2018) 035, [arXiv:1803.01038]. – 23 –

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  13. [13]

    J. M. Wachter, K. D. Olum, J. J. Blanco-Pillado, V. R. Gade, and K. Sivakumar,Numerical gravitational backreaction on cosmic string loops from simulations,Phys. Rev. D113(2026), no. 4 043521, [arXiv:2411.10366]

    work page arXiv 2026

  14. [14]

    J. M. Wachter, K. D. Olum, and J. J. Blanco-Pillado,More accurate gravitational wave backgrounds from cosmic strings,Phys. Rev. D113(2026), no. 8 083533, [arXiv:2411.16590]

    work page arXiv 2026

  15. [15]

    LISA Sensitivity to Gravitational Waves from Sound Waves

    K. Schmitz,LISA Sensitivity to Gravitational Waves from Sound Waves,Symmetry12 (2020), no. 9 1477, [arXiv:2005.10789]

    work page internal anchor Pith review Pith/arXiv arXiv 2020

  16. [16]

    The construction and use of LISA sensitivity curves

    T. Robson, N. J. Cornish, and C. Liu,The construction and use of LISA sensitivity curves, Class. Quant. Grav.36(2019), no. 10 105011, [arXiv:1803.01944]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  17. [17]

    M. C. Edwards, P. Maturana-Russel, R. Meyer, J. Gair, N. Korsakova, and N. Christensen, Identifying and Addressing Nonstationary LISA Noise,Phys. Rev. D102(2020), no. 8 084062, [arXiv:2004.07515]

    work page internal anchor Pith review Pith/arXiv arXiv 2020

  18. [18]

    Revisiting time delay interferometry for unequal-arm LISA and TAIJI

    G. Wang and W.-T. Ni,Revisiting time delay interferometry for unequal-arm LISA and TAIJI,Phys. Scripta98(2023), no. 7 075005, [arXiv:2008.05812]

    work page internal anchor Pith review Pith/arXiv arXiv 2023

  19. [19]

    TDI-infinity: time-delay interferometry without delays

    M. Vallisneri, J.-B. Bayle, S. Babak, and A. Petiteau,Time-delay interferometry without delays,Phys. Rev. D103(2021), no. 8 082001, [arXiv:2008.12343]

    work page internal anchor Pith review Pith/arXiv arXiv 2021

  20. [20]

    Construction of Sensitivity Curves for Dynamic LISA and Taiji

    H.-Y. Shi and Y. Tang,Construction of Sensitivity Curves for Dynamic LISA and Taiji, arXiv:2606.28071

    work page internal anchor Pith review Pith/arXiv arXiv

  21. [21]

    S. L. Larson, W. A. Hiscock, and R. W. Hellings,Sensitivity curves for spaceborne gravitational wave interferometers,Phys. Rev. D62(2000) 062001, [gr-qc/9909080]

    work page internal anchor Pith review Pith/arXiv arXiv 2000

  22. [22]

    Frequency response of space-based interferometric gravitational-wave detectors

    D. Liang, Y. Gong, A. J. Weinstein, C. Zhang, and C. Zhang,Frequency response of space-based interferometric gravitational-wave detectors,Phys. Rev. D99(2019), no. 10 104027, [arXiv:1901.09624]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  23. [23]

    Frequency response of time-delay interferometry for space-based gravitational wave antennas

    C. Zhang, Q. Gao, Y. Gong, D. Liang, A. J. Weinstein, and C. Zhang,Frequency response of time-delay interferometry for space-based gravitational wave antenna,Phys. Rev. D100 (2019), no. 6 064033, [arXiv:1906.10901]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  24. [24]

    Sensitivity functions for space-borne gravitational wave detectors

    X.-Y. Lu, Y.-J. Tan, and C.-G. Shao,Sensitivity functions for space-borne gravitational wave detectors,Phys. Rev. D100(2019), no. 4 044042, [arXiv:2007.03400]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  25. [25]

    Full analytical formulas for frequency response of space-based gravitational wave detectors

    C. Zhang, Q. Gao, Y. Gong, B. Wang, A. J. Weinstein, and C. Zhang,Full analytical formulas for frequency response of space-based gravitational wave detectors,Phys. Rev. D 101(2020), no. 12 124027, [arXiv:2003.01441]

    work page internal anchor Pith review Pith/arXiv arXiv 2020

  26. [26]

    Galactic binary science with the new LISA design

    N. Cornish and T. Robson,Galactic binary science with the new LISA design,J. Phys. Conf. Ser.840(2017), no. 1 012024, [arXiv:1703.09858]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  27. [27]

    Reconstructing the evolution of white dwarf binaries: further evidence for an alternative algorithm for the outcome of the common-envelope phase in close binaries

    G. Nelemans and C. A. Tout,Reconstructing the evolution of white dwarf binaries: Further evidence for an alternative algorithm for the outcome of the common-envelope phase in close binaries,Mon. Not. Roy. Astron. Soc.356(2005) 753–764, [astro-ph/0410301]

    work page internal anchor Pith review Pith/arXiv arXiv 2005

  28. [28]

    Supernova Type Ia progenitors from merging double white dwarfs: Using a new population synthesis model

    S. Toonen, G. Nelemans, and S. Portegies Zwart,Supernova Type Ia progenitors from merging double white dwarfs: Using a new population synthesis model,Astron. Astrophys. 546(2012) A70, [arXiv:1208.6446]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  29. [29]

    The stochastic background: scaling laws and time to detection for pulsar timing arrays

    X. Siemens, J. Ellis, F. Jenet, and J. D. Romano,The stochastic background: scaling laws – 24 – and time to detection for pulsar timing arrays,Class. Quant. Grav.30(2013) 224015, [arXiv:1305.3196]. – 25 –

    work page internal anchor Pith review Pith/arXiv arXiv 2013