Hopf bifurcations in multiple-parameter space of the Hodgkin-Huxley equations I. Global organization of bistable periodic solutions

 The Hodgkin-Huxley equations (HH) are parameterized by a number of parameters and shows a variety of qualitatively different behaviors depending on the parameter values. We explored the dynamics of the HH for a wide range of parameter values in the multiple-parameter space, that is, we examined the global structure of bifurcations of the HH. Results are summarized in various two-parameter bifurcation diagrams with I _ext (externally applied DC current) as the abscissa and one of the other parameters as the ordinate. In each diagram, the parameter plane was divided into several regions according to the qualitative behavior of the equations. In particular, we focused on periodic solutions emerging via Hopf bifurcations and identified parameter regions in which either two stable periodic solutions with different amplitudes and periods and a stable equilibrium point or two stable periodic solutions coexist. Global analysis of the bifurcation structure suggested that generation of these regions is associated with degenerate Hopf bifurcations. Received: 23 April 1999 / Accepted in revised form: 24 September 1999     

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     

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.     

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.     

具有三时滞Lotka-Volterra互惠系统的Hopf分支

建立了具有三个时滞的Lotka-Volterra互惠系统;获得了正平衡点和Hopf分支存在的条件等;并对所获得的结果进行了数值模拟．     

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).     

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.     

Emergence of oscillations in a model of weakly coupled two Bonhoeffer-van der Pol equations

Bifurcations of periodic solutions in a model of weakly coupled two Bonhoeffer-van der Pol equations are studied. The model realizes a half-center model with reciprocal inhibition, a typical model used in the field of neural motor control to account for the generation of alternating rhythmic bursts observed in motoneurons and spinal neural networks. Several oscillatory solutions such as in-phase, anti-phase as well as out-of-phase solutions emerge from the model's equilibrium as one of the parameters of the model changes. Among the variety of bifurcations exhibited by the model, we analyze Hopf bifurcations, by which several periodic solutions emerge, and illustrate generation mechanisms of alternating oscillations in the model.     

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.     

The effect of predator limitation on the dynamics of simple food chains

We investigate the influence of competition between predators on the dynamics of bitrophic predator–prey systems and of tritrophic food chains. Competition between predators is implemented either as interference competition, or as a density-dependent mortality rate. With interference competition, the paradox of enrichment is reduced or completely suppressed, but otherwise, the dynamical behavior of the systems is not fundamentally different from that of the Rosenzweig–MacArthur model, which contains no predator competition and shows only continuous transitions between fixed points or periodic oscillations. In contrast, with density-dependent predator mortality, the system shows a surprisingly rich dynamical behavior. In particular, decreasing the density regulation of the predator can induce catastrophic shifts from a stable fixed point to a large oscillation where the predator chases the prey through a cycle that brings both species close to the threshold of extinction. Other catastrophic bifurcations, such as subcritical Hopf bifurcations and saddle-node bifurcations of limit cycles, do also occur. In tritrophic food chains, we find again that fixed points in the model with predator interference become unstable only through Hopf bifurcations, which can also be subcritical, in contrast to the bitrophic situation. The model with a density limitation shows again catastrophic destabilization of fixed points and various nonlocal bifurcations. In addition, chaos occurs for both models in appropriate parameter ranges.     

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.     

Local and global bifurcations at infinity in models of glycolytic oscillations

 We investigate two models of glycolytic oscillations. Each model consists of two coupled nonlinear ordinary differential equations. Both models are found to have a saddle point at infinity and to exhibit a saddle-node bifurcation at infinity, giving rise to a second saddle and a stable node at infinity. Depending on model parameters, a stable limit cycle may blow up to infinite period and amplitude and disappear in the bifurcation, and after the bifurcation, the stable node at infinity then attracts all trajectories. Alternatively, the stable node at infinity may coexist with either a stable sink (not at infinity) or a stable limit cycle. This limit cycle may then disappear in a heteroclinic bifurcation at infinity in which the unstable manifold from one saddle at infinity joins the stable manifold of the other saddle at infinity. These results explain prior reports for one of the models concerning parameter values for which the system does not admit any physical (bounded) behavior. Analytic results on the scaling of amplitude and period close to the bifurcations are obtained and confirmed by numerical computations. Finally, we consider more realistic modified models where all solutions are bounded and show that some of the features stemming from the bifurcations at infinity are still present. Received 4 September 1995; received in revised form 18 September 1996     

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.     

一类具有时滞的传染病模型的稳定性分析

研究了一类具有时滞的传染病生物模型．首先研究了该模型的线性稳定性,并给出了一列Hopf分支值,然后利用中心流形定理和正规型方法,给出了确定分支周期解的分支方向与稳定性的计算公式．     

Bifurcation analysis of periodic SEIR and SIR epidemic models

