Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 830237, 11 pages

http://dx.doi.org/10.1155/2013/830237

## Disease Control of Delay SEIR Model with Nonlinear Incidence Rate and Vertical Transmission

^{1}School of Mathematical Sciences, Dalian University of Technology, Dalian 116024, China^{2}Department of Mathematics, Tonghua Normal University, Tonghua 136000, China^{3}School of Innovation Experiment, Dalian University of Technology, Dalian 116024, China

Received 19 July 2013; Revised 30 September 2013; Accepted 3 October 2013

Academic Editor: Seiya Imoto

Copyright © 2013 Yan Cheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The aim of this paper is to develop two delayed SEIR epidemic models with nonlinear incidence rate, continuous treatment, and impulsive vaccination for a class of epidemic with latent period and vertical transition. For continuous treatment, we obtain a basic reproductive number and prove the global stability by using the Lyapunov functional method. We obtain two thresholds and for impulsive vaccination and prove that if then the disease-free periodic solution is globally attractive and if then the disease is permanent by using the comparison theorem of impulsive differential equation. Numerical simulations indicate that pulse vaccination strategy or a longer latent period will make the population size infected by a disease decrease.

#### 1. Introduction

Mathematical models describing the population dynamics of infectious diseases have been playing an important role in understanding epidemiological patterns and disease control. Researchers have studied the epidemic models by ordinary differential equations [1–3] and the references cited therein. A customarily epidemic model is susceptible, infectious, and recovered model (SIR for short) [4–7]. But in real life, many diseases have a period of incubation time inside the hosts before the hosts become infectious; if we include incubation period of the hosts, the model is described as SEIR model. As tuberculosis (TB), measles and so on, a susceptible individual becomes exposed (infected but not infective) by adequate contact with an infectious individual. SEIR infections disease model has been studied by many authors for its important biological meaning [8–13]. In [13], the authors considered the following delayed SEIR epidemic model: where , , , and represent the number of individuals who are susceptible, exposed, infected, and removed, respectively. The parameters , , , and are positive constants, and here is the constant recruitment rate into the population, is the contact rate, is the birth and death rate, is the removal rate. represents a time delay describing the latent period of the disease and the term represents the individuals surviving in the latent period and becoming infective at time . The sufficient conditions are obtained for the global asymptotic stability of the endemic equilibrium.

In the study of epidemic model, the spread of an infectious disease is a crucial issue, which depends on both the population behavior and the infectivity of the disease. These two aspects are captured in the incidence rate of a disease. In many epidemiological models, the incidence rate is described as mass action incidence with bilinear interactions given by , where is the probability of transmission per contact and and represent the susceptible and infected populations, respectively. This contact law is more appropriate for a few of infected individuals; when the size of infected individuals is increasing, the underlying assumption of homogeneous mixing may not be valid. In fact, with the increase of infected populations, the susceptible individual will take measure to prevent unbounded contact rates. In [14], Anderson and May proposed a saturated incidence rate of the form in which measures the infection force of the disease and measures the inhibition effect from the behavioral change of the susceptible individuals. The same as the nonlinear incidence rates of the form were investigated by Liu et al. [15, 16].

In real life, some diseases may be transferred through horizontal transmission and vertical transmission (disease is the passing of an infection to offspring of infected parents). The offspring of infected parents may already be infected with the disease at birth, so many infections in nature transmit through both horizontal and vertical modes, such as tuberculosis (TB), rubella, hepatitis B, and AIDS [17–21].

Vaccination and treatment are important strategy for the elimination of infectious diseases. Recently, pulse vaccination has been confirmed as an effective method to prevent the spread of the disease [22–24]. Theoretical results show that the pulse vaccination strategy can be distinguished from the conventional strategies in leading to disease eradication at relatively low values of vaccination [25]. The study of vaccination, treatment, and associated behavioral changes related to disease transmission has been the subject of intense theoretical analysis.

