Classification of different Hepatitis B infected individuals with saturated incidence rate

The role of incidence rate is very important in the study of epidemiological models. In this article, the analysis of an epidemic problem for the transmission dynamic of HBV with saturated incidence rate is presented, which is more generalize than bilinear incidence rate. After formulating the new mathematical model, the threshold quantity reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0$$\end{document}R0 is investigated by using the well known approach i.e. next generation matrix and investigate the possible equilibriums such as disease free and endemic equilibria. Then for the local and global behavior of the proposed problem, the local asymptotic stability analysis as well as global asymptotic stability analysis are proved. To prove the global asymptotic stability at disease free equilibrium, the classic Lyapunov function theory is considered. Similarly to show global asymptotic stability at endemic equilibria, the geometrical approach is used, which is the generalization of Lyapunov theory. Finally, numeric of the proposed problem are carried out to show the feasibility of the obtained results and the role of saturated incidence rate.

and infected mother to her/his child during the time of birth, respectively. However Hepatitis B virus cannot be transmitted through water, food, hugging, kissing and causal contact such as in the work place, school etc. (McMahon 2005). The mode of transmission of Hepatitis B virus and HIV are the same, but Hepatitis B virus is fifty to hundreds time more infectious (Thornley et al. 2008;Ali and Zaman 2016).
Infection of Hepatitis B has different phases, such as acute Hepatitis, chronic and carrier Hepatitis. Acute infection of Hepatitis B denoted the first 6 months subsequently exposed some one to Hepatitis B virus. Herein the resistant organism is frequently capable to flawless the infection from the physique then the body must be recovered in a few months complectly for some one, but for the remaining infection grow and leads to more serious stage called chronic stage or life long sickness. Chronic stage of Hepatitis B raises to the infection, that happens when Hepatitis B virus remain in an individual body and over period, the contagion can develop stern healthiness problem. Individuals with carrier Hepatitis frequently have no antiquity of acute sickness, but it can cause liver damaging which become the cause liver damaging and could also grow to liver tumor Mann and Roberts (2011).
Numerous of authors have premeditated different widespread models, which label the dynamic of communicable diseases (Zou et al. 2010;Thornley et al. 2008;Mann and Roberts 2011). The incidence rate is one of the key concept and theatres an important role in the field of mathematical modeling. Bilinear occurrence rate βSI frequently recycled in various epidemic problems (Fan et al. 2001;Li and Ma 2002;Zou et al. 2010), where β represents the contact rate, S represents susceptible individuals and I represents the infectious individuals. For the first time Capasso and Serio introduced the saturated incidence rate βSI 1+αI , which is the generalized form of bilinear incidence rate (Capasso and Serio 1978). This occurrence rate is additional sensible than the bilinear occurrence rate, especially in the case of sexually transmittible diseases, since it embraces the negotiating alteration and swarming influence of the virulent individuals and inhibits the unboundedness of the interaction rate by indicating appropriate parameters, which was recycled in numerous of epidemic problems (Kar and Jana 2013;Gomes et al. 2005;Liu and Yang 2005).
In this article, we present an epidemic problem for the transmission dynamic of Hepatitis B virus with saturated incidence rate, which is the modified version of Zou et al. (2010). After developing the new model, we find the basic reproduction number by using the well known approach i.i. next generation matrix (NGM) approach. Furthermore, we investigate the possible equilibriums i.e. disease free and endemic equilibria and show the local asymptotic stability as well as global asymptotic stability at both equilibriums. For the local asymptotic stability, we use linearization and Routh-Hurwitz criteria, while to discuss the global stability, we then use the classic Lypanavo function theory and geometrical approach. Finally the numerical simulation are achieved by exhausting Runge-Kutta method of order fourth scheme to show the feasibility of the obtained results and the role of saturated incidence rate.
The organization of the paper is as follows. In the second section, we presents the proposed model and studied its different analysis, including basic reproduction number and equilibriums. In section "Stability analysis", we study the stability analysis and prove the local as well as global stability. Section "Numerical analysis" is devoted to numerical simulation and discussion. Finally a brief conclusion is presented in section "Conclusion".

