Pith. sign in

REVIEW 4 major objections 4 minor 18 references

Sub-horizon gravitational waves unavoidably seed large-scale white noise in cosmic curvature, and the quietness of the largest scales forbids nHz gravitational waves from sources earlier than redshift 10^8.

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 →

T0 review · deepseek-v4-flash

2026-08-01 09:11 UTC pith:ENM7KORN

load-bearing objection New and plausible GW-to-LSWN conversion; the headline constraint is real if the authors' unvalidated k_max bound holds, but the derivation itself is the value here. the 4 major comments →

arxiv 2607.27338 v1 pith:ENM7KORN submitted 2026-07-29 astro-ph.CO gr-qchep-ph

Gravitational Waves as a Source of Large-Scale White Noise: New Constraints

classification astro-ph.CO gr-qchep-ph PACS 04.30.-w98.80.-k
keywords gravitational waveslarge-scale white noisecurvature perturbationsradiation erapulsar timing arraysQCD phase transitioncosmological constraintsstochastic gravitational-wave background
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.

This paper tries to establish that a stochastic gravitational-wave background generated inside the horizon during the radiation era creates a measurable 'large-scale white noise' (LSWN) tail in the cosmic curvature field, and that the observed absence of such white noise yields a new, much stronger class of constraints on early-universe gravitational waves. The mechanism is universal: gravity waves shear the cosmic fluid, and the squared shear acts as a nonlinear source that accumulates curvature on the largest scales. The authors derive an amplitude k_BH^GW proportional to Ω_GW0*^2 (z_* H0)^4 / k_*^3 and combine it with the 99% confidence bound k_BH ≤ 1.80×10^-13 Mpc^-1 from their companion analysis of Planck 2018 data, obtaining z_*^2 Ω_GW0* < 5×10^7 (f_*/nHz)^{3/2}. A sympathetic reader would take away two things: pulsar-timing-array gravitational waves cannot come from before redshift roughly 10^8, and the QCD-scale gravitational-wave density today must be below about 5×10^-16, roughly seven orders of magnitude tighter than previous bounds. The paper stresses that this is a minimal constraint: acoustic waves and source granularity add more LSWN, so realistic models would be excluded even more strongly.

Core claim

At the paper's core is a conversion from gravitational waves to large-scale white noise. For sub-horizon tensor modes in the radiation era, the equal-time shear power spectrum is P_σ⊥ ≃ 6π^2 H^2 Ω_GW/k^3. Squaring the shear and integrating the unequal-time correlator over the radiation-era expansion history turns this into a curvature white-noise amplitude k_BH^GW ≃ 288π^3 Ω_GW0*^2 (z_* H0)^4 / k_*^3 (eq. 15). Requiring this not to exceed the observed total LSWN bound gives z_*^2 Ω_GW0* < 5×10^7 (f_*/nHz)^{3/2} — the paper's fundamental inequality. A covariant decomposition isolates a transverse part (gravitational-wave shear) from longitudinal parts (acoustic waves), showing the GW-only con

What carries the argument

The central object is k_BH, the large-scale white-noise amplitude in the comoving curvature perturbation, defined through ∆^2_R(k) ≈ ∆^2_iR(k) + k_BH/k; it is what the observations bound and what the theory predicts. The carrying mechanism is transverse shear: gravitational waves produce σ⊥ in the center-of-momentum frame, and the exact evolution equation for the covariant curvature K contains a θσ^2 source term, so squared shear integrates over time into curvature. The power spectrum of σ^2 is computed by convolving the GW shear power spectra, with a (1+μ^2)^2 kernel, and a double time integral with Hubble-weighted kernels turns it into white noise. The amplitude's Ω_GW^2 and (z_* H0)^4 sca

Load-bearing premise

The numerical constraints all rest on the companion claim that the total large-scale white-noise amplitude is k_BH ≤ 1.80×10^-13 Mpc^-1 at 99% confidence from Planck 2018 data; if that bound is wrong, too weak, or inapplicable to the gravitational-wave contribution, every inequality in this paper collapses.

