Mathematical models describing the population dynamics of infectious diseases have played an important role in better understanding epidemiological patterns and disease controls for a long time. The basis of modeling dynamics and evolution of infectious diseases are summarized by Anderson and May  . It is well known that the spread of a communicable disease involves disease-related factors such as infectious agent, mode of transmission, incubation periods, infectious periods, susceptibility, and resistance.
Continuous time deterministic epidemic models are traditionally elaborated as systems of ordinary differential equations. More realistic models should include some states of these systems, and ideally, a real system should be modeled by delay differential equation. Time delay plays an important part in propagation process of the epidemic, and we can simulate the exposed period of infectious diseases, the infections period of patients and the immunity period of recovery of the disease with time delay. However, compared with studies of the dynamical behaviors of the epidemic models with time delays on complex networks, only a few attentions about early studies have been paid to them on complex networks. In recent years, the dynamics of the SIS, SIR, SEI and SEIR epidemic models have received considerable attention      . In order to study the effects of disease latency or immunity in real life, the delay is incorporated in such models  -  . Liu and Deng et al. discussed epidemic SIS model with discrete time delay which represents the infectious period  , and they obtained the basic reproduction number and discussed the persistence of the disease, but they failed to give the proof in detail. Wang and Wang et al. discussed an epidemic SIR model with discrete time delay which represented the exposed period  . Based on the above, in this paper, we will present a suitable epidemic SEIRS model with time delay which represents the infectious period on complex networks using functional differential equations to investigate the epidemic spreading.
The rest of the paper is organized as follows: Section 2 derives a SEIRS model with time delay mechanism on scale-free networks. Then, Section 3 obtains two equilibriums and basic reproductive number. In Section 4, numerical simulations are performed. Finally, we conclude the paper in Section 5.
2. Model Formulation
One of the most effective interventions to contain the spread of epidemic diseases is the delay mechanism as discussed above. In order to investigate the efficiency of delay mechanism, we consider the new SEIRS model with delay mechanism on complex networks. On the complex networks, each individual is represented by a node of the network and the edges are the connections between individuals along which the infection may spread. Taking the connectivity among different nodes into consideration, let , , and be the relative densities of susceptible, exposed, infected and recovered nodes of degree k at time t respectively. In the course of disease transmission, a susceptible individual will be infected with probability if it connects to a exposed (infected) one. r is the state transition rate from exposed to infectious individual. The rate constant of recovery for infected individuals is denoted by . , for recovered nodes, is the effectively constant of immunity, and there are still recovered nodes, when moment from moment, which will lose the immunity ability to become the susceptible nodes again. Here, we assume that the birth rate equals the death rate, and the rate constant is . Thus, the dynamic mean-field reaction rate equations can be written as
where the probability denotes a link pointing to an exposed individual, and satisfies
The probability denotes a link pointing to an infected individual, and satisfies
Where describes the average degree and describes the connectivity distribution. is the total density of exposed individuals, and the is the total density of infected individuals in the whole network.
For simplicity, through this paper, let and . Thus, the system (1) is equivalent to the following model
With the normalization condition, there is
Obviously, the initial conditions for system satisfy
3. Analysis of the Novel SEIRS Model
In this part, we put forward the analytic solution of system (2) for describing the dynamic behavior.
Theorem 1. Denote
The disease-free equilibrium of system (2) always exists and there exists a positive endemic equilibrium of system (2) when .
Proof. To get the equilibrium solution , we need to make the right side of system (2) equal to zero. Then, the equilibrium should satisfy
Thus, we have
According to the normalization condition , we can get
For , we can find that satisfies Equation (5). Therefore, and is a disease-free equilibrium of system (2). Substituting and of Equation (5) into , we can obtain
so that, Equation (6) can be written
Due to and , the equation has a non-trivial solution if and only if .
So, we have
Let’s define , thus, we can compute the base reproduction number as follows
Namely, there is the unique nontrivial solution if and only if . From Equations ((4) and (5)), we are able to obtain
, , , .
So, we define the equilibrium . Hence, when , system (2) has a unique positive endemic equilibrium . The proof is completed.
Remark. The basic reproductive number is obtained by (7), which depends on different model parameters and fluctuations of the degree distribution. The delay parameter cannot change the basic reproductive number. Clearly, in the eco-limits of infinite size network the number of nodes grows to infinity, i.e.,
, then grows to infinity, so the increase of the basic reproductive number, i.e., , is obvious.
Theorem 2. If , the disease-free equilibrium is globally asymptotically stable. When , the epidemic disease is permanent, which means there exists , such that
where is any solution of system (2), satisfying or .
Proof. For , i.e., the values of , and are fixed, a unique corresponding exists. We can discuss the first three equations of system (2) with at all . System (2) satisfies
Thus, the jacobian matrix of the virus-free equilibrium for the system (2) can be obtained, the jacobian matrix is a matrix, such that
We can put forward the characteristic equation of the disease-free equilibrium from jacobian matrix, such that
where and .
Note that is equivalent to which
, so that
, which means that . Therefore, the
eigenvalue of are all negative if ; otherwise, if and only if , here exists a unique positive eigenvalue of . Using the Perron-Feobenius theorem, this suggests that the maximal real part of all eigenvalues of is positive only if . Thus, we can obtain the results of this theorem by using a theorem presented of Lajmanovich and York  . The proof is completed.
4. Numerical Simulations and Sensitivity Analysis
To support and explain our theoretical analysis results, we present several numerical simulations in this section. The system (2) is considered on complex networks, whose the degree distribution is , where the parameter satisfies , . We suppose , , and .
In Figure 1, the parameters are chosen as , , , , , thus the threshold value . According to , we can see that when , grows to zero, i.e., the infectious individuals will ultimately disappear. If , the equilibrium is globally attractive and the disease eventually disappear.
In Figure 2, the parameters are chosen as , , , , , thus the threshold value . We can see that when , grows to a constant, i.e., the epidemic disease is permanent and the number of infected individuals will converge to a positive constant.
In Figure 3, the parameters are chosen as , , , , , . We can see that the corresponding
Figure 1. With , the prevalence of and versus corresponding for
Figure 2. With , the prevalence of and versus corresponding for
Figure 3. With , the prevalence of and versus corresponding for and different parameter .
Figure 4. With , the prevalence of and versus corresponding for and different parameter .
decreases significantly as the delay parameter increases, i.e., a larger delay parameter can weaken the spreading of disease. The larger delay is, the lower corresponding level is.
In Figure 4, the parameters are chosen as , , , , , . We can see that the corresponding decreases significantly as the delay parameter increases, i.e., a larger delay parameter can reduce the endemic level.
In this paper, a novel SEIRS epidemic model of delay mechanism on complex networks has been presented. Using the field theory, we calculated two epidemiologically relevant quantities, i.e., the basic reproductive number for invasion and the endemic prevalence of epidemic diseases. We have proved that if the basic reproductive number , the disease-free equilibrium is globally asymptotically stable; while , the endemic equilibrium is permanent. Moreover, increasing delay parameters can result in the weakness of the diseases spreading and the decrease of population infected. Numerical simulations show that the endemic equilibrium is globally asymptotically stable when (as shown in Figure 2). It is interesting but challenging to discuss the stability of equilibrium . We leave it for our future work. This study has valuable guiding significance in effectively predicting epidemic spreading.
This work is supported in part by the National Natural Science Foundation of China under Grants 61672112, the Teaching Research Project of Hubei Province under Grant 2014260, and the Teaching Research Project of Yangtze University under Grant JY2015010.
 Wei, J. and Zou, X. (2006) Bifurcation Analysis of a Population Model and the Resulting SIS Epidemic Model with Delay. Journal of Computational and Applied Mathematics, 197, 169-187.
 Huang, G., Takeuchi, Y., Ma, W., et al. (2010) Global Stability for Delay SIR and SEIR Epidemic Models with Nonlinear Incidence Rate. Bulletin of Mathematical Biology, 72, 1192-1207.
 Xiao, Y., Chen, L. and ven den Bosch, F. (2002) Dynamical Behavior for a Stage-Structured SIR Infectious Disease Model. Nonlinear Analysis: Real World Applications, 3, 175-190.
 Gao, S., Chen, L. and Teng, Z. (2007) Impulsive Vaccination of an SEIRS Model with Time Delay and Varying Total Population Size. Bulletin of Mathematical Biology, 69, 731-745.
 Li, T., Liu, X.D., Wu, J., et al. (2016) An Epidemic Spreading Model on Adaptive Scale-Free Networks with Feedback Mechanism. Physica A: Statistical Mechanics and Its Applications, 450, 649-656.
 Enatsu, Y., Nakata, Y. and Muroya, Y. (2012) Global Stability of SIRS Epidemic Models with a Class of Nonlinear Incidence Rates and Distributed Delays. Acta Mathematica Scientia, 32, 851-865.
 Nakata, Y. (2011) On the Global Stability of a Delayed Epidemic Model with Transport-Related Infection. Nonlinear Analysis: Real World Applications, 12, 3028-3034.
 Zhang, T. and Teng, Z. (2008) Global Behavior and Permanence of SIRS Epidemic mOdel with Time Delay. Nonlinear Analysis: Real World Applications, 9, 1409-1424. https://doi.org/10.1016/j.nonrwa.2007.03.010
 Liu, Q., Jiang, D., Shi, N., et al. (2016) Asymptotic Behavior of a Stochastic Delayed SEIR Epidemic Model with Nonlinear Incidence. Physica A: Statistical Mechanics and Its Applications, 462, 870-882.
 Wang, J., Wang, J., Liu, M., et al. (2014) Global Stability Analysis of an SIR Epidemic Model with Demographics and Time Delay on Networks. Physica A: Statistical Mechanics and Its Applications, 410, 268-275.