Mathematical model and its analysis
Herein this section, we presents a Hepatitis B virus transmission epidemic model. For this, we split the entire populace into seven epidemiological subclasses, susceptible S(t), latent L(t), acute infected A(t), chronic infected individuals B(t), carrier individuals C(t), recovered with permanent immunity R(t) and vaccinated V(t). Thus the compartmental mathematical model can be represented by the succeeding system of seven ordinary differential equations, Here, b represents the birth rate, ξ represents the birth rate without successful vaccination, η represents the proportion of prenatally infected individuals, φ represents the rate of waning vaccine induced immunity, β represents the transmission rate from susceptible to infected, γ, ζ represent the reduced transmission rate of chronic and carrier Hepatitis B infected individuals, respectively. µ 0 represents the death rate, which occur naturally, v represents the vaccination rate, σ represents the moving rate from latent class to acute class, γ 1 represents the moving rate from acute to chronic and carrier, γ 2 represents the moving rate of chronic carrier to immune, γ 3 represents the moving rate of carrier to immune, µ 1 , µ 2 represents the death rate, which occur from the Hepatitis B and p represents the average probability of those individuals, who fails to recovers in acute stage and goes to chronic carrier.
Let "N (t)" represents the entire populace, such that, The initial conditions with the proposed model, make sure that, N (t) ≥ 0. Hence the entire populace "N (t)" is bounded and persist positive for t > 0. Now the time derivative of N(t) yields which implies that, for large time t, that is t → ∞, N (t) ≤ b µ 0 . Thus the feasible region for our proposed model is Therefore, we study the dynamic of our proposed model (1) in the feasible region , which is an attracting set. For obtaining the Jacobian matrix of our model, we take the reduced system, because R dose not appears explicitly in all others classes, so the Jacobian matrix of the reduced system (with out R) becomes where

Basic reproduction number
In epidemiological models the inception capacity R 0 is called the basic reproduction number is a key concept, which characterizes the anticipated average amount of new contaminations created directly and circuitously by a solitary infectious, when familiarized into a entirely susceptible populace (Anderson and May 1991;van den Driessche and Watmough 2002). To find the this quantity for our proposed model (1), we use the method of Van Den Driessche and Watmough (2008). Let χ = (L(t), A(t), B(t), C(t)), so from the model (1), we have In Eq. (5), F and V are define as . Now we find the Jacobian of F and V at diseases free equilibrium F 0 , thus we have where a 11 = σ + µ 0 , a 22 = µ 0 + γ 1 + ψ, a 32 = pγ 1 , a 33 = µ 0 + µ 1 + γ 2 , a 41 = (1− p)γ 1 and a 44 = µ 0 + µ 0 + γ 3 − bξη. Thus R 0 is the spectral radius of K = FV −1 , that . So the basic reproduction number R 0 for our proposed model (1)  becomes where

Equilibrium analysis
To study the dynamics of our proposed model (1), we need to find the equilibriums, thus the disease free equilibrium of the model (1) is denoted by F 0 and define as Similarly endemic equilibrium is denoted by F 1 and define as and w 1 = σβbξ ηγ 1 q 5 (q 4 + γ 1 pq 6 + γ 1 q 2 q 3 q 4 q 5 (αq 7 + βq 6 ) + αβbξ ηγ 2 1 q 4 q 6 ).

Stability analysis
Herein this section, we want to study asymptotic stability of our proposed model at both equilibria. So first, we prove the local stability of our model at illness free and endemic equilibrium, then the global asymptotic stability. For the local stability, we use the linearization method, but for global dynamic stability, we want to use both the Lypanavo function theory and geometrical approach.

