pith. sign in

arxiv: 2604.27561 · v1 · submitted 2026-04-30 · 🧮 math.AP

Finite-time blow-up in a class of chemotaxis systems with spatially heterogeneous diffusion sensitivity

Yashuang Zhao , Shijun Li , Shaopeng Xu This is my paper

Pith reviewed 2026-05-07 09:42 UTC · model grok-4.3

classification 🧮 math.AP
keywords chemotaxisKeller-Segelfinite-time blow-upheterogeneous diffusionradial symmetryparabolic-elliptic systemblow-up criterion
0
0 comments X

The pith

In this chemotaxis system with |x|^β diffusion, sufficiently concentrated radial initial mass produces finite-time blow-up.

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

The paper studies a parabolic-elliptic Keller-Segel model in which diffusion sensitivity varies with position through the coefficient |x|^β. It first constructs radially symmetric classical solutions for nonconstant nonnegative radial initial data and establishes their local existence, boundedness, and uniqueness under a compatibility condition when the data are radially decreasing. The central result shows that when the initial mass is large enough and sufficiently concentrated near the origin, these solutions cannot remain bounded and instead develop a singularity in finite time. This establishes blow-up for a heterogeneous variant of the standard chemotaxis system inside a ball.

Core claim

For the system u_t = ∇ · (|x|^β ∇u) − ∇ · (u^α ∇v) with Δv = μ − u under Neumann conditions in the ball Ω = B_R(0), any nonconstant nonnegative radial initial datum admits a radially symmetric classical solution in (Ω ∖ {0}) × (0,T); when the datum is C^{1+θ}, radially decreasing, and satisfies the compatibility criterion, the solution is unique and bounded up to T* < T, yet if the initial mass is sufficiently concentrated the solution blows up in finite time.

What carries the argument

The radially symmetric, radially decreasing initial data satisfying the compatibility criterion, which interacts with the singular diffusion coefficient |x|^β to control the evolution toward concentration at the origin.

If this is right

  • Blow-up can occur despite the spatially heterogeneous diffusion term.
  • The origin becomes the only possible location of the singularity because of radial symmetry and the form of the diffusion coefficient.
  • The maximal existence time T* is finite precisely when the initial mass meets the concentration criterion.
  • The result extends classical finite-time blow-up statements to models with position-dependent sensitivity.

Where Pith is reading between the lines

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

  • The same concentration mechanism might produce blow-up for initial data that are only approximately radial if the asymmetry is small.
  • Numerical integration of the system could locate the critical mass value separating global existence from blow-up.
  • Similar blow-up statements may hold when the domain is replaced by an annulus that avoids the origin singularity.

Load-bearing premise

The initial data must be radially symmetric and decreasing inside a ball, satisfy a compatibility criterion, and the diffusion coefficient may be singular at the origin.

What would settle it

A radially decreasing initial datum whose total mass exceeds the concentration threshold but whose solution remains bounded for all positive times would contradict the blow-up statement.

read the original abstract

