Nonlinear Oscillations in Models of Immune Responses to Persistent Viruses

This paper compares several population models for HIV dynamics. The effect of perfect or nonperfect drugs proposed by Ho, Perelson, and others are considered in the more complicated models of Nowak and Bangham for immune responses to persistent viruses. It is pointed out that two models proposed by Nowak and Bangham have different dynamical behaviors. One of them can have sustained oscillations via the Hopf bifurcation and the other does not. To overcome the difficulty in symbolic computations for complicated systems, we use parameter transforms and introduce a semisymbolic method. This method is very efficient in finding Hopf bifurcation parameters and can be used for other models.     

Multiparametric bifurcations of an epidemiological model with strong Allee effect

In this paper we completely study bifurcations of an epidemic model with five parameters introduced by Hilker et al. (Am Nat 173:72–88, 2009), which describes the joint interplay of a strong Allee effect and infectious diseases in a single population. Existence of multiple positive equilibria and all kinds of bifurcation are examined as well as related dynamical behavior. It is shown that the model undergoes a series of bifurcations such as saddle-node bifurcation, pitchfork bifurcation, Bogdanov–Takens bifurcation, degenerate Hopf bifurcation of codimension two and degenerate elliptic type Bogdanov–Takens bifurcation of codimension three. Respective bifurcation surfaces in five-dimensional parameter spaces and related dynamical behavior are obtained. These theoretical conclusions confirm their numerical simulations and conjectures by Hilker et al., and reveal some new bifurcation phenomena which are not observed in Hilker et al. (Am Nat 173:72–88, 2009). The rich and complicated dynamics exhibit that the model is very sensitive to parameter perturbations, which has important implications for disease control of endangered species.     

Bifurcations in a white-blood-cell production model

We study the dynamics of a model of white-blood-cell (WBC) production. The model consists of two compartmental differential equations with two discrete delays. We show that from normal to pathological parameter values, the system undergoes supercritical Hopf bifurcations and saddle-node bifurcations of limit cycles. We characterize the steady states of the system and perform a bifurcation analysis. Our results indicate that an increase in apoptosis rate of either hematopoietic stem cells or WBC precursors induces a Hopf bifurcation and an oscillatory regime takes place. These oscillations are seen in some hematological diseases.     

The role of transmissible diseases in the Holling-Tanner predator-prey model

Transmissible diseases are known to induce remarkable major behavioral changes in predator-prey systems. However, little attention has been paid to model such situations. The latter would allow to predict useful applications in both dynamics and control. Here the Holling-Tanner model is revisited to account for the influence of a transmissible disease, under the assumption that it spreads among the prey species only. We have found the equilibria and analyzed the behavior of the system around each one of them. A threshold result determining when the disease dies out has been identified. We also investigated the parametric space under which the system enters into Hopf and transcritical bifurcations, around the disease free equilibrium. The system is shown to experience neither saddle-node nor pitch-fork bifurcation. Global stability results are obtained by constructing suitable Lyapunov functions.     

Oscillatory viral dynamics in a delayed HIV pathogenesis model

We consider an HIV pathogenesis model incorporating antiretroviral therapy and HIV replication time. We investigate the existence and stability of equilibria, as well as Hopf bifurcations to sustained oscillations when drug efficacy is less than 100%. We derive sufficient conditions for the global asymptotic stability of the uninfected steady state. We show that time delay has no effect on the local asymptotic stability of the uninfected steady state, but can destabilize the infected steady state, leading to a Hopf bifurcation to periodic solutions in the realistic parameter ranges.     

The Dynamics of Growing Islets and Transmission of Schistosomiasis Japonica in the Yangtze River

We formulate and analyze a system of ordinary differential equations for the transmission of schistosomiasis japonica on the islets in the Yangtze River, China. The impact of growing islets on the spread of schistosomiasis is investigated by the bifurcation analysis. Using the projection technique developed by Hassard, Kazarinoff and Wan, the normal form of the cusp bifurcation of codimension 2 is derived to overcome the technical difficulties in studying the existence, stability, and bifurcation of the multiple endemic equilibria in high-dimensional phase space. We show that the model can also undergo transcritical bifurcations, saddle-node bifurcations, a pitchfork bifurcation, and Hopf bifurcations. The bifurcation diagrams and epidemiological interpretations are given. We conclude that when the islet reaches a critical size, the transmission cycle of the schistosomiasis japonica between wild rats Rattus norvegicus and snails Oncomelania hupensis could be established, which serves as a possible source of schistosomiasis transmission along the Yangtze River.     

