Next Article in Journal
The Finite Difference Method and Analysis for Simulating the Unsteady Generalized Maxwell Fluid with a Multi-Term Time Fractional Derivative
Previous Article in Journal
A Bearing Fault Diagnosis Method under Small Sample Conditions Based on the Fractional Order Siamese Deep Residual Shrinkage Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Minimal Wave Speed for a Nonlocal Viral Infection Dynamical Model

1
Key Laboratory of Eco-Environments in Three Gorges Reservoir Region (Ministry of Education), School of Mathematics and Statistics, Southwest University, Chongqing 400715, China
2
Shanxi Key Laboratory of Mathematical Techniques and Big Data Analysis on Disease Control and Prevention, Complex Systems Research Center, Shanxi University, Taiyuan 030006, China
*
Authors to whom correspondence should be addressed.
Submission received: 23 December 2023 / Revised: 9 February 2024 / Accepted: 13 February 2024 / Published: 26 February 2024
(This article belongs to the Topic Advances in Nonlinear Dynamics: Methods and Applications)

Abstract

:
To provide insights into the spreading speed and propagation dynamics of viruses within a host, in this paper, we investigate the traveling wave solutions and minimal wave speed for a degenerate viral infection dynamical model with a nonlocal dispersal operator and saturated incidence rate. It is found that the minimal wave speed c is the threshold that determines the existence of traveling wave solutions. The existence of traveling fronts connecting a virus-free steady state and a positive steady state with wave speed c c is established by using Schauder’s fixed-point theorem, limiting arguments, and the Lyapunov functional. The nonexistence of traveling fronts for c < c is proven by the Laplace transform. In particular, the lower-bound estimation of the traveling wave solutions is provided by adopting a rescaling method and the comparison principle, which is a crucial prerequisite for demonstrating that the traveling semifronts connect to the positive steady state at positive infinity by using the Lyapunov method and is a challenge for some nonlocal models. Moreover, simulations show that the asymptotic spreading speed may be larger than the minimal wave speed and the spread of the virus may be postponed if the diffusion ability or diffusion radius decreases. The spreading speed may be underestimated or overestimated if local dispersal is adopted.

1. Introduction

Although huge advances have been made in preventing and treating HIV and viral hepatitis, such as antiretroviral treatment for HIV and vaccination programs for the hepatitis B virus (HBV), HIV and HBV pandemics remain a major global public health problem. It is reported that there were 36.9 million people living with HIV worldwide in 2017 and 257 million people and 71 million people in 2015 were living with HBV and hepatitis C virus (HCV), respectively. Meanwhile, 0.94 million people in 2017 and 1.3 million people in 2015 died from AIDS-related and hepatitis-related illnesses, respectively [1,2]. Therefore, we have a long way to go to control and extinguish these viral infectious diseases.
To understand the pathogenesis of viruses within the host and then propose more effective control measures, many different methods have been developed. In particular, mathematical models have been verified as an effective method [3,4]. In 2000, Nowak and May [5] proposed the following basic viral infection model:
d w ( t ) d t = s b w ( t ) β w ( t ) v ( t ) , d u ( t ) d t = β w ( t ) v ( t ) μ u ( t ) , d v ( t ) d t = p u ( t ) γ v ( t ) ,
where w ( t ) , u ( t ) , and v ( t ) denote the concentrations of healthy target cells, infected cells, and free virions at time t, respectively. s is the recruitment rate of healthy target cells. b, μ , and γ represent the death rates of healthy target cells, infected cells, and free virions, respectively. The infectious incidence rate is β w ( t ) v ( t ) and p is the virus production rate. All the parameters in model (1) are positive. System (1) has a virus-free equilibrium point E 0 = ( w 0 , 0 , 0 ) , which is globally asymptotically stable when R 0 1 , where w 0 = s / b and R 0 = p β s / ( b μ γ ) = p β w 0 / ( μ γ ) . If R 0 > 1 , system (1) admits at a unique positive equilibrium point that is globally asymptotically stable [6]. Since then, many works concerning the impacts of various factors on within-host viral dynamics have been conducted using mathematical models [7,8]. Although the incidence rate in most of these viral models adopts a bilinear function response, this may be not so appropriate when the concentration of virions is high. In this case, the saturation effect may cause a viral response rate that is less than linear. Hence, it is more reasonable to adopt a saturation nonlinear incidence rate β w ( t ) v m ( t ) / ( 1 + α v n ( t ) ) , where m , n , α > 0 . The case where m = n = 1 has been studied in viral infection models by several researchers, including [9,10].
Note that many studies on viral infection models assume that the within-host environments are homogeneous, and ignore the impact of heterogeneous environments and the mobility of virions or cells. However, virions or cells may move within and between tissues and may face different environments in different locations within the host, which would consequently impact the dynamics of the virus [11,12]. Thus, it is more reasonable to incorporate spatial factors into the models, which have been studied by some researchers [13,14,15]. Strain et al. [13] introduced a lattice cellular automaton model to investigate the contribution of three-dimensional spatial correlations in viral propagation. Wang and Wang [14] proposed a degenerate HBV infection model with a local dispersal operator and investigated the existence of traveling wave solutions. Lai and Zou [15] established a reaction-diffusion viral infection model with a repulsion effect and investigated its spreading speed and the existence of traveling wave solutions. Most of these studies assume that the virions or cells diffuse in the form of local dispersal and follow Fickian diffusion, which can only be used to study situations where the density of the species is relatively low and the species diffuses in a small range [16]. However, the concentrations of virions and cells are relatively high within tissues, which suggests that nonlocal dispersal may be more reasonable in viral infection models. Moreover, the nonlocal dispersal operator can be viewed as an approximation of the local dispersal operator when the kernel function takes a special form [17]. Recently, Zhao and Ruan [18] assumed that the virions diffuse in the form of the nonlocal mode in domain Ω R n ( n 1 ) , and subsequently proposed and analyzed the following nonlocal viral infection model:
w ( x , t ) t = s b w ( x , t ) β w ( x , t ) v ( x , t ) , u ( x , t ) t = β w ( x , t ) v ( x , t ) μ u ( x , t ) , v ( x , t ) t = d v Ω J ( x y ) [ v ( y , t ) v ( x , t ) ] d y + p u ( x , t ) γ v ( x , t ) ,
where w ( x , t ) , u ( x , t ) , and v ( x , t ) are the concentrations of the target cells, infected cells, and free virions at time t and location x, respectively. d v represents the diffusion rate of the virions. Here, J ( x y ) can be viewed as the probability that virions jump from location y to location x and J ( x y ) = J ( y x ) . Thus, the nonlocal dispersal operator Ω J ( x y ) [ v ( y , t ) v ( x , t ) ] d y includes not only the rate that virions arrive at location x from other locations ( Ω J ( x y ) v ( y , t ) d y ), but also the leaving rate of virions at location x ( Ω J ( y x ) v ( x , t ) d y ). Other parameters have the same meanings as those in model (1). The authors in [18] investigated the threshold dynamics of model (2) and the impact of the dispersal rate on solutions of (2).
In the process of viral transmission, there is evidence exhibiting that virions can spread in a way like a traveling wave front [13]. Thus, if the virus diffusion takes the form of nonlocal dispersal in an unbounded domain, two interesting questions arise: (1) Can the model exhibit traveling wave solutions or not? (2) What is the spreading speed of the virus? Additionally, accurate estimates of the spreading speed, especially at the early stage of viral infection, can provide insights into how the virus propagates. From a mathematical point of view, estimates of the spreading speed can usually be obtained by studying the asymptotic spreading speed, which is relative to minimal wave speed. In this paper, inspired by the above-mentioned arguments, we intend to study the traveling wave solutions and minimal wave speed problems of the following viral infection model with the saturation incidence rate:
w ( x , t ) t = s b w ( x , t ) β w ( x , t ) v ( x , t ) 1 + α v ( x , t ) , u ( x , t ) t = β w ( x , t ) v ( x , t ) 1 + α v ( x , t ) μ u ( x , t ) , v ( x , t ) t = d v [ ( J v ) ( x , t ) v ( x , t ) ] + p u ( x , t ) γ v ( x , t ) ,
where ( J v ) ( x , t ) = R J ( x y ) v ( y , t ) d y . Here, the domain is R , and the incidence rate is in saturated mass action β w ( x , t ) v ( x , t ) / ( 1 + α v ( x , t ) ) . Other parameters are the same as those in (2). Throughout this paper, we assume that the dispersal kernel J satisfies
( H )
J C 1 ( R ) , J ( x ) = J ( x ) 0 , R J ( x ) d x = 1 , J is compactly supported and R J ( x ) e λ x d x < + for all λ > 0 .
Clearly, system (3) always admits in a virus-free steady state E 0 = ( w 0 , 0 , 0 ) , where w 0 = s / b . Moreover, the ODE system associated with system (3) admits a unique positive steady state E = ( w , u , v ) when R 0 > 1 , where
R 0 = p β s b μ γ = p β w 0 μ γ
is the basic reproduction number of the corresponding ODE system. In the rest of this paper, we always assume that R 0 > 1 holds.
The traveling wave solution of (3) is a positive solution ( w ( x , t ) , u ( x , t ) , v ( x , t ) ) of (3) which has the form
( w ( x , t ) , u ( x , t ) , v ( x , t ) ) = ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) , ξ = x + c t ,
where c > 0 is the wave speed. A positive traveling wave solution ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) is called a traveling semifront of (3) if it satisfies lim ξ ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) = E 0 , and it is called a traveling front if it satisfies
lim ξ ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) = E 0 , lim ξ + ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) = E .
It is clear that the traveling wave solution ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) satisfies
c W ( ξ ) = s b W ( ξ ) β W ( ξ ) V ( ξ ) 1 + α V ( ξ ) , c U ( ξ ) = β W ( ξ ) V ( ξ ) 1 + α V ( ξ ) μ U ( ξ ) , c V ( ξ ) = d v [ ( J V ) ( ξ ) V ( ξ ) ] + p U ( ξ ) γ V ( ξ ) ,
where ( J V ) ( ξ ) = R J ( ξ y ) V ( y ) d y .
The viral infection model is neither a cooperative system nor a competitive system, which together with the existence of the recruitment term of healthy target cells infers that the classic methods, such as the monotone semiflow method, the shooting method, and connection index theory, are all not valid. Meanwhile, as far as we know, few mathematical works have been performed to study the existence of traveling wave solutions and the minimal wave speed in viral infection models [14,19,20,21], especially for nonlocal systems. Furthermore, the nonlocal dispersal operator causes the solutions of system (3) to lack regularity and compactness, which may lead to new difficulties in analysis. Recently, Wang and Ma [22] investigated the traveling wave solution problem for a nonlocal HIV infection model with a Beddington–DeAngelis functional response, where they assumed that all cells and virions can nonlocally diffuse but have the same diffusion ability. They proved the existence of traveling wave solutions for c c , but there are some additional conditions for c = c . The existence of traveling wave solutions for c = c and the nonexistence for c < c were further studied in [23]. It is worth noting that using the Lyapunov function is an effective method to show the traveling wave solutions connect to the positive steady state. However, not only upper-bound estimations of the solution are required, but also lower-bound estimations, which is also a challenge for nonlocal systems. In particular, only free virions can diffuse in our model, which may also lead to some challenges. In this paper, we will overcome the aforementioned difficulties to obtain traveling wave solutions and the minimal wave speed of system (3) by utilizing Schauder’s fixed-point theorem, the rescaling method, the comparison principle, and so on.
The paper is organized as follows. In the next section, we establish the existence or nonexistence of traveling wave solutions with wave speed c > 0 , and give the minimal wave speed of system (3). In Section 3, we discuss the results. Some conclusions are presented in the final section.

2. Traveling Wave Solutions

In this section, we mainly focus on the traveling wave solutions and the minimal wave speed of system (3). Firstly, we establish the existence of traveling fronts for c > c . Secondly, we show the existence of traveling fronts for c = c . Finally, the nonexistence of traveling fronts is investigated for 0 < c < c .

