# Global Dynamics of an SIS Model on Metapopulation Networks with Demographics

Authors: Maoxing Liu; Xinjie Fu; Jie Zhang; Donghua Zhao
Journal: Complexity
 (2021)

---

## Abstract

In this paper, we propose a susceptible-infected-susceptible (SIS) epidemic model with demographics on heterogeneous metapopulation networks. We analytically derive the basic reproduction number, which determines not only the existence of endemic equilibrium but also the global dynamics of the model. The model always has the disease-free equilibrium, which is globally asymptotically stable when the basic reproduction number is less than unity and otherwise unstable. We also provide sufficient conditions on the global stability of the unique endemic equilibrium. Numerical simulations are performed to illustrate the theoretical results and the effects of the connectivity and diffusion. Furthermore, we find that diffusion rates play an active role in controlling the spread of infectious diseases.


---

## Body

## 1. Introduction

In real life, the spread of infectious diseases is always affected by the flow of population. In order to study the impact of demographics and the flow of population, we will focus on an epidemic model with demographics on metapopulation networks. As Masuda [1] mentioned, underlying contact networks are considered to be static on the time scale of epidemics, whereas underlying humans’ social networks of prevailing infectious diseases such as influenza are presumably dynamic even during one day. The dynamics of networks are induced by diffusion of individuals among residences, workplaces, places for social activities, and so on. As a simple framework within which the role of spatial processes in disease transmission can be examined, metapopulation models have been introduced to describe epidemics and ecological invasion in such a situation [1–3]. A node in such a model represents a metapopulation or a habitat and not an individual. A link represents a physical pathway connecting a pair of metapopulations. Individuals travel from one node to another. Simultaneously, interactions between individuals, such as infection, can occur in each metapopulation.In recent years, the effect of diffusion on disease spread has caught the attention of many researchers [4–8]. Luca et al. [9] considered that metapopulation epidemic models can describe the spatial spread of an infectious disease through a spatially structured host population [10, 11]. These models consist of patches or subpopulations of the system, connected through a coupling process generally characterizing hosts’ mobility [9, 12]. When the diffusion was first considered, the reaction-diffusion equation was regarded as a two-step process, and the reaction-diffusion was assumed to occur simultaneously to obtain the reaction-diffusion processes. Saldaña [13] derived continuous-time equations governing the limit dynamics of discrete-time reaction-diffusion processes defined on heterogeneous metapopulations, and the spread of infectious diseases under two different transmission mechanisms was considered. Juher et al. [14] proposed a system of continuous-time equations to analyze the spread of infectious diseases, which is based on the limited transmission and nonlimited transmission of diseases. Compared with reference [13], the authors studied the stability of endemic equilibrium in uncorrelated networks with nonlimited transmission. Therefore, a certain amount of diffusion is needed for the spread of diseases in metapopulation situation. After that, Masuda [1] investigated the effect of diffusion rate on the spread of epidemics in metapopulation networks. It can be seen from the above work that diffusion can change the epidemic threshold in many heterogeneous networks [1].Based on previous work, we will investigate an epidemic model with demographics in metapopulation networks. Here the structure of patches (nodes) is made up of the connectivity (degree) distribution. Moreover, each patch contains two types of individuals: susceptible and infected individuals. Within each patch, transmission or recovery (reaction processes) occurred between individuals of different types. The diffusion of individuals takes place at once among patches (diffusion process) at constant rates. In fact, reaction and diffusion processes can be considered to take place simultaneously. With this idea, many epidemic models have been proposed to describe the dynamics of disease spread among patches. The model of this paper is a description of an SIS metapopulation model at the node level with birth and death. Based on the model in [1], we consider the model with node-based birth and death. Wang et al. [15] made use of a mean-field approximation, and this model has birth and death, but its birth and death do not change with the change of degree. In some models, in order to facilitate calculation, the elements of the connectivity matrix A are given by aii′=ipi′/i. In this paper, we will analyze the model without using the mean-field approximation and transforming the elements of the adjacency matrix, and we study it in an arbitrary network.The remaining of this paper is organized as follows. In Section2, we firstly present the model and then study the existence of equilibria and their stability. There exist threshold dynamics that are determined by the basic reproduction number R0. The theoretical results and the effects of connectivity and diffusion are demonstrated by numerical simulations in Section 3. The paper concludes with a summary and a brief discussion in last section.


## 2. Main Results

### 2.1. The Model

