Next Article in Journal
An Analytical Expression for Magnet Shape Optimization in Surface-Mounted Permanent Magnet Machines
Next Article in Special Issue
Iterated Petrov–Galerkin Method with Regular Pairs for Solving Fredholm Integral Equations of the Second Kind
Previous Article in Journal
Differential Evolution Algorithm for Multilevel Assignment Problem: A Case Study in Chicken Transportation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Green’s Function of the Linearized Logarithmic Keller–Segel–Fisher/KPP System

Department of Mathematics, University of Alabama at Birmingham, Birmingham, AL 35294, USA
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2018, 23(4), 56; https://doi.org/10.3390/mca23040056
Submission received: 14 September 2018 / Revised: 1 October 2018 / Accepted: 1 October 2018 / Published: 3 October 2018

Abstract

:
We consider a Keller–Segel type chemotaxis model with logarithmic sensitivity and logistic growth. The logarithmic singularity in the system is removed via the inverse Hopf–Cole transformation. We then linearize the system around a constant equilibrium state, and obtain a detailed, pointwise description of the Green’s function. The result provides a complete solution picture for the linear problem. It also helps to shed light on small solutions of the nonlinear system.

1. Introduction

We consider a Keller–Segel type chemotaxis model with logarithmic sensitivity and logistic growth:
c t = ε c x x - μ u c - σ c , u t + χ [ u ( ln c ) x ] x = D u x x + a u ( 1 - u K ) , x R , t > 0 .
Here, the unknown functions c = c ( x , t ) and u = u ( x , t ) are the concentration of a chemical signal and the density of a cellular population, respectively. The system parameters are interpreted as follows.
  • ε 0 is the diffusion coefficient of chemical signal.
  • μ 0 is the coefficient of density-dependent production/degradation rate of chemical signal.
  • σ 0 is the natural degradation rate of chemical signal.
  • χ 0 is the coefficient of chemotactic sensitivity.
  • D 0 is the diffusion coefficient of cellular population.
  • a 0 is the natural growth rate of cellular population.
  • K > 0 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 ( χ > 0 ) or repulsive ( χ < 0 ). When χ > 0 and μ > 0 , Equation (1) describes the movement of cells that are attracted to and consume the chemical, say, for nutrition. When χ < 0 and μ < 0 , 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 χ μ > 0 throughout this paper. Mathematically, the non-diffusive part of the transformed system to be discussed below is hyperbolic in biologically relevant regimes when χ μ > 0 , while it may change type when χ μ < 0 [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]:
