Hopf bifurcations in multiple-parameter space of the Hodgkin-Huxley equations II. Singularity theoretic approach and highly degenerate bifurcations

In the Hodgkin-Huxley equations (HH), we have identified the parameter regions in which either two stable periodic solutions with different amplitudes and periods and an equilibrium point or two stable periodic solutions coexist. The global structure of bifurcations in the multiple-parameter space in the HH suggested that the bistabilities of the periodic solutions are associated with the degenerate Hopf bifurcation points by which several qualitatively different behaviors are organized. In this paper, we clarify this by analyzing the details of the degenerate Hopf bifurcations using the singularity theory approach which deals with local bifurcations near a highly degenerate fixed point. Received: 23 April 1999 / Accepted in revised form: 24 September 1999     

Stability of periodic solutions for a model of genetic repression with delays

A technique is discussed for locating the Hopf bifurcation of an n-dimensional system of delay differential equations which arises from a model for control of protein biosynthesis. Certain parameter values are shown to allow a Hopf bifurcation to periodic orbits. At the Hopf bifurcation the periodic orbits are shown to be stable either analytically or numerically depending on the parameter values.On leave from North Carolina State University.Supported in part by N.S.F. Grant # MCS 81-02828     

Mathematical modelling of dynamics and control in metabolic networks: VI. Dynamic bifurcations in single biochemical control loops

A dynamic stability analysis of an extended form of the Goodwin equations is presented. The Goodwin equations are extended to include Michaelis-Menten kinetics for the removal of the end-product. Inclusion of saturation kinetic behavior substantially increases the likelihood of dynamic instability in this model control loop. Oscillations are found for reaction chains of low order, as low as second order, and low degrees of co-operativity, as low as v = 2, simultaneously, thus indicating that dynamic instability in this system exists for physiologically realistic parameter values. The branches of bifurcated solutions are computed numerically and unstable Hopf bifurcations are found. Further, solution branches from stable Hopf bifurcation points are found to "fold back", i.e. have periodic limit points, producing situations where multiple stable limit cycles exist.     

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 analysis of a model for hormonal regulation of the menstrual cycle

A model for hormonal control of the menstrual cycle with 13 ordinary differential equations and 41 parameters is presented. Important changes in model behavior result from variations in two of the most sensitive parameters. One parameter represents the level of estradiol sufficient for significant synthesis of luteinizing hormone, which causes ovulation. By studying bifurcation diagrams in this parameter, an interval of parameter values is observed for which a unique stable periodic solution exists and it represent an ovulatory cycle. The other parameter measures mass transfer between the first two stages of ovarian development and is indicative of healthy follicular growth. Changes in this parameter affect the uniqueness interval defined with respect to the first parameter. Hopf, saddle-node and transcritical bifurcations are examined. To attain a normal ovulatory menstrual cycle in this model, a balance must be maintained between healthy development of the follicles and flexibility in estradiol levels needed to produce the surge in luteinizing hormone.     

Modelling the spatio-temporal dynamics of multi-species host-parasitoid interactions: heterogeneous patterns and ecological implications

A mathematical model of the spatio-temporal dynamics of a two host, two parasitoid system is presented. There is a coupling of the four species through parasitism of both hosts by one of the parasitoids. The model comprises a system of four reaction-diffusion equations. The underlying system of ordinary differential equations, modelling the host-parasitoid population dynamics, has a unique positive steady state and is shown to be capable of undergoing Hopf bifurcations, leading to limit cycle kinetics which give rise to oscillatory temporal dynamics. The stability of the positive steady state has a fundamental impact on the spatio-temporal dynamics: stable travelling waves of parasitoid invasion exhibit increasingly irregular periodic travelling wave behaviour when key parameter values are increased beyond their Hopf bifurcation point. These irregular periodic travelling waves give rise to heterogeneous spatio-temporal patterns of host and parasitoid abundance. The generation of heterogeneous patterns has ecological implications and the concepts of temporary host refuge and niche formation are considered.     

Ionic channel density of excitable membranes can act as a bifurcation parameter

As the maximal K⁺-conductance (or K⁺-channel density) of the Hodgkin-Huxley equations is reduced, the stable resting membrane potential bifurcates at a subcritical Hopf bifurcation into small amplitude unstable oscillations. These small amplitude solutions jump to large amplitude periodic solutions that correspond to a repetitive discharge of action potentials. Thus the specific channel density can act as a bifurcation parameter, and can control the excitability and autorhythmicity of excitable membranes.     

Phase resetting and bifurcation in the ventricular myocardium.