Bistability and oscillations in chemical reaction networks

Bifurcation theory is one of the most widely used approaches for analysis of dynamical behaviour of chemical and biochemical reaction networks. Some of the interesting qualitative behaviour that are analyzed are oscillations and bistability (a situation where a system has at least two coexisting stable equilibria). Both phenomena have been identified as central features of many biological and biochemical systems. This paper, using the theory of stoichiometric network analysis (SNA) and notions from algebraic geometry, presents sufficient conditions for a reaction network to display bifurcations associated with these phenomena. The advantage of these conditions is that they impose fewer algebraic conditions on model parameters than conditions associated with standard bifurcation theorems. To derive the new conditions, a coordinate transformation will be made that will guarantee the existence of branches of positive equilibria in the system. This is particularly useful in mathematical biology, where only positive variable values are considered to be meaningful. The first part of the paper will be an extended introduction to SNA and algebraic geometry-related methods which are used in the coordinate transformation and set up of the theorems. In the second part of the paper we will focus on the derivation of bifurcation conditions using SNA and algebraic geometry. Conditions will be derived for three bifurcations: the saddle-node bifurcation, a simple branching point, both linked to bistability, and a simple Hopf bifurcation. The latter is linked to oscillatory behaviour. The conditions derived are sufficient and they extend earlier results from stoichiometric network analysis as can be found in (Aguda and Clarke in J Chem Phys 87:3461–3470, 1987; Clarke and Jiang in J Chem Phys 99:4464–4476, 1993; Gatermann et al. in J Symb Comput 40:1361–1382, 2005). In these papers some necessary conditions for two of these bifurcations were given. A set of examples will illustrate that algebraic conditions arising from given sufficient bifurcation conditions are not more difficult to interpret nor harder to calculate than those arising from necessary bifurcation conditions. Hence an increasing amount of information is gained at no extra computational cost. The theory can also be used in a second step for a systematic bifurcation analysis of larger reaction networks. We have added a dedication of the paper to K. Gatermann.     

Hopf bifurcation and bursting synchronization in an excitable systems with chemical delayed coupling

In this paper we consider the Hopf bifurcation and synchronization in the two coupled Hindmarsh–Rose excitable systems with chemical coupling and time-delay. We surveyed the conditions for Hopf bifurcations by means of dynamical bifurcation analysis and numerical simulation. The results show that the coupled excitable systems with no delay have supercritical Hopf bifurcation, while the delayed system undergoes Hopf bifurcations at critical time delays when coupling strength lies in a particular region. We also investigated the effect of the delay on the transition of bursting synchronization in the coupled system. The results are helpful for us to better understand the dynamical properties of excitable systems and the biological mechanism of information encoding and cognitive activity.     

The subcritical collapse of predator populations in discrete-time predator-prey models.

Many discrete-time predator-prey models possess three equilibria, corresponding to (1) extinction of both species, (2) extinction of the predator and survival of the prey at its carrying capacity, or (3) coexistence of both species. For a variety of such models, the equilibrium corresponding to coexistence may lose stability via a Hopf bifurcation, in which case trajectories approach an invariant circle. Alternatively, the equilibrium may undergo a subcritical flip bifurcation with a concomitant crash in the predator's population. We review a technique for distinguishing between subcritical and supercritical flip bifurcations and provide examples of predator-prey systems with a subcritical flip bifurcation.     

The discrete dynamics of symmetric competition in the plane