2.1. The Existence of Traveling Fronts for c > c

In this subsection, we first give the definition of c and then study the existence of traveling semifronts for c > c .
Let Δ ( λ , c ) = d v R J ( y ) ( e λ y 1 ) d y c λ γ . For convenience, for any function ω ( ξ ) defined in R , we denote R J ( ξ y ) ω ( y ) d y as ( J ω ) ( ξ ) .
Consider the following linearized system:
c φ ( ξ ) = β w 0 ϕ ( ξ ) μ φ ( ξ ) , c ϕ ( ξ ) = d v [ ( J ϕ ) ( ξ ) ϕ ( ξ ) ] + p φ ( ξ ) γ ϕ ( ξ ) ,
where ( J ϕ ) ( ξ ) = R J ( ξ y ) ϕ ( y ) d y . Let ( φ ( ξ ) , ϕ ( ξ ) ) = ( φ 0 , ϕ 0 ) e λ 0 ξ be a solution of system (6), where λ 0 0 ,   φ 0 > 0 , ϕ 0 > 0 . Then, we have
β w 0 ϕ 0 = ( c λ 0 + μ ) φ 0 p φ 0 = Δ ( λ 0 , c ) ϕ 0 ,
which implies that Δ ( λ 0 , c ) = p φ 0 / ϕ 0 < 0 . Actually, when there exist λ 0 0 , φ 0 > 0 and ϕ 0 > 0 such that ( φ 0 , ϕ 0 ) e λ 0 ξ is a solution of system (6), we have Δ ( λ 0 , c ) < 0 ; otherwise, the expression of Δ ( λ , c ) implies that the sign of Δ ( λ , c ) is unclear.
In the following, we will study the existence of λ 0 0 , φ 0 > 0 , and ϕ 0 > 0 that satisfy (7). Firstly, we provide the range of λ such that Δ ( λ , c ) < 0 .
Lemma 1.
Δ ( λ , c ) = 0 admits a positive root λ + ( c ) such that
Δ ( λ , c ) < 0 f o r   a l l λ [ 0 , λ + ( c ) ) .
Proof. 
It is clear that for all c 0 ,
Δ ( 0 , c ) = γ < 0 , Δ ( λ , c ) λ = d v R J ( y ) ( y ) e λ y d y c ,
2 Δ ( λ , c ) λ 2 = d v R J ( y ) y 2 e λ y d y > 0 , λ R .
By assumption ( H ) , there exist r 1 and r 2 satisfying 0 < r 1 < r 2 such that
ϖ : = min x [ r 2 , r 1 ] J ( x ) > 0 .
It follows by Taylor’s formula that
d v ϖ e λ r 1 λ c λ = d v ϖ 1 λ + r 1 + λ r 1 2 2 ! + λ 3 r 1 4 4 ! + λ 4 r 1 5 5 ! + · · · + λ n 1 r 1 n n ! + · · · + λ d v ϖ λ r 1 3 3 ! c ,
where n ! = 1 × 2 × 3 × 4 · · · × n . Thus, for any c 0 , there exists λ 1 ( c ) 0 such that d v ϖ e λ r 1 λ > c λ for all λ > λ 1 ( c ) . Therefore, for c > 0 and λ > λ 1 ( c ) , one gets
Δ ( λ , c ) = d v R J ( y ) ( e λ y 1 ) d y c λ γ d v r 2 r 1 J ( y ) e λ y d y d v c λ γ d v ϖ e λ r 1 ( e λ ( r 2 r 1 ) 1 ) λ d v c λ γ . c λ ( e λ ( r 2 r 1 ) 2 ) d v γ + ,
as λ + , which indicates that Δ ( + , c ) = + , c > 0 . Similarly, it is easy to get Δ ( + , 0 ) = + .
Clearly, when c > 0 , Δ ( λ , c ) λ | ( 0 , c ) = c < 0 , which together with the above results implies that the conclusion is valid when c > 0 . In the case where c = 0 , the assumption ( H ) implies Δ ( λ , 0 ) is a limited function for λ > 0 . Moreover, we have Δ ( λ , c ) λ | ( 0 , 0 ) = 0 and Δ ( λ , c ) λ | ( λ , 0 ) = d v R J ( y ) ( y ) e λ y d y > 0 for λ > 0 , which together with Δ ( 0 , 0 ) < 0 , Δ ( + , 0 ) = + and 2 Δ ( λ , c ) λ 2 > 0 guarantees that the conclusion is also valid for c = 0 . This completes the proof.  □
For λ [ 0 , λ + ( c ) ) , we define
M ( λ , c ) = c λ + μ 0 0 Δ ( λ , c ) 1 0 β w 0 p 0 = 0 β w 0 c λ + μ p Δ ( λ , c ) 0 .
Then, by (7), we have M ( λ , c ) ( φ 0 , ϕ 0 ) T = ( φ 0 , ϕ 0 ) T . Let L ( λ , c ) be the maximum eigenvalue of M ( λ , c ) when λ [ 0 , λ + ( c ) ) . Clearly,
L 2 ( λ , c ) = p β w 0 ( c λ + μ ) Δ ( λ , c ) .
Denote
l ( λ , c ) : = L 2 ( λ , c ) , λ [ 0 , λ + ( c ) ) .
Lemma 2.
Suppose that R 0 > 1 . Then, there exist c > 0 and λ ( 0 , λ + ( c ) ) such that
l ( λ , c ) = 1 , l ( λ , c ) λ ( λ , c ) = 0 .
Furthermore, the following conclusions hold.
( i )
If c > c , the equation l ( λ , c ) = 1 admits two positive roots λ 1 and λ 2 satisfying λ 1 < λ < λ 2 < λ + ( c ) , l ( λ , c ) > 1 for λ ( 0 , λ 1 ) ( λ 2 , λ + ( c ) ) , l ( λ , c ) < 1 for λ ( λ 1 , λ 2 ) , and
l ( λ , c ) λ λ = λ 1 < 0 , l ( λ , c ) λ λ = λ 2 > 0 .
( i i )
If c = c , then λ 1 = λ 2 = λ < λ + ( c ) , and l ( λ , c ) 1 for all λ ( 0 , λ + ( c ) ) .
( i i i )
If 0 < c < c , then l ( λ , c ) > 1 for all λ ( 0 , λ + ( c ) ) .
Proof. 
If c = 0 , then it follows from the proof of Lemma 1 and
Δ ( λ , 0 ) λ λ = 0 = d v R J ( y ) ( y ) d y = 0
that Δ ( λ , 0 ) > γ for λ > 0 . Thus, l ( λ , 0 ) > 1 for all λ ( 0 , λ + ( 0 ) ) . By some calculations, it is easy to show that
l ( 0 , c ) = p β w 0 μ γ > 1 , l ( λ , c ) 0 as c + , l ( λ , c ) + as λ λ + ( c ) ,
l ( λ , c ) c = p β w 0 λ ( Δ ( λ , c ) ( c λ + μ ) ) ( c λ + μ ) 2 Δ ( λ , c ) 2 < 0 , for all c 0 , λ ( 0 , λ + ( c ) ) ,
2 l ( λ , c ) λ 2 = p β w 0 { ( c λ + μ ) 3 Δ λ λ ( λ , c ) Δ 2 ( λ , c ) ( c λ + μ ) Δ ( λ , c ) [ c Δ ( λ , c ) + ( c λ + μ ) Δ λ ( λ , c ) ] 2 ( c λ + μ ) Δ ( λ , c ) [ c 2 Δ 2 ( λ , c ) + ( c λ + μ ) 2 Δ λ 2 ( λ , c ) ] } / [ ( c λ + μ ) Δ ( λ , c ) ] 4 > 0 , for all c > 0 , λ ( 0 , λ + ( c ) ) ,
where Δ λ ( λ , c ) = Δ ( λ , c ) λ , Δ λ λ ( λ , c ) = 2 Δ ( λ , c ) λ 2 . Hence, by using the above results, the conclusions can be easily obtained.  □
Remark 1.
c obtained in Lemma 2 is the minimal wave speed of system (3), which will be proved later.
For any c > c , it follows from Lemma 2 that there exists λ c such that l ( λ c , c ) = 1 and l ( λ , c ) λ | λ = λ c < 0 . Thus, there exist q 1 > 0 and q 2 > 0 such that
M ( λ c + ϵ , c ) ( q 1 , q 2 ) T = L ( λ c + ϵ , c ) ( q 1 , q 2 ) T < ( q 1 , q 2 ) T ,
for ϵ > 0 small enough, where L ( λ c + ϵ , c ) is the maximum eigenvalue of M ( λ c + ϵ , c ) . Therefore,
( c ( λ c + ϵ ) + μ ) q 1 + β w 0 q 2 < 0 , p q 1 + Δ ( λ c + ϵ , c ) q 2 < 0 .
Recall that l ( λ c , c ) = 1 implies that 1 is the maximum eigenvalue of M ( λ c , c ) . Therefore, for c > c , there exist φ 0 > 0 and ϕ 0 > 0 such that ( φ 0 , ϕ 0 ) satisfies (7) and ( φ ( ξ ) , ϕ ( ξ ) ) = ( φ 0 , ϕ 0 ) e λ c ξ is a solution of (6).
In the following, we always assume that R 0 > 1 and c > c . Let
W ¯ ( ξ ) = w 0 , U ¯ ( ξ ) = min φ 0 e λ c ξ , β w 0 α μ , V ¯ ( ξ ) = min ϕ 0 e λ c ξ , p β w 0 α μ γ , and W ̲ ( ξ ) = max { 0 , w 0 σ e ϵ 1 ξ } , U ̲ ( ξ ) = max { 0 , e λ c ξ ( φ 0 M q 1 e ϵ ξ ) } , V ̲ ( ξ ) = max { 0 , e λ c ξ ( ϕ 0 M q 2 e ϵ ξ ) } ,
where σ , ϵ 1 , ϵ and M can be defined later. We always assume that
M > max φ 0 q 1 β w 0 α μ φ 0 ϵ λ c , ϕ 0 q 2 p β w 0 α μ γ ϕ 0 ϵ λ c ,
which can ensure U ¯ ( ξ ) > U ̲ ( ξ ) , V ¯ ( ξ ) > V ̲ ( ξ ) for all ξ R . Other restrictions on M can be found later.
For convenience, denote
x 1 = 1 λ c ln β w 0 α μ φ 0 , x 2 = 1 λ c ln p β w 0 α μ γ ϕ 0 , x 3 = 1 ϵ 1 ln w 0 σ , x 4 = 1 ϵ ln φ 0 M q 1 , x 5 = 1 ϵ ln ϕ 0 M q 2 .
Lemma 3.
The functions W ¯ ( ξ ) , U ¯ ( ξ ) and V ¯ ( ξ ) satisfy
c W ¯ ( ξ ) s b W ¯ ( ξ ) β W ¯ ( ξ ) V ̲ ( ξ ) 1 + α V ̲ ( ξ ) , f o r a n y ξ R ;
c U ¯ ( ξ ) β W ¯ ( ξ ) V ¯ ( ξ ) 1 + α V ¯ ( ξ ) μ U ¯ ( ξ ) , f o r ξ x 1 ;
c V ¯ ( ξ ) d v [ ( J V ¯ ) ( ξ ) V ¯ ( ξ ) ] + p U ¯ ( ξ ) γ V ¯ ( ξ ) , f o r ξ x 2 .
Proof. 
By the definitions of W ¯ ( ξ ) and V ̲ ( ξ ) , we have
s b W ¯ ( ξ ) β W ¯ ( ξ ) V ̲ ( ξ ) 1 + α V ̲ ( ξ ) c W ¯ ( ξ ) = β W ¯ ( ξ ) V ̲ ( ξ ) 1 + α V ̲ ( ξ ) 0 .
In the case of ξ < x 1 , we have U ¯ ( ξ ) = φ 0 e λ c ξ . It follows from V ¯ ( ξ ) ϕ 0 e λ c ξ that one has
β W ¯ ( ξ ) V ¯ ( ξ ) 1 + α V ¯ ( ξ ) μ U ¯ ( ξ ) c U ¯ ( ξ ) β W ¯ ( ξ ) V ¯ ( ξ ) μ U ¯ ( ξ ) c U ¯ ( ξ ) e λ c ξ ( β w 0 ϕ 0 ( μ + c λ c ) φ 0 ) = 0 .
In the case of ξ > x 1 , then U ¯ ( ξ ) = β w 0 / ( α μ ) . Following V ¯ ( ξ ) / ( 1 + α V ¯ ( ξ ) ) 1 / α , one gets
β W ¯ ( ξ ) V ¯ ( ξ ) 1 + α V ¯ ( ξ ) μ U ¯ ( ξ ) c U ¯ ( ξ ) β w 0 α μ U ¯ ( ξ ) = 0 .
Therefore, the above two cases yield that
c U ¯ ( ξ ) β W ¯ ( ξ ) V ¯ ( ξ ) 1 + α V ¯ ( ξ ) μ U ¯ ( ξ ) , for ξ x 1 .
Note that
( J V ¯ ) ( ξ ) x 2 J ( ξ y ) ϕ 0 e λ c y d y + x 2 + J ( ξ y ) p β w 0 α μ γ d y min p β w 0 α μ γ , ϕ 0 e λ c ξ R J ( y ) e λ c y d y .
If ξ > x 2 , then V ¯ ( ξ ) = p β w 0 / ( α μ γ ) . Since U ¯ ( ξ ) β w 0 / ( α μ ) ,
d v [ ( J V ¯ ) ( ξ ) V ¯ ( ξ ) ] + p U ¯ ( ξ ) γ V ¯ ( ξ ) c V ¯ ( ξ ) p U ¯ ( ξ ) γ V ¯ ( ξ ) 0 .
If ξ < x 2 , then V ¯ ( ξ ) = ϕ 0 e λ c ξ . Thus,
d v [ ( J V ¯ ) ( ξ ) V ¯ ( ξ ) ] + p U ¯ ( ξ ) γ V ¯ ( ξ ) c V ¯ ( ξ ) e λ c ξ ( Δ ( λ c , c ) ϕ 0 + p φ 0 ) = 0 ,
by the fact that U ¯ ( ξ ) φ 0 e λ c ξ . This completes the proof.  □
Lemma 4.
For ϵ 1 ( 0 , λ c ) and σ > max { w 0 , β w 0 ϕ 0 / b } , the function W ̲ ( ξ ) satisfies
c W ̲ ( ξ ) s b W ̲ ( ξ ) β W ̲ ( ξ ) V ¯ ( ξ ) 1 + α V ¯ ( ξ ) , f o r ξ x 3 .
Proof. 
If ξ > x 3 , then W ̲ ( ξ ) = 0 , and the conclusion clearly holds. Thus, it needs only to be shown that the conclusion is valid for ξ < x 3 . If ξ < x 3 < 0 , then W ̲ ( ξ ) = w 0 σ e ϵ 1 ξ . Therefore,
s b W ̲ ( ξ ) β W ̲ ( ξ ) V ¯ ( ξ ) 1 + α V ¯ ( ξ ) c W ̲ ( ξ ) = ( b + c ϵ 1 ) σ e ϵ 1 ξ β W ̲ ( ξ ) V ¯ ( ξ ) 1 + α V ¯ ( ξ ) ( b + c ϵ 1 ) σ e ϵ 1 ξ β w 0 ϕ 0 e λ c ξ e λ c ξ ( ( b + c ϵ 1 ) σ β w 0 ϕ 0 ) > 0 .
This completes the proof.  □
Lemma 5.
Let ϵ be small enough to satisfy 0 < ϵ < min { λ c , ϵ 1 } and M be large enough. Then, the functions U ̲ ( ξ ) and V ̲ ( ξ ) satisfy
c U ̲ ( ξ ) β W ̲ ( ξ ) V ̲ ( ξ ) 1 + α V ̲ ( ξ ) μ U ̲ ( ξ ) , f o r ξ x 4 ;
c V ̲ ( ξ ) d v [ ( J V ̲ ) ( ξ ) V ̲ ( ξ ) ] + p U ̲ ( ξ ) γ V ̲ ( ξ ) , f o r ξ x 5 .
Proof. 
We only consider the case in which x 4 x 5 ; the others can be considered similarly. Let
M > max φ 0 q 1 , ϕ 0 q 2 , ϕ 0 q 2 w 0 σ ϵ ϵ 1 , β w 0 α ϕ 0 2 + β σ ϕ 0 q 1 ( μ + c ( λ c + ϵ ) ) β w 0 q 2 .
Then, x 4 x 5 x 3 < 0 .
Obviously, the first inequality holds for ξ > x 4 , and the second one is valid for ξ > x 5 .
If x 4 ξ < x 5 , then U ̲ ( ξ ) = 0 , V ̲ ( ξ ) = e λ c ξ ( ϕ 0 M q 2 e ϵ ξ ) . Hence,
d v [ ( J V ̲ ) ( ξ ) V ̲ ( ξ ) ] + p U ̲ ( ξ ) γ V ̲ ( ξ ) c V ̲ ( ξ ) d v R J ( ξ y ) e λ c y ( ϕ 0 M q 2 e ϵ y ) d y d v V ̲ ( ξ ) γ V ̲ ( ξ ) c V ̲ ( ξ ) = e λ c ξ ϕ 0 Δ ( λ c , c ) M q 2 e ( λ c + ϵ ) ξ Δ ( λ c + ϵ , c ) p e λ c ξ ( φ 0 M q 1 e ϵ ξ ) 0 .
If ξ < x 4 , then W ̲ ( ξ ) = w 0 σ e ϵ 1 ξ , U ̲ ( ξ ) = e λ c ξ ( φ 0 M q 1 e ϵ ξ ) , V ̲ ( ξ ) = e λ c ξ ( ϕ 0 M q 2 e ϵ ξ ) . Thus, it is easy to show that
d v [ ( J V ̲ ) ( ξ ) V ̲ ( ξ ) ] + p U ̲ ( ξ ) γ V ̲ ( ξ ) c V ̲ ( ξ ) e ( λ c + ϵ ) ξ M ( q 2 Δ ( λ c + ϵ , c ) + p q 1 ) 0 .
Directly, one has
β W ̲ ( ξ ) V ̲ ( ξ ) 1 + α V ̲ ( ξ ) = β w 0 V ̲ ( ξ ) 1 + α V ̲ ( ξ ) β w 0 V ̲ ( ξ ) + β w 0 V ̲ ( ξ ) β σ e ϵ 1 ξ V ̲ ( ξ ) 1 + α V ̲ ( ξ ) = β w 0 V ̲ ( ξ ) β w 0 α V ̲ 2 ( ξ ) 1 + α V ̲ ( ξ ) β σ e ϵ 1 ξ V ̲ ( ξ ) 1 + α V ̲ ( ξ ) β w 0 V ̲ ( ξ ) β w 0 α ϕ 0 2 e 2 λ c ξ β σ ϕ 0 e ( λ c + ϵ 1 ) ξ ,
which together with ξ < 0 yields that
β W ̲ ( ξ ) V ̲ ( ξ ) 1 + α V ̲ ( ξ ) μ U ̲ ( ξ ) c U ̲ ( ξ ) M e ( λ c + ϵ ) ξ q 1 ( μ + c ( λ c + ϵ ) ) β w 0 q 2 β w 0 α ϕ 0 2 e 2 λ c ξ β σ ϕ 0 e ( λ c + ϵ 1 ) ξ e ( λ c + ϵ ) ξ ( M ( q 1 ( μ + c ( λ c + ϵ ) ) β w 0 q 2 ) β w 0 α ϕ 0 2 β σ ϕ 0 ) 0 .
The proof is complete.  □
For X > X ¯ : = max { | x 1 | , | x 2 | , x 3 , x 4 , x 5 } , define
Γ X = ( ψ ( · ) , φ ( · ) , ϕ ( · ) ) C ( [ X , X ] , R 3 ) W ̲ ( ξ ) ψ ( ξ ) W ¯ ( ξ ) , for ξ [ X , X ] ; U ̲ ( ξ ) φ ( ξ ) U ¯ ( ξ ) , for ξ [ X , X ] ; V ̲ ( ξ ) ϕ ( ξ ) V ¯ ( ξ ) , for ξ [ X , X ] ; ψ ( X ) = W ̲ ( X ) , φ ( X ) = U ̲ ( X ) , ϕ ( X ) = V ̲ ( X ) , .
For any ( ψ ( · ) , φ ( · ) , ϕ ( · ) ) Γ X , consider the following truncated problem
c W ( ξ ) = s b W ( ξ ) g ( ψ ( ξ ) , ϕ ( ξ ) ) , ξ ( X , X ] , c U ( ξ ) = g ( ψ ( ξ ) , ϕ ( ξ ) ) μ U ( ξ ) , ξ ( X , X ] , c V ( ξ ) = d v R J ( ξ y ) ϕ ^ ( y ) d y d v V ( ξ ) + p φ ( ξ ) γ V ( ξ ) , ξ ( X , X ] , W ( X ) = W ̲ ( X ) , U ( X ) = U ̲ ( X ) , V ( X ) = V ̲ ( X ) ,
where g ( ψ ( ξ ) , ϕ ( ξ ) ) = β ψ ( ξ ) ϕ ( ξ ) / ( 1 + α ϕ ( ξ ) ) , and
ϕ ^ ( ξ ) = ϕ ( X ) , ξ > X , ϕ ( ξ ) , | ξ | X , V ̲ ( ξ ) , ξ < X .
Then, it follows from the standard theory of ordinary differential equations that system (11) admits a unique solution ( W X ( ξ ) , U X ( ξ ) , V X ( ξ ) ) satisfying W X ( · ) , U X ( · ) , V X ( · ) C 1 ( [ X , X ] ) for any ( ψ ( · ) , φ ( · ) , ϕ ( · ) ) Γ X .
Define F = ( F 1 , F 2 , F 3 ) : Γ X C ( [ X , X ] , R 3 ) by
F 1 ( ψ ( ξ ) , φ ( ξ ) , ϕ ( ξ ) ) = W X ( ξ ) , F 2 ( ψ ( ξ ) , φ ( ξ ) , ϕ ( ξ ) ) = U X ( ξ ) , F 3 ( ψ ( ξ ) , φ ( ξ ) , ϕ ( ξ ) ) = V X ( ξ ) .
Lemma 6.
The operator F satisfies F ( Γ X ) Γ X . Moreover, operator F is completely continuous.
Proof. 
By using Lemmas 3–5 and similar arguments to those in ([24], Theorem 2.5), it is easy to show that F ( Γ X ) Γ X .
Now, we show that F is completely continuous. Let ( W X ( ξ ) , U X ( ξ ) , V X ( ξ ) ) be the unique solution of system (11) with ( ψ ( ξ ) , φ ( ξ ) , ϕ ( ξ ) ) Γ X . Then, we can obtain that
W X ( ξ ) = W ̲ ( X ) e b c ( ξ + X ) + 1 c X ξ ( s g ( ψ ( η ) , ϕ ( η ) ) ) e b c ( ξ η ) d η ,
U X ( ξ ) = U ̲ ( X ) e μ c ( ξ + X ) + 1 c X ξ e μ c ( ξ η ) g ( ψ ( η ) , ϕ ( η ) ) d η ,
V X ( ξ ) = V ̲ ( X ) e d v + γ c ( ξ + X ) + 1 c X ξ e d v + γ c ( ξ η ) p φ ( η ) + d v R J ( η y ) ϕ ^ ( y ) d y d η .
For any ( ψ 1 ( ξ ) , φ 1 ( ξ ) , ϕ 1 ( ξ ) ) Γ X and ( ψ 2 ( ξ ) , φ 2 ( ξ ) , ϕ 2 ( ξ ) ) Γ X , it is obvious that
R J ( η y ) ϕ ^ 1 ( y ) d y R J ( η y ) ϕ ^ 2 ( y ) d y = X X J ( η y ) ( ϕ 1 ( y ) ϕ 2 ( y ) ) d y + X + J ( η y ) ( ϕ 1 ( X ) ϕ 2 ( X ) ) d y .
Hence,
R J ( η y ) ϕ ^ 1 ( y ) d y R J ( η y ) ϕ ^ 2 ( y ) d y 2 max ξ [ X , X ] | ϕ 1 ( ξ ) ϕ 2 ( ξ ) | .
Then, by the definition of ( W X ( ξ ) , U X ( ξ ) , V X ( ξ ) ) , it is clear that F is continuous. Furthermore, since W X ( · ) , U X ( · ) , V X ( · ) C 1 ( [ X , X ] ) , we have that W X , U X and V X are uniformly bounded on [ X , X ] according to Equation (11). Therefore, we can get that operator F is compact on Γ X . This completes the proof.  □
It is obvious that Γ X is a closed and convex set. Then, it follows from Lemma 6 and Schauder’s fixed-point theorem that operator F admits a fixed point ( W X ( · ) , U X ( · ) , V X ( · ) ) Γ X ; that is,
( W X ( ξ ) , U X ( ξ ) , V X ( ξ ) ) = F ( W X ( ξ ) , U X ( ξ ) , V X ( ξ ) ) for any ξ [ X , X ] .
For simplicity, we drop the superscript ∗ and denote the fixed point as ( W X ( ξ ) , U X ( ξ ) , V X ( ξ ) ) in the following.
Define
C 1 , 1 ( [ X , X ] ) = u ( · ) C 1 ( [ X , X ] ) | u ( · ) and u ( · ) are Lipschitz continuous ,
with norm
u C 1 , 1 ( [ X , X ] ) = max ξ [ X , X ] | u ( ξ ) | + max ξ [ X , X ] | u ( ξ ) | + max ξ , η [ X , X ] , ξ η | u ( ξ ) u ( η ) | | ξ η | .
Now, we give some estimations of W X ( ξ ) , U X ( ξ ) , V X ( ξ ) in the space C 1 , 1 ( [ X , X ] ) .
Theorem 1.
There exists a positive constant C independent of X such that
W X C 1 , 1 ( [ X , X ] ) C , U X C 1 , 1 ( [ X , X ] ) C , V X C 1 , 1 ( [ X , X ] ) C ,
for any X > X ¯ .
Proof. 
Since ( W X ( · ) , U X ( · ) , V X ( · ) ) is a fixed point of F on Γ X , we have
W X ( ξ ) w 0 , U X ( ξ ) β w 0 α μ , V X ( ξ ) p β w 0 α μ γ , for any ξ [ X , X ] ,
and
c W X ( ξ ) = s b W X ( ξ ) β V X ( ξ ) 1 + α V X ( ξ ) W X ( ξ ) , c U X ( ξ ) = β V X ( ξ ) 1 + α V X ( ξ ) W X ( ξ ) μ U X ( ξ ) , c V X ( ξ ) = d v R J ( ξ y ) V ^ X ( y ) d y d v V X ( ξ ) + p U X ( ξ ) γ V X ( ξ ) .
where
V ^ X ( ξ ) = V X ( X ) , ξ > X , V X ( ξ ) , | ξ | X , V ̲ ( ξ ) , ξ < X .
Obviously, V ^ X ( ξ ) p β w 0 / ( α μ γ ) for any ξ R .
By simple calculations, we can obtain that
| W X ( ξ ) | L 1 = 1 c s + b w 0 + β w 0 α ,
| U X ( ξ ) | L 2 = 2 β w 0 c α ,
| V X ( ξ ) | L 3 = 2 p β w 0 c α μ 1 + d v γ .
Denote C 1 = max { L 1 , L 2 , L 3 } ; one has
| W X ( ξ ) W X ( η ) | C 1 | ξ η | , | U X ( ξ ) U X ( η ) | C 1 | ξ η | , | V X ( ξ ) V X ( η ) | C 1 | ξ η | ,
for any ξ , η [ X , X ] . Furthermore,
| W X ( ξ ) W X ( η ) | = 1 c b ( W X ( η ) W X ( ξ ) ) + β V X ( η ) ( W X ( η ) W X ( ξ ) ) 1 + α V X ( η ) + β W X ( ξ ) ( V X ( η ) V X ( ξ ) ) 1 + α V X ( ξ ) + β α V X ( η ) W X ( ξ ) ( V X ( ξ ) V X ( η ) ) ( 1 + α V X ( η ) ) ( 1 + α V X ( ξ ) ) 1 c b + β α | W X ( η ) W X ( ξ ) | + 2 β w 0 | V X ( η ) V X ( ξ ) | 1 c b + β α + 2 β w 0 C 1 | ξ η | .
Similarly, we can obtain
| U X ( ξ ) U X ( η ) | 1 c μ + β α + 2 β w 0 C 1 | ξ η | .
It follows from assumption ( H ) that the kernel function J is Lipschitz continuous. Let Q be its Lipschitz constant. Then, by similar arguments to the proofs in ([25], Theorem 2.8), it is easy to show
R ( J ( ξ y ) J ( η y ) ) V ^ X ( y ) d y ϕ 0 Q λ c + 3 p β ω 0 J L α μ γ + C 1 ξ η .
Thus,
| V X ( ξ ) V X ( η ) | d v c R ( J ( ξ y ) J ( η y ) ) V ^ X ( y ) d y + ( d v + γ + p ) C 1 c | ξ η | 1 c d v ϕ 0 Q λ c + 3 p β ω 0 J L α μ γ + C 1 + ( d v + γ + p ) C 1 | ξ η | .
Combining the above arguments, the conclusion follows. This completes the proof.  □
Lemma 7.
Suppose that R 0 > 1 . For any c > c , system (3) admits a positive traveling semifront ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) satisfying
0 < W ( ξ ) < w 0 , 0 < U ( ξ ) β w 0 α μ , 0 < V ( ξ ) p β w 0 α μ γ , f o r a n y ξ R .
Proof. 
Let { X n } n N + be an increasing sequence satisfying lim n + X n = + and X n > X ¯ for any n N + . Then, for any n N + , operator F has a fixed point ( W X n , U X n , V X n ) on Γ X n . Therefore, it follows from Theorem 1 and Arzela–Ascoli’s theorem that there exists a subsequence { X n k } k N + such that ( W X n k , U X n k , V X n k ) ( W , U , V ) in C l o c 1 ( R ) as k + , for some ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) [ C 1 ( R ) ] 3 . Furthermore, the Lebesgue’s dominated convergence theorem yields that
lim k + R J ( ξ y ) V ^ X n k ( y ) d y = R J ( ξ y ) V ( y ) d y = ( J V ) ( ξ ) .
Thus, ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) satisfies system (5), and for any ξ R ,
W ̲ ( ξ ) W ( ξ ) w 0 , U ̲ ( ξ ) U ( ξ ) U ¯ ( ξ ) β w 0 α μ , V ̲ ( ξ ) V ( ξ ) V ¯ ( ξ ) p β w 0 α μ γ ,
which can guarantee that lim ξ W ( ξ ) = w 0 , lim ξ U ( ξ ) = 0 , lim ξ V ( ξ ) = 0 .
Now, we show that 0 < W ( ξ ) < w 0 , U ( ξ ) > 0 , V ( ξ ) > 0 , for any ξ R . Suppose that there exists ξ 0 R such that W ( ξ 0 ) = 0 , and then W ( ξ 0 ) = 0 . Thus, by the first equation of (5), we have 0 = c W ( ξ 0 ) = s , which is a contradiction. Hence, W ( ξ ) > 0 for any ξ R . The second and third equations of (5) yield that
U ( ξ ) μ U ( ξ ) c , and V ( ξ ) d v + γ c V ( ξ ) for any ξ R .
By the comparison principle and the fact that U ( ξ ) U ̲ ( ξ ) > 0 for ξ < x 4 and V ( ξ ) V ̲ ( ξ ) > 0 for ξ < x 5 , it is easy to get that U ( ξ ) > 0 and V ( ξ ) > 0 for all ξ R . If there exists ξ 1 R such that W ( ξ 1 ) = w 0 , then W ( ξ 1 ) = 0 . Hence, the first equation of (5) implies that V ( ξ 1 ) = 0 , a contradiction. This completes the proof.  □
In order to show that the traveling semifront obtained in Lemma 7 is indeed a traveling front, we need to show that the following lemma holds.
Lemma 8.
Suppose that R 0 > 1 . For any c > c , denote ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) as the traveling semifront of system (3) obtained in Lemma 7. Then, lim inf ξ + V ( ξ ) > 0 .
Proof. 
Suppose by contradiction that there exists a nondecreasing sequence { ξ n } n N + satisfying ξ n + as n + and V ( ξ n ) 0 such that lim n + V ( ξ n ) = 0 . Define
W n ( ξ ) = W ( ξ n + ξ ) , U n ( ξ ) = U ( ξ n + ξ ) , V n ( ξ ) = V ( ξ n + ξ ) .
By similar arguments to those in Theorem 1, it is clear that ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) ) , ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) ) and W n C 1 , 1 ( R ) , U n C 1 , 1 ( R ) , V n C 1 , 1 ( R ) are all uniformly bounded. Thus, Arzela–Ascoli’s theorem implies that there exists a subsequence { n k } k N + such that
( W n k ( ξ ) , U n k ( ξ ) , V n k ( ξ ) ) ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) , in C l o c 1 ( R ) as k + ,
for some ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) [ C 1 ( R ) ] 3 . For simplicity, denote ( W n k ( ξ ) , U n k ( ξ ) , V n k ( ξ ) ) as ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) ) . Obviously, V ( 0 ) = 0 . Then, by similar arguments to those in ([25], Theorem 2.9), we can get that W = w 0 and V = 0 in R . Let
φ n ( ξ ) = U n ( ξ ) max { U n ( ξ ) , V n ( ξ ) } , ϕ n ( ξ ) = V n ( ξ ) max { U n ( ξ ) , V n ( ξ ) } .
Then, ( φ n ( ξ ) , ϕ n ( ξ ) ) satisfies φ n ( ξ ) 1 ,   ϕ n ( ξ ) 1 and
c φ n ( ξ ) = β W n ( ξ ) ϕ n ( ξ ) 1 + α V n ( ξ ) μ φ n ( ξ ) , c ϕ n ( ξ ) = d v [ ( J ϕ n ) ( ξ ) ϕ n ( ξ ) ] + p φ n ( ξ ) γ ϕ n ( ξ ) .
By a similar method to above, we can get that there exists a subsequence, denoted by ( φ n ( ξ ) , ϕ n ( ξ ) ) for simplicity, such that ( φ n ( ξ ) , ϕ n ( ξ ) ) ( φ ( ξ ) , ϕ ( ξ ) ) in C l o c 1 ( R ) as n + , where 0 φ ( ξ ) 1 ,   0 ϕ ( ξ ) 1 and
c φ ( ξ ) = β w 0 ϕ ( ξ ) μ φ ( ξ ) , c ϕ ( ξ ) = d v [ ( J ϕ ) ( ξ ) ϕ ( ξ ) ] + p φ ( ξ ) γ ϕ ( ξ ) .
is satisfied.
Clearly, max { φ , ϕ } = 1 . Note that if φ ( ξ 1 ) = 0 or ϕ ( ξ 1 ) = 0 for some ξ 1 R , then φ ( ξ ) 0 and ϕ ( ξ ) 0 for all ξ R . Thus, φ ( ξ ) > 0 and ϕ ( ξ ) > 0 for all ξ R .
Let φ ( x , t ) = φ ( x + c t ) ,   ϕ ( x , t ) = ϕ ( x + c t ) . Then, ( φ ( x , t ) ,   ϕ ( x , t ) ) satisfies
φ ( x , t ) t = β w 0 ϕ ( x , t ) μ φ ( x , t ) β w 0 ϕ ( x , t ) 1 + α 1 ϕ ( x , t ) μ φ ( x , t ) , t > 0 , x R , ϕ ( x , t ) t = d v [ ( J ϕ ) ( x , t ) ϕ ( x , t ) ] + p φ ( x , t ) γ ϕ ( x , t ) , t > 0 , x R , φ ( x , 0 ) = φ ( x ) , ϕ ( x , 0 ) = ϕ ( x ) , x R ,
where α 1 is small enough to satisfy ( p β w 0 μ γ ) / ( μ γ α 1 ) > 2 . Therefore, the comparison principle ([17], Lemma 2.3) implies that φ ( x , t ) φ ̲ ( x , t ) , ϕ ( x , t ) ϕ ̲ ( x , t ) for all t > 0 and x R , where ( φ ̲ ( x , t ) , ϕ ̲ ( x , t ) ) is the solution of the following system
φ ̲ ( x , t ) t = β w 0 ϕ ̲ ( x , t ) 1 + α 1 ϕ ̲ ( x , t ) μ φ ̲ ( x , t ) , t > 0 , x R , ϕ ̲ ( x , t ) t = d v [ ( J ϕ ̲ ) ( x , t ) ϕ ̲ ( x , t ) ] + p φ ̲ ( x , t ) γ ϕ ̲ ( x , t ) , t > 0 , x R , 0 < φ ̲ ( x , 0 ) < min { φ ( x ) , φ ̲ 0 } , 0 < ϕ ̲ ( x , 0 ) < min { ϕ ( x ) , ϕ ̲ 0 } , x R ,
with ( φ ̲ 0 , ϕ ̲ 0 ) = ( ( p β w 0 μ γ ) / ( μ p α 1 ) , ( p β w 0 μ γ ) / ( μ γ α 1 ) ) . Recall that R 0 > 1 . According to ([17], Theorem 3.2), we can get
lim t + ( φ ̲ ( 0 , t ) , ϕ ̲ ( 0 , t ) ) = ( φ ̲ 0 , ϕ ̲ 0 ) .
Thus, 1 < ϕ ̲ 0 / 2 lim inf t + ϕ ̲ ( 0 , t ) lim inf t + ϕ ( 0 , t ) = lim inf t + ϕ ( c t ) = ϕ ( ) , which contradicts ϕ ( ξ ) 1 . Therefore, lim inf ξ + V ( ξ ) > 0 holds.  □
With the aid of Lemmas 7 and 8, the existence of traveling fronts of system (3) for c > c can be obtained as follows.
Theorem 2.
Suppose that R 0 > 1 . For any c > c , system (3) admits a traveling front with wave speed c.
Proof. 
By Lemmas 7 and 8, we only need to prove that lim ξ + ( W ( ξ ) , U ( ξ ) , V ( ξ ) = E . Inspired by [22,26,27,28], we use the Lyapunov method to show that this conclusion holds. Let f ( z ) = z 1 ln z , χ 1 ( z ) = z + J ( y ) d y , χ 2 ( z ) = z J ( y ) d y . By the assumption ( H ) , without loss of generality, we assume that the compact support of J is [ r , r ] . Then, it is clear that
χ 1 ( y ) = 0 for y r , χ 2 ( y ) = 0 for y r , χ i ( 0 ) = 1 2 for i = 1 , 2 .
Define L ( W , U , V ) ( ξ ) = L 1 ( ξ ) + L 2 ( ξ ) , where
L 1 ( ξ ) = c w f W ( ξ ) w + c u f U ( ξ ) u + c μ v p f V ( ξ ) v ,
L 2 ( ξ ) = d v μ v p 0 + χ 1 ( y ) f V ( ξ y ) v d y d v μ v p 0 χ 2 ( y ) f V ( ξ y ) v d y ,
and ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) is the solution of system (5). It is clear that L ( W , U , V ) ( ξ ) is bounded from below by Lemmas 7 and 8. Following similar calculations to those in [22,27], one has
d L 2 ( ξ ) d ξ = d v μ v p f V ( ξ ) v R J ( y ) f V ( ξ y ) v d y .
Thus,
d L ( W , U , V ) ( ξ ) d ξ = d L 1 ( ξ ) d ξ + d L 2 ( ξ ) d ξ = b W ( W w ) 2 + μ u 4 w W W V u ( 1 + α v ) U w v ( 1 + α V ) U v u V 1 + α V 1 + α v μ α u ( V v ) 2 ( 1 + α V ) ( 1 + α v ) v + d v μ p 1 v V ( ( J V ) ( ξ ) V ( ξ ) ) + d L 2 ( ξ ) d ξ = b W ( W w ) 2 + μ u 4 w W W V u ( 1 + α v ) U w v ( 1 + α V ) U v u V 1 + α V 1 + α v μ α u ( V v ) 2 ( 1 + α V ) ( 1 + α v ) v d v μ v p R J ( y ) f V ( ξ y ) V ( ξ ) d y 0 .
Hence, L ( W , U , V ) ( ξ ) is non-increasing on ξ 0 . It is clear that d L ( W , U , V ) ( ξ ) d ξ = 0 if and only if W ( ξ ) = w , U ( ξ ) = u , V ( ξ ) = v . Then, by similar arguments to those in ([22], Theorem 2.1) or ([29], Theorem 2.3), the conclusion is valid.  □

