Transition from spirals to defect turbulence driven by a convective instability

Annihilation dynamics during spiral defect chaos revealed by particle models

Pair-annihilation events are ubiquitous in a variety of spatially extended systems and are often studied using computationally expensive simulations. Here, we develop an approach in which we simulate the pair-annihilation of spiral wave tips in cardiac models using a computationally efficient particle model. Spiral wave tips are represented as particles with dynamics governed by diffusive behavior and short-ranged attraction. The parameters for diffusion and attraction are obtained by comparing particle motion to the trajectories of spiral wave tips in cardiac models during spiral defect chaos. The particle model reproduces the annihilation rates of the cardiac models and can determine the statistics of spiral wave dynamics, including its mean termination time. We show that increasing the attraction coefficient sharply decreases the mean termination time, making it a possible target for pharmaceutical intervention.

Pulsating Active Matter

We reveal that the mechanical pulsation of locally synchronized particles is a generic route to propagate deformation waves. We consider a model of dense repulsive particles whose activity drives periodic change in size of each individual. The dynamics is inspired by biological tissues where cells consume fuel to sustain active deformation. We show that the competition between repulsion and synchronization triggers an instability which promotes a wealth of dynamical patterns, ranging from spiral waves to defect turbulence. We identify the mechanisms underlying the emergence of patterns, and characterize the corresponding transitions. By coarse-graining the dynamics, we propose a hydrodynamic description of an assembly of pulsating particles, and discuss an analogy with reaction-diffusion systems.

Composite spiral waves in discrete-time systems

Spiral waves are a type of typical pattern in open reaction-diffusion systems far from thermodynamic equilibrium. The study of spiral waves has attracted great interest because of its nonlinear characteristics and extensive applications. However, the study of spiral waves has been confined to continuous-time systems, while spiral waves in discrete-time systems have been rarely reported. In recent years, discrete-time models have been widely studied in ecology because of their appropriateness to systems with nonoverlapping generations and other factors. Therefore, spiral waves in discrete-time systems need to be studied. Here, we investigated a novel type of spiral wave called a composite spiral wave in a discrete-time predator-pest model, and we revealed the formation mechanism. To explain the observed phenomena, we defined and quantified a move state effect of multiperiod states caused by the coupling of adjacent stable multiperiod orbits, which is strictly consistent with the numerical results. The other move state effect is caused by an unstable focus, which is the state of the local points at the spiral center. The combined effect of these two influences can lead to rich dynamical behaviors of spiral waves, and the specific structure of the composite spiral waves is the result of the competition of the two effects in opposite directions. Our findings shed light on the dynamics of spiral waves in discrete-time systems, and they may guide the prediction and control of pests in deciduous forests.

Markov modeling of phase singularity interaction effects in human atrial and ventricular fibrillation

Atrial and ventricular fibrillation (AF/VF) are characterized by the repetitive regeneration of topological defects known as phase singularities (PSs). The effect of PS interactions has not been previously studied in human AF and VF. We hypothesized that PS population size would influence the rate of PS formation and destruction in human AF and VF, due to increased inter-defect interaction. PS population statistics were studied in computational simulations (Aliev-Panfilov), human AF and human VF. The influence of inter-PS interactions was evaluated by comparison between directly modeled discrete-time Markov chain (DTMC) transition matrices of the PS population changes, and M/M/∞ birth-death transition matrices of PS dynamics, which assumes that PS formations and destructions are effectively statistically independent events. Across all systems examined, PS population changes differed from those expected with M/M/∞. In human AF and VF, the formation rates decreased slightly with PS population when modeled with the DTMC, compared with the static formation rate expected through M/M/∞, suggesting new formations were being inhibited. In human AF and VF, the destruction rates increased with PS population for both models, with the DTMC rate increase exceeding the M/M/∞ estimates, indicating that PS were being destroyed faster as the PS population grew. In human AF and VF, the change in PS formation and destruction rates as the population increased differed between the two models. This indicates that the presence of additional PS influenced the likelihood of new PS formation and destruction, consistent with the notion of self-inhibitory inter-PS interactions.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Control of spiral waves in optogenetically modified cardiac tissue by periodic optical stimulation

