Existence and stability of periodic travelling wave solutions to Nagumo's nerve equation

Summary Nagumo's nerve conduction equation has a one-parameter family of spatially periodic travelling wave solutions. First, we prove the existence of these solutions by using a topological method. (There are some exceptional cases in which this method cannot be applied in showing the existence.) A periodic travelling wave solution corresponds to a closed orbit of a third-order dynamical system. The Poincaré index of the closed orbit is determined as a direct consequence of the proof of the existence. Second, we prove that the periodic travelling wave solution is unstable if the Poincaré index of the corresponding closed orbit is + 1. By using this result, together with the result of the author's previous paper, it is concluded that the slow periodic travelling wave solutions are always unstable. Third, we consider the stability of the fast periodic travelling wave solutions. We denote by L(c) the spatial period of the travelling wave solution with the propagation speed c. It is shown that the fast solution is unstable if its period is close to L_min, the minimum of L(c).     

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.     

Travelling wave phenomena in non-linear diffusion degenerate Nagumo equations

 In this paper we study the existence of one-dimensional travelling wave solutions u(x, t)=φ(x−ct) for the non-linear degenerate (at u=0) reaction-diffusion equation u _t =[D(u)u _x ] _x +g(u) where g is a generalisation of the Nagumo equation arising in nerve conduction theory, as well as describing the Allee effect. We use a dynamical systems approach to prove: 1. the global bifurcation of a heteroclinic cycle (two monotone stationary front solutions), for c=0, 2. The existence of a unique value c ^*>0 of c for which φ(x−c ^* t) is a travelling wave solution of sharp type and 3. A continuum of monotone and oscillatory fronts for c≠c ^*. We present some numerical simulations of the phase portrait in travelling wave coordinates and on the full partial differential equation. Received 15 December 1995; received in revised form 14 May 1996     

Calcium wave propagation by calcium-induced calcium release: An unusual excitable system

We discuss in detail the behaviour of a model, proposed by Goldbeteret al. (1990.Proc. natn. Acad. Sci. 87, 1461–1465), for intracellular calcium wave propagation by calcium-induced calcium release, focusing our attention on excitability and the propagation of waves in one spatial dimension. The model with no diffusion behaves like a generic excitable system, and threshold behaviour, excitability and oscillations can be understood within this general framework. However, when diffusion is included, the model no longer behaves like a generic excitable system; the fast and slow variables are not distinct and previous results on excitable systems do not necessarily apply. We consider a piecewise linear simplification of the model, and construct travelling pulse and periodic plane wave solutions to the simplified model. The analogous behaviour in the full model is studied numerically. Goldbeter's model for calciuminduced calcium release is an excitable system of a type not previously studied in detail.     

Interaction of peristaltic flow with pulsatile flow in a circular cylindrical tube

The effect of pulsatile flow on peristaltic transport in a circular cylindrical tube is analysed. The flow of a Newtonian viscous incompressible fluid in a flexible circular cylindrical tube on which an axisymmetric travelling sinusoidal wave is imposed, is considered. The initial flow in the tube is induced by an arbitrary periodic pressure gradient. A perturbation solution with amplitude ratio (wave amplitude/tube radius) as a parameter is obtained when the frequency of the travelling wave and that of the imposed pressure gradient are equal. The interaction effects of periodic wall induced flow and periodic pressure imposed flow are visualized through the presence of substantially different components of steady and higher harmonic oscillating flow in the first order flow solution. Numerical results show a strong variation of steady state velocity profiles with boundary wave number and Reynolds number and a strong phase shift behaviour of the flow in the radial direction.     

