Pith. sign in

REVIEW 2 major objections 4 minor 1 cited by

Stability of fluids in spacetimes with decelerated expansion

T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read The paper claims that homogeneous barotropic fluids in decelerating FLRW spacetimes are asymptotically stable under the explicit condition $K<1-2/(3\alpha)$, with shock formation for dust at $\alpha\le 1/2$ and for radiation at $\alpha\le…

desk verdict A credible new stability theorem for barotropic fluids on decelerating FLRW, with an explicit lowest-order proof and a load-bearing higher-order section that needs to be filled in before the result is fully established. read the letter →

arxiv 2501.12798 v1 pith:SFYUPCNK submitted 2025-01-22 gr-qc math-phmath.MP

David Fajman , Maximilian Ofner , Todd Oliynyk , Zoe Wyatt This is my paper · ORCID
classification gr-qcmath-phmath.MP MSC 35Q7583C0535B35
keywords relativisticEulerequationsFLRWspacetimesnonlinearstabilitydeceleratedexpansionbarotropicperfectfluidsshockformationenergyestimatesdustandradiation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to establish when a homogeneous, barotropic perfect fluid on a fixed expanding universe with scale factor $a(t)=t^\alpha$, $\alpha<1$, remains smooth for all future times instead of forming shocks. Its main theorem proves nonlinear stability of the quiet state $(L,\vec v)=(L_0,0)$ in expansion-normalized variables under the explicit condition $K<1-2/(3\alpha)$, where $K=c_s^2$ is the squared speed of sound and $\alpha$ the expansion exponent. The stability range is the slowest expansion in which the friction-like damping can overcome the sound-speed tendency to focus waves, and the proof exhibits a corrected energy functional that realizes this balance. The paper also proves shock formation from arbitrarily small data for dust with $\alpha\le 1/2$ and for radiation with $\alpha\le 1$, so the claimed phase boundary is flanked by genuine instability. A reader should care because this pins down, for the first time in the spatially compact decelerated regime, which cosmological fluids can be globally smooth perturbations of homogeneity.

What carries the argument

The load-bearing object is the corrected $L^2$ energy $E_{\ell,c}=E_\ell+\frac{1}{2}\frac{\alpha(1-3K)}{1+K}t^{\alpha-1}\int(D^\ell v,D^{\ell+1}L)$, whose correction term produces a definite decay for the spatial gradient of the rescaled density while the friction term gives decay for $\bar v$ and $Dv$. The optimal correction coefficient balances the extra positive contribution against the damping and is exactly where the condition $K<1-2/(3\alpha)$ enters; a hierarchical set of energy identities keeps the estimate closed in regularity, and separate Riccati and characteristic arguments produce the shock examples.

What would settle it

Check the estimate (4.15) numerically or analytically for the lowest allowed regularity $N=6$: construct initial data with $\|Dv\|_{H^6}+\|DL\|_{H^6}$ small and integrate the corrected energy inequality; if the corrected energy $E_{N,c}$ fails to remain equivalent to $E_N$ or the energy does not decay like $t^{-(\alpha(1-3K)-\mu)/2}$, the bootstrap in Theorem 5.1 collapses. A simpler witness would be a direct failure of the claimed bound $|J|\le(\ell+2)/2$ for a distribution of derivatives in a higher-order perturbative term.

Watch

Extended reading notes

Core claim

On the torus $\mathbb{T}^3$ with metric $-dt^2+a(t)^2\delta_{ij}dx^idx^j$ and $a(t)=t^\alpha$, the quiet solution $(L,\vec v)=(L_0,0)$ of the expansion-normalized relativistic Euler equations is asymptotically stable for $K\in(0,1/3)$ when $(1-K)\alpha>2/3$. Theorem 5.1 gives global smooth solutions for $N\ge 6$ and decay rates $|\bar v(t)|+\|Dv\|_{H^N}+\|DL\|_{H^N}\le C\varepsilon t^{-(\alpha(1-3K)-\mu)/2}$, with the same decay transferred to the physical variables in Corollary 5.2. The stability curve $K=1-2/(3\alpha)$ is claimed to be sharp: companion numerics indicate instability when it is violated, and the paper's own shock-formation results place the dust transition at $\alpha=1/2$ and the radiation transition at $\alpha=1$.

Load-bearing premise