Resonant drift of nonlinear spiral waves occurs in various types of excitable media under periodic stimulation. Recently a novel methodology of optogenetics has emerged, which allows to affect excitability of cardiac cells and tissues by optical stimuli. In this paper we study if resonant drift of spiral waves in the heart can be induced by a homogeneous weak periodic optical stimulation of cardiac tissue. We use a two-variable and a detailed model of cardiac tissue and add description of depolarizing and hyperpolarizing optogenetic ionic currents. We show that weak periodic optical stimulation at a frequency equal to the natural rotation frequency of the spiral wave induces resonant drift for both depolarizing and hyperpolarizing optogenetic currents. We quantify these effects and study how the speed of the drift and its direction depend on the initial conditions. We also derive analytical formulas based on the response function theory which correctly predict the drift velocity and its trajectory. We conclude that optogenetic methodology can be used for control of spiral waves in cardiac tissue and discuss its possible applications.

The inspection paradox: An important consideration in the evaluation of rotor lifetimes in cardiac fibrillation

Background and Objective: Renewal theory is a statistical approach to model the formation and destruction of phase singularities (PS), which occur at the pivots of spiral waves. A common issue arising during observation of renewal processes is an inspection paradox, due to oversampling of longer events. The objective of this study was to characterise the effect of a potential inspection paradox on the perception of PS lifetimes in cardiac fibrillation. Methods: A multisystem, multi-modality study was performed, examining computational simulations (Aliev-Panfilov (APV) model, Courtmanche-Nattel model), experimentally acquired optical mapping Atrial and Ventricular Fibrillation (AF/VF) data, and clinically acquired human AF and VF. Distributions of all PS lifetimes across full epochs of AF, VF, or computational simulations, were compared with distributions formed from lifetimes of PS existing at 10,000 simulated commencement timepoints. Results: In all systems, an inspection paradox led towards oversampling of PS with longer lifetimes. In APV computational simulations there was a mean PS lifetime shift of +84.9% (95% CI, ± 0.3%) (p < 0.001 for observed vs overall), in Courtmanche-Nattel simulations of AF +692.9% (95% CI, ±57.7%) (p < 0.001), in optically mapped rat AF +374.6% (95% CI, ± 88.5%) (p = 0.052), in human AF mapped with basket catheters +129.2% (95% CI, ±4.1%) (p < 0.05), human AF-HD grid catheters 150.8% (95% CI, ± 9.0%) (p < 0.001), in optically mapped rat VF +171.3% (95% CI, ±15.6%) (p < 0.001), in human epicardial VF 153.5% (95% CI, ±15.7%) (p < 0.001). Conclusion: Visual inspection of phase movies has the potential to systematically oversample longer lasting PS, due to an inspection paradox. An inspection paradox is minimised by consideration of the overall distribution of PS lifetimes.

Bulk-surface coupling identifies the mechanistic connection between Min-protein patterns in vivo and in vitro

Self-organisation of Min proteins is responsible for the spatial control of cell division in Escherichia coli, and has been studied both in vivo and in vitro. Intriguingly, the protein patterns observed in these settings differ qualitatively and quantitatively. This puzzling dichotomy has not been resolved to date. Using reconstituted proteins in laterally wide microchambers with a well-controlled height, we experimentally show that the Min protein dynamics on the membrane crucially depend on the micro chamber height due to bulk concentration gradients orthogonal to the membrane. A theoretical analysis shows that in vitro patterns at low microchamber height are driven by the same lateral oscillation mode as pole-to-pole oscillations in vivo. At larger microchamber height, additional vertical oscillation modes set in, marking the transition to a qualitatively different in vitro regime. Our work reveals the qualitatively different mechanisms of mass transport that govern Min protein-patterns for different bulk heights and thus shows that Min patterns in cells are governed by a different mechanism than those in vitro.

Arnold tongues in oscillator systems with nonuniform spatial driving