v = ( ln c ) x = c x c .
Under the variables v and u, Equation (1) is converted into
v t + ( μ u - ε v 2 ) x = ε v x x , u t + χ ( u v ) x = D u x x + a u ( 1 - u K ) .
Equation (3) can be further simplified by rescaling and/or non-dimensionalization:
t ˜ = χ μ K t , x ˜ = χ μ K x , v ˜ = sign ( χ ) χ μ K v , u ˜ = u K .
After dropping the tilde accent, we arrive at
v t + ( u - ε v 2 / χ ) x = ε v x x , u t + ( u v ) x = D u x x + r u ( 1 - u ) ,
where
r = a χ μ K 0 .
We consider the Cauchy problem of Equation (1):
( c , u ) ( x , 0 ) = ( c 0 , u 0 ) ( x ) ,
or equivalently, the Cauchy problem of Equation (5):
( v , u ) ( x , 0 ) = ( v 0 , u 0 ) ( x ) ,
where the Cauchy datum ( v 0 , u 0 ) is assumed to be a small perturbation of a constant equilibrium state ( v ¯ , u ¯ ) . To be an equilibrium state, we need u ¯ = 0 or u ¯ = 1 . It is clear that the former is unstable. Therefore, we set u ¯ = 1 . To discuss v ¯ , we apply Equation (2) to have
v 0 = c 0 c 0 , c 0 ( x ) = c 0 ( 0 ) e 0 x v 0 ( y ) d y   with   c 0 ( 0 ) > 0 ,
where for simplicity we have omitted the scaling constant sign ( χ ) χ / μ K from Equation (4). If v 0 - v ¯ L 1 ( R ) while v ¯ 0 , we have
0 v 0 ( y ) d y = ± , - 0 v 0 ( y ) d y = ± .
Therefore, from Equation (9) we have c 0 ( x ) either as x or as x - , depending on v ¯ > 0 or v ¯ < 0 . For physically interesting problems, we consider lim x ± c 0 ( x ) = c ± with 0 < c ± < .
Therefore, we take v ¯ = 0 . In summary,
lim x ± ( v 0 , u 0 ) = ( v ¯ , u ¯ ) = ( 0 , 1 ) .
Cauchy problem of Equations (5) and (8) has unique global-in-time small data solution, i.e., when ( v 0 , u 0 ) is a small perturbation of ( 0 , 1 ) , 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:
w 1 = v , w 2 = u - 1 .
Linearizing Equation (5) around ( 0 , 1 ) , we have
w 1 t + w 2 x = ε w 1 x x , w 2 t + w 1 x = D w 2 x x - r w 2 ,
where ε , D , r 0 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:
w t + A w x = B w x x + L w ,
where
w = w 1 w 2 , A = 0 1 1 0 , B = ε 0 0 D , L = 0 0 0 - r .
Here, ε , D , r 0 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 G ( x , t ) of
G t + A G x = B G x x + L G G ( x , 0 ) = δ ( x ) I 2 × 2 ,
where δ ( x ) is the Dirac δ -function, and I 2 × 2 is the 2 × 2 identity matrix. Our main results on G are the following theorems, concerning three different cases: r = 0 ; r > 0 while ε = D = 0 ; and r > 0 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 r = 0 , ε , D 0 , and at least one of ε and D be positive. Let l 0 be an integer. Then, for x R , t > 0 , the Green’s function G ( x , t ) of Equation (13) has the following estimates:
1. 
When ε , D > 0 ,
l x l G ( x , t ) = l x l 1 2 π ( ε + D ) t e - ( x + t ) 2 2 ( ε + D ) t 1 2 - 1 2 - 1 2 1 2 + 1 2 π ( ε + D ) t e - ( x - t ) 2 2 ( ε + D ) t 1 2 1 2 1 2 1 2 + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - ( x + t ) 2 C t + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - ( x - t ) 2 C t ,
where C > 0 is a constant.
2. 
When ε = 0 and D > 0 ,
l x l G ( x , t ) = l x l 1 2 π D t e - ( x + t ) 2 2 D t 1 2 - 1 2 - 1 2 1 2 + 1 2 π D t e - ( x - t ) 2 2 D t 1 2 1 2 1 2 1 2 + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - ( x + t ) 2 C t + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - ( x - t ) 2 C t + e - t / D j = 0 l δ ( l - j ) ( x ) Q j ,
where C > 0 is a constant, and Q j , 0 j l , is a 2 × 2 , symmetric, polynomial matrix in t with a degree not more than j. In particular,
Q 0 = 1 0 0 0 .
3. 
When ε > 0 and D = 0 ,
l x l G ( x , t ) = l x l 1 2 π ε t e - ( x + t ) 2 2 ε t 1 2 - 1 2 - 1 2 1 2 + 1 2 π ε t e - ( x - t ) 2 2 ε t 1 2 1 2 1 2 1 2 + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - ( x + t ) 2 C t + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - ( x - t ) 2 C t + e - t / ε j = 0 l δ ( l - j ) ( x ) Q j ,
where C > 0 is a constant, and Q j , 0 j l , is a 2 × 2 , symmetric, polynomial matrix in t with a degree not more than j. In particular,
Q 0 = 0 0 0 1 .
Under the assumption r = 0 , Equation (13) becomes
w t + A w x = B w x x .
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 l = 0 ). We also note that G is symmetric since A and B are, so are Q j in Equations (17) and (18).

