Pith. sign in

REVIEW 3 major objections 4 minor 28 references

This paper proves that a degenerate chemotaxis system with indirect signal production has a global weak solution that remains bounded, for the full subcritical range, at critical mass below an explicit threshold, and in a supercritical rang

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-01 08:12 UTC pith:6TBCCEB5

load-bearing objection Strong discrete estimates, but the paper delegates the central convergence step to a preprint and states a supercritical smallness condition the proof doesn't use. the 3 major comments →

arxiv 2607.21211 v1 pith:6TBCCEB5 submitted 2026-07-23 math.AP

Global existence and boundedness for a degenerate chemotaxis system with indirect signal production via minimizing movement schemes

Tatsuya Hosono (OMU) , Philippe Lauren{\cc}ot (LAMA) This is my paper
classification math.AP MSC 35K6535K4047J3035B33
keywords global existenceboundednessdegenerate chemotaxisindirect signal productionminimizing movement schemeWasserstein gradient flowdiscrete maximal regularityMoser iteration
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper proves global existence and boundedness of weak solutions for a degenerate chemotaxis system in which cells produce an intermediate substance that in turn produces the chemoattractant. It covers the full subcritical range (diffusion exponent m above the critical value 2−4/d in dimensions d≥5, and every m>1 in dimensions 1≤d≤4), the critical case m=2−4/d in d≥5 with total mass below an explicit threshold, and a supercritical range m between 2d/(d+4) and 2−4/d provided the L^m norm of the initial cell density is below a computed value. The construction uses a time-discrete minimizing-movement scheme: the cell-density step is a gradient flow in the 2-Wasserstein metric, while the signal and intermediate equations are L^2 variational steps. Uniform bounds for the discrete approximations come from a discrete Liapunov functional, flow interchange, discrete maximal regularity, and Moser iteration, and the approximations then converge to a weak solution. If correct, the result completes the subcritical and critical existence theory for this system and adds a new small-data supercritical existence and boundedness regime.

Core claim

The paper's central claim is Theorem 1.2: for initial data in suitable L^∞/W^{2,∞} classes, system (1.3) has a global weak solution satisfying u∈L^∞((0,∞)×R^d), v∈L^∞((0,∞),W^{2,∞}), and w∈L^∞((0,∞),W^{1,∞}) in three parameter regimes: (m1) d≥5 and m>2−4/d, or 1≤d≤4 and any m>1; (m2) d≥5, m=2−4/d, and mass M below the threshold M_* defined through the optimal constant in (2.14c); (m3) d≥5, 2d/(d+4)<m<2−4/d, and ||u0||_m^m below a quantity z_m depending only on d, m, and M. The engine is Theorem 2.9, which establishes uniform-in-time-step bounds for the minimizing-movement scheme (2.1), including L^m, L^∞, W^{1,∞}, and W^{2,∞} control of the discrete approximations, by combining a discrete Li

What carries the argument

The central object is the time-discrete minimizing-movement scheme (2.1): at each step, w and v are obtained as minimizers of L^2 variational problems tied to the linear heat equations, while u is updated by minimizing a functional consisting of E[u,v]=M^{m−1}/(m−1)||u||_m^m − M∫uv plus a 2-Wasserstein displacement cost. The key identity is the discrete Liapunov decrease L[u_n,v_n,w_n]+D_0+D_1 ≤ L[u_{n−1},v_{n−1},w_{n−1}] of Lemma 2.1, with L given by (2.3). A lower bound for L provides L^m control of u in the subcritical, critical, and supercritical regimes via the optimal constant K_m^2 and threshold mass M_* (Lemma 2.2); flow interchange upgrades this to L^ρ estimates, and discrete maxima

Load-bearing premise

The load-bearing premise is that the time-discrete scheme really converges to a weak solution of the original equations — a step handed off to a companion preprint — and, in the supercritical case, that the initial data also satisfy an extra energy-smallness condition that the theorem's stated assumptions do not guarantee.

What would settle it