The literature on SEIR model with nonlinear incidence, constant infectious period, impulsive vaccination, dealing with the analysis of disease that is vertically and horizontally transmitted is not extensive [17, 19]. But, in fact, under the situation of disease with vertical transmission, the continuous treatment should be considered for the infected, and impulsive vaccination to the susceptible, newborns of the susceptible, exposed and the removed, and newborns of infected which not be vertical infected.

Motivated by the literature above, we introduce delay epidemic models with nonlinear incidence rates of the form , and we also considered the constant latency period and vertically and horizontally in (2). The purpose of this paper is to study the nonlinear dynamics of system, and we consider two different strategies to the model which are constant treatment and pulse vaccination to the newborns and susceptible,

The basic assumptions are as follows.(i)The total population size at time (day) is denoted by . For , this shows that the total population has a constant size. Without loss of generality, we assume in this paper . The newborns of , and are susceptible individuals, and the newborns of who are not vertically infected are also susceptible individuals.(ii)The positive constant (per day) denotes the death rate and birth rate of susceptible, exposed, and recovered individuals. The positive constant (per day) denotes the death rate and birth rate of infective individuals. The positive constant (per day) is the natural recovery rate of infective individuals. The positive constant () (per day) is the vertical transmission rate, and note , (per day), and then . Fraction of all newborns with mothers in the susceptible, exposed, and recovered classes are vaccinated and appeared in the recovered class, while the remaining fraction, , appears in the susceptible class; suppose . is the proportion of those cured successfully.(iii)The incidence rate is described by a nonlinear function where (per day) is a positive constant describing the infection rate. is the length of the latent period and the term reflects the fact that an individual is surviving in the latent period and becoming infective at time .

The remaining part of this paper is organized as follows. In Section 2, we investigate the global stability of the endemic equilibrium of (2) by using Rouches theorem and Lyapunov-LaSalle type theorem. The global asymptotic stability of disease-free periodic solution and the conditions for the permanence of the disease by comparison techniques are described in Section 3. Numerical simulations are presented in Section 4. In Section 5, we conclude this paper with some remarks.

#### 2. Continuous Treatment Strategy of the SEIR Model

In this section, we consider a continuous treatment of SEIR model with constant latent period and nonlinear incidence rate. By using , notice that first and third equations of system (2) do not contain the variables and ; therefore, system (2) is equivalent to the following 2-dimensional system:

##### 2.1. Disease-Free Equilibrium and Its Stability

First, we define

Theorem 1. *If , the disease-free equilibrium of system (3) is locally asymptotically stable for all ; if , the disease-free equilibrium is unstable.*

*Proof. *Steady states of system satisfy the following system of equations:

Obviously, is the disease-free equilibrium of (3). In order to analyze the behavior of the system (3) near , we linearize the system about the equilibrium point; let , is one of the eigenvalues of the linearization of system (6) near the steady state , and the other eigenvalue is determined by equation:
Let
if , it is easy to show that, for real,
hence, has a positive real root. Therefore, if , the disease-free equilibrium is unstable.

If , we prove that the disease-free equilibrium is locally stable. Otherwise, . We note that
a contradiction. Hence, the disease-free equilibrium is locally asymptotically stable if .

Theorem 2. *If , the disease-free equilibrium of system (3) is globally asymptotically stable for all .*

To proof the global stability of the disease-free equilibrium , we choose Lyapunov function and it is easy to prove , ; it follows that , .

##### 2.2. Endemic Equilibrium and Its Stability

If , then system (3) has a unique positive equilibrium , where

Theorem 3. *If , conditions (17) and (22) are satisfied, then for the endemic equilibrium of system (3) is locally asymptotically stable.*

*Proof. *Let , ; the linearized system is obtained
From the linearized system we obtain the characteristic equation
where
For the characteristic equation becomes
and we can see that both roots are negative and real if and only if
Now for , if is a root of (16), we have
Separating the real and imaginary parts, we have
Adding both equations and regrouping by powers of , we obtain the following fourth degree polynomial
from which we have
It follows that if
are satisfied, (20) does not have positive solutions, and the characteristic equation (14) does not have purely imaginary roots. Inequalities in (17) and (22) guarantee that all roots of (14) have no positive roots. According to Rouche’s theorem, Theorem 3 is proved.