2.2. The Existence of a Traveling Front with Wave Speed c = c

Theorem 3.
Assume that R 0 > 1 and c = c . Then, system (3) admits a traveling front with wave speed c .
Proof. 
The proof is divided into the following three steps.
Step 1. System (3) admits a bounded traveling wave solution.
Choose any nonincreasing sequence { c n } n N + satisfying c n c as n + and c < c i + 1 < c i c 1 = c + 1 for any i N + . Then, for any c n , system (5) admits a positive solution ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) ) satisfying ( W n ( ) , U n ( ) , V n ( ) ) = ( w 0 , 0 , 0 ) and 0 < W n ( ξ ) < w 0 , 0 < U n ( ξ ) β w 0 / ( α μ ) , 0 < V n ( ξ ) p β w 0 / ( α μ γ ) for any ξ R by Theorem 2. Obviously, ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) ) is uniformly bounded. Moreover, it follows by similar arguments to those in Theorem 1 that ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) ) , W n C 1 , 1 ( R ) , U n C 1 , 1 ( R ) and V n C 1 , 1 ( R ) are all uniformly bounded. Then, Arzela–Ascoli’s theorem yields that there exists a subsequence { c n k } k N + such that
( W n k ( ξ ) , U n k ( ξ ) , V n k ( ξ ) ) ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) , in C l o c 1 ( R ) as k + ,
for some ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) [ C 1 ( R ) ] 3 . It is easy to see that ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) is a nonnegative bounded solution of system (5), and 0 W ( ξ ) w 0 , 0 U ( ξ ) β w 0 / ( α μ ) , 0 V ( ξ ) p β w 0 / ( α μ γ ) for any ξ R .
Step 2.  ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) is positive.
In the following, we still denote ( W n k ( ξ ) , U n k ( ξ ) , V n k ( ξ ) ) by ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) ) for simplicity. Suppose that
sup ξ R ( w 0 W n ( ξ ) ) 0 , sup ξ R U n ( ξ ) 0 , sup ξ R V n ( ξ ) 0 , as n + .
Then, for any 0 < ϵ < ( p β w 0 μ γ ) / ( p β + 2 μ p α ) small enough satisfying p β ( w 0 ϵ ) / ( μ γ ) > 1 , there exists a N > 0 large enough such that
w 0 ϵ < W n ( ξ ) w 0 , U n ( ξ ) < ϵ , V n ( ξ ) < ϵ , for n > N and ξ R .
Thus, for n > N , we have
c U n ( ξ ) > β ( w 0 ϵ ) V n ( ξ ) 1 + α V n ( ξ ) μ U n ( ξ ) , c V n ( ξ ) = d v [ ( J V n ) ( ξ ) V n ( ξ ) ] + p U n ( ξ ) γ V n ( ξ ) .
Let u n ( x , t ) = U n ( x + c t ) , v n ( x , t ) = V n ( x + c t ) . Clearly, u n ( x , t ) and v n ( x , t ) satisfy
u n ( x , t ) t β ( w 0 ϵ ) v n ( x , t ) 1 + α v n ( x , t ) μ u n ( x , t ) , t > 0 , x R , v n ( x , t ) t = d v [ ( J v n ) ( x , t ) v n ( x , t ) ] + p u n ( x , t ) γ v n ( x , t ) , t > 0 , x R , u n ( x , 0 ) = U n ( x ) , v n ( x , 0 ) = V n ( x ) , x R .
It then follows from the comparison principle ([17], Lemma 2.3) that u n ( x , t ) u ̲ ( x , t ) , v n ( x , t ) v ̲ ( x , t ) for all t > 0 and x R , where ( u ̲ ( x , t ) , v ̲ ( x , t ) ) is the solution of the following system
u ̲ ( x , t ) t = β ( w 0 ϵ ) v ̲ ( x , t ) 1 + α v ̲ ( x , t ) μ u ̲ ( x , t ) , t > 0 , x R , v ̲ ( x , t ) t = d v [ ( J v ̲ ) ( x , t ) v ̲ ( x , t ) ] + p u ̲ ( x , t ) γ v ̲ ( x , t ) , t > 0 , x R , 0 < u ̲ ( x , 0 ) < min { U n ( x ) , u ̲ 0 } , 0 < v ̲ ( x , 0 ) < min { V n ( x ) , v ̲ 0 } , x R ,
with ( u ̲ 0 , v ̲ 0 ) = ( ( p β ( w 0 ϵ ) μ γ ) / ( μ p α ) , ( p β ( w 0 ϵ ) μ γ ) / ( μ γ α ) ) . According to ([17], Theorem 3.2), we have:
lim t + ( u ̲ ( x , t ) , v ̲ ( x , t ) ) = ( u ̲ 0 , v ̲ 0 ) .
Therefore, ϵ lim inf t + U n ( c t ) = lim inf t + u n ( 0 , t ) lim inf t + u ̲ ( 0 , t ) u ̲ 0 / 2 > ϵ , a contradiction. Hence, (17) does not hold. Recall that lim ξ ( W n ( ξ ) , U n ( ξ ) , V n ( ξ ) )  =  ( w 0 , 0 , 0 ) . Therefore, we can assume by some transformation that for ϵ > 0 small enough,
w 0 ϵ W n ( ξ ) w 0 , U n ( ξ ) ϵ , V n ( ξ ) ϵ , for ξ 0 ,
and at least one of the following three equalities holds:
W n ( 0 ) = w 0 ϵ , U n ( 0 ) = ϵ , V n ( 0 ) = ϵ .
From the definition of ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) , one has
w 0 ϵ W ( ξ ) w 0 , U ( ξ ) ϵ , V ( ξ ) ϵ , for ξ 0 ,
and at least one of the following three equalities holds:
W ( 0 ) = w 0 ϵ , U ( 0 ) = ϵ , V ( 0 ) = ϵ .
By similar arguments to those in Lemma 7, it is clear that W ( ξ ) > 0 for all ξ R . Then, the second and third equations of (5) show that U ( ξ ) 0 and V ( ξ ) 0 for all ξ R if there exists ξ 1 R such that either U ( ξ 1 ) = 0 or V ( ξ 1 ) = 0 holds. Therefore, in the case that either U ( 0 ) = ϵ or V ( 0 ) = ϵ holds, U ( ξ ) > 0 and V ( ξ ) > 0 for all ξ R . In the case that W ( 0 ) = w 0 ϵ , suppose that there exists ξ 2 R such that U ( ξ 2 ) = 0 or V ( ξ 2 ) = 0 . Then, U ( ξ ) 0 and V ( ξ ) 0 for all ξ R , which yields that c W ( 0 ) = s b W ( 0 ) = b ϵ > 0 . However, w 0 ϵ W ( ξ ) w 0 for ξ 0 and W ( 0 ) = w 0 ϵ imply that W ( 0 ) 0 , a contradiction. Thus, U ( ξ ) > 0 and V ( ξ ) > 0 for all ξ R . Therefore, combining the above arguments, we can conclude that ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) is a positive solution of system (5).
Step 3.  ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) satisfies boundary conditions (4).
We now show that ( W ( ) , U ( ) , V ( ) ) = ( w 0 , 0 , 0 ) . By the second equation of (5) and U n ( ) = 0 , it is easy to show that
U n ( ξ ) = β c ξ e μ c ( ξ z ) W n ( z ) V n ( z ) 1 + α V n ( z ) d z , for ξ R .
Therefore, V n ( ξ ) satisfies
c V n ( ξ ) = d v [ ( J V n ) ( ξ ) V n ( ξ ) ] + h ( V n ( ξ ) , ξ ) , for ξ R ,
where
h ( V n ( ξ ) , ξ ) = γ V n ( ξ ) + p β c ξ e μ c ( ξ z ) W n ( z ) V n ( z ) 1 + α V n ( z ) d z .
Choosing an ϵ > 0 small enough to satisfy p β ( w 0 ϵ ) / μ γ > 4 ϵ p β α w 0 / μ and since W n ( ξ ) w 0 ϵ for ξ 0 , some simple calculations yield that, for ξ 0 ,
h V n ( 0 , ξ ) = γ + p β c ξ e μ c ( ξ z ) W n ( z ) d z γ + p β c ξ e μ c ( ξ z ) ( w 0 ϵ ) d z = p β ( w 0 ϵ ) μ γ > 0 ,
| h V n V n ( V n ( ξ ) , ξ ) | = p β c ξ e μ c ( ξ z ) 2 α W n ( z ) ( 1 + α V n ( z ) ) 3 d z 2 p β α w 0 μ .
Since V n ( ξ ) ϵ for ξ 0 and
h ( V n ( ξ ) , ξ ) h ( 0 , ξ ) h V n ( 0 , ξ ) V n ( ξ ) = h V n ( z 1 , ξ ) V n ( ξ ) h V n ( 0 , ξ ) V n ( ξ ) = h V n V n ( z 2 , ξ ) z 1 V n ( ξ ) ,
where z 1 , z 2 [ 0 , V n ( ξ ) ] , it is clear that h ( V n ( ξ ) , ξ ) h V n ( 0 , ξ ) V n ( ξ ) / 2 for any ξ 0 . Hence, for ξ 0 , we have
c V n ( ξ ) d v [ ( J V n ) ( ξ ) V n ( ξ ) ] + h V n ( 0 , ξ ) V n ( ξ ) / 2 .
Then, by similar arguments to those in ([30], Theorem 3.4), for ξ 0 , we have,
c V n ( ξ ) + d v R J ( y ) y 0 1 V n ( ξ θ y ) d θ d y ξ h V n ( 0 , z ) V n ( z ) 2 d z 1 2 p β ( w 0 ϵ ) μ γ ξ V n ( z ) d z .
Obviously, the left side of the above inequality is less than c ϵ + d v r ϵ when ξ r . Thus, for n large enough,
r V n ( ξ ) d ξ 2 μ γ ( c ϵ + d v r ϵ ) p β ( w 0 ϵ ) μ γ ,
which implies that
r V ( ξ ) d ξ 2 μ γ ( c ϵ + d v r ϵ ) p β ( w 0 ϵ ) μ γ .
Therefore, it follows from the boundedness of V ( ξ ) that lim ξ V ( ξ ) = 0 . Meanwhile, following ([31], Lemma 2.3), we can get lim ξ V ( ξ ) = 0 . Let { η n } be a nonincreasing sequence satisfying η n as n + . Then, by the Fatou Lemma, we have
V ( ) lim inf n + ( J V ) ( η n ) lim sup n + ( J V ) ( η n ) V ( ) ,
which yields that lim n + ( ( J V ) ( η n ) V ( η n ) ) = 0 . Thus, the third equation of (5) can guarantee that lim n + U ( η n ) = 0 . Then, it follows from the arbitrariness of sequence { η n } that lim ξ U ( ξ ) = 0 .
In the following, we show that lim ξ W ( ξ ) = w 0 . Let W ̲ = lim inf ξ W ( ξ ) and W ¯ = lim sup ξ W ( ξ ) . If W ̲ < W ¯ w 0 , then there exists a sequence { x n } satisfying x n as n + such that lim n + W ( x n ) = W ̲ and W ( x n ) = 0 . Since V ( ) = 0 and W ̲ < w 0 , for 0 < ϵ < b ( w 0 W ̲ ) / ( b + β w 0 ) , there exists N > 0 large enough such that V ( x n ) < ϵ and W ( x n ) < W ̲ + ϵ for any n > N . Therefore, for n > N , the first equation of (5) yields that
0 = c W ( x n ) = s b W ( x n ) β W ( x n ) V ( x n ) 1 + α V ( x n ) s b ( W ̲ + ϵ ) β w 0 ϵ > 0 ,
a contradiction. Hence, W ̲ = W ¯ . Thus, lim ξ W ( ξ ) exists, which implies that W ( ) = 0 . Then, by the first equation of (5), V ( ) = 0 and similar arguments to above, one has that W ( ) = w 0 . In addition, by similar arguments to those in Lemma 8 and Theorem 2, we can get lim ξ + ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) = E . The proof is complete.  □