Local stability analysis
For the local dynamic of the proposed model at disease free and endemic equilibrium points, we state and prove the following results.
Theorem 1 If R 0 < 1, then the model (1) is locally asymptotically stable at disease free equilibrium point F 0 and if R 0 > 1, then it is the unstable saddle point.
Proof The characteristic equation of the Jacobian matrix J (4) at disease free equilibrium F 0 becomes where If R 0 < 1, we have 0 <R i < 1, for i = 1, 2, 3. So a j > 0, for j = 1, 2, 3, 4, and also it is easy to show that, a 1 a 2 a 3 > a 2 3 + a 2 2 a 4 . Therefore Routh-Herwitz criteria is satisfied, that is all the roots of the characteristic polynomial P( ) have negative real parts, which ensure that the disease free equilibrium point F 0 is stable, while for R 0 > 1, the system have both negative and positive eigenvalues, which states that the disease free equilibrium point is a saddle point, always unstable.
Theorem 2 If R 0 > 1, then the model (1) is locally asymptotically stable at endemic equilibrium point F 1 and if R 0 < 1, then it is the unstable.

Global stability analysis
To show the global stability of the proposed model, using Lypanavo function theory and geometrical approach. Here, we want to prove the global stability of disease free equilibrium point F 0 by using Lypanavo function theory, while to prove the global stability at endemic equilibrium point F 1 , we use geometrical approach. Before to presents the stability results, the following important lemma due to Li and Muldowney (1996) is quoted. So, we obtain the condition at which the proposed model is globally asymptotically stable.
Lemma 1 If the system ẏ = g(y), g : D → R n containing a unique equilibrium of the form y * and there is a compact absorbing set, then this system is globally asymptotically stable around that equilibrium, if there exist a function D(y) and a Lozinskii measure ℓ, such that lim t→∞ supsup 1 t t 0 ℓ(B)dt < 0.
Thus regarding the global stability of the proposed model at disease free and endemic equilibrium points, we have the following results.
Theorem 3 If R 0 < 1, then the model (1) is globally asymptotically stable at disease free equilibrium point F 0 and unstable otherwise.
Theorem 4 If R 0 > 1, then the model (1) is globally asymptotically stable at endemic equilibrium point F 1 .
Proof Let J 2 and J 3 be the Jacobian matrix and second additive compound matrix containing only the first three equation of the model (1), then we have Let us consider the function P(χ) = P(S, L, A) = diag S L , S L , S L , which implies that P −1 (χ) = diag L S , L S , L S , then taking the time derivative, that is P f (χ), we get Thus we take B = P f P −1 + PJ 3 P −1 , which can be written as (16) In the above matrix ) be a vector in R 3 and its norm . defined by Now we take the Lozinski measure ℓ(B) with respect to the above norm described by Martin (1974), that is ℓ(B) ≤ sup{g 1 , g 2 } = sup{ℓ(B 11 ) + �B 12 �, ℓ(B 22 ) + �B 21 �}, where g i = ℓ(B ii ) + �B ij � for i = 1, 2 and i � = j, which implies that {v, σ } and �B 21 � = max{σ , 0} = σ . Therefore g 1 and g 2 becomes Thus, we can write Hence ℓ(B) ≤Ṡ S − 2µ 0 . Now integrating the Lozinski measure ℓ(B) with respect to t in the interval [0, t] and taking lim t→∞ , we obtain So finally, we can write Thus the subsystem, which containing first three equations of the model (1) is globally asymptotically stable around its interior equilibrium (S 1 , L 1 , A 1 ). Now consider the subsystem of the model (1), such that Taking the limit system of the model (22), we get Solving system (38) and using the initial conditions B(0), C(0), R(0) and V(0). So for large time t that is t → ∞, B(t) → B 1 , C(t) → C 1 , R(t) → R 1 and V (t) → V 1 , which is sufficient to prove that the endemic equilibrium point E 1 is globally asymptotically stable.