With the dynamic differential equations of Beeler, G. W., and H. Reuter (1977, J. Physiol. [Lond.]. 268:177-210), we have studied the oscillatory behavior of the ventricular muscle fiber stimulated by a depolarizing applied current I app. The dynamic solutions of BR equations revealed that as I app increases, a periodic repetitive spiking mode appears above the subthreshold I app, which transforms to a periodic spiking-bursting mode of oscillations, and finally to chaos near the suprathreshold I app (i.e., near the termination of the periodic state). Phase resetting and annihilation of repetitive firing in the ventricular myocardium were demonstrated by a brief current pulse of the proper magnitude applied at the proper phase. These phenomena were further examined by a bifurcation analysis. A bifurcation diagram constructed as a function of I app revealed the existence of a stable periodic solution for a certain range of current values. Two Hopf bifurcation points exist in the solution, one just above the lower periodic limit point and the other substantially below the upper periodic limit point. Between each periodic limit point and the Hopf bifurcation, the cell exhibited the coexistence of two different stable modes of operation; the oscillatory repetitive firing state and the time-independent steady state. As in the Hodgkin-Huxley case, there was a low amplitude unstable periodic state, which separates the domain of the stable periodic state from the stable steady state. Thus, in support of the dynamic perturbation methods, the bifurcation diagram of the BR equation predicts the region where instantaneous perturbations, such as brief current pulses, can send the stable repetitive rhythmic state into the stable steady state.     

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.     

Isolated periodic solutions of the Hodgkin-Huxley equations

Periodic solutions of the current clamped Hodgkin-Huxley equations (Hodgkin & Huxley, 1952 J. Physiol. 117, 500) that arise by degenerate Hopf bifurcation were studied recently by Labouriau (1985 SIAM J. Math. Anal. 16, 1121, 1987 Degenerate Hopf Bifurcation and Nerve Impulse (Part II), in press). Two parameters, temperature T and sodium conductance gNa were varied from the original values obtained by Hodgkin & Huxley. Labouriau's work proved the existence of small amplitude periodic solution branches that do not connect locally to the stationary solution branch, and had not been previously computed. In this paper we compute these solution branches globally. We find families of isolas of periodic solutions (i.e. branches not connected to the stationary branch). For values of gNa in the range measured by Hodgkin & Huxley, and for physically reasonable temperatures, there are isolas containing orbitally asymptotically stable solutions. The presence of isolas of periodic solutions suggests that in certain current space clamped membrane experiments, action potentials could be observed even though the stationary state is stable for all current stimuli. Once produced, such action potentials will disappear suddenly if the current stimulus is either increased or decreased past certain values. Under some conditions, "jumping" between action potentials of different amplitudes might be observed.     

Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models

When the traditional assumption that the incidence rate is proportional to the product of the numbers of infectives and susceptibles is dropped, the SIRS model can exhibit qualitatively different dynamical behaviors, including Hopf bifurcations, saddle-node bifurcations, and homoclinic loop bifurcations. These may be important epidemiologically in that they demonstrate the possibility of infection outbreak and collapse, or autonomous periodic coexistence of disease and host. The possible mechanisms leading to nonlinear incidence rates are discussed. Finally, a modified general criterion for supercritical or subcritical Hopf bifurcation of 2-dimensional systems is presented.     

The pyrite iron cycle catalyzed by Acidithiobacillus ferrooxidans

In this paper, we study a model of the biotic pyrite iron cycle catalyzed by bacteria Acidithiobacillus ferrooxidans, in mining activity sites waste dumps. Chemical reactions, reaction rates and the population growth model are mostly taken from the existing literature. Analysis of the corresponding dynamical system shows the existence of up to four non-trivial steady state solutions. The stability of the equilibria solutions is determined, finding up to two coexisting stable solutions. Two Hopf bifurcations and a region of parameter space in which there are stable periodic solutions are found. In addition, we find a homoclinic bifurcation which gives rise to a stable periodic orbit of large period. The existence of these stable oscillatory solutions gives a possible explanation for the oscillations of bacteria concentration and pH for the iron cycle, described in Jaynes et al. (Water Resour Res 20:233–242, 1984).     

Degenerate Hopf bifurcation and isolated periodic solutions of the Hodgkin-Huxley model with varying sodium ion concentration

Points of degenerate Hopf bifurcation in the Hodgkin-Huxley model are found as parameters temperature T and voltage level of sodium VNa are varied. Local techniques of degenerate Hopf bifurcation analysis are used to show the existence of families of periodic solutions of the model: isolated branches of periodic solutions (i.e. branches not connected to the stationary branch) are found in addition to Hopf branches. Purely numerical techniques are used to show that the isolas persist for VNa up to a value slightly greater than 114 mV. Under some conditions there are multiple stable periodic solutions, so "jumping" between action potentials of different amplitudes might be observed.     

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.     

Periodic metabolic systems: Oscillations in multiple-loop negative feedback biochemical control networks