Nonlinear oscillator systems are ubiquitous in biology and physics, and their control is a practical problem in many experimental systems. Here we study this problem in the context of the two models of spatially coupled oscillators: the complex Ginzburg-Landau equation (CGLE) and a generalization of the CGLE in which oscillators are coupled through an external medium (emCGLE). We focus on external control drives that vary in both space and time. We find that the spatial distribution of the drive signal controls the frequency ranges over which oscillators synchronize to the drive and that boundary conditions strongly influence synchronization to external drives for the CGLE. Our calculations also show that the emCGLE has a low density regime in which a broad range of frequencies can be synchronized for low drive amplitudes. We study the bifurcation structure of these models and find that they are very similar to results for the driven Kuramoto model, a system with no spatial structure. We conclude by discussing qualitative implications of our results for controlling coupled oscillator systems such as the social amoebae Dictyostelium and populations of Belousov Zhabotinsky (BZ) catalytic particles using spatially structured external drives.

Turbulent pattern in the 1,4-cyclohexanedione Belousov-Zhabotinsky reaction

Chemical turbulence was observed experimentally in the 1,4-cyclohexanedione Belousov-Zhabotinsky (CHD-BZ) reaction in a double layer consisting of a catalyst-loaded gel and uncatalyzed liquid on a Petri dish. The chemical patterns in the CHD-BZ reaction occur spontaneously in various forms as follows: the initial, regular, transient, and turbulent patterns, subsequently. These four patterns are characterized by using the two-dimensional Fourier transform (2D-FT). Mechanism of the onset of the turbulence in the CHD-BZ reaction is proposed. Turbulence in the CHD-BZ reaction is reproducible under a well defined protocol and it exists for a period of time of about 50 minutes, which is sufficiently long to offer a good opportunity to study and control the turbulence in the future. Two models of the BZ reaction were used to simulate the spiral breakup. Both are capable of producing spiral turbulence from initially regular patterns in each layer and reflect certain features of dynamics observed in experiments.

Single to multiple acoustic vortex excitations in the transition to defect-mediated dust acoustic wave turbulence

We experimentally investigate the cooperative excitations in the transition from a self-excited three-dimensional ordered plane wave to a defect-mediated turbulence (DMT) state with multiple unstable defect filaments in a dusty plasma system. It is found that, with increasing effective driving, a single acoustic vortex (AV) with positive or negative helicity winding around a long straight defect filament with small wiggling in the 2+1D (dimensional) space-time space starts to emerge along the center axis of the small dust cluster. The sequential ruptures of the crest surfaces from the cluster boundary followed by their reconnection with adjacent ruptured crest surfaces, or repelling one of the pairwise generated defects out of the boundary, are the key for the single AV generation. Further increasing driving makes the single defect filament exhibit helical excursion in the 2+1D space. The system eventually enters the state with a few short-lived AVs and the DMT state with multiple AVs. The gradual increasing defect filament fluctuations and defect number in the transition to the DMT more strongly distort the nearby waveforms, which leads to the transition from the emergence of distinct sideband peaks to the broadened peaks in the power spectra of temporal dust density fluctuation. For the system with a larger cluster size, the single AV states are skipped in the transition to the DMT state.

Extinction dynamics of spiral defect chaos

Spatially extended excitable systems can exhibit spiral defect chaos (SDC) during which spiral waves continuously form and disappear. To address how this dynamical state terminates using simulations can be computationally challenging, especially for large systems. To circumvent this limitation, we treat the number of spiral waves as a stochastic population with a corresponding birth-death equation and use techniques from statistical physics to determine the mean episode duration of SDC. Motivated by cardiac fibrillation, during which the heart's electrical activity becomes disorganized and shows fragmenting spiral waves, we use generic models of cardiac electrophysiology. We show that the duration can be computed in minimal computational time and that it depends exponentially on domain size. Therefore, the approach can result in efficient and accurate predictions of mean episode duration which may be extended to more complex geometries and models.

Spiral wave initiation in excitable media