We consider the generalized Lotka-Volterra two-species system xn + 1 = xn exp(r1(1 - xn) - s1yn) yn + 1 = yn exp(r2(1 - yn) - s2xn) originally proposed by R. M. May as a model for competitive interaction. In the symmetric case that r1 = r2 and s1 = s2, a region of ultimate confinement is found and the dynamics therein are described in some detail. The bifurcations of periodic points of low period are studied, and a cascade of period-doubling bifurcations is indicated. Within the confinement region, a parameter region is determined for the stable Hopf bifurcation of a pair of symmetrically placed period-two points, which imposes a second component of oscillation near the stable cycles. It is suggested that the symmetric competitive model contains much of the dynamical complexity to be expected in any discrete two-dimensional competitive model.     

Effects of quarantine in six endemic models for infectious diseases

Thresholds, equilibria, and their stability are found for SIQS and SIQR epidemiology models with three forms of the incidence. For most of these models, the endemic equilibrium is asymptotically stable, but for the SIQR model with the quarantine-adjusted incidence, the endemic equilibrium is an unstable spiral for some parameter values and periodic solutions arise by Hopf bifurcation. The Hopf bifurcation surface and stable periodic solutions are found numerically.     

Bifurcation control of the Morris–Lecar neuron model via a dynamic state-feedback control

In this paper, we address the control problem of bifurcations in the Morris–Lecar (ML) neuron model. With the use of a dynamic state-feedback control, two Hopf bifurcation points in the ML neuron model with Type II excitability can be relocated to new desired locations simultaneously. Also, with the proposed control law, the neuronal excitability characteristics can be transformed from Type I excitability to Type II excitability by changing the type of bifurcation, in which the neuron goes from quiescence to periodic spiking from a saddle node on an invariant circle bifurcation to a Hopf bifurcation. Simulation results are provided.     

Equilibrium analysis and phase synchronization of two coupled HR neurons with gap junction

The properties of equilibria and phase synchronization involving burst synchronization and spike synchronization of two electrically coupled HR neurons are studied in this paper. The findings reveal that in the non-delayed system the existence of equilibria can be turned into intersection of two odd functions, and two types of equilibria with symmetry and non-symmetry can be found. With the stability and bifurcation analysis, the bifurcations of equilibria are investigated. For the delayed system, the equilibria remain unchanged. However, the Hopf bifurcation point is drastically affected by time delay. For the phase synchronization, we focus on the synchronization transition from burst synchronization to spike synchronization in the non-delayed system and the effect of coupling strength and time delay on spike synchronization in delayed system. In addition, corresponding firing rhythms and spike synchronized regions are obtained in the two parameters plane. The results allow us to better understand the properties of equilibria, multi-time-scale properties of synchronization and temporal encoding scheme in neuronal systems.     

Hopf bifurcation in three-species food chain models with group defense.

Three-species food-chain models, in which the prey population exhibits group defense, are considered. Using the carrying capacity of the environment as the bifurcation parameter, it is shown that the model without delay undergoes a sequence of Hopf bifurcations. In the model with delay it is shown that using a delay as a bifurcation parameter, a Hopf bifurcation can also occur in this case. These occurrences may be interpreted as showing that a region of local stability (survival) may exist even though the positive steady states are unstable. A computer code BIFDD is used to determine the stability of the bifurcation solutions of a delay model.     

Equilibrium properties of a multi-locus, haploid-selection, symmetric-viability model

Under haploid selection, a multi-locus, diallelic, two-niche Levene (1953) model is studied. Viability coefficients with symmetrically opposing directional selection in each niche are assumed, and with a further simplification that the most and least favored haplotype in each niche shares no alleles in common, and that the selection coefficients monotonically increase or decrease with the number of alleles shared. This model always admits a fully polymorphic symmetric equilibrium, which may or may not be stable.We show that a stable symmetric equilibrium can become unstable via either a supercritical or subcritical pitchfork bifurcation. In the supercritical bifurcation, the symmetric equilibrium bifurcates to a pair of stable fully polymorphic asymmetric equilibria; in the subcritical bifurcation, the symmetric equilibrium bifurcates to a pair of unstable fully polymorphic asymmetric equilibria, which then connect to either another pair of stable fully polymorphic asymmetric equilibria through saddle-node bifurcations, or to a pair of monomorphic equilibria through transcritical bifurcations. As many as three fully polymorphic stable equilibria can coexist, and jump bifurcations can occur between these equilibria when model parameters are varied.In our Levene model, increasing recombination can act to either increase or decrease the genetic diversity of a population. By generating more hybrid offspring from the mating of purebreds, recombination can act to increase genetic diversity provided the symmetric equilibrium remains stable. But by destabilizing the symmetric equilibrium, recombination can ultimately act to decrease genetic diversity.     