Numerical analysis
In this section, we want to observe the dynamical behavior of our proposed model. In order to do this, we purpose numerical results by using Runge-Kutta of order 4th scheme which have used several authors for a wide range of problems consisting of ordinary differential equations. For the simulation purpose, we use different value of parameters used in the proposed model are given in the Table 1. In the set of parameters some are taken from published articles, while some of the parameters are taken in such a way, α Saturation 0-0.9 Assumed that would be much more biologically feasible. Furthermore the time interval is taken 50 months with initial population for susceptible, latent, acute infected, chronically infected, carrier, recovered and vaccinated individuals as 100, 80, 60, 40, 20, 0 and 20, respectively. In the study of biological dynamic especially in the transmission dynamic of infectious disease sensitivity analysis is very important. Because by sensitive analysis, we are able to judge the role of every parameter and thus easily, we are in the position to develop a strategy for eradicating the disease from the community. Here we carry out a sensitivity analysis of our proposed model (1) with respect to the parameter α, which are presented in Fig. 1.
In Fig. 1, the sensitivity analysis of susceptible, latent, acute infected with Hepatitis B, chronic infected with Hepatitis B, carrier, recovered and vaccinated are presented, which shows that, susceptible, recovered and vaccinated population are directly proportional to saturation α, while latent, acute infected, chronically infected and carrier population are inversely proportional to saturation α. So from this brief analysis, it is observed that, saturation is also one of the factor to prevent the spreading of Hepatitis B in the world, which is one of the top three infectious diseases. Almost one third of the population are infected with Hepatitis B. Figure 2 represents the phase space diagram of susceptible, latent, recovered and susceptible, latent, vaccinated with stable endemic equilibrium and respective set of parameters given in Table 1. Figure 3 represents the phase space diagram of latent, acute, chronically infected and latent, acute, recovered. Similarly Figs. 4 and 5 represents the phase space diagram of carrier, recovered, vaccinated, time, recovered, vaccinated, time, latent, carrier, time, latent and recovered respectively with stable endemic equilibrium and respective set of parameters given in Table 1.

Conclusion
The contagious disease of Hepatitis B is one of the topmost three infectious diseases reported by WHO and the incidence rate is precise significant in the study of epidemic models. In this article, we developed an epidemic problem for the spreading dynamic of Hepatitis B virus by classification of different subclasses with saturated occurrence rate. This occurrence rate is more generalized from than bilinear incidence rate and seems more realistic. Therefore we divided the host population into seven epidemiological subclasses, namely, susceptible, latent, acute infected with Hepatitis B, chronically infected with Hepatitis B, carrier, recovered and vaccinated, then formulated the model with this new features. First, we find the basic reproduction number R 0 . As in epidemiological models, the mathematical models has two steady states, infected and noninfected steady states. So we obtained the two non-negative equilibria that is disease free F 0 , which is always exist and locally stable for R 0 < 1 and endemic equilibria F 1 , which exist only,  Fig. 1 The plot shows the dynamical behavior of the proposed model (1) with saturation and without saturation if R 0 > 1, therefore locally stable, if R 0 > 1. This portents occur, since totally the flows along the axes of susceptible, vaccinated and recovered are always attractor to the point F 0 , but the flow along the axes of the infected compartment that is latent, infected with acute, chronic Hepatitis B and carrier depend on the value of R 0 . So if R 0 < 1, then the axes of infected compartments are attractor towards the illness free equilibrium F 0 . But when the R 0 acroses one that is R 0 > 1, the axes of these compartment do not attract the disease free equilibrium F 0 and repels from it. Further more to show the global stability at F 0 , we established the Lyapnovo function, while to to prove the global stability at F 1 , we used the geometrical approach. Finally numerics of the proposed model are retrieved out to show the feasibility of the model.  Fig. 2 The plot shows the phase space diagram of susceptible, latent, recovered and susceptible, latent and vaccinated with stable endemic equilibrium F 1 and respective set of parameters given in Table 1  0 Fig. 3 The plot shows the phase space diagram of latent, acute infected, chronically infected and latent, acute infected, recovered with stable endemic equilibrium F 1 and respective set of parameters given in Table 1 0  Fig. 4 The plot shows the phase space diagram of carrier, recovered and vaccinated with stable endemic equilibrium F 1 and respective set of parameters given in Table 1 Fig. 5 The plot shows the phase space diagram of time, susceptible, vaccinated, time, latent, carrier and time, latent, recovered with stable endemic equilibrium F 1 and respective set of parameters given in Table 1