Find initial data satisfying (1.4)–(1.5) in d≥5 with ||u0||_m^m < z_m but L[u0,v0,w0] > f_m(θz_m) for every θ∈(0,1); for such data, Corollary 2.4(m3) cannot be applied, so Theorem 1.2(m3) would be unsupported as stated. More generally, a counterexample to the companion preprint's convergence step would invalidate the existence conclusion.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • For dimensions 1≤d≤4, global bounded weak solutions exist for every m>1 and every admissible initial data, fully covering the subcritical regime.
  • For d≥5, global bounded weak solutions exist for m>2−4/d with no size restriction, and at the critical exponent when the total mass is below the explicit threshold M_*.
  • In the supercritical range 2d/(d+4)<m<2−4/d, smallness of ||u0||_m^m below z_m guarantees global boundedness and the L^m norm of u is non-increasing in time.
  • The bounds are uniform with respect to the time step, so the same estimates hold for the discrete approximations and pass to the limit, yielding extra regularity such as ∇u^m∈L^2 and u∈L^∞.
  • The thresholds are explicit in terms of d, m, M, and optimal constants, giving computable criteria for global existence in each regime.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The supercritical case as printed may overreach: the proof's Corollary 2.4(m3) requires the additional condition L[u0,v0,w0] ≤ f_m(θz_m) for some θ∈(0,1), which is not shown to follow from the theorem's assumptions; unless that condition is implied, Theorem 1.2(m3) needs an extra hypothesis.
  • The convergence of the discrete scheme to a weak solution is deferred to a companion preprint, and Section 3 only sketches the argument; any gap there would affect Theorem 1.2, though the discrete bounds of Theorem 2.9 are self-contained.
  • The same combination of a discrete Liapunov functional, flow interchange, and discrete maximal regularity may extend to other indirect-production chemotaxis variants, such as bounded domains or nonlocal signal production, whenever a suitable lower bound for the Liapunov functional is available.
  • Computing or estimating z_m explicitly for concrete dimensions would allow comparison with known critical-mass thresholds for direct Keller–Segel systems and could predict a sharp smallness condition in the supercritical range.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies the fully parabolic degenerate chemotaxis system with indirect signal production (1.3) on R^d. The authors introduce a time-discrete minimizing movement scheme (2.1), prove a discrete Liapunov functional inequality (Lemma 2.1), derive lower bounds on the Liapunov functional in the subcritical, critical, and supercritical ranges (Lemmas 2.2 and 2.3), obtain L^m bounds (Corollary 2.4), L^\rho bounds via flow interchange and discrete maximal regularity (Lemma 2.5), and then L^\infty bounds via discrete heat-kernel smoothing and Moser iteration (Lemmas 2.7 and 2.8). These estimates are collected in Theorem 2.9. Section 3 states compactness results and then asserts convergence to a weak solution by reference to [25, Sections 6--7], without writing the limit identification or the Euler--Lagrange equation for the first discrete equation. The main theorem (Theorem 1.2) claims global bounded weak solutions in three regimes: subcritical, critical with small mass, and supercritical with small initial L^m norm.

Significance. If the convergence step were fully proved, the results would materially extend the existing theory: they cover the full subcritical range 1\le d\le 4 and d\ge 5 with m>m_*, settle the critical case m=m_* for d\ge 5 with small mass, and add a supercritical small-initial-data regime m\in(m_*,m^*). The discrete estimates are nontrivial and appear internally consistent; in particular, the Liapunov dissipation algebra, the sharp-constant identity (2.12), and the absorption/Moser exponents in Lemmas 2.5 and 2.8 are presented in detail and spot-checks support them. The paper also provides explicit, parameter-free thresholds M_* and z_m. However, the central existence claim is not actually established in the manuscript: the final passage to a weak solution is delegated to a preprint by another author, and the supercritical smallness hypothesis in Theorem 1.2(m3) is not the hypothesis used in Corollary 2.4(m3). These are load-bearing gaps, not cosmetic issues.

