Global stability of delay-distributed viral infection model with two modes of viral transmission and B-cell impairment

This paper formulates a virus dynamics model with impairment of B-cell functions. The model incorporates two modes of viral transmission: cell-free and cell-to-cell. The cell-free and cell-cell incidence rates are modeled by general functions. The model incorporates both, latently and actively, infected cells as well as three distributed time delays. Nonnegativity and boundedness properties of the solutions are proven to show the well-posedness of the model. The model admits two equilibria that are determined by the basic reproduction number R₀. The global stability of each equilibrium is proven by utilizing Lyapunov function and LaSalle's invariance principle. The theoretical results are illustrated by numerical simulations. The effect of impairment of B-cell functions and time delays on the virus dynamics are studied. We have shown that if the functions of B-cell is impaired, then the concentration of viruses is increased in the plasma. Moreover, we have observed that increasing the time delay will suppress the viral replication.     

Stability of a general delay‐distributed virus dynamics model with multi‐staged infected progression and immune response

In this paper, we formulate a (n + 3)‐dimensional nonlinear virus dynamics model that considers n‐stages of the infected cells and n + 1 distributed time delays. The model incorporates humoral immune response and general nonlinear forms for the incidence rate of infection, the generation and removal rates of the cells and viruses. Under a set of conditions on the general functions, the basic infection reproduction number and the humoral immune response activation number are derived. Utilizing Lyapunov functionals and LaSalle's invariance principle, the global asymptotic stability of all steady states of the model are proven. Numerical simulations are carried out to confirm the theoretical results. Copyright © 2016 John Wiley & Sons, Ltd.     

Global stability of humoral immunity virus dynamics models with nonlinear infection rate and removal

In this paper, we investigate the dynamical behavior of two nonlinear models for viral infection with humoral immune response. The first model contains four compartments; uninfected target cells, actively infected cells, free virus particles and B cells. The intrinsic growth rate of uninfected cells, incidence rate of infection, removal rate of infected cells, production rate of viruses, neutralization rate of viruses, activation rate of B cells and removal rate of B cells are given by more general nonlinear functions. The second model is a modification of the first one by including an eclipse stage of infected cells. We assume that the latent-to-active conversion rate is also given by a more general nonlinear function. For each model we derive two threshold parameters and establish a set of conditions on the general functions which are sufficient to determine the global dynamics of the models. By using suitable Lyapunov functions and LaSalle’s invariance principle, we prove the global asymptotic stability of the all equilibria of the models. We perform some numerical simulations for the models with specific forms of the general functions and show that the numerical results are consistent with the theoretical results.     

Stability of general virus dynamics models with both cellular and viral infections and delays

We consider general virus dynamics model with virus‐to‐target and infected‐to‐target infections. The model is incorporated by intracellular discrete or distributed time delays. We assume that the virus‐target and infected‐target incidences, the production, and clearance rates of all compartments are modeled by general nonlinear functions that satisfy a set of reasonable conditions. The non‐negativity and boundedness of the solutions are studied. The existence and stability of the equilibria are determined by a threshold parameter. We use suitable Lyapunov functionals and apply LaSalle's invariance principle to prove the global asymptotic stability of the all equilibria of the model. We confirm the theoretical results by numerical simulations. Copyright © 2017 John Wiley & Sons, Ltd.     

NONLINEAR MATRIX MODELS AND POPULATION DYNAMICS

Nonlinear matrix difference equations are studied as models for the discrete time dynamics of a population whose individual members have been categorized into a finite number of classes. The equations are treated with sufficient generality so as to include virtually any type of structuring of the population (the sole constraint is that all newborns lie in the same class) and any types of nonlinearities which arise from the density dependence of fertility rates, survival rates and transition probabilities between classes. The existence and stability of equilibrium class distribution vectors are studied by means of bifurcation theory techniques using a single composite, biologically meaningful quantity as a bifurcation parameter, namely the inherent net reproductive rate r. It is shown that, just as in the case of linear matrix equations, a global continuum of positive equilibria exists which bifurcates as a function of r from the zero equilibrium state at and only at r = 1. Furthermore the zero equilibrium loses stability as r is increased through 1. Unlike the linear case however, for which the bifurcation is “vertical” (i.e., equilibria exist only for r = 1), the nonlinear equation in general has positive equilibria for an interval of r values. Methods for studying the geometry of the continuum based upon the density dependence of the net reproductive rate at equilibrium are developed. With regard to stability, it is shown that in general the positive equilibria near the bifurcation point are stable if the bifurcation is to the right and unstable if it is to the left. Some further results and conjectures concerning stability are also given. The methods are illustrated by several examples involving nonlinear models of various types taken from the literature.     