The energy argument depends on every non-explicit, higher-order error term in the commuted energy identities satisfying the stated bounds (4.5)/(4.14) and being absorbed together with the sign-indefinite correction term by the decay terms, a bookkeeping step that the paper presents only as a sketch for $N\ge 6$.

Editorial extensions

If this is right

  • Under the theorem, small perturbations of the quiet fluid on a decelerating torus universe stay smooth and the rescaled velocity and density gradients decay polynomially in $t$, with the density perturbation itself bounded by a constant times the initial data size.
  • The decay exponent $\alpha(1-3K)-\mu$ is positive exactly when $K<1-2/(3\alpha)$, so the theorem cannot be extended to slower expansion without either a smaller sound speed or a stronger decay mechanism.
  • Setting $\alpha=1$ in the corrected energy yields a streamlined proof of the known linear-expansion stability theorem for $K\in(0,1/3)$.
  • For dust, stability above $\alpha=1/2$ and shock formation below or at that exponent combine to locate the dust phase transition at $\alpha=1/2$; for radiation the transition sits at $\alpha=1$.
  • If the companion numerics are correct, the inequality (1.3) is not merely sufficient but marks the true phase boundary between stabilization and shock formation.

Reading between the lines

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

  • A natural next test is to probe the curve $K=1-2/(3\alpha)$ from the unstable side: the numerics of the companion paper predict blow-up for any $K>0$ when $\alpha$ lies strictly between $1/2$ and $2/3$, which would make the $K\to 0$ limit of the stability diagram discontinuous.
  • The same correction mechanism may extend to inhomogeneous or self-gravitating settings, where the density-gradient correction term would have to compete with curvature errors; locating its analogue could yield stability thresholds for Einstein-Euler systems in decelerated cosmologies.
  • Because the shock-formation examples are arbitrarily small, the linearized decay rates cannot be the whole story: the nonlinear focusing term is what breaks the solution, so a purely spectral stability analysis would miss the true boundary.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The paper studies the relativistic Euler equations with a linear barotropic equation of state p = Kρ on fixed FLRW spacetimes with scale factor a(t) = t^α in the decelerated regime α < 1. In expansion-normalized variables (L, v), the authors construct a corrected L^2-based energy functional and prove nonlinear stability of the quiet fluid solution (L_0, 0) under the condition K < 1 - 2/(3α), with explicit power-law decay rates for the mean velocity and higher Sobolev norms. The same energy functional is adapted to the linear expansion case α = 1, yielding a new proof of previously known stability for 0 < K < 1/3. The paper also contains two shock-formation results: radiation fluids form shocks for arbitrarily small perturbations when α ≤ 1, and dust forms shocks for arbitrarily small perturbations when α ≤ 1/2. The stability proof relies on a sequence of energy identities for the uncorrected and corrected energies, followed by a bootstrap argument in Theorem 5.1.

Significance. If the main theorem is correct, it is the first nonlinear stability result in the spatially compact decelerated regime for fluids with non-vanishing speed of sound, and it identifies a concrete threshold curve K = 1 - 2/(3α). The paper gives an explicit mechanism for the decay: a correction term generates a negative contribution proportional to (1 - 3Kα)t^{α-2}∫v·∂L, which is balanced against error terms. The algebraic cancellations in Proposition 4.8 and the sign structure in Lemma 4.11 are convincing and are among the paper's strengths. The adaptation to linear expansion in Section 6 is elegant, and the shock-formation results in Section 7 complement the stability theorem by showing that the dust threshold α = 1/2 is sharp and that radiation is always unstable in the decelerated range. The main weakness is that the higher-order error bookkeeping in Lemma 4.17, which is the load-bearing step for the full Sobolev estimate, is only sketched.

major comments (2)
  1. [§4.4.3, Lemma 4.17] Lemma 4.17 is the step that converts the lowest-order cancellation into a uniform estimate at all derivative levels, but its proof is presented as a sketch: only the representative term (4.15) is estimated, with a single derivative split, and the text states that the remaining terms follow similarly. In particular, the borderline case |I| = |J| = (ℓ+2)/2 at ℓ = N = 6 is not treated, and the terms in which an undifferentiated factor v is decomposed as v = \bar v + (v - \bar v) require a separate argument because \bar v is not pointwise small. Since Lemma 4.17 supplies the absorption of error terms needed to close the bootstrap in Theorem 5.1, a complete verification of these cases is necessary for the decay estimate (5.2) and hence for the threshold (1.3).
  2. [Definition 4.16 and Theorem 5.1 proof] Definition 4.16 defines a higher-order perturbative term by the bound (4.14), which contains no time weight, while the bootstrap in Theorem 5.1 uses error terms of size C(1+t^{1-α})t^{-1}E_{N,c}^{3/2} and C t^{α-1}t^{-1}E_{N,c}. The paper does not explain how these time weights follow from (4.14) after commuting with D^I, nor does it quantify the smallness of t0 and ε required to absorb the sign-indefinite term -c(1-3Kα)t^{α-2}∫(D^ℓ v, D^{ℓ+1}L). If any omitted term were quadratic in E_{N,c} rather than E_{N,c}^{3/2}, or if the absorption required t0^{α-1} to be small in a way not guaranteed by the stated hypotheses, the decay (5.2) and the threshold (1.3) would not follow.