What would settle it

Re-run the companion white-noise analysis on the public Planck 2018 temperature and polarization maps with independent foreground and covariance modeling; if the recovered 99% confidence upper limit on k_BH is more than a few times 1.80×10^-13 Mpc^-1, the QCD and redshift-horizon bounds weaken by the same factor and may stop excluding PTA waves from z_* ~ 10^8. A positive detection of curvature white noise above that amplitude in CMB lensing or large-scale-structure data would directly contradict the premise.

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

If this is right

  • Pulsar-timing-array gravitational waves (f_* near nHz, Ω_GW0 near 10^-8) cannot have been generated at redshift z_* ≳ 10^8; if they were, their shear would have left LSWN exceeding the Planck-based bound.
  • At the QCD phase-transition horizon scale, the present-day gravitational-wave density must satisfy Ω_GW0 ≲ 5×10^-16, about seven orders of magnitude tighter than direct-detection or BBN/CMB limits at that frequency.
  • For any proposed gravitational-wave observatory, eq. (20) defines a 'redshift horizon': the maximum source redshift at which its sensitivity could detect primordial waves without violating LSWN; under optimistic assumptions some proposed instruments could still reach early times, but foregrounds and extra LSWN shrink that horizon.
  • The constraints tighten, not loosen, in realistic models: acoustic waves from the same sources produce additional LSWN (eq. 40 suggests ϵ ≪ 1), so the early universe must have been dynamically quiet at small scales.
  • Strongly first-order phase transitions at high redshift with significant bulk kinetic energy are effectively ruled out as gravitational-wave sources because the associated flows would generate observable white noise even if the GWs themselves did not.

Where Pith is reading between the lines

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

  • If the mechanism is generic, LSWN becomes a universal back-reaction bound: any early-universe process with significant anisotropic stress or bulk flows — phase transitions, cosmic defects, primordial magnetic fields — is constrained by large-scale quietness, not just by direct GW searches.
  • A concrete testable extension is to compute ϵ for a specific phase-transition model; if ϵ is as small as the paper suggests, the real exclusion region could extend to lower frequencies and smaller Ω_GW0 than the minimal curve, potentially ruling out detection of TeV-scale transitions by proposed space-based interferometers.
  • An independent measurement of k_BH from CMB lensing or galaxy surveys would turn the inequality into a measurement: a positive white-noise tail would constrain the redshift distribution of the gravitational-wave background, information direct detectors cannot provide.
  • One could also search for the same white-noise signature in future CMB or 21-cm data, where the larger small-scale nonlinearities of the later universe might allow a positive detection of the shear-squared mechanism in action.

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

4 major / 4 minor

Summary. The paper derives the amplitude of large-scale white noise (LSWN) generated by a sub-horizon stochastic gravitational wave background in the radiation era, obtaining k_BH^GW ≃ 288π³ Ω_GW0*² (z_*H0)^4 / k_*³ (eq. 15), and combines it with the 99% CL LSWN bound k_max^BH = 1.80×10^-13 Mpc^-1 from the companion paper [2] to derive the constraint z_*² Ω_GW0* ≲ 5×10^7 (f_*/nHz)^{3/2}. This is used to argue that PTA-scale GWs cannot originate from z_* ≳ 10^8, and to set Ω_GW0 ≲ 5×10^-16 at the QCD scale. The paper also develops a covariant longitudinal/transverse decomposition of curvature and a discrete-event model for simultaneous GW and LSWN generation.

Significance. If correct, the paper provides a potentially powerful new constraint on high-redshift gravitational wave production, orders of magnitude stronger than existing bounds. Strengths include a self-contained derivation in outline from the covariant 1+3 formalism, an explicit minimal model with a conservative bound, and a useful decomposition separating GW-induced from matter-induced LSWN. The central claim is not circular: the input bound from [2] does not assume the GW result. However, the numerical constraints are conditional on the unpublished companion analysis [2], and the manuscript contains internal sign and numerical inconsistencies that must be fixed before the headline numbers can be taken at face value.

