1. Introduction
We consider a Keller–Segel type chemotaxis model with logarithmic sensitivity and logistic growth:
Here, the unknown functions and are the concentration of a chemical signal and the density of a cellular population, respectively. The system parameters are interpreted as follows.
is the diffusion coefficient of chemical signal.
is the coefficient of density-dependent production/degradation rate of chemical signal.
is the natural degradation rate of chemical signal.
is the coefficient of chemotactic sensitivity.
is the diffusion coefficient of cellular population.
is the natural growth rate of cellular population.
is the typical carrying capacity of cellular population.
The system describes the dynamics when certain biological organism releases or consumes a chemical signal in the local environment while both entities are naturally diffusing and reacting. It includes logarithmic chemotactic response of cells to the signal, and some or all of the following mechanisms: random walk/diffusion, consumption/deposition of the chemical by cells, natural degradation of the chemical, and the logistic growth of cells.
Biologically, the sign of
indicates whether the chemotactic movement is attractive (
) or repulsive (
). When
and
, Equation (
1) describes the movement of cells that are attracted to and consume the chemical, say, for nutrition. When
and
, as adopted in [
1] for the non-growth model, it describes the movement of cells that deposit a chemical signal to modify the local environment for succeeding passages. Such a scenario has found applications in cancer research [
2]. Since there is no difference in the analysis of these two scenarios, we assume
throughout this paper. Mathematically, the non-diffusive part of the transformed system to be discussed below is hyperbolic in biologically relevant regimes when
, while it may change type when
[
3].
The logarithmic singularity in Equation (
1) accounts for Fechner’s law, which states that subjective sensation is proportional to the logarithm of the stimulus intensity [
4]. It can be removed via the inverse Hopf–Cole transformation [
5]:
Under the variables
v and
u, Equation (
1) is converted into
Equation (
3) can be further simplified by rescaling and/or non-dimensionalization:
After dropping the tilde accent, we arrive at
where
We consider the Cauchy problem of Equation (
1):
or equivalently, the Cauchy problem of Equation (
5):
where the Cauchy datum
is assumed to be a small perturbation of a constant equilibrium state
. To be an equilibrium state, we need
or
. It is clear that the former is unstable. Therefore, we set
. To discuss
, we apply Equation (
2) to have
where for simplicity we have omitted the scaling constant
from Equation (
4). If
while
, we have
Therefore, from Equation (
9) we have
either as
or as
, depending on
or
. For physically interesting problems, we consider
with
.
Therefore, we take
. In summary,
Cauchy problem of Equations (
5) and (
8) has unique global-in-time small data solution, i.e., when
is a small perturbation of
, see [
6,
7]. To study small data solutions, especially their long time behavior, one needs to study the corresponding linear system, linearized around the constant equilibrium state. For this, we introduce new variables for the perturbation:
Linearizing Equation (
5) around
, we have
where
are constant parameters.
The goal of this paper is to obtain an accurate and detailed pointwise description, both in
x and in
t, of the Green’s function of Equation (
12). The Green’s function provides a complete solution picture to Equation (
12) and is significant in the linear theory. As discussed above, it also sheds light on the behavior of small data solutions for Equations (
5) and (
8), which will be studied in a future work.
2. Main Results and Discussion
To obtain the Green’s function, we write Equation (
12) in vector form:
where
Here,
are constants. We assume that at least one of them is positive. Otherwise, Equation (
13) has no dissipation, and its Green’s function consists of
-functions along the characteristic lines, a different scenario to what we discuss below.
The Green’s Function of Equation (
13) is the solution matrix
of
where
is the Dirac
-function, and
is the
identity matrix. Our main results on
G are the following theorems, concerning three different cases:
;
while
; and
while at least one of
and
D is positive. The cases correspond to different types of systems: hyperbolic–parabolic conservation laws, hyperbolic balance laws, and hyperbolic–parabolic balance laws.
2.1. Hyperbolic–Parabolic Conservation Laws
Theorem 1. Let , , and at least one of ε and D be positive. Let be an integer. Then, for , , the Green’s function of Equation (
13)
has the following estimates: - 1.
When ,where is a constant. - 2.
When and ,where is a constant, and , , is a , symmetric, polynomial matrix in t with a degree not more than j. In particular, - 3.
When and ,where is a constant, and , , is a , symmetric, polynomial matrix in t with a degree not more than j. In particular,
Under the assumption
, Equation (
13) becomes
Green’s function estimates on a general system in the form of Equation (
19) are detailed in [
8] (see Theorems 6.2 and 6.15 therein). It is straightforward to verify that the assumptions of those theorems are satisfied when
A and
B are given in Equation (
14). Therefore, by direct calculation and straightforward application of those theorems, we obtain Theorem 1. We note that Equations (
16)–(
18) are precise and explicit in the leading terms (and in the singular terms if
). We also note that
G is symmetric since
A and
B are, so are
in Equations (
17) and (
18).
2.2. Hyperbolic Balance Laws
Theorem 2. Let , , and be an integer. Then, for , , the Green’s function of Equation (13) has the following estimate:where is a constant, and and , , are , symmetric, polynomial matrices in t whose degrees are not more than j. In particular, Under the assumptions of Theorem 2, Equation (
13) becomes
Green’s function estimates on a general system in the form of Equation (
21) are detailed in [
9] (see Theorem 3.6 therein). It is straightforward to verify that the assumptions of that theorem are satisfied when
A and
L are given in Equation (
14). Therefore, direct application of that theorem would gives us an estimate similar to Equation (
20). Here, our result (Equation (
20)) has slightly more details in the higher order terms, the second and third terms on the righthand side of Equation (
20). This is due to the special structure of
A and
L in Equation (
14), and is justified in
Section 3.
2.3. Hyperbolic–Parabolic Balance Laws
Theorem 3. Let , , and at least one of ε and D be positive. Let be an integer. Then, for , , the Green’s function of Equation (
13)
has the following estimates: - 1.
When ,where is a constant. - 2.
When and ,where is a constant, and , , is a , symmetric, polynomial matrix in t with a degree not more than . In particular, - 3.
When and ,where is a constant, and , , is a , symmetric, polynomial matrix in t with a degree not more than . In particular,
Comparing Theorems 1–3 we observe that the solution behavior for is very different to that for . When , from Theorem 1, we see that the leading term in time decay is two heat kernels along the characteristics of A, while, for , from Theorems 2 and 3, it is a heat kernel along t-axis. Therefore, the logistic growth of cells completely changes the solution picture.
From all three theorems, we also observe that the regularity of solution depends solely on the number of nonzero diffusion coefficients and D. If both are positive, there is no -functions in the Green’s function (see Theorems 1 and 3, Case 1). If one of them is zero, then there is a -function (and its derivatives as appropriate) (see Theorems 1 and 3, Cases 2 and 3). If both are zero, then there are two -functions (see Theorem 2).
The last comment is on the role of D. If there is no logistic growth of cells, the two diffusion coefficients and D play the same role (see Theorem 1). However, if there is logistic growth, , then only r and but not D appear in the leading heat kernel (see Theorem 3). That is, logistic growth of cells overwhelms their diffusion.
In next section, we prove Theorem 3 and justify Theorem 2 to finish this paper.