\indent In this paper, we study a class of parabolic-elliptic Keller-Segel systems with diffusion sensitivity dependent on spatial position, given by type \begin{equation} \left\{ \begin{array}{ll} u_{t} = \bigtriangledown\cdot(|x|^{\beta} \bigtriangledown u)-\bigtriangledown\cdot(u^{\alpha} \bigtriangledown v), 0=\bigtriangleup v-\mu +u, \qquad \mu:=\frac{1}{|\Omega|}\int_{\Omega}udx,\end{array}\right. \end{equation} under homogeneous Neumann conditions in a ball $\Omega=B_{R}(0)\subset \mathbb{R}^{n}$ with $\alpha \ge 1$, $\beta>0$ and $n\ge 2$.\par \indent It is proved that any nonconstant nonnegative radial initial data $u_{0}\in C^{\theta}(\overline{\Omega})$, where $\theta \in (0,1)$, there exists a radially symmetric classical solution of the system (0.1) in $(\Omega \setminus \{ 0 \})\times (0,T)$ for some $T>0$; moreover, if the initial values $u_{0}\in C^{1+\theta}(\overline{\Omega})$ for some $\theta \in (0,1)$ and satisfy a certain compatibility criterion and are radially decreasing, then this solution is bounded and unique in $(\Omega \setminus \{ 0 \})\times (0,T^{*})$ with $T^{*}<T$.\par Finally, it is found that the initial mass corresponding to this parabolic-elliptic problem (0.1) is sufficiently concentrated to allow the solution to blow up in finite time.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript analyzes a parabolic-elliptic Keller-Segel system with spatially heterogeneous diffusion sensitivity |x|^β (β>0) in a ball Ω=B_R(0)⊂R^n (n≥2), subject to homogeneous Neumann boundary conditions. It proves local-in-time existence of radially symmetric classical solutions in the punctured domain (Ω∖{0})×(0,T) for any nonconstant nonnegative radial initial datum u_0∈C^θ(Ω̄) with θ∈(0,1). Under the additional assumptions that u_0∈C^{1+θ}(Ω̄) is radially decreasing and satisfies a compatibility criterion, the solution is asserted to be bounded and unique on (Ω∖{0})×(0,T*) for some T*<T. The paper further claims that when the initial mass is sufficiently concentrated, the solution blows up in finite time.

Significance. If the technical details are completed without gaps, the results would extend finite-time blow-up criteria for chemotaxis models to the setting of singular diffusion coefficients at the origin. The radial symmetry and compatibility assumptions permit explicit moment or ODE comparisons that are unavailable in the nonradial case. The work is of moderate interest for mathematical biology applications involving heterogeneous media, though the restriction to balls and radial data narrows its scope. The paper does not appear to contain machine-checked proofs or fully parameter-free derivations.

major comments (2)
  1. [Abstract] Abstract: The claim that the solution 'is bounded and unique in (Ω∖{0})×(0,T*) with T*<T' while simultaneously asserting finite-time blow-up for concentrated mass is internally inconsistent with standard local-existence/continuation theory for quasilinear parabolic systems. If the solution remains bounded on [0,T*), the maximal existence time can be extended beyond T*, contradicting the asserted maximality of T*. The manuscript must clarify the precise definition of T* (e.g., whether it is the first time a norm diverges) and explain how boundedness on the whole interval up to a finite T* is compatible with blow-up; this logical step is load-bearing for the central finite-time blow-up conclusion.
  2. [Abstract and §2 (local existence)] The local existence statement in the punctured domain (Ω∖{0})×(0,T) and the subsequent boundedness result rely on radial symmetry and a compatibility criterion whose explicit form is not stated in the abstract. In the main text, the criterion must be written out (e.g., as a boundary or regularity condition at x=0) and shown to be sufficient to control the singular diffusion term |x|^β near the origin; without this, the handling of the degeneracy at x=0 cannot be verified.
minor comments (2)
  1. [Abstract] The abstract sentence 'it is found that the initial mass corresponding to this parabolic-elliptic problem (0.1) is sufficiently concentrated to allow the solution to blow up in finite time' is grammatically incomplete and should be rephrased to state the precise mass threshold or concentration condition.
  2. [Abstract] Replace the nonstandard symbol 'bigtriangledown' with the conventional '∇' throughout the displayed system (0.1).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and insightful comments on our manuscript. We address the major comments below and will make the necessary revisions to clarify the abstract and strengthen the presentation of the local existence and compatibility results.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The claim that the solution 'is bounded and unique in (Ω∖{0})×(0,T*) with T*<T' while simultaneously asserting finite-time blow-up for concentrated mass is internally inconsistent with standard local-existence/continuation theory for quasilinear parabolic systems. If the solution remains bounded on [0,T*), the maximal existence time can be extended beyond T*, contradicting the asserted maximality of T*. The manuscript must clarify the precise definition of T* (e.g., whether it is the first time a norm diverges) and explain how boundedness on the whole interval up to a finite T* is compatible with blow-up; this logical step is load-bearing for the central finite-time blow-up conclusion.

    Authors: We acknowledge that the abstract wording is ambiguous and potentially misleading. In the revised manuscript we will explicitly distinguish T (a time of local existence from the local existence theorem) from T* (the maximal existence time). Under the radial decreasing assumption and compatibility criterion the solution remains classical, hence bounded on every compact subinterval [0,t] with t<T*. We prove finite-time blow-up by showing that for sufficiently concentrated initial data the maximal time T* must be finite, with ||u(·,t)||_∞ → ∞ as t → T*−. This is compatible with standard theory because the L^∞ bound is not uniform on [0,T*); the continuation criterion fails precisely when the norm diverges. We will rewrite the abstract to state the definition of T* and the precise meaning of blow-up. revision: yes

  2. Referee: [Abstract and §2 (local existence)] The local existence statement in the punctured domain (Ω∖{0})×(0,T) and the subsequent boundedness result rely on radial symmetry and a compatibility criterion whose explicit form is not stated in the abstract. In the main text, the criterion must be written out (e.g., as a boundary or regularity condition at x=0) and shown to be sufficient to control the singular diffusion term |x|^β near the origin; without this, the handling of the degeneracy at x=0 cannot be verified.

    Authors: We agree that the compatibility criterion needs to be stated explicitly and its role in controlling the degeneracy made transparent. In the revised version we will include the precise statement of the compatibility criterion already in the abstract and restate it verbatim at the beginning of Section 2. We will add a dedicated paragraph (and the corresponding estimates) showing that the criterion, together with radial symmetry, yields sufficient Hölder regularity at the origin to absorb the singular factor |x|^β into the diffusion term, thereby justifying the classical solution concept in the punctured domain. These estimates will be written out in full detail. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper establishes local existence of radially symmetric classical solutions for nonconstant radial initial data, proves boundedness and uniqueness on (0,T*) under radial-decreasing C^{1+θ} data satisfying a compatibility criterion, and concludes finite-time blow-up when initial mass is sufficiently concentrated. These steps rely on direct PDE analysis (local existence theory, continuation criteria, and presumably moment or energy methods for blow-up) without reducing any claim to a fitted parameter, self-definition, or load-bearing self-citation. The abstract and structure present independent mathematical arguments that do not collapse to tautological inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work relies on classical local existence theory for parabolic-elliptic systems and radial symmetry preservation; no new free parameters or invented entities are introduced.

axioms (2)
  • standard math Local existence and continuation theory for quasilinear parabolic-elliptic systems with Neumann boundary conditions
    Invoked to obtain short-time classical solutions away from the origin.
  • domain assumption Radial symmetry of solutions is preserved when initial data are radial
    Used throughout to reduce the problem to a one-dimensional radial setting.

pith-pipeline@v0.9.0 · 5640 in / 1274 out tokens · 32299 ms · 2026-05-07T09:42:37.208969+00:00 · methodology

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

32 extracted references · 32 canonical work pages

  1. [1]

    Clifford S. Patlak. Random walk with persistence and external bias.The Bulletin of Math- ematical Biophysics, 15(3):311–338, 1953

    work page 1953

  2. [2]

    Keller and Lee A

    Evelyn F. Keller and Lee A. Segel. Model for chemotaxis.Journal of Theoretical Biology, 30(2):225–234, 1971

    work page 1971

  3. [3]

    Signaling mechanisms for regulation of chemotaxis.Cell Research, 15(1):52– 56, 2005

    Dianqing Wu. Signaling mechanisms for regulation of chemotaxis.Cell Research, 15(1):52– 56, 2005

    work page 2005

  4. [4]

    Free-flight responses of drosophila melanogaster to attractive odors.The Journal of experimental biology, 209(15):3001–3017, 2006

    Seth A Budick and Michael H Dickinson. Free-flight responses of drosophila melanogaster to attractive odors.The Journal of experimental biology, 209(15):3001–3017, 2006

    work page 2006

  5. [5]

    J. S. Kennedy and David Marsh. Pheromone-regulated anemotaxis in flying moths.Science, 184(4140):1001–999, 1974

    work page 1974

  6. [6]

    Signaling pathways induced by vascular endothelial growth factor

    B Larriv´ ee and A Karsan. Signaling pathways induced by vascular endothelial growth factor. International journal of molecular medicine, 5:447–56, 06 2000

    work page 2000

  7. [7]

    Finite dimensional attractor for one-dimensional keller-segel equations.Funkcialaj Ekvacioj, 44(3):441–469, 2001

    Koichi Osaki and Atsushi Yagi. Finite dimensional attractor for one-dimensional keller-segel equations.Funkcialaj Ekvacioj, 44(3):441–469, 2001

    work page 2001

  8. [8]

    Nagai, T

    T. Nagai, T. Senba, and K. Yoshida. Application of the trudinger-moser inequality to a parabolic system of chemotaxis.Funkc Ekvacioj, 40(3):411–433, 1997

    work page 1997

  9. [9]

    Aggregation vs

    Michael Winkler. Aggregation vs. global diffusive behavior in the higher-dimensional keller–segel model.Journal of Differential Equations, 248(12):2889–2905, 2010

    work page 2010

  10. [10]

    Toward a mathematical theory of keller–segel models of pattern formation in biological tissues.Mathematical Models and Methods in Applied Sciences, 25(9):1663–1763, 2015

    Nicola Bellomo, Nicola Bellomo, Abdelghani Bellouquid, Youshan Tao, and Michael Win- kler. Toward a mathematical theory of keller–segel models of pattern formation in biological tissues.Mathematical Models and Methods in Applied Sciences, 25(9):1663–1763, 2015

    work page 2015

  11. [11]

    Thomas Hillen and Kevin J. Painter. A user’s guide to pde models for chemotaxis.Journal of Mathematical Biology, 58(1-2):183–217, 2009

    work page 2009

  12. [12]

    Horstmann

    D. Horstmann. From 1970 until present: The keller-segel model in chemotaxis and its consequences i.Jahresbericht der Deutschen Mathematiker-Vereinigung, 105(3):103–165, 2003

    work page 1970

  13. [13]

    Blowup of nonradial solutions to parabolic–elliptic systems model- ing chemotaxis in two-dimensional domains.Journal of Inequalities and Applications, 2001(1):37–55, 2001

    Toshitaka Nagai. Blowup of nonradial solutions to parabolic–elliptic systems model- ing chemotaxis in two-dimensional domains.Journal of Inequalities and Applications, 2001(1):37–55, 2001

    work page 2001

  14. [14]

    A blow-up mechanism for a chemotaxis model

    Miguel A Herrero and Juan JL Vel´ azquez. A blow-up mechanism for a chemotaxis model. Annali della Scuola Normale Superiore di Pisa-Classe di Scienze, 24(4):633–683, 1997

    work page 1997

  15. [15]

    Blow-up in a chemotaxis model without symmetry assumptions.European Journal of Applied Mathematics, 12(2):159–177, 2001

    Dirk Horstmann and Guofang Wang. Blow-up in a chemotaxis model without symmetry assumptions.European Journal of Applied Mathematics, 12(2):159–177, 2001

    work page 2001

  16. [16]

    Does a’volume-filling effect’always prevent chemotactic collapse?Mathe- matical Methods in the Applied Sciences, 33(1):12–24, 2010

    Michael Winkler. Does a’volume-filling effect’always prevent chemotactic collapse?Mathe- matical Methods in the Applied Sciences, 33(1):12–24, 2010

    work page 2010

  17. [17]

    Boundedness vs

    Dirk Horstmann and Michael Winkler. Boundedness vs. blow-up in a chemotaxis system. Journal of Differential Equations, 215(1):52–107, 2005. 25

    work page 2005

  18. [18]

    On explosions of solutions to a system of partial differ- ential equations modelling chemotaxis.Transactions of the american mathematical society, 329(2):819–824, 1992

    Willi J¨ ager and Stephan Luckhaus. On explosions of solutions to a system of partial differ- ential equations modelling chemotaxis.Transactions of the american mathematical society, 329(2):819–824, 1992

    work page 1992

  19. [19]

    Michael Winkler. How unstable is spatial homogeneity in keller-segel systems? a new critical mass phenomenon in two-and higher-dimensional parabolic-elliptic cases.Mathematische Annalen, 373(10):1237–1282, 2019

    work page 2019

  20. [20]

    Boundedness and finite-time collapse in a chemotaxis system with volume-filling effect.Nonlinear Analysis: Theory, Methods & Applications, 72(2):1044–1064, 2010

    Michael Winkler and Kianhwa C Djie. Boundedness and finite-time collapse in a chemotaxis system with volume-filling effect.Nonlinear Analysis: Theory, Methods & Applications, 72(2):1044–1064, 2010

    work page 2010

  21. [21]

    Slow grow-up in a quasilinear keller–segel system.Journal of Dynamics and Differential Equations, 36(2):1677–1702, 2024

    Michael Winkler. Slow grow-up in a quasilinear keller–segel system.Journal of Dynamics and Differential Equations, 36(2):1677–1702, 2024

    work page 2024

  22. [22]

    Complete infinite-time mass aggregation in a quasilinear keller–segel sys- tem.Israel Journal of Mathematics, 263(1):93–127, 2024

    Michael Winkler. Complete infinite-time mass aggregation in a quasilinear keller–segel sys- tem.Israel Journal of Mathematics, 263(1):93–127, 2024

    work page 2024

  23. [23]

    Xinyu Tu, Chunlai Mu, and Pan Zheng. On effects of the nonlinear signal production to the boundedness and finite-time blow-up in a flux-limited chemotaxis model.Mathematical Models and Methods in Applied Sciences, 32(4):647–711, 2022

    work page 2022

  24. [24]

    Finite-time blow-up in a degenerate chemotaxis system with flux limitation.Transactions of the American mathematical society, Series B, 4(2):31– 67, 2017

    Nicola Bellomo and Michael Winkler. Finite-time blow-up in a degenerate chemotaxis system with flux limitation.Transactions of the American mathematical society, Series B, 4(2):31– 67, 2017

    work page 2017

  25. [25]

    Finite-time blow-up in a quasi- linear degenerate chemotaxis system with flux limitation.Acta Applicandae Mathematicae, 167(1):231–259, 2020

    Yuka Chiyoda, Masaaki Mizukami, and Tomomi Yokota. Finite-time blow-up in a quasi- linear degenerate chemotaxis system with flux limitation.Acta Applicandae Mathematicae, 167(1):231–259, 2020

    work page 2020

  26. [26]

    Finite time blow-up for a one-dimensional quasilin- ear parabolic–parabolic chemotaxis system.Annales de l’Institut Henri Poincar´ e C, Analyse non lin´ eaire, 27(1):437–446, 2010

    Tomasz Cie´ slak and Philippe Lauren¸ cot. Finite time blow-up for a one-dimensional quasilin- ear parabolic–parabolic chemotaxis system.Annales de l’Institut Henri Poincar´ e C, Analyse non lin´ eaire, 27(1):437–446, 2010

    work page 2010

  27. [27]

    Solutions to a chemotaxis system with spatially heterogeneous diffusion sensitivity.Journal of Differential Equations, 424:1–29, 2025

    Gregor Fl¨ uchter. Solutions to a chemotaxis system with spatially heterogeneous diffusion sensitivity.Journal of Differential Equations, 424:1–29, 2025

    work page 2025

  28. [28]

    Blow-up in a higher-dimensional chemotaxis system despite logistic growth restriction.Journal of Mathematical Analysis and Applications, 384(2):261–272, 2011

    Michael Winkler. Blow-up in a higher-dimensional chemotaxis system despite logistic growth restriction.Journal of Mathematical Analysis and Applications, 384(2):261–272, 2011

    work page 2011

  29. [29]

    Finite-time blow-up in a quasilinear system of chemo- taxis.Nonlinearity, 21(5):1057–1076, 2008

    Tomasz Cie´ slak and Michael Winkler. Finite-time blow-up in a quasilinear system of chemo- taxis.Nonlinearity, 21(5):1057–1076, 2008

    work page 2008

  30. [30]

    A blow-up result for a quasilinear chemotaxis system with logistic source in higher dimensions.Journal of Mathematical Analysis and Applica- tions, 464(1):435–455, 2018

    Ke Lin, Chunlai Mu, and Hua Zhong. A blow-up result for a quasilinear chemotaxis system with logistic source in higher dimensions.Journal of Mathematical Analysis and Applica- tions, 464(1):435–455, 2018

    work page 2018

  31. [31]

    Michael Winkler. Approaching logarithmic singularities in quasilinear chemotaxis- consumption systems with signal-dependent sensitivities.Discrete & Continuous Dynamical Systems-Series B, 27(11):6565–6587, 2022

    work page 2022

  32. [32]

    L ∞ bounds in a two-dimensional doubly degenerate nutrient taxis system with general cross-diffusive flux.Journal of Differential Equations, 400(10):423–456, 2024

    Michael Winkler. L ∞ bounds in a two-dimensional doubly degenerate nutrient taxis system with general cross-diffusive flux.Journal of Differential Equations, 400(10):423–456, 2024. 26

    work page 2024