2.3. The Nonexistence of Traveling Fronts for 0 < c < c

In order to show the nonexistence of traveling fronts for 0 < c < c , we first need to prove the following lemma.
Lemma 9.
Suppose that R 0 > 1 . Let ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) be a bounded positive solution satisfying boundary conditions (4) and system (5). Then, there exists ρ > 0 such that sup ξ R { U ( ξ ) e ρ ξ } < + , sup ξ R { V ( ξ ) e ρ ξ } < + , sup ξ R { U ( ξ ) e ρ ξ } < + , and sup ξ R { V ( ξ ) e ρ ξ } < + .
Proof. 
Obviously, W ( ξ ) w 0 for all ξ R . In fact, if there exits ξ 0 R such that W ( ξ 0 ) > w 0 , then the first equation of (5) implies that W ( ξ 0 ) < 0 , which induces that W ( ) W ( ξ 0 ) > w 0 , a contradiction.
Since R 0 > 1 , there exists ϵ > 0 small enough such that
R ˜ 0 = p β ( w 0 ϵ ) μ γ ( 1 + α ϵ ) > 1 .
Note that lim ξ ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) = ( w 0 , 0 , 0 ) . Then, for a fixed ϵ > 0 small enough, there exists ξ 1 < 0 such that
w 0 ϵ < W ( ξ ) w 0 , U ( ξ ) < ϵ , V ( ξ ) < ϵ , for any ξ < ξ 1 .
Thus, the second and third equations of system (5) yield that
c U ( ξ ) > β ( w 0 ϵ ) V ( ξ ) 1 + α ϵ μ U ( ξ ) , ξ < ξ 1 , c V ( ξ ) = d v [ ( J V ) ( ξ ) V ( ξ ) ] + p U ( ξ ) γ V ( ξ ) , ξ < ξ 1 .
Denote
Q 1 ( ξ ) = ξ U ( s ) d s , Q 2 ( ξ ) = ξ V ( s ) d s , K 1 ( ξ ) = ξ Q 1 ( s ) d s , K 2 ( ξ ) = ξ Q 2 ( s ) d s .
Clearly, ξ ( J V ) ( s ) d s = R J ( y ) ξ y V ( z ) d z d y = R J ( y ) Q 2 ( ξ y ) d y = ( J Q 2 ) ( ξ ) . Therefore, integrating the both sides of (22) from to ξ with ξ < ξ 1 yields that
c U ( ξ ) > β ( w 0 ϵ ) Q 2 ( ξ ) 1 + α ϵ μ Q 1 ( ξ ) , c V ( ξ ) = d v [ ( J Q 2 ) ( ξ ) Q 2 ( ξ ) ] + p Q 1 ( ξ ) γ Q 2 ( ξ ) .
Again, integrating the both sides of (23) from to ξ with ξ < ξ 1 , we have
c Q 1 ( ξ ) > β ( w 0 ϵ ) 1 + α ϵ K 2 ( ξ ) μ K 1 ( ξ ) , c Q 2 ( ξ ) d v ξ [ ( J Q 2 ) ( s ) Q 2 ( s ) ] d s = p K 1 ( ξ ) γ K 2 ( ξ ) .
Moreover, according to that fact that y Q 2 ( ξ θ y ) is nonincreasing with respect to θ [ 0 , 1 ] , we have
d v ξ [ ( J Q 2 ) ( s ) Q 2 ( s ) ] d s = d v R J ( y ) ξ y ξ Q 2 ( s ) d s d y = d v R J ( y ) y 0 1 Q 2 ( ξ θ y ) d θ d y d v R J ( y ) y 0 1 Q 2 ( ξ ) d θ d y = 0 .
Hence, for ξ < ξ 1 , (24) yields that
c Q 1 ( ξ ) > β ( w 0 ϵ ) 1 + α ϵ K 2 ( ξ ) μ K 1 ( ξ ) , c Q 2 ( ξ ) p K 1 ( ξ ) γ K 2 ( ξ ) .
Obviously, for any τ > 0 ,
K 1 ( ξ ) = ξ Q 1 ( s ) d s = 0 + Q 1 ( ξ z ) d z 0 τ Q 1 ( ξ z ) d z τ Q 1 ( ξ τ ) ,
K 2 ( ξ ) = ξ Q 2 ( s ) d s = 0 + Q 2 ( ξ z ) d z 0 τ Q 2 ( ξ z ) d z τ Q 2 ( ξ τ ) .
Then, by the inequalities of (25), it is easy to see that
c ( p Q 1 ( ξ ) + μ Q 2 ( ξ ) ) p β ( w 0 ϵ ) 1 + α ϵ μ γ K 2 ( ξ ) p β ( w 0 ϵ ) 1 + α ϵ μ γ τ Q 2 ( ξ τ ) ,
c γ Q 1 ( ξ ) + β ( w 0 ϵ ) 1 + α ϵ Q 2 ( ξ ) p β ( w 0 ϵ ) 1 + α ϵ μ γ K 1 ( ξ ) p β ( w 0 ϵ ) 1 + α ϵ μ γ τ Q 1 ( ξ τ ) .
Therefore,
c γ + β ( w 0 ϵ ) 1 + α ϵ + p + μ ( Q 1 ( ξ ) + Q 2 ( ξ ) ) c ( p + γ ) Q 1 ( ξ ) + μ + β ( w 0 ϵ ) 1 + α ϵ Q 2 ( ξ ) p β ( w 0 ϵ ) 1 + α ϵ μ γ τ ( Q 1 ( ξ τ ) + Q 2 ( ξ τ ) ) .
Hence, there exists τ 1 > 0 such that Q 1 ( ξ τ 1 ) + Q 2 ( ξ τ 1 ) < ( Q 1 ( ξ ) + Q 2 ( ξ ) ) / 2 .
Let H ( ξ ) = ( Q 1 ( ξ ) + Q 2 ( ξ ) ) e ρ ξ , where ρ = ln 2 / τ 1 > 0 . It is clear that H ( ξ τ 1 ) < H ( ξ ) for all ξ < ξ 1 . Denote by [ r , r ] the compact support of J. Then, for ξ < ξ 1 r , we have
| ( J Q 2 ) ( ξ ) Q 2 ( ξ ) | = r r J ( y ) ξ y ξ V ( s ) d s d y ϵ r .
By the equations of (23), for ξ < ξ 1 r , we have
p c U ( ξ ) + c μ V ( ξ ) d v μ [ ( J Q 2 ) ( ξ ) Q 2 ( ξ ) ] p β ( w 0 ϵ ) 1 + α ϵ μ γ Q 2 ( ξ ) ,
γ c U ( ξ ) + c β ( w 0 ϵ ) 1 + α ϵ V ( ξ ) d v β ( w 0 ϵ ) 1 + α ϵ [ ( J Q 2 ) ( ξ ) Q 2 ( ξ ) ] p β ( w 0 ϵ ) 1 + α ϵ μ γ Q 1 ( ξ ) .
Thus, it follows from the boundedness of U ( ξ ) and V ( ξ ) that Q 1 ( ξ ) and Q 2 ( ξ ) are bounded for ξ < ξ 1 r , which together with H ( ξ τ 1 ) < H ( ξ ) indicates that H ( ξ ) is bounded as ξ . Thus, Q 1 ( ξ ) e ρ ξ and Q 2 ( ξ ) e ρ ξ are bounded as ξ .
On the one hand,
( J Q 2 ) ( ξ ) e ρ ξ = r r J ( y ) Q 2 ( ξ y ) e ρ ξ d y = ξ r ξ + r J ( ξ z ) e ρ ( ξ z ) Q 2 ( z ) e ρ z d z r r J ( y ) e ρ y d y sup z [ ξ r , ξ + r ] { Q 2 ( z ) e ρ z } ,
which yields that ( J Q 2 ) ( ξ ) e ρ ξ is bounded as ξ . On the other hand, the second and third equations of system (5) yield for ξ < ξ 1 ,
c U ( ξ ) < β ( w 0 + ϵ ) V ( ξ ) μ U ( ξ ) , c V ( ξ ) = d v [ ( J V ) ( ξ ) V ( ξ ) ] + p U ( ξ ) γ V ( ξ ) , c U ( ξ ) < β ( w 0 + ϵ ) Q 2 ( ξ ) μ Q 1 ( ξ ) , c V ( ξ ) = d v [ ( J Q 2 ) ( ξ ) Q 2 ( ξ ) ] + p Q 1 ( ξ ) γ Q 2 ( ξ ) .
Thus, one can get that U ( ξ ) e ρ ξ , V ( ξ ) e ρ ξ , U ( ξ ) e ρ ξ and V ( ξ ) e ρ ξ are all bounded as ξ . Then, since ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) satisfies (4), it is obvious that sup ξ R { U ( ξ ) e ρ ξ } < + , sup ξ R { V ( ξ ) e ρ ξ } < + , sup ξ R { U ( ξ ) e ρ ξ } < + , and sup ξ R { V ( ξ ) e ρ ξ } < + . This completes the proof.  □
Theorem 4.
Suppose that R 0 > 1 and 0 < c < c , then system (3) does not have a traveling front with wave speed c.
Proof. 
Suppose, by contradiction, that there exists ( W ( ξ ) , U ( ξ ) , V ( ξ ) ) satisfying system (5) and boundary conditions (4) with W ( ξ ) > 0 , U ( ξ ) > 0 , V ( ξ ) > 0 for all ξ R . Then, it follows from Lemma 9 that W ( ξ ) w 0 for all ξ R and there exists ρ > 0 such that the results in Lemma 9 hold.
Now, we show that sup ξ R { ( w 0 W ( ξ ) ) e ρ ξ } < + . Let T ( ξ ) = w 0 W ( ξ ) . Then, T ( ) = 0 . By the first equation of (5), we can obtain that
c T ( ξ ) + b T ( ξ ) = β W ( ξ ) V ( ξ ) 1 + α V ( ξ ) ,
which yields that
T ( ξ ) = 1 c ξ e b c ( y ξ ) β W ( y ) V ( y ) 1 + α V ( y ) d y β w 0 c ξ e b c ( y ξ ) V ( y ) d y .
Thus, sup ξ R { T ( ξ ) e ρ ξ } < + holds by using sup ξ R { V ( ξ ) e ρ ξ } < + .
For a bounded function φ ( ξ ) 0 , define its two-sided Laplace transform as
L φ ( λ ) = R φ ( ξ ) e λ ξ d ξ , λ 0 , φ ( ξ ) 0 .
Clearly, L φ ( λ ) is defined in [ 0 , λ φ + ) such that λ φ + + satisfies lim λ λ φ + L φ ( λ ) = + or λ φ + = + . Denote the two-sided Laplace transform of U ( ξ ) and V ( ξ ) by L U ( λ ) and L V ( λ ) , respectively. Obviously, λ U + ρ and λ V + ρ .
Taking the two-sided Laplace transform on both sides of the second and third equations of system (5), we get that
( c λ + μ ) L U ( λ ) = R β W ( ξ ) V ( ξ ) 1 + α V ( ξ ) e λ ξ d ξ , Δ ( λ , c ) L V ( λ ) = p L U ( λ ) .
Since
R β W ( ξ ) V ( ξ ) 1 + α V ( ξ ) e λ ξ d ξ β w 0 L V ( λ ) ,
the first equation of (27) implies that λ U + λ V + . Then, for λ ( 0 , min { λ + ( c ) , λ V + } ) , where λ + ( c ) is defined in Lemma 1, the first equation of (27) yields that
β w 0 L V ( λ ) ( c λ + μ ) L U ( λ ) = β R w 0 W ( ξ ) 1 + α V ( ξ ) V ( ξ ) e λ ξ d ξ ,
and the second equation of (27) implies that
β w 0 L V ( λ ) ( c λ + μ ) L U ( λ ) = p β w 0 + ( c λ + μ ) Δ ( λ , c ) p L V ( λ ) .
Hence,
p β w 0 + ( c λ + μ ) Δ ( λ , c ) p L V ( λ ) = β R w 0 W ( ξ ) 1 + α V ( ξ ) V ( ξ ) e λ ξ d ξ .
Since 0 < c < c , Lemmas 1 and 2 infer that p β w 0 + ( c λ + μ ) Δ ( λ , c ) > 0 for λ ( 0 , λ + ( c ) ] . Furthermore,
R w 0 W ( ξ ) 1 + α V ( ξ ) V ( ξ ) e λ ξ d ξ = R w 0 W ( ξ ) + α W ( ξ ) V ( ξ ) 1 + α V ( ξ ) V ( ξ ) e λ ξ d ξ R w 0 W ( ξ ) V ( ξ ) e λ ξ d ξ + α w 0 R V 2 ( ξ ) e λ ξ d ξ ,
which together with sup ξ R { ( w 0 W ( ξ ) ) e ρ ξ } < + and sup ξ R { V ( ξ ) e ρ ξ } < + infers that λ V + > λ + ( c ) . Thus, L V ( λ ) is well defined for all λ ( 0 , λ + ( c ) ] . It follows from Lemma 1 that Δ ( λ , c ) 0 as λ λ + ( c ) . Therefore, it follows by (28) that
0 = p β w 0 + ( c λ + μ ) Δ ( λ , c ) p L V ( λ ) β R w 0 W ( ξ ) 1 + α V ( ξ ) V ( ξ ) e λ ξ d ξ β R W ( ξ ) 1 + α V ( ξ ) V ( ξ ) e λ + ( c ) ξ d ξ > 0 , as λ λ + ( c ) ,
a contradiction. This completes the proof.  □
Remark 2.
Theorems 2, 3 and 4 imply that c defined in Lemma 2 is the minimal wave speed of system (3).