Early development and quorum sensing in bacterial biofilms

 We develop mathematical models to examine the formation, growth and quorum sensing activity of bacterial biofilms. The growth aspects of the model are based on the assumption of a continuum of bacterial cells whose growth generates movement, within the developing biofilm, described by a velocity field. A model proposed in Ward et al. (2001) to describe quorum sensing, a process by which bacteria monitor their own population density by the use of quorum sensing molecules (QSMs), is coupled with the growth model. The resulting system of nonlinear partial differential equations is solved numerically, revealing results which are qualitatively consistent with experimental ones. Analytical solutions derived by assuming uniform initial conditions demonstrate that, for large time, a biofilm grows algebraically with time; criteria for linear growth of the biofilm biomass, consistent with experimental data, are established. The analysis reveals, for a biologically realistic limit, the existence of a bifurcation between non-active and active quorum sensing in the biofilm. The model also predicts that travelling waves of quorum sensing behaviour can occur within a certain time frame; while the travelling wave analysis reveals a range of possible travelling wave speeds, numerical solutions suggest that the minimum wave speed, determined by linearisation, is realised for a wide class of initial conditions. Received: 10 February 2002 / Revised version: 29 October 2002 / Published online: 19 March 2003 Key words or phrases: Bacterial biofilm – Quorum sensing – Mathematical modelling – Numerical solution – Asymptotic analysis – Travelling wave analysis     

Effects of noise on models of spiny dendrites

We study the effects of noise in two models of spiny dendrites. Through the introduction of different types of noise to both the Spike-diffuse-spike (SDS) and Baer–Rinzel (BR) models we investigate the change in behaviour of the travelling wave solution present in both deterministic systems, as noise intensity increases. We show that the speed of wave propagation in both the SDS and BR models respectively differs as the noise intensity in the spine heads increases. In contrast the cable is very robust to noise and as such the speed shows very little variation from the deterministic system. We introduce a space-dependent spine density, ρ(x), to the original Baer–Rinzel model and show how this modified model can mimic behaviour (under influence of noise) of both original systems, through variation of one parameter. We also show that the correlation time and length scales of the noise can enhance propagation of travelling wave solutions where the white noise dominates the underlying signal and produces noise induced phenomena.     

Interaction of Peristaltic Flow with Pulsatile Couple-Stress Fluid

1IntroductionPeristalsisIsnow－wellknowntothephyslologlststobeoneofthem8JormechanismforfluidtransportInmanybiologicalsystems．Inpatlcular，peristaltlcmechanismmaybeInvolvedInswal－lowing恤throughtheesophagus，urinetransPOrtfromkidneytobladderthoughuner．Inaddl－tion，perlstaltlcpumpingoccursInmanypracticalapplicationsInvolvingbio－mechanicalsystems．Thestudyofthemechanismofperistalsis，Inbothmechanicalandphysiologicalsituations，hasre－centlybecometheoNectofs。；ent；f；crese。roh．S；nce…     

The competitive dynamics between tumor cells, a replication-competent virus and an immune response

Replication-competent viruses have been used as an alternative therapeutic approach for cancer treatment. However, new clinical data revealed an innate immune response to virus that may mitigate the effects of treatment. Recently, Wein, Wu and Kirn have established a model which describes the interaction between tumor cells, a replication-competent virus and an immune response (Cancer Research 63 (2003):1317–1324). The purpose of this paper is to extend their model from the viewpoints of mathematics and biology and then prove global existence and uniqueness of solution to this new model, to study the dynamics of this novel therapy for cancers, and to explore a explicit threshold of the intensity of the immune response for controlling the tumor. We also study a time-delayed version of the model. We analytically prove that there exists a critical value τ₀ of the time-delay τ such that the system has a periodic solution if τ>τ₀. Numerical simulations are given to verify the analytical results. Furthermore, we numerically study the spatio-temporal dynamics of the model. The effects of the diffusivity of the immune response on the tumor growth are also discussed.Mathematics Subject Classification (2000): 35R35, 92A15     

Lateralization of Travelling Wave Response in the Hearing Organ of Bushcrickets