2.2. Hyperbolic Balance Laws

Theorem 2.
Let r > 0 , ε = D = 0 , and l 0 be an integer. Then, for x R , t > 0 , the Green’s function G ( x , t ) of Equation (13) has the following estimate:
l x l G ( x , t ) = l x l 1 4 π t / r e - r x 2 4 t 1 0 0 0 + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - x 2 C t 0 1 1 0 + ( t + 1 ) - 1 t - l + 1 2 e - x 2 C t O ( 1 ) 0 0 O ( 1 ) + e - r 2 t j = 0 l δ ( l - j ) ( x + t ) Q 1 j + e - r 2 t j = 0 l δ ( l - j ) ( x - t ) Q 2 j ,
where C > 0 is a constant, and Q 1 j and Q 2 j , 0 j l , are 2 × 2 , symmetric, polynomial matrices in t whose degrees are not more than j. In particular,
Q 10 = 1 2 - 1 2 - 1 2 1 2 , Q 20 = 1 2 1 2 1 2 1 2 .
Under the assumptions of Theorem 2, Equation (13) becomes
w t + A w x = L w .
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 r > 0 , ε , D 0 , and at least one of ε and D be positive. Let l > 0 be an integer. Then, for x R , t > 0 , the Green’s function G ( x , t ) of Equation (13) has the following estimates:
1. 
When ε , D > 0 ,
l x l G ( x , t ) = l x l 1 4 π ( ε + 1 / r ) t e - x 2 4 ( ε + 1 / r ) t 1 0 0 0 + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - x 2 C t 0 1 1 0 + ( t + 1 ) - 1 t - l + 1 2 e - x 2 C t O ( 1 ) 0 0 O ( 1 ) ,
where C > 0 is a constant.
2. 
When ε = 0 and D > 0 ,
l x l G ( x , t ) = l x l 1 4 π t / r e - r x 2 4 t 1 0 0 0 + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - x 2 C t 0 1 1 0 + ( t + 1 ) - 1 t - l + 1 2 e - x 2 C t O ( 1 ) 0 0 O ( 1 ) + e - t D j = 0 l δ ( l - j ) ( x ) Q j ,
where C > 0 is a constant, and Q j , 0 j l , is a 2 × 2 , symmetric, polynomial matrix in t with a degree not more than j / 2 . In particular,
Q 0 = 1 0 0 0 .
3. 
When ε > 0 and D = 0 ,
l x l G ( x , t ) = l x l 1 4 π ( ε + 1 / r ) t e - x 2 4 ( ε + 1 / r ) t 1 0 0 0 + O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - x 2 C t 0 1 1 0 + ( t + 1 ) - 1 t - l + 1 2 e - x 2 C t O ( 1 ) 0 0 O ( 1 ) + e - ( r + 1 ε ) t j = 0 l δ ( l - j ) ( x ) Q j ,
where C > 0 is a constant, and Q j , 0 j l , is a 2 × 2 , symmetric, polynomial matrix in t with a degree not more than j / 2 . In particular,
Q 0 = 0 0 0 1 .
Comparing Theorems 1–3 we observe that the solution behavior for r > 0 is very different to that for r = 0 . When r = 0 , from Theorem 1, we see that the leading term in time decay is two heat kernels along the characteristics of A, while, for r > 0 , 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, r > 0 , then only r and ε > 0 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.

3. Green’s Function Estimates

Notation 1.
Throughout this paper, C denotes a universal positive constant, whose value may vary line by line according to the context.
To study a linear system, we perform Fourier transform with respect to x:
w ^ ( ξ , t ) = R w ( x , t ) e - i x ξ d x , w ( x , t ) = 1 2 π R w ^ ( ξ , t ) e i x ξ d ξ .
Taking Fourier transform of Equation (15), we have
G ^ t = E G ^ , G ^ ( ξ , 0 ) = I 2 × 2 ,
where
E = E ( i ξ ) = L - i ξ A - ξ 2 B .
Solving Equation (26) gives us
G ^ ( ξ , t ) = e t E ( i ξ ) .
Applying the inverse transform, we arrive at
l x l G ( x , t ) = 1 2 π R ( i ξ ) l G ^ ( ξ , t ) e i x ξ d ξ
for an integer l 0 . Our goal in this section is to estimate the righthand side of Equation (29) to obtain the results in Theorems 2 and 3 under the assumption r > 0 .