major comments (4)
  1. [Eqs. (12), (14), (15)] Substituting eq. (14) into eq. (12) gives k_BH^GW ≃ 288π³ Ω_GW0*² (z_*H0)^4 e^{+2τW} / k_*³, since the h_* factors cancel. Eq. (15) as printed has e^{-2τW}, the opposite sign. For τW ≈ 0.2 this is a 49% error in the amplitude and about a 22% error in the derived Ω_GW0 bounds. Please correct the sign and propagate the correction through eqs. (17)-(20).
  2. [Eqs. (17)-(18) and abstract] Using eq. (15) with z_* = 10^12, k_* = 1 pc^-1 ≃ 10^6 Mpc^-1 (or the nHz comoving wavenumber 6.5×10^5 Mpc^-1) and H0 ≃ 2.3×10^-4 Mpc^-1, I obtain k_BH^GW-QCD ∼ 10^4 (Ω_GW0/10^-8)^2 Mpc^-1, not the value 100 in eq. (17). The resulting 99% bound is Ω_GW0^QCD ≲ 5×10^-17 for f_* = 1 nHz, consistent with the abstract inequality but about an order of magnitude stronger than the 5×10^-16 quoted in eq. (18). Please reconcile the numerical estimates and use consistent values throughout.
  3. [Section 1, eq. (2) and all numerical results] The entire constraint chain scales with the 99% CL upper limit k_max^BH = 1.80×10^-13 Mpc^-1, taken from the authors' companion preprint [2]. The manuscript does not re-derive this bound, provide an error budget, or compare it with published LSWN limits. Since [2] is an unpublished Planck re-analysis, the headline constraints are conditional on its correctness and applicability. Please state this dependence explicitly, quantify the sensitivity (the derived Ω bounds scale as sqrt(k_max)), and ideally provide an independent check or a derivation of k_max.
  4. [Section 2, eq. (7)] The step from the equal-time result, eq. (6), to the unequal-time correlator, eq. (7), is asserted rather than derived. The text's justification that the phase average ⟨e^{ik(η-η′)}⟩ equals 1/2 is not correct as written; the relevant average appears to be over cos²(k Δη), and it must be shown that this factor survives the double time integral in eq. (8). Since this factor directly enters the amplitude in eq. (15), please provide a derivation or a controlled approximation.
minor comments (4)
  1. [Throughout] The word 'kurvature' is a recurring typo for 'curvature'; also 'arelic' in the Introduction should be 'a relic'.
  2. [Eq. (20)] The typesetting '1p6 Ω' is ambiguous; it should be rendered as 1/(√6 √Ω_GW0*) or similar.
  3. [Eq. (16)] The cutoff at k = 3√3 k_* is a factor of 5.2 in wavenumber, which is not 'only slightly' different from k_*. Please clarify the intended meaning.
  4. [Section 3.3 / Fig. 1] The exclusion plot would be easier to read if the LSWN bound curve included the τW damping factor or stated explicitly that damping is ignored.

Circularity Check

1 steps flagged

Headline GW constraints are a constant rescaling of the companion-paper bound k_max^BH from [2]; the GW-to-LSWN conversion is derived in this paper, but the numerical limits are not self-contained.

specific steps
  1. self citation load bearing [Section 1, Eq. (2); Section 3.3, Eqs. (17)-(19); abstract inequality]
    "The observational constraint we use is the 99% confidence level limit on LSWN [2] k_BH ≤ k_max^BH = 1.80×10−13 Mpc−1 (2) based on Planck 2018 data ([5]). ... which must not exceed the observational constraint of eq. (2). This leads to the constraint Ω_QCD_GW0 <5×10−16 and Ω_QCD_GW* <8×10−12 (18)"

    Every numerical bound in the Letter (Eqs. 17-20 and the abstract inequality) is obtained by algebraically inverting the derived relation k_GW_BH ≃ 288π³ Ω_GW0*² (z* H0)⁴ / k*³ against this single k_max value. The limit k_max^BH = 1.80×10−13 Mpc−1 is not derived or cross-checked here; it is imported from the same-authors companion paper [2]. Thus the reported Ω_GW0 constraints and redshift horizons are a constant rescaling of [2]'s input bound. The GW-to-LSWN conversion is independent, but the numerical content of the headline result is load-bearing on a self-citation that is not independently validated within this manuscript.