major comments (3)
  1. [Section 3, end of proof of Theorem 1.2] The proof of the main existence theorem stops with compactness and the sentence "we proceed as in [25, Section 7], to which we refer." The weak formulation in Definition 1.1 is never verified. In particular, the Euler--Lagrange equation for the Wasserstein minimization (2.1c) is not written down, and the passage to the limit in the degenerate flux M^{m-1}\nabla u^m and in the transport term u\nabla v is not shown. This is not a routine corollary of the uniform bounds: it requires identifying the limit of the discrete optimality conditions and handling the nonlinear flux compactness. Since Theorem 1.2 is exactly an existence theorem, this delegation to a preprint is a major gap that must be filled in the manuscript.
  2. [Theorem 1.2(m3) vs. Corollary 2.4(m3)] Theorem 1.2(m3) assumes only \|u_0\|_m^m < z_m. However, Corollary 2.4(m3) requires, in addition, L[u_0,v_0,w_0] \le f_m(\theta z_m) for some \theta\in(0,1). That extra condition is not implied by (1.4)--(1.5) together with the smallness of \|u_0\|_m^m, because v_0 and w_0 can be arbitrarily large while still satisfying the stated assumptions, and L[u_0,v_0,w_0] is not controlled by \|u_0\|_m alone. Consequently Theorem 2.9(m3) and Theorem 1.2(m3) are not supported by the proof as written. The statement of the theorem must be strengthened (or the proof must show that the L-condition follows from the stated hypotheses, which it does not).
  3. [Lemma 2.2(m2) and Theorem 1.2(m2)] If the interval in (m2) is intended to include the endpoint M=M_*, the proof does not cover it: the strict positivity of \delta_{m_*} in (2.14c) holds only for M<M_* after the equivalence using m_*-2<0. At M=M_* the lower bound (2.10b) degenerates. The statement and the proof must be reconciled: either exclude the endpoint or provide a separate approximation argument. If the intended interval is M<M_*, the notation should be made unambiguous.
minor comments (4)
  1. [Lemma 2.5, Step 1] The phrase "Classical properties of (1.5)" should refer to the porous medium equation (2.21), not to the initial-data assumption (1.5).
  2. [Throughout Section 2] The notation for the mass interval in the critical case is inconsistent: Theorem 1.2(m2), Corollary 2.4(m2), and Lemma 2.2(m2) should all use the same open/closed bracket convention, and the endpoint M=M_* should be clearly excluded if it is not covered.
  3. [References] Since the proof relies on [25] for the essential convergence argument, the preprint status of [25] should be prominently flagged in the text and the specific arguments borrowed should be stated, not merely referenced.
  4. [Display (2.32)] The formula has a minor typo: the term involving \|v_0-\Delta v_0\|_m is written with an extra factor and a missing power in the displayed line. This does not affect the subsequent estimates but should be corrected.

Circularity Check

0 steps flagged

No circular derivation: constants are explicit variational thresholds, estimates come from external tools, and the delegated convergence argument is an incompleteness rather than a self-referential reduction.

full rationale

The paper's central claims are not obtained by fitting parameters to the quantities they purport to predict, nor by importing a self-citation as the load-bearing premise. The thresholds M* and z_m are explicitly defined in (2.14c)-(2.14d) in terms of K_m, a well-defined supremum in (2.12); they are not tuned to the constructed solution. The lower bounds on the Liapunov functional and the resulting a priori estimates in Corollary 2.4 and Theorem 2.9 are derived from functional inequalities with stated constants, and the flow interchange, discrete maximal regularity, and Moser iteration arguments rely on external sources [2, 3, 19, 23, 28] and standard auxiliary lemmas [5, 20]. The self-citations that occur are not load-bearing in a circular way: [20, Lemma A.1] and [5, Proposition 8] are technical iteration/limiting lemmas from published work and do not contain the theorem's conclusion. The most serious issue flagged in the manuscript is that the convergence of the time-discrete scheme to a weak solution is not proved in detail; Section 3 states “we proceed as in [25, Section 7], to which we refer.” This is an omission or external delegation, not a circular reduction: the missing limit identification is an independent step, and relying on another author's preprint does not make the conclusion equivalent to its own inputs. Similarly, the mismatch in case (m3), where Corollary 2.4(m3) requires the additional condition L[u0,v0,w0] ≤ f_m(θ z_m) while Theorem 1.2 states only ‖u0‖_m^m < z_m, is a possible overreach or missing hypothesis, but it is not a fitted-parameter or self-definitional circularity. No construction in the paper reduces a claimed prediction to the data that define it.