Summary For a general multiple loop feedback inhibition system in which the end product can inhibit any or all of the intermediate reactions it is shown that biologically significant behaviour is always confined to a bounded region of reaction space containing a unique equilibrium. By explicit construction of a Liapunov function for the general n dimensional differential equation it is shown that some values of reaction parameters cause the concentration vector to approach the equilibrium asymptotically for all physically realizable initial conditions. As the parameter values change, periodic solutions can appear within the bounded region. Some information about these periodic solutions can be obtained from the Hopf bifurcation theorem. Alternatively, if specific parameter values are known a numerical method can be used to find periodic solutions and determine their stability by locating a zero of the displacement map. The single loop Goodwin oscillator is analysed in detail. The methods are then used to treat an oscillator with two feedback loops and it is found that oscillations are possible even if both Hill coefficients are equal to one.     

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.     

Stability switches,oscillatory multistability,and spatio-temporal patterns of nonlinear oscillations in recurrently delay coupled neural networks

A model of time-delay recurrently coupled spatially segregated neural assemblies is here proposed. We show that it operates like some of the hierarchical architectures of the brain. Each assembly is a neural network with no delay in the local couplings between the units. The delay appears in the long range feedforward and feedback inter-assemblies communications. Bifurcation analysis of a simple four-units system in the autonomous case shows the richness of the dynamical behaviors in a biophysically plausible parameter region. We find oscillatory multistability, hysteresis, and stability switches of the rest state provoked by the time delay. Then we investigate the spatio-temporal patterns of bifurcating periodic solutions by using the symmetric local Hopf bifurcation theory of delay differential equations and derive the equation describing the flow on the center manifold that enables us determining the direction of Hopf bifurcations and stability of the bifurcating periodic orbits. We also discuss computational properties of the system due to the delay when an external drive of the network mimicks external sensory input.     

Periodic orbits arising from Hopf bifurcations in a Volterra prey-predator model

 In this paper we derive a formula which enables the stability of periodic solutions to a Volterra integro-differential system to be determined. This system which has been studied by Cushing [1], models a predator-prey interaction with distributed delays. The results are obtained by using the algorithm developed by Kazarinoff, Wan and van den Driessche [2] based on the centre manifold formulas of Hassard and Wan [3]. We discuss an example of the formula for the case of weak kernels and show that under certain conditions stable periodic solutions arising from Hopf bifurcations at different critical values of the parameters can exist together. Received 30 December 1994; received in revised form 12 December 1995     

Transient oscillations induced by delayed growth response in the chemostat

In this paper, in order to try to account for the transient oscillations observed in chemostat experiments, we consider a model of single species growth in a chemostat that involves delayed growth response. The time delay models the lag involved in the nutrient conversion process. Both monotone response functions and nonmonotone response functions are considered. The nonmonotone response function models the inhibitory effects of growth response of certain nutrients when concentrations are too high. By applying local and global Hopf bifurcation theorems, we prove that the model has unstable periodic solutions that bifurcate from unstable nonnegative equilibria as the parameter measuring the delay passes through certain critical values and that these local periodic solutions can persist, even if the delay parameter moves far from the critical (local) bifurcation values.When there are two positive equilibria, then positive periodic solutions can exist. When there is a unique positive equilibrium, the model does not have positive periodic oscillations and the unique positive equilibrium is globally asymptotically stable. However, the model can have periodic solutions that change sign. Although these solutions are not biologically meaningful, provided the initial data starts close enough to the unstable manifold of one of these periodic solutions they may still help to account for the transient oscillations that have been frequently observed in chemostat experiments. Numerical simulations are provided to illustrate that the model has varying degrees of transient oscillatory behaviour that can be controlled by the choice of the initial data.Mathematics Subject Classification: 34D20, 34K20, 92D25Research was partially supported by NSERC of Canada.This work was partly done while this author was a postdoc at McMaster.     

Coexistence of Multiple Periodic and Chaotic Regimes in Biochemical Oscillations with Phase Shifts

The numerical study of a glycolytic model formed by a system of three delay differential equations reveals a multiplicity of stable coexisting states: birhythmicity, trirhythmicity, hard excitation and quasiperiodic with chaotic regimes. For different initial functions in the phase space one may observe the coexistence of two different quasiperiodic motions, the existence of a stable steady state with a stable torus, and the existence of a strange attractor with different stable regimes (chaos with torus, chaos with bursting motion, and chaos with different periodic regimes). For a single range of the control parameter values our system may exhibit different bifurcation diagrams: in one case a Feigenbaum route to chaos coexists with a finite number of successive periodic bifurcations, in other conditions it is possible to observe the coexistence of two quasiperiodicity routes to chaos. These studies were obtained both at constant input flux and under forcing conditions.