Stability of HIV/HTLV-I co-infection model with delays

In this paper, we formulate a within-host dynamics model for HIV/HTLV-I co-infection under the influence of cytotoxic T lymphocytes (CTLs). The model incorporates silent HIV-infected CD4⁺T cells and silent HTLV-infected CD4⁺T cells. The model includes two routes of HIV transmission, virus to cell (VTC) and cell to cell (CTC). It also incorporates two modes of HTLV-I transmission, horizontal transmission via direct CTC contact and vertical transmission through mitotic division of Tax-expressing HTLV-infected cells. The model takes into account five types of distributed-time delays. We analyze the model by proving the nonnegativity and boundedness of the solutions, calculating all possible equilibria, deriving a set of key threshold parameters, and proving the global stability of all equilibria. The global asymptotic stability of all equilibria is established by utilizing Lyapunov function and LaSalle's invariance principle. We present numerical simulations to justify the applicability and effectiveness of the theoretical results. In addition, we discuss the effect of HTLV-I infection on the HIV dynamics and vice versa.     

Global dynamics of delay‐distributed HIV infection models with differential drug efficacy in cocirculating target cells

In this paper, we investigate the dynamical behaviors of three human immunodeficiency virus infection models with two types of cocirculating target cells and distributed intracellular delay. The models take into account both short‐lived infected cells and long‐lived chronically infected cells. In the two types of target cells, the drug efficacy is assumed to be different. The incidence rate of infection is given by bilinear and saturation functional responses in the first and second models, respectively, while it is given by a general function in the third model. Lyapunov functionals are constructed and LaSalle invariance principle is applied to prove the global asymptotic stability of all equilibria of the models. We have derived the basic reproduction number R₀ for the three models. For the first two models, we have proven that the disease‐free equilibrium is globally asymptotically stable (GAS) when R₀≤1, and the endemic equilibrium is GAS when R₀>1. For the third model, we have established a set of sufficient conditions for global stability of both equilibria of the model. We have checked our theoretical results with numerical simulations. Copyright © 2015 John Wiley & Sons, Ltd.     

TWO DIFFERENTIAL INFECTIVITY EPIDEMIC MODELS WITH NONLINEAR INCIDENCE RATE

This paper considers two differential infectivity（DI） epidemic models with a nonlinear incidence rate and constant or varying population size. The models exhibits two equilibria, namely., a disease-free equilibrium O and a unique endemic equilibrium. If the basic reproductive number σ is below unity,O is globally stable and the disease always dies out. If σ〉1, O is unstable and the sufficient conditions for global stability of endemic equilibrium are derived. Moreover,when σ〈 1 ,the local or global asymptotical stability of endemic equilibrium for DI model with constant population size in n-dimensional or two-dimensional space is obtained.     

Stability of Equilibria for the $\mathfrak{so}(4)$ Free Rigid Body

The stability for all generic equilibria of the Lie–Poisson dynamics of the \mathfrakso(4)\mathfrak{so}(4) rigid body dynamics is completely determined. It is shown that for the generalized rigid body certain Cartan subalgebras (called of coordinate type) of \mathfrakso(n)\mathfrak{so}(n) are equilibrium points for the rigid body dynamics. In the case of \mathfrakso(4)\mathfrak{so}(4) there are three coordinate type Cartan subalgebras whose intersection with a regular adjoint orbit gives three Weyl group orbits of equilibria. These coordinate type Cartan subalgebras are the analogues of the three axes of equilibria for the classical rigid body in \mathfrakso(3)\mathfrak{so}(3). In addition to these coordinate type Cartan equilibria there are others that come in curves.     

Global bifurcation of travelling waves in discrete nonlinear Schrödinger equations

The discrete nonlinear Schrödinger equations of n sites are studied with periodic boundary conditions. These equations have n branches of standing waves that bifurcate from zero. Travelling waves appear as a symmetry-breaking from the standing waves for different amplitudes. The bifurcation is proved using the global Rabinowitz alternative in subspaces of symmetric functions. Applications to the Schrödinger and Saturable lattices are presented.     

Stability preserving NSFD scheme for a multi‐group SVIR epidemic model

A discrete multi‐group SVIR epidemic model with general nonlinear incidence rate and vaccination is investigated by utilizing Mickens' nonstandard finite difference scheme to a corresponding continuous model. Mathematical analysis shows that the global asymptotic stability of the equilibria is fully determined by the basic reproduction number by constructing Lyapunov functions. The results imply that the discretization scheme can efficiently preserves the global asymptotic stability of the equilibria for corresponding continuous model, and numerical simulations are carried out to illustrate the theoretical results. Copyright © 2017 John Wiley & Sons, Ltd.     

Global optimization of mixed-integer bilevel programming problems