For a metapopulation network, each node stands for a metapopulation or patch. LetN be the number of nodes; the N×N adjacency matrix, denoted by A, is defined by(1)Aji=0,ifj=iorjandiare not adjacent whenj≠i,1,ifj≠iandjandiare adjacent.Clearly,A is symmetric. The model is an extension of the following one mentioned in Masuda [1]:(2)dρS,idt=−βρS,iρI,i+γρI,i−DSρS,i+DS∑j=1NAjikjρS,j,dρI,idt=βρS,iρI,i−γρI,i−DIρI,i+DI∑j=1NAjikjρI,j,where ρS,it and ρI,it represent the number of susceptible and infected individuals at node i and at time t, respectively, DS and DI are the diffusion rates for the susceptible and infected individuals, respectively, β is the transmission rate, γ is the recovery rate, and ki is the degree of the node i. For (2), it is shown that diffusion increases the epidemic threshold of the SIS dynamics in arbitrary heterogeneous networks.In this paper, we incorporate demographics into (2), and the model is as follows:(3)dρS,idt=bi−βρS,iρI,i+γρI,i−DSρS,i+DS∑j=1NAjikjρS,j−μρS,i,dρI,idt=βρS,iρI,i−γρI,i−DIρI,i+DI∑j=1NAjikjρI,j−μρI,i,where bi is the birth rate at node i and μ is the natural death rate of individuals (whose meaning is different from that in (2)).It is easy to see that (3) is well posed, that is, for any given nonnegative initial condition, it has a unique and global solution through it. Moreover, the solution is also nonnegative. Let b¯=∑i=1Nbi and H=∑i=1NρS,i+ρI,i. Adding up all equations of (3) yields(4)dHtdt=b¯−μHt,which implies that limsupt⟶∞Ht≤b¯/μ. Thus, the feasible region of (3) can be chosen as(5)Γ=ρS,1,…,ρS,N,ρI,1,…,ρI,N∈R+2N|∑i=1NρS,i+ρI,i≤b¯μ.It can be verified thatΓ is attractive and positively invariant with respect to (3) (see [16]). Let Γ∘ denote the interior of Γ and ∂Γ denote the boundary of Γ.


### 2.2. Existence of Disease-Free Equilibria

We start the analysis of (3) with the existence of disease-free equilibria. At a disease-free equilibrium of (3), we have the following linear system [17]:(6)bi−DSρS,i+DS∑j=1NAjikjρS,j−μρS,i=0,i=1,2,…,N,which can be written in the matrix form(7)L¯S=B¯,where(8)L¯=DS+μ−DSA21k2⋯−DSAN1kN−DSA12k1DS+μ⋯−DSAN2kN⋮⋮⋱⋮−DSA1Nk1−DSA2Nk2⋯DS+μ,S=ρS,1,ρS,2,…,ρS,NT,B¯=b1,b2,…,bNT.AsL¯ is a column diagonal dominant matrix with positive diagonal elements and negative off-diagonal elements, L¯ is nonsingular and any element of L¯−1 is positive [18]. Thus, (6) has a unique positive solution, which has the following result.Proposition 1.
System (3) always has a unique disease-free equilibrium, denoted by E0=ρS,10,…,ρS,N0,0,…,0.Since the adjacency matrixAji is randomly given, the disease-free equilibrium E0 given in Proposition 1 does not have an explicit expression. Thus, we calculate it in specific network: globally coupled network. In a globally coupled network,(9)A=01⋯110⋯1⋮⋮⋱⋮11⋯0.Then,(10)L¯=DS+μ−DSN−1⋯−DSN−1−DSN−1DS+μ⋯−DSN−1⋮⋮⋱⋮−DSN−1−DSN−1⋯DS+μ.By complex calculations, we get(11)L¯−1=N−1μ+DSNμDS+μ2N−1DSNμDS+μ2N−1⋯DSNμDS+μ2N−1DSNμDS+μ2N−1N−1μ+DSNμDS+μ2N−1⋯DSNμDS+μ2N−1⋮⋮⋱⋮DSNμDS+μ2N−1DSNμDS+μ2N−1⋯N−1μ+DSNμDS+μ2N−1,and hence(12)ρS,i0=N−1μ+DSbi+DS−bi+∑j=1NbjNμDS+μ2N−1,fori=1,2,…,N.In particular, ifb1=b2=⋯=bN=b, then ρS,i0=μ/b and the disease-free equilibrium is E0=μ/b,…,μ/b︸N,0,…,0.


### 2.3. Stability of the Disease-Free Equilibrium