3. Discussion of Results

It was found that when the kernel function takes a special form, the model with a nonlocal dispersal operator exhibits similar wave propagation properties to the model with a fractional Laplacian operator [32]. In fact, fractional Laplacian and fractional derivatives are special cases of nonlocal dispersal operators [33,34]. As far as we know, there are few results on the propagation dynamics of the degenerate viral dynamical model with fractional diffusion or a nonlocal dispersal operator. Thus, the results obtained in this paper can not only provide some insights into the spreading speed and the propagation dynamics of a virus but also provide a basis for the propagation properties of a viral dynamical model with fractional diffusion.
Recall that system (3) is neither a cooperative system nor a competitive system. At present, there are still some difficulties in giving an exact expression for the asymptotic spreading speed of system (3) and in elucidating the relationship between the minimal wave speed and the asymptotic spreading speed. In the following, we show some numerical arguments by using MATLAB R2016a. We divide the simulation into two steps.
  • Choose an appropriate spatial domain and then discretize it. We take the domain to be [ 500 , 500 ] . The discretization step size is 0.2 , which results in 5001 ordinary differential equations. Under our specified parameters and initial values, the viruses are always away from the boundaries of the domain during our simulation.
  • Let 0.05 be the time step. We use the ode45 function in Matlab to solve the ordinary differential equations for numerical simulation.
