The vertex arboricity a(G) of a graph G is the minimum number of colors required to color the vertices of G such that no cycle is monochromatic. The list vertex arboricity al(G) is the list version of this concept. In this paper, we prove that if G is a planar graph without 5-cycles intersecting with 6-cycles, then al(G)≤2.
Citation: Yanping Yang, Yang Wang, Juan Liu. List vertex arboricity of planar graphs without 5-cycles intersecting with 6-cycles[J]. AIMS Mathematics, 2021, 6(9): 9757-9769. doi: 10.3934/math.2021567
[1] | Da Shi, Zhao Li, Dan Chen . New traveling wave solutions, phase portrait and chaotic patterns for the dispersive concatenation model with spatio-temporal dispersion having multiplicative white noise. AIMS Mathematics, 2024, 9(9): 25732-25751. doi: 10.3934/math.20241257 |
[2] | Yujuan Jiao, Jinmiao Yang, Hang Zhang . Traveling wave solutions to a cubic predator-prey diffusion model with stage structure for the prey. AIMS Mathematics, 2022, 7(9): 16261-16277. doi: 10.3934/math.2022888 |
[3] | Baojian Hong, Jinghan Wang, Chen Li . Analytical solutions to a class of fractional coupled nonlinear Schrödinger equations via Laplace-HPM technique. AIMS Mathematics, 2023, 8(7): 15670-15688. doi: 10.3934/math.2023800 |
[4] | Mahmoud Soliman, Hamdy M. Ahmed, Niveen Badra, Taher A. Nofal, Islam Samir . Highly dispersive gap solitons for conformable fractional model in optical fibers with dispersive reflectivity solutions using the modified extended direct algebraic method. AIMS Mathematics, 2024, 9(9): 25205-25222. doi: 10.3934/math.20241229 |
[5] | Bruce Brewer, Jake Daniels, Nghiem V. Nguyen . Exact Jacobi elliptic solutions of some models for the interaction of long and short waves. AIMS Mathematics, 2024, 9(2): 2854-2873. doi: 10.3934/math.2024141 |
[6] | Khalid Khan, Amir Ali, Manuel De la Sen, Muhammad Irfan . Localized modes in time-fractional modified coupled Korteweg-de Vries equation with singular and non-singular kernels. AIMS Mathematics, 2022, 7(2): 1580-1602. doi: 10.3934/math.2022092 |
[7] | Jin Hyuk Choi, Hyunsoo Kim . Exact traveling wave solutions of the stochastic Wick-type fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation. AIMS Mathematics, 2021, 6(4): 4053-4072. doi: 10.3934/math.2021240 |
[8] | Xijun Deng, Aiyong Chen . Traveling wave fronts in a single species model with cannibalism and strongly nonlocal effect. AIMS Mathematics, 2024, 9(10): 26688-26701. doi: 10.3934/math.20241298 |
[9] | M. TarikulIslam, M. AliAkbar, M. Abul Kalam Azad . Traveling wave solutions in closed form for some nonlinear fractional evolution equations related to conformable fractional derivative. AIMS Mathematics, 2018, 3(4): 625-646. doi: 10.3934/Math.2018.4.625 |
[10] | Ahmed Yahya, Adam Zakria, Ibrahim-Elkhalil Ahmed, Shams A. Ahmed, Husam E. Dargail, Abdelgabar Adam Hassan, Eshraga Salih . Viscoelastic Pasternak foundation analysis of a thermoelastic microbeam using Moore–Gibson–Thompson heat conduction under Klein–Gordon (KG) nonlocality. AIMS Mathematics, 2025, 10(5): 10340-10358. doi: 10.3934/math.2025471 |
The vertex arboricity a(G) of a graph G is the minimum number of colors required to color the vertices of G such that no cycle is monochromatic. The list vertex arboricity al(G) is the list version of this concept. In this paper, we prove that if G is a planar graph without 5-cycles intersecting with 6-cycles, then al(G)≤2.
Population dynamic models with age structure have a long history from articles of Lotka [16] in 1907 and Sharp and Lotka [24] in 1911. The first nonlinear age-structured model was proposed by Gurtin and MacCamy [18] in 1974. The theory of traveling wave solutions of reaction-diffusion systems has attracted much attention due to its significant nature in biology, chemistry, epidemiology and physics. Ducrot and Magal [3] researched the existence of traveling wave solutions for infection-age structured model with local diffusion by the methods of upper-lower solutions and Schauder's fixed point theorem. More recently, the traveling wave solutions of age-structured models received a lot of interest in the literature and we refer to [3,4,5,26,29] for more results about this topic.
The article is devoted to the study of traveling wave solutions of the nonlocal dispersal SIRS epidemic model with age structure
{∂∂tS(t,x)=d[J∗S(t,x)−S(t,x)]−S(t,x)∫A0β(a)i(t,a,x)da+δR(t,x),∂∂ti(t,a,x)+∂∂ai(t,a,x)=d[J∗i(t,a,x)−i(t,a,x)]−γi(t,a,x),∂∂tR(t,x)=d[J∗R(t,x)−R(t,x)]+γ∫A0i(t,a,x)da−δR(t,x),i(t,0,x)=S(t,x)∫A0β(a)i(t,a,x)da, | (1.1) |
where (t,a,x)∈D≜R+×(0,A]×Rand
(J∗S−S)(t,x)≜∫RJ(x−y)S(t,y)dy−S(t,x), |
(J∗i−i)(t,a,x)≜∫RJ(x−y)i(t,a,y)dy−i(t,a,x), |
(J∗R−R)(t,x)≜∫RJ(x−y)R(t,y)dy−R(t,x). |
S(t,x) and R(t,x) represent the densities of susceptible and removed individuals located at time t and position x and i(t,a,x) represents the density of infected individuals located at time t and location x with age a. And a⩾0 is the time since individuals were infected. Based on the actual situation, we assume that the age of infection is a constant. A∈(0,+∞) means the maximum age of infection. Therefore, we assume that the maximum infection age A is a sufficiently large constant and i(t,A,x)=0 throughout this paper. The parameter d>0 represents the diffusion rate of individuals. The function β(a) denotes the infection age-specific transmission rate. γ>0 is the recovery rate of the infected individuals and δ>0 is the loss of immunity rate. Since there are no birth or death rates in this model, N(t,x)=S(t,x)+∫A0i(t,a,x)da+R(t,x) is a constant for all t⩾0. Therefore, we assume that the initial value of (1.1) is S(0,x)+∫A0i(0,a,x)da+R(0,x)=S0. And in this whole paper we assume δ>γ and we mainly apply this assumption in Lemma 3.3. Moreover, J(x−y) denotes the probability of jumping from position y to position x and J∗S(t,x) denotes the total number of susceptible individuals located at time t, moving from the whole space to position x. In this paper, we give the following assumptions.
(A1) J∈C1(R), J(x)=J(−x)≥0, ∫RJ(x)dx=1 and J is compactly supported.
(A2) The map a→β(a) is almost everywhere bounded and belongs to L∞(0,A].
Kermack and McKendrick [14] first proposed a compartmental model to describe the spread of infectious diseases in 1927. Since then, the SIR epidemic model and its various adaptations have been widely studied in pathology[22,26,30,32,34]. In epidemiology, traveling wave solutions have attracted a great deal of attention, which denote propagation through space at a constant velocity. In order to control and prevent diseases, it is crucial to determine whether traveling wave solutions exist. Therefore, many studies utilize spatially and temporally correlated models to investigate the existence of traveling wave solutions[15,21,29,31,33].
Hosono and Ilyas [7] investigated the reaction-diffusion equation with the nonlinear reaction term. They verified the existence of a noncritical traveling wave solution. However, it is worth mentioning that the model involves the classical Laplacian diffusion in [7], which is also called the local diffusion. The local diffusion denotes that a population at position x will only be affected by variations in the population near position x, which is deficient in describing diffusion[1,6,10,12,17,21]. To overcome the limitations of diffusion of individuals, the Laplacian operator can be replaced by the nonlocal diffusion J∗u−u. Traveling wave solutions of infectious disease models with nonlocal diffusion terms have drawn great attention in recent years. For example, Yang et al. [33] studied the existence and nonexistence of traveling wave solutions in a Kermack-McKendrick epidemic model with nonlocal diffusion. They established a non-empty closed convex cone on a bounded closed interval and constructed appropriate upper and lower solutions and then obtained a nonlinear operator of an ordinary differential system. By applying Schauder's fixed point theorem, they proved the existence of nontrivial traveling wave solutions. Furthermore, Yang et al. [31,32] also investigated the traveling wave solutions of the SIR epidemic model with critical wave speed. After that, Ma and Yuan [17] studied the traveling wave solutions of the nonlocal dispersal SIRS model with spatio-temporal delay. Qiao et al. [21] researched the traveling wave solutions of the nonlocal dispersal SEIR model with standard incidence.
In addition, age is one of the important features in describing epidemics due to the fact that individuals of different ages may have different survivability and behaviors in their natural conditions. For example, hand-foot-mouth disease, chickenpox, measles, and influenza are all susceptible to occur in childhood. Novel coronavirus infections are universally susceptible to the population and can occur in all age groups, but are less common in children under three years of age. The disease is mainly concentrated in older people over the age of sixty-five and the virus is easily exacerbated. Therefore, age structure is a significant factor in studying the epidemic patterns of infectious diseases. In recent years, lots of age-structure models have been proposed in the epidemic spread (see [2,8,9,13,19,20,29]). For instance, Ducrot et al. [5] studied the existence of traveling wave solutions for multigroup age-structure models. At the same time, Ducrot and Magal [3] researched the existence of traveling wave solutions for the infection-age structured model with diffusion. In addition, Ducrot and Magal [4] also proposed the SI infectious disease model with external supplies and age structure. And age is the time since the individuals were infected. They proved the existence of traveling wave solutions by constructing suitable upper and lower solutions and applying Schauder's fixed point theorem. Meanwhile, they studied the convergence behavior at positive infinity by constructing a suitable Lyapunov functional. In addition, Tian and Guo [26] obtained the existence of traveling wave solutions of nonlocal dispersal Fisher-KPP model with age structure. In recent years, Kang and Ruan [11] have proposed an age-structure SIS infectious disease model with nonlocal diffusion
{(∂∂t+∂∂a)S(t,a,x)=d[LS](t,a,x)−(λ(t,a,x)+μ(a,x))S(t,a,x)+γ(a,x)I(t,a,x),(∂∂t+∂∂a)I(t,a,x)=d[LI](t,a,x)+λ(t,a,x)S(t,a,x)−(μ(a,x)+γ(a,x))I(t,a,x),S(t,0,x)=∫a+0∫Ωβ(a,x,y)S(t,a,y)dyda,t>0, | (1.2) |
with
I(t,0,x)=0,S(0,a,x)=S0(a,x),I(0,a,x)=I0(a,x), |
where a≥0, t≥0, x∈Ω⊂RN and the nonlocal operator L is defined by
[Lu](t,a,x):=(J∗u−u)(t,a,x)=∫ΩJ(x−y)u(t,a,y)dy−u(t,a,x). |
And μ(a,x) and γ(a,x) represent the mortality and recovery rates, respectively, for individuals of age a at location x. λ(t,a,x) is the infectivity of infected individuals to susceptible individuals of age a at time t in position x. β(a,x,y) indicates the birth rate of individuals of age a on position y giving birth to newborns at position x. They prove the existence of traveling wave solutions by constructing upper and lower solutions in [11]. We can refer to [15,22,23,25,28,34] for the relevant conclusions on the research of traveling wave solutions for many infectious disease models. As far as we know, there are few researches on the traveling wave solutions of epidemic models with age structure [4,29].
This paper is structured as follows. In Section 2, we give some key preliminaries of this paper. Then, we establish the existence of traveling wave solutions for c>c∗ in Section 3. Section 4 obtains the nonexistence of traveling wave solutions for 0<c<c∗ by using the Laplace transform. We draw the conclusions in Section 5.
In this section, we present some preliminaries and the characteristic equation.
Let N(t,x)=S(t,x)+∫A0i(t,a,x)da+R(t,x). The system (1.1) can be rewritten as the following system:
{∂∂tN(t,x)=d[J∗N(t,x)−N(t,x)],∂∂ti(t,a,x)+∂∂ai(t,a,x)=d[J∗i(t,a,x)−i(t,a,x)]−γi(t,a,x),∂∂tR(t,x)=d[J∗R(t,x)−R(t,x)]+γ∫A0i(t,a,x)da−δR(t,x),i(t,0,x)=(N(t,x)−∫A0i(t,a,x)da−R(t,x))∫A0β(a)i(t,a,x)da. | (2.1) |
The homogeneous system for (2.1) always exists a disease-free equilibrium (S0,0,0), where S0 denotes the density of susceptible individuals at the start of the infection. In addition, if the basic reproduction number R0=S0∫A0β(a)e−γada>1, the homogeneous system for (2.1) exists a unique endemic equilibrium (S0,i∗(a),R∗), where
i∗(a)=δγδ+γ(S0−1∫A0β(a)e−γada)e−γa1−e−γA and R∗=γδ+γ(S0−1∫A0β(a)e−γada). |
Our purpose in this paper is to discuss the traveling wave solutions connecting the disease-free equilibrium and the endemic equilibrium. More precisely, we aim to study the traveling wave solutions of system (2.1) with the form
(N(t,x),i(t,a,x),R(t,x))=(ˆN(z),ˆi(a,z),ˆR(z)), z=x+ct, | (2.2) |
where the parameter c>0 indicates the wave speed. Substituting (2.2) into (2.1) and replacing ˆN, ˆi, ˆR with N, i, R, we deduce the following system:
{cN′(z)=d[J∗N(z)−N(z)],c∂∂zi(a,z)+∂∂ai(a,z)=d[J∗i(a,z)−i(a,z)]−γi(a,z),cR′(z)=d[J∗R(z)−R(z)]+γ∫A0i(a,z)da−δR(z),i(0,z)=(N(z)−∫A0i(a,z)da−R(z))∫A0β(a)i(a,z)da. | (2.3) |
Thus, the solution of system (2.3) connecting (S0,0,0) and (S0,i∗(a),R∗) is a special solution, which satisfies the asymptotic boundary conditions
(N,i,R)(−∞)=(S0,0,0),(N,i,R)(+∞)=(S0,i∗(a),R∗). | (2.4) |
We intend to obtain the existence of traveling waves of system (2.1), which satisfies N(±∞)=S0. Since N(z) is a constant, system (2.3) can be simplified by removing the first equation to the following system:
{c∂∂zi(a,z)+∂∂ai(a,z)=d[J∗i(a,z)−i(a,z)]−γi(a,z),cR′(z)=d[J∗R(z)−R(z)]+γ∫A0i(a,z)da−δR(z),i(0,z)=(S0−∫A0i(a,z)da−R(z))∫A0β(a)i(a,z)da, | (2.5) |
which satisfies the asymptotic boundary conditions
(i,R)(−∞)=(0,0)=E0,(i,R)(+∞)=(i∗(a),R∗)=E∗. | (2.6) |
Assume that R0=S0∫A0β(a)e−γada>1. By linearizing the third equation of system (2.5) around the disease-free equilibrium (0, 0) and letting i(a,z)=eλzϕ(a), we can obtain the characteristic equation as follows:
F1(λ,c)=S0∫A0β(a)eg(λ,c)ada−1, | (2.7) |
where g(λ,c)=d(∫RJ(x)eλxdx−1)−γ−cλ and the function ϕ satisfies ϕ′(a)=g(λ,c)ϕ(a). Notice that
F1(0,c)=S0∫A0β(a)e−γada−1=R0−1>0,∂F1(λ,c)∂c=−λS0∫A0β(a)eg(λ,c)ada<0,∂F1(λ,c)∂λ|λ=0=−cS0∫A0β(a)e−γaada<0,∂2F1(λ,c)∂λ2=S0∫A0β(a)eg(λ,c)aa(a(gλ(λ,c))2+gλλ(λ,c))da>0. | (2.8) |
In order to construct R+(z), we consider the following function and study its properties:
F2(λ,c)=d∫RJ(x)e−λxdx−d−δ−cλ+γR∗M, |
where M is a positive constant to be determined and M can be taken to be e(λ1−λ3)z3 in Lemma 3.3. And λ1,λ3 can be seen in Lemma 2.1 and z3 can be obtained from R+(z). Notice that
F2(0,c)=−δ+γR∗M,∂F2(λ,c)∂λ|λ=0=−c<0,∂F2(λ,c)∂c=−λ<0,∂2F2(λ,c)∂λ2=d∫RJ(x)x2eλxdx>0. | (2.9) |
Based on (2.8) and (2.9), we can obtain the following lemma.
Lemma 2.1. If R0=S0∫A0β(a)e−γada>1, then there exist positive constants c∗i and λ∗i such that
Fi(λ∗i,c∗i)=0and∂Fi(λ,c)∂λ|(λ∗i,c∗i)=0,i=1,2. |
Furthermore,
(i) ● when c∈(0,c∗1), it is known that F1(λ,c)>0 holds for all λ≥0;
● when c>c∗1, F1(λ,c)=0 has two positive roots λ1(c) and λ2(c) satisfying
0<λ1(c)<λ∗1<λ2(c)<+∞. |
In addition, when c>c∗1, F1(λ,c) is less than zero in λ∈(λ1(c),λ2(c)) and greater than zero beyond [λ1(c),λ2(c)].
(ii) ● When c∈(0,c∗2), it is known that F2(λ,c)>0 holds for all λ≥0;
● when c>c∗2, F2(λ,c)=0 has two positive roots λ3(c) and λ4(c) satisfying
0<λ3(c)<λ∗2<λ4(c). |
In addition, when c>c∗2, F2(λ,c) is less than zero in λ∈(λ3(c),λ4(c)) and greater than zero beyond [λ3(c),λ4(c)].
Remark 2.1. When δ=0 and γ=0, we obtain c∗2=0, λ3=0. Therefore, we get c∗2<c∗1 and λ3<λ1 by choosing δ and γ small enough. Taking c∗=c∗1, we consider c∗ as the critical wave speed in this paper.
Remark 2.2. The function ϕ(a) satisfies the equation ϕ′(a)=g(λ,c)ϕ(a) for c>c∗. In the process of deriving the characteristic equation, it can be obtained that ϕ(0)=S0∫+∞0β(a)ϕ(a)da.
In this section, we assume c>c∗ and give the definitions of the upper solution and lower solution. To construct the suitable upper and lower solutions, we assume that δ and γ are small enough throughout this paper (see Remark 2.1 and (3.3)). Then, the existence of traveling wave solutions is proved by using Schauder's fixed point theorem.
Definition 3.1. A pair of the continuous functions Φ+=(i+(a,z),R+(z)) and Φ−=(i−(a,z),R−(z)) are called the upper solution and lower solution of system (2.5), if there exists a finite set S={Zi∈R:i=1,2,...,m} such that Φ′+,Φ′− exist and are bounded for z∈RnS, and satisfy
{c∂∂zi+(a,z)+∂∂ai+(a,z)⩾d[J∗i+(a,z)−i+(a,z)]−γi+(a,z),cR′+(z)⩾d[J∗R+(z)−R+(z)]+γ∫A0i+(a,z)da−δR+(z),i+(0,z)⩾(S0−∫A0i+(a,z)da−R−(z))∫A0β(a)i+(a,z)da, | (3.1) |
and
{c∂∂zi−(a,z)+∂∂ai−(a,z)⩽d[J∗i−(a,z)−i−(a,z)]−γi−(a,z),cR′−(z)⩽d[J∗R−(z)−R−(z)]+γ∫A0i−(a,z)da−δR−(z),i−(0,z)⩽(S0−∫A0i−(a,z)da−R+(z))∫A0β(a)i−(a,z)da, | (3.2) |
for z∈RnS, respectively.
Let ε>0 be small enough and ε1=ε+ε2. Next, we define two continuous functions Φ+=(i+(a,z),R+(z)) and Φ−=(i−(a,z),R−(z)) as follows:
i+(a,z)={eλ1zϕ(a), z⩽z1,i∗(a)+ε1e−λzϕ(a), z>z1,i−(a,z)={0,z⩽z2,i∗(a)−ε1e−λzϕ(a),z>z2,R+(z)={R∗eλ3z, z⩽z3,R∗+εe−λz, z>z3,R−(z)={0, z⩽z4,R∗−εe−λz, z>z4, | (3.3) |
where z1,z3>0, z2<z4<0 and λ>0 is a small enough constant.
Remark 3.1. By the expressions for i−(a,z) and R−(z), it follows that z2=−1λlnR∗ε−1λlnδe−γa(1+ε)ϕ(a) and z4=−1λlnR∗ε. And it is easy to hold that z2<z4, e.g., when γ=1A, where A is a sufficiently large constant.
Lemma 3.1. The function i+(a,z) satisfies
{c∂∂zi+(a,z)+∂∂ai+(a,z)⩾d[J∗i+(a,z)−i+(a,z)]−γi+(a,z),i+(0,z)⩾(S0−∫A0i+(a,z)da−R−(z))∫A0β(a)i+(a,z)da, | (3.4) |
for z≠z1.
Proof. For the first equation of system (3.4), when z<z1, i+(a,z)=eλ1zϕ(a). Thus, we just need to prove
cλ1eλ1zϕ(a)+eλ1zϕ′(a)⩾d(∫RJ(z−y)eλ1yϕ(a)dy)−(γ+d)eλ1zϕ(a). |
It is sufficient to prove
ϕ′(a)ϕ(a)⩾d∫RJ(x)eλ1xdx−γ−d−cλ1. |
Recalling Remark 2.2, the above inequality holds.
When z>z1, i+(a,z)=i∗(a)+ε1e−λzϕ(a), we need to prove
−cλε1e−λzϕ(a)+i∗′(a)+ε1e−λzϕ′(a)⩾d∫RJ(z−y)(i∗(a)+ε1e−λyϕ(a))dy−(d+γ)(i∗(a)+ε1e−λzϕ(a)). |
It is enough to verify
ϕ′(a)ϕ(a)⩾d∫RJ(x)e−λxdx−γ−d+cλ. |
That is to say,
g(λ1)⩾g(λ)+2cλ. | (3.5) |
When λ>0 is small enough, (3.5) is established.
For the second equation of the system (3.4), when z<z1, i+(a,z)=eλ1zϕ(a). Therefore, it is sufficient to demonstrate
eλ1zϕ(0)⩾(S0−∫A0i+(a,z)da−R−(z))eλ1z∫A0β(a)ϕ(a)da. | (3.6) |
Recalling Remark 2.2, inequality (3.6) holds.
When z>z1, i+(a,z)=i∗(a)+ε1e−λzϕ(a). We need to show
i∗(0)+ε1e−λzϕ(0)⩾(S0−∫A0i+(a,z)da−R−(z))∫A0β(a)(i∗(a)+ε1e−λzϕ(a))da. | (3.7) |
Due to ϕ(0)=S0∫A0β(a)ϕ(a)da and i∗(0)>(S0−∫A0i+(a,z)da−R−(z))∫A0β(a)i∗(a)da, it is evident that (3.7) is true.
Lemma 3.2. The function i−(a,z) satisfies
{c∂∂zi−(a,z)+∂∂ai−(a,z)⩽d[J∗i−(a,z)−i−(a,z)]−γi−(a,z),i−(0,z)⩽(S0−∫A0i−(a,z)da−R+(z))∫A0β(a)i−(a,z)da, | (3.8) |
for z≠z2.
Proof. For the first equation of system (3.8), when z<z2, i−(a,z)=0. Therefore, it is obviously true.
When z>z2, i−(a,z)=i∗(a)−ε1e−λzϕ(a). We need to show
cλε1e−λzϕ(a)+i∗′(a)−ε1e−λzϕ′(a)⩽d∫RJ(z−y)(i∗(a)−ε1e−λyϕ(a))dy−(d+γ)(i∗(a)−ε1e−λzϕ(a)). | (3.9) |
(3.9) can be simplified to
ϕ′(a)ϕ(a)⩾d∫RJ(x)e−λxdx−γ−d+cλ. |
That is to say,
g(λ1)⩾g(λ)+2cλ. | (3.10) |
When λ>0 is small enough, (3.10) holds.
For the second equation of system (3.8), when z<z2, i−(a,z)=0. It is clearly true.
When z>z2, i−(a,z)=i∗(a)−ε1e−λzϕ(a). We need to show
i∗(0)−ε1e−λzϕ(0)⩽(S0−∫A0(i∗(a)−ε1e−λzϕ(a))da−R+(z))∫A0β(a)(i∗(a)−ε1e−λzϕ(a))da. | (3.11) |
Next, we need to discuss the two cases of R+(z) in (3.11).
When z2<z<z3, R+(z)=R∗eλ3z. It is sufficient to ensure that
i∗(0)−ε1e−λzϕ(0)≤(S0−∫A0i∗(a)da+ε1e−λz∫A0ϕ(a)da−R∗eλ3z)(∫A0β(a)i∗(a)da−ε1e−λz∫+∞0β(a)ϕ(a)da). | (3.12) |
(3.12) can be simplified to
(ε1e−λz∫A0ϕ(a)da−R∗eλ3z+R∗)∫A0β(a)i∗(a)da+(∫A0i∗(a)da−ε1e−λz∫A0ϕ(a)da+R∗eλ3z)ε1e−λz∫A0β(a)ϕ(a)da⩾0. | (3.13) |
Under the conditions that
ε1e−λz3∫A0ϕ(a)da−R∗eλ3z3+R∗⩾0 | (3.14) |
and
∫A0i∗(a)da−ε1e−λz∫A0ϕ(a)da+R∗eλ3z⩾0 | (3.15) |
the inequality (3.13) holds. In addition, (3.14) is equal to
(1+ε)∫A0ϕ(a)da⩾1. | (3.16) |
We can choose ϕ(0)=g(λ1)eg(λ1)A−1 and ε1 small enough to guarantee (3.15) and (3.16) hold.
When z>z3, R+(z)=R∗+εe−λz. It is to prove that
i∗(0)−ε1e−λzϕ(0)≤(S0−∫A0i∗(a)da−R∗+ε1e−λz∫A0ϕ(a)da−εe−λz)(∫A0β(a)i∗(a)da−ε1e−λz∫A0β(a)ϕ(a)da). | (3.17) |
(3.17) can be simplified to
εe−λz((1+ε)∫A0ϕ(a)da−1)∫A0β(a)i∗(a)da+(∫A0i∗(a)da+R∗−ε1e−λz∫A0ϕ(a)da+εe−λz)(ε1e−λz∫A0β(a)ϕ(a)da)⩾0. | (3.18) |
Due to ϕ(0)=g(λ1)eg(λ1)A−1, then (1+ε)∫A0ϕ(a)da⩾1 holds. Since ε1 is small enough, then
∫A0i∗(a)da+R∗−ε1e−λz∫A0ϕ(a)da+εe−λz⩾0 | (3.19) |
is true. Therefore, (3.17) holds.
Lemma 3.3. The function R+(z) satisfies
cR′+(z)⩾d[J∗R+(z)−R+(z)]+γ∫A0i+(a,z)da−δR+(z), | (3.20) |
for z≠z3.
Proof. When z<z3, R+(z)=R∗eλ3z. It is necessary to demonstrate
cλ3R∗eλ3z⩾d∫RJ(z−y)R∗eλ3ydy−(d+δ)R∗eλ3z+γ∫A0i+(a,z)da. |
That is to say, we only need to prove this simplified expression as follows:
d∫RJ(z−y)eλ3(y−z)dy−d−δ+γR∗e−λ3z∫A0i+(a,z)da−cλ3⩽0. | (3.21) |
At this point, (3.21) can be organized in the following inequalities:
F2(λ3,c)−γR∗(M−e−λ3z∫A0i+(a,z)da)⩽0. |
Since F2(λ3,c)=0, we need to demonstrate
M⩾e−λ3z∫A0i+(a,z)da. | (3.22) |
Next, we will consider different cases for i+(a,z).
If z<z1, then i+(a,z)=eλ1zϕ(a). Noticing that ∫A0ϕ(a)da=1, we can choose M=e(λ1−λ3)z3 to ensure that (3.22) holds.
If z>z1, then i+(a,z)=i∗(a)+ε1e−λzϕ(a). We just need to verify
e(λ1−λ3)z3⩾e−λ3z∫A0(i∗(a)+ε1e−λzϕ(a))da. | (3.23) |
Simplifying (3.23), we obtain
e(λ1−λ3)z3eλ3z1⩾∫A0i∗(a)da+ε1e−λz1. | (3.24) |
Since λ3<λ1 and z1<z<z3, it follows that the left side of (3.24) is greater than 1. When ε1>0 is sufficiently small, ε1e−λz1 tends to zero. By taking S0⩽1 so that ∫A0i∗(a)da⩽1, we have that (3.24) holds.
When z>z3, R+(z)=R∗+εe−λz, i+(a,z)⩽i∗(a)+ε1e−λ(z∨z1)ϕ(a), where z∨z1:=max{z,z1}. It is to prove that
−cλεe−λz⩾d∫RJ(z−y)(R∗+εe−λy)dy−(d+δ)(R∗+εe−λz)+γ∫A0(i∗(a)+ε1e−λ(z∨z1)ϕ(a))da. | (3.25) |
Due to
|d∫RJ(x)e−λxdx−d+cλ|→0asλ→0, |
then there exists δ1>0, for any λ∈(0,δ1),
|d∫RJ(x)e−λxdx−d+cλ|<η1. | (3.26) |
Thus, (3.25) can be simplified to
d∫RJ(x)e−λxdx−d+cλ−δ+γ(1+ε)e−λ(z∨z1−z)∫A0ϕ(a)da≤η1−δ+γ(1+ε). | (3.27) |
We can choose η1=δ−γ(1+ε), then (3.25) holds.
Lemma 3.4. The function R−(z) satisfies
cR′−(z)⩽d[J∗R−(z)−R−(z)]+γ∫A0i−(a,z)da−δR−(z), | (3.28) |
for z≠z4.
Proof. When z<z4, R−(z)=0. Thus, it is obviously true.
When z>z4, R−(z)=R∗−εe−λz. We need to show
cελe−λz≤d∫RJ(z−y)(R∗−εe−λy)dy−(d+δ)(R∗−εe−λz)+γ∫A0i−(a,z)da. | (3.29) |
Since z2<z4 when A is sufficiently large, (3.29) reduces to
d∫RJ(x)e−λxdx−d+cλ⩽δ−γ(1+ε). | (3.30) |
Due to
|d∫RJ(x)e−λxdx−d+cλ|→0asλ→0, |
then there exists δ1>0, for any λ∈(0,δ1),
|d∫RJ(x)e−λxdx−d+cλ|<η1. | (3.31) |
We choose
η1=δ−γ(1+ε), | (3.32) |
then (3.29) holds true.
For any b>0, define the bounded closed convex set
Cb={i0∈C([−b,b]):i−(0,z)⩽i0(z)⩽i+(0,z)}. |
Next, we study the operator T:Cb→C([−b,b]) given by
T(i0)(z)=(S0−∫A0i(a,z)da−R(z))∫A0β(a)i(a,z)da, |
where i is the solution of the problem
{c∂∂zi(a,z)+∂∂ai(a,z)=d[J∗i(a,z)−i(a,z)]−γi(a,z),a>0,z∈(−b,b),i(0,z)=i0(z),z∈(−b,b),i(a,±b)=i−(a,±b),a>0, | (3.33) |
while R is the solution of the problem
{cR′(z)=d[J∗R(z)−R(z)]+γ∫A0i(a,z)da−δR(z),z∈(−b,b),R(±b)=R−(±b). | (3.34) |
Therefore, the existence of solutions to systems (3.33) and (3.34) can be converted to the existence of a fixed point for the operator T.
Now, we can present the main result of the existence of traveling waves in this section.
Theorem 3.1. (Existence of traveling waves) We assume that R0>1, for any c>c∗, system (2.5) has a solution that connect the disease-free equilibrium E0 and the endemic equilibrium E∗.
Proof. From the definition of Cb, it can be shown that Cb is closed, convex and bounded. Moreover, by the similar argument in [17], we have operator T is completely continuous and T(Cb)⊂Cb. Using the Schauder's fixed theorem, T has a fixed point i0. Let (ib(a,z),Rb(z))(z∈(−b,b)) be the solution of systems (3.33) and (3.34) for any fixed b>0. To obtain the existence of traveling wave solutions in R, we choose an increasing sequence {bn}+∞n=1 such that bn>max{z1,z3} and limn→+∞bn=+∞. By similar arguments in [3, Section 4.4] and [4, Proposition 2.5], for the sequence (ibn,Rbn), we can extract a subsequence by a standard diagonal extract argument, denoted by {ibnk}k∈N, {Rbnk}k∈N, which tend towards functions (i,R) in the following topologies ibnk→i and Rbnk→R as k→+∞ uniformly on every bounded closed interval and pointwise on R for any given a>0. Due to the fact that J is compactly supported, applying the Lebesgue dominated convergence theorem, we obtain the following results:
limk→+∞∫+∞−∞J(y)ibnk(a,z−y)dy=∫+∞−∞J(y)i(a,z−y)dy=J∗i(a,z) |
and
limk→+∞∫+∞−∞J(y)Rbnk(z−y)dy=∫+∞−∞J(y)R(z−y)dy=J∗R(z) |
for any z∈R. Therefore, we have that i(a,z) and R(z) satisfy system (2.5). Note the fact that
i−(a,z)⩽i(a,z)⩽i+(a,z), R−(z)⩽R(z)⩽R+(z),∀z∈R, | (3.35) |
and
limz→−∞(i−(a,z),R−(z))=(0,0),limz→+∞(i−(a,z),R−(z))=(i∗,R∗),limz→−∞(i+(a,z),R+(z))=(0,0),limz→+∞(i+(a,z),R+(z))=(i∗,R∗), | (3.36) |
we obtain
limz→−∞(i(a,z),R(z))=(0,0),limz→+∞(i(a,z),R(z))=(i∗,R∗). | (3.37) |
Thus, (i(a,z),R(z)) satisfies the asymptotic boundary conditions (2.6). The proof is completed.
In this section, we mainly focus on the nonexistence of traveling waves when 0<c<c∗ with R0>1 by using the Laplace transform.
Theorem 4.1. (Nonexistence of traveling waves) Assume that R0>1, for any speed c∈(0,c∗), there exist no nontrivial traveling wave solutions (i(a,z),R(z)) of system (2.5) satisfying (2.6).
Proof. By contradiction, we assume that there exists a nontrivial traveling wave solution (i(a,z),R(z)) of system (2.5) that satisfies
(i,R)(−∞)=(0,0),(i,R)(+∞)=(i∗,R∗). |
Due to R0=S0∫A0β(a)e−γada>1 and (i,R)(−∞)=(0,0), there exists ˆz<0 such that S0−I(z)−R(z)>S02+12∫A0β(a)e−γada for any z<ˆz, where I(z)=∫A0i(a,z)da. Integrating the first equation of (2.5) with respect to a from 0 to A, we have
cI′(z)=d(J∗I(z)−I(z))−γI(z)+(S0−I(z)−R(z))∫A0β(a)i(a,z)da⩾d(J∗I(z)−I(z))−γI(z)+(S02+12∫A0β(a)e−γada)β1I(z), | (4.1) |
where β1=infa∈(0,A]β(a). For any z<ˆz, let H(z)=∫z−∞I(s)ds. Integrating two sides of inequality (4.1) from −∞ to z, we obtain
cI(z)⩾d∫z−∞(J∗I(s)−I(s))ds+[(S02+12∫A0β(a)e−γada)β1−γ]H(z). | (4.2) |
By applying Fubini theorem, it holds that
∫z−∞J∗I(s)ds=∫+∞−∞∫z−∞J(x)I(s−x)dsdx=∫+∞−∞J(x)∫z−x−∞I(s)dsdx=J∗H(z). | (4.3) |
Substituting (4.3) into (4.2), we get
cI(z)⩾d(J∗H(z)−H(z))+[(S02+12∫A0β(a)e−γada)β1−γ]H(z). | (4.4) |
Thanks to
∫z−∞(J∗H(s)−H(s))ds=∫z−∞∫+∞−∞J(x)(H(s−x)−H(s))dxds=∫z−∞∫+∞−∞(−x)J(x)∫10H′(s−θx)dθdxds=∫10∫+∞−∞(−x)J(x)∫z−∞I(s−θx)dsdxdθ=∫+∞−∞(−x)J(x)∫10H(z−θx)dθdx, |
we have J∗H(z)−H(z) is integrable on (−∞,z] for any z∈R. From Eq (4.4), we obtain that H(z) is integrable on (−∞,z] for any z∈R. Then integrating both sides of inequality (4.4) from −∞ to z with z⩽ˆz, we have
[(S02+12∫A0β(a)e−γada)β1−γ]∫z−∞H(s)ds⩽cH(z)+d∫+∞−∞xJ(x)∫10H(z−θx)dθdx. |
Due to xH(z−θx) is non-increasing for θ∈[0,1] with any fixed z∈R. We obtain
[(S02+12∫A0β(a)e−γada)β1−γ]∫z−∞H(s)ds⩽(c+d∫+∞−∞xJ(x)dx)H(z). |
By the property that J is an even function, we have ∫+∞−∞xJ(x)dx=0 holds. Thus, for any z⩽ˆz, one gets
[(S02+12∫A0β(a)e−γada)β1−γ]∫+∞0H(z−s)ds⩽cH(z). |
Since H(z) is increasing with respect to z, there exists some τ>0 such that
[(S02+12∫A0β(a)e−γada)β1−γ]τH(z−τ)⩽cH(z). |
Hence, there exists a constant τ0>0 large enough and some ν∈(0,1) such that H(z−τ0)⩽νH(z) for each z⩽ˆz. Set Q(x)=H(x)e−μ1x and μ1=1τ0ln1ν, then
Q(z−τ0)=H(z−τ0)e−μ1(z−τ0)⩽νH(z)e−μ1(z−τ0)=Q(z). |
Therefore, there exists Q0>0 such that Q(z)⩽Q0 for any z⩽ˆz, which implies
H(z)⩽Q0eμ1zfor any z⩽ˆz. |
It is noticed that
cI′(z)⩽d(J∗I(z)−I(z))+(S0‖ |
then there exists P_1 > 0 , such that I(z) \leqslant P_1 \mathrm{e}^{\mu_1 z} for any z \leqslant \hat{z} . Due to I(z) is bounded, it is possible to obtain
\begin{equation} \sup _{z \in \mathbb{R}}\left\{I(z)\mathrm{e}^{-\mu_1 z}\right\} < +\infty \quad \text{and} \quad \sup _{z \in \mathbb{R}}\left\{I^{\prime}(z) \mathrm{e}^{-\mu_1 z}\right\} < +\infty. \end{equation} | (4.5) |
By the same process, we have
\begin{equation} \sup _{z \in \mathbb{R}}\left\{R(z)\mathrm{e}^{-\mu_2 z}\right\} < +\infty \quad \text{and} \quad \sup _{z \in \mathbb{R}}\left\{R^{\prime}(z) \mathrm{e}^{-\mu_2 z}\right\} < +\infty \quad \text{for some} \quad \mu_2 > 0. \end{equation} | (4.6) |
According to (4.5) and (4.6), it can be obtained that
\begin{equation*} \sup _{z \in \mathbb{R}} \left\{\mathrm{e}^{-\mu_0 z}(I(z)+R(z))\right\} < +\infty \end{equation*} |
and
\begin{equation*} \sup _{z \in \mathbb{R}}\left\{\mathrm{e}^{-\mu_0 z}\int_0^A \beta(a) i(a, z) \mathrm{d}a\right\} < \sup _{z \in \mathbb{R}}\left\{\|\beta\|_{\infty} \mathrm{e}^{-\mu_0z} I(z) \right\} < +\infty, \end{equation*} |
where \mu_0 = \min \left\{\mu_1, \mu_2\right\}. Therefore, we obtain
\begin{equation} \sup _{z \in \mathbb{R}}\left\{\mathrm{e}^{-2\mu_0 z } \left(I(z)+R(z)\right)\int_0^A \beta(a) i(a, z) \mathrm{d}a\right\} < +\infty. \end{equation} | (4.7) |
Next, we define
\begin{equation*} \tilde{I}(z) = \int_{0}^{A} \beta(a)i(a, z)\mathrm{d}a \end{equation*} |
and a two-sided Laplace transform of I(\cdot) by
\begin{equation*} L_{I}(\lambda) = \int_{-\infty}^{+\infty} \mathrm{e}^{-\lambda z} I(z) \mathrm{d}z \end{equation*} |
and a two-sided Laplace transform of \tilde{I}(\cdot) by
\begin{equation*} L_{\tilde{I}}(\lambda) = \int_{-\infty}^{+\infty} \mathrm{e}^{-\lambda z} \int_{0}^{A} \beta(a)i(a, z)\mathrm{d}a \mathrm{d}z \end{equation*} |
for \lambda \in \mathbb{C} with 0 < \operatorname{Re} \lambda < \mu_0 .
Then, in view of
\begin{equation} d(J * I(z)-I(z))-cI^{\prime}(z)-\gamma I(z)+S_0 \int_0^A \beta(a) i(a, z)\mathrm{d}a = \left(I(z)+R(z)\right) \int_0^A \beta(a) i(a, z)\mathrm{d}a \end{equation} | (4.8) |
and
\begin{equation*} \int_{-\infty}^{+\infty} \mathrm{e}^{-\lambda z} J * I(z) \mathrm{d} z = \int_{-\infty}^{+\infty} J(x) \mathrm{e}^{-\lambda x} L_{I}(\lambda)\mathrm{d}x, \end{equation*} |
we have
\begin{equation} \begin{aligned} g(\lambda, c)L_{I}(\lambda)+S_0 L_{\tilde{I}}(\lambda) = &\int_{-\infty}^{+\infty} \mathrm{e}^{-\lambda z} \left(I(z)+R(z)\right)\int_0^A \beta(a) i(a, z) \mathrm{d}a \mathrm{d} z. \end{aligned} \end{equation} | (4.9) |
By the property of the Laplace transform (see[27]), it follows that either there exists a real number \hat{\lambda} such that L_{I}(\lambda) and L_{\tilde{I}}(\lambda) are analytic for \lambda \in \mathbb{C} with 0 < \operatorname{Re} \lambda < \hat{\lambda} and \lambda = \hat{\lambda} is singular point of L_{I}(\lambda) and L_{\tilde{I}}(\lambda) , or L_{I}(\lambda) and L_{\tilde{I}}(\lambda) are well defined for \lambda \in \mathbb{C} with \operatorname{Re} \lambda > 0 . Indeed, motivated by Zhou, Xu, Wei et al. [35, Section 3], we can obtain that the right-hand side of Eq (4.9) is well defined for \lambda \in \mathbb{C} with 0 < \operatorname{Re} \lambda < 2 \mu_0 according to (4.7). One can obtain that L_{I}(\lambda) and L_{\tilde{I}}(\lambda) are well defined with \operatorname{Re} \lambda > 0 . Nevertheless, (4.9) can be rewritten as
\begin{equation} \int_{-\infty}^{+\infty} \mathrm{e}^{-\lambda z} \left[g(\lambda, c)I(z)+S_0 \int_0^A \beta(a) i(a, z)\mathrm{d}a-\left(I(z)+R(z)\right)\int_0^A \beta(a) i(a, z) \mathrm{d}a\right]\mathrm{d}z = 0. \end{equation} | (4.10) |
It is evident from the definition of g(\lambda, c) and Lemma 2.1 that g(\lambda, c) \rightarrow +\infty as \lambda \rightarrow +\infty . This is a contradiction to Eq (4.10) and we complete the proof.
In this paper, we propose and consider a nonlocal dispersal SIRS infectious disease model with age structure that has practical significance. Due to the consideration of age structure, the investigation of traveling wave problem becomes more complicated. In addition, the epidemic system is non-monotone such that the theory for monotone semiflow can not be applied. To overcome these difficulties, we obtain the existence of traveling wave solutions by constructing appropriate upper and lower solutions and applying Schauder fixed point theorem. In fact, the solution of system (2.5) can be derived when the basic reproduction number R_0 > 1 and the wave speed c > c^* . That is to say, infectious diseases can spread among populations when R_0 > 1 and c > c^* . Furthermore, we prove the nonexistence of traveling wave solutions when 0 < c < c^{*} with R_0 > 1 by using the Laplace transform.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors would like to thank the anonymous referees for their careful reading and valuable comments. This work is supported by the National Natural Science Foundation of China (Grant No. 12001502) and the 2023 Graduate Innovation Fund Project of China University of Geosciences, Beijing (Grant No. YB2023YC007).
All authors declare no conflict of interest that could affect the publication of this paper.
[1] |
O. V. Borodin, A. O. Ivanova, Planar graphs without 4-cycles adjacent to 3-cycles are list vertex 2-arborable, J. Graph Theory, 62 (2009), 234–240. doi: 10.1002/jgt.20394
![]() |
[2] |
G. Chartrand, H. V. Kronk, C. E. Wall, The point-arboricity of a graph, Israel J. Math., 6 (1968), 169–175. doi: 10.1007/BF02760181
![]() |
[3] |
M. Chen, A. Raspaud, W. F. Wang, Vertex-arboricity of planar graphs without intersecting triangles, Eur. J. Combin., 33 (2012), 905–923. doi: 10.1016/j.ejc.2011.09.017
![]() |
[4] |
M. Chen, L. Huang, W. F. Wang, List vertex-arboricity of toroidal graphs without 4-cycles adjacent to 3-cycles, Discrete Math., 339 (2016), 2526–2535. doi: 10.1016/j.disc.2016.04.010
![]() |
[5] |
I. Choi, H. H. Zhang, Vertex arboricity of toroidal graphs with a forbidden cycle, Discrete Math., 333 (2014), 101–105. doi: 10.1016/j.disc.2014.06.011
![]() |
[6] |
S. L. Hakimi, E. F. Schmeichel, A note on the vertex arboricity of a graph, SIAM J. Discrete Math., 2 (1989), 64–67. doi: 10.1137/0402007
![]() |
[7] |
D. J. Huang, W. C. Shiu, W. F. Wang, On the vertex-arboricity of planar graphs without 7-cycles, Discrete Math., 312 (2012), 2304–2315. doi: 10.1016/j.disc.2012.03.035
![]() |
[8] | L. Huang, M. Chen, W. F. Wang, Toroidal graphs without 3-cycles adjacent to 5-cycles have list vertex-arboricity at most 2, Int. J. Math. Stat., 16 (2015), 97–105. |
[9] |
W. F. Wang, L. Huang, M. Chen, List vertex-arboricity of planar graphs without intersecting 5-cycles, Acta Math. Appl. Sin. Engl. Ser., 36 (2020), 439–447. doi: 10.1007/s10255-020-0936-1
![]() |
[10] |
A. Raspaud, W. F. Wang, On the vertex-arboricty of planar graphs, Eur. J. Combin., 29 (2008), 1064–1075. doi: 10.1016/j.ejc.2007.11.022
![]() |
[11] |
Y. Q. Wang, M. Chen, W. F. Wang, A note on the list vertex arboricity of toroidal graphs, Discrete Math., 341 (2018), 3344–3347. doi: 10.1016/j.disc.2018.08.021
![]() |
[12] | Y. P. Yang, Y. Q. Wang, P. Wang, W. F. Wang, List vertex arboricity of planar graphs of diameter two, Adv. Math. (China), 50 (2021), 335–344. |
[13] |
A. F. Yang, J. J. Yuan, On the vertex arboricity of planar graphs of diameter two, Discrete Math., 307 (2007), 2438–2447. doi: 10.1016/j.disc.2006.10.017
![]() |
[14] | H. Zhang, On list vertex 2-arboricity of toroidal graphs without cycles of specific length, Bull. Iranian Math. Soc., 42 (2016), 1293–1303. |