Before considering the stability of the disease-free equilibrium of (3), we introduce the basic reproduction number.Linearizing (3) at the disease-free equilibrium E0, we can get(13)dItdt=MIt,where It=ρI,1t,ρI,2,…,ρI,NtT and(14)M=βρS,10−γ−μ−DIDIA21k2⋯DIAN1kNDIA12k1βρS,20−γ−μ−DI⋯DIAN2kN⋮⋮⋱⋮DIA1Nk1DIA2Nk2⋯βρS,N0−γ−μ−DI.Then,sM is a simple eigenvalue of M with a positive eigenvector. Let FiI denote the rate of appearance of new infections in node i, Vi+I be the transfer rate of infected individuals into node i by all other means, and Vi−I be the transfer rate of infected individuals out of node i. Then, we can choose(15)F=βρS,100⋯00βρS,20⋯0⋮⋮⋱⋮00⋯βρS,N0,(16)V=V1−I−V1+IV2−I−V2+I⋮VN−I−VN+I=μ+DI+γ−DIA21k2⋯−DIAN1kN−DIA12k1μ+DI+γ⋯−DIAN2kN⋮⋮⋱⋮−DIA1Nk1−DIA2Nk2⋯μ+DI+γ.We callFV−1 the next generation matrix for (3), and the basic reproduction number R0 is defined as(17)R0=ρFV−1,where ρFV−1 denotes the spectral radius of the matrix FV−1 (see [19] for details). Note that R0 is simply the product of the infection rate and the mean duration of the infection [15]. In a globally coupled network, the basic reproduction number is R0=bβ/μγ+μ. Refer to [19] for Theorem 1.Theorem 1.
The following statements on (3) are true.(i)
R0>1 if and only if sM>0 and R0<1 if and only if sM<0.(ii)
IfR0<1, then the disease-free equilibrium E0 is locally asymptotically stable while it is unstable if R0>1.
Now, we establish the global stability of the disease-free equilibrium.Theorem 2.
AssumeR0<1. Suppose that A is irreducible. Then, the disease-free equilibrium E0 of (3) is globally asymptotically stable in R+2N.Proof.
AsΓ is an attracting and positively invariant set, we only need to show the global stability of E0 in Γ.
LetF and V be given by (15) and (16), respectively. Since all off-diagonal entries of V are nonpositive and the sum of the entries in each column of V is positive, V is a nonsingular M-matrix. Also, V−1 is irreducible. By the Perron–Frobenius theorem ([18], p. 27), the nonnegative irreducible matrix V−1F has a positive left eigenvector w1,w2,…,wN corresponding to the eigenvalue ρV−1F. Since F is a diagonal matrix, ρV−1F=ρFV−1. Consequently, we have(18)w1,w2,…,wNV−1F=R0w1,w2,…,wN,or(19)1R0w1,w2,…,wN=w1,w2,…,wNF−1V.
Fori=1,…,N, denote ci=wi/βρS,i0>0 and I=ρI,1,ρI,2,…,ρI,NT. Consider the Lyapunov function(20)L=∑i=1NciρI,i,and the derivative of L along the solutions of (3) is(21)dLdt|2.2=∑i=1NciβρS,iρI,i−γρI,i−DIρI,i+DI∑j=1NAjikjρI,j−μρI,i≤∑i=1NciβρS,i0−γ−μ−DIρI,i+DI∑j≠iNAjikjρI,j=w1βρS,10,w2βρS,20,…,wNβρS,N0F−VI=w1,w2,…,wN1−F−1VI=w1,w2,…,wN1−1R0I≤0.
Sinceci>0 for all i, dL/dt|2.2=0 implies that either ρS,i=ρS,i0 or ρI,i=0 for any 1≤i≤N. When ρS,i=ρS,i0, it follows from the first equation of (3) that(22)0=dρS,i0dt=bi−βρS,i0ρI,i+γρI,i−DSρS,i0+DS∑j=1NAjikjρS,j0−μρS,i0.
Comparing this with (6) gives ρI,i=0. Thus, we have shown that dL/dt|2.2=0 which implies that ρI,i=0 for all i. By LaSalle’s invariance principle, we see that E0 is globally asymptotically stable in Γ and hence in R+2N.
The following result tells us that whenR0>1, the disease is persistent.Proposition 2 (see [17]).
Suppose thatA is irreducible. If R0>1, then system (3) is uniformly persistent and there exists an endemic equilibrium E∗ in Γ∘.The proof of Proposition2 is similar to the proof of Proposition 3.3 in [20]. Using uniform persistence of (3) and uniform boundedness of solutions in Γ∘ we can obtain the existence of an equilibrium of (3) in Γ∘(see Theorem 2.8.6 in [21]). In Proposition 2, the assumption that the adjacency matrix A=Aij is irreducible is necessary. If A=0, then system (3) has an asymptotically stable boundary equilibrium when R0>1, and thus the system is not persistent [22].


### 2.4. The Global Stability of the Endemic Equilibrium

