Discrete Dynamics in Nature and Society

Volume 2015, Article ID 563127, 11 pages

http://dx.doi.org/10.1155/2015/563127

## Stability Analysis for a Fractional HIV Infection Model with Nonlinear Incidence

^{1}Department of Basic Course, Haikou College of Economics, Haikou 571127, China^{2}School of Mathematics and Physics, China University of Geosciences, Wuhan 430074, China^{3}Department of Mathematics and Statistics, University of Helsinki, 00014 Helsinki, Finland

Received 30 September 2014; Accepted 16 December 2014

Academic Editor: Bai-Lian Larry Li

Copyright © 2015 Linli Zhang 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

We introduce the fractional-order derivatives into an HIV infection model with nonlinear incidence and show that the established model in this paper possesses nonnegative solution, as desired in any population dynamics. We also deal with the stability of the infection-free equilibrium, the immune-absence equilibrium, and the immune-presence equilibrium. Numerical simulations are carried out to illustrate the results.

#### 1. Introduction

##### 1.1. The History of Fractional Calculus

Fractional calculus which is a branch of mathematical analysis extends derivatives and integrals to an arbitrary order (real, even, or complex order). The study on fractional calculus started at the end of seventeenth century and the first reference was proposed by Leibniz and L’Hospital in 1695, in which half-order derivative was mentioned. The emergence of fractional calculus aroused interest of some famous mathematicians such as Lacroix, Lagrange, Fourier, and Laplace. J. L. Lagrange developed the law of exponents for differential operators in 1772, which promoted indirectly the development of fractional calculus. P. S. Laplace defined a fractional derivative by means of integral in 1812 and S. F. Lacroix and J. B. J. Fourier in succession mentioned the arbitrary order derivative, respectively, in 1819 and in 1822. However, the basic theory of fractional calculus was established with the studies of Liouville, Grunwald, Letnikov, and Riemann until the end of the nineteenth century.

At the initial stage, the development of fractional calculus had been restricted to the pure mathematics research. Until the nineteen nineties, the theory of fractional calculus was applied to nature and social sciences. In the field of anomalous diffusion, the researchers applied fractional calculus to describe the diffusion processes which do not conform to the brown motion and proposed the fractional anomalous diffusion model. According to the theoretical analysis and the experimental data or results, they found that the fractional model is more reasonable to describe these processes [1]. Afterwards, many mathematicians and applied researchers also have tried to demonstrate applications of fractional differentials in the areas of non-Newtonian fluids [2], signal processing [3–5], viscoelasticity [6, 7], fluid-dynamic traffic model [8], colored noise [9], bioengineering [10], solid mechanics [11], continuum and statistical mechanics [12], and economics [13] and brought new research view for those fields.

##### 1.2. Mathematical Modeling

In biology, it has been deduced that the membranes of cells of biological organism have fractional-order electrical conductance [14]; hence some mathematical models which describe cells behavior are classified into groups of non-integer-order models. In the field of rheology, fractional-order derivatives embody essential features of cell rheological behavior and have enjoyed greatest success [15]. Some researchers also found that fractional ordinary differential equations are naturally related to systems with memory which exists in most biological systems. And fractional-order equations are also closely related to fractals, which are abundant in biological systems.

Mathematical models have been proven valuable in understanding the dynamics of viral infection. Although a large number of works on modeling the dynamics of viral infection have been done [16–20], it has been restricted to integer-order (delay) differential equations. In recent years, it has turned out that many phenomena in virus infection can be described very successfully by the models using fractional-order differential equations [21, 22]. Motivated by these references, in this paper, we will consider a fractional-order HIV model.

Among a large number of virus infection models described with integer-order or delay differential equations, a typical model is given as the following simple three-dimensional system:Here represents the concentration of uninfected cells at time ; represents the concentration of infected cells that produce virus at time ; represents the concentration of viruses at time . and are the death rate of uninfected cells and the rate constant characterizing infection of the cells, respectively. is the death rate of infected cells due to either the virus or the immune system. Free virus is produced from infected cells at the rate and removed at rate .