full rationale

The paper's central derivation is not logically circular: Eq. (15) follows from the shear power spectrum and the shear-driven kurvature equation, and it does not assume the non-detection of LSWN. The observational input k_max^BH from [2] is an external data-based bound (Planck 2018), not the GW result itself. However, the Letter is not numerically self-contained: all quoted constraints scale directly with that single 99% CL number, which comes from the authors' own companion preprint, and no independent re-analysis or error budget is provided. This is a load-bearing self-citation and a validation risk, but the new physical content — the translation from GW shear to white-noise curvature — is derived in the present paper. The paper also repeatedly labels its bounds 'minimal' and 'conservative', which is a caveat on strength, not a circularity. An independent CMB re-analysis of k_max would resolve the concern; hence the moderate score of 4 rather than a higher constructional-circularity score.

Axiom & Free-Parameter Ledger

2 free parameters · 7 axioms · 0 invented entities

The constraint calculation rests on the authors' prior LSWN framework for both the physical mechanism and the key observational bound, plus several explicit approximations (Gaussian random phases, factor 1/2, instantaneous production). The only fitted number introduced externally is k_max^BH; the GW density and redshift are the targets being constrained, not fit parameters.

free parameters (2)
  • k_max^BH (observational LSWN bound) = 1.80e-13 Mpc^-1
    99% CL upper limit on total LSWN amplitude fitted from Planck 2018 in the authors' companion paper [2]; every numerical constraint in this paper scales with this value.
  • phase-averaging factor for unequal-time correlator = 1/2
    Introduced by hand in eq. (7) to represent time averaging of rapidly oscillating GW modes; it is load-bearing for the time-integrated LSWN amplitude but not derived from the spectrum.
axioms (7)
  • domain assumption The exact evolution equation for kurvature K, eq. (23), and its solution eq. (3), in which shear squared sources large-scale white noise.
    Taken from Papers I/II and covariant GR; the kurvature/LSWN framework is the authors' own construction, not independently established.
  • domain assumption Gravitational waves in the cosmic fluid generate shear with P_σ⊥ = 6π²H²Ω_GW/k³ for sub-horizon radiation-era modes.
    Standard linear tensor-mode relation plus the GW energy density definition; central input to Step 1.
  • ad hoc to paper The unequal-time correlator of σ² factorizes as H(η)²H(η')²Ω(η)Ω(η') with a constant phase factor 1/2.
    Eq. (7) asserts the 1/2 factor without a complete derivation; this choice affects the double time integral in Step 3 and hence the final amplitude.
  • ad hoc to paper The minimal model assumes instantaneous, Gaussian, statistically homogeneous GW production at z*, with a truncated white-noise spectrum.
    Defines the fiducial scenario and the characteristic k* in eq. (13); the authors state realistic sources would give more LSWN, so this is conservative.
  • domain assumption Vorticity of the cosmic fluid is negligible, justifying the longitudinal/transverse decomposition.
    Stated in §4.2 and §A.2; if vorticity is not small, the clean separation of GW-induced kurvature breaks down.
  • domain assumption The observational non-detection of LSWN (eq. 2) is a valid bound on the total white-noise curvature.
    External input from the authors' previous paper [2]; not independently verified in this manuscript.
  • standard math Radiation-dominated background H(z)=H0√Ωr0(1+z)² and standard degrees-of-freedom scaling h(z).
    Used throughout to evaluate redshift integrals and to set the horizon wavenumber.