In this section, we provide sufficient conditions on the global stability of the endemic equilibrium. The method of proof is the graph-theoretical approach developed in [23–25].Theorem 3.
Assume thatR0>1. Let E∗=ρS,1∗,…,ρS,N∗,ρI,1∗,…,ρI,N∗ be an endemic equilibrium of (3). If A is irreducible and there exists λ>0 such that DS∑j=1NAji/kjρS,j∗=λDI∑j=1NAji/kjρI,j∗ for all 1≤i,j≤N, then the endemic equilibrium E∗ is globally asymptotically stable in Γ∘, and hence (3) only has a unique endemic equilibrium.Proof.
For each1≤i≤N, set(23)ViρS,i,ρI,i=ρS,i−ρS,i∗−ρS,i∗lnρS,iρS,i∗+ρI,i−ρI,i∗−ρI,i∗lnρI,iρI,i∗.
Noting(24)DSρS,i∗=bi−βρS,i∗ρI,i∗+γρI,i∗+DS∑j=1NAjikjρS,j∗−μρS,i∗,DI+μρI,i∗=βρS,i∗ρI,i∗+DI∑j=1NAjikjρI,j∗−γρI,i∗,we calculate the derivative of Vi along solutions of (3) as(25)dVidt|2.2=bi−DSρS,i+DS∑j=1NAjikjρS,j−μρS,i−biρS,i∗ρS,i+βρI,iρS,i∗−γρI,iρS,i∗ρS,i+DSρS,i∗−DS∑j=1NAjikjρS,jρS,i∗ρS,i+μρS,i∗−DIρI,i+DI∑j=1NAjikjρI,j−μρI,i−βρS,iρI,i∗+γρI,i∗+DIρI,i∗+μρI,i∗−DI∑j=1NAjikjρI,jρI,i∗ρI,i=bi1−ρS,iρS,i∗+lnρS,iρS,i∗+1−ρS,i∗ρS,i+lnρS,i∗ρS,i+DS∑j=1NAjikjρS,j∗1−ρS,i∗ρS,jρS,iρS,j∗+lnρS,i∗ρS,jρS,iρS,j∗+DS∑j=1NAjikjρS,j∗ρS,jρS,j∗+lnρS,j∗ρS,j−ρS,iρS,i∗−lnρS,i∗ρS,i+DI∑j=1NAjikjρI,j∗1−ρI,i∗ρI,jρI,iρI,j∗+lnρI,i∗ρI,jρI,iρI,j∗+DI∑j=1NAjikjρI,j∗ρI,jρI,j∗+lnρI,j∗ρI,j−ρI,iρI,i∗−lnρI,i∗ρI,i+γρI,i∗1−ρS,i∗ρI,iρS,iρI,i∗+lnρS,i∗ρI,iρS,iρI,i∗+γρI,i∗lnρI,i∗ρI,i−lnρS,i∗ρS,i+γρI,i1−ρI,i∗ρS,iρI,iρS,i∗+lnρI,i∗ρS,iρI,iρS,i∗+γρI,ilnρI,iρI,i∗−lnρS,iρS,i∗≤DS∑j=1NAjikjρS,j∗ρS,jρS,j∗+lnρS,j∗ρS,j−ρS,iρS,i∗−lnρS,i∗ρS,i+DI∑j=1NAjikjρI,j∗ρI,jρI,j∗+lnρI,j∗ρI,j−ρI,iρI,i∗−lnρI,i∗ρI,i+γρI,i∗lnρI,i∗ρI,i−lnρS,i∗ρS,i+γρI,ilnρI,iρI,i∗−lnρS,iρS,i∗=DS∑j=1NAjikjρS,j∗ρS,jρS,j∗+lnρS,j∗ρS,j−ρS,iρS,i∗−lnρS,i∗ρS,i+DI∑j=1NAjikjρI,j∗ρI,jρI,j∗+lnρI,j∗ρI,j−ρI,iρI,i∗−lnρI,i∗ρI,i=DI∑j=1NAjikjρI,j∗λρS,jρS,j∗+λlnρS,j∗ρS,j+ρI,jρI,j∗+lnρI,j∗ρI,j−DI∑j=1NAjikjρI,j∗λρS,iρS,i∗+λlnρS,i∗ρS,i+ρI,iρI,i∗+lnρI,i∗ρI,i=DI∑j=1NAjikjρI,j∗GjρS,j,ρI,j−GiρS,i,ρI,i,where(26)GiρS,i,ρI,i=λρS,iρS,i∗+λlnρS,i∗ρS,i+ρI,iρI,i∗+lnρI,i∗ρI,i.
Here we have used the fact that1−x+lnx≤0 for x>0, and the equality holds if and only if x=1.
Consider a weight matrixW=wij with entries wij=DIAji/kjρI,j∗ and denote the corresponding weighted digraph as G,W. Let li=∑T∈TiwT≥0 be the same as that given in (A.1) in the Appendix associated with G,W. Then, by (A.2), the following identity holds:(27)∑i=1Nli∑j=1NDIAjikjρI,j∗GjρS,j,ρI,j−GiρS,i,ρI,i=0.
Set(28)VρS,1,ρI,1,…,ρS,N,ρI,N=∑i=1NliViρS,i,ρI,i.
Then, we obtain(29)dVdt|2.2=∑i=1NlidVidt|2.2≤∑i=1Nli∑j=1NDIAjikjρI,j∗GjρS,j,ρI,j−GiρS,i,ρI,i=0,for all ρS,1,ρI,1,…,ρS,N,ρI,N∈Γ∘. Since A is irreducible, we know that li>0 for all i (see the Appendix), and thus dV/dt|2.2=0 implies that ρS,i=ρS,i∗ for all i. From the first equation of (3), we obtain(30)0=bi−βρS,i∗ρI,i+γρI,i−DSρS,i∗+DS∑j=1NAjikjρS,j∗−μρS,i∗,for i=1,…,N, which implies that ρI,i=ρI,i∗. This means that the largest invariant set in dV/dt|2.2=0 is the singleton E∗. By LaSalle’s invariance principle, E∗ is globally asymptotically stable in Γ∘.