Subsequently we discuss the sufficient conditions under which the endemic equilibrium is globally asymptotically stable for the system (3). For , hence, the dynamics of system (3) in the first octant of is equivalent to that of the following system:

The initial conditions for system (23) take the form where , the space of continuous functions mapping the interval into , where .

For continuity of the initial conditions, we require

It is well known by the fundamental theory of functional differential equations [26], the system (23) has a unique solution satisfying the initial conditions. It is easy to show that all solutions of system (23) with initial conditions are defined on and remain positive for all .

Lemma 4 (see [25]). *Let the initial condition be , and , for all . Then .*

Theorem 5. *Let the initial condition be , and , for all . Further suppose ; then for any infectious period satisfying
**
where , the endemic equilibrium is globally asymptotically stable.*

*Proof. *Let , , ; the linearized system is:

The trivial solution of system (27) is globally asymptotically stable and is equivalent to the fact that the endemic equilibrium of system (23) is globally asymptotically stable. We will employ Lyapunov functional technique to prove it.

Now let us introduce the following functions:
where is an arbitrary real constant. Choosing , the derivative of is
and applying Cauchy-Chwartz inequality to all product terms, we obtain the following expression:
We choose Lyapunov function to be the form , and we get
Substituting this in the inequality for , we get
The right-hand expression of the above inequality is always negative provided that (26) holds. A direct application of the Lyapunov-LaSalle type theorem shows that , , . The proof is complete.

#### 3. Continuous Treatment and Pulse Vaccination Strategies

When continuous treatment and pulse vaccination strategies are included in the SEIR epidemic model with the nonlinear infectious force and vertical transmission, it can be written as follows: and , , , and are the number of susceptible, exposed, infectious, and recovered at time , respectively. is the proportion of those vaccinated successfully at , which is called pulse vaccination rate.

We also consider the following reduced systems: Let be the following subset of , . From biological considerations, we discuss system (34) in the closed set . It can be verified that is positively invariant with respect to system (34).

We first state two important lemmas which are useful in our following discussions.

Lemma 6 (see [12]). *Consider the following impulsive differential equation
**
where , , . Then the above system has a unique positive periodic solution given by
**
which is globally asymptotically stable, where
*

Lemma 7 (see [27, 28]). *Consider the following equation:
**
where , for . One has *(1)*if , then ;*(2)*if , then .*

##### 3.1. Global Stability of the Disease-Free Periodic Solution

Now we will prove the disease-free periodic solution is global attractively. We first demonstrate the existence of the disease-free periodic solution, in which infectious individuals are entirely absent from the population permanently, that is, for all . Under this condition, the growth of susceptible individuals must satisfy By Lemma 6, we obtain the periodic solution of system (39), and this solution is globally asymptotically stable. Hence, the system (34) has a disease-free periodic solution .

Theorem 8. *Let be any solution of (34); then the disease-free periodic solution is globally asymptotically stable provided that
*

*Proof. *Since , we can choose sufficiently small such that
From the first equation of system (34), we have , and then we consider the following comparison system with pulse:
In view of Lemma 6, we obtain and
There exists an integer , such that
That is,
Furthermore, from the second equation, we have
and then we consider the following comparison equation:
then, from (41) we have . In view of Lemma 7, we have , . So there must exist an integer , such that for all .

When , from the first equation of system (34), we have
Consider the following comparison impulsive differential equation for all
By Lemma 6, we have the unique periodic solution of system (50) given by
By the comparison theorem, there exists an integer such that
Because and are sufficiently small, it follows from (46) and (52) that . Therefore, the disease-free solution of system (34) is globally attractive. The proof is completed.

Denote

Corollary 9. *(i) If , then the infection-free periodic solution is globally attractive.**(ii) If , then the infection-free periodic solution is globally attractive provided that .*

Corollary 10. *(i) If , then the infection-free periodic solution is globally attractive.**(ii) If , then the disease will be endemic and system (34) is permanent provided that .*

