{"id":"64b8c294-ed6e-4522-afc6-0ce311a62c86","arxiv_id":"2506.11757","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the chemotaxis-Navier-Stokes system with porous medium diffusion, as m→∞ the solutions converge to a Hele-Shaw free boundary problem whose pressure satisfies the complementarity relation P∞(ΔP∞ − ∇·(χ(c∞)∇c∞))=0.","lead":"This paper proves global weak solutions for a chemotaxis-Navier-Stokes system with nonlinear porous medium diffusion for all exponents m≥3, and then shows that as m tends to infinity these solutions converge to a Hele-Shaw type free boundary problem. The result connects models of bacterial aggregation in fluids to the classical incompressible limit framework.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main-text proof of (1−n∞)∇P∞=0 uses false identity (4.37); the gap is real but a one-line Stampacchia argument repairs it, so the central claim is conditionally supported.","rationale":"The paper proves a substantial new Hele-Shaw limit for a chemotaxis-Navier-Stokes system, and the main body of uniform estimates (Lemmas 4.1–4.4) appears sound under the stated assumptions. The most direct proof-level defect is in the derivation of the third graph relation: equation (4.37) is mathematically incoherent, since it replaces the vector u∞ by ∇P∞ and then asserts a convergence that does not follow from the displayed line. This is a real gap in the proof of Theorem 2.2 as written. Nevertheless, the gap is immediately repairable: the relation (1−n∞)P∞=0 together with P∞≥0 and P∞(t)∈H^1 implies, via Stampacchia's lemma, that ∇P∞=0 almost everywhere on the set where P∞=0, and hence almost everywhere on {n∞<1}. Thus the third graph relation is not in doubt, and the complementarity relation proof in (4.38)–(4.40) goes through once this replacement is made. The sign error in equation (2.5) is another genuine defect, but it affects only the formal motivation and is not used in the rigorous convergence proof. The reader's weakest assumption, the uniform high-L^p initial bounds in H2/H3, is a limitation on the class of admissible data rather than an inconsistency: for data outside these bounds the proof's key estimates fail, but Theorem 2.2 is explicitly conditional on H2/H3. Because all identified issues are addressable without changing the theorem's statement, I agree with the conditional verdict and recommend no change to the reader's assessment.","tokens_in":32670,"tokens_out":38059,"duration_ms":336625,"concrete_test":"Analytically verify the replacement step: for a.e. t, take P∞(t)∈H^1(R^d) with P∞≥0 and (1−n∞)P∞=0; check that Stampacchia's lemma yields ∇P∞=0 a.e. on {n∞<1}, and confirm that the derivation of (4.40) from the weak form of (2.6)_1 uses only (1−n∞)P∞=0, (1−n∞)∇P∞=0, ∇·u∞=0, and P∞∈L^2(H^1). If this verification succeeds, the third graph relation and the complementarity relation (2.9) are established without (4.37).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of the third Hele-Shaw graph relation (1−n∞)∇P∞=0 in Section 4.2 relies on the displayed identity (4.37), which is incorrect as written: it mixes u∞ with ∇P∞, and the assertion that ∇P∞^{1+α} ⇀ ∇P∞ in L^{3/2}_{loc} as α→0+ does not follow from the preceding line. This is a genuine gap in the proof of Theorem 2.2 as written. However, the gap is trivial to close: from (1−n∞)P∞=0 and P∞≥0, the set {n∞<1} is contained, modulo null sets, in the zero set of P∞. Since P∞(t)∈H^1(R^d) for a.e. t and is nonnegative, Stampacchia's lemma gives ∇P∞=0 a.e. on {P∞=0}, hence a.e. on {n∞<1}. Thus (1−n∞)∇P∞=0 follows directly, and the subsequent complementarity proof (4.38)–(4.40) only needs that relation. A second, separate defect is the sign error in the pressure equation (2.5): the chemotaxis terms should carry a minus sign (∂tP+u·∇P = (m−1)P(ΔP−∇·(χ∇c)) + ∇P·(∇P−χ∇c)), otherwise the formal derivation of (2.9) in the introduction is inconsistent. This equation is not used in the rigorous proof of Theorem 2.2, so it does not by itself threaten the result. The reader's emphasis on the uniform L^{m−1}/L^{m+1} initial bounds is a scope limitation, not an internal inconsistency: without those bounds the estimates in Lemmas 3.2 and 4.3 break down, but the theorem is stated under exactly those hypotheses.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Cauchy problem for a chemotaxis-Navier-Stokes system with porous-medium diffusion, ∂t n + u·∇n = Δn^m − ∇·(χ(c)n∇c), coupled to an oxygen equation and the incompressible Navier-Stokes equations. It first proves global existence of weak solutions uniformly for m ≥ 3 under hypotheses (H1)-(H2), with estimates independent of m. It then establishes that as m → ∞ subsequences converge to a Hele-Shaw-type limit, consisting of the system (2.6), the graph relations 0 ≤ n∞ ≤ 1, (1−n∞)P∞ = 0, (1−n∞)∇P∞ = 0, and the complementarity relation P∞(ΔP∞ − ∇·(χ(c∞)∇c∞)) = 0 in the distributional sense. The proof is based on the pressure formulation P_m = (m/(m−1))n_m^{m−1}, uniform-in-m energy and compactness estimates, and a special test-function argument for the complementarity relation; Appendix A gives an alternative proof via strong compactness of ∇n_m^m.","tokens_in":33100,"tokens_out":15127,"duration_ms":152734,"significance":"If the result is accepted after revision, this is a significant contribution: it provides a uniform global-existence theorem for m ≥ 3 and the first rigorous Hele-Shaw/stiff-pressure limit for the chemotaxis-Navier-Stokes system, including fluid convection and buoyancy effects. The proof is largely self-contained and is built from explicit, parameter-free estimates; no fitted constants or assumed limiting equations are used. The alternative proof of the complementarity relation in Appendix A is a useful consistency check. However, two technical defects, one in the proof of the third graph relation and one in the displayed pressure equation, need to be repaired. Both are local and do not undermine the overall strategy, so the central claim is conditionally supported.","major_comments":[{"comment":"The identity (4.37) is false as written: it asserts u∞·∇P∞ = ∇P∞, which is not justified by the preceding arguments, and the claimed weak convergence ∇P∞^{1+α} ⇀ ∇P∞ in L^{3/2}_loc as α → 0+ does not follow from the displayed line. Since this is the only step offered for the third Hele-Shaw graph relation (1−n∞)∇P∞ = 0, the proof of Theorem 2.2 has a genuine gap at this point. The gap is easily repaired: from (1−n∞)P∞ = 0 and P∞ ≥ 0, the set {n∞ < 1} is contained, up to null sets, in {P∞ = 0}; applying Stampacchia's lemma to P∞(t) ∈ H^1(R^d) gives ∇P∞ = 0 a.e. on {P∞ = 0}, and hence (1−n∞)∇P∞ = 0. The paragraph around (4.37) should be replaced by this direct argument.","section":"Section 4.2, Eq. (4.37)"},{"comment":"Equation (2.5) has the wrong signs in the chemotaxis terms. Direct differentiation of P_m = (m/(m−1))n_m^{m−1} along (2.4) gives ∂tP_m + u_m·∇P_m = (m−1)P_m(ΔP_m − ∇·(χ(c_m)∇c_m)) + ∇P_m·(∇P_m − χ(c_m)∇c_m), not the formula with plus signs displayed in (2.5). This makes the formal derivation of the complementarity relation (2.9) inconsistent and should be corrected. The rigorous proof of Theorem 2.2 does not rely on (2.5), so this is a repairable error rather than a fatal flaw.","section":"Section 2, Eq. (2.5)"}],"minor_comments":[{"comment":"The constants in (H2) and (H3) are required to be independent of m, which forces the initial densities n_{m,0} to be essentially bounded with a uniform bound; since the L^p norm converges to the L^∞ norm as p → ∞, concentrated or genuinely unbounded initial data are excluded. This is a scope restriction of Theorem 2.2 and should be stated explicitly, for example in Remark 2.2 or immediately after the assumptions.","section":"Assumptions (H2)-(H3)"},{"comment":"The weak-* convergence statement for n_m in L∞(0,T;L^q(R^d)) is stated for q ∈ [1,∞); for q = 1, the weak-* notation is nonstandard because L^1 is not the dual of a Banach space in the usual sense. The statement should restrict q to (1,∞), or the q = 1 case should be formulated separately, for instance as weak convergence of measures.","section":"Theorem 2.2, Eq. (2.10)"},{"comment":"There are several typographical issues, including 'Texting' instead of 'Testing' in Lemma 3.3 and inconsistent typesetting of n_m^m as 'nm m'; a careful proofreading pass is recommended.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the main results are potentially valuable. The two identified issues are repairable without changing the overall framework: the proof of the third graph relation should use the Stampacchia argument, and the sign error in (2.5) should be corrected. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe headline: this paper proves a real new result and the main argument is mostly in good shape, but one displayed identity in the proof of the third Hele-Shaw graph relation is wrong as written. The error is easily repairable, so I would not reject over it.\n\nWhat is new: global-in-time weak solutions for the chemotaxis-Navier-Stokes system with porous medium diffusion for any m≥3, uniform in m, and the first rigorous Hele-Shaw limit (m→∞) for the full system with fluid convection, leading to a free boundary problem with the complementarity relation (2.9). The proof strategy follows the weak-solution framework of He, Li, Perthame [22,23,25], but the extension to the Navier-Stokes coupled system is genuinely nontrivial: uniform L^{m+1} estimates, control of ∂t n_m, and a clean test-function argument for the complementarity relation that avoids needing strong L2 compactness of ∇ n_m^m (though Appendix A supplies that route under extra hypotheses). The novelty is real.\n\nThe soft spot: in Section 4.2, equation (4.37) claims u∞·∇P∞^{1+α} = (1+α)P∞^α ∇P∞, which drops u∞ and would imply u∞·∇P∞ = ∇P∞. That is false. The subsequent conclusion (1−n∞)∇P∞=0 is therefore not justified by the text as printed. The gap is trivial to close: from (1−n∞)P∞=0 and P∞≥0, the set {n∞<1} is contained, modulo null sets, in {P∞=0}; since P∞∈H^1 and nonnegative, Stampacchia's lemma gives ∇P∞=0 a.e. on that set, which is exactly the desired relation. So the central claim is conditionally supported. There is also a sign error in the formal pressure equation (2.5) — the chemotaxis terms should carry minuses — but that equation is not used in the rigorous proof, so it is a typo-level issue, though it should be fixed because the formal derivation of (2.9) depends on it.\n\nThe uniform L^{m−1}/L^{m+1} initial bounds in (H2)–(H3) are strong: as m→∞ they force the initial density to be essentially bounded by a constant independent of m. That is a scope limitation, not an internal inconsistency — the theorem is stated under exactly those hypotheses. It should be highlighted as a limitation.\n\nOverall: this deserves a serious referee. The errors are real but repairable, and the main theorem is new and important for the PDE analysis of chemotaxis. I would recommend acceptance after a revision that fixes (4.37), or replaces it with the Stampacchia argument, and corrects the sign in (2.5).","headline":"A real new result with a repairable gap: the proof of the third Hele-Shaw graph relation uses a false identity, but a one-line Stampacchia argument fixes it.","tokens_in":33597,"tokens_out":5209,"would_cite":true,"duration_ms":45748,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35A01","35B40","35B44","35K55","76D27","92C17"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that, as the diffusion exponent m grows, weak solutions of the chemotaxis–Navier–Stokes system with porous-medium diffusion converge to a Hele-Shaw type free boundary problem in which the cell density is capped at one…","keywords":["Hele-Shaw limit","chemotaxis-Navier-Stokes system","porous medium diffusion","stiff pressure law","free boundary problem","complementarity relation","global weak solutions","degenerate elliptic equation"],"falsifier":"Take a sequence of initial data with fixed $L^1$ mass but $\\|n_{m,0}\\|_{L^{m+1}(\\mathbb{R}^d)}\\to\\infty$ as $m\\to\\infty$, so that (H3) fails, and check whether any limit still satisfies $0\\le n_\\infty\\le1$ with $P_\\infty$ supported on $\\{n_\\infty=1\\}$; if the conclusions persist, the uniform high-integrability premise is not necessary, while if they fail the premise is confirmed.","tokens_in":32487,"feed_emoji":"🦠","tokens_out":6966,"duration_ms":68109,"temperature":0.7,"pith_summary":"This paper proves two related results for the chemotaxis–Navier–Stokes system with porous-medium diffusion, a model of bacteria swimming in a viscous fluid and consuming oxygen. First, for every diffusion exponent $m\\ge3$, global weak solutions exist in $\\mathbb{R}^d$, $d\\ge2$, with regularity bounds that do not depend on $m$. Second, as $m\\to\\infty$ these solutions converge to a Hele-Shaw type free boundary problem: the limiting cell density is capped at $1$, the limiting bacterium pressure obeys the stiff pressure law, and the pressure satisfies the degenerate elliptic complementarity relation $P_\\infty(\\Delta P_\\infty-\\nabla\\cdot(\\chi(c_\\infty)\\nabla c_\\infty))=0$ in the sense of distributions. If correct, the paper gives the first rigorous derivation of the Hele-Shaw limit for a chemotaxis-fluid system and identifies the stiff pressure law as the mechanism that prevents the cell density from exceeding one.","feed_headline":"Bacteria-fluid flows harden into Hele-Shaw fronts","feed_subtitle":"Proof that as diffusion stiffens, cell density saturates at 1 and pressure solves a degenerate elliptic equation.","key_machinery":"The argument is carried by the effective bacterium pressure $P_m=\\frac{m}{m-1}n_m^{m-1}$, whose equation is close to a porous-medium equation with chemotactic drift, together with the energy functional $E(t)=\\int_{\\mathbb{R}^d}(\\frac{1}{m-2}P_m+\\frac12|\\nabla c_m|^2+\\frac12|u_m|^2)\\,dx$. Uniform-in-$m$ estimates from this energy, an $L^{m+1}$ estimate for $n_m$, the bound $m\\|(n_m-1)_+\\|^2_{L^2}\\le C$ obtained through a Newtonian-potential test function, and the identity $n_m^m=\\frac{m-1}{m}n_mP_m$ identify both weak limits as $P_\\infty$. The complementarity relation is then verified without strong compactness of gradients, by testing the limiting system with $\\varphi P_\\infty$ and using a difference-quotient argument to show $\\int\\partial_t n_\\infty\\,\\varphi P_\\infty=0$.","core_discovery":"Theorem 2.2 is the central claim: under assumptions (H1)–(H3), for $m\\ge\\max\\{2d+1,5\\}$, the weak solutions $(n_m,c_m,u_m)$ constructed in Theorem 2.1 converge, up to subsequences, to a limit $(n_\\infty,c_\\infty,u_\\infty,P_\\infty,\\Pi_\\infty)$. The pressure $P_m=\\frac{m}{m-1}n_m^{m-1}$ and the scaled density $n_m^m/m$ have the same weak limit $P_\\infty$ in $L^2(0,T;H^1)$, and the limit satisfies the Hele-Shaw type system (2.6)–(2.7) together with the graph relations $0\\le n_\\infty\\le1$, $(1-n_\\infty)P_\\infty=0$, and $(1-n_\\infty)\\nabla P_\\infty=0$ almost everywhere. The complementarity relation $P_\\infty(\\Delta P_\\infty-\\nabla\\cdot(\\chi(c_\\infty)\\nabla c_\\infty))=0$ holds distributionally; in the saturation region $\\{n_\\infty=1\\}$, where the pressure is supported, this is a degenerate elliptic equation for the limiting pressure.","pith_inferences":["Not spelled out in the paper, the estimate $m\\|(n_m-1)_+\\|^2_{L^2}\\le C$ implies an $O(m^{-1/2})$ $L^2$ decay rate for the positive part of $n_m-1$, which could be checked numerically.","Because the special test-function route to the complementarity relation bypasses the Aronson–Bénilan estimates and strong-gradient compactness used in earlier Keller–Segel Hele-Shaw proofs, the same mechanism may verify the complementarity relation in settings where strong compactness of $\\nabla P_m$ is unavailable.","The proof uses only the porous-medium-type structure of the density equation, so it likely extends to chemotaxis-fluid systems with volume-filling effects, logistic growth, or more general chemotactic sensitivities; a concrete test would be the same limit with the reaction term $nf(c)$ replaced by $nf(c)+g(n)$.","The uniform bounds in (H3) grow with $m$, so the theorem as stated covers only essentially bounded initial cell densities; whether the Hele-Shaw limit holds for unbounded or concentrated data is a natural open extension."],"forward_implications":["Global weak solutions exist for the full chemotaxis–Navier–Stokes system in any dimension $d\\ge2$ for every $m\\ge3$, without the structural conditions (1.6) or (1.8) and without the spatial-weight assumption on $n_0$.","The Hele-Shaw limit is justified: as $m\\to\\infty$, $u_m\\to u_\\infty$ strongly in $L^2(0,T;L^2_{\\mathrm{loc}})$ and $c_m\\to c_\\infty$ strongly in $L^2(0,T;W^{1,p}_{\\mathrm{loc}})$, while the pressure converges weakly and $n_m$ converges in $\\dot{H}^{-1}_{\\mathrm{loc}}$.","The limiting cell density is saturated: $0\\le n_\\infty\\le1$ almost everywhere, and the pressure $P_\\infty$ is supported exactly on the set where $n_\\infty=1$.","In the saturation region the pressure solves the degenerate elliptic equation $\\Delta P_\\infty=\\nabla\\cdot(\\chi(c_\\infty)\\nabla c_\\infty)$ where $P_\\infty>0$.","If the initial cell mass is finite, the saturation region is bounded and the pressure is compactly supported, so the limit carries a genuine free boundary."],"supporting_citations":[{"why":"Introduces the chemotaxis–Navier–Stokes system (1.1) as the model of bacteria swimming in a viscous fluid.","marker":"[48]"},{"why":"Establishes the Hele-Shaw asymptotics for porous-medium tumor growth models, the template for the stiff pressure limit.","marker":"[42]"},{"why":"Gives the first Hele-Shaw limit for the Patlak-Keller-Segel model, the baseline that this paper extends to a fluid-coupled system.","marker":"[5]"},{"why":"Provides the weak-solution framework for the incompressible limit of Patlak-Keller-Segel with complementarity relation.","marker":"[23]"},{"why":"Develops the porous-medium-with-chemotaxis incompressible limit whose compactness method Appendix A follows.","marker":"[22]"},{"why":"Supplies the special test-function and difference-quotient device used here to verify the complementarity relation.","marker":"[6]"},{"why":"Prior global existence result for coupled chemotaxis-fluid equations with nonlinear diffusion whose conditions Theorem 2.1 relaxes.","marker":"[34]"}],"fun_headline_variants":["Chemotaxis flows stiffen into Hele-Shaw limit","Bacterial pressure solves degenerate elliptic in Hele-Shaw limit","As diffusion hardens, cell flow converges to Hele-Shaw","Hele-Shaw limit proven for chemotaxis-Navier-Stokes","Stiff pressure law emerges in cell-flow Hele-Shaw limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The uniform-in-$m$ bounds on the initial cell density in $L^{m-1}$ and $L^{m+1}$ (assumptions (H2) and (H3)) are load-bearing, since these norms grow with $m$ and force the initial density to be essentially bounded by a constant independent of $m$, excluding genuinely unbounded or strongly concentrated initial data.","fun_headline_variants_meta":{"raw":{"variants":["Chemotaxis flows stiffen into Hele-Shaw limit","Bacterial pressure solves degenerate elliptic in Hele-Shaw limit","As diffusion hardens, cell flow converges to Hele-Shaw","Hele-Shaw limit proven for chemotaxis-Navier-Stokes","Stiff pressure law emerges in cell-flow Hele-Shaw limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000723,"raw_usage":{"total_tokens":3247,"prompt_tokens":951,"completion_tokens":2296,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":2216}},"tokens_in":567,"tokens_out":2296,"duration_ms":17589,"temperature":1.0,"reasoning_tokens":2216,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:07:18.438142+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a sequence of initial data with fixed $L^1$ mass but $\\|n_{m,0}\\|_{L^{m+1}(\\mathbb{R}^d)}\\to\\infty$ as $m\\to\\infty$, so that (H3) fails, and check whether any limit still satisfies $0\\le n_\\infty\\le1$ with $P_\\infty$ supported on $\\{n_\\infty=1\\}$; if the conclusions persist, the uniform high-integrability premise is not necessary, while if they fail the premise is confirmed.","supporting_citations":[{"cited_title":"Tuval, L","cited_arxiv_id":null,"evidence_quote":"Introduces the chemotaxis–Navier–Stokes system (1.1) as the model of bacteria swimming in a viscous fluid."},{"cited_title":"Perthame, F","cited_arxiv_id":null,"evidence_quote":"Establishes the Hele-Shaw asymptotics for porous-medium tumor growth models, the template for the stiff pressure limit."},{"cited_title":"Craig, I","cited_arxiv_id":null,"evidence_quote":"Gives the first Hele-Shaw limit for the Patlak-Keller-Segel model, the baseline that this paper extends to a fluid-coupled system."},{"cited_title":"He, H.-L","cited_arxiv_id":null,"evidence_quote":"Provides the weak-solution framework for the incompressible limit of Patlak-Keller-Segel with complementarity relation."},{"cited_title":"Incompressible limit of porous media equation with chemotaxis and growth","cited_arxiv_id":"2312.16869","evidence_quote":"Develops the porous-medium-with-chemotaxis incompressible limit whose compactness method Appendix A follows."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the special test-function and difference-quotient device used here to verify the complementarity relation."},{"cited_title":"Liu and A","cited_arxiv_id":null,"evidence_quote":"Prior global existence result for coupled chemotaxis-fluid equations with nonlinear diffusion whose conditions Theorem 2.1 relaxes."}],"review_version":1}