## 2.1. The Model

For a metapopulation network, each node stands for a metapopulation or patch. LetN be the number of nodes; the N×N adjacency matrix, denoted by A, is defined by(1)Aji=0,ifj=iorjandiare not adjacent whenj≠i,1,ifj≠iandjandiare adjacent.Clearly,A is symmetric. The model is an extension of the following one mentioned in Masuda [1]:(2)dρS,idt=−βρS,iρI,i+γρI,i−DSρS,i+DS∑j=1NAjikjρS,j,dρI,idt=βρS,iρI,i−γρI,i−DIρI,i+DI∑j=1NAjikjρI,j,where ρS,it and ρI,it represent the number of susceptible and infected individuals at node i and at time t, respectively, DS and DI are the diffusion rates for the susceptible and infected individuals, respectively, β is the transmission rate, γ is the recovery rate, and ki is the degree of the node i. For (2), it is shown that diffusion increases the epidemic threshold of the SIS dynamics in arbitrary heterogeneous networks.In this paper, we incorporate demographics into (2), and the model is as follows:(3)dρS,idt=bi−βρS,iρI,i+γρI,i−DSρS,i+DS∑j=1NAjikjρS,j−μρS,i,dρI,idt=βρS,iρI,i−γρI,i−DIρI,i+DI∑j=1NAjikjρI,j−μρI,i,where bi is the birth rate at node i and μ is the natural death rate of individuals (whose meaning is different from that in (2)).It is easy to see that (3) is well posed, that is, for any given nonnegative initial condition, it has a unique and global solution through it. Moreover, the solution is also nonnegative. Let b¯=∑i=1Nbi and H=∑i=1NρS,i+ρI,i. Adding up all equations of (3) yields(4)dHtdt=b¯−μHt,which implies that limsupt⟶∞Ht≤b¯/μ. Thus, the feasible region of (3) can be chosen as(5)Γ=ρS,1,…,ρS,N,ρI,1,…,ρI,N∈R+2N|∑i=1NρS,i+ρI,i≤b¯μ.It can be verified thatΓ is attractive and positively invariant with respect to (3) (see [16]). Let Γ∘ denote the interior of Γ and ∂Γ denote the boundary of Γ.


## 2.2. Existence of Disease-Free Equilibria

We start the analysis of (3) with the existence of disease-free equilibria. At a disease-free equilibrium of (3), we have the following linear system [17]:(6)bi−DSρS,i+DS∑j=1NAjikjρS,j−μρS,i=0,i=1,2,…,N,which can be written in the matrix form(7)L¯S=B¯,where(8)L¯=DS+μ−DSA21k2⋯−DSAN1kN−DSA12k1DS+μ⋯−DSAN2kN⋮⋮⋱⋮−DSA1Nk1−DSA2Nk2⋯DS+μ,S=ρS,1,ρS,2,…,ρS,NT,B¯=b1,b2,…,bNT.AsL¯ is a column diagonal dominant matrix with positive diagonal elements and negative off-diagonal elements, L¯ is nonsingular and any element of L¯−1 is positive [18]. Thus, (6) has a unique positive solution, which has the following result.Proposition 1.
System (3) always has a unique disease-free equilibrium, denoted by E0=ρS,10,…,ρS,N0,0,…,0.Since the adjacency matrixAji is randomly given, the disease-free equilibrium E0 given in Proposition 1 does not have an explicit expression. Thus, we calculate it in specific network: globally coupled network. In a globally coupled network,(9)A=01⋯110⋯1⋮⋮⋱⋮11⋯0.Then,(10)L¯=DS+μ−DSN−1⋯−DSN−1−DSN−1DS+μ⋯−DSN−1⋮⋮⋱⋮−DSN−1−DSN−1⋯DS+μ.By complex calculations, we get(11)L¯−1=N−1μ+DSNμDS+μ2N−1DSNμDS+μ2N−1⋯DSNμDS+μ2N−1DSNμDS+μ2N−1N−1μ+DSNμDS+μ2N−1⋯DSNμDS+μ2N−1⋮⋮⋱⋮DSNμDS+μ2N−1DSNμDS+μ2N−1⋯N−1μ+DSNμDS+μ2N−1,and hence(12)ρS,i0=N−1μ+DSbi+DS−bi+∑j=1NbjNμDS+μ2N−1,fori=1,2,…,N.In particular, ifb1=b2=⋯=bN=b, then ρS,i0=μ/b and the disease-free equilibrium is E0=μ/b,…,μ/b︸N,0,…,0.


## 2.3. Stability of the Disease-Free Equilibrium

