REVIEW 2 major objections 4 minor 25 references
Pathological solutions of Navier-Stokes equations on $\mathbb{T}^2$ with gradients in Hardy spaces
T0 review · 2 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read The paper proves that p=1 is the sharp boundary between uniqueness and non-uniqueness for very weak 2D Navier-Stokes solutions when vorticity is measured in Hardy spaces, and constructs the p<1 pathology for every smooth initial datum.
desk verdict Sharp p=1 threshold for 2D NSE non-uniqueness in C(H^p) vorticity path spaces — structurally sound, but the new periodic Hardy-space estimates are load-bearing and need careful referee scrutiny. 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 engine is a single iteration step for the Navier-Stokes form of the Reynolds system, where a solution with a smooth stress error R0 is replaced by a new solution whose stress error is arbitrarily small in L1 and whose perturbation gradient is small in H^p. The perturbation is built from fast oscillating, temporally concentrated building blocks, and a crucial point is that the principal corrector is written with enough derivatives so that, after applying the periodic Hardy-space atomic estimates, its gradient lands in H^p for any p<1. Around this step the paper develops the needed periodic Hardy-space theory: maximal-function definition, atomic decomposition, the fact that on the torus at
What would settle it
Compute the H^p norm (with p<1) of the gradient of a concentrated, fast-oscillating corrector of the form used in Section 4, either analytically or numerically on the torus, and check the bound (20): a counterexample where the big-atom estimate of Corollary 3.6 or the projection estimate of Proposition A.17 gives a strictly worse power of the frequency would refute the iteration. More directly, one could search for a periodic function supported on a large ball that is an H^p atom under Definition 3.5 but whose H^p norm is unbounded in the ball's size.
Extended reading notes
Core claim
On the two-dimensional torus, take any smooth divergence-free initial datum. The paper's central claim is a dichotomy for very weak solutions with vorticity in C([0,T],H^p): for p≥1 the Cauchy problem is unique there, while for every 0<p<1 there are infinitely many distinct very weak solutions, all starting from the same datum, with gradients in C([0,T],H^p). It proves an even stronger flexibility statement: any two smooth solutions can be joined by a single nonzero very weak solution that equals the first on [0,1/8] and the second on [7/8,1]. The threshold is not about the smoothness of the initial value alone: the same smooth datum sits in every H^p, so the deciding object is the path spac
Load-bearing premise
The load-bearing premise is that the newly developed periodic Hardy-space estimates are correct as used—in particular, that on the torus an atom supported on a large ball needs no vanishing moments, and that the divergence-free projection is bounded on H^p; if either estimate fails, the p<1 construction collapses.
Editorial extensions
If this is right
- The index p=1 becomes a sharp threshold: measurable vorticity in H^1 is enough for uniqueness, while any p<1 admits non-unique solutions.
- Every smooth initial datum supports infinitely many pathological solutions, so non-uniqueness cannot be blamed on rough data.
- Any two smooth solutions can be interpolated by a single weak solution that follows one early and the other late—a strong form of flexibility.
- The constructed solutions have gradient controlled uniformly in time in H^p and lie in W^{σ,1} for every σ<1, improving the temporal control of previous non-unique 2D Navier-Stokes solutions.
- Well-posedness statements for these equations should specify the path space for vorticity, not merely the initial-data space.
Reading between the lines
- If the periodic Hardy-space estimates carry over to higher dimensions, the same error-reduction step might yield non-unique 3D Navier-Stokes or Euler solutions with vorticity in H^p for p<1; the paper does not claim this.
- The H^1 threshold suggests a natural open question: does initial vorticity in H^1 produce a mild solution whose vorticity remains in C([0,T],H^1)? The paper notes this is not addressed.
- If the big-atom property is robust, similar constructions could work on other compact manifolds, where atoms on large balls also need no cancellation.
- A direct check of gradient estimate (20) on simpler oscillatory test functions would allow the Hardy-space bounds to be tested numerically before relying on the full iteration.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper establishes a sharp dichotomy for very weak solutions of the 2D Navier-Stokes equations on T^2, measured through the path space of vorticity. For p ≥ 1, if two very weak solutions have vorticities in C([0,T], H^p) and agree at time zero, then they coincide (Theorem 1.2(A)). For p ∈ (0,1), the paper claims that for any smooth initial datum there are infinitely many very weak solutions with gradients in C([0,T], H^p), and more generally that any two smooth solutions can be connected by a nonzero solution agreeing with the first on [0,1/8] and with the second on [7/8,1] (Theorems 1.2(B) and 1.3). The construction follows the convex integration scheme for the NS-Reynolds system, with a key iteration Proposition 2.2. The new ingredient is a theory of real Hardy spaces on T^n, developed in the appendix, including atomic decompositions and the boundedness of the Leray projector on H^p(T^n).
Significance. If the periodic Hardy-space toolbox is fully justified, the paper gives a clean and striking threshold: p = 1 is exactly the boundary between uniqueness and non-uniqueness when the vorticity path space is C([0,T], H^p). The convex integration part is standard in structure, but the explicit parameter hierarchy in Section 6 is transparent and checkable. The appendix is a valuable contribution on its own, especially the periodic atomic decomposition and the observation that periodic atoms need not have vanishing moments when supported on large balls. The paper also makes explicit that the ill-posedness phenomenon is a property of the path space, not merely of the initial-data space. No machine-checked proofs are supplied, but the main computation is written in sufficient detail to be audited.
major comments (2)
- [Appendix A.5, Proposition A.17] Proposition A.17 is load-bearing: it is used in §4.3.2 to control ∇w_t and hence enters the gradient estimate (20). The transference proof has two gaps. First, the displayed chain after (35) bounds ∥sup_{0<ζ≤R}|Ψζ*Pa|∥_{L^p(T^n)} by N^{-n}∫_{[-N/2,N/2]^n} sup |Ξ_{1/N}(Ψζ*Pa)_ext|^p dx. Since |Ξ_{1/N}|≤1, this integrand is pointwise smaller than |Ψζ*Pa|; to obtain an upper bound one must use periodicity to show that the average of |Ξ_{1/N}|^p over integer shifts is bounded below by a positive constant independent of N. No such argument is supplied. Second, E_N(x)→0 is asserted uniformly for ζ∈(0,R], but the final estimate requires uniformity as R→∞; the large-ζ regime is not treated. The proof of Proposition 3.8 is therefore incomplete, and the p<1 construction depends on it.
- [Appendix A.4, Proposition 3.9] Proposition 3.9 is the mechanism used in Theorem 1.2(A) to pass from C_tH^1 vorticity to u∈C_tL^2. The proof invokes the statement 'In two dimensions H^1 ,→(H^1)^*' to control the distribution pairing. As written this is not a standard fact: H^1 denotes both the Hardy space and the Sobolev space, and the asserted embedding is not justified. The desired bound can be obtained directly from H^1⊂L^1, div f=0, the elliptic estimate giving ∇f∈L^1, and the Sobolev embedding W^{1,1}(T^2)↪L^2. The conclusion is correct, but the present proof needs to be replaced or substantially clarified.
minor comments (4)
- [§4.3.2, final estimate] The text says 'Remark 3.6 and (19)' but no Remark 3.6 exists; the relevant result is the embedding L^s⊂H^p for s>1 (Proposition 3.3 together with Lemma A.4). Please cite the correct statement.
- [Theorem 1.2(B) proof] The proof that there are 'infinitely many' distinct solutions is only implicit. Since Theorem 1.3 gives one solution for each choice of v2, one should explicitly argue that different choices of v2 (e.g., different final data) produce different solutions on [7/8,1] and hence distinct u. Please make this step explicit.
- [Appendix A.4, notation] In the proof of Proposition 3.9, the same symbol H^1 is used for the Hardy space and for the Sobolev space H^1(T^2). This is confusing and should be resolved by notation or a sentence.
- [Throughout] Minor typos: 'a ill-posedness space' should be 'an ill-posedness space'; §4.3.2 refers to 'Corollary 3.8' when Proposition 3.8 is meant; the proof of Proposition A.17 overloads P for both periodic and full-space Leray projectors.
Circularity Check
No significant circularity: the p=1 threshold is assembled from an external uniqueness theorem and an external convex-integration construction; the periodic Hardy-space appendix is new, not self-referential.
full rationale
The paper's derivation chain does not reduce any of its claims to its own assumptions. The well-posedness side (Theorem 1.2(A)) depends on the external Ladyzhenskaya-Prodi-Serrin uniqueness criterion as stated in Cheskidov-Luo [12], after a short reduction using Proposition 3.9. Proposition 3.9 is proved in Appendix A.4 from a Helmholtz-style decomposition; regardless of whether that proof is fully detailed, it does not presuppose the uniqueness being proved. The ill-posedness side (Theorem 1.2(B), Theorem 1.3) is a convex-integration iteration built on Proposition 2.2. The gradient-in-H^p estimate (20) is derived by applying the periodic Hardy-space tools of the appendix (Corollary 3.6, Proposition 3.7, Proposition 3.8/A.17) to explicit correctors. Those appendix results are proved from the maximal-function and atomic definitions in the paper, not assumed from the target theorem. There is no fitted parameter renamed as a prediction: the iteration parameters are chosen explicitly and all error terms are shown to be small by negative powers of λ. The paper does cite the first author's earlier work [10] for the antidivergence Lemma 3.1, but that is a published, independent technical lemma and not the source of the threshold dichotomy. The Skeptic's concern about the correctness of Corollary 3.6, Proposition 3.7, and Proposition A.17 is a verification risk, not circularity: the estimates are not equivalent to the theorem's conclusion by construction, and no displayed equation in the appendix assumes the p<1 non-uniqueness result. The paper's own footnote 2 noting that existence of mild solutions for p=1 is open is a limitation, not a circular step. Therefore no circularity step can be exhibited, and the circularity score is 0.
Assumptions & free parameters
free parameters (1)
- Convex integration parameter hierarchy (exponents alpha, a, b, beta, gamma and integrability s) =
a=5, b=11, beta=4, gamma=1, s close to 1, alpha large (alpha > 2/p, plus Section 6 inequalities)
assumptions (7)
- domain assumption LPS-type uniqueness theorem for very weak 2D NSE solutions (Theorem 1.3 of Cheskidov-Luo [12])
- domain assumption Global smooth solvability of 2D NSE for smooth periodic data
- domain assumption Building-block estimates (7), corrector splitting (Proposition 6.1), and decorrelation (Lemma 7.2, 7.3) from Buck-Modena [5]
- standard math Classical H^p(R^n) theory of Stein [23] and Folland-Stein [19] (maximal function equivalence, atomic decomposition)
- domain assumption Antidivergence operators of Burczak-Modena-Szekelyhidi [10, Proposition 4] (Lemma 3.1)
- standard math Nash-type decomposition of near-identity symmetric matrices into rank-one pieces (Lemma 3.2)
- standard math Hodge decomposition with W^{1,2} potentials and the 2D embedding W^{1,2} into BMO
Cite this review
Pith. "Pith review of Pathological solutions of Navier-Stokes equations on $\mathbb{T}^2$ with gradients in Hardy spaces." pith.science (2026). https://pith.science/paper/ALS7KGBY
@misc{pith2026250908168,
author = {Pith},
title = {Pith review of: Pathological solutions of Navier-Stokes equations on $\mathbbT^2$ with gradients in Hardy spaces},
year = {2026},
howpublished = {\url{https://pith.science/paper/ALS7KGBY}},
note = {Machine review of arXiv:2509.08168}
}
abstract
For an arbitrary smooth initial datum, we construct multiple nonzero solutions to the $2$d Navier-Stokes equations, with their gradients in the Hardy space $\mathcal{H}^p$ with any $p \in (0,1)$. Thus, in terms of the path space $C(\mathcal{H}^p)$ for vorticity, $p=1$ is the threshold value distinguishing between non-uniqueness and uniqueness regimes. In order to obtain our result, we develop the needed theory of Hardy spaces on periodic domains.
Figures
Reference graph
Works this paper leans on
-
[5]
Compactly supported anomalous weak solutions for 2D Euler equations with vorticity in Hardy spaces
M. Buck and S. Modena. Compactly supported anomalous weak solutions for 2D Euler equations with vorticity in Hardy spaces.ArXiv 2405.19214, 2024
work page Pith review arXiv 2024
-
[1]
Albritton, E
D. Albritton, E. Bru´ e, and M. Colombo. Non-uniqueness of Leray solutions of the forced Navier-Stokes equations.Ann. of Math. (2), 196(1):415–455, 2022
2022
-
[2]
D. Albritton, E. Bru´ e, M. Colombo, C. De Lellis, V. Giri, M. Janisch, and H. Kwon. Instability and non-uniqueness for the 2D Euler equations, after M. Vishik, volume 219 of Annals of Mathematics Studies. Princeton University Press, Princeton, NJ, 2024
work page 2024
-
[3]
E. Bru´ e and M. Colombo. Nonuniqueness of solutions to the Euler equations with vorticity in a Lorentz space.Communications in Mathematical Physics, 403(2):1171–1192, 2023
work page 2023
- [4]
-
[6]
M. Buck and S. Modena. Non-uniqueness and energy dissipation for 2D Euler equations with vorticity in Hardy spaces.Journal of Mathematical Fluid Mechanics, 26(2):1–39, 2024
work page 2024
-
[7]
Buckmaster, M
T. Buckmaster, M. Colombo, and V. Vicol. Wild solutions of the Navier-Stokes equations whose singular sets in time have Hausdorff dimension strictly less than 1.J. Eur. Math. Soc. (JEMS), 24(9):3333–3378, 2022
2022
-
[8]
T. Buckmaster and V. Vicol. Nonuniqueness of weak solutions to the Navier-Stokes equa- tion.Ann. of Math. (2), 189(1):101–144, 2019
work page 2019
Show all 25 references
-
[9]
Buckmaster and V
T. Buckmaster and V. Vicol. Convex integration constructions in hydrodynamics.Bull. Amer. Math. Soc. (N.S.), 58(1):1–44, 2021
2021
-
[10]
Burczak, S
J. Burczak, S. Modena, and L. Sz´ ekelyhidi. Non uniqueness of power-law flows.Commu- nications in Mathematical Physics, 388:199–243, 2021
2021
-
[11]
Castro, D
A. Castro, D. Faraco, F. Mengual, and M. Solera. A proof of Vishik’s nonuniqueness theorem for the forced 2d Euler equation.arXiv:2404.15995
-
[12]
Cheskidov and X
A. Cheskidov and X. Luo. Sharp nonuniqueness for the Navier-Stokes equations.Inventiones mathematicae, 229(3):987–1054, 2022
2022
-
[13]
Cheskidov and X
A. Cheskidov and X. Luo.L 2-critical nonuniqueness for the 2D Navier-Stokes equations. Annals of PDE, 9(2):13, 2023
2023
-
[14]
M. P. Coiculescu and S. Palasek. Non-uniqueness of smooth solutions of the navier-stokes equations from critical data.arXiv:2503.14699
-
[15]
De Lellis and L
C. De Lellis and L. Sz´ ekelyhidi, Jr. The Euler equations as a differential inclusion.Ann. of Math. (2), 170(3):1417–1436, 2009
2009
-
[16]
De Lellis and L
C. De Lellis and L. Sz´ ekelyhidi, Jr. Weak stability and closure in turbulence.Philos. Trans. Roy. Soc. A, 380(2218):Paper No. 20210091, 16, 2022
2022
-
[17]
T. M. Elgindi, R. W. Murray, and A. R. Said. Wellposedness and singularity formation beyond the Yudovich class.arXiv:2312.17610, 2023
2023 arXiv
-
[18]
Fan and S
D. Fan and S. Sato. Transference on certain multilinear multiplier operators.Journal of the Australian Mathematical Society, 70(1):37–55, 2001
2001
-
[19]
G. B. Folland and E. M. Stein.Hardy Spaces on Homogeneous Groups.(MN-28), Volume 28, volume 28. Princeton University Press, 2020
2020
-
[20]
Goldberg
D. Goldberg. A local version of real Hardy spaces.Duke Mathematical Journal, 46(1):27–42, 1979
1979
-
[21]
P. Isett. A proof of Onsager’s conjecture.Ann. of Math. (2), 188(3):871–963, 2018
2018
-
[22]
Jia and V
H. Jia and V. Sverak. Are the incompressible 3d navier–stokes equations locally ill-posed in the natural energy space?Journal of Functional Analysis, 268(12):3734–3766, 2015
2015
-
[23]
E. Stein. Harmonic analysis: real-variable methods, orthogonality, and oscillatory integrals. Princeton Mathematical Series, 43, 1993. 30
1993
-
[24]
M. Vishik. Instability and non-uniqueness in the Cauchy problem for the Euler equations of an ideal incompressible fluid. Part I.arXiv:1805.09426, 2018
2018 arXiv
-
[25]
M. Vishik. Instability and non-uniqueness in the Cauchy problem for the Euler equations of an ideal incompressible fluid. Part II.arXiv:1805.09440, 2018. 31
2018 arXiv
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.