pith-pipeline@v1.3.0-daily-deepseek · 15776 in / 32695 out tokens · 296556 ms · 2026-08-01T09:11:40.907186+00:00 · methodology

0 comments
read the original abstract

A stochastic gravitational wave (GW) background sources a shear in the flow of cosmic fluid which, through non-linear mode coupling, generates large-scale white noise (LSWN) in the kurvature density field. Building on the LSWN framework of our previous work, we derive the amplitude of this GW-induced LSWN and translate the observational non-detection of LSWN into bounds on the production redshift and density of gravity waves. In particular, a minimal constraint on gravity waves with $z=0$ density parameter $\Omega_\mathrm{GW0}^*$ in frequency band $f_*$ generated at redshift $z_*$ must satisfy ${z_*}^2\,\Omega_\mathrm{GW0}^*<5\times10^7\,(f_*/\mathrm{nHz})^{3/2}$. This, for example, precludes the gravity waves recently detected by pulsar timing arrays \cite{NANOGrav:2023hvm} from being present before $z_*\sim10^8$, long after the quark hadron phase transition. While orders of magnitude stronger than other constraints on gravity waves, this is a minimal LSWN constraint as realistic modeling of early universe gravity wave production, including the granularity of the gravity-wave sources and the LSWN produced by the associated acoustic waves would probably tighten this constraint by orders of magnitude.

Figures

Figures reproduced from arXiv: 2607.27338 by Albert Stebbins, Gabriela Barenboim.