In addition, inspired by [35], we use the slope of the boundaries of the virus’s spreading domain to estimate the asymptotic spreading speed of the virus.
We now give the estimation of the asymptotic spreading speed and show the relationship between the minimal wave speed and the asymptotic spreading speed of system (3) by simulations. Let the parameters values be
s = 2.6 × 10 4 cells mL 1 day 1 , μ = 0.26 day 1 , p = 2.9 virions day 1 cells 1 ,
β = 2.25 × 10 7 mL day 1 virions 1 , γ = 6.0 day 1 , b = 0.0026 day 1 ,
which were used for HCV infectious transmission [36]. Then, the basic reproduction number R 0 = 4.1827 > 1 . Additionally, we assume that α = 1 × 10 7 mL virions 1 , d = 0.1 mm 2 day 1 , J ( x ) with compact support [ r , r ] satisfies
J ( x ) = e 1 x 2 r 2 r r e 1 x 2 r 2 d x | x | < r , 0 | x | r ,
and the initial data w ( x , 0 ) = 1 × 10 7 for x R , u ( 0 , 0 ) = 200 ,   v ( 0 , 0 ) = 1500 , u ( x , 0 ) = v ( x , 0 ) = 0 for all x R { 0 } . Setting the radius of compact support as r = 2 , we can get that the minimal wave speed c = 0.3177 by Lemma 2 and find that system (3) admits a non-monotonic traveling front which has a hump in the profile (see Figure 1a). Let v = 0.0001 be the threshold value above which the virus can be detected. It is found that the asymptotic spreading speed is approximately equal to 0.32 > c (see Figure 1b), which implies that the asymptotic spreading speed may be larger than the minimal wave speed.
Next, we studied the influences of the diffusion ability dv and the radius r of compact support on the minimal wave speed c . Figure 2 shows that c increases as dv or r increases (the parameter values are fixed to those in Figure 1 except for dv or r). Hence, decreasing the diffusion ability or diffusion radius may postpone the spread of the virus.
Finally, we investigated the influences of the diffusion mode on the spreading speed. Assume that the virions can move either in the form of nonlocal dispersal or in the form of local dispersal (Laplace diffusion). Let the parameter values and initial data be the same as those in Figure 1 except for the radius r of compact support. Figure 3 shows that the solutions have a large hump for both local and nonlocal dispersals. It also shows that the virus with nonlocal dispersal spreads faster than the virus with local dispersal when the radius r is larger, while the inverse is true when the radius r is smaller. Thus, there may exist a threshold value r such that a virus with nonlocal dispersal and a virus with local dispersal have the same asymptotic spreading speed when r = r and a virus with nonlocal dispersal spreads faster (slower) than a virus with local dispersal when r > r ( r < r ). Hence, nonlocal dispersal can postpone the spread of a virus when the diffusion radius is smaller and accelerate the spread of a virus when the diffusion radius is larger. In fact, it is found that the minimal wave speed for nonlocal dispersal is smaller than the minimal wave speed for local dispersal when the diffusion radius is small enough, and it can surpass the minimal wave speed for local dispersal when the radius increases (see Figure 2b), where the minimal wave speed for local dispersal can be defined by similar arguments to those in Section 2.