Periodicity in an epidemic model with a generalized non-linear incidence

We develop and analyze a simple SIV epidemic model including susceptible, infected and perfectly vaccinated classes, with a generalized non-linear incidence rate subject only to a few general conditions. These conditions are satisfied by many models appearing in the literature. The detailed dynamics analysis of the model, using the Poincaré index theory, shows that non-linearity of the incidence rate leads to vital dynamics, such as bistability and periodicity, without seasonal forcing or being cyclic. Furthermore, it is shown that the basic reproductive number is independent of the functional form of the non-linear incidence rate. Under certain, well-defined conditions, the model undergoes a Hopf bifurcation. Using the normal form of the model, the first Lyapunov coefficient is computed to determine the various types of Hopf bifurcation the model undergoes. These general results are applied to two examples: unbounded and saturated contact rates; in both cases, forward or backward Hopf bifurcations occur for two distinct values of the contact parameter. It is also shown that the model may undergo a subcritical Hopf bifurcation leading to the appearance of two concentric limit cycles. The results are illustrated by numerical simulations with realistic model parameters estimated for some infectious diseases of childhood.     

The computer-aided bifurcation analysis of predator-prey models

Mathematical analysis of dynamical systems can often benefit from accompanying numerical computations. This is particularly true if one has software (e.g. AUTO [6, 7]) capable of providing an automatic bifurcation analysis of such systems. Computer programs of this type now exist. We describe the application of such software to a predator-prey model. Phenomena that arise in this analysis include stationary bifurcations, limit points, Hopf bifurcations and secondary periodic bifurcations. A two-parameter numerical analysis leads quite naturally to the detection of higher order singularities.Supported in part by NSERC Canada (#4274) and FCAC Québec (#EQ1438)     

一个离散单种群动力系统的超临界Flip分支和Hoph分支(英文)

文中考虑一个非线性的具有年龄阶段结构的单种群离散模型,致力于揭示该系统的动力学行为,说明了系统在分支阈值附近会出现超临界Flip分支和Hoph分支.与一些文献不同的是文中用数学工具给出证明过程而不是用数值模拟的结果说明.     

Biomass and biodiversity in species-rich tritrophic communities

We consider a tritrophic system with one basal and one top species and a large number of primary consumers, and derive upper and lower bounds for the total biomass of the middle trophic level. These estimates do not depend on dynamical regime, holding for fixed point, periodic, or chaotic dynamics. We have two kinds of estimates, depending on whether the predator abundance is zero. All these results are uniform in a self-limitation parameter, which regulates prey diversity in the system. For strong self-limitation, diversity is large; for weak self-limitation, it is small. Diversity depends on the variance of species’ parameter values. The larger this variance, the lower the diversity, and vice versa. Moreover, variation in the parameters of the Holling type II functional response changes the bifurcation character, with the equilibrium state with nonzero predator abundance losing stability. If that variation is small then the bifurcation can lead to oscillations (the Hopf bifurcation). Under certain conditions, there exists a supercritical Hopf bifurcation. We then find a connection between diversity and Hopf bifurcations. We also show that the system exhibits top-down regulation and a hump-shaped diversity-productivity curve.We then extend the model by allowing species to experience self-regulation. For this extended model, explicit estimates of prey diversity are obtained. We study the dynamics of this system and find the following. First, diversity and system dynamics crucially depend on variation in species parameters. We show that under certain conditions, the system undergoes a supercritical Hopf bifurcation. We also establish a connection between diversity and Hopf bifurcations. For strong self-limitation, diversity is large and complex dynamics are absent. For weak self-limitation, diversity is small and the equilibrium with non-zero predator abundance is unstable.