Spiral waves represent an important example of dissipative structures observed in many distributed systems in chemistry, biology and physics. By definition, excitable media occupy a stationary resting state in the absence of external perturbations. However, a perturbation exceeding a threshold results in the initiation of an excitation wave propagating through the medium. These waves, in contrast to acoustic and optical ones, disappear at the medium's boundary or after a mutual collision, and the medium returns to the resting state. Nevertheless, an initiation of a rotating spiral wave results in a self-sustained activity. Such activity unexpectedly appearing in cardiac or neuronal tissues usually destroys their dynamics which results in life-threatening diseases. In this context, an understanding of possible scenarios of spiral wave initiation is of great theoretical importance with many practical applications.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 2)'.

Spiral wave initiation in excitable media

Spiral waves represent an important example of dissipative structures observed in many distributed systems in chemistry, biology and physics. By definition, excitable media occupy a stationary resting state in the absence of external perturbations. However, a perturbation exceeding a threshold results in the initiation of an excitation wave propagating through the medium. These waves, in contrast to acoustic and optical ones, disappear at the medium's boundary or after a mutual collision, and the medium returns to the resting state. Nevertheless, an initiation of a rotating spiral wave results in a self-sustained activity. Such activity unexpectedly appearing in cardiac or neuronal tissues usually destroys their dynamics which results in life-threatening diseases. In this context, an understanding of possible scenarios of spiral wave initiation is of great theoretical importance with many practical applications.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 2)'.

Chimera States in Continuous Media: Existence and Distinctness

The defining property of chimera states is the coexistence of coherent and incoherent domains in systems that are structurally and spatially homogeneous. The recent realization that such states might be common in oscillator networks raises the question of whether an analogous phenomenon can occur in continuous media. Here, we show that chimera states can exist in continuous systems even when the coupling is strictly local, as in many fluid and pattern forming media. Using the complex Ginzburg-Landau equation as a model system, we characterize chimera states consisting of a coherent domain of a frozen spiral structure and an incoherent domain of amplitude turbulence. We show that in this case, in contrast with discrete network systems, fluctuations in the local coupling field play a crucial role in limiting the coherent regions. We suggest these findings shed light on new possible forms of coexisting order and disorder in fluid systems.

Dynamics of spiral waves rotating around an obstacle and the existence of a minimal obstacle

Pinning of vortices by obstacles plays an important role in various systems. In the heart, anatomical reentry is created when a vortex, also known as the spiral wave, is pinned to an anatomical obstacle, leading to a class of physiologically very important arrhythmias. Previous analyses of its dynamics and instability provide fine estimates in some special circumstances, such as large obstacles or weak excitabilities. Here, to expand theoretical analyses to all circumstances, we propose a general theory whose results quantitatively agree with direct numerical simulations. In particular, when obstacles are small and pinned spiral waves are destabilized, an accurate explanation of the instability in two-dimensional media is provided by the usage of a mapping rule and dimension reduction. The implications of our results are to better understand the mechanism of arrhythmia and thus improve its early prevention.

Weakly and strongly coupled Belousov-Zhabotinsky patterns

We investigate experimentally and numerically the synchronization of two-dimensional spiral wave patterns in the Belousov-Zhabotinsky reaction due to point-to-point coupling of two separate domains. Different synchronization modalities appear depending on the coupling strength and the initial patterns in each domain. The behavior as a function of the coupling strength falls into two qualitatively different regimes. The weakly coupled regime is characterized by inter-domain interactions that distorted but do not break wave fronts. Under weak coupling, spiral cores are pushed around by wave fronts in the other domain, resulting in an effective interaction between cores in opposite domains. In the case where each domain initially contains a single spiral, the cores form a bound pair and orbit each other at quantized distances. When the starting patterns consist of multiple randomly positioned spiral cores, the number of cores decreases with time until all that remains are a few cores that are synchronized with a partner in the other domain. The strongly coupled regime is characterized by interdomain interactions that break wave fronts. As a result, the wave patterns in both domains become identical.

Spatiotemporal chaos involving wave instability