Before considering the stability of the disease-free equilibrium of (3), we introduce the basic reproduction number.Linearizing (3) at the disease-free equilibrium E0, we can get(13)dItdt=MIt,where It=ρI,1t,ρI,2,…,ρI,NtT and(14)M=βρS,10−γ−μ−DIDIA21k2⋯DIAN1kNDIA12k1βρS,20−γ−μ−DI⋯DIAN2kN⋮⋮⋱⋮DIA1Nk1DIA2Nk2⋯βρS,N0−γ−μ−DI.Then,sM is a simple eigenvalue of M with a positive eigenvector. Let FiI denote the rate of appearance of new infections in node i, Vi+I be the transfer rate of infected individuals into node i by all other means, and Vi−I be the transfer rate of infected individuals out of node i. Then, we can choose(15)F=βρS,100⋯00βρS,20⋯0⋮⋮⋱⋮00⋯βρS,N0,(16)V=V1−I−V1+IV2−I−V2+I⋮VN−I−VN+I=μ+DI+γ−DIA21k2⋯−DIAN1kN−DIA12k1μ+DI+γ⋯−DIAN2kN⋮⋮⋱⋮−DIA1Nk1−DIA2Nk2⋯μ+DI+γ.We callFV−1 the next generation matrix for (3), and the basic reproduction number R0 is defined as(17)R0=ρFV−1,where ρFV−1 denotes the spectral radius of the matrix FV−1 (see [19] for details). Note that R0 is simply the product of the infection rate and the mean duration of the infection [15]. In a globally coupled network, the basic reproduction number is R0=bβ/μγ+μ. Refer to [19] for Theorem 1.Theorem 1.
The following statements on (3) are true.(i)
R0>1 if and only if sM>0 and R0<1 if and only if sM<0.(ii)
IfR0<1, then the disease-free equilibrium E0 is locally asymptotically stable while it is unstable if R0>1.
Now, we establish the global stability of the disease-free equilibrium.Theorem 2.
AssumeR0<1. Suppose that A is irreducible. Then, the disease-free equilibrium E0 of (3) is globally asymptotically stable in R+2N.Proof.
AsΓ is an attracting and positively invariant set, we only need to show the global stability of E0 in Γ.
LetF and V be given by (15) and (16), respectively. Since all off-diagonal entries of V are nonpositive and the sum of the entries in each column of V is positive, V is a nonsingular M-matrix. Also, V−1 is irreducible. By the Perron–Frobenius theorem ([18], p. 27), the nonnegative irreducible matrix V−1F has a positive left eigenvector w1,w2,…,wN corresponding to the eigenvalue ρV−1F. Since F is a diagonal matrix, ρV−1F=ρFV−1. Consequently, we have(18)w1,w2,…,wNV−1F=R0w1,w2,…,wN,or(19)1R0w1,w2,…,wN=w1,w2,…,wNF−1V.
Fori=1,…,N, denote ci=wi/βρS,i0>0 and I=ρI,1,ρI,2,…,ρI,NT. Consider the Lyapunov function(20)L=∑i=1NciρI,i,and the derivative of L along the solutions of (3) is(21)dLdt|2.2=∑i=1NciβρS,iρI,i−γρI,i−DIρI,i+DI∑j=1NAjikjρI,j−μρI,i≤∑i=1NciβρS,i0−γ−μ−DIρI,i+DI∑j≠iNAjikjρI,j=w1βρS,10,w2βρS,20,…,wNβρS,N0F−VI=w1,w2,…,wN1−F−1VI=w1,w2,…,wN1−1R0I≤0.
Sinceci>0 for all i, dL/dt|2.2=0 implies that either ρS,i=ρS,i0 or ρI,i=0 for any 1≤i≤N. When ρS,i=ρS,i0, it follows from the first equation of (3) that(22)0=dρS,i0dt=bi−βρS,i0ρI,i+γρI,i−DSρS,i0+DS∑j=1NAjikjρS,j0−μρS,i0.
Comparing this with (6) gives ρI,i=0. Thus, we have shown that dL/dt|2.2=0 which implies that ρI,i=0 for all i. By LaSalle’s invariance principle, we see that E0 is globally asymptotically stable in Γ and hence in R+2N.
The following result tells us that whenR0>1, the disease is persistent.Proposition 2 (see [17]).
Suppose thatA is irreducible. If R0>1, then system (3) is uniformly persistent and there exists an endemic equilibrium E∗ in Γ∘.The proof of Proposition2 is similar to the proof of Proposition 3.3 in [20]. Using uniform persistence of (3) and uniform boundedness of solutions in Γ∘ we can obtain the existence of an equilibrium of (3) in Γ∘(see Theorem 2.8.6 in [21]). In Proposition 2, the assumption that the adjacency matrix A=Aij is irreducible is necessary. If A=0, then system (3) has an asymptotically stable boundary equilibrium when R0>1, and thus the system is not persistent [22].