minor comments (4)
  1. [§1.1] There is a typo: 'consistute' should be 'constitute'.
  2. [Lemma 4.9 proof] The proof of Lemma 4.9 is very brief; given that the lemma contains explicit coefficients, it would help to list the terms that are declared perturbative and to indicate which cancellations remove the remaining explicit contributions.
  3. [§7.3] The sentence 'nothing about the intersections of the characteristics of (7.15) can be said at this point' is immediately followed by the enveloping construction; the logical transition would be clearer if the role of the envelopes as upper and lower bounds were stated before drawing the intersection conclusion.
  4. [§4.4.3, proof of Lemma 4.17] The phrase 'Sobolov-embedding' contains a typo; it should be 'Sobolev embedding'.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the stability threshold (1.3) is obtained from the corrected-energy balance, and author self-citations are used only as supporting or contextual evidence.

full rationale

The central derivation is self-contained. The stability threshold (1.3) is produced inside the paper by the balance between the decay term -alpha(1-3K)/t E and the error terms after the correction term with coefficient c = alpha(1-3K)/(2(1+K)) is added; Theorem 5.1 then requires alpha(1-3K)-mu > 2(1-alpha), which is algebraically equivalent to (1-K)alpha > 2/3, i.e. K < 1 - 2/(3alpha). No parameter in this chain is fitted to the companion numerics, and the companion paper [4] is cited only for the hedged claim that the condition is 'presumably sharp', not as an input to the proof of Theorem 5.1. Other author self-citations [5,6,7] provide prior results and context; Section 6 independently reproves the known linear-expansion stability theorem using the same energy, which is a consistency check rather than a circular input. The shock-formation sections use the characteristic method and cite external results [20,23] for comparison, not in-paper fitted values. Thus no prediction reduces by construction to a fitted parameter or to a self-citation chain. Two non-circular proof gaps are worth noting for completeness: Definition 4.16 omits the time-weight factor (1+t^{1-alpha})/t that appears in the bootstrap, and Lemma 4.17's proof checks only the representative term (4.15), leaving the higher-order bookkeeping sketched; these affect verifiability, not circularity.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The stability proof is self-contained apart from standard analytic tools: the Poincaré inequality on T^3 and Sobolev embedding estimates are invoked with references (Lemma 2.3 from Hebey [12]), and local existence with the continuation principle is imported from Speck [23]. The domain of the theorem is restricted to the polynomial scale factor a(t) = t^α with α < 1 and the linear equation of state p = Kρ, both stated in Section 1. There are no parameters fitted to data: the correction coefficient c = α(1-3K)/(2(1+K)) is fixed by the energy balance in Lemma 4.11, and the small margin µ in Theorem 5.1 is any positive number below α(1-3K) - 2(1-α). No new physical entities are introduced; the corrected energy functional is a proof device. The only external input used outside the proof is numerical evidence for sharpness from the companion paper [4], cited as 'suggesting' rather than proving the converse.