In this paper, we investigate pattern formation in a model of a reaction confined in a microemulsion, in a regime where both Turing and wave instability occur. In one-dimensional systems, the pattern corresponds to spatiotemporal intermittency where the behavior of the systems alternates in both time and space between stationary Turing patterns and traveling waves. In two-dimensional systems, the behavior initially may correspond to Turing patterns, which then turn into wave patterns. The resulting pattern also corresponds to a chaotic state, where the system alternates in both space and time between standing wave patterns and traveling waves, and the local dynamics may show vanishing amplitude of the variables.

Pattern transitions in spatial epidemics: Mechanisms and emergent properties

Infectious diseases are a threat to human health and a hindrance to societal development. Consequently, the spread of diseases in both time and space has been widely studied, revealing the different types of spatial patterns. Transitions between patterns are an emergent property in spatial epidemics that can serve as a potential trend indicator of disease spread. Despite the usefulness of such an indicator, attempts to systematize the topic of pattern transitions have been few and far between. We present a mini-review on pattern transitions in spatial epidemics, describing the types of transitions and their underlying mechanisms. We show that pattern transitions relate to the complexity of spatial epidemics by, for example, being accompanied with phenomena such as coherence resonance and cyclic evolution. The results presented herein provide valuable insights into disease prevention and control, and may even be applicable outside epidemiology, including other branches of medical science, ecology, quantitative finance, and elsewhere.

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that is triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media with applications to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact for cardiac arrhythmias.

Optical textures: characterizing spatiotemporal chaos

Macroscopic systems subjected to injection and dissipation of energy can exhibit complex spatiotemporal behaviors as result of dissipative self-organization. Here, we report a one- and two-dimensional pattern forming setup, which exhibits a transition from stationary patterns to spatiotemporal chaotic textures, based on a nematic liquid crystal layer with spatially modulated input beam and optical feedback. Using an adequate projection of spatiotemporal diagrams, we determine the largest Lyapunov exponent. Jointly, this exponent and Fourier transform allow us to distinguish between spatiotemporal chaos and amplitude turbulence concepts, which are usually merged.

Optical Rogue Waves in Vortex Turbulence

We present a spatiotemporal mechanism for producing 2D optical rogue waves in the presence of a turbulent state with creation, interaction, and annihilation of optical vortices. Spatially periodic structures with bound phase lose stability to phase unbound turbulent states in complex Ginzburg-Landau and Swift-Hohenberg models with external driving. When the pumping is high and the external driving is low, synchronized oscillations are unstable and lead to spatiotemporal vortex-mediated turbulence with high excursions in amplitude. Nonlinear amplification leads to rogue waves close to turbulent optical vortices, where the amplitude tends to zero, and to probability density functions (PDFs) with long tails typical of extreme optical events.

Propagation Properties of the Precipitation Band in an AlCl₃/NaOH System

When inherently immobile solid particles collectively form precipitates in a reaction-diffusion system involving a redissolution reaction, a propagation phenomenon may occur in which a dynamic pattern of precipitation bands forms. This propagating precipitation phenomenon has been studied by many researchers. However, two completely different processes-i.e., the reaction-diffusion of reactants and the crystal growth of products-progress simultaneously in the system, thereby rendering the phenomenon complex. There are no well-established experimental laws for this propagating precipitation phenomenon, such as the spacing, time, and width laws associated with the well-known Liesegang phenomenon, which is static in the sense that precipitation bands form and remain at the same position. In fact, it has not been clarified which of the processes controls the propagation phenomenon. Accordingly, we have investigated the apparent diffusion coefficient associated with the dynamics of propagating precipitation band in an AlCl3/NaOH system for the case in which a large excess of outer electrolytes (i.e., OH(-)) diffuses into gel in which inner electrolytes (i.e.,Al(3+)) are homogeneously distributed. An isolated precipitation band of Al(OH)3 was formed horizontally in a test tube and propagated vertically in proportion to the square root of time. In our experimental results, we found that the apparent diffusion coefficient, D(p), possesses an exponential dependence on the initial concentrations of the outer electrolyte, and the inner electrolyte; the measured relation was D(p) = D[Al(3+)](-0.6)[OH(-)](0.6), where D = (0.63 ± 0.04) × 10(5) cm(2)/s. From our model equations based on the prenucleation theory, which take into account a redissolution reaction, we found that the dynamics of the reaction front of the outer and the inner electrolytes was an important factor in controlling the propagation of the precipitation band. In our simulation results, we obtained a similar dependence of the apparent diffusion coefficient on the electrolyte concentrations.

