Bifurcation analysis of nonlinear reaction-diffusion equations—I. Evolution equations and the steady state solutions

A model nonlinear network involving chemical reactions and diffusion is studied. The time evolution and bounds on the steady state solutions are analyzed. Spatially ordered solutions of the equations of the dissipative structure type are found by bifurcation theory. These solutions are calculated analytically and their qualitative properties are discussed.     

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.     

Exact solutions of Lotka-Volterra equations

The Lotka-Volterra nonlinear differential equations for two competing species P and Q contain six independent parameters. Their general analytic solutions, valid for arbitrary values of the parameters, are at present unknown. However, when two or more of these parameters are interrelated, it is possible to obtain the exact solutions in the P, Q phase plane, and six cases of solvability are given in this paper. The dependence of the solutions on the parameters and the initial conditions can thus be readily investigated.     

Non-existence of wave solutions for the class of reaction--diffusion equations given by the Volterra interacting-population equations with diffusion.

The system of equations is reduced to a single nonlinear parabolic equation on which a maximum principle can be used. It is then shown that the effect of uniform diffusion on the Volterra equations for any even number of interacting populations which have non-zero equilibrium values, is to damp out all spatial variations. The inclusion of population saturation terms is shown to enhance the damping process, as would be expected. The main consequence of the results is that such reaction-diffusion equations (given in section 5) cannot have physically realistic wave-like solutions, that is stable solutions, with non-negative values of the concentrations, which evolve from a time dependent solution.     

Asymptotic solutions of population dynamic equations

An approach is offered to construct the asymptomatics of the solutions on the small parameter in the close neighborhood of the equilibrium condition of the well-known Volterra-Lotka "prey-predator" system and one of its modifications which takes into account the intraspecies competition of preys and limitation of food resources of a predator. On the basis of the formulas obtained, possible dynamic modes of the size of populations of both kinds are analyzed.     

Asymptotic behaviour of solutions of the diffusive Lotka-Volterra equations

A spatially discrete version of the diffusive Lotka-Volterra equations is considered. Asymptotical spatial homogeneity of solutions of the equations with equilibrium, periodic or zero flux boundary conditions is proved without regard to crowding effects. The proof does not require the assumption of equal diffusion coefficients and the restrictions on the dimension of space and on the initial data, which are necessary in the spatially continuous model.     

Periodic solutions to nonautonomous difference equations

A technique is presented for determining when periodic solutions to nonautonomous periodic difference equations exist. Under certain constraints, stable periodic solutions can be guaranteed to exist, and this is used to compare the analogous behavior of a nonautonomous periodic hyperbolic difference equation to that of the nonautonomous periodic Pearl-Verhulst logistic differential equation.     

Travelling wave solutions of diffusive Lotka-Volterra equations

We establish the existence of travelling wave solutions for two reaction diffusion systems based on the Lotka-Volterra model for predator and prey interactions. For simplicity, we consider only 1 space dimension. The waves are of transition front type, analogous to the travelling wave solutions discussed by Fisher and Kolmogorov et al. for a scalar reaction diffusion equation. The waves discussed here are not necessarily monotone. For any speed c there is a travelling wave solution of transition front type. For one of the systems discussed here, there is a distinguished speed c^* dividing the waves into two types, waves of speed c < c^* being one type, waves of speed c ? c^* being of the other type. We present numerical evidence that for this system the wave of speed c^* is stable, and that c^* is an asymptotic speed of propagation in some sense. For the other system, waves of all speeds are in some sense stable. The proof of existence uses a shooting argument and a Lyapunov function. We also discuss some possible biological implications of the existence of these waves.     

Generalization of the Spiegler-Kedem-Katchalsky frictional model equations of the transmembrane transport for multicomponent non-electrolyte solutions.

The Spiegler-Kedem-Katchalsky frictional model equations of the transmembrane transport for systems containing n-component, non-ionic solutions is presented. The frictional interpretation of the phenomenological coefficients of membrane and the expressions connecting the practical coefficients (Lp, sigma i, omega ij) with frictional coefficients (fij) are presented.     

Properties of the solutions to certain equations arising in the central core model of the renal medulla.

We discuss the mathematical properties of the solutions to two of the differential equations that arise in the central core model of the renal medulla. Explicit solutions are obtained for a particular functional form of the flux of solute.     

Asymptotic behaviour of solutions to abstract logistic equations

We analyze the asymptotic behaviour of solutions of the abstract differential equation u'(t)=Au(t)-F(u(t))u(t)+f. Our results are applicable to models of structured population dynamics in which the state space consists of population densities with respect to the structure variables. In the equation the linear term A corresponds to internal processes independent of crowding, the nonlinear logistic term F corresponds to the influence of crowding, and the source term f corresponds to external effects. We analyze three separate cases and show that for each case the solutions stabilize in a way governed by the linear term. We illustrate the results with examples of models of structured population dynamics -- a model for the proliferation of cell lines with telomere shortening, a model of proliferating and quiescent cell populations, and a model for the growth of tumour cord cell populations.     

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     

Multiple solutions of a model describing cancerous growth

Multiple solutions of a model describing immune surveillance against cancer are studied. The model was proposed by R. Lefever and co-workers. A limit cycle solution and two kinds of wave front solutions are worked out by means of a perturbation method. The problem of nucleation concerning this model is studied by a singular-perturbation technique. In addition, the pulse wave solution is also found by numerical simulations.     

Evolution of discontinuous solutions to transport equations in stellarators

It was shown earlier that, in the range of rare collisions, transport equations for stellarators allow steady discontinuous solutions for the ambipolar electric field and for the plasma density and temperature gradients. Moreover, such solutions are non-single-valued; that is, their explicit form depends on the initial values of the ambipolar electric field. The time-independent transport equations are derived under the conventional quasineutrality condition; i.e., it is assumed that the electron and ion densities, N _e and N _i , are related by the relationship N _e = ZN _i (where Z is the ion charge number). In other words, the plasma charge density is assumed to be much less than the product e _i N _i . Under typical conditions, the corresponding inequality is satisfied by a large margin. However, if the electric field E has discontinuities, then it can be seen from the equation ▿·E = 4πρ that, at the discontinuity points, the charge density becomes infinite and the relationship N _e = ZN _i fails to hold, so it is necessary to replace it with N _e = ZN _i + ρ/e _e . In the transport equations, this latter replacement produces additional terms, proportional to the second radial derivative of the field E. With these additional terms, the steady solutions are modified substantially. First, the ambipolar field and the derivatives of the density and temperatures all become continuous functions of the coordinates, a result that seems to be quite obvious. The second, not-so-obvious result is that the steady solutions become single-valued, i.e., independent of the initial values of the ambipolar electric field. It turns out that, in this case, two regimes are possible, depending on the values of the plasma parameters. In the first regime, the solution is unique and is independent of the initial conditions. In the second regime, two steady solutions can exist, depending on the initial conditions. One of the solution is similar to that obtained in the first regime, and the other differs from the first one both in the ambipolar field profile and in the dependence of the density and temperatures on the minor plasma radius. It cannot be excluded that different plasma confinement modes revealed in experiments are associated with the existence of such solutions.