Wave trains in an excitable FitzHugh-Nagumo model: bistable dispersion relation and formation of isolas

Double-wave reentry in excitable media

A monolayer of chick embryo cardiac cells grown in an annular geometry supports two simultaneous reentrant excitation waves that circulate as a doublet. We propose a mechanism that can lead to such behavior. The velocity restitution gives the instantaneous velocity of a wave as a function of the time since the passage of the previous wave at a given point in space. Nonmonotonic restitution relationships will lead to situations in which various spacings between circulating waves are possible. In cardiology, the situation in which two waves travel in an anatomically defined circuit is referred to as double-wave reentry. Since double-wave reentry may arise as a consequence of pacing during cardiac arrhythmias, understanding the dynamic features of double-wave reentry may be helpful in understanding the physiological properties of cardiac tissue and in the design of therapy.

Origin of finite pulse trains: Homoclinic snaking in excitable media

Many physical, chemical, and biological systems exhibit traveling waves as a result of either an oscillatory instability or excitability. In the latter case a large multiplicity of stable spatially localized wavetrains consisting of different numbers of traveling pulses may be present. The existence of these states is related here to the presence of homoclinic snaking in the vicinity of a subcritical, finite wavenumber Hopf bifurcation. The pulses are organized in a slanted snaking structure resulting from the presence of a heteroclinic cycle between small and large amplitude traveling waves. Connections of this type require a multivalued dispersion relation. This dispersion relation is computed numerically and used to interpret the profile of the pulse group. The different spatially localized pulse trains can be accessed by appropriately customized initial stimuli, thereby blurring the traditional distinction between oscillatory and excitable systems. The results reveal a new class of phenomena relevant to spatiotemporal dynamics of excitable media, particularly in chemical and biological systems with multiple activators and inhibitors.

Travelling waves in a neural field model with refractoriness

At one level of abstraction neural tissue can be regarded as a medium for turning local synaptic activity into output signals that propagate over large distances via axons to generate further synaptic activity that can cause reverberant activity in networks that possess a mixture of excitatory and inhibitory connections. This output is often taken to be a firing rate, and the mathematical form for the evolution equation of activity depends upon a spatial convolution of this rate with a fixed anatomical connectivity pattern. Such formulations often neglect the metabolic processes that would ultimately limit synaptic activity. Here we reinstate such a process, in the spirit of an original prescription by Wilson and Cowan (Biophys J 12:1–24, 1972), using a term that multiplies the usual spatial convolution with a moving time average of local activity over some refractory time-scale. This modulation can substantially affect network behaviour, and in particular give rise to periodic travelling waves in a purely excitatory network (with exponentially decaying anatomical connectivity), which in the absence of refractoriness would only support travelling fronts. We construct these solutions numerically as stationary periodic solutions in a co-moving frame (of both an equivalent delay differential model as well as the original delay integro-differential model). Continuation methods are used to obtain the dispersion curve for periodic travelling waves (speed as a function of period), and found to be reminiscent of those for spatially extended models of excitable tissue. A kinematic analysis (based on the dispersion curve) predicts the onset of wave instabilities, which are confirmed numerically.

Drifting solitary waves in a reaction-diffusion medium with differential advection

A distinct propagation of solitary waves in the presence of autocatalysis, diffusion, and symmetry-breaking (differential) advection, is being studied. These pulses emerge at lower reaction rates of the autocatalytic activator, i.e., when the advective flow overcomes the fast excitation and induces a fluid type "drifting" behavior, making the phenomenon unique to reaction-diffusion-advection class systems. Using the spatial dynamics analysis of a canonical model, we present the properties and the organization of such drifting pulses. The insights underly a general understanding of localized transport in simple reaction-diffusion-advection models and thus provide a background to potential chemical and biological applications.

Feedback control of subcritical Turing instability with zero mode

A global feedback control of a system that exhibits a subcritical monotonic instability at a nonzero wave number (short-wave or Turing instability) in the presence of a zero mode is investigated using a Ginzburg-Landau equation coupled to an equation for the zero mode. This system is studied analytically and numerically. It is shown that feedback control, based on measuring the maximum of the pattern amplitude over the domain, can stabilize the system and lead to the formation of localized unipulse stationary states or traveling solitary waves. It is found that the unipulse traveling structures result from an instability of the stationary unipulse structures when one of the parameters characterizing the coupling between the periodic pattern and the zero mode exceeds a critical value that is determined by the zero mode damping coefficient.

Generation of finite wave trains in excitable media

Spatiotemporal control of excitable media is of paramount importance in the development of new applications, ranging from biology to physics. To this end, we identify and describe a qualitative property of excitable media that enables us to generate a sequence of traveling pulses of any desired length, using a one-time initial stimulus. The wave trains are produced by a transient pacemaker generated by a one-time suitably tailored spatially localized finite amplitude stimulus, and belong to a family of fast pulse trains. A second family, of slow pulse trains, is also present. The latter are created through a clumping instability of a traveling wave state (in an excitable regime) and are inaccessible to single localized stimuli of the type we use. The results indicate that the presence of a large multiplicity of stable, accessible, multi-pulse states is a general property of simple models of excitable media.

Anomalous pulse interaction in dissipative media

We review a number of phenomena occurring in one-dimensional excitable media due to modified decay behind propagating pulses. Those phenomena can be grouped in two categories depending on whether the wake of a solitary pulse is oscillatory or not. Oscillatory decay leads to nonannihilative head-on collision of pulses and oscillatory dispersion relation of periodic pulse trains. Stronger wake oscillations can even result in a bistable dispersion relation. Those effects are illustrated with the help of the Oregonator and FitzHugh-Nagumo models for excitable media. For a monotonic wake, we show that it is possible to induce bound states of solitary pulses and anomalous dispersion of periodic pulse trains by introducing nonlocal spatial coupling to the excitable medium.

Wave-pinned filaments of scroll waves

Scroll waves are three-dimensional excitation patterns that rotate around one-dimensional space curves. Typically these filaments are closed loops or end at the system boundary. However, in excitable media with anomalous dispersion, filaments can be pinned to the wake of traveling wave pulses. This pinning is studied in experiments with the 1,4-cyclohexanedione Belousov-Zhabotinsky reaction and a three-variable reaction-diffusion model. We show that wave-pinned filaments are related to the coexistence of rotating and translating wave defects in two dimensions. Filament pinning causes a continuous expansion of the total filament length. It can be ended by annihilating the pinning pulse in a frontal wave collision. Following such an annihilation, the filament connects itself to the system boundary. Its postannihilation shape that is initially the exposed rim of the scroll wave unwinds continuously over numerous rotation periods.