## 2.4. The Global Stability of the Endemic Equilibrium

In this section, we provide sufficient conditions on the global stability of the endemic equilibrium. The method of proof is the graph-theoretical approach developed in [23–25].Theorem 3.
Assume thatR0>1. Let E∗=ρS,1∗,…,ρS,N∗,ρI,1∗,…,ρI,N∗ be an endemic equilibrium of (3). If A is irreducible and there exists λ>0 such that DS∑j=1NAji/kjρS,j∗=λDI∑j=1NAji/kjρI,j∗ for all 1≤i,j≤N, then the endemic equilibrium E∗ is globally asymptotically stable in Γ∘, and hence (3) only has a unique endemic equilibrium.Proof.
For each1≤i≤N, set(23)ViρS,i,ρI,i=ρS,i−ρS,i∗−ρS,i∗lnρS,iρS,i∗+ρI,i−ρI,i∗−ρI,i∗lnρI,iρI,i∗.
Noting(24)DSρS,i∗=bi−βρS,i∗ρI,i∗+γρI,i∗+DS∑j=1NAjikjρS,j∗−μρS,i∗,DI+μρI,i∗=βρS,i∗ρI,i∗+DI∑j=1NAjikjρI,j∗−γρI,i∗,we calculate the derivative of Vi along solutions of (3) as(25)dVidt|2.2=bi−DSρS,i+DS∑j=1NAjikjρS,j−μρS,i−biρS,i∗ρS,i+βρI,iρS,i∗−γρI,iρS,i∗ρS,i+DSρS,i∗−DS∑j=1NAjikjρS,jρS,i∗ρS,i+μρS,i∗−DIρI,i+DI∑j=1NAjikjρI,j−μρI,i−βρS,iρI,i∗+γρI,i∗+DIρI,i∗+μρI,i∗−DI∑j=1NAjikjρI,jρI,i∗ρI,i=bi1−ρS,iρS,i∗+lnρS,iρS,i∗+1−ρS,i∗ρS,i+lnρS,i∗ρS,i+DS∑j=1NAjikjρS,j∗1−ρS,i∗ρS,jρS,iρS,j∗+lnρS,i∗ρS,jρS,iρS,j∗+DS∑j=1NAjikjρS,j∗ρS,jρS,j∗+lnρS,j∗ρS,j−ρS,iρS,i∗−lnρS,i∗ρS,i+DI∑j=1NAjikjρI,j∗1−ρI,i∗ρI,jρI,iρI,j∗+lnρI,i∗ρI,jρI,iρI,j∗+DI∑j=1NAjikjρI,j∗ρI,jρI,j∗+lnρI,j∗ρI,j−ρI,iρI,i∗−lnρI,i∗ρI,i+γρI,i∗1−ρS,i∗ρI,iρS,iρI,i∗+lnρS,i∗ρI,iρS,iρI,i∗+γρI,i∗lnρI,i∗ρI,i−lnρS,i∗ρS,i+γρI,i1−ρI,i∗ρS,iρI,iρS,i∗+lnρI,i∗ρS,iρI,iρS,i∗+γρI,ilnρI,iρI,i∗−lnρS,iρS,i∗≤DS∑j=1NAjikjρS,j∗ρS,jρS,j∗+lnρS,j∗ρS,j−ρS,iρS,i∗−lnρS,i∗ρS,i+DI∑j=1NAjikjρI,j∗ρI,jρI,j∗+lnρI,j∗ρI,j−ρI,iρI,i∗−lnρI,i∗ρI,i+γρI,i∗lnρI,i∗ρI,i−lnρS,i∗ρS,i+γρI,ilnρI,iρI,i∗−lnρS,iρS,i∗=DS∑j=1NAjikjρS,j∗ρS,jρS,j∗+lnρS,j∗ρS,j−ρS,iρS,i∗−lnρS,i∗ρS,i+DI∑j=1NAjikjρI,j∗ρI,jρI,j∗+lnρI,j∗ρI,j−ρI,iρI,i∗−lnρI,i∗ρI,i=DI∑j=1NAjikjρI,j∗λρS,jρS,j∗+λlnρS,j∗ρS,j+ρI,jρI,j∗+lnρI,j∗ρI,j−DI∑j=1NAjikjρI,j∗λρS,iρS,i∗+λlnρS,i∗ρS,i+ρI,iρI,i∗+lnρI,i∗ρI,i=DI∑j=1NAjikjρI,j∗GjρS,j,ρI,j−GiρS,i,ρI,i,where(26)GiρS,i,ρI,i=λρS,iρS,i∗+λlnρS,i∗ρS,i+ρI,iρI,i∗+lnρI,i∗ρI,i.
Here we have used the fact that1−x+lnx≤0 for x>0, and the equality holds if and only if x=1.
Consider a weight matrixW=wij with entries wij=DIAji/kjρI,j∗ and denote the corresponding weighted digraph as G,W. Let li=∑T∈TiwT≥0 be the same as that given in (A.1) in the Appendix associated with G,W. Then, by (A.2), the following identity holds:(27)∑i=1Nli∑j=1NDIAjikjρI,j∗GjρS,j,ρI,j−GiρS,i,ρI,i=0.
Set(28)VρS,1,ρI,1,…,ρS,N,ρI,N=∑i=1NliViρS,i,ρI,i.
Then, we obtain(29)dVdt|2.2=∑i=1NlidVidt|2.2≤∑i=1Nli∑j=1NDIAjikjρI,j∗GjρS,j,ρI,j−GiρS,i,ρI,i=0,for all ρS,1,ρI,1,…,ρS,N,ρI,N∈Γ∘. Since A is irreducible, we know that li>0 for all i (see the Appendix), and thus dV/dt|2.2=0 implies that ρS,i=ρS,i∗ for all i. From the first equation of (3), we obtain(30)0=bi−βρS,i∗ρI,i+γρI,i−DSρS,i∗+DS∑j=1NAjikjρS,j∗−μρS,i∗,for i=1,…,N, which implies that ρI,i=ρI,i∗. This means that the largest invariant set in dV/dt|2.2=0 is the singleton E∗. By LaSalle’s invariance principle, E∗ is globally asymptotically stable in Γ∘.