3.1. Spectral Analysis

We carry out spectral analysis of E ( i ξ ) defined in Equation (27). By straightforward calculation, the eigenvalues of E ( i ξ ) are
λ 1 , 2 ( i ξ ) = - 1 2 [ r + ( ε + D ) ξ 2 ] ± 1 4 [ r + ( ε + D ) ξ 2 ] 2 - ξ 2 [ ( ε r + 1 ) + D ε ξ 2 ] ,
and the corresponding eigenprojections are
P 1 , 2 ( i ξ ) = 1 - ξ 2 + ( λ 1 , 2 + ε ξ 2 ) 2 - ξ 2 - i ξ ( λ 1 , 2 + ε ξ 2 ) - i ξ ( λ 1 , 2 + ε ξ 2 ) ( λ 1 , 2 + ε ξ 2 ) 2 .
Therefore, with Equation (28), we have
G ^ ( ξ , t ) = e t E ( i ξ ) = e λ 1 ( i ξ ) t P 1 ( i ξ ) + e λ 2 ( i ξ ) t P 2 ( i ξ ) ,
with λ 1 , 2 and P 1 , 2 given in Equations (30) and (31).
The leading term in e t E ( i ξ ) comes from small ξ . Thus, we consider Taylor expansions for | ξ | 1 :
λ 1 ( i ξ ) = - ( ε + 1 r ) ξ 2 + O ( ξ 4 ) , λ 2 ( i ξ ) = - r + O ( ξ 2 ) , P 1 ( i ξ ) = 1 + O ( ξ 2 ) - i ξ r + O ( ξ 3 ) - i ξ r + O ( ξ 3 ) O ( ξ 2 ) , P 2 ( i ξ ) = O ( ξ 2 ) O ( ξ ) O ( ξ ) 1 + O ( ξ 2 ) .
Similarly, the regularity of G and its derivatives comes from the expansions as | ξ | . To simplify our formulation, λ 1 takes the positive square root in Equation (30) if ε > D , and the negative one if ε < D , while λ 2 is the other one. For ε D , we have
λ 1 ( i ξ ) = - D ξ 2 + λ 10 ( i ξ ) , λ 2 ( i ξ ) = - ε ξ 2 + λ 20 ( i ξ ) , P 1 ( i ξ ) = p 1 p 2 ε - D ( i ξ ) - 1 p 2 ε - D ( i ξ ) - 1 1 - p 1 , P 2 ( i ξ ) = 1 - p 1 - p 2 ε - D ( i ξ ) - 1 - p 2 ε - D ( i ξ ) - 1 p 1 ,
where
λ 10 ( i ξ ) = - ( r + 1 ε - D ) + j = 1 c 1 j ξ - 2 j , λ 20 ( i ξ ) = 1 ε - D + j = 1 c 2 j ξ - 2 j , p 1 = p 1 ( i ξ ) = j = 1 c ˜ 1 j ξ - 2 j , p 2 = p 2 ( i ξ ) = 1 + j = 1 c ˜ 2 j ξ - 2 j
are analytic at ∞, with real coefficients c 1 j , c 2 j , c ˜ 1 j and c ˜ 2 j in Equation (35).
If ε = D , on the other hand, we have
λ 1 , 2 ( i ξ ) = - D ξ 2 ± i ξ - 1 2 r ± λ - 1 ( i ξ ) , P 1 , 2 ( i ξ ) = 1 2 ± p 3 1 2 ± p 4 1 2 ± p 4 1 2 p 3 ,
where
λ - 1 ( i ξ ) = r 2 8 ( i ξ ) - 1 + ( i ξ ) - 1 j = 1 c 3 j ξ - 2 j , p 3 ( i ξ ) = r 4 ( i ξ ) - 1 + j = 2 c ˜ 3 j ( i ξ ) - j , p 4 ( i ξ ) = j = 2 c ˜ 4 j ( i ξ ) - j
are analytic at ∞, and c 3 j , c ˜ 3 j and c ˜ 4 j are real coefficients.