Standing wave-like patterns in the Gray-Scott model

Standing wave-like patterns are obtained in the Gray-Scott model when the dynamics that correspond to defect-mediated turbulence for equal diffusivities interact with a Turing instability. The Turing instability can be caused by either differential or cross-diffusion. We compare results with the Oregonator model, for which standing wave-like patterns are also observed under similar conditions.

Chiral selection and frequency response of spiral waves in reaction-diffusion systems under a chiral electric field

Chirality is one of the most fundamental properties of many physical, chemical, and biological systems. However, the mechanisms underlying the onset and control of chiral symmetry are largely understudied. We investigate possibility of chirality control in a chemical excitable system (the Belousov-Zhabotinsky reaction) by application of a chiral (rotating) electric field using the Oregonator model. We find that unlike previous findings, we can achieve the chirality control not only in the field rotation direction, but also opposite to it, depending on the field rotation frequency. To unravel the mechanism, we further develop a comprehensive theory of frequency synchronization based on the response function approach. We find that this problem can be described by the Adler equation and show phase-locking phenomena, known as the Arnold tongue. Our theoretical predictions are in good quantitative agreement with the numerical simulations and provide a solid basis for chirality control in excitable media.

Defect-mediated turbulence and transition to spatiotemporal intermittency in the Gray-Scott model

In this paper, we show that the Gray-Scott model is able to produce defect-mediated turbulence. This regime emerges from the limit cycle, close or far from the Hopf bifurcation, but always right before the Andronov homoclinic bifurcation of the homogeneous system. After this bifurcation, as the control parameter is further changed, the system starts visiting more and more frequently the stable node of the model. Consequently, the defect-mediated turbulence gradually turns into spatiotemporal intermittency.

Phase-bistable pattern formation in oscillatory systems via rocking: application to nonlinear optical systems

We present a review, together with new results, of a universal forcing of oscillatory systems, termed 'rocking', which leads to the emergence of a phase bistability and to the kind of pattern formation associated with it, characterized by the presence of phase domains, phase spatial solitons and phase-bistable extended patterns. The effects of rocking are thus similar to those observed in the classic 2 : 1 resonance (the parametric resonance) of spatially extended systems of oscillators, which occurs under a spatially uniform, time-periodic forcing at twice the oscillations' frequency. The rocking, however, has a frequency close to that of the oscillations (it is a 1 : 1 resonant forcing) and hence is a good alternative to the parametric forcing when the latter is inefficient (e.g. in optics). The key ingredient is that the rocking amplitude is modulated either in time or in space, such that its sign alternates (exhibits π-phase jumps). We present new results concerning a paradigmatic nonlinear optical system (the two-level laser) and show that phase domains and dark-ring (phase) solitons replace the ubiquitous vortices that characterize the emission of free-running, broad area lasers.

Transient spatiotemporal chaos in the Morris-Lecar neuronal ring network

Transient behavior is thought to play an integral role in brain functionality. Numerical simulations of the firing activity of diffusively coupled, excitable Morris-Lecar neurons reveal transient spatiotemporal chaos in the parameter regime below the saddle-node on invariant circle bifurcation point. The neighborhood of the chaotic saddle is reached through perturbations of the rest state, in which few initially active neurons at an effective spatial distance can initiate spatiotemporal chaos. The system escapes from the neighborhood of the chaotic saddle to either the rest state or to a state of pulse propagation. The lifetime of the chaotic transients is manipulated in a statistical sense through a singular application of a synchronous perturbation to a group of neurons.

Core solutions of rigidly rotating spiral waves in highly excitable media

Analytical spiral wave solutions for reaction-diffusion equations play an important role in studying spiral wave dynamics. In this paper, we focus on such analytical solutions in the case of highly excitable media. We present numerical evidence that, for rigidly rotating spiral waves in highly excitable media, the species values in the spiral core region do harmonic oscillations but not relaxation ones, and their amplitudes grow linearly with the distance from the rotation center. An analytical solution is proposed to describe such spiral wave dynamics, and the quantitative comparisons between the numerical results and the analytical solutions show that the proposed spiral core solution works well in highly excitable media.