Theorem 8 determines the global attractivity of (34) in for the case . Its epidemiological implication is that the infectious population vanishes in time so the disease dies out. Corollaries 9 and 10 imply that the disease will disappear if the vaccination rate or the length of latent period of the disease is large enough.

##### 3.2. Persistent

In this section we say the disease is endemic if the infectious population persists above a certain positive level for sufficiently large time. The endemicity of the disease can be well captured and studied through the notion of uniform persistence.

*Definition 11. *System (34) is said to be uniformly persistent if there exist positive constants , , (both are independent of the initial values), such that every solution with positive initial conditions of system (34) satisfies

Denote

Theorem 12. *If , then there is a positive constant such that each positive solution of system (34) satisfies for all t sufficiently large.*

*Proof. *Let be any solution with initial values of system (34), and then it is obvious that , for all . We are left to prove there exist positive constants , and ( is sufficiently large) such that , for all .

Firstly, from the first equation of system (34), we have
Consider the following comparison equations:

By Lemma 6 and the comparison theorem [29], we know that for any sufficiently small , there exists a ( is sufficiently large) such that

Now, we will prove that there exist and a sufficiently large such that for all . Since the proof is rather long, it will be convenient to divide it into two steps.*Step 1.* Since , there exist , sufficiently small such that
where .

We claim that for any , it is impossible that for all . Suppose that the claim is not valid. There exists a such that for all , and then follows from the first equation of system (34) that for ,
Consider the comparison impulsive system for ,
According to Lemma 6, there exists such that
for all . The second equation of system (34) can be translated into the following form:
Define a function such that
then the derivative of along the solution of system (34) is
From (55), we obtain , , which implies that , . This is contrary to . Hence, there exists a such that .*Step 2.* According to Step 1, for any positive solution of system (34), we are left to consider two cases. First, if for all , then our aim is obtained. Second oscillates about for all large . In this case, setting , there are two possible cases for .

Define
We hope to show that for all large . The conclusion is evident in the first case. For the second case, let and satisfy , and , for . Therefore, it is certain that there exists a () such that
In this case, we will discuss three possible cases in terms of the sizes of , , and .*Case 1.* If , then for .*Case 2.* If , then from the second equation of system (34), we can deduce for and , and it is obvious that for .*Case 3.* If , we will consider the following two cases, respectively.*Subcase 3.1.* For , it is easy to obtain .*Subcase 3.2.* For , it is easy to obtain . Then, proceeding exactly as the proof for the above claim, we see that for . Since this kind of interval is chosen in an arbitrary way (we only need to be large), we conclude that for all large in the second case. In view of our above discussions, the choices of are independent of the positive solution, and we have proved that any positive solution of (34) satisfies for all large . The proof is completed.

Set

From Theorems 12, we also easily obtain the following results.

Corollary 13. *(i) If , then the disease will be endemic and system (34) is permanent provided that .**(ii) If , then the disease will be endemic and system (34) is permanent provided that .*

#### 4. Numerical Simulations

In this section, we present some numerical simulations to demonstrate our theoretical results established in this paper.

*Example 1*. Letting , , , , , , , , , , , we consider the pulse vaccination strategy (see Figure 1), the result shows that the disease fades away when the proportion of those vaccinated successfully ; but the disease will exist everlasting when the proportion of those vaccinated successfully . So, this verifies the results in Corollaries 9 and 13, for the epidemic disease with vertical transition periodical vaccination is an effective method to prevent the disease. Also, it can be seen that with the increase of which typically causes oscillation bigger of the susceptible.

*Example 2*. We use the same parameters as in Example 1 except choosing and , respectively. Compare with is obviously that the longer of the latent period the lower of the infective number (see Figure 2), and this shows disease with long latent period disadvantage to the spread of the disease.

*Example 3*. For the same parameters as in Example 1 and take and , we obtain the number of infected individual (see Figure 3), when . With the increase of and the number of infected individual is decreasing, so these results indicate that it will be helpful to control the disease with vertical transition for bigger and .

#### 5. Conclusion

