REVIEW 3 major objections 4 minor 26 references
The Fate of Large Scale White Noise in Second-Order Cosmological Perturbation Theory
T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Second-order cosmology kills the claimed infrared enhancement of the curvature power spectrum.
desk verdict A clear, honest self-correction that fixes the second-order dictionary and kills the claimed IR pole in R2, with the caveat that the vanishing of the pole coefficient is conditional on assuming no primordial 1/k^2 second-order mode and the proof is limited to a perfect fluid in radiation domination. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing identity is the exact second-order relation $\Delta\rho_2 = \frac{1}{4\pi G a^2}\!\left(\nabla^2 R_2 + Q - \nabla^2 C\right)$, which replaces the linear Poisson relation $\nabla^2 R = 4\pi G a^2 \Delta\rho$. Here $Q$ is a quadratic, locally constructed source built from shear and vorticity, and $C$ is the quadratic completion needed to make $R_2$ gauge invariant; in the soft limit $k\to0$, $\nabla^2 C$ vanishes while $Q$ stays finite, so the white noise is carried entirely by $Q$. The second load-bearing object is the conserved pole coefficient $F = P'/H + 3P + \tilde{\Xi}_0$, defined from the $1/k^2$ part $P$ of the second-order metric perturbation and the soft part of the source $\tilde{\Xi}_0$. The paper proves $F' = 0$ identically from the second-order equation of motion, so $F$ can only be set by initial data; requiring no independent second-order super-horizon mode sets $F=0$.
What would settle it
A direct check would be to compute $F$ or $\lim_{k\to0}(k^2 R_2)$ in an explicit second-order radiation-domination setup that includes a nonvanishing homogeneous solution $P_{\rm hom}=G+D/\eta^3$; if a nonzero constant mode $G$ can be physically generated rather than inserted by hand, the conclusion $F=0$ fails. A second, complementary check is to repeat the soft-limit analysis for an imperfect fluid or during matter domination: if the source $\tilde{\Xi}$ develops an $O(\epsilon^{-2})$ component or the conservation law $F'=0$ is broken, then the infrared protection of $R_2$ does not hold in those regimes.
Extended reading notes
Core claim
The central discovery is that the white noise carried by the second-order kurvature density $\Delta\rho_2$ resides entirely in a quadratic extrinsic-curvature source $Q$, and that this source is not inherited by the comoving curvature perturbation. Working in Poisson gauge and specializing to a perfect fluid in radiation domination, the paper derives the exact second-order relation $\Delta\rho_2 = \frac{1}{4\pi G a^2}(\nabla^2 R_2 + Q - \nabla^2 C)$, where $C$ is the quadratic completion that converts the metric combination $\psi_2 - H v_2$ into the gauge-invariant $R_2$. The earlier linear Poisson relation $\nabla^2 R = 4\pi G a^2 \Delta\rho$ omitted the term $Q - \nabla^2 C$; since $\nabla^2 C\to 0$ in the soft limit while $Q$ tends to a constant, the omitted term is not a small correction but the entire white-noise content. Solving the second-order Einstein equations, the paper proves that the coefficient $F$ of the potential $O(\epsilon^{-2})$ pole in $R_2$ is an exact constant of motion during radiation domination and, with no primordial second-order homogeneous mode, $F=0$. Hence $\lim_{k\to0}(k^2 R_2)=0$, the claimed infrared enhancement is an artifact of mixing first- and second-order equations, and the $k_{\rm BH}$ bound derived from it does not follow.
Load-bearing premise
The proof assumes there is no pre-existing second-order curvature mode with a $1/k^2$ component on super-horizon scales—the homogeneous soft-mode coefficients $G$ and $D$ are set to zero—and it is restricted to a perfect fluid in radiation domination, so the protection of $R_2$ is not established for imperfect fluids or other epochs.
Editorial extensions
If this is right
- The $k_{\rm BH}$ bound proposed in Ref. [2], and its inferred constraint on the small-scale primordial power spectrum, do not apply.
- The second-order kurvature density $\Delta\rho_2$ still has a genuine, large-scale white-noise spectrum, now computed explicitly during radiation domination including the exact weight function.
- The comoving curvature perturbation $R_2$ is protected: $\lim_{k\to0}(k^2 R_2)=0$, so no $1/k$ contribution to the dimensionless curvature power spectrum arises from this mechanism.
- Because $C$ is $O(\epsilon^0)$, $R_2$ itself acquires a finite, infrared-convergent white-noise contribution at second order, distinct from the claimed enhancement.
- The proof holds for a perfect fluid in radiation domination; extending it to imperfect fluids, neutrino free streaming, or other eras is left as an open problem.
- In the soft limit, the conserved pole coefficient $F$ may be related to a Langlois-Vernizzi-type conserved charge, a connection the paper flags but does not establish.
Reading between the lines
- Going beyond the paper: the finite white-noise piece in $R_2$ from the completion $C$ could still leave a subdominant, scale-independent large-scale signature in the curvature power spectrum, and characterizing its amplitude would give a concrete target for future lattice or numerical Boltzmann calculations.
- Going beyond the paper: because the perfect-fluid radiation source $Q$ grows with internal wavenumber for a scale-invariant primordial spectrum, the same UV sensitivity that motivated the $k_{\rm BH}$ bound is still present in $\Delta\rho$, so realistic imperfect-fluid effects such as viscosity or Silk damping may act as the physical cutoff.
- Going beyond the paper: the no-pole result depends on the absence of a primordial second-order $1/k^2$ curvature mode; a testable extension would be to check whether non-Gaussian initial conditions or vector/tensor sources can generate a nonzero $F$ at later times.
- Going beyond the paper: the mismatch $Q \neq \nabla^2 C$ is the precise sense in which the linear Poisson relation fails as a second-order dictionary, and the same failure could affect other attempts to convert covariant nonlinear curvature measures into statements about gauge-invariant perturbation variables.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper revisits the Large Scale White Noise (LSWN) proposal of Refs. [1,2], working to second order in cosmological perturbation theory in Poisson gauge. It confirms that the covariant kurvature density Δρ does acquire a white-noise power spectrum on large scales, and it explicitly computes the quadratic source Q for a perfect radiation-dominated fluid. The central new claim is that this white noise is not inherited by the second-order comoving curvature perturbation: lim_{k→0} k²R₂ = 0. The argument proceeds by deriving a conserved coefficient F for the would-be 1/k² pole in ψ₂−Hv₂, showing that, under the chosen homogeneous-mode initial conditions G=D=0, this coefficient is zero. The paper concludes that the k_BH bound on the small-scale primordial power spectrum proposed in Ref. [2] is invalid.
Significance. If the central claim were established without qualification, the paper would resolve an important inconsistency in the recent LSWN literature: the white noise in Δρ is real, but it does not translate into an infrared enhancement of the curvature power spectrum. The paper's positive contributions include a clean second-order derivation of the relation Δρ₂ = (4πGa²)⁻¹(∇²R₂ + Q − ∇²C), an explicit computation of the weight function W(η,q), and the observation that the soft-limit pole coefficient obeys a conservation law. These are genuine analytic results and are presented transparently. The main weakness is that the headline protection claim is proved only under the initial-condition assumption G=D=0 and only for a perfect fluid in radiation domination, yet the abstract and introduction state it unconditionally. The structural derivation is not circular, but the absence of a primordial second-order superhorizon mode is an extra physical assumption that needs to be stated and justified.
major comments (3)
- [§IV.C, Eqs. (64)–(70)] The proof that F=0, and hence that lim_{k→0} k²(ψ̃₂−Hṽ₂)=0, does not follow from the second-order Einstein equations alone. Equation (66) establishes only that F is time-independent; its value is fixed by the initial data through the homogeneous solution P_hom = G + D/η³. The choice G=D=0 is an additional physical assumption, namely that there is no primordial second-order superhorizon curvature mode. The manuscript itself acknowledges that G≠0 would give F≠0. Since Eq. (74) and the central claim depend on this choice, the unqualified statement in the abstract and Introduction that R₂ is protected is not supported by the derivation as presented. The claim must be made conditional on the vanishing of such a mode, or a physical argument for G=0 must be supplied.
- [§V and §IV.C] The proof is restricted to a perfect fluid in radiation domination, and the early-time limit η→0 is used to discard the homogeneous mode. In a realistic cosmology the radiation era begins at a finite reheating time η_i>0, so a constant second-order soft mode present at η_i is not excluded by regularity at η→0; matching from an inflationary phase is not discussed. Section V explicitly limits the proof to a perfect fluid in radiation domination, but the abstract and Introduction state the protection of R₂ without these qualifications. The scope of the central result should be stated consistently throughout, or the generality of the conclusion should be established.
- [Abstract and Introduction] The abstract says that “R₂ is protected from developing any such IR enhancement” and that the constraint proposed in Ref. [2] is therefore invalid. Given Major Comments 1 and 2, the second statement is stronger than what is proved: the refutation of the k_BH bound is conditional on the absence of a primordial second-order homogeneous mode and on the perfect-fluid/radiation-domination setting. While the paper later notes these limitations, the headline claims should carry the same qualifications; otherwise the reader is presented with an overstatement of the result.
minor comments (4)
- [Eq. (44) and preceding text] After Eq. (36) the “prim” superscript on R₁ is dropped, so it should be made explicit that P_R1(q) in Eq. (44) denotes the primordial power spectrum, not the full time-dependent one.
- [Fig. 1] The left panel caption says “Behavior of the kernel (48),” but Eq. (48) defines k_BH as an integral and the weight function W; please clarify that the plotted object is W(η,q), not the integrand f(η,k,q) of Eq. (43).
- [§III, after Eq. (47)] The sentence “the apparent IR divergence in Eq. (47) is an artifact of using the first-order dictionary” is accurate, but it would be helpful to restate explicitly that Eq. (47) also assumes P_R1 is the first-order dimensionless spectrum and that the second-order contribution R₂ is being identified with the full R; this is clear from context but deserves a one-sentence reminder.
- [§V, future directions] The remark that R₂ itself develops a white-noise contribution from C is interesting, but it is presented in a footnote-like aside; consider placing it in the main text with a precise statement of the expected k→0 limit of the R₂ power spectrum, since this is a natural follow-up for readers.
Circularity Check
No significant circularity; the derivations are self-contained, and the main limitation is an openly stated initial-condition assumption rather than a circular reduction.
full rationale
The paper's derivation chain is self-contained. It obtains the exact second-order relation Δρ2 = (1/4πGa^2)(∇²R2 + Q − ∇²C) by expanding the covariant kurvature density in metric variables and using the second-order Einstein equations; it then computes Q explicitly during radiation domination (Eqs. 38–44), derives the soft-limit conservation law F' = 0 (Eq. 66) from the equation of motion for the pole coefficient, and checks that C is O(ε^0) so that ∇²C vanishes in the soft limit (Eq. 73). The final result lim_{k→0} k²R2 = 0 follows from Eq. (65) together with the conservation law and the explicitly stated initial condition G = D = 0 in the homogeneous solution; the paper itself notes that G ≠ 0 would give F ≠ 0. This is a transparent assumption that no independent primordial second-order homogeneous mode is present, not a concealed reuse of the target result: the conservation law and the Q ≠ ∇²C identification are independent of the value of F and are not imported from the author's prior papers. The self-citations are to the work being critiqued or to the definition of Δρ, whose relevant properties are re-derived here. The abstract's unconditional phrasing is stronger than the conditional proof, and Section V restricts the proof to a perfect fluid in radiation domination, but those are scope and correctness caveats, not circularity.
Assumptions & free parameters
free parameters (2)
- G (constant homogeneous soft mode) =
0 (assumed)
- D (decaying homogeneous soft mode) =
0 (regularity)
assumptions (6)
- domain assumption Standard second-order cosmological perturbation theory and gauge transformations in the Malik-Wands formalism
- domain assumption Perfect fluid with vanishing anisotropic stress (π_μν=0)
- domain assumption Radiation domination with w=c_s^2=1/3 and a∝η
- domain assumption Gaussian initial conditions for the primordial curvature perturbation
- ad hoc to paper No primordial second-order homogeneous mode (G=D=0)
- domain assumption Scalar perturbations only, with vorticity contribution neglected
Cite this review
Pith. "Pith review of The Fate of Large Scale White Noise in Second-Order Cosmological Perturbation Theory." pith.science (2026). https://pith.science/paper/H4GLB4A6
@misc{pith2026260809897,
author = {Pith},
title = {Pith review of: The Fate of Large Scale White Noise in Second-Order Cosmological Perturbation Theory},
year = {2026},
howpublished = {\url{https://pith.science/paper/H4GLB4A6}},
note = {Machine review of arXiv:2608.09897}
}
abstract
Working to second order in cosmological perturbation theory, we reconsider the Large Scale White Noise (LSWN) effect proposed in [2511.13866] and [2511.15803]. We demonstrate that the kurvature density variable $\Delta \rho$ does indeed develop LSWN at second order. However, contrary to the claims of [2511.15803], we show that this is not inherited as an infrared (IR) pole in the second-order comoving curvature perturbation $R_2$. This apparent enhancement was an artifact of using a linear Poisson equation to relate two genuinely second order quantities. Further, we show that $R_2$ is protected from developing any such IR enhancement: $\lim_{k \rightarrow 0} \left( k^2 R_2 \right) = 0$, which follows from a conservation law in the soft limit. The constraint proposed in [2511.15803] and its implications for the small-scale primordial power spectrum are therefore invalid.
Figures
Reference graph
Works this paper leans on
-
[2]
+∇ 2(ϕ2−ψ 2) + 3 2a′′ a −H 2 ϕ2 +Q (3) = 12πGa2c2 eff,2δρ2,(49c) δρ′ 2 + 3H(1 +c 2 eff,2)δρ2 + (1 +w)ρ 0(∇2v2−3ψ′
-
[1]
(B6) For future convenience, note that the trace of the latter reads Ci i =−3ψ+∇ 2E
Metric Perturbations We write the covariant metric as gµν =g (0) µν +δg (1) µν + 1 2δg(2) µν , (B3) with components g(0) 00 =−a 2, δg (1) 00 =−2a 2ϕ1, δg (2) 00 =−2a 2ϕ2,(B4a) g(0) 0i = 0, δg (1) 0i =a 2B1i, δg (2) 0i =a 2B2i,(B4b) g(0) ij =a 2δij, δg (1) ij = 2a2C1ij, δg (2) ij = 2a2C2ij,(B4c) where at each order the 0iandijcomponents can be decomposed i...
-
[3]
+Q (4) = 0, (49d) (1 +w)ρ 0 ρ′ 0 ρ0 ∂iv2 +∂iv′ 2 +∂iϕ2 + 4H∂iv2 +c 2 eff,2∂iδρ2 +Q (5) i = 0,(49e) wherec 2 eff,2≡δP 2/δρ2, and the quadratic first-order sourcesQ (1),Q (2) i ,Q (3),Q (4),Q (5) i are collected in Appendix E. Note that the Bianchi identity relates the 00, 0i,ijtrace, and energy conservation equations, and so of these 5 equations, only 4 ar...
-
[4]
(61) 11 Now, recall from Eq. (58) that the combination that matters isψ 2−Hv 2. From Eq. (49b), we have ˜ψ2−H˜v2 = 3 2 ˜ψ2 + ˜ψ′ 2 2H + ˜q(2) 2H + ˜Ξ 2k2 , (62) where we have also used Eq. (54) to substitute ˜ϕ2 = ˜ψ2 + ˜Ξ/k2. As already established,q (2) is pole-free, so theO(ϵ−2) piece vanishes provided P′ + 3HP=−H ˜Ξ0. (63) Notice that this combination...
-
[5]
Connection Coefficients The connection coefficients can be expanded as Γα βγ = (0)Γα βγ + (1)δΓα βγ + 1 2 (2)δΓα βγ, (B10) with the following entries: 000: (0)Γ0 00 =H, (B11a) (1)δΓ0 00 =ϕ′ 1, (B11b) (2)δΓ0 00 =ϕ′ 2−4ϕ 1ϕ′ 1 + 2Bi 1 B′ 1i +HB 1i +∂iϕ1 , (B11c) 00i: (0)Γ0 0i = 0, (B12a) 16 (1)δΓ0 0i =∂ iϕ1 +HB 1i, (B12b) (2)δΓ0 0i =∂ iϕ2 +HB 2i−4ϕ 1∂iϕ1−4H...
-
[6]
Matter Perturbations We define the fluid 4-velocityu µ to be the tangent vector field to the worldlines of fluid elements. Lettingx µ(τ) be a fluid element’s worldline, withτthe proper time comoving with the fluid, one has uµ = dxµ dτ . (B17) From this definition along with the definition of proper timedτ 2 =−g µνdxµdxν, it follows that the fluid 4-veloci...
-
[7]
Geometric Quantities Here we derive and list the components of the spatial projection tensor, expansion, shear, vorticity, and acceleration corresponding to the matter 4-velocityu µ. Spatial Projection T ensor: Expanding Eq. (A1) as Pµν =P (0) µν +δP (1) µν + 1 2δP(2) µν , (B27) and using Eqs. (B4) and (B25), the components of the spatial projector at eac...
-
[8]
Background At the background level, the momentum constraint is trivially satisfied and the energy constraint is simply the first Friedmann equation, H2 = 8πG 3 a2ρ0, (C3) whereH=aHis the conformal Hubble rate. Similarly, there is only one independent evolution equation — the second Friedmann equation, H′ =− 4πG 3 a2 ρ0 + 3p0 , (C4) where ′ denotes a deriv...
Show all 26 references
-
[9]
For the scalar perturbations, the energy and momentum constraints are 3H ψ′ 1 +Hϕ 1 −∇ 2ψ1−H∇ 2 E′ 1−B 1 =−4πGa 2δρ1, (C6) ψ′ 1 +Hϕ 1 =−4πGa 2 ρ0 +p 0 v1 +B 1
First Order At first order, scalar, vector, and tensor perturbations decouple, and so one can derive the constraint and dynamical equations for the various sectors independently. For the scalar perturbations, the energy and momentum constraints are 3H ψ′ 1 +Hϕ 1 −∇ 2ψ1−H∇ 2 E′...
-
[10]
Second Order Here we collect energy-momentum conservation and Einstein equations at second order. From the former, we have the energy conservation equation δρ′ 2 + 3H(δρ 2 +δP 2) + (1 +w)ρ0 Ci 2i ′ +∂ iv2i + 2vi 1∂i (δρ1 +δP 1) + 2 (δρ1 +δP 1) Ci 1i ′ +∂ iv1i + 2(1 +w)ρ 0 2vi ...
-
[11]
(v1i +B 1i) + 2 (δρ1 +δP 1) (v′ 1i +B′ 1i) + 2 (δρ1 +δP 1) (∂iϕ1 + 4H(v1i +B 1i)) −2(1 +w)ρ ′ 0 (v1i + 2B1i)ϕ 1−2C 1ijvj 1 + 2(1 +w)ρ 0 (v1i +B 1i) Cj 1j ′ +∂jvj 1 −B 1i (ϕ′ 1 + 8Hϕ1) + 2C′ 1ijvj 1 + 2Cj 1iv′ 1j + 2(1 +w)ρ 0 vj 1 (∂j(v1i +B 1i)−∂ iB1j + 8HC1ij)−ϕ 1 (v′ 1i + 2B...
-
[12]
(C17) For the Einstein equations, Eq
+ 8H(∂iv1)(∂iv1) + 3ψ′ 1ψ1 +ψ 1∇2v1−(∂ iv1)(∂iψ1) = 0, (C16) and (C13) simplifies to (1 +w)ρ 0 ρ′ 0 ρ0 ∂iv2 +∂iv′ 2 +∂iϕ2 + 4H∂iv2 +∂iδP2 + 2 (δρ′ 1 +δP ′ 1)∂iv1 + 2 (δρ1 +δP 1) (∂iv′ 1 +∂iϕ1 + 4H∂iv1)−6(1 +w)ρ ′ 0ψ1∂iv1 + 2(1 +w)ρ 0 (∂iv1)∇2v1 + (∂jv1)∂i∂jv1−5ψ′ 1∂iv1−3ψ 1∂iv...
-
[13]
+∇ 2(ϕ2−ψ 2) + 3 2a′′ a −H 2 ϕ2 −6(∂ iψ1)(∂iψ1)−24Hψ 1ψ′ 1−8ψ 1∇2ψ1−3ψ′ 1ψ′ 1 + 12ψ2 1 H2−2 a′′ a = 4πGa2 3δP2 + 2ρ0(1 +w)(∂iv1)(∂iv1) . (C20) Appendix D: Second Order Comoving Curvature Perturbation Here we review the construction of the gauge invariant comoving curvature per...
-
[14]
(E1e) 25
+ 8H(∂iv1)(∂iv1) + 3ψ′ 1ψ1 +ψ 1∇2v1−(∂ iv1)(∂iψ1) , (E1d) Q(5) i =2(1 +c 2 s)δρ1 δρ′ 1 δρ1 ∂iv1 +∂iv′ 1 +∂iϕ1 + 4H∂iv1 −6(1 +w)ρ ′ 0ψ1∂iv1 + 2(1 +w)ρ 0 (∂iv1)∇2v1 + (∂jv1)∂i∂jv1−5ψ′ 1∂iv1−3ψ 1∂iv′ 1−2ψ 1∂iψ1−12Hψ 1∂iv1 . (E1e) 25
-
[15]
Barenboim, A
G. Barenboim, A. Ireland, and A. Stebbins, Large Scale White Noise and Cosmology, (2025), arXiv:2511.13866 [astro- ph.CO]
2025 arXiv
-
[16]
Barenboim, A
G. Barenboim, A. Ireland, and A. Stebbins, The Noisy Universe, (2025), arXiv:2511.15803 [astro-ph.CO]
2025 arXiv
-
[17]
P. J. E. Peebles, The Gravitational-Instability Picture and the Nature of the Distribution of Galaxies, apjl189, L51 (1974)
1974
-
[18]
P. J. E. Peebles,The large-scale structure of the universe(1980)
1980
-
[19]
K. A. Malik and D. Wands, Cosmological perturbations, Phys. Rept.475, 1 (2009), arXiv:0809.4944 [astro-ph]
2009 arXiv
-
[20]
A. J. Christopherson,Applications of Cosmological Perturbation Theory, Ph.D. thesis, Queen Mary, U. of London, Math. Sci. (2011), arXiv:1106.0446 [astro-ph.CO]
2011 arXiv
-
[21]
D. H. Lyth, Large Scale Energy Density Perturbations and Inflation, Phys. Rev. D31, 1792 (1985)
1985
-
[22]
Bartolo, S
N. Bartolo, S. Matarrese, and A. Riotto, CMB Anisotropies at Second Order I, JCAP06, 024, arXiv:astro-ph/0604416
-
[23]
Bartolo, S
N. Bartolo, S. Matarrese, and A. Riotto, CMB Anisotropies at Second-Order. 2. Analytical Approach, JCAP01, 019, arXiv:astro-ph/0610110
-
[24]
Langlois and F
D. Langlois and F. Vernizzi, Evolution of non-linear cosmological perturbations, Phys. Rev. Lett.95, 091303 (2005), arXiv:astro-ph/0503416
2005 arXiv
-
[25]
Langlois and F
D. Langlois and F. Vernizzi, Conserved non-linear quantities in cosmology, Phys. Rev. D72, 103501 (2005), arXiv:astro- ph/0509078
2005
-
[26]
R. M. Wald,General Relativity(Chicago Univ. Pr., Chicago, USA, 1984)
1984
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.