Quasiperiodicity route to spatiotemporal chaos in one-dimensional pattern-forming systems

We propose a route to spatiotemporal chaos for one-dimensional stationary patterns, which is a natural extension of the quasiperiodicity route for low-dimensional chaos to extended systems. This route is studied through a universal model of pattern formation. The model exhibits a scenario where stationary patterns become spatiotemporally chaotic through two successive bifurcations. First, the pattern undergoes a subcritical Andronov-Hopf bifurcation leading to an oscillatory pattern. Subsequently, a secondary bifurcation gives rise to an oscillation with an incommensurable frequency with respect to the former one. This last bifurcation is responsible for the spatiotemporally chaotic behavior. The Lyapunov spectrum enables us to identify the complex behavior observed as spatiotemporal chaos, and also from the larger Lyapunov exponents characterize the above instabilities.

Cross-diffusion in the two-variable Oregonator model

We explore the effect of cross-diffusion on pattern formation in the two-variable Oregonator model of the Belousov-Zhabotinsky reaction. For high negative cross-diffusion of the activator (the activator being attracted towards regions of increased inhibitor concentration) we find, depending on the values of the parameters, Turing patterns, standing waves, oscillatory Turing patterns, and quasi-standing waves. For the inhibitor, we find that positive cross-diffusion (the inhibitor being repelled by increasing concentrations of the activator) can induce Turing patterns, jumping waves and spatially modulated bulk oscillations. We qualitatively explain the formation of these patterns. With one model we can explain Turing patterns, standing waves and jumping waves, which previously was done with three different models.

Chiral symmetry breaking in a reaction-diffusion system

The emergence of order in far-from-equilibrium systems is often accompanied by the formation of spatially asymmetric patterns. About 30 years ago, a general mechanism to select a chiral solution by coupling a reaction-diffusion system to an external chiral electric field was proposed by Nicolis and Prigogine [Proc. Natl. Acad. Sci. USA 78, 659 (1981)]. However, no experimental or even numerical evidence in reaction-diffusion systems has been reported yet. Here we report a chiral symmetry-breaking phenomenon in a reaction-diffusion system coupled to a circularly polarized electric field (CPEF). Specifically, we show that the CPEF breaks the zero-rotation chiral symmetry between clockwise and counterclockwise spiral defects and that ordered spiral waves with preferred chirality arise from defect-mediated turbulence. The occurrence of such chiral symmetry breaking can be understood by the competition between spiral defects with opposite chirality.

Noise-induced spatiotemporal patterns in Hodgkin-Huxley neuronal network

The effect of noise on the pattern selection in a regular network of Hodgkin-Huxley neurons is investigated, and the transition of pattern in the network is measured from subexcitable to excitable media. Extensive numerical results confirm that kinds of travelling wave such as spiral wave, circle wave and target wave could be developed and kept alive in the subexcitable network due to the noise. In the case of excitable media under noise, the developed spiral wave and target wave could coexist and new target-like wave is induced near to the border of media. The averaged membrane potentials over all neurons in the network are calculated to detect the periodicity of the time series and the generated traveling wave. Furthermore, the firing probabilities of neurons in networks are also calculated to analyze the collective behavior of networks.

Type of spiral wave with trapped ions

Pattern formation in ultracold quantum systems has recently received a great deal of attention. In this work, we investigate a two-dimensional model system simulating the dynamics of trapped ions. We find a spiral wave that is rigidly rotating, but with a peculiar core region in which adjacent ions oscillate in antiphase. The formation of this spiral wave is ascribed to the excitability previously reported by Lee and Cross. The breakup of the spiral wave is probed and, especially, an extraordinary scenario of the disappearance of the spiral wave, caused by spontaneous expansion of the antiphase core, is unveiled.

Two staged pattern formation and spontaneous wave breakup in the ferroin-bromate-pyrocatechol reaction