Travelling waves are the physical basis of frequency discrimination in many vertebrate and invertebrate taxa, including mammals, birds, and some insects. In bushcrickets (Tettigoniidae), the crista acustica is the hearing organ that has been shown to use sound-induced travelling waves. Up to now, data on mechanical characteristics of sound-induced travelling waves were only available along the longitudinal (proximal-distal) direction. In this study, we use laser Doppler vibrometry to investigate in-vivo radial (anterior-posterior) features of travelling waves in the tropical bushcricket Mecopoda elongata. Our results demonstrate that the maximum of sound-induced travelling wave amplitude response is always shifted towards the anterior part of the crista acustica. This lateralization of the travelling wave response induces a tilt in the motion of the crista acustica, which presumably optimizes sensory transduction by exerting a shear motion on the sensory cilia in this hearing organ.     

Existence of traveling wave solutions in a diffusive predator-prey model

 We establish the existence of traveling front solutions and small amplitude traveling wave train solutions for a reaction-diffusion system based on a predator-prey model with Holling type-II functional response. The traveling front solutions are equivalent to heteroclinic orbits in R ⁴ and the small amplitude traveling wave train solutions are equivalent to small amplitude periodic orbits in R ⁴ . The methods used to prove the results are the shooting argument and the Hopf bifurcation theorem. Received: 25 May 2001 / Revised version: 5 August 2002 / Published online: 19 November 2002 RID="*" ID="*" Research was supported by the National Natural Science Foundations (NNSF) of China. RID="*" ID="*" Research was partially supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada. On leave from the Department of Mathematics and Statistics, Dalhousie University, Halifax, Nova Scotia B3H 3J5, Canada. Mathematics Subject Classification (2000): 34C35, 35K57 Key words or phrases: Traveling wave solution – Wazewski set – Shooting argument – Hopf bifurcation Acknowledgements. We would like to thank the two referees for their careful reading and helpful comments.     

Partition of the Hodgkin-Huxley type model parameter space into the regions of qualitatively different solutions

We have examined the problem of obtaining relationships between the type of stable solutions of the Hodgkin-Huxley type system, the values of its parameters and a constant applied current (I). As variable parameters of the system the maximal Na⁺(g_Na),K⁺(g_K) conductances and shifts (Gm, Gh, Gn) of the voltage-dependences have been chosen. To solve this problem it is sufficient to find points belonging to the boundary, partitioning the parameter space of the system into the regions of the qualitatively different types of stable solutions (steady states and stable periodic oscillations). Almost all over the physiological range of I, a type of stable solution is determined by the type of steady state (stable or unstable). Using this fact, the approximate solution of this problem could be obtained by analyzing the spectrum of eigenvalues of the Jacobian matrix for the linearized system. The families of the plane sections of the boundary have been constructed in the three-parameter spaces (I, g_Na,g_K), (I, Gm, Gh), (I, Gm, Gn).     

Criterion for stable reentry in a ring of cardiac tissue

We model electrical wave propagation in a ring of cardiac tissue using an mth-order difference equation, where m denotes the number of cells in the ring. Under physiologically reasonable assumptions, the difference equation has a unique equilibrium solution. Applying Jury’s stability test, we prove a theorem concerning the local asymptotic stability of this equilibrium solution. Our results yield conditions for sustained reentrant tachycardia, a type of cardiac arrhythmia.      

From periodic travelling waves to travelling fronts in the spike-diffuse-spike model of dendritic waves

In the vertebrate brain excitatory synaptic contacts typically occur on the tiny evaginations of neuron dendritic surface known as dendritic spines. There is clear evidence that spine heads are endowed with voltage-dependent excitable channels and that action potentials invade spines. Computational models are being increasingly used to gain insight into the functional significance of a spine with an excitable membrane. The spike-diffuse-spike (SDS) model is one such model that admits to a relatively straightforward mathematical analysis. In this paper we demonstrate that not only can the SDS model support solitary travelling pulses, already observed numerically in more detailed biophysical models, but that it has periodic travelling wave solutions. The exact mathematical treatment of periodic travelling waves in the SDS model is used, within a kinematic framework, to predict the existence of connections between two periodic spike trains of different interspike interval. The associated wave front in the sequence of interspike intervals travels with a constant velocity without degradation of shape, and might therefore be used for the robust encoding of information.     

