Chemotaxis models with signal-dependent sensitivity and a logistic-type source, II: Persistence and stabilization
Pith reviewed 2026-05-13 19:57 UTC · model grok-4.3
The pith
A threshold χ^*(u^*) on the sensitivity parameter separates stable and unstable regimes for the positive equilibrium.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper identifies a threshold χ^*(u^*) such that the unique positive equilibrium (u^*, v^*) is linearly stable when χ0 < χ^*(u^*), producing local exponential decay, and linearly unstable when χ0 > χ^*(u^*). It supplies conditions under which every bounded solution converges exponentially to this equilibrium. For the mass-constrained case a = b = 0 the same linearization technique produces a stability threshold for constant states, and the Lyapunov and rectangle/ODE comparison methods are extended to the general parameter regime m ≥ 1 and β > 0.
What carries the argument
The threshold function χ^*(u^*) obtained by linearizing the PDE system at the positive equilibrium and determining the sign of the principal eigenvalue; it directly controls whether small perturbations decay or grow.
If this is right
- When χ0 is below χ^*(u^*) small perturbations of the equilibrium decay exponentially to zero.
- Uniform persistence holds for all positive initial data whenever m ≥ 1.
- Large values of the saturation exponent β or negative χ0 suppress aggregation and favor relaxation to the constant state.
- Under the stated conditions on a, b, α, γ, β the entire orbit of any bounded solution approaches the equilibrium exponentially fast.
Where Pith is reading between the lines
- The explicit threshold may be used to locate the onset of spatial patterns by tracking when χ0 crosses χ^*(u^*).
- The stability criterion supplies a concrete parameter regime in which biological aggregation models are guaranteed to relax rather than form clusters.
- The same linearization technique could be tested on related systems that replace the logistic source with other growth functions.
Load-bearing premise
Global existence and uniform boundedness of solutions, which were proved in Part I of the series.
What would settle it
Explicit computation of χ^*(u^*) for concrete parameter values followed by numerical integration of the system showing convergence when χ0 lies below the threshold and divergence or pattern formation when χ0 lies above it.
read the original abstract
This paper is Part II of a series on global existence and asymptotic behavior of positive solutions to \begin{equation*} \begin{cases} \displaystyle u_t=\Delta u-\chi_0\nabla\cdot\left(\frac{u^m}{(1+v)^\beta}\nabla v\right)+au-bu^{1+\alpha}, & x\in\Omega, \cr \displaystyle 0=\Delta v-\mu v+\nu u^\gamma, & x\in\Omega, \cr \displaystyle \frac{\partial u}{\partial n}=\frac{\partial v}{\partial n}=0, & x\in\partial\Omega, \end{cases} \end{equation*} where $\Omega\subset\mathbb{R}^N$ is a bounded and smooth domain. The parameters $\alpha,\gamma,m,\mu,\nu$ are positive, $\chi_0$ is real, and $a,b,\beta$ are nonnegative. In Part I, we established boundedness and global existence. Here, we study persistence and stabilization, quantifying how $\beta$ and $\chi_0$ influence long-time dynamics. First, we prove uniform persistence if $m\ge 1$. Next, for $a,b>0$, the unique positive equilibrium is $(u^*,v^*) = \left((\tfrac{a}{b})^{1/\alpha},(\tfrac{\nu}{\mu})(\tfrac{a}{b})^{\gamma/\alpha}\right)$. We identify a threshold $\chi^*(u^*)$: $(u^*,v^*)$ is linearly stable if $\chi_0<\chi^*(u^*)$, with local exponential decay, unstable if $\chi_0>\chi^*(u^*)$. We also give conditions ensuring every bounded solution converges exponentially to $(u^*,v^*)$. For $a=b=0$, we study stability of the constant equilibria under mass constraint, obtaining a linear stability threshold and global stabilization. We extend the Lyapunov method from $m=1$ to $m>1$ and the rectangle/ODE method from $\beta=0$ to $\beta>0$. For $m\ge 1$, signal saturation (large $\beta$) or repulsion ($\chi_0<0$) prevents aggregation and promotes relaxation. In Part III, we study bifurcation and pattern formation when $\chi_0$ passes through critical thresholds.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper (Part II) studies persistence and stabilization for the parabolic-elliptic chemotaxis system with sensitivity χ(u,v)=χ₀ u^m/(1+v)^β and logistic source au-bu^{1+α}. Building on global existence/boundedness from Part I, it proves uniform persistence when m≥1, identifies an explicit linear stability threshold χ^*(u*) for the unique positive equilibrium (u^*,v^*)=( (a/b)^{1/α}, (ν/μ)(a/b)^{γ/α} ), establishes local exponential decay for χ₀<χ^*(u^*) and instability for χ₀>χ^*(u^*), and gives conditions under which every bounded solution converges exponentially to the equilibrium. For a=b=0 it treats mass-constrained stability of constants. The proofs extend the Lyapunov functional from the m=1 case and rectangle/ODE comparison from the β=0 case.
Significance. If the claims hold, the work supplies concrete parameter thresholds quantifying how signal saturation (large β) and repulsion (χ₀<0) suppress aggregation and promote relaxation to the homogeneous state. The explicit threshold χ^*(u*) and the direct extensions of Lyapunov and comparison techniques constitute clear technical progress over the special cases treated earlier in the literature. The results are conditional on the boundedness already proved in Part I and prepare the ground for the bifurcation analysis announced for Part III.
major comments (2)
- [§4] §4 (linearization around (u^*,v^*)): the sign of the principal eigenvalue that defines χ^*(u*) is asserted to change at the threshold, but the manuscript does not display the explicit characteristic equation or the verification that the zero eigenvalue is simple; this step is load-bearing for the stability claim and should be written out in full.
- [Theorem 5.3] Theorem 5.3 (global exponential convergence): the stated conditions on m and β are sufficient but the proof sketch invokes an extension of the Lyapunov functional whose dissipation estimate is only indicated for m>1; the case m=1 requires a separate (but presumably simpler) argument that is not supplied.
minor comments (3)
- [Introduction] The definition of χ^*(u*) is given only after the statement of the main stability theorem; move the explicit formula to the introduction or to a preliminary lemma so that the threshold is visible before its use.
- [§3] Notation: the symbol χ^*(u^*) is used both for the threshold value and for the function; a clearer distinction (e.g., χ_crit(u^*)) would avoid confusion in the statements.
- [§2] Add a short remark recalling the precise boundedness result from Part I that is invoked as a standing assumption for all convergence statements.
Simulated Author's Rebuttal
We thank the referee for the thorough review and the recommendation for minor revision. We are pleased that the significance of the results is recognized. Below we address the major comments point by point, and we will make the necessary revisions to the manuscript.
read point-by-point responses
-
Referee: [§4] §4 (linearization around (u^*,v^*)): the sign of the principal eigenvalue that defines χ^*(u*) is asserted to change at the threshold, but the manuscript does not display the explicit characteristic equation or the verification that the zero eigenvalue is simple; this step is load-bearing for the stability claim and should be written out in full.
Authors: We thank the referee for this comment. Upon review, we acknowledge that the explicit form of the characteristic equation and the simplicity of the zero eigenvalue were not fully detailed in Section 4. In the revised version, we will include the derivation of the characteristic equation from the linearized system and provide a proof that the eigenvalue zero is simple at the critical value χ^*(u^*), relying on the positivity and the structure of the steady-state equations. This will strengthen the stability analysis. revision: yes
-
Referee: [Theorem 5.3] Theorem 5.3 (global exponential convergence): the stated conditions on m and β are sufficient but the proof sketch invokes an extension of the Lyapunov functional whose dissipation estimate is only indicated for m>1; the case m=1 requires a separate (but presumably simpler) argument that is not supplied.
Authors: We appreciate the referee's careful examination of the proof of Theorem 5.3. The extension of the Lyapunov functional is presented for m > 1, but for m = 1 the functional simplifies, and the dissipation estimate follows from integrating by parts and using the logistic source terms. We will add a dedicated subsection or paragraph detailing the m = 1 case to ensure the proof is complete for all m ≥ 1 as stated. revision: yes
Circularity Check
Minor self-citation in method extensions; core linearization and equilibrium analysis independent
full rationale
The paper computes the unique positive equilibrium explicitly as (u^*,v^*) = ((a/b)^{1/α}, (ν/μ)(a/b)^{γ/α}) and obtains the stability threshold χ^*(u^*) via standard linearization of the parabolic-elliptic system, yielding the Jacobian eigenvalues without any fitted parameters or self-referential definitions. Extensions of the Lyapunov functional (from m=1) and rectangle/ODE comparison (from β=0) are described as direct adaptations of prior techniques, but these are not load-bearing for the threshold itself and do not reduce the central claims to the paper's own inputs by construction. Global boundedness is conditioned on Part I, which is external to the present derivation. No self-definitional, fitted-prediction, or ansatz-smuggling patterns appear in the quoted claims.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Global existence and boundedness of positive solutions established in Part I
- standard math Standard linearization and spectral analysis for parabolic systems on bounded domains
Reference graph
Works this paper leans on
-
[1]
Real World Appl.49(2019), 312–330,https://doi.org/10.1016/j.nonrwa.2019.03.012
Jaewook Ahn, Kyungkeun Kang, and Jihoon Lee,Eventual smoothness and stabilization of global weak solutions in parabolic-elliptic chemotaxis systems with logarithmic sensitivity, Nonlinear Anal. Real World Appl.49(2019), 312–330,https://doi.org/10.1016/j.nonrwa.2019.03.012. MR3937416
-
[2]
Khadijeh Baghaei,Global boundedness of solutions to a chemotaxis consumption model with signal dependent motility and logistic source, C. R. Math. Acad. Sci. Paris362(2024), 1131–1146. MR4824914
work page 2024
-
[3]
Piotr Biler,Global solutions to some parabolic-elliptic systems of chemotaxis, Adv. Math. Sci. Appl.9(1999), no. 1, 347–359. MR1690388
work page 1999
-
[4]
Tobias Black,Global generalized solutions to a parabolic-elliptic Keller-Segel system with singular sensi- tivity, Discrete Contin. Dyn. Syst. Ser. S13(2020), no. 2, 119–137,https://doi.org/10.3934/dcdss. 2020007. MR4043685
-
[5]
Tobias Black, Johannes Lankeit, and Masaaki Mizukami,Stabilization in the Keller-Segel system with signal- dependent sensitivity, Appl. Anal.99(2020), no. 16, 2877–2891,https://doi.org/10.1080/00036811.2019. 1585534. MR4174428
-
[6]
Junhong Cao, Wei Wang, and Hao Yu,Asymptotic behavior of solutions to two-dimensional chemotaxis system with logistic source and singular sensitivity, J. Math. Anal. Appl.436(2016), no. 1, 382–392,https://doi.org/ 10.1016/j.jmaa.2015.11.058. MR3440100
-
[7]
Differential Equations (accepted pending revision), preprint arXiv:2512.14858 (2025)
Le Chen, Ian Ruau, and Wenxian Shen,Chemotaxis models with signal-dependent sensitivity and a logistic- type source, I: Boundedness and global existence, J. Differential Equations (accepted pending revision), preprint arXiv:2512.14858 (2025)
-
[8]
,Chemotaxis models with signal-dependent sensitivity and a logistic-type source, III: Bifurcation, 2026. In preparation
work page 2026
-
[9]
Mourad Choulli, El Maati Ouhabaz, and Masahiro Yamamoto,Stable determination of a semilinear term in a parabolic equation, Commun. Pure Appl. Anal.5(2006), no. 3, 447–462,https://doi.org/10.3934/cpaa.2006. 5.447. MR2217590
-
[10]
Real World Appl.49(2019), 286–311,https: //doi.org/10.1016/j.nonrwa.2019.03.009
Mengyao Ding, Wei Wang, and Shulin Zhou,Global existence of solutions to a fully parabolic chemotaxis system with singular sensitivity and logistic source, Nonlinear Anal. Real World Appl.49(2019), 286–311,https: //doi.org/10.1016/j.nonrwa.2019.03.009. MR3936798
- [11]
-
[12]
,Partial differential equations of parabolic type, Prentice-Hall, Inc., Englewood Cliffs, N.J.,
-
[13]
Kentarou Fujie,Boundedness in a fully parabolic chemotaxis system with singular sensitivity, J. Math. Anal. Appl.424(2015), no. 1, 675–684,https://doi.org/10.1016/j.jmaa.2014.11.045. MR3286587
-
[14]
Kentarou Fujie and Takasi Senba,Global existence and boundedness in a parabolic-elliptic Keller-Segel system with general sensitivity, Discrete Contin. Dyn. Syst. Ser. B21(2016), no. 1, 81–102,https://doi.org/10.3934/ dcdsb.2016.21.81. MR3426833
work page 2016
-
[15]
Kentarou Fujie, Michael Winkler, and Tomomi Yokota,Blow-up prevention by logistic sources in a parabolic- elliptic Keller-Segel system with singular sensitivity, Nonlinear Anal.109(2014), 56–71,https://doi.org/10. 1016/j.na.2014.06.017. MR3247293
work page 2014
-
[16]
,Boundedness of solutions to parabolic-elliptic Keller-Segel systems with signal-dependent sensitivity, Math. Methods Appl. Sci.38(2015), no. 6, 1212–1224,https://doi.org/10.1002/mma.3149. MR3338145
-
[17]
Ignacio Tello,On a parabolic-elliptic system with chemotaxis and logistic type growth, J
Evgeny Galakhov, Olga Salieva, and J. Ignacio Tello,On a parabolic-elliptic system with chemotaxis and logistic type growth, J. Differential Equations261(2016), no. 8, 4631–4647,https://doi.org/10.1016/j.jde.2016. 07.008. MR3537839
-
[18]
David Gilbarg and Neil S. Trudinger,Elliptic partial differential equations of second order, Classics in Mathe- matics, Springer-Verlag, Berlin, 2001. Reprint of the 1998 edition. MR1814364
work page 2001
-
[19]
Xiao He and Sining Zheng,Convergence rate estimates of solutions in a higher dimensional chemotaxis system with logistic source, J. Math. Anal. Appl.436(2016), no. 2, 970–982,https://doi.org/10.1016/j.jmaa.2015. 12.058. MR3446989
-
[20]
840, Springer-Verlag, Berlin-New York, 1981
Daniel Henry,Geometric theory of semilinear parabolic equations, Lecture Notes in Mathematics, vol. 840, Springer-Verlag, Berlin-New York, 1981. MR610244. 49
work page 1981
-
[21]
M. A. Herrero, E. Medina, and J. J. L. Vel´ azquez,Finite-time aggregation into a single point in a reaction- diffusion system, Nonlinearity10(1997), no. 6, 1739–1754,https://doi.org/10.1088/0951-7715/10/6/
-
[22]
Miguel A. Herrero and Juan J. L. Vel´ azquez,Singularity patterns in a chemotaxis model, Math. Ann.306(1996), no. 3, 583–623,https://doi.org/10.1007/BF01445268. MR1415081
-
[23]
Liang Hong, Miaoqing Tian, and Sining Zheng,An attraction-repulsion chemotaxis system with nonlinear productions, J. Math. Anal. Appl.484(2020), no. 1, 123703, 8,https://doi.org/10.1016/j.jmaa.2019. 123703. MR4042624
-
[24]
Bingran Hu and Youshan Tao,Boundedness in a parabolic-elliptic chemotaxis-growth system under a crit- ical parameter condition, Appl. Math. Lett.64(2017), 1–7,https://doi.org/10.1016/j.aml.2016.08
-
[25]
Evelyn F. Keller and Lee A. Segel,Initiation of slime mold aggregation viewed as an instability, J. Theoret. Biol. 26(1970), no. 3, 399–415,https://doi.org/10.1016/0022-5193(70)90092-5. MR3925816
-
[26]
Evelyn F Keller and Lee A Segel,Model for chemotaxis, J. Theoret. Biol.30(1971), no. 2, 225–234
work page 1971
-
[27]
,Traveling bands of chemotactic bacteria: a theoretical analysis, J. Theoret. Biol.30(1971), no. 2, 235– 248
work page 1971
-
[28]
Real World Appl.69(2023), Paper No
Halil Ibrahim Kurt and Wenxian Shen,Chemotaxis systems with singular sensitivity and logistic source: bound- edness, persistence, absorbing set, and entire solutions, Nonlinear Anal. Real World Appl.69(2023), Paper No. 103762, 27,https://doi.org/10.1016/j.nonrwa.2022.103762. MR4489287
-
[29]
Halil ibrahim Kurt and Wenxian Shen,Stabilization in two-species chemotaxis systems with singular sensitivity and Lotka-Volterra competitive kinetics, Discrete Contin. Dyn. Syst.44(2024), no. 4, 882–904,https://doi. org/10.3934/dcds.2023130. MR4696144
-
[30]
Halil Ibrahim Kurt, Wenxian Shen, and Shuwen Xue,Stability, bifurcation and spikes of stationary solutions in a chemotaxis system with singular sensitivity and logistic source, Math. Models Methods Appl. Sci.34(2024), no. 9, 1797–1845,https://doi.org/10.1142/S0218202524500386
-
[31]
O. A. Ladyˇ zenskaja, V. A. Solonnikov, and N. N. Ural′ceva,Linear and quasilinear equations of parabolic type, Translations of Mathematical Monographs, Vol. 23, American Mathematical Society, Providence, R.I., 1968. Translated from the Russian by S. Smith. MR0241822
work page 1968
-
[33]
Gary M. Lieberman,Second order parabolic differential equations, World Scientific Publishing Co., Inc., River Edge, NJ, 1996. MR1465184
work page 1996
-
[34]
Masaaki Mizukami and Tomomi Yokota,A unified method for boundedness in fully parabolic chemotaxis systems with signal-dependent sensitivity, Math. Nachr.290(2017), no. 16, 2648–2660,https://doi.org/10.1002/mana. 201600399. MR3722501
-
[35]
Toshitaka Nagai,Blowup of nonradial solutions to parabolic-elliptic systems modeling chemotaxis in two-dimensional domains, J. Inequal. Appl.6(2001), no. 1, 37–55,https://doi.org/10.1155/ S1025583401000042. MR1887324
work page 2001
-
[36]
Toshitaka Nagai and Takasi Senba,Global existence and blow-up of radial solutions to a parabolic-elliptic system of chemotaxis, Adv. Math. Sci. Appl.8(1998), no. 1, 145–156. MR1623326
work page 1998
-
[37]
A. Pazy,Semigroups of linear operators and applications to partial differential equations, Applied Mathematical Sciences, vol. 44, Springer-Verlag, New York, 1983. MR710486
work page 1983
-
[38]
Ignacio Tello and Michael Winkler,A chemotaxis system with logistic source, Comm
J. Ignacio Tello and Michael Winkler,A chemotaxis system with logistic source, Comm. Partial Differential Equations32(2007), no. 4-6, 849–877,https://doi.org/10.1080/03605300701319003. MR2334836
-
[39]
Giuseppe Viglialoro and Thomas E. Woolley,Boundedness in a parabolic-elliptic chemotaxis system with non- linear diffusion and sensitivity and logistic source, Math. Methods Appl. Sci.41(2018), no. 5, 1809–1824, https://doi.org/10.1002/mma.4707. MR3778089
-
[40]
Chang-Jian Wang, Ya-Jie Zhu, and Xin-Cai Zhu,Long time behavior of the solution to a chemotaxis system with nonlinear indirect signal production and logistic source, Electron. J. Qual. Theory Differ. Equ. (2023), Paper No. 11, 21. MR4577527
work page 2023
-
[41]
Michael Winkler,Absence of collapse in a parabolic chemotaxis system with signal-dependent sensitivity, Math. Nachr.283(2010), no. 11, 1664–1673,https://doi.org/10.1002/mana.200810838. MR2759803. 50
-
[42]
global diffusive behavior in the higher-dimensional Keller-Segel model, J
,Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model, J. Differential Equations248(2010), no. 12, 2889–2905,https://doi.org/10.1016/j.jde.2010.02.008. MR2644137
-
[43]
,Global solutions in a fully parabolic chemotaxis system with singular sensitivity, Math. Methods Appl. Sci.34(2011), no. 2, 176–190,https://doi.org/10.1002/mma.1346. MR2778870
-
[44]
11 Under review as a conference paper at ICLR 2026 Ricky T
,Global asymptotic stability of constant equilibria in a fully parabolic chemotaxis system with strong logistic dampening, J. Differential Equations257(2014), no. 4, 1056–1077,https://doi.org/10.1016/j.jde. 2014.04.023. MR3210023
-
[45]
Michael Winkler and Tomomi Yokota,Stabilization in the logarithmic Keller-Segel system, Nonlinear Anal.170 (2018), 123–141,https://doi.org/10.1016/j.na.2018.01.002. MR3765558
-
[46]
Tian Xiang,Dynamics in a parabolic-elliptic chemotaxis system with growth source and nonlinear secretion, Commun. Pure Appl. Anal.18(2019), no. 1, 255–284,https://doi.org/10.3934/cpaa.2019014. MR3845565
-
[47]
Real World Appl.42(2018), 120–139,https: //doi.org/10.1016/j.nonrwa.2017.12.007
Xiangdong Zhao and Sining Zheng,Global existence and asymptotic behavior to a chemotaxis-consumption system with singular sensitivity and logistic source, Nonlinear Anal. Real World Appl.42(2018), 120–139,https: //doi.org/10.1016/j.nonrwa.2017.12.007. MR3773354
-
[48]
Pan Zheng,On a parabolic-elliptic Keller-Segel system with nonlinear signal production and nonlocal growth term, Dyn. Partial Differ. Equ.21(2024), no. 1, 61–76,https://doi.org/10.4310/dpde.2024.v21.n1. a3. MR4672992
-
[49]
Pan Zheng, Chunlai Mu, Robert Willie, and Xuegang Hu,Global asymptotic stability of steady states in a chemotaxis-growth system with singular sensitivity, Comput. Math. Appl.75(2018), no. 5, 1667–1675,https: //doi.org/10.1016/j.camwa.2017.11.032. MR3766542. 51
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.