assumptions (6)
  • domain assumption Local-in-time existence and the continuation principle for (3.4): a solution extends beyond any t1 > t0 whenever its H^{N+1} Sobolev norm stays bounded.
    Invoked in Section 3 after (3.6) ('can be extended beyond any t1 > t0 if its Sobolev norm of order N+1 remains bounded [23]'); it converts the uniform energy bound from Section 5 into global existence.
  • standard math Poincaré inequality on the compact torus (Lemma 2.3, cited from Hebey, Theorem 2.10).
    Used throughout Section 4 to control the mean velocity bar v and the fluctuations v - bar v, e.g., in Lemma 4.1 and Lemma 4.17.
  • domain assumption Restriction to the polynomial scale factor a(t) = t^α in the decelerated range α < 1; the authors state more general scale factors are expected but not treated.
    Section 1.1: 'we restrict to a polynomial scale factor a(t) = t^α but we expect more general results could be derived in future work.' The theorem's domain of validity is exactly this class.
  • ad hoc to paper Bootstrap assumptions t > 1 and |v| < 1/10 are propagated by the energy estimates for sufficiently small data.
    Section 3 ('We make the following choices and bootstrap assumptions for convenience: t > 1, |v| < 1/10'); closure is the content of the proof of Theorem 5.1, and Lemma 4.19 additionally needs t large and ||v|| small.
  • domain assumption Linear barotropic equation of state p = Kρ with 0 < K < 1/3 for the stability theorem; dust (K = 0) and radiation (K = 1/3) are handled as separate cases in the shock-formation section.
    The EOS is imposed in (1.1) and Theorem 5.1; the K → 0 discontinuity noted in Section 1.6 is a consequence of treating dust separately rather than as a limit of the stability theorem.
  • standard math Sobolev embedding and Cauchy-Schwarz estimates for the perturbative terms (Definitions 4.3 and 4.16), including the borderline term estimate (4.15).
    The whole energy bookkeeping rests on these bounds; the N ≥ 6 condition is chosen so that ||D^J L||_{L∞} is controlled, as shown in the proof of Lemma 4.17.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Stability of fluids in spacetimes with decelerated expansion." pith.science (2026). https://pith.science/paper/SFYUPCNK

@misc{pith2026250112798,
  author       = {Pith},
  title        = {Pith review of: Stability of fluids in spacetimes with decelerated expansion},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SFYUPCNK}},
  note         = {Machine review of arXiv:2501.12798}
}
read the original abstract

We prove the nonlinear stability of homogeneous barotropic perfect fluid solutions in fixed cosmological spacetimes undergoing decelerated expansion. The results hold provided a specific inequality between the speed of sound of the fluid and the expansion rate of spacetime is valid. Numerical studies in our earlier complementary paper provide strong evidence that the aforementioned condition is sharp, i.e. that instabilities occur when the inequality is violated. In this regard, our present result covers the regime of slowest possible expansion which allows for fluids to stabilize, depending on their speed of sound. Our proof relies on an energy functional which is universal in the sense that it also applies to the case of linear expansion and enables a significantly simplified proof of bounds for fluids on linearly expanding spacetimes. Finally, we consider the special cases of dust and radiation fluids in the decelerated regime and prove shock formation for arbitrarily small perturbations of homogeneous solutions.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Boundedness and decay of waves on spatially flat decelerated FLRW spacetimes

    gr-qc 2025-05 conditional novelty 6.0 of 10

    A twisted vector-field method yields energy boundedness, local energy decay, r^p-weighted estimates, and energy and pointwise decay for waves on all spatially flat decelerated FLRW backgrounds with scale factor t^q, 0<q<1.

Reference graph

Works this paper leans on

23 extracted references · 22 canonical work pages · cited by 1 Pith paper

  1. [1]

    Beyer, E

    F. Beyer, E. Marshall, and T. A. Oliynyk. Future instabil ity of flrw fluid solutions for linear equations of state p = kρ with 1 / 3 < k < 1. Phys. Rev. D , 107:104030, May 2023

    work page 2023

  2. [2]

    Brauer, A

    U. Brauer, A. Rendall, and O. Reula. The cosmic no-hair th eorem and the non-linear stability of homogeneous Newtonian cosmological models. Classical Quantum Gravity , 11(9):2283–2296, 1994

    work page 1994

  3. [3]

    Christodoulou

    D. Christodoulou. The Formation of Shocks in 3-Dimensional Fluids . EMS, 2007

    work page 2007

  4. [4]

    Fajman, M

    D. Fajman, M. Maliborski, M. Ofner, T. Oliynyk, and Z. Wya tt. Phase transition between shock formation and stability in cosmological fluids. 2024. preprint [arXiv:24 05.03431]

    work page 2024

  5. [5]

    Fajman, M

    D. Fajman, M. Ofner, T. A. Oliynyk, and Zoe Wyatt. The stab ility of relativistic fluids in linearly expanding cosmologies. Int. Math. Res. Not. IMRN , (5):4328–4383, 2024

    work page 2024

  6. [6]

    Fajman, M

    D. Fajman, M. Ofner, and Z. Wyatt. Slowly expanding stabl e dust spacetimes. Arch. Rat. Mach. Anal., 248(5):66, 2024

    work page 2024

  7. [7]

    Fajman, T.A

    D. Fajman, T.A. Oliynyk, and Zoe Wyatt. Stabilizing rela tivistic fluids on spacetimes with non-accelerated expansion. Commun. Math. Phys. , 383:401–426, 2021

    work page 2021

  8. [8]

    Fournodavlos, E

    G. Fournodavlos, E. Marshall, and T. A. Oliynyk. Future s tability of perfect fluids with extreme tilt and linear equation of state p = c2 sρ for the Einstein-Euler system with positive cosmological c onstant: The range 1 3 < c 2 s < 3 7 , 2024

    work page 2024