4. Conclusions

Inspired by the phenomenon of viruses spreading like traveling waves [13], and considering the actual situation of virus transmission, we established a degenerate viral infection dynamical model with a nonlocal dispersal operator and analyzed the existence of traveling wave solutions of the model. We proved the existence of traveling wave solutions connecting the virus-free steady state and the positive steady state with wave speed c c , as well as the nonexistence of traveling wave solutions with 0 < c < c . Thus, we can conclude that c defined in Lemma 2 is the minimal wave speed of system (3). It is worth mentioning that the lower-bound estimation of the traveling wave solutions was achieved by adopting rescaling methods and the comparison principle, which is a challenge for some nonlocal models. While other methods may exist, our method is much simpler and can be easily adapted for application to other models with nonlocal dispersal.
Furthermore, the relationship between the minimal wave speed and the asymptotic spreading speed and the influences of the diffusion mode and diffusion ability on the minimal wave speed or the asymptotic spreading speed were investigated via simulations. Both the theoretical and numerical simulation results indicate the existence of traveling wave solutions of system (3), which is consistent with the evidence presented in [13]. Based on the simulations, we conclude that the asymptotic spreading speed may be larger than the minimal wave speed, and decreasing the diffusion ability or diffusion radius may postpone the spread of the virus. Nonlocal dispersal can postpone the spread of the virus when the diffusion radius is smaller and accelerate the spread when the diffusion radius is larger.
For the proposed model in this study, there is a typical characteristic, i.e., the target cell cannot move freely within the host, which is suitable for HBV or HCV infections. However, due to the diversity of viruses, there also exists some viruses, such as HIV or HTLV, for which their susceptible target cells and infected cells can move freely within the host and may have different mobilities. Therefore, if we consider nonlocal dispersal and different mobilities in both the target cells and virions, two interesting questions naturally arise that are worth further study: can the virions propagate as a traveling wave front, and what is its minimal wave speed? Moreover, during our analysis, we assumed that the kernel function K ( x ) is symmetric. However, the actual environment is very complex, and the virions may diffuse asymmetrically within the host. The traveling wave solution and minimal wave speed of a model with an asymmetric dispersal kernel function should be further studied.