## 3. Numerical Simulations

The results in the previous section presented the stability of the disease-free and endemic equilibria. Through numerical simulations of infection dynamics, it is shown that the connectivity between patches has a large effect on disease dynamics. In this section, an Erdős–Rényi network with 100 nodes is established to study the influence of diffusion, and the connection probability is 0.4. All parameters are positive. It should be noted that we simulated the numerical solution of the model equation, not Monte Carlo (MC) simulations of the metapopulation dynamics [14].Figure1 shows the average number of susceptible and infected individuals with different initial values. The parameters are used as N=100, b=1, μ=0.00312, β=0.125, DS=0.2, DI=0.2, and γ=0.6. We select three different initial conditions for simulation, and the three initial values were C1: the ratio of susceptible number to infected number is 1 : 9; C2: the ratio of susceptible number to infected number is 1 : 1; and C3: the ratio of susceptible number to infected number is 9 : 1. Susceptible and infected individuals obey Poisson distribution in nodes. It can be seen in Figure 1 that under different initial conditions, both susceptible individuals and infected individuals converge to a unique globally asymptotically stable endemic equilibrium.Figure 1
The average number of susceptible and infected individuals with different initial values.
(a)(b)In Figure2, we consider the impact of diffusion on the number of infected individuals at the endemic equilibrium. For this purpose, we take N=100, b=6, μ=0.00512, β=0.325, and γ=0.6. In Figure 2(a), we fix DI=1. We observe that ρI,i∗ increases as DS increases and for the same DS, when i is larger, ρI,i∗ is also larger. In Figure 2(b), we take DS=2. It is observed that the same phenomena happen with respect to DI and i with fixed DI.Figure 2
The impact of diffusion on the number of infected individuals at the endemic equilibria, and each line corresponds to the nodes with degreesk of 1, 10, 26, and 66, respectively (from bottom to top). (a) Relationship between DS and ρI,i∗. (b) Relationship between DI and ρI,i∗.
(a)(b)The relationship betweenβ, DI, and R0 is simulated in Figure 3(a), and the relationship between μ, DI, and R0 is simulated in Figure 3(b). It can be seen from Figure 3 that the influence of DI on R0 is greater than that of β and μ, so the influence of the diffusion rates for the infected individuals on disease outbreak cannot be ignored.Figure 3
The impact of diffusion and other parameters onR0. (a) Relationship between β, DI, and R0. (b) Relationship between μ, DI, and R0.
(a)(b)


## 4. Conclusion

In this paper, we formulated and analyzed an SIS model with demographics and diffusion on metapopulation networks. In this model, we solved the stability of the disease-free equilibrium and deduced the basic reproduction number. We concluded that whenR0<1, the disease-free equilibrium is globally asymptotically stable; otherwise, it is unstable. Next, we proved the existence and stability of endemic equilibria. The endemic equilibrium is globally asymptotically stable when R0>1. Finally, we gave some numerical simulations to illustrate the main results and demonstrate the effects of connectivity and diffusion on the number of infected individuals at the endemic equilibrium. As the connectivity or diffusion is increasing, the number of infected individuals is also increasing. This means that when an epidemic comes, the more people there are and the greater the degree is, the more likely the disease will break out. Therefore, we should pay attention to the prevention and control of diseases and make plans in advance.The research of this paper can be extended to the following aspects: in this paper, the Erdős–Rényi network is used to simulate, and the models based on other networks can be analyzed; secondly, the infectious disease model with heterogeneous diffusion rate can also be analyzed. It is necessary to study the mathematical model under a specific network in the future, so as to make the model more realistic.



---

*Source: 8884236-2021-09-20.xml*