The bifurcations of the periodic solutions of SEIR and SIR epidemic models with sinusoidally varying contact rate are investigated. The analysis is carried out with respect to two parameters: the mean value and the degree of seasonality of the contact rate. The corresponding portraits in the two-parameter space are obtained by means of a numerical continuation method. Codimension two bifurcations (degenerate flips and cusps) are detected, and multiple stable modes of behavior are identified in various regions of the parameter space. Finally, it is shown how the parametric portrait of the SEIR model tends to that of the SIR model when the latent period tends to zero.     

Dynamics and bifurcation of a model for hormonal control of the menstrual cycle with inhibin delay

A system of 13 ordinary differential equations with 42 parameters is presented to model hormonal regulation of the menstrual cycle. For an excellent fit to clinical data, the model requires a 36 h time delay for the effect of inhibin on the synthesis of follicle stimulating hormone. Biological and mathematical reasons for this delay are discussed. Bifurcations with respect to changes in three important parameters are examined. One parameter represents the level of estradiol adequate for significant synthesis of luteinizing hormone. Bifurcation diagrams with respect to this parameter reveal an interval of parameter values for which a unique stable periodic solution exists and this solution represents a menstrual cycle during which ovulation occurs. The second parameter measures mass transfer between the first two stages of ovarian development and is indicative of healthy follicular growth. The third parameter is the time delay. Changes in the second parameter and the time delay affect the size of the uniqueness interval defined with respect to the first parameter. Saddle-node, transcritical and degenerate Hopf bifurcations are studied.     

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.     

Bifurcation analysis of the Nowak-Bangham model in CTL dynamics

This paper investigates the local bifurcations of a CTL response model published by Nowak and Bangham [M.A. Nowak, C.R.M. Bangham, Population dynamics of immune responses to persistent viruses, Science 272 (1996) 74]. The Nowak-Bangham model can have three equilibria depending on the basic reproduction number, and generates a Hopf bifurcation through two bifurcations of equilibria. The main result shows a sufficient condition for the interior equilibrium to have a unique bifurcation point at which a simple Hopf bifurcation occurs. For this proof, some new techniques are developed in order to apply the method established by Liu [W.M. Liu, Criterion of Hopf bifurcations without using eigenvalues, J. Math. Anal. Appl. 182 (1) (1994) 250]. In addition, to demonstrate the result obtained theoretically, some bifurcation diagrams are presented with numerical examples.     

Instabilities in spatially extended predator-prey systems: spatio-temporal patterns in the neighborhood of Turing-Hopf bifurcations

We investigate the emergence of spatio-temporal patterns in ecological systems. In particular, we study a generalized predator-prey system on a spatial domain. On this domain diffusion is considered as the principal process of motion. We derive the conditions for Hopf and Turing instabilities without specifying the predator-prey functional responses and discuss their biological implications. Furthermore, we identify the codimension-2 Turing-Hopf bifurcation and the codimension-3 Turing-Takens-Bogdanov bifurcation. These bifurcations give rise to complex pattern formation processes in their neighborhood. Our theoretical findings are illustrated with a specific model. In simulations a large variety of different types of long-term behavior, including homogenous distributions, stationary spatial patterns and complex spatio-temporal patterns, are observed.     

Hopf and torus bifurcations,torus destruction and chaos in population biology

One of the simplest population biological models displaying a Hopf bifurcation is the Rosenzweig–MacArthur model with Holling type II response function as essential ingredient. In seasonally forced versions the fixed point on one side of the Hopf bifurcation becomes a limit cycle and the Hopf limit cycle on the other hand becomes a torus, hence the Hopf bifurcation becomes a torus bifurcation, and via torus destruction by further increasing relevant parameters can follow deterministic chaos. We investigate this route to chaos also in view of stochastic versions, since in real world systems only such stochastic processes would be observed.However, the Holling type II response function is not directly related to a transition from one to another population class which would allow a stochastic version straight away. Instead, a time scale separation argument leads from a more complex model to the simple 2 dimensional Rosenzweig–MacArthur model, via additional classes of food handling and predators searching for prey. This extended model allows a stochastic generalization with the stochastic version of a Hopf bifurcation, and ultimately also with additional seasonality allowing a torus bifurcation under stochasticity.Our study shows that the torus destruction into chaos with positive Lyapunov exponents can occur in parameter regions where also the time scale separation and hence stochastic versions of the model are possible. The chaotic motion is observed inside Arnol’d tongues of rational ratio of the forcing frequency and the eigenfrequency of the unforced Hopf limit cycle.Such torus bifurcations and torus destruction into chaos are also observed in other population biological systems, and were for example found in extended multi-strain epidemiological models on dengue fever. To understand such dynamical scenarios better also under noise the present low dimensional system can serve as a good study case.