Figure 1
Figure 1. Figure 1: Exclusion plot in the (k⋆, ΩGW) plane of gravity waves produced at the horizon scale. The solid blue curve shows the LSWN bound from this work, Eq. (19), with the blue shaded region excluded. The gray dashed line and shading show the BBN/CMB/Neff integrated ceiling [10]. The orange dashed segment shows the direct GW search bound over its relevant k band [11]. The purple dashed vertical line marks the QCD h… view at source ↗
Figure 2
Figure 2. Figure 2: Plotted is the redshift horizon of various surveys. This is the maximum redshift from which [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

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

18 extracted references · 12 linked inside Pith

  1. [1]

    Large Scale White Noise and Cos- mology

    Gabriela Barenboim, Aurora Ireland, and Albert Stebbins. “Large Scale White Noise and Cos- mology”. In: (Nov. 2025). arXiv:2511.13866 [astro-ph.CO]

  2. [2]

    The Noisy Universe

    Gabriela Barenboim, Aurora Ireland, and Albert Stebbins. “The Noisy Universe”. In: (Nov. 2025). arXiv:2511.15803 [astro-ph.CO]

  3. [3]

    The NANOGrav 15-Year Data Set: Evidence for a Gravitational- Wave Background

    NANOGrav Collaboration. “The NANOGrav 15-Year Data Set: Evidence for a Gravitational- Wave Background”. In:Astrophys. J. Lett.951.1 (2023), p. L8.doi:10 . 3847 / 2041 - 8213 / acdac6. arXiv:2306.16213 [astro-ph.HE]

  4. [4]

    Daniel Baumann.TASI Lectures on Inflation. 2012. arXiv:0907.5424 [hep-th].url:https: //arxiv.org/abs/0907.5424

  5. [5]

    Planck 2018 results. VI. Cosmological parameters

    N. Aghanim et al. “Planck 2018 results. VI. Cosmological parameters”. In:Astron. Astrophys.641 (2020). [Erratum: Astron.Astrophys. 652, C4 (2021)], A6.doi:10.1051/0004-6361/201833910. arXiv:1807.06209 [astro-ph.CO]

  6. [6]

    Damping of Tensor Modes in Cosmology

    Steven Weinberg. “Damping of Tensor Modes in Cosmology”. In:Phys. Rev. D69 (2004), p. 023503.doi:10.1103/PhysRevD.69.023503. arXiv:astro-ph/0306304 [astro-ph]

  7. [7]

    The Atacama Cosmology Telescope: DR6 Constraints on Extended Cosmological Models

    Erminia Calabrese et al. “The Atacama Cosmology Telescope: DR6 Constraints on Extended Cosmological Models”. In:Journal of Cosmology and Astroparticle Physics2025.11 (2025), p. 063.doi:10.1088/1475-7516/2025/11/063. arXiv:2503.14454 [astro-ph.CO]

  8. [8]

    Cosmological Backgrounds of Gravitational Waves

    Chiara Caprini and Daniel G. Figueroa. “Cosmological Backgrounds of Gravitational Waves”. In:Class. Quant. Grav.35.16 (2018), p. 163001.doi:10 . 1088 / 1361 - 6382 / aac608. arXiv: 1801.04268 [astro-ph.CO]

  9. [9]

    Probing Physics Beyond the Standard Model: Limits from BBN and the CMB Independently and Combined

    Tsung-Han Yeh et al. “Probing Physics Beyond the Standard Model: Limits from BBN and the CMB Independently and Combined”. In: (2022). arXiv:2207.13133 [astro-ph.CO]

  10. [10]

    Constraints on primordial grav- itational waves from the Cosmic Microwave Background

    Thomas J. Clarke, Edmund J. Copeland, and Adam Moss. “Constraints on primordial grav- itational waves from the Cosmic Microwave Background”. In:JCAP10 (2020), p. 002.doi: 10.1088/1475-7516/2020/10/002. arXiv:2004.11396 [astro-ph.CO]

  11. [11]

    Upper Limits on the Isotropic Gravitational-Wave Background from the first part of LIGO, Virgo, and KAGRA’s fourth Observing Run

    R. Abbott et al. “Upper Limits on the Isotropic Gravitational-Wave Background from the first part of LIGO, Virgo, and KAGRA’s fourth Observing Run”. In: (2025). arXiv:2508 . 20721 [gr-qc]

  12. [12]

    Measuring the spectrum of primordial gravitational waves with CMB, PTA and laser interferometers

    Paolo Campeti et al. “Measuring the spectrum of primordial gravitational waves with CMB, PTA and laser interferometers”. In:Journal of Cosmology and Astroparticle Physics2021.1, 012 (Jan. 2021), p. 012.doi:10.1088/1475-7516/2021/01/012. arXiv:2007.04241 [astro-ph.CO]

  13. [13]

    A Space-Time Fluid (Unabridged)

    Albert Stebbins. “A Space-Time Fluid (Unabridged)”. In: (Jan. 2026). arXiv:2601 . 16996 [physics.gen-ph]

  14. [14]

    Relativistic cosmology and large- scale structure

    Christos G. Tsagas, Anthony Challinor, and Roy Maartens. “Relativistic cosmology and large- scale structure”. In:Phys. Rept.465 (2008), pp. 61–147.doi:10.1016/j.physrep.2008.03.003. arXiv:0705.4397 [astro-ph]

  15. [15]

    George F. R. Ellis, Roy Maartens, and Malcolm A. H. MacCallum.Relativistic Cosmology. Cam- bridge: Cambridge University Press, 2012.isbn: 9780521381154

  16. [16]

    Theodore Frankel.The Geometry of Physics: An Introduction. 3rd. Cambridge: Cambridge Uni- versity Press, 2011.doi:10.1017/CBO9781139061377

  17. [17]

    Conformally invariant orthogonal decomposition of symmetric tensors on Rie- mannian manifolds and the initial-value problem of general relativity

    James W. York. “Conformally invariant orthogonal decomposition of symmetric tensors on Rie- mannian manifolds and the initial-value problem of general relativity”. In:J. Math. Phys.14 (1973), pp. 456–464.doi:10.1063/1.1666338

  18. [18]

    Eric Gourgoulhon.3+1 Formalism and Bases of Numerical Relativity. 2007. arXiv:gr - qc / 0703035 [gr-qc]. 16