Author Contributions

Conceptualization, X.R. and T.Z.; methodology, X.R., T.Z. and X.L.; software, X.R. validation, X.R., T.Z., L.L. and X.L.; formal analysis, X.R.; investigation, X.R.; writing—original draft preparation, X.R.; writing—review and editing, L.L., T.Z. and X.L.; supervision, T.Z. and X.L.; project administration and funding acquisition, X.R., T.Z. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (12071382, 11901477, 12371503, 12126349, 11601293), the Project funded by China Postdoctoral Science Foundation (2019M653816XB), and the Natural Science Foundation of Shanxi Province (202303021211003).

Data Availability Statement

Data are contained within the article.

Acknowledgments

We are grateful to the editor and reviewers for their valuable comments and suggestions that greatly improved the presentation of this paper.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. World Health Organization. World health Statistics 2017: Monitoring health for the SDGs, Sustainable Development Goals; World Health Organization: Geneva, Switzerland, 2017. [Google Scholar]
  2. Jonit United Nations Programme on HIV/AIDS (UNAIDS). Miles to Go: Closing Gaps, Breaking Barriers, Righting Injustices; UNAIDS: Geneva, Switzerland, 2018. [Google Scholar]
  3. Perelson, A.; Neumann, A.; Markowitz, M.; Leonard, J.; Ho, D. HIV-1 dynamics in vivo: Virion clearance rate, infected cell life-span, and viral generation time. Science 1996, 271, 1582–1586. [Google Scholar] [CrossRef]
  4. Perelson, A.S.; Nelson, P.W. Mathematical analysis of HIV-1 dynamics in vivo. SIAM Rev. 1999, 41, 3–44. [Google Scholar] [CrossRef]
  5. Nowak, M.A.; May, R.M. Virus Dynamics: Mathematical Principles of Immunology and Virology; Oxford University Press: Oxford, UK, 2000. [Google Scholar]
  6. Korobeinikov, A. Global properties of basic virus dynamics models. Bull. Math. Biol. 2004, 66, 879–883. [Google Scholar] [CrossRef]
  7. Manna, K.; Chakrabarty, S.P. Chronic hepatitis B infection and HBV DNA-containing capsids: Modeling and analysis. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 383–395. [Google Scholar] [CrossRef]
  8. Shu, H.; Wang, L.; Watmough, J. Sustained and transient oscillations and chaos induced by delayed antiviral immune response in an immunosuppressive infection model. J. Math. Biol. 2014, 68, 477–503. [Google Scholar] [CrossRef]
  9. Song, X.; Neumann, A.U. Global stability and periodic solution of the viral dynamics. J. Math. Anal. Appl. 2007, 329, 281–297. [Google Scholar] [CrossRef]
  10. Xu, R.; Ma, Z. An HBV model with diffusion and time delay. J. Theor. Biol. 2009, 257, 499–509. [Google Scholar] [CrossRef] [PubMed]
  11. Murooka, T.T.; Deruaz, M.; Marangoni, F.; Vrbanac, V.D.; Seung, E.; von Andrian, U.H.; Tager, A.M.; Luster, A.D.; Mempel, T.R. HIV-infected T cells are migratory vehicles for viral dissemination. Nature 2012, 490, 283–287. [Google Scholar] [CrossRef]
  12. Fackler, O.T.; Murooka, T.T.; Imle, A.; Mempel, T.R. Adding new dimensions: Towards an integrative understanding of HIV-1 spread. Nat. Rev. Microbiol. 2014, 12, 563–574. [Google Scholar] [CrossRef] [PubMed]
  13. Strain, M.C.; Richman, D.D.; Wong, J.K.; Levine, H. Spatiotemporal dynamics of HIV propagation. J. Theor. Biol. 2002, 218, 85–96. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, K.; Wang, W. Propagation of HBV with spatial dependence. Math. Biosci. 2007, 210, 78–95. [Google Scholar] [CrossRef] [PubMed]
  15. Lai, X.; Zou, X. Repulsion effect on superinfecting virions by infected cells. Bull. Math. Biol. 2014, 76, 2806–2833. [Google Scholar] [CrossRef] [PubMed]
  16. Murray, J.D. Mathematical biology: Spatial models and biomedical applications II. In Interdisciplinary Applied Mathematics, 3rd ed.; Springer: New York, NY, USA, 2003; Volume 18. [Google Scholar]
  17. Wang, J.B.; Li, W.T.; Sun, J.W. Global dynamics and spreading speeds for a partially degenerate system with non-local dispersal in periodic habitats. Proc. R. Soc. Edinb. Sect. 2018, 148, 849–880. [Google Scholar] [CrossRef]
  18. Zhao, G.; Ruan, S. Spatial and temporal dynamics of a nonlocal viral infection model. SIAM J. Appl. Math. 2018, 78, 1954–1980. [Google Scholar] [CrossRef]
  19. Bocharov, G.; Meyerhans, A.; Bessonov, N.; Trofimchuk, S.; Volpert, V. Spatiotemporal dynamics of virus infection spreading in tissues. PLoS ONE 2016, 11, e0168576. [Google Scholar] [CrossRef]
  20. Ren, X.; Tian, Y.; Liu, L.; Liu, X. A reaction-diffusion within-host HIV model with cell-to-cell transmission. J. Math. Biol. 2018, 76, 1831–1872. [Google Scholar] [CrossRef]
  21. Wang, W.; Ma, W.; Lai, X. Repulsion effect on superinfecting virions by infected cells for virus infection dynamic model with absorption effect and chemotaxis. Nonlinear Anal. Real World Appl. 2017, 33, 253–283. [Google Scholar] [CrossRef]
  22. Wang, W.; Ma, W. Travelling wave solutions for a nonlocal dispersal HIV infection dynamical model. J. Math. Anal. Appl. 2018, 457, 868–889. [Google Scholar] [CrossRef]
  23. Yang, Y.; Hsu, C.H.; Zou, L.; Zhou, J. A note on the propagation dynamics in a nonlocal dispersal HIV infection model. Proc. Am. Math. Soc. 2022, 150, 4867–4877. [Google Scholar] [CrossRef]
  24. Yang, F.Y.; Li, Y.; Li, W.T.; Wang, Z.C. Traveling waves in a nonlocal dispersal Kermack-McKendrick epidemic model. Discret. Contin. Dyn. Syst. B 2013, 18, 1969–1993. [Google Scholar] [CrossRef]
  25. Zhu, C.C.; Li, W.T.; Yang, F.Y. Traveling waves in a nonlocal dispersal SIRH model with relapse. Comput. Math. Appl. 2017, 73, 1707–1723. [Google Scholar] [CrossRef]
  26. Ducrot, A.; Magal, P. Travelling wave solutions for an infection-age structured epidemic model with external supplies. Nonlinearity 2011, 24, 2891–2911. [Google Scholar] [CrossRef]
  27. Li, Y.; Li, W.T.; Yang, F.Y. Traveling waves for a nonlocal dispersal SIR model with delay and external supplies. Appl. Math. Comput. 2014, 247, 723–740. [Google Scholar] [CrossRef]
  28. Wu, C.; Yang, Y.; Zhao, Q.; Tian, Y.; Xu, Z. Epidemic waves of a spatial SIR model in combination with random dispersal and non-local dispersal. Appl. Math. Comput. 2017, 313, 122–143. [Google Scholar] [CrossRef]
  29. Zhang, R.; Liu, S. Traveling waves for SVIR epidemic model with nonlocal dispersal. Math. Biosci. Eng. 2019, 16, 1654–1682. [Google Scholar] [CrossRef] [PubMed]
  30. Zhang, L.; Li, W.T.; Wu, S.L. Multi-type entire solutions in a nonlocal dispersal epidemic model. J. Dynam. Differ. Equations 2016, 28, 189–224. [Google Scholar] [CrossRef]
  31. Wu, J.; Zou, X. Traveling wave fronts of reaction-diffusion systems with delay. J. Dynam. Differ. Equations 2001, 13, 651–687. [Google Scholar] [CrossRef]
  32. Hou, R.; Xu, W.B. The traveling wave solutions in a mixed-diffusion epidemic model. Fractal Fract. 2022, 6, 217. [Google Scholar] [CrossRef]
  33. D’Elia, M.; Gunzburger, M. The fractional Laplacian operator on bounded domains as a special case of the nonlocal diffusion operator. Comput. Math. Appl. 2013, 66, 1245–1260. [Google Scholar] [CrossRef]
  34. Defterli, O.; D’Elia, M.; Du, Q.; Gunzburger, M.; Lehoucq, R.; Meerschaert, M.M. Fractional diffusion on bounded domains. Fract. Calc. Appl. Anal. 2015, 18, 342–360. [Google Scholar] [CrossRef]
  35. Neubert, M.G.; Parker, I.M. Projecting rates of spread for invasive species. Risk Anal. 2004, 24, 817–831. [Google Scholar] [CrossRef] [PubMed]
  36. Dahari, H.; Lo, A.; Ribeiro, R.M.; Perelson, A.S. Modeling hepatitis C virus dynamics: Liver regeneration and critical drug efficacy. J. Theor. Biol. 2007, 247, 371–381. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Solutions of system (3). (a) Evolution of virus population. (b) Evolution of the virus spreading domain.
Figure 1. Solutions of system (3). (a) Evolution of virus population. (b) Evolution of the virus spreading domain.
Fractalfract 08 00135 g001
Figure 2. The influence of parameters on minimal wave speed. (a) The influence of diffusion ability on c (nonlocal dispersal). (b) The influence of radius r of compact support.
Figure 2. The influence of parameters on minimal wave speed. (a) The influence of diffusion ability on c (nonlocal dispersal). (b) The influence of radius r of compact support.
Fractalfract 08 00135 g002
Figure 3. The concentration of v ( x , 1300 ) for local or nonlocal operator.
Figure 3. The concentration of v ( x , 1300 ) for local or nonlocal operator.
Fractalfract 08 00135 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ren, X.; Liu, L.; Zhang, T.; Liu, X. Minimal Wave Speed for a Nonlocal Viral Infection Dynamical Model. Fractal Fract. 2024, 8, 135. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract8030135

AMA Style

Ren X, Liu L, Zhang T, Liu X. Minimal Wave Speed for a Nonlocal Viral Infection Dynamical Model. Fractal and Fractional. 2024; 8(3):135. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract8030135

Chicago/Turabian Style

Ren, Xinzhi, Lili Liu, Tianran Zhang, and Xianning Liu. 2024. "Minimal Wave Speed for a Nonlocal Viral Infection Dynamical Model" Fractal and Fractional 8, no. 3: 135. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract8030135

Article Metrics

Back to TopTop