Axiom & Free-Parameter Ledger

0 free parameters · 7 axioms · 0 invented entities

No constants are fitted to data or chosen ad hoc: M* (2.14c), z_m (2.14d), δ_m and δ_{m,θ} (2.14a),(2.14e) are explicit functions of the sharp variational constant K_m (2.12); the initial mass M is part of the data. The burden is entirely in the axioms above: (i) discrete-scheme well-posedness from [25]; (ii) W_2 gradient-flow/flow-interchange theory from [2] and [5]; (iii) discrete maximal regularity from [3,19]; (iv) standard interpolation lemmas from [28] and the author's [20]; (v) the transfer of [25]'s convergence argument. Most of these are external and legitimate; (v) is the fragile one — a soundness/reproducibility concern, not circularity.

axioms (7)
  • domain assumption The time-discrete scheme (2.1) is well-posed: minimizers exist for each step and enjoy u_n ∈ L^2, v_n ∈ W^{4,2}, w_n ∈ W^{2,2}, with v_n,w_n nonnegative.
    Invoked at the start of Section 2 ('Arguing as in [25, Proposition 2.2]') and used throughout (e.g., nonnegativity of v_n in (2.24), regularity in Lemmas 2.5–2.8). The proof is not reproduced; it is taken from a preprint.
  • domain assumption Flow-interchange inequality (2.23) from [2, Theorem 11.1.4] applied to the porous-medium flow (2.21), including the t→0 passage in (2.25) deferred to [5, Proposition 8, Step 3].
    Load-bearing for the entire L^ρ estimate (2.20): without (2.25) there is no Lyapunov-type control of ||u_n||_ρ. Requires the PME to be a W_2-gradient flow of F[z]=||z||_ρ^ρ/(ρ−1) and suitable differentiability at t=0.
  • domain assumption Discrete maximal regularity estimates (A.3)/(A.8) from [3,19] hold on R^d for the fully implicit scheme z_n = z_{n−1} + τ(∆z_n − z_n + f_n).
    Used in Lemma 2.5 Steps 3–4 and Lemma A.1 to convert L^s bounds on f_n into L^q bounds on z_n −∆z_n with q>1; on the whole space this rests on maximal regularity for sectorial operators and A-stable discretizations, cited from [3,19].
  • standard math Gagliardo–Nirenberg interpolation lemma [28, Lemma 2.4] holds with the stated exponents (2.26), (2.35), (2.49).
    Used repeatedly in Lemmas 2.5 and 2.8; standard for these porous-medium-type interpolations, but the paper does not reproduce the lemma.
  • standard math Bessel potential bound ||(I−∆)^{−2}u||_{m/(m−1)} ≤ C||u||_1 (Lemma 2.3) and the continuous embeddings of W^{2,2}, W^{2,s} used in Lemma 2.5 Step 4.
    Needed for the Liapunov lower bound in 1≤d≤4 and for the r′ estimates in Step 4; standard harmonic-analysis facts, stated with the required exponent ranges.
  • standard math Moser-iteration lemma [20, Lemma A.1] used to close Lemma 2.8.
    An elementary discrete iteration lemma from one of the authors' 1994 paper; not proven here, no independent verification offered.
  • domain assumption The convergence argument of [25, Sections 6–7] transfers verbatim to the present estimates (Propositions 3.1–3.3), including existence of the limit triple and identification of the limit as a weak solution of (1.3) in the sense of Definition 1.1.
    Paper states convergence 'follows very closely' [25] and completes Theorem 1.2 'as in [25, Section 7]'; the discrete Euler–Lagrange equation for u_n is not written down and the passage to the limit in the u-equation is not performed here.

pith-pipeline@v1.3.0-alltime-deepseek · 36874 in / 48185 out tokens · 435353 ms · 2026-08-01T08:12:47.787890+00:00 · methodology

0 comments
read the original abstract