Two approaches that solve the mixed-integer nonlinear bilevel programming problem to global optimality are introduced. The first addresses problems mixed-integer nonlinear in outer variables and C²-nonlinear in inner variables. The second adresses problems with general mixed-integer nonlinear functions in outer level. Inner level functions may be mixed-integer nonlinear in outer variables, linear, polynomial, or multilinear in inner integer variables, and linear in inner continuous variables. This second approach is based on reformulating the mixed-integer inner problem as continuous via its vertex polyheral convex hull representation and solving the resulting nonlinear bilevel optimization problem by a novel deterministic global optimization framework. Computational studies illustrate proposed approaches.     

Modeling diseases with latency and nonlinear incidence rates: global dynamics of a multi‐group model

In this paper, we perform global stability analysis of a multi‐group SEIR epidemic model in which we can consider the heterogeneity of host population and the effects of latency and nonlinear incidence rates. For a simpler version that assumes an identical natural death rate for all groups, and with a gamma distribution for the latency, the basic reproduction number is defined by the theory of the next generation operator and proved to be a sharp threshold determining whether or not disease spread. Under certain assumptions, the disease‐free equilibrium is globally asymptotically stable if R₀≤1 and there exists a unique endemic equilibrium which is globally asymptotically stable if R₀>1. The proofs of global stability of equilibria exploit a matrix‐theoretic method using Perron eigenvetor, a graph‐theoretic method based on Kirchhoff's matrix tree theorem and Lyapunov functionals. Copyright © 2015 John Wiley & Sons, Ltd.     

The existence of asymptotically stable solutions for a nonlinear functional integral equation with values in a general Banach space

Motivated by the recent known results about the solvability and existence of asymptotically stable solutions for nonlinear functional integral equations in spaces of functions defined on unbounded intervals with values in the n-dimensional real space, we establish asymptotically stable solutions for a nonlinear functional integral equation in the space of all continuous functions on R₊ with values in a general Banach space, via a fixed point theorem of Krasnosel’skii type. In order to illustrate the result obtained here, an example is given.     

Piecewise-Linear Approximations of Multidimensional Functions

We develop explicit, piecewise-linear formulations of functions f(x):ℝ ⁿ ↦ℝ, n≤3, that are defined on an orthogonal grid of vertex points. If mixed-integer linear optimization problems (MILPs) involving multidimensional piecewise-linear functions can be easily and efficiently solved to global optimality, then non-analytic functions can be used as an objective or constraint function for large optimization problems. Linear interpolation between fixed gridpoints can also be used to approximate generic, nonlinear functions, allowing us to approximately solve problems using mixed-integer linear optimization methods. Toward this end, we develop two different explicit formulations of piecewise-linear functions and discuss the consequences of integrating the formulations into an optimization problem.     

Canonical dual least square method for solving general nonlinear systems of quadratic equations

This paper presents a canonical dual approach for solving general nonlinear algebraic systems. By using least square method, the nonlinear system of m-quadratic equations in n-dimensional space is first formulated as a nonconvex optimization problem. We then proved that, by the canonical duality theory developed by the second author, this nonconvex problem is equivalent to a concave maximization problem in ℝ ^m , which can be solved easily by well-developed convex optimization techniques. Both existence and uniqueness of global optimal solutions are discussed, and several illustrative examples are presented.     

Global dynamics of a multistage SIR model with distributed delays and nonlinear incidence rate

In this paper, a multistage susceptible‐infectious‐recovered model with distributed delays and nonlinear incidence rate is investigated, which extends the model considered by Guo et al. [H. Guo, M. Y. Li and Z. Shuai, Global dynamics of a general class of multistage models for infectious diseases, SIAM J. Appl. Math., 72 (2012), 261–279]. Under some appropriate and realistic conditions, the global dynamics is completely determined by the basic reproduction number R₀. If R₀≤1, then the infection‐free equilibrium is globally asymptotically stable and the disease dies out in all stages. If R₀>1, then a unique endemic equilibrium exists, and it is globally asymptotically stable, and hence the disease persists in all stages. The results are proved by utilizing the theory of non‐negative matrices, Lyapunov functionals, and the graph‐theoretical approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Global stability of an SIS epidemic model with delays

In this paper, we consider the global dynamics of the S(E)IS model with delays denoting an incubation time. By constructing a Lyapunov functional, we prove stability of a disease‐free equilibrium E₀ under a condition different from that in the recent paper. Then we claim that R₀≤1 is a necessary and sufficient condition under which E₀ is globally asymptotically stable. We also propose a discrete model preserving positivity and global stability of the same equilibria as the continuous model with distributed delays, by means of discrete analogs of the Lyapunov functional.