Creating bound states in excitable media by means of nonlocal coupling

Invariant manifolds and global bifurcations

Invariant manifolds are key objects in describing how trajectories partition the phase spaces of a dynamical system. Examples include stable, unstable, and center manifolds of equilibria and periodic orbits, quasiperiodic invariant tori, and slow manifolds of systems with multiple timescales. Changes in these objects and their intersections with variation of system parameters give rise to global bifurcations. Bifurcation manifolds in the parameter spaces of multi-parameter families of dynamical systems also play a prominent role in dynamical systems theory. Much progress has been made in developing theory and computational methods for invariant manifolds during the past 25 years. This article highlights some of these achievements and remaining open problems.

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.

Analyzing critical propagation in a reaction-diffusion-advection model using unstable slow waves

The effect of advection on the propagation and in particular on the critical minimal speed of traveling waves in a reaction-diffusion model is studied. Previous theoretical studies estimated this effect on the velocity of stable fast waves and predicted the existence of a critical advection strength below which propagating waves are not supported anymore. In this paper, an analytical expression for the advection-velocity relation of the unstable slow wave is derived. In addition, the critical advection strength is calculated taking into account the unstable slow wave solution. We also analyze a two-variable reaction-diffusion-advection model numerically in a wide parameter range. Due to the new control parameter (advection) we can find stable wave propagation in the otherwise non-excitable parameter regime, if the advection strength exceeds a critical value. Comparing theoretical predictions to numerical results, we find that they are in good agreement. Theory provides an explanation for the observed behaviour.

Dynamics of reaction-diffusion patterns controlled by asymmetric nonlocal coupling as a limiting case of differential advection

A one-component bistable reaction-diffusion system with asymmetric nonlocal coupling is derived as a limiting case of a two-component activator-inhibitor reaction-diffusion model with differential advection. The effects of asymmetric nonlocal couplings in such a bistable reaction-diffusion system are then compared to the previously studied case of a system with symmetric nonlocal coupling. We carry out a linear stability analysis of the spatially homogeneous steady states of the model and numerical simulations of the model to show how the asymmetric nonlocal coupling controls and alters the steady states and the front dynamics in the system. In a second step, a third fast reaction-diffusion equation is included which induces the formation of more complex patterns. A linear stability analysis predicts traveling waves for asymmetric nonlocal coupling, in contrast to a stationary Turing patterns for a system with symmetric nonlocal coupling. These findings are verified by direct numerical integration of the full equations with nonlocal coupling.

Formation of localized structures in bistable systems through nonlocal spatial coupling. I. General framework

The present work studies the influence of nonlocal spatial coupling on the existence of localized structures in one-dimensional extended systems. We consider systems described by a real field with a nonlocal coupling that has a linear dependence on the field. Leveraging spatial dynamics we provide a general framework to understand the effect of the nonlocality on the shape of the fronts connecting two stable states. In particular we show that nonlocal terms can induce spatial oscillations in the front tails, allowing for the creation of localized structures, that emerge from pinning between two fronts. In parameter space the region where fronts are oscillatory is limited by three transitions: the modulational instability of the homogeneous state, the Belyakov-Devaney transition in which monotonic fronts acquire spatial oscillations with infinite wavelength, and a crossover in which monotonically decaying fronts develop spatial oscillations with a finite wavelength. We show how these transitions are organized by codimension 2 and 3 points and illustrate how by changing the parameters of the nonlocal coupling it is possible to bring the system into the region where localized structures can be formed.

Formation of localized structures in bistable systems through nonlocal spatial coupling. II. The nonlocal Ginzburg-Landau equation

We study the influence of a linear nonlocal spatial coupling on the interaction of fronts connecting two equivalent stable states in the prototypical 1-dimensional real Ginzburg-Landau equation. While for local coupling the fronts are always monotonic and therefore the dynamical behavior leads to coarsening and the annihilation of pairs of fronts, nonlocal terms can induce spatial oscillations in the front, allowing for the creation of localized structures, emerging from pinning between two fronts. We show this for three different nonlocal influence kernels. The first two, mod-exponential and Gaussian, are positive definite and decay exponentially or faster, while the third one, a Mexican-hat kernel, is not positive definite.

Self-organized pacemakers and bistability of pulses in an excitable medium

Pattern formation in an excitable medium described by a three-component reaction-diffusion system is investigated. Our focus is on stable self-organized pacemakers which give rise to spatially extended target patterns. Bistability of pulse solutions in the excitable regime is also reported, and interactions of the different pulses with each other and the pacemaker are studied. Self-organized pacemakers are created by a suitable perturbation from the steady state or through interaction of pulses. Bound states of one-dimensional pacemakers and phase flips are also observed.

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.