Global existence and boundedness of weak solutions for a fully parabolic degenerate chemotaxis system with indirect signal production are proved for any initial data in the subcritical case and under smallness conditions in the critical and supercritical cases. To construct weak solutions, a time discrete scheme is set up, for which the first equation has a gradient flow structure with respect to the 2-Wasserstein distance, while the other two equations feature an L\textsuperscript{2}-variational structure. The proof relies in particular on the flow interchange method and discrete maximal regularity.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 2 linked inside Pith

  1. [1]

    N. D. Alikakos , L p \ bounds of solutions of reaction-diffusion equations , Comm. Partial Differential Equations, 4 (1979), pp. 827--868

    1979

  2. [2]

    Ambrosio, N

    L. Ambrosio, N. Gigli, and G. Savar \'e , Gradient flows in metric spaces and in the space of probability measures , Basel: Birkh \"a user, 2005

    2005

  3. [3]

    Ashyralyev, S

    A. Ashyralyev, S. Piskarev, and L. Weis , On well-posedness of difference schemes for abstract parabolic equations in L^p([0,T];E) spaces , Numer. Funct. Anal. Optim., 23 (2002), pp. 669--693

    2002

  4. [4]

    Blanchet, J

    A. Blanchet, J. A. Carrillo, D. Kinderlehrer, M. Kowalczyk, Ph . Lauren c ot, and S. Lisini , A hybrid variational principle for the K eller- S egel system in R^2 , ESAIM Math. Model. Numer. Anal., 49 (2015), pp. 1553--1576

    2015

  5. [5]

    Blanchet and Ph

    A. Blanchet and Ph . Lauren c ot , The parabolic-parabolic K eller- S egel system with critical diffusion as a gradient flow in R^d,\ d 3 , Comm. Partial Differential Equations, 38 (2013), pp. 658--686

    2013

  6. [6]

    Chen, J.-G

    L. Chen, J.-G. Liu, and J. Wang , Multidimensional degenerate K eller- S egel system with critical diffusion exponent 2n/(n+2) , SIAM J. Math. Anal., 44 (2012), pp. 1077--1102

    2012

  7. [7]

    Chen and J

    L. Chen and J. Wang , Exact criterion for global existence and blow up to a degenerate K eller- S egel system , Doc. Math., 19 (2014), pp. 103--120

    2014

  8. [8]

    Ding and W

    M. Ding and W. Wang , Global boundedness in a quasilinear fully parabolic chemotaxis system with indirect signal production , Discrete Contin. Dyn. Syst. Ser. B, 24 (2019), pp. 4665--4684

    2019

  9. [9]

    Engel and R

    K.-J. Engel and R. Nagel , One-parameter semigroups for linear evolution equations , vol. 194 of Graduate Texts in Mathematics, Springer-Verlag, New York, 2000. With contributions by S. Brendle, M. Campiti, T. Hahn, G. Metafune, G. Nickel, D. Pallara, C. Perazzoli, A. Rhandi, S. Romanelli and R. Schnaubelt

    2000

  10. [10]

    Fujie and T

    K. Fujie and T. Senba , Application of an A dams type inequality to a two-chemical substances chemotaxis system , J. Differential Equations, 263 (2017), pp. 88--148

    2017

  11. [11]

    Differential Equations, 266 (2019), pp

    height 2pt depth -1.6pt width 23pt, Blowup of solutions to a two-chemical substances chemotaxis system in the critical dimension , J. Differential Equations, 266 (2019), pp. 942--976

    2019

  12. [12]

    Herda, A

    M. Herda, A. Trescases, and A. Zurek , A finite volume scheme for the local sensing chemotaxis model , SMAI J. Comput. Math., 11 (2025), pp. 637--676

    2025

  13. [13]

    Hosono , On the C auchy problem of a chemotaxis system with indirect signal production , in Mathematical A nalysis and A pproximation of PDE - C hemotaxis M odels, vol

    T. Hosono , On the C auchy problem of a chemotaxis system with indirect signal production , in Mathematical A nalysis and A pproximation of PDE - C hemotaxis M odels, vol. 43 of SEMA SIMAI Springer Ser., Springer, Cham, 2026, pp. 61--82

    2026

  14. [14]

    Hosono and Ph

    T. Hosono and Ph . Lauren c ot , Global existence and boundedness of solutions to a fully parabolic chemotaxis system with indirect signal production in R^4 , J. Differential Equations, 416 (2025), pp. 2085--2133

    2025

  15. [15]

    Ishida and T

    S. Ishida and T. Yokota , Boundedness in a quasilinear fully parabolic K eller- S egel system via maximal S obolev regularity , Discrete Contin. Dyn. Syst. Ser. S, 13 (2020), pp. 212--232

    2020

  16. [16]

    Jordan, D

    R. Jordan, D. Kinderlehrer, and F. Otto , The variational formulation of the F okker- P lanck equation , SIAM J. Math. Anal., 29 (1998), pp. 1--17

    1998

  17. [17]

    E. F. Keller and L. A. Segel , Model for chemotaxis , J. Theoret. Biol., 30 (1971), pp. 225--234

    1971

  18. [18]

    Kimijima, K

    A. Kimijima, K. Nakagawa, and T. Ogawa , Threshold of global behavior of solutions to a degenerate drift-diffusion system in between two critical exponents , Calc. Var. Partial Differential Equations, 53 (2015), pp. 441--472

    2015

  19. [19]

    Kov \'a cs, B

    B. Kov \'a cs, B. Li, and C. Lubich , A-stable time discretizations preserve maximal parabolic regularity , SIAM J. Numer. Anal., 54 (2016), pp. 3600--3624

    2016

  20. [20]

    Lauren c ot , Solutions to a P enrose- F ife model of phase-field type , J

    Ph . Lauren c ot , Solutions to a P enrose- F ife model of phase-field type , J. Math. Anal. Appl., 185 (1994), pp. 262--274

    1994

  21. [21]

    Lin and S

    K. Lin and S. Wang , Sharp critical mass criterion for the fully parabolic K eller- S egel system with the intermediate exponent , Nonlinearity, 39 (2026), pp. Paper No. 015008, 25

    2026

  22. [22]

    X. Mao, M. Liu, and Y. Li , Finite-time blowup in a fully parabolic chemotaxis model involving indirect signal production . Preprint, arXiv :2503.12439 [math. AP ] (2025), 2025

    Pith/arXiv arXiv 2025

  23. [23]

    Matthes, R

    D. Matthes, R. J. McCann, and G. Savar \'e , A family of nonlinear fourth order equations of gradient flow type , Commun. Partial Differ. Equations, 34 (2009), pp. 1352--1397

    2009

  24. [24]

    Mimura , The variational formulation of the fully parabolic K eller- S egel system with degenerate diffusion , J

    Y. Mimura , The variational formulation of the fully parabolic K eller- S egel system with degenerate diffusion , J. Differential Equations, 263 (2017), pp. 1477--1521

    2017

  25. [25]

    Preprint, arXiv :2406.14536 [math

    height 2pt depth -1.6pt width 23pt, Formulation of C himera gradient flows for chemotaxis systems with indirect signal production and degenerate diffusion . Preprint, arXiv :2406.14536 [math. AP ] (2024), 2024

    Pith/arXiv arXiv 2024

  26. [26]

    height 2pt depth -1.6pt width 23pt, Global existence of solutions to parabolic-parabolic Keller - Segel system in between two critical exponents , Adv. Math. Sci. Appl., 33 (2024), pp. 77--96

    2024

  27. [27]

    Ogawa , The degenerate drift-diffusion system with the S obolev critical exponent , Discrete Contin

    T. Ogawa , The degenerate drift-diffusion system with the S obolev critical exponent , Discrete Contin. Dyn. Syst. Ser. S, 4 (2011), pp. 875--886

    2011

  28. [28]

    Sugiyama , Global existence in sub-critical cases and finite time blow-up in super-critical cases to degenerate Keller - Segel systems

    Y. Sugiyama , Global existence in sub-critical cases and finite time blow-up in super-critical cases to degenerate Keller - Segel systems. , Differ. Integral Equ., 19 (2006), pp. 841--876

    2006