Show all 23 references
  1. [9]

    Friedman

    A. Friedman. ¨Uber die Kr¨ ummung des Raumes. Zeitschrift f¨ ur Physik, 10(1):377–386, 1922

    1922

  2. [10]

    Friedrich

    H. Friedrich. Sharp asymptotics for Einstein- λ -dust flows. Comm. Math. Phys. , 350:803 – 844, 2017

    2017

  3. [11]

    Hadˇ zi´ c and J

    M. Hadˇ zi´ c and J. Speck. The global future stability of the FLR W solutions to the Dust-Einstein system with a positive cosmological constant. J. Hyper. Differential Equations , 12:87–188, 2015

    2015

  4. [12]

    E. Hebey. Nonlinear analysis on manifolds: Sobolev spaces and inequa lities, volume 5 of Courant Lecture Notes in Mathematics . New York University, Courant Institute of Mathematical Sc iences, New York; American Math- ematical Society, Providence, RI, 1999

    1999

  5. [13]

    P. G. LeFloch and C. Wei. Nonlinear stability of self-gr avitating irrotational Chaplygin fluids in a FLR W geom- etry. Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 38(3):787–814, 2021

    2021

  6. [14]

    Liu and C

    C. Liu and C. Wei. Future stability of the FLR W spacetime for a large class of perfect fluids. Ann. Henri Poincar´ e, 22(3):715–770, 2021

    2021

  7. [15]

    L¨ ubbe and J

    C. L¨ ubbe and J. A. Valiente Kroon. A conformal approach for the analysis of the non-linear stability of radiation cosmologies. Annals of Physics , 328:1–25, 2013

    2013

  8. [16]

    Marshall and T

    E. Marshall and T. A. Oliynyk. On the stability of relati vistic perfect fluids with linear equations of state p = kρ where 1 / 3 < k < 1. Letters in Mathematical Physics , 113(5):102, 2023

    2023

  9. [17]

    P. Mondal. The linear stability of the n+1 dimensional FLR W spacetimes. Classical Quantum Gravity, 38(22):Pa- per No. 225009, 51, 2021

    2021

  10. [18]

    T. A. Oliynyk. Future stability of the FLR W fluid solutio ns in the presence of a positive cosmological constant. Commun. Math. Phys. , 346:293–312; see the preprint [arXiv:1505.00857] for a co rrected version, 2016

    2016 arXiv

  11. [19]

    Todd A. Oliynyk. Future global stability for relativis tic perfect fluids with linear equations of state p = Kρ where 1 / 3 < K < 1/ 2. SIAM J. Math. Anal. , 53(4):4118–4141, 2021

    2021

  12. [20]

    Rendall and F

    A.D. Rendall and F. Stahl. Shock waves in plane symmetri c spacetimes. Communications in Partial Differential Equations, 33:2020–2039, 2008

    2020

  13. [21]

    Rodnianski and J

    I. Rodnianski and J. Speck. The stability of the irrotat ional Euler-Einstein system with a positive cosmological constant. J. Eur. Math. Soc. , 15:2369–2462, 2013

    2013

  14. [22]

    J. Speck. The nonlinear future-stability of the FLR W fa mily of solutions to the Euler-Einstein system with a positive cosmological constant. Selecta Mathematica, 18:633–715, 2012

    2012

  15. [23]

    J. Speck. The stabilizing effect of spacetime expansion on relativistic fluids with sharp results for the radiation equation of state. Arch. Rat. Mech. , 210:535–579, 2013. David F ajman & Maximilian Ofner, Todd Oliynyk Zoe Wyatt F aculty of Physics, School of Mathematics Depart...

    2013

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.