REVIEW 3 major objections 4 minor 36 references
On the sharp critical mass threshold for the 3D Patlak-Keller-Segel-Navier-Stokes system via Couette flow
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A strong enough Couette flow suppresses 3D chemotaxis blow-up whenever the total cell mass stays below $16\pi^{2}$, a threshold the paper argues is sharp.
desk verdict A credible global-existence theorem for 3D PKS-NS with mass below 16π², but the 'sharp' threshold in the title is heuristic, not proved; the paper deserves review after softening that claim. 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 object is the zero mode $n_{0}(t,y,z)$, the $x$-average of the cell density, whose dynamics (1.8) is a 2D Keller–Segel-type diffusion–aggregation equation with extra 'good terms' from the velocity and the non-zero modes. The proof decomposes every field into the zero mode (depending on $(y,z)$) and the non-zero mode (depending on $(x,y,z)$): non-zero modes feel the mixing of the Couette flow and decay at the enhanced rate $e^{-aA^{-1/3}t}$, while the non-constant part of the zero modes obeys the heat-type dissipation of Lemmas 4.2–4.3. The boundedness of the density is obtained from a free-energy estimate for $n_{0}$ using the logarithmic Hardy–Littlewood–Sobolev inequality, which yields a finite $\|n_{0}\log^{+}n_{0}\|_{L^{1}}$ bound precisely when the rescaled mass $m = M/(2\pi)$ is below $8\pi$; from there, a Moser–Alikakos iteration gives $\|n\|_{L^{\infty}}$. The 3D lift-up effect on $u_{1,0}$ is tamed by splitting it into a good part $U_{1}$ and a bad part $U_{2}$, and the velocity estimates are closed with the quasi-linear scheme based on the modified operator $L_{V}$ and the 'good derivative' $\partial_{z} - \kappa\partial_{y}$, with $\kappa = \partial_{z}V/\partial_{y}V$, $V = y + U_{2}/A$.
What would settle it
Take initial data satisfying (1.6) but with rescaled zero-mode mass $m = M/(2\pi)$ strictly above $8\pi$, i.e., $M > 16\pi^{2}$, with the zero-mode density $n_{0}$ close to a supercritical 2D Keller–Segel profile, and run (1.5) with arbitrarily large $A$; the 2D critical-mass theory predicts finite-time blow-up. A numerical simulation showing global regularity for such supercritical $m$, or a proof of finite-time blow-up, would settle whether the claimed threshold is genuinely sharp.
Extended reading notes
Core claim
The central claim, Theorem 1.1, is that for initial data $0 < n_{\mathrm{in}} \in H^{2}(\mathbb{T}^{3})$, $u_{\mathrm{in}} \in H^{2}(\mathbb{T}^{3})$, there exist constants $\epsilon$ and $A_{1}$ (depending on the data) such that if $A \geq A_{1}$, $\|(u_{2,\mathrm{in}})^{0}\|_{H^{2}} + \|(u_{3,\mathrm{in}})^{0}\|_{H^{1}} \leq \epsilon$, and $M = \int_{\mathbb{T}^{3}} n_{\mathrm{in}}\,dx\,dy\,dz < 16\pi^{2}$, then the solution of the rescaled perturbation system (1.5) is global in time. The threshold is proposed as sharp: when the density is independent of $x$ and the velocity perturbation vanishes, the zero-mode equation (1.8)–(1.9) becomes the parabolic–elliptic 2D Keller–Segel equation, whose critical mass $\int_{\mathbb{T}^{2}} n_{0}\,dy\,dz = 8\pi$ corresponds exactly to $M = 16\pi^{2}$. What is new here is the observation that the non-constant parts of the zero modes $u_{2,0}$, $u_{3,0}$ decay dissipatively (Lemma 4.3), which is what allows the logarithmic Hardy–Littlewood–Sobolev argument to close.
Load-bearing premise
The argument's bootstrap starts from the requirement that the initial $x$-averages of the transverse velocity components $(u_{2,\mathrm{in}})^{0}$ and $(u_{3,\mathrm{in}})^{0}$ be $\epsilon$-small in $H^{2} \times H^{1}$ (condition (1.6)); choosing $A$ large does not create this smallness, and every zero-mode estimate in the proof, including the new decay of Lemma 4.3, depends on it, so if it fails the proof of Theorem 1.1 collapses.
Editorial extensions
If this is right
- For every initial density with total mass below $16\pi^{2}$ and sufficiently small transverse zero-mode velocity, a strong enough Couette flow makes the coupled system globally regular, so no chemotactic singularity forms at any finite time.
- The proof reaches the threshold $16\pi^{2} = 2\pi \cdot 8\pi$ in the 3D setting, improving earlier thresholds of the same type, namely $\tfrac{8\pi}{9}$ and $\tfrac{24}{5}\pi^{2}$.
- The mechanism that carries the argument is the dissipative decay of the non-constant parts of the transverse velocity zero modes (Lemma 4.3), which supplies the weak damping needed for the logarithmic Hardy–Littlewood–Sobolev bound on the density zero mode.
- If the claim is correct, the critical-mass picture for chemotaxis acquires a shear-flow analogue: the Couette flow effectively renormalizes the dimension seen by the zero-mode density to $2$, preserving the $8\pi$ value measured on the $2$D torus.
Reading between the lines
- The sharpness of $16\pi^{2}$ is inferred from the reduction of the zero-mode equation to 2D Keller–Segel; the paper does not construct a blow-up example for $M \geq 16\pi^{2}$, so the optimality of the threshold is a heuristic consequence of the 2D theory rather than a proved failure above critical mass.
- Condition (1.6) is a genuine smallness restriction on the $x$-averages of the transverse velocity, independent of $A$; large $A$ cannot create this smallness, so Theorem 1.1 covers only data whose transverse zero modes start very small. A natural test is whether the same or a similar threshold can be obtained without this restriction, for instance using a time-dependent shear flow.
- The structure of the proof suggests the $16\pi^{2}$ threshold should also appear for other shear or relaxation-enhancing flows in $\mathbb{T}^{3}$, provided the zero-mode transverse velocity is damped by an analogous heat-type dissipation; checking this for a Poiseuille flow or for a sequence of alternating shear flows would be a direct test of the mechanism.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the 3D parabolic-elliptic Patlak-Keller-Segel system coupled with the incompressible Navier-Stokes equations, written in perturbation variables around the Couette flow (Ay,0,0). After rescaling time by A, the system takes the form (1.5). The main theorem (Theorem 1.1) asserts that for positive H^2 initial density, H^2 initial velocity perturbation, A sufficiently large, and the smallness condition (1.6) on the x-average of the transverse velocity components, the solution is global in time whenever the total mass M is below 16π^2. The proof is a bootstrap over five energy functionals: zero-mode velocity energies E_{1,1} and E_{1,2}, nonzero-mode density and vorticity energies E_{2,1} and E_{2,2}, the L^∞ density norm E_3, the lift-up-related energy E_4, and two auxiliary energies E_{5,1} and E_{5,2}. The closure uses enhanced dissipation for nonzero x-frequencies, heat dissipation for the non-average zero mode, the logarithmic Hardy-Littlewood-Sobolev inequality for the free energy of the x-averaged density, and a quasi-linear change of variables to control the bad component of the first velocity component. The title and abstract additionally claim that 16π^2 is the sharp critical mass threshold, based on a heuristic reduction of the zero-mode density to the 2D Keller-Segel equation.
Significance. The global-existence statement, if fully verified, is a substantial quantitative improvement over the earlier thresholds M < 8π/9 and M < 24π^2/5 in related work, and it demonstrates that a strong Couette flow suppresses chemotactic blow-up in 3D for a natural class of data. The energy structure is elaborate, and the proof contains a genuinely new dissipative-decay estimate for the non-average parts of the transverse zero-mode velocities (Lemma 4.3), combined coherently with the logarithmic Hardy-Littlewood-Sobolev inequality. The paper also makes the role of the zero-mode density explicit through the free energy argument in Section 6. However, the paper does not prove that 16π^2 is sharp: no blow-up or unboundedness result for M > 16π^2 is supplied, and the heuristic in Remark 1.1 is not a valid reduction of the full system. The reliable contribution is the conditional global-existence theorem under the stated hypotheses, not the sharpness assertion.
major comments (3)
- [§1, Remark 1.1, and the abstract] The assertion that 16π² is the sharp threshold is not established. The reduction 'u=0 and n=n(t,y,z)' is not a solution of the perturbation system (1.3)/(1.5): for u=0 with x-independent positive n, the first velocity equation becomes ∂_t u_1 + ∂_x P = n/A after rescaling, and no periodic pressure can balance an x-independent positive buoyancy force. The zero-mode equation in Section 2.2 gives the same obstruction: at u=0 one has ∂_t u_{1,0} = n_0/A, so u_1 is immediately forced to grow. Hence the exact 2D Keller-Segel critical mass 8π is not a rigorous consequence of the zero-mode dynamics, and the manuscript supplies no finite-time blow-up or unboundedness argument for M > 16π². The title, the abstract, and Remark 1.1 should be revised to claim only global existence below 16π² under the stated smallness assumptions, unless a rigorous sharpness proof is added.
- [Lemmas 3.6, A.2, A.3, Proposition A.2, and Lemma 8.1] Several load-bearing estimates are imported without proof from the authors' companion preprint [8] or are stated without proof in this manuscript. Lemma 3.6, Lemma A.2, and Proposition A.2 are used in the energy closures of Sections 4, 5, and 8 and are justified only by 'the proof is omitted' or 'can be found in [8]'. Lemma A.3 is only partially proved in the appendix. Lemma 8.1 is stated without any proof despite being used in the crucial estimates of Lemmas 8.2–8.5. Since [8] is an unpublished preprint and the current paper advertises an improved threshold, the verification is not self-contained. Please include proofs of these statements, or provide precise hypotheses and enough detail for a reader to verify the adaptation to the present setting; at minimum, Lemma 3.6 and Proposition A.2 need this treatment.
- [Theorem 1.1, condition (1.6), and the abstract] Condition (1.6) is a separate smallness assumption on the x-average of the transverse velocity components, not a consequence of large A. Lemma 4.1 and all later zero-mode controls depend on it, and without it the lift-up term in the u_1 equation cannot be controlled. The abstract states that if the Couette flow is sufficiently strong, then global existence follows from M < 16π² alone; this omits condition (1.6) and is stronger than what Theorem 1.1 proves. Please state the smallness condition explicitly in the abstract and in the introduction, and clarify that large A alone is not sufficient under the current proof.
minor comments (4)
- [Abstract and Remark 1.1] The word 'complication' in the abstract should be 'combination'; the phrase 'whose critical mass in 2D is 8π' should also be adjusted so that the heuristic discussion is explicitly labeled as non-rigorous.
- [§2.1 and Notations] The notation f^0 and f_0 is overloaded: in (2.1), f^0 denotes the average and f_0 denotes the zero mode, but later expressions such as 'n_0' sometimes mean the zero mode and sometimes the mean value. Standardizing the notation (for example, using an overline for the spatial average) would remove a recurring source of confusion.
- [Appendix A.1] Lemma A.3 states twelve estimates but only (A.2)_1 and (A.2)_2 are proved, with the comment that the rest are similar to Lemma 3.2 in [8]. Since the remaining estimates are used repeatedly in the paper, at least a sketch of the different proof structures would help the reader, especially because the case distinctions in the norms are not all identical.
- [Throughout] There are numerous typographical issues, including 'Hölder' rendered as 'H¨ older' and several inconsistent exponential factors in the proof of Lemma 4.3. A careful proofreading pass is needed before publication.
Circularity Check
No circularity: the global-existence proof is a bootstrap with assumptions distinct from the conclusion, and the sharpness heuristic is not a derived prediction.
full rationale
The main theorem is proved by a bootstrap on explicit energy functionals; no parameter is fitted to the target inequality, and the threshold M<16π² is not defined in terms of the solution. The condition enters through Lemma 6.3 as the standard 2D logarithmic Hardy-Littlewood-Sobolev threshold m=M/(2π)<8π, i.e. an external known result scaled by the x-period, not a construction that presupposes the theorem. Several technical estimates are quoted from the same authors' prior work [8], including Lemma 3.6 and Proposition A.2, but these are auxiliary lemmas with stated hypotheses that do not include the conclusion of Theorem 1.1; relying on them is a verification risk, not a circular reduction. Remark 1.1's sharpness discussion assumes u=0 and x-independent n, which is not a solution of the forced system (1.5), so the word 'sharp' is not justified by the provided argument; that is a correctness concern about an extra heuristic claim, not a circular step in the existence proof. The central derivation chain is therefore free of self-definitional, fitted-input, or self-citation circularity.
Assumptions & free parameters
free parameters (2)
- α (interpolation exponent) =
any value in (1/2, 3/4)
- a, b (exponential decay weights in X_a and X_b norms) =
0 < a < b < 2a
assumptions (6)
- standard math Poincaré inequality and Gagliardo-Nirenberg-Sobolev inequalities on T^n
- standard math Logarithmic Hardy-Littlewood-Sobolev inequality on compact 2D manifolds
- domain assumption Initial data satisfy n_in > 0, n_in ∈ H^2, u_in ∈ H^2, on T^3
- domain assumption Zero-mode smallness (1.6) and mass bound M < 16π²
- ad hoc to paper Correctness of Lemma 3.6, Lemma A.2, and Proposition A.2, imported from the companion preprint [8]
- standard math Local well-posedness of system (1.5)
Cite this review
Pith. "Pith review of On the sharp critical mass threshold for the 3D Patlak-Keller-Segel-Navier-Stokes system via Couette flow." pith.science (2026). https://pith.science/paper/CJXUYCG4
@misc{pith2026250610578,
author = {Pith},
title = {Pith review of: On the sharp critical mass threshold for the 3D Patlak-Keller-Segel-Navier-Stokes system via Couette flow},
year = {2026},
howpublished = {\url{https://pith.science/paper/CJXUYCG4}},
note = {Machine review of arXiv:2506.10578}
}
abstract
As is well-known, the solution of the Patlak-Keller-Segel system in 3D may blow up in finite time regardless of any initial cell mass. In this paper, we are interested in the suppression of blow-up and the critical mass threshold for the 3D Patlak-Keller-Segel-Navier-Stokes system via the Couette flow $(Ay, 0, 0)$. It is proved that if the Couette flow is sufficiently strong ($A$ is large enough), then the solutions for the system are global in time in the periodic domain $(x,y,z)\in\mathbb{T}^{3}$ as long as the initial cell mass is less than $16\pi^{2}$. This result seems to be sharp, since the zero-mode function (the mean value in $x-$direction) of the three dimensional density is a complication of the two-dimensional Keller-Segel equations, whose critical mass in 2D is $8\pi$. One new observation is the dissipative decay of $(\widetilde{u}_{2,0},\widetilde{u}_{3,0})$ (see Lemma 4.3 for more details), then we combine the quasi-linear method proposed by Wei-Zhang (Comm. Pure Appl. Math., 2021) with the zero-mode estimate of the density by the logarithmic Hardy-Littlewood-Sobolev inequality as Bedrossian-He (SIAM J. Math. Anal., 2017) or He (Nonlinearity, 2025) to obtain the bounded-ness of the density and the velocity.
Reference graph
Works this paper leans on
-
[8]
Cui S., Wang L. and Wang W. Suppression of blow-up for the 3D Patlak-Keller-Segel-Navier- Stokes system via the Couette flow. arXiv: 2412.19197
- [1]
-
[2]
Blanchet A., Carrillo J. A. and Masmoudi N. (2008). Infinite time aggregation for the critical Patlak–Keller–Segel model inR 2. Communications on Pure and Applied Mathematics, 61, 1449-1481
work page 2008
-
[3]
Blanchet A., Dolbeault J. and Perthame B. (2006). Two-dimensional Keller-Segel model: optimal critical mass and qualitative properties of the solutions. Electron. J. Differential Equations, No. 44, 32 pp
work page 2006
-
[4]
Calvez V. and Corrias L. (2008). The parabolic-parabolic Keller-Segel model inR 2. Commu- nications in Mathematical Sciences, 6(2), 417-447
work page 2008
-
[5]
Chen Q., Wei D. and Zhang Z. (2024). Transition threshold for the 3-D Couette flow in a finite channel. Memoirs of the American Mathematical Society, 296, no. 1478, v+178 pp
work page 2024
-
[6]
Cui S. and Wang W. (2024). Suppression of blow-up in multi-species Patlak-Keller-Segel- Navier-Stokes system via the Poiseuille flow in a finite channel. SIAM Journal on Mathemat- ical Analysis, 56(6) 7683-7712
work page 2024
-
[7]
Cui S., Wang L. and Wang W. (2025). Suppression of blow-up in Patlak-Keller-Segel system coupled with linearized Navier-Stokes equations via the 3D Couette flow. Journal of Differ- ential Equations, 432, 113196
work page 2025
Show all 36 references
-
[9]
and Wei J
Davila J., Pino M., Dolbeault J., Musso M. and Wei J. (2024). Existence and stability of infinite time blow-up in the Keller-Segel system. Archive for Rational Mechanics and Analy- sis, 248(4), 61
2024
-
[10]
and Wang W
Deng S., Shi B. and Wang W. (2025). Suppression of blow-up in 3-D Keller-Segel model via the Couette flow in whole space. J. Differential Equations, 432, 113265
2025
-
[11]
and Wang W
Feng Y., Shi B. and Wang W. (2022). Dissipation enhancement of planar helical flows and applications to three-dimensional Kuramoto-Sivashinsky and Keller-Segel equations. Journal of Differential Equations, 313, 420-449
2022
-
[12]
He S. (2018). Suppression of blow-up in parabolic-parabolic Patlak-Keller-Segel via strictly monotone shear flows. Nonlinearity, 31(8), 3651. SUPPRESSION OF BLOW-UP FOR PKS-NS SYSTEM 69
2018
-
[13]
He S. (2023). Enhanced dissipation and blow-up suppression in a chemotaxis-fluid system. SIAM Journal on Mathematical Analysis, 55(4), 2615-2643
2023
-
[14]
(2025) Time-dependent Flows and Their Applications in Parabolic-parabolic Patlak- Keller-Segel Systems Part I: Alternating Flows
He S. (2025) Time-dependent Flows and Their Applications in Parabolic-parabolic Patlak- Keller-Segel Systems Part I: Alternating Flows. Journal of Functional Analysis, 288(5), 110786
2025
-
[15]
He S. (2025). Time-dependent shear flows and their applications in parabolic-parabolic Patlak- Keller-Segel systems. Nonlinearity, 38(3), 035029
2025
-
[16]
and Tadmor E
He S. and Tadmor E. (2021). Multi-species patlak-keller-segel system. Indiana University Mathematics Journal, 70(4), 1577-1624
2021
-
[17]
and Painter K.J
Hillen T. and Painter K.J. (2009). A user’s guide to PDE models for chemotaxis. Journal of mathematical biology, 58(1), 183-217
2009
-
[18]
(2023) Suppression of Chemotactic Singularity via Viscous Flow with Large Buoyancy
Hu Z. (2023) Suppression of Chemotactic Singularity via Viscous Flow with Large Buoyancy. SIAM Journal on Mathematical Analysis, 56(6), 7866-7902
2023
-
[19]
and Kiselev A
Hu Z. and Kiselev A. (2023). Suppression of chemotactic blow up by strong buoyancy in stokes-boussinesq flow with cold boundary. Journal of Functional Analysis, 287(7), 110541
2023
-
[20]
and Yao Y
Hu Z., Kiselev A. and Yao Y. Suppression of chemotactic singularity by buoyancy. arXiv: 2305.01036
-
[21]
and Segel L
Keller E. and Segel L. (1970). Initiation of slime mold aggregation viewed as an instability. Journal of Theoretical Biology, 26(3), 399
1970
-
[22]
and Xu X
Kiselev A. and Xu X. (2016). Suppression of chemotactic explosion by mixing. Archive for Rational Mechanics and Analysis, 222, 1077-1112
2016
-
[23]
and Xu X
Li H., Xiang Z. and Xu X. (2025). Suppression of blow-up in Patlak-Keller-Segel-Navier-Stokes system via the Poiseuille flow. J. Differential Equations, 434, 113301
2025
-
[24]
Nagai T. (2000). Behavior of solutions to a parabolic-elliptic system modelling chemotaxis. Journal of the Korean Mathematical Society, 37, 721-732
2000
-
[25]
Patlak C. (1953). Random walk with persistence and external bias. The bulletin of mathe- matical biophysics, 15, 311-338
1953
-
[26]
Schweyer R. (2014). Stable blow-up dynamic for the parabolic-parabolic Patlak-Keller-Segel model. ariv:1403.4975
2014 arXiv
-
[27]
and Wolansky G
Shafrirand I. and Wolansky G. (2005). The logarithmic HLS inequality for systems on compact manifolds, J.Funct.Anal., 227 pp. 200-226
2005
-
[28]
and Wang W
Shi B. and Wang W. (2024). Enhanced dissipation and blow-up suppression for the three di- mensional Keller-Segel equation with the plane Couette-Poiseuille flow. Journal of Differential Equations, 403, 368-405
2024
-
[29]
and Winkler M
Souplet P. and Winkler M. (2019). Blow-up profiles for the parabolic-elliptic Keller-Segel system in dimensionsn≥3. Communications in Mathematical Physics, 367(2), 665-681
2019
-
[30]
and Winkler M
Tao Y. and Winkler M. (2016). Blow-up prevention by quadratic degradation in a two- dimensional Keller-Segel-Navier-Stokes system. Z. Angew. Math. Phys., 67, 1-23
2016
-
[31]
and Winkler M
Tao Y. and Winkler M. (2017). Critical mass for infinite-time aggregation in a chemotaxis model with indirect signal production. Journal of the European Mathematical Society, 19(12), 3641-3678
2017
-
[32]
and Zhang Y
Wang L., Wang W. and Zhang Y. Stability of a class of supercritical volume-filling chemotaxis- fluid model near Couette flow. arXiv:2410.02214. 70 SHIKUN CUI, LILI W ANG, WENDONG W ANG, AND JUNCHENG WEI
-
[33]
Wei D. (2018). Global well-posedness and blow-up for the 2-D Patlak-Keller-Segel equation. Journal of Functional Analysis, 274(2), 388-401
2018
-
[34]
and Zhang Z
Wei D. and Zhang Z. (2021). Transition threshold for the 3D Couette flow in Sobolev space. Communications on Pure and Applied Mathematics, 74(11), 2398-2479
2021
-
[35]
Winkler M. (2013). Finite-time blow-up in the higher-dimensional parabolic-parabolic Keller- Segel system. Journal de Mathematiques Pures et Appliquees, 100(5), 748-767
2013
-
[36]
and Zi R
Zeng L., Zhang Z. and Zi R. (2021). Suppression of blow-up in Patlak-Keller-Segel-Navier- Stokes system via the Couette flow. Journal of Functional Analysis, 280(10), 108-967. (Shikun Cui)School of Mathematical Sciences, Dalian University of Technology, Dalian, 116024, China E...
2021
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.