The infection rate in model (1) is usually assumed to be bilinear with respect to virus and uninfected target cells . However, there are some evidences which show that a bilinear infection might not be an effective assumption when the number of target cells is large enough [23–25]. And Huang et al. [17] incorporated Holling type II functional response infection rate into basic virus dynamics (1), which was expressed as where () is a constant. Hence, system (1) can be modified into the following system:

The immune response following viral infection is universal and necessary to eliminate or control the disease. Antibodies, cytokines, natural killer cells, B cells, and T cells are all essential components of a normal immune response to viral infection. Since cytotoxic T lymphocytes (CTLs) play a critical role in antiviral defense by attacking the virus-infected cells, Huang et al. [17] considered the immune response to model (3) with .

System (3) with can be further simplified if we take into consideration the fact that an average life span of viral particles is usually significantly shorter than one of the infected cells. It can be assumed, therefore, that, compared with a “slow” variation of the infected cells level, the virus load relatively quickly reaches a quasiequilibrium level. The equality holds in the quasiequilibrium state and hence . This assumption is referred to as “separation of time scales” and is in common use in the virus dynamics [18]. We have to stress that this assumption does not imply that the virus concentration remains constant; on the contrary, it is assumed to be proportional to the varying concentration of infected cells . Accordingly, system (3) with can now be reformulated as a system of three differential equations:

Further we introduce the fractional order into model (4) and a new system can be described by the following set of FODEs of order :where the immune response is assumed to get stronger at a rate , which is proportional to the number of infected cells and their current concentration; the immune response decays exponentially at a rate , which is proportional to their current concentration; the parameter expresses the efficacy of nonlytic component.

The paper is organized as follows. In the next section, we give two definitions and two lemmas about fractional calculus. In Section 3, we show the nonnegativeness of the solutions of system (5) with initial condition (6). In Section 4, we give a detailed stability analysis for three equilibria. Finally numerical simulations are presented to illustrate the obtained results in Section 5.

#### 2. Fractional Calculus

*Definition 1 (see [26]). *The Riemann-Liouville (R-L) fractional integral operator of order of a function is defined as Here is the Euler gamma function which is defined as This function is generalization of a factorial in the following form:

*Definition 2 (see [26]). *The Caputo (C) fractional derivative of order , , , is defined as where the function has absolutely continuous derivatives up to order . In particular, when , one has

In this paper we use Caputo fractional derivative definition. The main advantage of Caputo’s definition is that the initial conditions for fractional differential equations with Caputo derivatives take the same form as that for integer-order differential equations.

Lemma 3 (see [26]). *Consider the following commensurate fractional-order system:**with and . The equilibrium points of system (12) are calculated by solving the following equation: . These points are locally asymptotically stable if all eigenvalues of Jacobian matrix evaluated at the equilibrium points satisfy*

*Definition 4 (see [27]). *The discriminant of a polynomial is defined by , where is the derivative of . If , is the determinant of the corresponding Sylvester matrix. The Sylvester matrix is formed by filling the matrix beginning with the upper left corner with the coefficients of and then shifting down one row and one column to the right and filling in the coefficients starting there until they hit the right side. The process is then repeated for the coefficients of .

*Lemma 5 (see [27]). For the polynomial equation,the conditions which make all the roots of (15) satisfy (13) are displayed as follows: (i)for , the condition is ;(ii)for , the conditions are either Routh-Hurwitz conditions or (iii)for ,(a)if the discriminant of , is positive, then Routh-Hurwitz conditions are the necessary and sufficient conditions; that is, , , and if ;(b)if , , , and , then (13) for (15) holds when ;(c)if , , and , then (13) for (15) holds when ;(d)if , , , and , then (13) for (15) holds for all .*

*3. Nonnegative Solutions*

*3. Nonnegative Solutions*

*Let and . For the proof of the theorem about nonnegative solutions, we would need the following lemma.*

*Lemma 6 (see [28] generalized mean value theorem). Let and for . Then one has with , , where .*

*Remark 7. *Suppose that and , for . It is clear from Lemma 6 that if , , then is nondecreasing for each . If , , then is nonincreasing for each .

*We now prove the main theorem.*

*Theorem 8. There is a unique solution for the initial value problem (5) with (6) and the solution remains in .*

