![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
In this paper, we consider a flux-limited Keller–Segel model derived in [Citation16, Citation18] in a one-dimensional bounded domain and give a refined asymptotic result of the spiky steady state by using the Sturm oscillation theorem in a more meticulous way based on the existence result of spiky steady state in [Citation4], showing a more accurate characterization of the cell aggregation phenomenon.
1. Introduction
Chemotaxis is a fundamental cellular reaction for survival, which describes the directional movements towards or away from the chemical stimuli. It is a fundamental mechanism in many important biological phenomenon such as the healing of wound [Citation17], pattern formation [Citation2, Citation14] and so on. Mathematical models of chemotaxis were widely studied in recent years and the model proposed by Keller and Segel [Citation9, Citation10] was well known and received widespread acceptance. A lot of outstanding work has been done on the Keller–Segel system and its variants to study the aggregation and wave propagation of cells, two important phenomena in chemotaxis, and we refer the reader to [Citation5, Citation7, Citation8] for details. About the wave propagation phenomenon, we refer the reader to [Citation12] and [Citation19], which are the papers for the traveling waves and the periodic waves generated from the Keller–Segel model.
In this paper, we focus on the flux-limited Keller–Segel (FLKS) system, which reads as follows: (1)
(1) where
represents the density of cells and
the concentration of chemical substance at time t and position x; ψ is the function to describe the chemical response mechanism, which is so-called the chemotactic sensitivity function with the chemotactic coefficient
. ψ is assumed to be a smooth function and a typical form of ψ is
. The assumption of limited flux here means that there exists a positive constant
such that
and
. So when the cells are in the large gradient environment (
is large), the chemotactic flux,
, is bounded to model velocity saturation.
The earliest relevant work about the FLKS system is in [Citation15], where Patlak first used the kinetic theory to express the chemotactic velocity in term of the average of velocities and run times of individual cells. Later, this approach was boosted and developed by Alt et al. [Citation1, Citation13] and so on. The derivation of the kinetic chemotaxis model in previous work gives the motivation to study the FLKS system. One can see details of the derivation of the FLKS system in [Citation16, Citation18].
In the FLKS system, since the global existence in time t of the solutions of (Equation1(1)
(1) ) with flux limited assumption is proven in [Citation4, Citation6], blow-up will not happen in finite time. In [Citation16], the large-time behaviour of the FLKS system in
was studied with a fixed initial mass
, denoted as m, for
. From [Citation16], we know that the existence of positive non-constant radial steady states depends on the dimension n and the initial mass m, that is, for n = 2, radial symmetric solutions exist if and only if
while there is no positive non-constant radial steady state for any m>0 when n>2. For any bounded domain
,
, the FLKS system with Neumann boundary conditions and mass constraint has positive non-constant radial steady states when the chemotactic coefficient χ is large enough (see [Citation20]), where the mass constraint means that
which comes from the integration of (Equation1
(1)
(1) ) with Neumann boundary conditions. In [Citation20], the existence of radial spiky steady states is proved, that is, when χ tends to be infinity, the non-constant radial steady states found in [Citation20] tends to be like spikes, where finding spiky steady states is another way in mathematics to describe the aggregation phenomenon of chemotaxis models except blow-up.
In this paper, we consider the one-dimensional steady states of the FLKS system (Equation1(1)
(1) ) with the mass constraint and Neumann boundary conditions in a bounded domain as follows:
(2)
(2) where α, M and L are fixed positive constants. We denote ML as the total mass of cells and denote
for convenience. Hence
satisfying
(3)
(3) The flux limited situation
and the smoothness of the function ϕ imply that there must exist a constant
such that
. Without loss of generality, we assume
in the following for simplicity.
The existence of monotone solution of (Equation2(2)
(2) ) is proved in [Citation4] for
where
is a constant depends on L, α, M and ϕ. As χ goes to infinity, it is proved that the decreasing solution ρ concentrates at x = 0, that is
in the sense of distribution, and the decreasing solution S converges to an explicit function uniformly on
, by using Helly's compactness theorem and the Sturm oscillation theorem (see Theorem 5.4 in [Citation4]), which implies the existence of the spiky steady state and hence the aggregation of cells. In this paper, we use the Sturm oscillation theorem in a more meticulous way to give a refined asymptotic result, that is, to give a more accurate characterization of the cell aggregation.
2. A priori estimate
Here we give the priori estimate first.
Lemma 2.1
If is the positive classical solution of (Equation2
(2)
(2) ), we have
(4)
(4)
(5)
(5) and
(6)
(6)
Proof.
We use the Green's function to give the lower bound of S. For any fixed , let
be the Green's function satisfying
Observe that G can be explicitly given by
Since
,
, we have that for any
,
Now we give the upper bound of v. Recall
Since
and
, we can rewrite
as
for
. Because the solution is positive, that is,
on
, we have
where the first inequality comes from the second equality in (Equation1
(1)
(1) ). The mean value theorem of integration implies that there exists some
such that
. Similarly, S can be estimated by
By the equation in the first line of (Equation2
(2)
(2) ) and the boundary condition, we have
which implies that
Take
for some
. Then
that is
3. Refined asymptotic
Recalling in [Citation4], the results that as ,
in the sense of distribution and
uniformly on
were proved. These results are true in our paper, and we give a refined asymptotic result below with a method coming from [Citation3]. Although both the functions ρ and S depend on the parameter χ, for convenience, in the following we still use the notations
and
other than
and
.
Theorem 3.1
For the solutions of (Equation2(2)
(2) ) in Theorem 5.4 in [Citation4], we have a further result, which is
(7)
(7) and
(8)
(8)
Remark 3.1
Comparing Equation (53) in [4] and Equation (Equation2(2)
(2) ) in our paper, we denote
. Observe that if the assumptions in (54) in [Citation4] are satisfied, that is,
,
and
, then our assumptions in (Equation3
(3)
(3) ) are satisfied by the facts
and
. So the assumptions of function ϕ in our paper give more relaxed restrictions of the FLKS system (Equation2
(2)
(2) ) than the one, system (53) and (56), in [Citation4]. Actually, using the same methods in [Citation4] and our paper, we can also obtain the results that as
,
in the sense of distribution and
uniformly on
and the result in Theorem 3.1 for system (Equation2
(2)
(2) ) with ϕ defined in (Equation3
(3)
(3) ).
Remark 3.2
From (Equation7(7)
(7) ), we can see that as
, the function
tends to be 0 with the order
on
and tends to be infinity with the order χ at the point 0. So we have
for
and
as
and Theorem 3.1 gives more detailed description of the profile of
when χ is large, that is, how the profile of ρ tends to be like a spike at the point 0 as
. Denote
and
, it is easy to check that
, and we can see the change of the profile of
clearly by (Equation8
(8)
(8) ) when χ is large.
Proof.
We will use Sturm's oscillation theorem in a more meticulous way to get our result. Recall for
. Using Sturm's oscillation theorem by a comparison between
and
, and using the fact that
on
for fixed
, we obtain that
for
. By (Equation5
(5)
(5) ) and the fact that ρ is decreasing, we have
for
and hence
where
is well defined because
is
smooth. Thus we have
which implies that
(9)
(9) Recall that
. Thus with Lemma 2.1 we have
Integrating the above equation over
, we obtain
which is
Since
is uniformly bounded for all large χ, we can assume that there exists some constant B such that
after passing to a subsequence as
. Then we take
and
Thus
and
Since S and
are uniformly bounded for all
,
is uniformly bounded. So we can assume
locally uniformly on
after further passing to a subsequence as
. For any fixed point
, since
as
,
as
, which implies that
By the integration of
, we have
which is
Now we need to find the value of B. Observe that
By (Equation9
(9)
(9) ), we know
which is,
.
Since the limit is unique, we have
By the monotonicity of
, we have
(10)
(10) Then the Lebesgue's dominated convergence theorem gives us
Recall that
for some positive constant C and
. Then we have
by (Equation10
(10)
(10) ) and the fact that
converges to
as
. Thus we prove (Equation7
(7)
(7) ).
By the representation of S via Green's function, which is
then we have
With a same discussion of the estimate of
in the proof of Theorem 3.2 in [Citation11], we have (Equation8
(8)
(8) ) and we finish this proof.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
References
- W. Alt, Biased random walk model for chemotaxis and related diffusion approximation, J. Math. Biol. 9 (1980), pp. 147–177.
- E. Budrene and H. Berg, Complex patterns formed by motile cells of Escherichia coli, Nature 349 (1991), pp. 630–633.
- X. Chen, J. Hao, X. Wang, Y. Wu, and Y. Zhang, Stability of spiky solution of Keller–Segel's minimal chemotaxis model, J. Differ. Equ. 257 (2014), pp. 3102–3134.
- A. Chertock, A. Kurganov, X. Wang, and Y. Wu, On a chemotaxis model with saturated chemotactic flux, Kinet. Relat. Models 5 (2012), pp. 51–95.
- T. Hillen and K. Painter, A user's guide to PDE models for chemotaxis, J. Math. Biol. 58 (2009), pp. 183–217.
- T. Hillen, K. Painter, and C. Schmeiser, Global existence for chemotaxis with finite sampling radius, Discrete Contin. Dyn. Syst. Ser. B. 7 (2007), pp. 125–144.
- D. Horstmann, From 1970 until present: The Keller–Segel model in chemotaxis and its consequences I, Jahresber. DMV 105 (2003), pp. 103–165.
- D. Horstmann, From 1970 until present: The Keller–Segel model in chemotaxis and its consequences II, Jahresber. DMV 106 (2004), pp. 51–69.
- E. Keller and L. Segel, Initiation of slime mold aggregation viewed as an instability, J. Theor. Biol.26 (1970), pp. 399–415.
- E. Keller and L. Segel, Models for chemotaxis, J. Theor. Biol. 30 (1971), pp. 225–234.
- H. Li, Spiky steady states of a chemotaxis system with singular sensitivity, J. Dyn. Differ. Equ. 30 (2018), pp. 1775–1795.
- P. Liu, J.P. Shi, and Z.A. Wang, Pattern formation of the attraction-repulsion Keller–Segel system, Discrete Contin. Dyn. Syst. Ser. B 18 (2013), pp. 2597–2625.
- H. Othmer, S.R. Dunbar, and W. Alt, Models of dispersal in biological systems, J. Math. Biol. 26 (1988), pp. 263–298.
- K.J. Painter, P.K. Maini, and H.G. Othmer, Stripe formation in juvenile pomacanthus explained by a generalized turing mechanism with chemotaxis, Proc. Natl. Acad. Sci. 96 (1999), pp. 5549–5554.
- C.S. Patlak, Random walk with persistence and external bias, Bull. Math. Biophys. 15 (1953), pp. 311–338.
- B. Perthame, N. Vauchelet, and Z. Wang, The flux limited Keller–Segel system; properties and derivation from kinetic equations, Rev. Mat. Iberoam. 36 (2020), pp. 357–386.
- G.J. Petter, H.M. Byrne, D.L.S. Mcelwain, and J. Norbury, A model of wound healing and angiogenesis in soft tissue, Math. Biosci. 136 (1996), pp. 35–63.
- J. Saragosti, V. Calvez, N. Bournaveas, A. Buguin, P. Silberzan, and B. Perthame, Mathematical description of bacterial traveling pulses, PLoS Comput. Biol. 6 (8) (2010), Article ID e1000890.
- Z.A. Wang, Mathematics of traveling waves in chemotaxis, Discrete Contin. Dyn. Syst. Ser. B. 18 (2013), pp. 601–641.
- Z.A. Wang and X. Xu, Radial spiky steady states of a flux-limited Keller–Segel model: Existence, asymptotics and stability, Stud. Appl. Math. 148 (2022), pp. 1251–1273.