In this paper, for a class of epidemic disease with latent period and vertical transition, we present two delayed SEIR epidemic models with nonlinear incidence rate based on the spread characters of the disease (such as tuberculosis). Our model is more approach to the realistic problem which is different from [17, 19]. Moreover, the methods in our model are different from the existing results because more factors are considered. When only considering constant treatment, we obtain basic reproductive number and prove the global stability by using the Lyapunov functional method. For the SEIR model with pulse vaccination we also get the theoretical result, if , the disease-free periodic solution is globally attractive; and if , the disease is permanent by using the comparison theorem of impulsive differential equation. By some simulation experiments, it clearly shows that the larger of the proportion of those vaccinated successfully the lower of the infective individuals, and the longer of the latent period the lower of the infective individuals. So these results demonstrate that it will be helpful to control the disease with vertical transition for bigger and .

#### Acknowledgment

The work is supported by the National Natural Science Foundation of China (no. 1124319).

#### References

- M. E. Alexander and S. M. Moghadas, “Periodicity in an epidemic model with a generalized non-linear incidence,”
*Mathematical Biosciences*, vol. 189, no. 1, pp. 75–96, 2004. View at Publisher · View at Google Scholar · View at Scopus - R. M. Anderson and R. M. May, Eds.,
*Infectious Diseases of Humans: Dynamics and Control*, Oxford University Press, New York, NY, USA, 1991. - Y. Jin, W. Wang, and S. Xiao, “An SIRS model with a nonlinear incidence rate,”
*Chaos, Solitons and Fractals*, vol. 34, no. 5, pp. 1482–1497, 2007. View at Publisher · View at Google Scholar · View at Scopus - A. Korobeinikov, “Lyapunov functions and global stability for SIR and SIRS epidemiological models with non-linear transmission,”
*Bulletin of Mathematical Biology*, vol. 68, no. 3, pp. 615–626, 2006. View at Publisher · View at Google Scholar · View at Scopus - L. Nie, Z. Teng, and A. Torres, “Dynamic analysis of an SIR epidemic model with state dependent pulse vaccination,”
*Nonlinear Analysis: Real World Applications*, vol. 13, no. 4, pp. 1621–1629, 2012. View at Publisher · View at Google Scholar · View at Scopus - X. Meng and L. S. Chen, “The dynamics of a new SIR epidemic model concerning pulse vaccination strategy,”
*Applied Mathematics and Computation*, vol. 197, no. 2, pp. 582–597, 2008. View at Publisher · View at Google Scholar · View at Scopus - D. Xiao and S. Ruan, “Global analysis of an epidemic model with nonmonotone incidence rate,”
*Mathematical Biosciences*, vol. 208, no. 2, pp. 419–429, 2007. View at Publisher · View at Google Scholar · View at Scopus - A. d'Onofrio, “Stability properties of pulse vaccination strategy in SEIR epidemic model,”
*Mathematical Biosciences*, vol. 179, no. 1, pp. 57–72, 2002. View at Publisher · View at Google Scholar · View at Scopus - M. A. Safi and S. M. Garba, “Global stability analysis of SEIR model with holling type II incidence function,”
*Computational and Mathematical Methods in Medicine*, vol. 2012, Article ID 826052, 8 pages, 2012. View at Publisher · View at Google Scholar - S. Gao, Z. Teng, and D. Xie, “The effects of pulse vaccination on SEIR model with two time delays,”
*Applied Mathematics and Computation*, vol. 201, no. 1-2, pp. 282–292, 2008. View at Publisher · View at Google Scholar · View at Scopus - G. Li and Z. Jin, “Global stability of a SEIR epidemic model with infectious force in latent, infected and immune period,”
*Chaos, Solitons and Fractals*, vol. 25, no. 5, pp. 1177–1184, 2005. View at Publisher · View at Google Scholar · View at Scopus - S. Gao, L. S. Chen, and Z. Teng, “Impulsive vaccination of an SEIRS model with time delay and varying total population size,”
*Bulletin of Mathematical Biology*, vol. 69, no. 2, pp. 731–745, 2007. View at Publisher · View at Google Scholar · View at Scopus - R. Xu and Z. Ma, “Global stability of a delayed SEIRS epidemic model with saturation incidence rate,”
*Nonlinear Dynamics*, vol. 61, no. 1-2, pp. 229–239, 2010. View at Publisher · View at Google Scholar · View at Scopus - R. M. Anderson and R. M. May, “Regulation and stability of host-parasite population interactions. I. Regulatory processes,”
*Journal of Animal Ecology*, vol. 47, no. 1, pp. 219–247, 1978. View at Publisher · View at Google Scholar - W. M. Liu, H. W. Hethcote, and S. A. Levin, “Dynamical behavior of epidemiological models with nonlinear incidence rates,”
*Journal of Mathematical Biology*, vol. 25, no. 4, pp. 359–380, 1987. View at Publisher · View at Google Scholar · View at Scopus - W. M. Liu, S. A. Levin, and Y. Iwasa, “Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models,”
*Journal of Mathematical Biology*, vol. 23, no. 2, pp. 187–204, 1985. View at Google Scholar · View at Scopus - M. Herrera et al., “Modeling the spread of tuberculosis in semiclosed communities,”
*Computational and Mathematical Methods in Medicine*, vol. 2013, Article ID 648291, 19 pages, 2013. View at Publisher · View at Google Scholar - A. d'Onofrio, “On pulse vaccination strategy in the SIR epidemic model with vertical transmission,”
*Applied Mathematics Letters*, vol. 18, no. 7, pp. 729–732, 2005. View at Publisher · View at Google Scholar · View at Scopus - Z. Lu, X. Chi, and L. S. Chen, “The effect of constant and pulse vaccination on SIR epidemic model with horizontal and vertical transmission,”
*Mathematical and Computer Modelling*, vol. 36, no. 9-10, pp. 1039–1057, 2002. View at Google Scholar · View at Scopus - M. Kgosimore and E. M. Lungu, “The effects of vaccination and treatment on the spread of HIV/AIDS,”
*Journal of Biological Systems*, vol. 12, no. 4, pp. 399–417, 2004. View at Publisher · View at Google Scholar · View at Scopus - M. C. Boily and R. M. Anderson, “Sexual contact patterns between men and women and the spread of HIV-1 in urban centres in Africa,”
*IMA Journal of Mathematics Applied in Medicine and Biology*, vol. 8, no. 4, pp. 221–247, 1991. View at Google Scholar · View at Scopus - A. B. Sabin, “Measles, killer of millions in developing countries: strategy for rapid elimination and continuing control,”
*European Journal of Epidemiology*, vol. 7, no. 1, pp. 1–22, 1991. View at Google Scholar · View at Scopus - C. A. de Quadros, J. K. Andrus, J. Olive et al., “Eradication of poliomyelitis: progress in the Americas,”
*Pediatric Infectious Disease Journal*, vol. 10, no. 3, pp. 222–229, 1991. View at Google Scholar · View at Scopus - M. Ramsay, N. Gay, E. Miller et al., “The epidemiology of measles in England and Wales: rationale for the 1994 national vaccination campaign,”
*Communicable Disease Report*, vol. 4, no. 12, pp. R141–R146, 1994. View at Google Scholar · View at Scopus - P. Yongzhen, L. Shuping, L. Changguo, and S. Chen, “The effect of constant and pulse vaccination on an SIR epidemic model with infectious period,”
*Applied Mathematical Modelling*, vol. 35, no. 8, pp. 3866–3878, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. Hale,
*Theory of Functional Differential Equations*, Springer, Heidelberg, Germany, 1977. - Y. Kuang,
*Delay Differential Equation with Application in Population Dynamics*, Academic Press, NewYork, NY, USA, 1993. - Y. N. Xiao and L. S. Chen, “Modeling and analysis of a predator-prey model with disease in the prey,”
*Mathematical Biosciences*, vol. 171, no. 1, pp. 59–82, 2001. View at Publisher · View at Google Scholar · View at Scopus - V. Lakshmikantham, D. Bainov, and P. S. Simeonov,
*Theory of Impulsive Differential Equations*, World Scientific, Singapore, 1989.