GLOBAL ANALYSIS OF WITHIN HOST VIRUS MODELS WITH CELL-TO-CELL VIRAL TRANSMISSION. Hossein Pourbashash and Sergei S. Pilyugin and Patrick De Leenheer
|
|
- Mildred Nichols
- 5 years ago
- Views:
Transcription
1 Manuscript submitted to AIMS Journals olume X, Number0X, XX200X doi:0.3934/xx.xx.xx.xx pp. X XX GLOBAL ANALYSIS OF WIHIN HOS IRUS MODELS WIH CELL-O-CELL IRAL RANSMISSION Abstract. Recent experimental studies have shown that HI can be transmitted directly from cell to cell when structures called virological synapses form during interactions between cells. In this article, we describe a new withinhost model of HI infection that incorporates two mechanisms: infection by free virions and the direct cell-to-cell transmission. We conduct the local and global stability analysis of the model. We show that if the basic reproduction number R 0, the virus is cleared and the disease dies out; if R 0 >, the virus persists in the host. We also prove that the unique positive equilibrium attracts all positive solutions under additional assumptions on the parameters. Finally, a multi strain model incorporating cell-to-cell viral transmission is proposed and shown to exhibit a competitive exclusion principle. Hossein Pourbashash and Sergei S. Pilyugin and Patrick De Leenheer Department of Mathematics, University of Florida 400 Stadium Road Gainesville, FL 326, USA Connell McCluskey Department of Mathematics, Wilfrid Laurier University 75 University Avenue West Waterloo, ON, Canada, N2L 3C5 his paper is dedicated to Chris Cosner, on the occasion of his 60th birthday. (Communicated by the associate editor name). Introduction. Mathematical modeling of within host virus models has flourished over the past few decades. hese models have been used to describe the dynamics inside the host of various infectious diseases such as HI, HC, HL, as well as the flu or even the malaria parasite. esting specific hypotheses based on clinical data is often difficult since samples cannot always be taken too frequently from patients, or because detection techniques of the virus may not be accurate. his justifies the central role played by mathematical models in this area of research. One of the early models, known as the standard model and reviewed in the next section, was used by Nowak and May [6], and by Perelson and Nelson [7] to model HI. It was successful in numerically reproducing the dynamics of the early stages of HI and its target, the CD4 + cells, following an infection event. he global behavior of the standard model was first investigated analytically in [2]. Here, the global stability of the disease steady state was proved using the powerful second compound matrix methods developed by Muldowney [5] and by Li and Muldowney [0]. Among the first successes of this geometrical method was the long 200 Mathematics Subject Classification. Primary: 37N25, 97M60; Secondary: 34A34, 34D23. Key words and phrases. within-host virus model; cell-to-cell transmission; global stability; Lyapunov function, second compound matrices.
2 2 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER sought-after proof of global stability in the classical epidemiological SEIR model []. Since epidemic models and within host virus models are structurally similar, it should not come as a surprise that the geometrical methods of [0, 5] are at the basis of the global stability proof of the standard model in [2]. It is well-known from the converse Lyapunov theorem [7], that whenever a nonlinear system of ODEs has a globally stable steady state, that then there exists a global Lyapunov function. Korobeinikov in [8] was the first to discover a Lyapunov function in the context of SEIR and SEIS models. Korobeinikov s Lyapunov function was exploited further in [3] to establish global stability for the standard model under relaxed conditions for the model parameters, as well as for within host virus models with multiple, possibly mutating, strains. Although our discussion so far was focused on the global stability of certain steady states, it has been demonstrated that more complex dynamical behavior, in particular taking the form of sustained oscillations, may arise in the standard model for specific proliferation rates of healthy cells [2, 24], or under the assumption that infected cells proliferate as well [23], a feature which is not encoded in the standard model. Based on more recent findings about additional infection processes, it is the main goal of this paper to extend the standard model in two ways:. In the standard model it is assumed that healthy cells can only be infected by viruses. However, recent work shows that cell-to-cell transmission of viruses also occurs when a healthy cell comes into contact with an infected cell [8]. In fact, in vitro experiments reported in [2] have shown that in shaken cultures, viral transfer via cell-to-cell contact is much more rapid and efficient than infection by free virus because it avoids several biophysical and kinetic barriers. In vivo, direct cell-to-cell transmission is also more potent and more efficient according to [4]. he details of the infection mechanism via cellto-cell contact are described in [3, 4, 9] and have been attributed to the formation of virological synapses (Ss), filopodia and nanotubes. he model proposed here deviates from, and in fact, expands the standard model, by including the additional cell-to-cell infection mechanism. 2. According to the standard model, viruses do not enter previously infected cells, although most HI infected splenocytes for example, are found to be multiply infected [4]. Moreover, the loss of virus upon entry of a healthy cell, has often been neglected in the literature. he model we propose here, accounts for both processes. 2. he mathematical model. he standard model of a within-host virus infection [6, 7] is = f( ) k, = k β, = Nβ ( k ), where the term in the parentheses may or may not be present. Our mathematical model has the following form: = f( ) k k 2, = k + k 2 β, = Nβ k 3 k 4. ()
3 WIHIN HOS IRUS MODELS 3 where, and denote the concentrations of uninfected host cells, infected host cells and free virus particles, respectively. Parameters k,k 2,β,N and are all positive constants. k is the contact rate between uninfected cells and viruses. k 2 is the contact rate between uninfected cells and infected cells. he parameters β and represent the death rate of infected cells and virus particles, respectively. In the case of budding viruses, N is the average number of virus particles produced by an infected cell during its lifetime; in case of lytic viruses, N represents the average burst size of an infected cell. he parameter k 3 models the rate of absorption of free virions into healthy cells during the infection process. In some instances, this term has been neglected [6, 7]. he parameter k 4 models the absorption of free virions into already infected cells, a process which may be biologically relevant for some viruses, yet in each case it will be determined by the specific mechanism of the viral cell entry [4]. he growth rate of the uninfected cell population is modeled by the smooth function f : R + R, which is assumed to satisfy the following: 0 > 0:f( )( 0 ) < 0, = 0, and f ( ) < 0 [0, 0 ]. (2) he continuity of f implies that f( 0 ) = 0, and hence E 0 =( 0, 0, 0) is an equilibrium point of system (). Biologically, E 0 represents the disease-free equilibrium. he Jacobian matrix of () at E 0 is J 0 = f ( 0 ) k 2 0 k 0 0 k 2 0 β k 0. 0 Nβ k 3 0 he submatrix of J 0 corresponding to the infectious compartments can be split as follows: k2 J 0 = 0 β k 0 k2 = 0 k 0 β 0 =: F. Nβ k Nβ + k 3 0 A quick calculation shows that the next generation matrix is given by F = k2 0 β + Nk0 k 0 +k 3 0 +k he basic reproductive number, R 0, is defined as the spectral radius of F [22], hence R 0 := ρ(f )= k 2 0 β + Nk 0. (3) + k 3 0 One eigenvalue of J 0 is given by f ( 0 ) < 0, and the remaining two are also eigenvalues of J 0. he determinant of J 0 is given by det J 0 = (β k 2 0 )( + k 3 0 ) Nβk 0 = β( + k 3 0 ) k 2 0 β Nk 0 + k 3 0 = β( + k 3 0 )( R 0 ). Consequently, if R 0 >, then det J 0 < 0, and J 0 has a positive eigenvalue. If R 0 <, then det J 0 > 0, and k 2 0 <βwhich implies that trj 0 < 0, so that all eigenvalues of J 0 have negative real parts. We obtain our first result. heorem 2.. If R 0 <, then E 0 is locally asymptotically stable; if R 0 >, E 0 is unstable.
4 4 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER A positive equilibrium E =(,, ) must satisfy the following relations: β = f( ), = Nf( ) + k 3 + k4, =k 2 β + Nk + k 3 + k 4f( ) β := G( ). he function G is monotonically increasing, since the first term is linear in, and the reciprocal of the second term + k 3 + k 4f( ) Nk Nk Nk β is strictly decreasing in when is in (0, 0 ), by (2). A positive equilibrium corresponds to a solution of G( )=with0< < 0. Since G(0) = 0 and G( 0 )=R 0, we conclude that a unique positive equilibrium exists if and only if R 0 >. Lemma 2.2. If R 0 the equilibrium E 0 is the only equilibrium of (), and if R 0 > then E 0 and E =(,, ) are the only two equilibria of (). 3. Preliminary results for the general model. Lemma 3.. Every forward solution of () starting in R 3 +, is defined in R 3 + for all t 0, and bounded. Furthermore, system () is uniformly dissipative in R 3 +. Proof. Positive invariance of R 3 + follows from the standard argument that the vector field of system () is inward pointing on the boundary of R 3 +. o show boundedness, fix a forward solution in R 3 +. Observe that the -equation in () implies that (t) M := max( 0,(0)) for all t 0 as long as the solution is defined. Since f( ) is bounded on [0,M], there exist A, B > 0 such that f( ) A B for all [0,M]. Choose a constant 0 <α < /N and consider an auxiliary function U = + + α.dueto(), we have U f( ) β( α N) α A B β( α N) α A α 2 U, where α 2 =min(b,β( α N),) > 0. Hence, U(t) max U(0), Aα2 for all t 0. his implies that the corresponding forward solution of () is bounded in R 3 + and exists for all t 0. From the above argument, it is clear that any solution of () inr 3 + eventually satisfies (t) M 0 := 0 +, and therefore one can choose uniform constants A 0,B 0 > 0 so that the function U(t) eventually satisfies the inequality U A 0 α 0 U where α 0 =min(b 0,β( α N),) > 0. Hence, any solution of () inr 3 + eventually enters a forward invariant simplex D = (,,) R 3 + : + + α A 0 +. (4) α 0 Hence, system () is uniformly dissipative in R 3 +. At this point, we make the following (possibly trivial) remark: Consider any solution of () inr 3 +. If there exists a t 0 such that (t 0 ) < 0 then (t) < 0 for all t t 0. Our next result provides the necessary and sufficient condition for the global stability of E 0. heorem 3.2. If R 0 then the infection free equilibrium E 0 attracts all solutions in R 3 +.
5 WIHIN HOS IRUS MODELS 5 Proof. If (0) = 0, then (0) > 0, so (t) > 0 for all sufficiently small t>0. If (t) = (t) = 0, then it is clear that (t) 0 as t. If exactly one of (0) or (0) is zero, then its derivative at t = 0 is positive, and hence this variable is positive for all small t>0. Furthermore, none of (t), (t) or (t) can decrease to zero in finite time, hence it suffices to consider only solutions with both (t) and (t) positive for all t 0. We previously remarked that either (t) 0 for all t 0, or there exists a t 0 0 such that (t) < 0 for all t>t 0. Consider the former case, that is, (t) 0 for all t 0. hen (t) 0 for all t 0, and (t) 0 as t. Barbalat s Lemma implies that lim t (t) = 0. Noticing that by (2), f( ) f( 0 ) = 0, we use the -equation in (), to conclude that = 0, and lim t (t) =lim t (t) = 0. Now, suppose that there exists a t 0 0 such that (t) < 0 for all t>t 0. Without loss of generality, we may assume that (t) < 0 for all t 0. We introduce an auxiliary function: and evaluate Ẇ along a solution of (): Ẇ = ( + k 3 0 ) + k 0 W =( + k 3 0 ) + k 0, (5) = ( + k 3 0 )(k + k 2 β )+k 0 (Nβ k 3 k 4 ) = (k 0 Nβ +( + k 3 0 )(k 2 β)) + (k ( + k 3 0 ) k 0 ( + k 3 ) k 0 k 4 ) β( + k 3 0 ) (R 0 ) + k ( 0 ) k k 4 0 0, where all three terms are non-positive. herefore, there exists lim t W (t) = W 0. Barbalat s Lemma then implies that lim t Ẇ (t) = 0. If W = 0, it is clear that lim t (t) =lim t (t) = 0 and lim t (t) = 0. If W > 0, then there exists > 0 such that k + k 2 > 0 for all sufficiently large t. he -equation in () implies that f( ), and thus there exists 2 > 0 such that (t) 0 2 for all sufficiently large t. Since Ẇ k ( 0 ) k 2,lim t Ẇ (t) = 0 implies that lim t (t) = 0, and by Barbalat s Lemma, we have that lim t (t) = 0, and the -equation in () implieslim t (t) = 0, hence W = 0, a contradiction. herefore, W = 0, and the proof is complete. he next heorem shows that the positive equilibrium E is locally stable whenever it exists. heorem 3.3. If R 0 >, then E is locally asymptotically stable. Proof. If R 0 >, a unique positive equilibrium E exists by Lemma 2.2. For notational convenience, we will drop all bars so that all variables are at their equilbrium values. he following relations hold at E: and β = f( ), k + k 2 = f( )/, + k 3 + k 4 = Nβ β k 2 = k. = Nf( ),
6 6 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER he Jacobian at E has the form f ( ) k k 2 k 2 k J(E) = k + k 2 k 2 β k. k 3 Nβ k 4 k 3 k 4 Using the equilibrium relations, J(E) simplifiesto f f/ k 2 k J(E) = f/ k k, k 3 ( + k 3 ) Nf where we have suppressed the arguments of f and f. he characteristic polynomial p(λ) =λ 3 + a λ 2 + a 2 λ + a 3 is given by λ f + f/ k 2 k p(λ) =det f/ λ + k k. k 3 ( + k 3 ) λ + Nf A direct computation using the equilibrium relations above to simplify a 2 and a 3, we obtain a = f/ f + k + Nf > (f/ f )+ Nf > 0, a 2 = k k 4 +( f k + f/ ) + Nf + k 2 f k k 3, a 3 = ( f + f/ )k k 4 + f/ Since f > 0, we have that Since kf a 2 >k k 4 + k f = k β, Nf k 2 Nf + k ( + k 3 ) + Nf2 + k 2f k k 3. >k 3, and f/ > k, we have that a 2 >k k 4 + k β + k 2 f>0. Simplifying a 3 in a similar fashion, we find that k k 3 β. a 3 =( f + f/ )k k 4 + k 2Nf 2 + k β > 0. Finally, since Nf >, we have that a a 2 >a 3, and J(E) satisfies the Routh- Hurwitz stability conditions. Lemma 3.4. System () is uniformly persistent in D, the interior of D, ifand only if R 0 >. Here, the set D is defined in (4). Proof. If R 0, then by heorem 3.2, E 0 attracts all positive solutions. Hence R 0 > is necessary for persistence. he sufficiency of the condition R 0 > follows from a uniform persistence result, heorem 4.3 in [6]. o demonstrate that system () satisfies all the conditions of heorem 4.3 in [6] whenr 0 >, choose X = R 3 + and E = D. he maximal invariant set N on the boundary D is the singleton {E 0 } and it is isolated. hus, the hypothesis (H) of [6] holds for (). We conclude the proof by observing that, in the setting of (), the necessary and sufficient condition for uniform persistence in heorem 4.3 of [6] is equivalent to E 0 being unstable. In the remainder of this paper, we prove that the positive equilibrium E is globally asymptotically stable in several special cases of system ().
7 WIHIN HOS IRUS MODELS 7 4. he case k 4 =0. We begin this section by introducing the method of demonstrating global stability through the use of compound matrices; see [0, 5]. Let Ω be an open subset of R n, and let f : Ω R n be a vector field that generates the flow φ(t, x). hat is, for a given x Ω, z(t) =φ(t, x) istheunique solution of the system ż = f(z) (6) with z(0) = x. A compact set K Ω is called absorbing if for each compact set F Ωthereexistsa>0suchthat φ(t, F ) K for all t. It is clear that an absorbing set must contain all equilibria of f in Ω. We say that an open set Ω R n is simply connected if each closed curve in Ω can be contracted to a singleton within Ω. Let denote the exterior product on R n R n, that is, a surjective antisymmetric bilinear transformation : R n R n R d,whered = n 2. Consider the basis of R d consisting of pairwise products e i e j, i<j n ordered lexicographically, where the vectors e i, i n are the elements of the canonical basis of R n. For a real n n matrix M, its second additive compound matrix, denoted by M [2],is defined by its action on the chosen basis as follows: M [2] (e i e j ):=(Me i ) e j + e i (Me j ), i<j n. According to this definition, M [2] is a d d matrix. For instance, if M =(m ij )is a3 3 matrix, then its second additive compound is given by M [2] = m + m 22 m 23 m 3 m 32 m + m 33 m 2. (7) m 3 m 2 m 22 + m 33 For more details on compound matrices, we refer the reader to [5, 5]. Let denote a vector norm on R d and let be the induced matrix norm on R d d, the space of all real d d matrices. For a given matrix A R d d,wedefine its Lozinskii measure with respect to the norm as µ(a) = I + ha lim. h 0 + h More details on the properties of Lozinskii measures can be found in []. For future reference we recall without proof the Lozinskii measure, associated to the l vector norm: Lemma 4.. he Lozinskii measure of a real n n matrix M with respect to the matrix norm induced by the l vector norm, is given by µ (M) = max i m ii + j=i m ji Let the map A :Ω R d d be continuously differentiable and nonsingular for all x Ω, and let µ be a Lozinskii measure on R d d,whered = n 2. Let K be a compact absorbing set for the flow φ(t, x) in Ω. Define a quantity q 2 as where q 2 =limsupsup x K t t t 0 µ(m(φ(s, x)))ds, (8) M(x) =A f (x)a (x)+a(x)j [2] (x)a (x), (9)
8 8 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER with A f (x) := d A(φ(t, x)), dt We have the following result [0, 24]. J(x) := f x (x). heorem 4.2. Consider system (6), and suppose that Ω is simply connected and contains a compact absorbing set K. Let (6) have a unique equilibrium point x K. If A :Ω R d d is C and everywhere nonsingular, and if µ is a Lozinskii measure such that q 2 < 0, then x is globally asymptotically stable in Ω. For the proof of the global stability of E, we will need the following Lemma. Lemma 4.3. Every forward solution of () in R 3 + satisfies lim sup (t) Nβ f(0) t β + 0 =: max. Proof. Consider a nonnegative forward solution of () and let >0. he -equation in () implies that (t) 0 + for all sufficiently large t. Hence, since f ( ) < 0 on [0, 0 ]by(2), the total number of cells, tot := +, for all large t satisfies: tot = f( ) β f(0) + β β tot f(0) + β( 0 + ) β tot. It follows that (t) tot (t) f(0) β for all sufficiently large t. he -equation in () now implies that f(0) Nβ β for all sufficiently large t, thus (t) Nβ f(0) β for all sufficiently large t. Since >0 is arbitrary, the claim follows. heorem 4.4. Suppose that k 4 =0and k 3 0. Let D be the set defined in (4), and let Ω= D. SupposethatR 0 >, 2k 2 0 <β,andk 3 (f(0) + β 0 ) <β, then E is globally asymptotically stable for system () in Ω. Proof. Consider a positive solution of (). As in the proof of heorem 3.2, any solution with (t) 0 for all t 0 converges to E 0. Since R 0 >, Lemma 3.4 implies persistence, hence any positive solution must satisfy (t) < 0 for all sufficiently large t. In addition, Lemma 4.3 implies that lim sup t (t) max. he condition k 3 (f(0) + β 0 ) <βimpliesthat k 3 max <Nβ,hencek 3 (t) <Nβ for all sufficiently large t. By shifting time (which does not affect the value of the quantity q 2 ), if necessary, we will assume that (t) < 0 and k 3 (t) <Nβfor all t 0. Now, consider the Jacobian matrix of system (): f ( ) k k 2 k 2 k J = k + k 2 k 2 β k. (0) k 3 Nβ k 3 he second compound matrix corresponding to J is given by J [2] = a β k k Nβ b k 2, () k 3 k + k 2 c β
9 WIHIN HOS IRUS MODELS 9 where a = f ( ) k k 2 + k 2, b = f ( ) k k 2 k 3, c = k 2 k 3. (2) Define an auxiliary matrix function A on Ω as A := diag,,. Since, > 0everywhereinΩ,A is smooth and nonsingular. Furthermore, we find that A f = diag ( ) 2, 2, 2, A f A = diag,,, hence a β AJ [2] A = Nβ M = A f A + AJ [2] A = where M = a β M 22 = k k b k 2 k 3 k + k 2 c β = a β M 2 Nβ k, (3) k b k 2 k 3 k + k 2 c β M M 2, (4) M 22, M 2 = k b k + k 2 k, M2 = k 2 c β. Nβ k 3, Let (u, v, w) denote a vector in R 3, and define a vector norm on R 3 as (u, v, w) = max( u, v + w ) and let µ be the corresponding Lozinskii measure. hen from [, 2], we have that where µ(m) max(g,g 2 ), (5) g = µ (M )+M 2, g 2 = µ (M 22 )+M 2. M 2, M 2 are the operator norms associated to the linear mappings M 2 : R 2 R, and M 2 : R R 2 respectively, where R 2 is endowed with the l vector norm in both cases.
10 0 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER We have the following µ (M ) = a β, M 2 = k, M 2 = max Nβ,k 3, µ (M 22 ) = + max (k + k 2 + b, k 2 + c β), where Lemma 4. has been used to calculate µ (M 22 ). Since k 3 < Nβ, it follows that Nβ k 3, and M 2 = Nβ. Using (), we find that = k = Nβ Recalling the expressions for a, b and c in (2), we have that + k 2 β, (6) k 3. (7) µ (M 22 )= Nβ + max(f ( ), 2k 2 β), and therefore that g = a β + k = f ( ) k k 2, g 2 = max(f ( ), 2k 2 β). Since f ( ) k k 2 max(f ( ), 2k 2 β), we find that g g 2, and then (5) implies that µ(m) g 2. Let p = max [0, 0] f ( ), and p 2 =2k 2 0 β. he assumptions in 2 imply that p < 0, and the assumptions of this heorem imply that p 2 < 0. Hence, µ(m) p, p := max(p,p 2 ) < 0. Integrating the above inequality, we find that herefore, t lim sup t t t 0 t 0 µ(m) ds p. µ(m) ds p<0. his clearly holds for all solutions in Ω, thus q 2 p<0. he uniform persistence established in Lemma 3.4 implies the existence of a compact absorbing set K Ω. Observing that Ω is simply connected, we apply heorem 4.2 to complete the proof. 5. he case k 3 =0. Although this case is biologically relevant only if k 4 = 0, the following global stability result holds for all k 4 0. heorem 5.. Suppose that k 3 =0and k 4 0. Let D and Ω be as defined in heorem 4.4. IfR 0 >, then the equilibrium E is globally asymptotically stable for system () in Ω.
11 WIHIN HOS IRUS MODELS Proof. Local stability of E follows from heorem 3.3. equilibrium relations hold: k + k2 Combining (9) and (20), we find that Consider the following function on R 3 +: W = dτ + τ τ So, dw dt = d dt + he first term, A,inẆ A = = = Recall that the following = β, (8) = k, β k 2 (9) Nβ = + k 4. (20) Nβk = β k2 + k k 4. (2) d dt + k can be rewritten as (f( ) k k 2 ) (f( ) f( )) + (f( ) f( )) + dτ + k τ dτ. d dt := A + A 2 + A 3. f( ) (k + k 2 )+ (k + k 2 ) (k + k2 ) (k + k 2 ) + (k + k 2 ) = (f( ) f( )) + (k + k2 ) 2 (k + k 2 ) (k + k 2 )+ (k + k 2 ). Due to (8), the second term, A 2,inẆ takes the form A 2 = (k + k 2 β )= k + k 2 β k he third term, A 3,inẆ is A 3 = k (Nβ k 4 )= k Nβ k k k 4 Using (2), A 3 can be written as k Nβ k 2 + k + k2. + k + k k 4. A 3 = β k 2 + k k 4 k k 4 k β + k 2 k k k + k k 4.
12 2 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER Adding A + A 2 + A 3, and using (8), we obtain Ẇ = (f( ) f( )) + k + k k k he first term is always non-positive due to our assumptions on f. he second, third and forth terms are non-positive as well due to the arithmetic-geometric mean (AM- GM) inequality. Hence, Ẇ 0in R 3 +. o characterize the subset of R 3 + where Ẇ equals zero, we distinguish two cases: k 4 > 0 and k 4 = 0. If k 4 > 0, then Ẇ equals zero if and only if (,,)=E. Since W is a proper Lyapunov function in R 3 +, it follows that E is globally asymptotically stable in R 3 +.Ifk 4 = 0, then Ẇ =0if and only if = and =. Since all solutions of () in R3 + are bounded by Lemma 3., the LaSalle s invariance principle implies that any ω-limit set in R 3 + is a subset of the largest invariant set in M = {(,,) R 3 + =, = }. Any such invariant set in M must satisfy = 0, hence 0=f( ) k + k2 = f( ), which implies that = and =. herefore, the largest invariant set in M is the singleton {E}, and since W is a proper function in R 3 +, E attracts all solutions in R Competitive exclusion. he global stability result of heorem 5. can be readily extended to the model with several viral strains competing for the same pool of susceptible cells. he multistrain competition model takes the form = f( ) i ki i i ki 2i, i = k i i + k2 i i β ii, i =,...,n, (22) i = N i β i i i i, i =,...,n, on R 2n+ +, where we neglect the loss of virions due to cellular absorption. For each strain, we define the quantities i 0 = ki 2 0 β i + N ik i 0 i, i = We have the following extinction result. i β i k2 i i + N i k i β = 0 i i. heorem 5.2. Suppose that 0 > max{, 2 0,..., n 0 }. hen the equilibrium E = (,,, 0,...,0) attracts all solutions of (22) with (0) + (0) > 0. Here, := f( )/β and := N f( )/. Proof. We begin by observing that 0 > 0 i is equivalent to < i for i =2,...,n. It is also easy to check that E is indeed an equilibrium of (22). Define W as follows W = τ dτ + + i 2( i + B i i ), τ dτ + k 0 τ dτ
13 WIHIN HOS IRUS MODELS 3 where B i = ki i for i =2,...,n.Evaluating Ẇ, we find that Ẇ = A +A A n, where A = ( )(f( ) f( )) +k2 2 + k 3 is the same as Ẇ in the proof of heorem 5., and A i = (k i i + k i 2 i ) β i i + B i (N i β i i i ) = i ( k2 i β i + B i N i β i )+ i ( k i B i i ) = i k2 i β i + ki N i β i i = ki 2 i + N i k i β i i Clearly, all A i s are non-positive, hence Ẇ = 0 is given by i ( i )= β i i 0 0 i ( i ). Ẇ 0. Furthermore, the set M where M = { =, =, 2 =... = n =0}, thus the largest invariant set of (22) inm is E. An application of the LaSalle s invariance principle concludes the proof. 6. he case 0 <k 3 << k 4. In this section, we use a global perturbation result [20] to show that the positive equilibrium E is globally asymptotically stable for system () in R 3 +, provided that k i > 0 for i =,...,4, and k 3 is sufficiently small. We begin with the following result. Lemma 6.. Suppose that k,k 2,k 4 > 0, k 3 =0,andR 0 >. Let W be the Lyapunov function defined in the proof of heorem 5.. hen there exist, 2 > 0 such that Ẇ (x) for all x R 3 + such that x E 2. In other words, Ẇ is bounded away from zero outside of any compact neighborhood of E in R 3 +. Proof. From the proof of heorem 5., we have that Ẇ k 3 + k k 4 o simplify the notation, let 2 + k 2. 2 A = k, A2 = k 2, A 3 = k k 4, u =, v =, w =, so that Ẇ A 3 u u v v A 2 (u ) 2 u A 3 v(w ) 2. For each v>0, the first term is maximized when u = v,hence 3 u u v v 3 2 v v = 2+ v v ( v ) 2 = 2+ v v( v + ) 2 (v )2.
14 4 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER herefore, 2+ v Ẇ A (v (u ) 2 v( v + ) 2 )2 A 2 A 3 v(w ) 2. u Now, let >0 be sufficiently small, and suppose that Ẇ ( 2, 0). hen we have (u ) 2 A 2 < 2 u <κ u <κ, for some κ > 0 which is independent of. Similarly, A 2+ v (v v( v + ) 2 )2 < 2 v <κ 2, for some κ 2 > 0 which is independent of. his implies that w A 3 ( κ 2 ) <κ 3 <κ 3, for an appropriate κ 3 > 0. Lastly, we have that = w w v (κ 2 + κ 3 ) <κ 4, κ 3 for some κ 4 > 0. Combining the above inequalities, we conclude that there exists κ 0 > 0 such that for all sufficiently small >0, the inequality 2 < Ẇ (x) < 0 implies that x E <κ 0, and the claim follows. heorem 6.2. Suppose that k,k 2,k 4 > 0, k 3 =0,andR 0 >. hen there exist δ>0 and a smooth mapping E :[0,δ) R 3 +, such that for every k 3 [0,δ) the point E(k 3 ) is a globally asymptotically stable equilibrium of () in R 3 +. Proof. First, we observe that R 0 is a continuous function of k 3 and for k 3 = 0, R 0 >. Hence, there exists δ > 0 such that R 0 > as long as k 3 [0,δ ). By Lemma 2.2, for each k 3 [0,δ ) there exists a unique positive equilibrium E(k 3 ) of () which is locally asymptotically stable by heorem 3.3. By the Implicit Function heorem, the mapping E :[0,δ ) R 3 + is smooth. Now, let, 2 > 0 be as defined in Lemma 6.. Let W be the Lyapunov function defined in the proof of heorem 5.. SinceW is a proper Lyapunov function, there exists a sufficiently small 0 > 0 such that W (x) 0 implies x E(0) 2.Define the compact set K = W ([0, 0 ]). Let δ 2 := 2k (0 + ), and suppose that k 3 [0,δ 2 ). hen for any solution of () in R 3 + it holds that (t) < 0 + for all sufficiently large t. We have that Ẇ = (f( ) f( )) + k 3 + k k k 4 2 k k 3 ( ). As long as x K and (t) < 0 +, Lemma 6. implies that Ẇ + k δ 2 (0 + ) 2 < 0, hence every positive solution enters the set K in finite time, so that K is an absorbing forward invariant compact set for () for all k 3 [0,δ 2 ). Now, letting
15 WIHIN HOS IRUS MODELS Log0 3 2 Uninfected cells Infected cells irus ime Figure. ime series of a solution ( (t), (t),(t)) of system () with parameter values as described in Example in the text. δ =min(δ,δ 2 ), we apply the global perturbation result (Proposition 2.3 in [20]) to conclude that E(k 3 ) is globally stable under () for all k 3 [0,δ). 7. Numerical examples. In this Section we numerically explore the behavior of our model when some of the conditions of the main results, heorem 4.4 and heorem 6.2, fail. In both cases, it appears that the disease steady state attracts the simulated solution, providing evidence for the conjecture that this steady state is globally stable when it exists. Example : In this example, we set k 4 = 0 and we have chosen the parameters such that both 2k 2 0 >βand k 3 (f(0) + β 0 ) >β, which implies that heorem 4.4 is not applicable. We used f( )=a b with a = 0 4 ml day and b =0.0day (which implies that 0 = 0 6 ml ), k = ml day,k 2 = 0 6 ml day,k 3 = 0 4 ml day,n = 3000, = 23day,β =day. Notice that with these assignments we have that R 0 =.44 >, 2k 2 0 β = > 0, and k 3 (f(0) + β 0 ) β = 78 > 0. he time series of the components of a simulated solution with initial condition ( (0), (0),(0)) = (0 6, 0, 0 6 ) are displayed in Figure. hey suggest convergence to a positive steady state. Example 2: In this example, we have chosen k 3 >> k 4 > 0 (in fact, k 3 = 0k 4 ), which implies that heorem 6.2 may not be applicable. Once again, we have set f( )=a b, and for a, b, k,k 2,k 3,N, and β we used the same values as in Example. On the other hand, here, k 4 = 0 5 ml day. Notice that these assignments yield the same value for R 0 =.44 as in Example because this value is independent of k 4 ;see(3). he time series of the components of a simulated solution with the same initial condition as in Example are displayed in Figure 2, and they suggest convergence to a positive steady state as before. 8. Conclusions. In this paper, we have analyzed a within-host virus model which incorporated the mechanisms of direct cell-to-cell viral transmission and the viral coinfection (absorption of free virions into already infected cells). We obtained a complete analytic description of equilibria and their local stability. We presented the basic reproductive number R 0 and proved that the infection persists when R 0 > and becomes extinct when R 0. We also obtained sufficient conditions for the global stability of the endemic equilibrium in several particular cases. We were unable to establish the global stability in the most general cases where all coefficients
16 6 H. POURBASHASH, S.S. PILYUGIN, C. MCCLUSKEY AND P. DE LEENHEER Log Uninfected cells Infected cells irus ime Figure 2. ime series of a solution ( (t), (t),(t)) of system () with parameter values as described in Example 2 in the text. are positive. In one instance, namely when the absorption of free virions into infected cells can be neglected (the case k 3 = 0), we used the geometric approach developed by Li and Muldowney to prove the global stability [0]. In another instance, when the absorption of free virions into uninfected cells can be neglected (admittedly, not a biologically plausible scenario), we employed a Lyapunov function approach with a specific function similar to that used by Li et al in [9]. If both k 3 = k 4 = 0, we extended the global stability result to include a competition of several viral strains, and proved that the competitive exclusion is the typical outcome. Moreover, as it is common in competition models, it turns out that the fittest competitor is determined by the highest value of the basic reproduction number, R 0. Another implication of this result is that a direct cell-to-cell viral transmission cannot induce coexistence of viral strains. Using the same Lyapunov function, we we were able to invoke a global perturbation result of Smith and Waltman [20] to prove the global stability of the endemic equilibrium when k 3 > 0 is a small parameter. Finally, we presented two numerical examples to illustrate the global stability of the endemic equilibrium. We conjecture that such equilibrium is globally stable whenever it exists, but the proof of this general result remains an open question. REFERENCES [] W.A. Coppel, Stability and asymptotic behavior of differential equations, Heath and Co., Boston, 965. [2] P. De Leenheer and H.L. Smith, irus dynamics: A global analysis, SIAM J. Appl. Math. 63 (2003), [3] P. De Leenheer and S.S. Pilyugin, Multistrain virus dynamics with mutations: a global analysis, Math.Med.Biol.25 (2008), [4] N. Dixit and A. Perelson, Multiplicity of human immunodeficiency virus infections in lymphoid tissue, J.irol.78 (2004), [5] M. Fiedler, Additive compound matrices and inequality for eigenvalues of stochastic matrices, Czech. Math. J. (974), [6] H.I. Freedman, M.X. ang and S.G. Ruan, Uniform persistence and flows near a closed positively invariant set, J.Dynam.DifferentialEquations6 (994), [7] H.K. Khalil, Nonlinear Systems, 3rd Edition, Prentice Hall, [8] A. Korobeinikov, Lyapunov functions and global properties for SEIR and SEIS epidemic models, Math.Med.Biol.2 (2005), [9] M.Y. Li, J.R. Graef, L. Wang and J. Karsai, Global dynamics of a SEIR model with varying total population size, Math.Biosci.60 (999), [0] M.Y. Li and J.S. Muldowney, A geometric approach to the global-stability problems, SIAMJ. Math. Anal. 27 (996),
17 WIHIN HOS IRUS MODELS 7 [] M.Y. Li and J.S. Muldowney, Global stability for the SEIR model in epidemiology, Math. Biosci. 25 (995), [2] R.H. Jr. Martin, Logarithmic norms and projections applied to linear differential systems, J. Math. Anal. Appl. 45 (974), [3] D. Mazurov, A. Ilinskaya, G.Heidecker, P. Lloyd and D. Derse, Quantitative Comparison of HL- and HI- Cell-to- Cell Infection with New Replication Dependent ectors, PLoS Pathogens 6 (200) e [4] B. Monel, E. Beaumont, D. endrame, O. Schwartz, D. Brand and F. Mammano, HI cellto-cell transmission requires the production of infectious virus particles and does not proceed through Env-mediated fusion pores, J.irol.86 (202), [5] J.S. Muldowney, Compound matrices and ordinary differential equations, RockyMount.J. Math. 20 (990), [6] M.A. Nowak and R.M. May, irus dynamics, Oxford University press, New York, [7] A.S. Perelson and P.W. Nelson, Mathematical analysis of HI- dynamics in vivo, SIAM Rev. 4 (999), [8]. Piguet and Q. Sattentau, Dangerous liaisons at the virological synapse, J.Clin.Invest. 4 (2004), [9] O. Schwartz, Immunological and virological aspects of HI cell-to-cell transfer, Retrovirology 6 (2009), I6. [20] H.L. Smith and P. Waltman, Perturbation of a globally stable steady state, Proc.Am.Math. Soc. 27 (999) [2] M. Sourisseau, N. Sol-Foulon, F. Porrot, F. Blanchet and O. Schwartz, Inefficient human immunodeficiency virus replication in mobile lymphocytes, J.irol.8 (2007), [22] P. van den Driessche and J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math.Biosci.80 (2002), [23] L. Wang and S. Ellermeyer, HI infection and CD4 + celldynamics,discretecontin.dyn. Syst. Ser. B 6 (2006), [24] L. Wang and M.Y. Li, Mathematical analysis of the global dynamics of a model for HI infection of CD4 + cells,math.biosci.200 (2006), Received xxxx 20xx; revised xxxx 20xx. address: pourbashash@ufl.edu address: pilyugin@ufl.edu address: ccmcc8@gmail.com address: deleenhe@ufl.edu
Global Stability of SEIRS Models in Epidemiology
Global Stability of SRS Models in pidemiology M. Y. Li, J. S. Muldowney, and P. van den Driessche Department of Mathematics and Statistics Mississippi State University, Mississippi State, MS 39762 Department
More informationGLOBAL STABILITY OF SIR MODELS WITH NONLINEAR INCIDENCE AND DISCONTINUOUS TREATMENT
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 304, pp. 1 8. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY
More informationResearch Article Modeling Computer Virus and Its Dynamics
Mathematical Problems in Engineering Volume 213, Article ID 842614, 5 pages http://dx.doi.org/1.1155/213/842614 Research Article Modeling Computer Virus and Its Dynamics Mei Peng, 1 Xing He, 2 Junjian
More informationMathematical models on Malaria with multiple strains of pathogens
Mathematical models on Malaria with multiple strains of pathogens Yanyu Xiao Department of Mathematics University of Miami CTW: From Within Host Dynamics to the Epidemiology of Infectious Disease MBI,
More informationGlobal stability of an SEIS epidemic model with recruitment and a varying total population size
Mathematical Biosciences 170 2001) 199±208 www.elsevier.com/locate/mbs Global stability of an SEIS epidemic model with recruitment and a varying total population size Meng Fan a, Michael Y. Li b, *, Ke
More informationA Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and Cure of Infected Cells in Eclipse Stage
Applied Mathematical Sciences, Vol. 1, 216, no. 43, 2121-213 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.216.63128 A Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and
More informationThursday. Threshold and Sensitivity Analysis
Thursday Threshold and Sensitivity Analysis SIR Model without Demography ds dt di dt dr dt = βsi (2.1) = βsi γi (2.2) = γi (2.3) With initial conditions S(0) > 0, I(0) > 0, and R(0) = 0. This model can
More informationGLOBAL DYNAMICS OF A MATHEMATICAL MODEL OF TUBERCULOSIS
CANADIAN APPIED MATHEMATICS QUARTERY Volume 13, Number 4, Winter 2005 GOBA DYNAMICS OF A MATHEMATICA MODE OF TUBERCUOSIS HONGBIN GUO ABSTRACT. Mathematical analysis is carried out for a mathematical model
More informationGlobal Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response
Global Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response Gang Huang 1,2, Wanbiao Ma 2, Yasuhiro Takeuchi 1 1,Graduate School of Science and Technology, Shizuoka University,
More informationFeedback control for a chemostat with two organisms
Feedback control for a chemostat with two organisms Patrick De Leenheer and Hal Smith Arizona State University Department of Mathematics and Statistics Tempe, AZ 85287 email: leenheer@math.la.asu.edu,
More informationFeedback-mediated oscillatory coexistence in the chemostat
Feedback-mediated oscillatory coexistence in the chemostat Patrick De Leenheer and Sergei S. Pilyugin Department of Mathematics, University of Florida deleenhe,pilyugin@math.ufl.edu 1 Introduction We study
More informationGlobal Qualitative Analysis for a Ratio-Dependent Predator Prey Model with Delay 1
Journal of Mathematical Analysis and Applications 266, 401 419 (2002 doi:10.1006/jmaa.2001.7751, available online at http://www.idealibrary.com on Global Qualitative Analysis for a Ratio-Dependent Predator
More informationGlobal Stability of a Computer Virus Model with Cure and Vertical Transmission
International Journal of Research Studies in Computer Science and Engineering (IJRSCSE) Volume 3, Issue 1, January 016, PP 16-4 ISSN 349-4840 (Print) & ISSN 349-4859 (Online) www.arcjournals.org Global
More informationGLOBAL STABILITY OF A VACCINATION MODEL WITH IMMIGRATION
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 92, pp. 1 10. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY
More informationA REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 A REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES KING-YEUNG LAM
More informationStability of SEIR Model of Infectious Diseases with Human Immunity
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 6 (2017), pp. 1811 1819 Research India Publications http://www.ripublication.com/gjpam.htm Stability of SEIR Model of Infectious
More informationMathematical Analysis of Epidemiological Models III
Intro Computing R Complex models Mathematical Analysis of Epidemiological Models III Jan Medlock Clemson University Department of Mathematical Sciences 27 July 29 Intro Computing R Complex models What
More informationarxiv: v2 [q-bio.pe] 3 Oct 2018
Journal of Mathematical Biology manuscript No. (will be inserted by the editor Global stability properties of renewal epidemic models Michael T. Meehan Daniel G. Cocks Johannes Müller Emma S. McBryde arxiv:177.3489v2
More informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationCompetitive Exclusion in a Discrete-time, Size-structured Chemostat Model
Competitive Exclusion in a Discrete-time, Size-structured Chemostat Model Hal L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 1804, USA E-mail: halsmith@asu.edu Xiao-Qiang Zhao
More informationViral evolution model with several time scales
Viral evolution model with several time scales AA Archibasov Samara National Research University, 4 Moskovskoe Shosse, 4486, Samara, Russia Abstract In this paper a viral evolution model with specific
More informationGlobal Analysis of an SEIRS Model with Saturating Contact Rate 1
Applied Mathematical Sciences, Vol. 6, 2012, no. 80, 3991-4003 Global Analysis of an SEIRS Model with Saturating Contact Rate 1 Shulin Sun a, Cuihua Guo b, and Chengmin Li a a School of Mathematics and
More informationHETEROGENEOUS MIXING IN EPIDEMIC MODELS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number 1, Spring 212 HETEROGENEOUS MIXING IN EPIDEMIC MODELS FRED BRAUER ABSTRACT. We extend the relation between the basic reproduction number and the
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationAustralian Journal of Basic and Applied Sciences. Effect of Personal Hygiene Campaign on the Transmission Model of Hepatitis A
Australian Journal of Basic and Applied Sciences, 9(13) Special 15, Pages: 67-73 ISSN:1991-8178 Australian Journal of Basic and Applied Sciences Journal home page: wwwajbaswebcom Effect of Personal Hygiene
More informationThe effect of population dispersal on the spread of a disease
J. Math. Anal. Appl. 308 (2005) 343 364 www.elsevier.com/locate/jmaa The effect of population dispersal on the spread of a disease Yu Jin, Wendi Wang Department of Mathematics, Southwest China Normal University,
More informationEpidemics in Two Competing Species
Epidemics in Two Competing Species Litao Han 1 School of Information, Renmin University of China, Beijing, 100872 P. R. China Andrea Pugliese 2 Department of Mathematics, University of Trento, Trento,
More informationarxiv: v2 [q-bio.pe] 22 Feb 2018
Noname manuscript No. (will be inserted by the editor) Dynamics of Virus and Immune Response in Multi-Epitope Network Cameron J. Browne Hal L. Smith arxiv:75.88v2 [q-bio.pe] 22 Feb 28 the date of receipt
More informationThe death of an epidemic
LECTURE 2 Equilibrium Stability Analysis & Next Generation Method The death of an epidemic In SIR equations, let s divide equation for dx/dt by dz/ dt:!! dx/dz = - (β X Y/N)/(γY)!!! = - R 0 X/N Integrate
More informationSIR Epidemic Model with total Population size
Advances in Applied Mathematical Biosciences. ISSN 2248-9983 Volume 7, Number 1 (2016), pp. 33-39 International Research Publication House http://www.irphouse.com SIR Epidemic Model with total Population
More informationDynamics of Disease Spread. in a Predator-Prey System
Advanced Studies in Biology, vol. 6, 2014, no. 4, 169-179 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/asb.2014.4845 Dynamics of Disease Spread in a Predator-Prey System Asrul Sani 1, Edi Cahyono
More informationMathematical Analysis of Epidemiological Models: Introduction
Mathematical Analysis of Epidemiological Models: Introduction Jan Medlock Clemson University Department of Mathematical Sciences 8 February 2010 1. Introduction. The effectiveness of improved sanitation,
More informationDelay SIR Model with Nonlinear Incident Rate and Varying Total Population
Delay SIR Model with Nonlinear Incident Rate Varying Total Population Rujira Ouncharoen, Salinthip Daengkongkho, Thongchai Dumrongpokaphan, Yongwimon Lenbury Abstract Recently, models describing the behavior
More informationStability Analysis of a SIS Epidemic Model with Standard Incidence
tability Analysis of a I Epidemic Model with tandard Incidence Cruz Vargas-De-León Received 19 April 2011; Accepted 19 Octuber 2011 leoncruz82@yahoo.com.mx Abstract In this paper, we study the global properties
More informationPreprint núm May Stability of a stochastic model for HIV-1 dynamics within a host. L. Shaikhet, A. Korobeinikov
Preprint núm. 91 May 2014 Stability of a stochastic model for HIV-1 dynamics within a host L. Shaikhet, A. Korobeinikov STABILITY OF A STOCHASTIC MODEL FOR HIV-1 DYNAMICS WITHIN A HOST LEONID SHAIKHET
More informationGlobal attractors, stability, and population persistence
Global attractors, stability, and population persistence Horst R. Thieme (partly joint work with Hal L. Smith*) School of Mathematical and Statistical Sciences Arizona State University partially supported
More informationGlobal stability for a four dimensional epidemic model
Note di Matematica SSN 1123-2536, e-ssn 1590-0932 Note Mat. 30 (2010) no. 2, 83 95. doi:10.1285/i15900932v30n2p83 Global stability for a four dimensional epidemic model Bruno Buonomo Department of Mathematics
More informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
More informationSimple Mathematical Model for Malaria Transmission
Journal of Advances in Mathematics and Computer Science 25(6): 1-24, 217; Article no.jamcs.37843 ISSN: 2456-9968 (Past name: British Journal of Mathematics & Computer Science, Past ISSN: 2231-851) Simple
More informationSensitivity and Stability Analysis of Hepatitis B Virus Model with Non-Cytolytic Cure Process and Logistic Hepatocyte Growth
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 3 2016), pp. 2297 2312 Research India Publications http://www.ripublication.com/gjpam.htm Sensitivity and Stability Analysis
More informationSTUDY OF THE DYNAMICAL MODEL OF HIV
STUDY OF THE DYNAMICAL MODEL OF HIV M.A. Lapshova, E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper is devoted to the study of the dynamical model of HIV. An application
More informationApplied Mathematics Letters
Applied athematics Letters 25 (212) 156 16 Contents lists available at SciVerse ScienceDirect Applied athematics Letters journal homepage: www.elsevier.com/locate/aml Globally stable endemicity for infectious
More informationDYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS. Wei Feng
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 7 pp. 36 37 DYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS Wei Feng Mathematics and Statistics Department
More informationAnalysis of SIR Mathematical Model for Malaria disease with the inclusion of Infected Immigrants
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 14, Issue 5 Ver. I (Sep - Oct 218), PP 1-21 www.iosrjournals.org Analysis of SIR Mathematical Model for Malaria disease
More information6. Age structure. for a, t IR +, subject to the boundary condition. (6.3) p(0; t) = and to the initial condition
6. Age structure In this section we introduce a dependence of the force of infection upon the chronological age of individuals participating in the epidemic. Age has been recognized as an important factor
More informationThe dynamics of disease transmission in a Prey Predator System with harvesting of prey
ISSN: 78 Volume, Issue, April The dynamics of disease transmission in a Prey Predator System with harvesting of prey, Kul Bhushan Agnihotri* Department of Applied Sciences and Humanties Shaheed Bhagat
More informationGlobal compact attractors and their tripartition under persistence
Global compact attractors and their tripartition under persistence Horst R. Thieme (joint work with Hal L. Smith) School of Mathematical and Statistical Science Arizona State University GCOE, September
More informationPermanence Implies the Existence of Interior Periodic Solutions for FDEs
International Journal of Qualitative Theory of Differential Equations and Applications Vol. 2, No. 1 (2008), pp. 125 137 Permanence Implies the Existence of Interior Periodic Solutions for FDEs Xiao-Qiang
More informationCan multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species
Can multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species Xingfu Zou Department of Applied Mathematics University of Western Ontario London, Ontario, Canada (Joint work
More informationSpectrally arbitrary star sign patterns
Linear Algebra and its Applications 400 (2005) 99 119 wwwelseviercom/locate/laa Spectrally arbitrary star sign patterns G MacGillivray, RM Tifenbach, P van den Driessche Department of Mathematics and Statistics,
More informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationGLOBAL STABILITY OF THE ENDEMIC EQUILIBRIUM OF A TUBERCULOSIS MODEL WITH IMMIGRATION AND TREATMENT
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 19, Number 1, Spring 2011 GLOBAL STABILITY OF THE ENDEMIC EQUILIBRIUM OF A TUBERCULOSIS MODEL WITH IMMIGRATION AND TREATMENT HONGBIN GUO AND MICHAEL Y. LI
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationStability Analysis of an HIV/AIDS Epidemic Model with Screening
International Mathematical Forum, Vol. 6, 11, no. 66, 351-373 Stability Analysis of an HIV/AIDS Epidemic Model with Screening Sarah Al-Sheikh Department of Mathematics King Abdulaziz University Jeddah,
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More informationMathematical Modeling and Analysis of Infectious Disease Dynamics
Mathematical Modeling and Analysis of Infectious Disease Dynamics V. A. Bokil Department of Mathematics Oregon State University Corvallis, OR MTH 323: Mathematical Modeling May 22, 2017 V. A. Bokil (OSU-Math)
More informationHAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
Georgian Mathematical Journal Volume 9 (2002), Number 3, 591 600 NONEXPANSIVE MAPPINGS AND ITERATIVE METHODS IN UNIFORMLY CONVEX BANACH SPACES HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
More informationDynamical Systems in Biology
Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems
More informationGlobal Analysis of an Epidemic Model with Nonmonotone Incidence Rate
Global Analysis of an Epidemic Model with Nonmonotone Incidence Rate Dongmei Xiao Department of Mathematics, Shanghai Jiaotong University, Shanghai 00030, China E-mail: xiaodm@sjtu.edu.cn and Shigui Ruan
More informationProofs for Large Sample Properties of Generalized Method of Moments Estimators
Proofs for Large Sample Properties of Generalized Method of Moments Estimators Lars Peter Hansen University of Chicago March 8, 2012 1 Introduction Econometrica did not publish many of the proofs in my
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationAustralian Journal of Basic and Applied Sciences. Effect of Education Campaign on Transmission Model of Conjunctivitis
ISSN:99-878 Australian Journal of Basic and Applied Sciences Journal home page: www.ajbasweb.com ffect of ducation Campaign on Transmission Model of Conjunctivitis Suratchata Sangthongjeen, Anake Sudchumnong
More informationStability and bifurcation analysis of epidemic models with saturated incidence rates: an application to a nonmonotone incidence rate
Published in Mathematical Biosciences and Engineering 4 785-85 DOI:.3934/mbe.4..785 Stability and bifurcation analysis of epidemic models with saturated incidence rates: an application to a nonmonotone
More informationESTIMATES ON THE NUMBER OF LIMIT CYCLES OF A GENERALIZED ABEL EQUATION
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 ESTIMATES ON THE NUMBER OF LIMIT CYCLES OF A GENERALIZED ABEL EQUATION NAEEM M.H. ALKOUMI
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationApparent paradoxes in disease models with horizontal and vertical transmission
Apparent paradoxes in disease models with horizontal and vertical transmission Thanate Dhirasakdanon, Stanley H.Faeth, Karl P. Hadeler*, Horst R. Thieme School of Life Sciences School of Mathematical and
More informationSTABILITY ANALYSIS OF A GENERAL SIR EPIDEMIC MODEL
VFAST Transactions on Mathematics http://vfast.org/index.php/vtm@ 2013 ISSN: 2309-0022 Volume 1, Number 1, May-June, 2013 pp. 16 20 STABILITY ANALYSIS OF A GENERAL SIR EPIDEMIC MODEL Roman Ullah 1, Gul
More informationQualitative Analysis of a Discrete SIR Epidemic Model
ISSN (e): 2250 3005 Volume, 05 Issue, 03 March 2015 International Journal of Computational Engineering Research (IJCER) Qualitative Analysis of a Discrete SIR Epidemic Model A. George Maria Selvam 1, D.
More informationWeighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 401 412 (2013) http://campus.mst.edu/adsa Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes
More informationDynamics of a Hepatitis B Viral Infection Model with Logistic Hepatocyte Growth and Cytotoxic T-Lymphocyte Response
Nonlinear Analysis and Differential Equations, Vol. 4, 16, no. 3, 19-1 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/nade.16.5164 Dynamics of a Hepatitis B Viral Infection Model with Logistic Hepatocyte
More informationNOTES ON LINEAR ODES
NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary
More informationSUBTHRESHOLD AND SUPERTHRESHOLD COEXISTENCE OF PATHOGEN VARIANTS: THE IMPACT OF HOST AGE-STRUCTURE
SUBTHRESHOLD AND SUPERTHRESHOLD COEXISTENCE OF PATHOGEN VARIANTS: THE IMPACT OF HOST AGE-STRUCTURE MAIA MARTCHEVA, SERGEI S. PILYUGIN, AND ROBERT D. HOLT Abstract. It is well known that in the most general
More informationGlobal Stability Analysis on a Predator-Prey Model with Omnivores
Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani
More informationOn Stability of Linear Complementarity Systems
On Stability of Linear Complementarity Systems Alen Agheksanterian University of Maryland Baltimore County Math 710A: Special Topics in Optimization On Stability of Linear Complementarity Systems Alen
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationarxiv: v1 [math.ds] 14 Dec 2017
Global properties of an age-structured virus model with saturated antibody immune response, multi-target cells and general incidence rate Ángel Cervantes-Pérez, Eric Ávila-Vales arxiv:1712.564v1 [math.ds]
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationA Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems
53rd IEEE Conference on Decision and Control December 15-17, 2014. Los Angeles, California, USA A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems Seyed Hossein Mousavi 1,
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a ournal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationIntroduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium
Introduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium point and the model reproduction number Prove the stability
More informationATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.
ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS Emerson A. M. de Abreu Alexandre N. Carvalho Abstract Under fairly general conditions one can prove that
More informationDynamics of a Networked Connectivity Model of Waterborne Disease Epidemics
Dynamics of a Networked Connectivity Model of Waterborne Disease Epidemics A. Edwards, D. Mercadante, C. Retamoza REU Final Presentation July 31, 214 Overview Background on waterborne diseases Introduction
More informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationCompetitive and Cooperative Differential Equations
Chapter 3 Competitive and Cooperative Differential Equations 0. Introduction This chapter and the next one focus on ordinary differential equations in IR n. A natural partial ordering on IR n is generated
More informationStochastic Viral Dynamics with Beddington-DeAngelis Functional Response
Stochastic Viral Dynamics with Beddington-DeAngelis Functional Response Junyi Tu, Yuncheng You University of South Florida, USA you@mail.usf.edu IMA Workshop in Memory of George R. Sell June 016 Outline
More informationCHAPTER 10 Shape Preserving Properties of B-splines
CHAPTER 10 Shape Preserving Properties of B-splines In earlier chapters we have seen a number of examples of the close relationship between a spline function and its B-spline coefficients This is especially
More informationGLOBAL ASYMPTOTIC STABILITY OF A DIFFUSIVE SVIR EPIDEMIC MODEL WITH IMMIGRATION OF INDIVIDUALS
Electronic Journal of Differential Equations, Vol. 16 (16), No. 8, pp. 1 1. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu GLOBAL ASYMPTOTIC STABILITY OF A DIFFUSIVE SVIR
More informationOn the Dynamics of a Rational Difference Equation, Part 1
International Journal of Difference Equations (IJDE). ISSN 0973-6069 Volume 3 Number 1 (2008), pp. 1 35 Research India Publications http://www.ripublication.com/ijde.htm On the Dynamics of a Rational Difference
More informationGlobal Analysis for an HIV Infection Model with CTL Immune Response and Infected Cells in Eclipse Phase
Article Global Analysis for an HIV Infection Model with CTL Immune Response and Infected Cells in Eclipse Phase Karam Allali 1 ID Jaouad Danane 1 ID and Yang Kuang 2 * ID 1 Laboratory of Mathematics and
More informationAARMS Homework Exercises
1 For the gamma distribution, AARMS Homework Exercises (a) Show that the mgf is M(t) = (1 βt) α for t < 1/β (b) Use the mgf to find the mean and variance of the gamma distribution 2 A well-known inequality
More informationAsymptotic of Enumerative Invariants in CP 2
Peking Mathematical Journal https://doi.org/.7/s4543-8-4-4 ORIGINAL ARTICLE Asymptotic of Enumerative Invariants in CP Gang Tian Dongyi Wei Received: 8 March 8 / Revised: 3 July 8 / Accepted: 5 July 8
More informationDiscrete Time Markov Chain of a Dynamical System with a Rest Phase
Discrete Time Markov Chain of a Dynamical System with a Rest Phase Abstract A stochastic model, in the form of a discrete time Markov chain, is constructed to describe the dynamics of a population that
More information