{"id":"5995ef4c-d1df-408c-b56d-e34971d5bb08","arxiv_id":"2505.23545","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Global well-posedness, extinction criteria, and, for the quasi-steady system, unique stable equilibria are proved rigorously for a one-dimensional moving-boundary biofilm model.","lead":"This mathematics paper proves that a simplified one-dimensional model of biofilm growth has a unique well-behaved solution at all times, and describes when the biofilm shrinks to nothing or settles at a stable nonzero thickness. It is read because it offers rigorous guarantees and long-time behavior predictions for a moving-boundary model that biofilm researchers use but had not analyzed.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Step 4 of Theorem 2.1 asserts global existence from a linear bound without the required regularizing estimate; the conclusion needs an explicit variation-of-constants/Grönwall argument.","rationale":"The reader's weakest_assumption flags both the import of Amann's theory and the Step 4 continuation argument. I agree with the second as the genuinely load-bearing gap. The Amann import in Steps 1–3 is standard: the space-verification sketches in (2.3)–(2.6) are adequate, and the chosen interpolation range makes the abstract spaces independent of the boundary operator, so I do not see a concrete failure there. Step 4, however, contains an actual omitted argument: linear growth of the semilinear term in the weaker norm E_γ does not by itself give the E_α bound required by the continuation criterion. In classical quasilinear parabolic theory this is supplied by a regularizing estimate for the evolution family, and the paper neither derives nor cites it. Since this is a standard tool, the gap is fillable rather than fatal; the proper disposition remains CONDITIONAL. The quasi-steady section is largely self-contained and correct; the only blemishes I found are sign inconsistencies in the comparison argument of Lemma 3.4 that appear to be typos and do not affect the stated upper bound. Thus the verdict should stay unchanged.","tokens_in":18490,"tokens_out":25754,"duration_ms":271405,"concrete_test":"Complete the missing argument: for T<t^+(w0), let U(t,s) be the evolution family of \\hat A(t); write v(t)=U(t,0)v0+∫_0^t U(t,s)\\hat f(s,v(s))ds, apply the regularizing estimate ||U(t,s)||_{L(E_γ,E_α)} ≤ C(t-s)^{-θ} with θ=α-γ, combine with (2.17), and derive a Grönwall inequality for ||v(t)||_{E_α}. If this derivation succeeds, Step 4 of Theorem 2.1 is valid; if it yields only a bound in a strictly weaker norm, global existence is not established. Independently, check whether the needed regularizing estimate is explicitly stated in Amann [1] or in Lunardi's monograph; a precise citation would also resolve the omission.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Step 4 of the proof of Theorem 2.1 (Section 2.1) is the load-bearing point. After deriving 0≤v≤c* and h∈BUC^{1-}, the proof defines the semilinear problem v_t = \\hat A(t)v + \\hat f(t,v) and states the linear bound (2.17) in E_γ = L_p: ||\\hat f(t,v)||_{E_γ} ≤ m0(1+||v||_{E_α}). It then concludes without further argument that ||v(t)||_{W^1_p} = ||v(t)||_{E_α} is bounded on the whole maximal interval. That conclusion does not follow from the two displayed inequalities alone. The continuation criterion (2.9) requires control of ||v(t)||_{E_α}; (2.17) only controls the forcing in the strictly weaker space E_γ. The missing ingredient is a regularizing estimate for the evolution system generated by \\hat A(t), e.g. ||U(t,s)||_{L(E_γ,E_α)} ≤ C(t-s)^{-(α-γ)}, followed by a Grönwall argument. No such estimate is stated or cited. If the estimate cannot be supplied, Theorem 2.1's global existence assertion is not proven, and every later evolutionary result (Propositions 2.2 and 2.3) inherits this gap. The gap is fixable by standard theory, but as written it is an unproved assertion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a one-dimensional moving-boundary biofilm model (1.1) and its reduced forms. After a fixed-domain change of variables, the evolution problem (2.2) is treated as a quasilinear parabolic system in Sobolev spaces; Theorem 2.1 claims unique global strong well-posedness. The paper also studies the quasi-steady system (3.2): existence and uniqueness of solutions, equilibria for general nonlinearities, and for the affine growth rate g(s)=alpha(s-b), uniqueness of the nontrivial equilibrium and convergence to it when r(c*)>b, or extinction when r(c*)<=b.","tokens_in":18574,"tokens_out":20513,"duration_ms":214508,"significance":"If the main theorem is fully proved, the paper would provide a rigorous well-posedness and long-time behavior framework for a widely used biofilm model, with the quasi-steady reduction giving a fairly complete qualitative picture. Strengths include the self-contained fixed-point and shooting arguments, the use of standard abstract quasilinear theory rather than ad hoc assumptions, and the candid discussion of the open stability question in Remark 3.11. The claims involve no fitted parameters and are concretely falsifiable. However, the global-existence proof in Theorem 2.1 contains an unproved smoothing step, so the central claim is not yet established as written.","major_comments":[{"comment":"The global-existence proof contains a load-bearing gap. After deriving 0<=v<=c* and h in BUC^{1-}, the proof only obtains the linear bound ||hat f(t,v)||_{E_gamma} <= m0(1+||v||_{E_alpha}) in (2.17), with E_gamma=L_p and E_alpha=W^1_p, and then asserts that ||v(t)||_{W^1_p}=||v(t)||_{E_alpha}<=C on the whole maximal interval. The continuation criterion (2.9) requires a bound in E_alpha, but (2.17) is not a differential inequality for ||v||_{E_alpha}. One needs a regularizing estimate for the evolution family generated by hat A(t), such as ||U(t,s)||_{L(E_gamma,E_alpha)}<=C(t-s)^{-(alpha-gamma)}, followed by a Gronwall argument. No such estimate is stated or cited. Since Theorem 2.1 and the subsequent evolutionary results depend on global existence, this step must be repaired; it is likely fixable by standard theory, but as written it is an unproved assertion.","section":"Section 2.1, Step 4 (after (2.17))"},{"comment":"The comparison argument for the upper bound u=c*-v is not stated correctly. From (2.2c) one obtains the boundary condition kappa_L/(L h(t)) u(t,1)+kappa/h(t)^2 u_y(t,1)=kappa_L c*/(L h(t)), not B(t)u(t)=c*, if B(t) is the operator defined in Step 1. If a normalized boundary operator is intended, this should be defined explicitly; otherwise the invocation of [1, Theorem 15.1] with boundary value c* is not literally justified. The positivity conclusion is plausible also with the positive time-dependent boundary value, but the reduction needs to be made explicit.","section":"Section 2.1, Step 3"}],"minor_comments":[{"comment":"The inference that h(t) tends to 0 needs the uniform negativity bound g(r(s))<=-delta<0 on [0,c*], which follows from the assumptions and compactness; without stating it, the exponential decay of h is not immediate from h' = h G(v(t)).","section":"Section 2.2, Proposition 2.2"},{"comment":"The sentence about the monotonicity of h is terse; please spell out why h is monotone, for instance by using the constant sign of f on each side of h_e, before taking the limit h_*.","section":"Section 3.2.2, Proposition 3.10"},{"comment":"There are a few minor typos (for example 'susbtrate' in the introduction), and the equivalence between (2.2c) and the original Robin condition (1.3) should be shown explicitly because the normalization of B is used repeatedly.","section":"Sections 1 and 2"},{"comment":"Please check the constants and missing factors in the displayed boundary values, in particular the expression for d(t,h(t)) in Proposition 2.3 and the boundary value for w in the proof of Lemma 3.4; as printed the scaling is not consistent.","section":"Proposition 2.3 and Lemma 3.4"}],"recommendation":"major_revision","confidential_remarks":"The decisive issue is the missing regularizing estimate in Step 4 of Theorem 2.1. If the authors supply that estimate by a standard variation-of-constants/Gronwall argument or by citing the appropriate theorem from Amann's paper, the manuscript is likely acceptable. The remaining comments are local and should not require substantial reworking. The paper is a good fit for math.AP."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a competent paper that proves genuinely new results for a specific reduced 1-D biofilm model. The main theorems are Theorem 2.1 (global strong well-posedness) and the quasi-steady analysis (equilibria, uniqueness, attraction/extinction). The methods are standard—Amann's quasilinear theory, fixed-point and shooting arguments—but applied carefully and with clear statements. The paper is honest about its limitations, including that stability of the nonzero equilibrium for the full time-dependent problem is open.\n\nThe proofs I checked line up: the fixed-domain transform, the elliptic subproblem in Proposition 3.2, the comparison and decay arguments, and the shooting argument for unique equilibrium are coherent. The reader's conditional verdict matches my reading. The one real soft spot is Step 4 of Theorem 2.1. After deriving the linear bound (2.17) in L_p, the authors assert without further argument that ||v(t)||_{W^1_p} stays bounded on the maximal interval. That step needs a regularizing estimate for the evolutionary system generated by \\hat A(t), e.g., ||U(t,s)||_{L(E_gamma,E_alpha)} <= C(t-s)^{-(alpha-gamma)}, plus a Gronwall argument. The paper neither states nor cites such an estimate. I think it is fixable with standard theory, but as written it is a missing argument, and the global existence theorem leans on it. Same issue propagates into Propositions 2.2 and 2.3.\n\nMinor things: the numerics in Remark 3.11 are illustrative without code, which is fine for a theorem paper but means the simulation-based suggestion is uncheckable. The model is a toy, but the authors say so and the analysis is rigorous.\n\nWho is this for? People working on free-boundary biofilm models or on quasilinear parabolic theory with a moving boundary can use the theorems and the techniques. I would not desk-reject it; it deserves a serious referee. I'd recommend sending it to review with a request to fill the Step 4 gap, either by adding the regularizing estimate or citing a theorem that directly implies it.","headline":"Solid, careful analysis of a toy biofilm model; new well-posedness and long-time results, but Theorem 2.1's global-existence proof skips a needed regularizing estimate.","tokens_in":19308,"tokens_out":1873,"would_cite":true,"duration_ms":20403,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q92","35R35"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a one-dimensional biofilm model with a moving boundary, the paper proves unique global strong solutions exist and that the reduced dynamics converge either to extinction or to a unique nontrivial equilibrium.","keywords":["biofilm","free boundary problem","moving boundary","global well-posedness","quasi-steady approximation","equilibria","parabolic PDE","Sobolev spaces"],"falsifier":"Solve (2.2) numerically with $p=2$, a concrete admissible pair such as $r(s)=s$, $g(s)=s-b$ with $r(c_*)>b$, and initial data $v_0$ in $[0,c_*]$ with large $W^1_2$ norm. Theorem 2.1 predicts that $h(t)$ remains positive and finite and that $\\|v(t)\\|_{W^1_2}$ stays bounded on every finite time interval; any computed finite-time blow-up or $h(t)$ hitting zero before that time would refute the global-existence claim. The decisive point to inspect is the Step 4 argument: without the asserted non-autonomous parabolic smoothing estimate, the uniform bound on the maximal interval does not follow.","tokens_in":18082,"feed_emoji":"🦠","tokens_out":8727,"duration_ms":83569,"temperature":0.7,"pith_summary":"The paper studies a reduced one-dimensional model of a biofilm growing in a diffusive boundary layer, where the biofilm height is a moving boundary driven by substrate consumption. It proves that this moving-boundary parabolic system is globally well-posed in strong Sobolev spaces: for every initial height and every initial substrate profile in $W^1_p$ taking values between $0$ and the bulk concentration $c_*$, there is a unique global solution, and the substrate stays within the physical range throughout the evolution. For the quasi-steady approximation, the paper shows that equilibria exist when the growth function changes sign between starvation and saturated conditions, and that for a linear growth law a unique nonzero equilibrium attracts all solutions when the saturated consumption rate exceeds the maintenance level, with extinction otherwise. These results put a widely used toy model of biofilm dynamics on a rigorous footing and exhibit the full range of long-time behaviour the model can produce.","feed_headline":"Global solutions proven for moving-boundary biofilm model","feed_subtitle":"Full evolution has unique strong solutions; quasi-steady dynamics end in extinction or a unique nontrivial equilibrium.","key_machinery":"The load-bearing object is the dimensionless fixed-domain transformation (2.1) mapping the moving interface $z=h(t)$ to the fixed interval $y\\in(0,1)$, converting the free-boundary problem into a quasilinear parabolic initial-boundary-value problem (2.2) for the pair $(v,h)$. Existence of a unique maximal solution follows from an abstract quasilinear parabolic framework; comparison arguments force the substrate deficit $v$ to stay in $[0,c_*]$, and a closing uniform bound on $\\|v(t)\\|_{W^1_p}$ rules out finite-time blow-up, yielding global existence. For the quasi-steady approximation, the key object is the elliptic profile $u[h]$, the unique solution of a fixed-$h$ boundary-value problem; a Schäfer fixed-point argument gives existence and uniqueness of $u[h]$, and a shooting argument using the monotone solution family of Lemma 3.7 selects the unique equilibrium height for affine growth.","core_discovery":"The central result is Theorem 2.1: after the change of variables $y=z/h(t)$, $v=c_*-c$, the moving-boundary problem becomes a fixed-domain quasilinear parabolic system for the pair $(v,h)$, and the paper proves—by importing an abstract quasilinear parabolic theory—that for every $p\\in(1,\\infty)$ and every initial $(h_0,v_0)$ with $0\\le v_0\\le c_*$ there is a unique global strong solution with the stated regularity, preserving the invariant interval $0\\le v\\le c_*$. The companion quasi-steady results show that the equilibrium set is governed by the sign of $g(0)$ and $g(r(c_*))$: for monotone $g$, equilibria exist exactly when the growth function changes sign, and for the affine choice $g(s)=\\alpha(s-b)$ the nonzero equilibrium is unique and globally attracting when $r(c_*)>b$, while extinction is the only long-time outcome when $r(c_*)\\le b$.","pith_inferences":["Editorial extension: the same fixed-domain transformation plus an elliptic-profile reduction should apply to other one-dimensional moving-boundary models with a diffusing nutrient, so the paper's method is a template as much as a result.","Editorial extension: the paper notes in Remark 3.11 that stability of the nontrivial equilibrium for the full evolutionary system remains open, since the linearization of the height equation vanishes; a center-manifold or higher-order argument would be needed, and the observed oscillatory convergence suggests the attraction is not exponential.","Editorial extension: the quasi-steady approximation is justified formally by taking the diffusion time scales to zero; a rigorous singular-limit proof connecting Theorem 2.1 to Propositions 3.9 and 3.10 would test whether the simplified dynamics faithfully represent the full model."],"forward_implications":["For every admissible initial configuration, the full moving-boundary system has a unique global strong solution, and the substrate concentration remains within the physical range $[0,c_*]$ for all time.","If the biofilm growth rate is always negative, the biofilm height decays monotonically to zero and the substrate concentration becomes spatially uniform at the bulk level $c_*$.","If the initial substrate profile is monotone in space, it stays monotone, and the substrate gradient is controlled by the initial gradient and the ratio of diffusivities.","In the quasi-steady approximation, equilibria exist whenever the growth function is negative at zero substrate and positive at saturated substrate; with an increasing growth function that never becomes positive, extinction is the only outcome, with the profile converging to $c_*$.","For affine growth $g(s)=\\alpha(s-b)$, the quasi-steady dynamics are sharp: a unique nonzero equilibrium exists and attracts every trajectory when $r(c_*)>b$, while extinction holds when $r(c_*)\\le b$."],"supporting_citations":[{"why":"Supplies the abstract quasilinear parabolic well-posedness, maximal continuation, comparison, and regularity theory on which Theorem 2.1 and the global-existence argument rest.","marker":"[1]"},{"why":"Introduces the one-dimensional biofilm toy model and its moving-boundary formulation that the paper analyzes.","marker":"[4]"},{"why":"Derives the kinetic condition (1.1f) for the biofilm height from incompressibility and Darcy's law.","marker":"[2]"},{"why":"Provides the Schäfer fixed-point theorem used to construct the elliptic profile u[h] in Proposition 3.2.","marker":"[3]"}],"fun_headline_variants":["Moving-boundary biofilm: global strong well-posedness proven","Biofilm model: extinction or unique stable equilibrium","Parabolic free boundary problem for biofilm: global solvability","Biofilm equilibria: uniqueness and stability in special case","One-dimensional biofilm: global solutions and equilibrium analysis"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the abstract quasilinear parabolic theory imported in Step 1 of the proof of Theorem 2.1 applies with all hypotheses satisfied, and that the Step 4 closing estimate—that a linear growth bound forces the solution's gradient to stay uniformly bounded on the whole maximal time interval—is valid; if either ingredient fails, the global-existence theorem collapses.","fun_headline_variants_meta":{"raw":{"variants":["Moving-boundary biofilm: global strong well-posedness proven","Biofilm model: extinction or unique stable equilibrium","Parabolic free boundary problem for biofilm: global solvability","Biofilm equilibria: uniqueness and stability in special case","One-dimensional biofilm: global solutions and equilibrium analysis"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000598,"raw_usage":{"total_tokens":2740,"prompt_tokens":832,"completion_tokens":1908,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":448,"completion_tokens_details":{"reasoning_tokens":1829}},"tokens_in":448,"tokens_out":1908,"duration_ms":15205,"temperature":1.0,"reasoning_tokens":1829,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:47:50.163543+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve (2.2) numerically with $p=2$, a concrete admissible pair such as $r(s)=s$, $g(s)=s-b$ with $r(c_*)>b$, and initial data $v_0$ in $[0,c_*]$ with large $W^1_2$ norm. Theorem 2.1 predicts that $h(t)$ remains positive and finite and that $\\|v(t)\\|_{W^1_2}$ stays bounded on every finite time interval; any computed finite-time blow-up or $h(t)$ hitting zero before that time would refute the global-existence claim. The decisive point to inspect is the Step 4 argument: without the asserted non-autonomous parabolic smoothing estimate, the uniform bound on the maximal interval does not follow.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the abstract quasilinear parabolic well-posedness, maximal continuation, comparison, and regularity theory on which Theorem 2.1 and the global-existence argument rest."},{"cited_title":"Klapper and J","cited_arxiv_id":null,"evidence_quote":"Introduces the one-dimensional biofilm toy model and its moving-boundary formulation that the paper analyzes."},{"cited_title":"Dockery and I","cited_arxiv_id":null,"evidence_quote":"Derives the kinetic condition (1.1f) for the biofilm height from incompressibility and Darcy's law."}],"review_version":1}