{"id":"311ac817-1e20-4b38-8908-513ef76d64f0","arxiv_id":"2506.10578","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For the 3D Patlak-Keller-Segel-Navier-Stokes system near a strong Couette flow, global regularity is proved for initial cell mass below 16π².","lead":"A new proof shows that in a 3D chemotaxis-fluid model, a strong Couette shear flow can prevent blow-up of the cell density whenever the total initial mass is below 16π², matching the known 2D Keller-Segel critical mass heuristic. The paper is a technical PDE stability analysis that could interest applied mathematicians working on mixing, enhanced dissipation, and chemotaxis.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sharpness claim rests on an invalid u=0 reduction: in the forced system (1.3), u=0 with x-independent n is not a solution because the buoyancy term drives u1; no blow-up above 16π² is proved.","rationale":"The reader's weakest assumption, condition (1.6), is indeed restrictive and load-bearing for the zero-mode estimates, but it is an explicit hypothesis of Theorem 1.1 rather than an invalid step. My concern targets the separate 'sharp threshold' assertion in the title and abstract: the paper's heuristic reduction to 2D Keller-Segel via u=0 is not a valid trajectory of the forced PKS-NS system because the buoyancy term n e_x cannot be balanced by a periodic pressure when n is positive and x-independent. This does not undermine the global-existence proof below 16π², so the reader's CONDITIONAL verdict remains appropriate. I would not change the verdict, but the sharpness claim should be either removed or backed by a genuine blow-up theorem for M>16π².","tokens_in":74019,"tokens_out":34436,"duration_ms":403848,"concrete_test":"Substitute the proposed reduction u=0, n=n0(t,y,z), c=c0(t,y,z) into (1.3). The u1 equation becomes ∂t u1 + ∂x P = n0(t,y,z). Taking the x-average of this identity gives ∂t (u1)^0 = n0(t,y,z), which is nonzero for positive n0; hence the velocity cannot stay zero. This settles that the heuristic is not an exact solution. To test sharpness directly, one would need a rigorous blow-up construction for M>16π² or a positive result that no global solution exists for some supercritical data; neither is supplied.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The main bootstrap for Theorem 1.1 is extensive and I did not find a concrete error in its core estimates. The load-bearing weakness is the advertised 'sharp' threshold in the title and abstract. Remark 1.1 justifies sharpness by saying that when u vanishes and n=n(t,y,z), equation (1.8) reduces to the 2D Keller-Segel equation with critical mass 8π, hence M<16π². This reduction is not a solution of the perturbation system. In (1.3)/(1.5), the u1 equation contains the explicit buoyancy forcing n (after rescaling, n/A). Setting u=0 makes the u1 equation read ∂t u1 + ∂x P = n(t,y,z); for x-independent positive n, no periodic pressure can balance n e_x, so u cannot remain zero. The zero-mode equation in Section 2.2 gives the same obstruction: at u=0, ∂t u1,0 = n0/A, forcing immediate growth of u1. Thus the exact 2D Keller-Segel threshold is not a rigorous consequence of the reduced zero-mode equation, and the paper supplies no other finite-time blow-up argument for M>16π². Theorem 1.1 remains a credible global-existence statement below the threshold, but the word 'sharp' is unsupported by the provided argument.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":74308,"tokens_out":9789,"duration_ms":109924,"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":[{"comment":"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.","section":"§1, Remark 1.1, and the abstract"},{"comment":"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.","section":"Lemmas 3.6, A.2, A.3, Proposition A.2, and Lemma 8.1"},{"comment":"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.","section":"Theorem 1.1, condition (1.6), and the abstract"}],"minor_comments":[{"comment":"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.","section":"Abstract and Remark 1.1"},{"comment":"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.","section":"§2.1 and Notations"},{"comment":"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.","section":"Appendix A.1"},{"comment":"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.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The main global-existence theorem appears credible based on the overall structure and the absence of an evident concrete error in the bootstrap estimates. However, the advertised sharpness of the threshold is not supported by any rigorous argument, and the load-bearing reliance on the unpublished companion preprint [8] makes verification difficult. The paper would be acceptable after the sharpness claim is removed or proved, the omitted proofs or precise statements are supplied, and the abstract is aligned with the actual hypotheses of Theorem 1.1. I recommend major revision rather than rejection because the global-existence result itself is a meaningful advance over the previously known thresholds."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main theorem is a real step forward: global existence for the 3D PKS-NS system near Couette flow when the initial cell mass is below 16π², improving the earlier 8π/9 and 24/5π² thresholds. The proof is a long bootstrap, and I did not find a concrete error in the core estimates. The genuinely new piece is the dissipative decay for the transverse velocity zero modes in Lemma 4.3, which lets the authors apply the logarithmic Hardy-Littlewood-Sobolev inequality to control the density zero mode in L² uniformly in time. That part is coherent.\n\nThe main problem is the word 'sharp' in the title. The paper proves no blow-up for masses above 16π². The heuristic in Remark 1.1—that with u=0 and x-independent density the zero-mode equation reduces to 2D Keller-Segel with critical mass 8π—does not hold as a statement about the full system. If u=0 and n=n(t,y,z)>0, the u1 equation contains the buoyancy forcing n, which cannot be balanced by a periodic pressure; u=0 is not a solution. The reduction is an analogy, not a rigorous argument. At best the threshold is conjecturally optimal. This needs to be fixed in the title and abstract.\n\nSecond, several lemmas are imported from the authors' own preprint [8], including the Fourier-relationship Lemma 3.1, the interaction Lemma 3.6, and parts of the space-time estimates in Appendix A. Some are standard, but a referee cannot verify the proof without chasing an unpublished arXiv paper. The authors should either include full proofs or cite a published version of [8] and confirm the statements match.\n\nThird, the smallness condition on (u2,in)^0 and (u3,in)^0 is a separate restriction on the data; large A does not produce it. The paper is clear about this, but it narrows the class of flows covered.\n\nI did not see evidence of circular reasoning: the mass condition enters through the log-HLS bound, and the bootstrap constants are chosen independently of the threshold.\n\nThis is a paper for specialists in chemotaxis-fluid equations and hydrodynamic stability. It deserves a serious referee. My recommendation is to send it to peer review, with the clear expectation of major revision: prove or explicitly soften the sharpness claim, and make the deferred lemmas available.","headline":"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.","tokens_in":74872,"tokens_out":4511,"would_cite":true,"duration_ms":51322,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q35","35Q92","35B44","76D05"],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["suppression of blow-up","Patlak-Keller-Segel-Navier-Stokes system","Couette flow","critical mass threshold","16 pi^2","enhanced dissipation","zero mode","chemotaxis-fluid coupling"],"falsifier":"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.","tokens_in":73797,"feed_emoji":"🦠","tokens_out":10107,"duration_ms":106633,"temperature":0.7,"pith_summary":"Alone, the three-dimensional Patlak–Keller–Segel system may blow up in finite time for arbitrarily small initial cell mass, so any threshold for global regularity has to come from additional mechanisms. This paper proves that a sufficiently strong Couette flow $(Ay,0,0)$ provides such a mechanism: if the initial $x$-averages of the two transverse velocity components are sufficiently small and the total initial cell mass satisfies $M < 16\\pi^{2}$, the coupled chemotaxis–Navier–Stokes system around the Couette flow has a global-in-time solution on $\\mathbb{T}^{3}$. The paper argues that $16\\pi^{2}$ is the sharp threshold, because the $x$-averaged density obeys an equation that reduces to the two-dimensional Keller–Segel equation, whose critical mass is $8\\pi$; with the torus normalization this is exactly $M = 2\\pi \\cdot 8\\pi = 16\\pi^{2}$. A new dissipative decay estimate for the non-constant part of the transverse velocity zero modes is what keeps the density below the critical value.","feed_headline":"Couette flow tames 3D chemotaxis blow-up below mass 16π²","feed_subtitle":"Global regularity for the fluid-coupled chemotaxis system when total cell mass stays under the 16π² threshold.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the zero-mode density estimate via the logarithmic Hardy–Littlewood–Sobolev inequality for shear flows, the template for the present blow-up suppression argument.","marker":"[1]"},{"why":"Provides the quasi-linear method and the space-time estimates (including the coordinate-transform propositions) used to control nonlinear interactions and the 3D lift-up effect.","marker":"[34]"},{"why":"Establishes the 2D Keller–Segel critical mass $8\\pi$ for $M < 8\\pi$, which is the basis for the heuristic that $16\\pi^{2}$ is sharp.","marker":"[3]"},{"why":"Gives the 2D dichotomy global well-posedness iff $M \\leq 8\\pi$, underpinning the sharpness statement for the reduced zero-mode equation.","marker":"[33]"},{"why":"States the logarithmic Hardy–Littlewood–Sobolev inequality used in Lemma A.6 to bound $\\|n_{0}\\log n_{0}\\|_{L^{1}}$ below the critical mass.","marker":"[27]"},{"why":"Previous threshold $\\tfrac{24}{5}\\pi^{2}$ obtained by the same authors in the same Couette setting; the present result improves it to $16\\pi^{2}$.","marker":"[8]"}],"fun_headline_variants":["Strong Couette flow halts 3D chemotaxis blow-up below 16π²","Sharp 16π² mass threshold for global 3D chemotaxis via Couette flow","Global regularity for chemotaxis with Couette flow and mass < 16π²","Couette flow suppresses 3D chemotaxis blow-up for mass under 16π²","Mass < 16π² and strong Couette flow give global chemotaxis solutions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Strong Couette flow halts 3D chemotaxis blow-up below 16π²","Sharp 16π² mass threshold for global 3D chemotaxis via Couette flow","Global regularity for chemotaxis with Couette flow and mass < 16π²","Couette flow suppresses 3D chemotaxis blow-up for mass under 16π²","Mass < 16π² and strong Couette flow give global chemotaxis solutions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000674,"raw_usage":{"total_tokens":3181,"prompt_tokens":1173,"completion_tokens":2008,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":789,"completion_tokens_details":{"reasoning_tokens":1891}},"tokens_in":789,"tokens_out":2008,"duration_ms":18162,"temperature":1.0,"reasoning_tokens":1891,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:24:06.101822+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"and He S","cited_arxiv_id":null,"evidence_quote":"Supplies the zero-mode density estimate via the logarithmic Hardy–Littlewood–Sobolev inequality for shear flows, the template for the present blow-up suppression argument."},{"cited_title":"and Zhang Z","cited_arxiv_id":null,"evidence_quote":"Provides the quasi-linear method and the space-time estimates (including the coordinate-transform propositions) used to control nonlinear interactions and the 3D lift-up effect."},{"cited_title":"and Perthame B","cited_arxiv_id":null,"evidence_quote":"Establishes the 2D Keller–Segel critical mass $8\\pi$ for $M < 8\\pi$, which is the basis for the heuristic that $16\\pi^{2}$ is sharp."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the 2D dichotomy global well-posedness iff $M \\leq 8\\pi$, underpinning the sharpness statement for the reduced zero-mode equation."},{"cited_title":"and Wolansky G","cited_arxiv_id":null,"evidence_quote":"States the logarithmic Hardy–Littlewood–Sobolev inequality used in Lemma A.6 to bound $\\|n_{0}\\log n_{0}\\|_{L^{1}}$ below the critical mass."}],"review_version":1}