3.2. Estimates on Inverse Transform

To estimate Equation (29), we focus on the case r > 0 , ε = 0 and D > 0 . All other cases are similar, and are discussed at the end of the section. Our goal is to obtain Equation (23). For this, we apply Equations (34) and (35) to have
( i ξ ) l e λ 1 ( i ξ ) t P 1 ( i ξ ) = O ( 1 ) ( i ξ ) l e - D ξ 2 t + O ( 1 ) t , ( i ξ ) l e λ 2 ( i ξ ) t P 2 ( i ξ ) = e - t D [ j = 0 l ( i ξ ) l - j Q j + ( i ξ ) - 1 Q l + 1 + O ( 1 ) ( i ξ ) - 2 ( 1 + t + + t l + l 2 ) + O ( 1 ) ( i ξ ) - 2 - l t l + l + 2 2 e O ( 1 ) ξ - 2 t ]
as | ξ | . Here, l 0 is an integer, l = 1 if l is odd, and l = 0 if l is even. On the other hand, Q j is a 2 × 2 polynomial matrix in t with a degree not more than j / 2 , 0 j l + 1 . In particular,
Q 0 = 1 0 0 0 .
With the same Q j in Equation (38), we define
R ( l ) ( x , t ) = l x l G ( x , t ) - l x l 1 4 π t / r e - r x 2 4 t 1 0 0 0 - e - t D j = 0 l δ ( l - j ) ( x ) Q j .
To obtain Equation (23), we need to prove
R ( l ) ( x , t ) = O ( 1 ) ( t + 1 ) - 1 2 t - l + 1 2 e - x 2 C t 0 1 1 0 + ( t + 1 ) - 1 t - l + 1 2 e - x 2 C t O ( 1 ) 0 0 O ( 1 )
for a constant C > 0 . Using the inverse Fourier transform in Equation (25), we have
R ( l ) ( x , t ) = 1 2 π R R ^ ( l ) ( ξ , t ) e i x ξ d ξ , R ^ ( l ) ( ξ , t ) = ( i ξ ) l G ^ ( ξ , t ) - e - 1 r ξ 2 t 1 0 0 0 - e - t D j = 0 l ( i ξ ) - j Q j .
Lemma 1.
Let r > 0 , ε = 0 and D > 0 . For x R and t > 0 , we have
| R 12 ( l ) ( x , t ) | = | R 21 ( l ) ( x , t ) | C ( t + 1 ) - 1 2 t - l + 1 2 e - x 2 C t + C t - l + 1 2 e - t C ,
| R 11 ( l ) ( x , t ) | + | R 22 ( l ) ( x , t ) | C ( t + 1 ) - 1 t - l + 1 2 e - x 2 C t + C t - l + 1 2 e - t C ,
where R j k ( l ) denotes the ( j , k ) entry of R ( l ) , 1 j , k 2 .
Proof. 
Let n > 0 be small such that Equation (33) applies for | ξ | 2 n , and N > 0 be large such that Equation (38) applies for | ξ | N . Denote the ( k , k ) entries of G ^ and Q j as G ^ k k and Q k k ( j ) , respectively. We write
R 11 ( l ) ( x , t ) = I 1 + I 2 + I 3 , I 1 = 1 2 π | ξ | n ( i ξ ) l [ G ^ 11 ( ξ , t ) - e - 1 r ξ 2 t ] e i x ξ d ξ , I 2 = 1 2 π | ξ | N ( i ξ ) l G ^ 11 ( ξ , t ) - e - t D j = 0 l ( i ξ ) - j Q 11 ( j ) e i x ξ d ξ , I 3 = 1 2 π [ n | ξ | N ( i ξ ) l G ^ 11 ( ξ , t ) e i x ξ d ξ - | ξ | n ( i ξ ) l e - 1 r ξ 2 t e i x ξ d ξ - | ξ | N e - t D j = 0 l ( i ξ ) l - j Q 11 ( j ) e i x ξ d ξ ] ,
where all integrals are over subsets of R .
For I 1 , we apply Equations (32) and (33) to have
I 1 = 1 2 π | ξ | n ( i ξ ) l [ e λ 1 ( i ξ ) t P 11 1 ( ξ , t ) - e - 1 r ξ 2 t ] e i x ξ d ξ + O ( 1 ) e - r 2 t ,
where P 11 1 is the ( 1 , 1 ) entry of P 1 . From Equations (30) and (31), we note that the integrand in Equation (45) is holomorphic in ξ (as a complex variable) in a neighborhood of the origin. Taking n small, we apply Cauchy integral theorem to replace the domain of integration [ - n , n ] by a path Γ ( α ) { - n + i η | η  is from  α } { ζ + i α | - n ζ n } { n + i η | η  is from  α   to   0 } . Here, α can be positive or negative, but | α | n . With Equation (33) we have
I 1 = Γ ( α ) ( i ξ ) l [ e - 1 r ξ 2 t + O ( | ξ | 4 ) t [ O ( | ξ | 4 ) t + O ( | ξ | 2 ) ] e i x ξ d ξ + O ( 1 ) e - r 2 t .
If r | x | / t n , we set α = r x / ( 2 t ) . Integrating over each pieces of Γ ( α ) , we have
| I 1 | C - n n e - r x 2 8 t - ζ 2 t 2 r ( | ζ | + | x | t ) l + 2 d ζ + C e - n 2 2 r t + C e - r 2 t C ( t + 1 ) - l + 3 2 e - x 2 C t + C e - t C .
If r | x | / t > n , we set α = n 2 sign ( x ) . The straightforward calculation yields
| I 1 | C e - t C .
To estimate I 2 , we apply Equation (38) to Equation (32) to have
I 2 = O ( 1 ) | ξ | N | ξ | l e - D 2 ξ 2 t d ξ + 1 2 π e - t D | ξ | N ( i ξ ) - 1 e i x ξ d ξ Q 11 ( l + 1 ) + O ( 1 ) e - t 2 D | ξ | N ξ - 2 d ξ = O ( 1 ) t - l + 1 2 e - t C ,
noting the second integral on the right-hand side of Equation (48) is the principal value.
To estimate I 3 , we write Equation (44)
I 3 = 1 2 π n | ξ | N ( i ξ ) l G ^ 11 ( ξ , t ) e i x ξ d ξ + O ( 1 ) t - l + 1 2 e - t C .
From Equation (30), we note that, for ξ R { 0 } , the real parts of λ 1 , 2 , ( λ 1 , 2 ) , are negative. As a continuous function on a compact set,
( λ 1 , 2 ) ( i ξ ) - β , ξ [ - N , - n ] [ n , N ]
for a constant β > 0 .
From the characteristic equation of E ( i ξ ) , it is straightforward to verify
- ξ 2 + λ 1 , 2 2 = ± λ 1 , 2 ( λ 1 - λ 2 ) .
Noting ε = 0 , P 1 , 2 ( i ξ ) in Equation (31) is analytic in any domain of C where λ 1 and λ 2 are distinct. From Equation (30), there are ξ j C , 1 j 4 , such that λ 1 = λ 2 . If a ξ j is on ( - N , - n ) ( n , N ) for n > 0 small and N > 0 large, we replace ( ξ j - α , ξ j + α ) , with α > 0 small, by a semi-circle centered at ξ j with radius α . In this way, we replace [ - N , - n ] [ n , N ] by a union of two paths, denoted as Γ. Noting Equation (50) and the continuity of ( λ 1 , 2 ) , we have
( λ 1 , 2 ) ( i ξ ) - β 2 , | P 1 , 2 ( i ξ ) | C , ξ Γ ,
by choosing α small. Since the integrand of Equation (49) is an entire function, applying Cauchy theorem and Equation (32), and substituting Equation (51) into Equation (49), we arrive at
I 3 = 1 2 π Γ ( i ξ ) l G ^ 11 ( ξ , t ) e i x ξ d ξ + O ( 1 ) t - l + 1 2 e - t C = O ( 1 ) t - l + 1 2 e - t C .
Here, we have chosen the small semi-circles in Γ on the upper-half complex plane if x 0 , and lower half-plane if x < 0 , so that ( i x ξ ) 0 .
Combining Equations (44), (46)–(48) and (52), we obtain the estimate for R 11 ( l ) in Equation (43). The estimates for R 22 ( l ) and R 12 ( l ) are obtained in the same way. In particular, the slower decay rate in R 12 ( l ) comes from - i ξ / r in the ( 1 , 2 ) entry of P 1 ( i ξ ) in Equation (33), comparing O ( ξ 2 ) in the ( 2 , 2 ) entry. Finally, Equations (14) and (27)–(29) imply that G ( x , t ) is symmetric. Therefore, Q j in Equation (38) hence in Equation (23) and R ( l ) in Equation (39) are symmetric. This gives us R 12 ( l ) = R 21 ( l ) . ☐
Lemma 2.
Let K > 0 be large. Under the assumptions of Lemma 1, for ξ R , t > 0 , and | x | / t K , we have
| R ( l ) ( x , t ) | C t - l + 1 2 e - x 2 C t ,
where R ( l ) is defined in Equation (39).
Proof. 
From Equations (41) and (28), we have
R ( l ) ( x , t ) = I 4 - l x l 1 4 π t / r e - r x 2 4 t 1 0 0 0 = I 4 + O ( 1 ) t - l + 1 2 e - x 2 C t , I 4 = 1 2 π R ( i ξ ) l e t E ( i ξ ) - e - t D j = 0 l ( i ξ ) - j Q j e i x ξ d ξ .
Since the integrand of I 4 is an entire function, by a standard argument, we apply Cauchy theorem to replace the integral path R by Γ x / t = { ζ + i x 2 D t | ζ R } . Taking K large and applying Equations (32) and (38) gives us
I 4 = 1 2 π Γ x / t ( i ξ ) l e t E ( i ξ ) - e - t D j = 0 l ( i ξ ) - j Q j e i x ξ d ξ = 1 2 π Γ x / t { O ( 1 ) ( i ξ ) l e - D ξ 2 t + O ( 1 ) t + e - t D [ ( i ξ ) - 1 Q l + 1 + O ( 1 ) | ξ | - 2 ( 1 + t + + t l + l 2 ) + O ( 1 ) | ξ | - 2 - l t l + l + 2 2 e O ( 1 ) | ξ | - 2 t ] } e i x ξ d ξ = O ( 1 ) R ( | ζ | l + | x t | l ) e - D ζ 2 t - x 2 8 D t d ζ + e - t D 2 π R e i x ζ i ζ - x 2 D t d ζ e - x 2 2 D t Q l + 1 + O ( 1 ) e - t 2 D - x 2 2 D t R d ζ ζ 2 + K 2 4 D 2 .
Note that
1 i ζ - x 2 D t = 1 i ζ - x 2 D t 1 ζ 2 + ( x 2 D t ) 2 + 1 i ζ x 2 D t ζ 2 + ( x 2 D t ) 2 ,
which implies
R e i x ζ i ζ - x 2 D t d ζ = R e i x ζ i ζ d ζ + O ( 1 ) | x t | + | x t | 2 | x | .
Therefore,
I 4 = O ( 1 ) t - l + 1 2 e - x 2 C t .
Substituting Equation (55) into Equation (54) gives us Equation (53). ☐
Combining Lemmas 1 and 2, and noting that G is symmetric, we arrive at Equation (23). The proof of Equation (24) is parallel. The proof of Equation (22) is simpler since Equations (34) and (36) imply that G and its derivatives contain no δ -functions if ε , D > 0 . This settles Theorem 3.
Theorem 2 can be either proved as Theorem 3, or derived from the general framework in [9], noting that G is symmetric, and that those O ( ξ 2 ) on the diagonal of P 1 ( i ξ ) in Equation (33) give an extra ( t + 1 ) - 1 2 , comparing to O ( ξ ) in the general framework.

Author Contributions

J.R. and Y.Z. contributed equally in the investigation of the project; Y.Z. wrote the paper.

Acknowledgments

The research of Yanni Zeng was partially supported by a grant from the Simons Foundation (#244905 to Yanni Zeng).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Othmer, H.; Stevens, A. Aggregation, blowup and collapse: The ABC’s of taxis in reinforced random walks. SIAM J. Appl. Math. 1997, 57, 1044–1081. [Google Scholar]
  2. Levine, H.A.; Sleeman, B.D.; Nilsen-Hamilton, M. A mathematical model for the roles of pericytes and macrophages in the initiation of angiogenesis. I. The role of protease inhibitors in preventing angiogenesis. Math. Biosci. 2000, 168, 71–115. [Google Scholar] [CrossRef]
  3. Zeng, Y.; Zhao, K. On the logarithmic Keller–Segel–Fisher/KPP model. 2018; submitted. [Google Scholar]
  4. Keller, E.F.; Segel, L.A. Traveling bands of chemotactic bacteria: A theoretical analysis. J. Theor. Biol. 1971, 30, 235–248. [Google Scholar] [CrossRef]
  5. Levine, H.A.; Sleeman, B.D. A system of reaction diffusion equations arising in the theory of reinforced random walks. SIAM J. Appl. Math. 1997, 57, 683–730. [Google Scholar]
  6. Zeng, Y. Lp time asymptotic decay for general hyperbolic-parabolic balance laws with applications. 2018; submitted. [Google Scholar]
  7. Zeng, Y. Hyperbolic-Parabolic Balance Laws: Asymptotic Behavior and a Chemotaxis Model. 2018; submitted. [Google Scholar]
  8. Liu, T.-P.; Zeng, Y. Large Time Behavior of Solutions for General Quasilinear Hyperbolic–Parabolic Systems of Conservation Laws; American Mathematical Society: Providence, RI, USA, 1997. [Google Scholar]
  9. Zeng, Y. Gas dynamics in thermal nonequilibrium and general hyperbolic systems with relaxation. Arch. Rational Mech. Anal. 1999, 150, 225–279. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Rugamba, J.; Zeng, Y. Green’s Function of the Linearized Logarithmic Keller–Segel–Fisher/KPP System. Math. Comput. Appl. 2018, 23, 56. https://doi.org/10.3390/mca23040056

AMA Style

Rugamba J, Zeng Y. Green’s Function of the Linearized Logarithmic Keller–Segel–Fisher/KPP System. Mathematical and Computational Applications. 2018; 23(4):56. https://doi.org/10.3390/mca23040056

Chicago/Turabian Style

Rugamba, Jean, and Yanni Zeng. 2018. "Green’s Function of the Linearized Logarithmic Keller–Segel–Fisher/KPP System" Mathematical and Computational Applications 23, no. 4: 56. https://doi.org/10.3390/mca23040056

APA Style

Rugamba, J., & Zeng, Y. (2018). Green’s Function of the Linearized Logarithmic Keller–Segel–Fisher/KPP System. Mathematical and Computational Applications, 23(4), 56. https://doi.org/10.3390/mca23040056

Article Metrics

Back to TopTop