This report investigated nonlinear spatiotemporal behavior in the ferroin-bromate-pyrocatechol reaction, in which two stages of wave formation, separated by several hours of quiescent period, were observed. In addition to its great photosensitivity, the second stage wave activity could undergo spontaneous breakups at broad reaction conditions. Analysis based on one-dimensional space-time plot suggests that the breakup was caused by propagation slowdown of the leading wave. Due to the presence of coupled autocatalytic reactions, the propagation of the initial and the second stage waves exhibited different and subtle responses to the variation of the concentration of each reagent.

Length scale of interaction in spatiotemporal chaos

Extensive systems have no long scale correlations and behave as a sum of their parts. Various techniques are introduced to determine a characteristic length scale of interaction beyond which spatiotemporal chaos is extensive in reaction-diffusion networks. Information about network size, boundary condition, or abnormalities in network topology gets scrambled in spatiotemporal chaos, and the attenuation of information provides such characteristic length scales. Space-time information flow associated with the recovery of spatiotemporal chaos from finite perturbations, a concept somewhat opposite to the paradigm of Lyapunov exponents, defines another characteristic length scale. High-precision computational studies of asymptotic spatiotemporal chaos in the complex Ginzburg-Landau system and transient spatiotemporal chaos in the Gray-Scott network show that these different length scales are comparable and thus suitable to define a length scale of interaction. Preliminary studies demonstrate the relevance of these length scales for stable chaos.

Novel type of amplitude spiral wave in a two-layer system

Interaction of spiral waves in a two-layer system described by a model of coupled complex Ginzburg-Landau equations with negative-feedback couplings ε(1) and ε(2) is studied. Synchronization of two spiral waves can be broadly found if ε(1)+ε(2) is sufficiently large. Prior to the synchronization, under the condition of strongly asymmetric coupling (∣ε(1)-ε(2)∣≫0), a novel type of spiral wave, amplitude spiral wave, exists in the driven system. The pattern of amplitude spiral wave shows the spiral in the amplitude and without a singularity point (tip), compared to usual spiral waves known for phase with amplitude uniform far away from tips and rotating around tips.

Pattern formation through phase bistability in oscillatory systems with space-modulated forcing

We propose a novel forcing technique of spatially extended self-oscillatory systems able to excite phase bistability and the dissipative structures associated with it. The forcing is time periodic at a frequency close to the oscillators' frequency and is spatially modulated. The effects of this type of forcing are demonstrated analytically and numerically in a directly driven complex Ginzburg-Landau equation. Both spatially periodic and spatially random drives prove to be effective.

Spontaneous suppression of spiral turbulence based on feedback strategy

Control of defect-mediated turbulence or spiral turbulence has been a subject of pattern dynamics study in the last two decades. As a result, many control strategies has been proposed in theory and tested in computer simulation. Here we propose a self-regulation system that can spontaneously eliminate topological defects in the reaction-diffusion system that suffers from the Doppler instability. The core of the self-regulation scheme is a feedback introduction of local inhomogeneities around spiral tips. We used modified FitzHugh-Nagumo model in computer simulations and light-sensitive Belousov-Zhabotinsky reaction in experiments to test the functionality of the proposed system. We found that if the inhomogeneities were induced on the spiral tips once they were identified, the spiral can be exterminated thus the spiral turbulence can be controlled. Compared with the studies on anatomic obstacles in heart, this self-regulation model may provide a possible mechanism for spontaneous termination of cardiac fibrillation.

Theoretical analysis of defect-mediated turbulence in a catalytic surface reaction

We present a statistical analysis of defect-mediated turbulence in a kinetic model of catalytic CO oxidation on Pt(110). A probabilistic description based on the gain and loss rates of defects is derived. For low values of the CO partial pressure the statistics of topological defects agree with earlier results for the complex Ginzburg-Landau equation. For high values of the CO partial pressure, we observe an additional autocatalytic reproduction of defects that results in a linear dependence of the defect creation rate on the number of defects in the system. The role of correlations between defects of opposite topological charge was found to be weaker than in the experimental system.