A net reproductive number for periodic matrix models

We give a definition of a net reproductive number R ₀ for periodic matrix models of the type used to describe the dynamics of a structured population with periodic parameters. The definition is based on the familiar method of studying a periodic map by means of its (period-length) composite. This composite has an additive decomposition that permits a generalization of the Cushing–Zhou definition of R ₀ in the autonomous case. The value of R ₀ determines whether the population goes extinct (R ₀<1) or persists (R ₀>1). We discuss the biological interpretation of this definition and derive formulas for R ₀ for two cases: scalar periodic maps of arbitrary period and periodic Leslie models of period 2. We illustrate the use of the definition by means of several examples and by applications to case studies found in the literature. We also make some comparisons of this definition of R ₀ with another definition given recently by Bacaër.     

Modeling Spatial Spread of Infectious Diseases with a Fixed Latent Period in a Spatially Continuous Domain

In this paper, with the assumptions that an infectious disease in a population has a fixed latent period and the latent individuals of the population may diffuse, we formulate an SIR model with a simple demographic structure for the population living in a spatially continuous environment. The model is given by a system of reaction-diffusion equations with a discrete delay accounting for the latency and a spatially non-local term caused by the mobility of the individuals during the latent period. We address the existence, uniqueness, and positivity of solution to the initial-value problem for this type of system. Moreover, we investigate the traveling wave fronts of the system and obtain a critical value c ^* which is a lower bound for the wave speed of the traveling wave fronts. Although we can not prove that this value is exactly the minimal wave speed, numeric simulations seem to suggest that it is. Furthermore, the simulations on the PDE model also suggest that the spread speed of the disease indeed coincides with c ^*. We also discuss how the model parameters affect c ^*.     

Monotone travelling fronts in an age-structured reaction-diffusion model of a single species

 We consider a partially coupled diffusive population model in which the state variables represent the densities of the immature and mature population of a single species. The equation for the mature population can be considered on its own, and is a delay differential equation with a delay-dependent coefficient. For the case when the immatures are immobile, we prove that travelling wavefront solutions exist connecting the zero solution of the equation for the matures with the delay-dependent positive equilibrium state. As a perturbation of this case we then consider the case of low immature diffusivity showing that the travelling front solutions continue to persist. Our findings are contrasted with recent studies of the delayed Fisher equation. Travelling fronts of the latter are known to lose monotonicity for sufficiently large delays. In contrast, travelling fronts of our equation appear to remain monotone for all values of the delay. Received: 1 November 2001 / Revised version: 10 May 2002 / Published online: 23 August 2002 Mathematics Subject Classification (2000): 35K57, 92D25 Key words or phrases: Age-structure – Time-delay – Travelling Fronts – Reaction-diffusion     

An approximate solution to the periodic bidomain equations in one dimension

An approximate, computationally tractable solution is proposed for the potentials in the bidomain model with periodic intracellular junctions (the periodic bidomain model). This new approach is based on the one-dimensional rigorous spectral method described previously by Trayanova and Pilkington (IEEE Trans. Biomed. Eng., May 1993). The total solution to the one-dimensional periodic bidomain problem is decomposed in the spectral domain into solutions to (1) the single-fiber classical bidomain problem in which the intracellular conductivity value incorporates the average contribution from cytoplasm and junction and (2) the “junctional” potential problem due to the presence of junctions at discrete locations alone. Solving for the junctional term rigorously requires most of the numerical effort in the solution for the periodic bidomain potentials. Here the junctional potential is found approximately with little numerical effort. A comparison between the rigorous and the approximate solutions serves as a justification for the proposed approximate solution procedure. The procedure outlined in this paper is applicable to higher spatial dimensions where both tissue anisotropy and junctional inhomogeneities play a role in establishing the transmembrane potential distribution.