*Proof. *From Theorem 3.1 and Remark 3.2 of [29], we obtain the solution on solving the initial value problem (5) with (6) which is not only existent but also unique. Next, we will show the nonnegative orthant is a positively invariant region. What is needed for this is to show that, on each hyperplane bounding the nonnegative orthant, the vector field points into . From (5), we findBy Remark 7, the solution will remain in .

*4. Equilibrium States and Their Stability*

*4. Equilibrium States and Their Stability*

*In this section, we investigate the stability of the fractional-order model of HIV infection of CD4 + T cells, that is, system (5) with (6). Consider the initial value problem (5) with (6) with satisfying .*

*In order to obtain the equilibria of system (5), we set , , and and we haveSystem (5) has three types of relevant nonnegative equilibrium states. System (5) always has an infection-free equilibrium , where which means that the infected cells are cleared. The basic reproductive number of the viruses for system (5) is given byThis number describes the average number of newly generated infected cells from one infected cell at the beginning of the infection process. When , in addition to the infection-free equilibrium , there is an immune-absence equilibrium , where Note that implies . When the HIV infection is in the immune-absence equilibrium , the infected cells and virus exist but the immune response is not activated yet. Further, we denoteNote that , which implies that is equivalent to . The latter is seen as the immune reproductive number, which expresses the average number of activated CTLs generated from one CTL during its life time through the stimulation of the infected cells . It is reasonable that immune response is activated in the case where . When , in addition to the infection-free equilibrium and the immune-absence equilibrium , there is an interior immune-presence equilibrium , where expresses the state where CTLs immune response is present.*

*Next, we establish the local asymptotic stability of model (5) by the characteristic equation.*

*Theorem 9. Consider system (5). (1)The infection-free equilibrium is locally asymptotically stable if .(2)If , the equilibrium is unstable, and if , it is a critical case.*

*Proof. *The characteristic equation for the infection-free equilibrium is given as follows: It is reduced toIt is clear that (26) has the characteristic roots which means , which means , and . Since the imaginary part of characteristic root is zero, which means is necessary and sufficient to ensure the local asymptotic stability of . If , ; hence is unstable. If , , which is a critical case.

*Theorem 10. Consider system (5). (1)The immune-free equilibrium is locally asymptotically stable if .(2)If , the immune-free equilibrium is unstable; if , it is a critical state.*

*Proof. *The characteristic equation for the immune-free equilibrium is given as follows: Using , it is reduced toThe root of (28) is negative and when , positive and when , and zero when , which is a critical case.

Now, we consider the equationDenoteSince and , (29) has two negative real roots and we denote them by and . It is easy to see and . Hence, when , the immune-free equilibrium is locally asymptotically stable, when , is unstable, and when , it is a critical case.

*To discuss the local stability of the immune-present equilibrium state , we consider the linearized system of (5) at . The Jacobian matrix at is given by The characteristic equation of the linearized system iswhereBased on Definition 4, we obtain the discriminant of (32)Using the result (iii) of Lemmas 5 and 3, we have the following theorem.*

*Theorem 11. Consider system (5). Under the condition of , (1)if the discriminant of , is positive, namely, , then the immune-present equilibrium is locally asymptotically stable for .(2)if , then the immune-present equilibrium is locally asymptotically stable for .*

*5. Numerical Method*

*5. Numerical Method**Atanackovic and Stankovic introduced a numerical method to solve the single linear FDE in 2004 [1]. A few years later, they developed again a method to solve the nonlinear FDE [30]. It was shown that the fractional derivative of a function with order satisfying may be expressed aswherewith the following properties:We approximate by using terms in sums appearing in (35) as follows:We can rewrite (38) as follows:where We set for . We can rewrite system (19) as the following form:whereNow we can rewrite (39) and (42) as the following form:with the following initial conditions:Now we consider the numerical solution of system of ordinary differential equations (44) with the initial conditions (45) by using the well-known Runge-Kutta method of fourth order.*

*6. Numerical Simulation and Discussion*

*6. Numerical Simulation and Discussion**Firstly, by using GEM (generalized Euler method) [21], we simulate system (5) with the parameter values as shown in Table 1.*