Anomalous Dispersion of Chemical Waves in a Homogeneously Catalyzed Reaction System

Propagating wave merging in a precipitation reaction

Propagating precipitation waves are a remarkable form of spatiotemporal behavior that arise through the coupling of reaction, diffusion, and precipitation. We study a system with a sodium hydroxide outer electrolyte and an aluminum hydroxide inner electrolyte. In a redissolution Liesegang system, a single propagating precipitation band moves down through the gel, with precipitate formed at the band front and precipitate dissolved at the band back. Complex spatiotemporal waves occur within the propagating precipitation band, including counter-rotating spiral waves, target patterns, and annihilation of waves on collision. We have also carried out experiments in thin slices of gel, which have revealed propagating waves of a diagonal precipitation feature within the primary precipitation band. These waves display a wave merging phenomenon in which two horizontally propagating waves merge into a single wave. Computational modeling permits the development of a detailed understanding of the complex dynamical behavior.

Temporally quasiperiodic data, propagating in the laboratory frame, can be rendered periodic by Galilean transformation

For a broad class of distributions of temperature, concentration, or another quantity propagating rectilinearly, we show that temporally quasiperiodic behavior in the laboratory frame can be rendered periodic by Galilean transformation. The approach is illustrated analytically and numerically using as an example a closed-form model distribution generated from a one-dimensional partial differential equation, and a detailed process is developed to determine frame speed from more general quasiperiodic, one-dimensional, temporally- and spatially-discretized data. The approach is extended to two- and three-dimensional rectilinear propagation, and its application to nonrectilinear propagation, along with implications for interpreting noise-corrupted data, are also discussed.

Chemo-hydrodynamic pulsations in simple batch A + B → C systems

Spatio-temporal oscillations can be induced under batch conditions with ubiquitous bimolecular reactions in the absence of any nonlinear chemical feedback, thanks to an active interplay between the chemical process and chemically driven hydrodynamic flows. When two reactants A and B, initially separated in space, react upon diffusive contact, they can power convective flows by inducing a localized variation of surface tension and density at the mixing interface. These flows feedback with the reaction-diffusion dynamics, bearing damped or sustained spatio-temporal oscillations of the concentrations and flow field. By means of numerical simulations, we detail the mechanism underlying these chemohydrodynamic oscillations and classify the main dynamical scenarios in the relevant space drawn by parameters ΔM and ΔR, which rule the surface tension- and buoyancy-driven contributions to convection, respectively. The reactor height is found to play a critical role in the control of the dynamics. The analysis reveals the intimate nature of these oscillatory phenomena and the hierarchy among the different phenomena at play: oscillations are essentially hydrodynamic and the chemical process features the localized trigger for Marangoni flows unstable toward oscillatory instabilities. The characteristic size of Marangoni convective rolls mainly determines the critical conditions and properties of the oscillations, which can be further tuned or suppressed by the buoyancy competition. We finally discuss the possible experimental implementation of such a class of chemo-hydrodynamic oscillator and its implications in fundamental and applied terms.

Long-Lasting Complex Reaction Behavior in a Closed Ferroin-Bromate-Hydroxybenzenesulfonate System

The bromate-phenolsulfonate reaction was found to exhibit spontaneous oscillations in a batch reactor, where the addition of small amounts of ferroin would result in nonoscillatory behavior. As the ferroin concentration was increased, the system produced very rich nonlinear behavior, including three isolated oscillatory regimes that were separated by as long as 48 h nonoscillatory period. The long-lasting nonlinear behavior may be attributed to the slow desulfonation of phenolsulfonate in an acidic solution, forming phenol-like intermediates. However, unlike the bromate-phenol oscillator, oxygen was found to greatly influence the reaction, and various complex oscillations could be observed by tuning the oxygen concentration. Mechanistic studies performed through employing ¹H NMR spectroscopy and mass spectrometry to measure intermediate species at different stages of the reaction were able to identify 1,4-benzoquinone, 2-bromo-1,4-benzoquinone, 2,6-dibromo-1,4-benzoquinone, and 2,4,6-tribromophenol as major components during the reaction.

Unusual spiral wave dynamics in the Kessler-Levine model of an excitable medium

The Kessler-Levine model is a two-component reaction-diffusion system that describes spatiotemporal dynamics of the messenger molecules in a cell-to-cell signaling process during the aggregation of social amoeba cells. An excitation wave arising in the model has a phase wave at the wave back, which simply follows the wave front after a fixed time interval with the same propagation velocity. Generally speaking, the medium excitability and the refractoriness are two important factors which determine the spiral wave dynamics in any excitable media. The model allows us to separate these two factors relatively easily since the medium refractoriness can be changed independently of the medium excitability. For rigidly rotating waves, the universal relationship has been established by using a modified free-boundary approach, which assumes that the front and the back of a propagating wave are thin in comparison to the wave plateau. By taking a finite thickness of the domain boundary into consideration, the validity of the proposed excitability measure has been essentially improved. A novel method of numerical simulation to suppress the spiral wave instabilities is introduced. The trajectories of the spiral tip observed for a long refractory period have been investigated under a systematic variation of the medium refractoriness.

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.

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.

Gas production in the bromate-pyrocatechol oscillator

A significant amount of gas production has been observed in the bromate-pyrocatechol oscillator under high concentrations of bromate and pyrocatechol. The observation is in contrast to the general perception that aromatic compounds can form bromate-based oscillators that are free of gas bubbles, which is a desired property in investigating pattern formation. Analysis with (1)H NMR, (13)C NMR, mass spectrometry, and X-ray crystallography illustrate the production of 5-(dibromomethylene)-2(5H)-furanone from pyrocatechol, where the loss of one carbon atom from the aromatic ring causes the formation of gas bubbles. Possible mechanisms have been proposed to explain the observed phenomenon.

Propagation failures, breathing fronts, and nonannihilation collisions in the ferroin-bromate-pyrocatechol system

The emergence of propagating pulses was investigated with the photosensitive ferroin-bromate-pyrocatechol reaction in capillary tubes, in which various interesting spatiotemporal behaviors such as propagation failure, breathing fronts, and transitions between propagating pulses and fronts have been observed. Rather than a mutual annihilation, the collision of a propagating pulse and a growing front forces the front to recede gradually. A phase diagram in the pyrocatechol-bromate concentration space shows that the pulse instabilities take place throughout the conditions at which the system generates wave activities, suggesting that the presence of coupled autocatalytic feedbacks may facilitate the onset of pulse instabilities.

Bifurcation analysis of solitary and synchronized pulses and formation of reentrant waves in laterally coupled excitable fibers

We study the dynamics of a reaction-diffusion system comprising two mutually coupled excitable fibers. We consider a case in which the dynamical properties of the two fibers are nonidentical due to the parameter mismatch between them. By using the spatially one-dimensional FitzHugh-Nagumo equations as a model of a single excitable fiber, synchronized pulses are found to be stable in some parameter regime. Furthermore, there exists a critical coupling strength beyond which the synchronized pulses are stable for any amount of parameter mismatch. We show the bifurcation structures of the synchronized and solitary pulses and identify a codimension-2 cusp singularity as the source of the destabilization of synchronized pulses. When stable solitary pulses in both fibers disappear via a saddle-node bifurcation on increasing the coupling strength, a reentrant wave is formed. The parameter region, where a stable reentrant wave is observed in direct numerical simulation, is consistent with that obtained by bifurcation analysis.

Three-dimensional spiral waves in an excitable reaction system: initiation and dynamics of scroll rings and scroll ring pairs

We report experimental results on spiral and scroll waves in the 1,4-cyclohexanedione Belousov-Zhabotinsky reaction. The propagating concentration waves are detected by two-dimensional photometry and optical tomography. Wave pulses can disappear in front-to-front and front-to-back collisions. This anomaly causes the nucleation of vortices from collisions of three nonrotating waves. In three-dimensional systems, these vortices are scroll rings that rotate around initially circular filaments. Depending on reactant concentrations, the filaments shrink or expand indicating positive and negative filament tensions, respectively. Shrinkage results in vortex annihilation. Expansion is accompanied by filament buckling and bending, which is interpreted as developing Winfree turbulence. We also describe the initiation of scroll ring pairs in four-wave collisions. The two filaments are stacked on top of each other and their motion suggests filament repulsion.

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.

Chemical oscillations and waves in the catalyzed bromate-pyrocatechol reaction

Long time series of temporal oscillations and wave formation are observed in the catalyzed bromate-pyrocatechol reaction conducted in a batch reactor, in which the induction time is insensitive to the presence of ferroin but is greatly shortened by Ce(III) or Mn(II). On the other hand, the number of oscillations is significantly increased by ferroin, while it is less sensitive to Ce(III) and Mn(II). The ferroin-catalyzed system also exhibits strong photosensitivity, in which illumination could quench the oscillatory behavior. A phase diagram illustrates that the oscillatory behavior of the studied system is more sensitive to the ratio of [pyrocatechol]/[bromate] than their absolute concentrations. Reactions conducted in a spatially extended medium show that the ferroin-catalyzed system supports a two-stage pattern formation with the wave activity surviving for up to 10 h.

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.

Propagation failures, breathing pulses, and backfiring in an excitable reaction-diffusion system

We report results from experiments with a pseudo-one-dimensional Belousov-Zhabotinsky reaction that employs 1,4-cyclohexanedione as its organic substrate. This excitable system shows traveling oxidation pulses and pulse trains that can undergo complex sequences of propagation failures. Moreover, we present examples for (i) breathing pulses that undergo periodic changes in speed and size and (ii) backfiring pulses that near their back repeatedly generate new pulses propagating in opposite direction.

Creating bound states in excitable media by means of nonlocal coupling

We consider pulses of excitation in reaction-diffusion systems subjected to nonlocal coupling. This coupling represents long-range connections between the elements of the medium; the connection strength decays exponentially with the distance. Without coupling, pulses interact only repulsively and bound states with two or more pulses propagating at the same velocity are impossible. Upon switching on nonlocal coupling, pulses begin to interact attractively and form bound states. First we present numerical results on the emergence of bound states in the excitable Oregonator model for the photosensitive Belousov-Zhabotinsky reaction with nonlocal coupling. Then we show that the appearance of bound states is provided solely by the exponential decay of nonlocal coupling and thus can be found in a wide class of excitable systems, regardless of the particular kinetics. The theoretical explanation of the emergence of bound states is based on the bifurcation analysis of the profile equations that describe the spatial shape of pulses. The central object is a codimension-4 homoclinic orbit which exists for zero coupling strength. The emergence of bound states is described by the bifurcation to 2-homoclinic solutions from the codimension-4 homoclinic orbit upon switching on nonlocal coupling. We stress that the high codimension of the bifurcation to bound states is generic, provided that the coupling range is sufficiently large.

Propagation failure dynamics of wave trains in excitable systems

We report experimental and numerical results on temporal patterns of propagation failures in reaction-diffusion systems. Experiments employ the 1,4-cyclohexanedione Belousov-Zhabotinsky reaction. The propagation failures occur in the frontier region of the wave train and can profoundly affect its expansion speed. The specific rhythms observed vary from simple periodic to highly complex and possibly chaotic sequences. All but the period-1 sequences are found in the transition region between "merging" and "tracking" dynamics, which correspond to wave behavior caused by two qualitatively different types of anomalous dispersion relations.

Coexistence of Two Bifurcation Regimes in a Closed Ferroin-Catalyzed Belousov−Zhabotinsky Reaction

The ferroin-catalyzed Belousov-Zhabotinsky (BZ) reaction was studied in a batch reactor under anaerobic conditions and was found to evolve through two separated regimes of complex oscillations. Significantly, the two bifurcation regimes exhibited qualitatively different dependence on compositions of the reaction mixture, i.e., initial concentrations of bromate, sulfuric acid, malonic acid, and ferroin. The reaction temperature also showed opposite effects on the two bifurcation regimes, in which complexities of the first bifurcation regime were enhanced while oscillations in the second bifurcation regime became simpler as a result of decreasing temperature. Numerical simulations with a 12-variable model developed specifically for the ferroin-BZ system were able to reproduce transient complex oscillations observed in experiments. These calculations further illustrated that reactions such as ferroin and HOBr, ferroin and HBrO2, and ferriin and Br- were not essential in describing complex dynamics of the ferroin-BZ reaction.

Dynamics of excitation pulses with attractive interaction: kinematic analysis and chemical wave experiments

We present a theoretical analysis of stacking and destacking wave trains in excitable reaction-diffusion systems with anomalous velocity-wavelength dependence. For linearized dispersion relations, kinematic analysis yields an analytical function that rigorously describes front trajectories. The corresponding accelerations have exactly one extremum that slowly decays with increasing pulse number. For subsequent pulses these maxima occur with a lag time equal to the inverse slope of the linearized dispersion curve. These findings are reproduced in experiments with chemical waves in the 1,4-cyclohexanedione Belousov-Zhabotinsky reaction but should be also applicable to step bunching on crystal surfaces and certain traffic phenomena.

Tracking waves and vortex nucleation in excitable systems with anomalous dispersion

We report experimental results obtained from a chemical reaction-diffusion system in which wave propagation is limited to a finite band of wavelengths and in which no solitary pulses exist. Wave patterns increase their size through repeated annihilation events of the frontier pulse that allow the succeeding pulses to advance farther. A related type of wave dynamics involves a stable but slow frontier pulse that annihilates subsequent waves in front-to-back collisions. These so-called merging dynamics give rise to an unexpected form of spiral wave nucleation. All of these phenomena are reproduced by a simple, three-species reaction-diffusion model that reveals the importance of the underlying anomalous dispersion relation.

Multiple mechanisms of spiral wave breakup in a model of cardiac electrical activity

It has become widely accepted that the most dangerous cardiac arrhythmias are due to reentrant waves, i.e., electrical wave(s) that recirculate repeatedly throughout the tissue at a higher frequency than the waves produced by the heart's natural pacemaker (sinoatrial node). However, the complicated structure of cardiac tissue, as well as the complex ionic currents in the cell, have made it extremely difficult to pinpoint the detailed dynamics of these life-threatening reentrant arrhythmias. A simplified ionic model of the cardiac action potential (AP), which can be fitted to a wide variety of experimentally and numerically obtained mesoscopic characteristics of cardiac tissue such as AP shape and restitution of AP duration and conduction velocity, is used to explain many different mechanisms of spiral wave breakup which in principle can occur in cardiac tissue. Some, but not all, of these mechanisms have been observed before using other models; therefore, the purpose of this paper is to demonstrate them using just one framework model and to explain the different parameter regimes or physiological properties necessary for each mechanism (such as high or low excitability, corresponding to normal or ischemic tissue, spiral tip trajectory types, and tissue structures such as rotational anisotropy and periodic boundary conditions). Each mechanism is compared with data from other ionic models or experiments to illustrate that they are not model-specific phenomena. Movies showing all the breakup mechanisms are available at http://arrhythmia.hofstra.edu/breakup and at ftp://ftp.aip.org/epaps/chaos/E-CHAOEH-12-039203/ INDEX.html. The fact that many different breakup mechanisms exist has important implications for antiarrhythmic drug design and for comparisons of fibrillation experiments using different species, electromechanical uncoupling drugs, and initiation protocols. (c) 2002 American Institute of Physics.

Excitation fronts in a spatially modulated light-sensitive Belousov-Zhabotinsky system

The evolution of excitation wave fronts in a spatially modulated light-sensitive Belousov-Zhabotinsky system is investigated experimentally and theoretically. The excitation wave propagates in a thin, quasi-two-dimensional reaction layer, which is illuminated through a periodical gray level mask. The light-induced differences in excitability and velocity give rise to a temporal and spatial modulation of the initially flat fronts. The experimental front evolution is described in the framework of a kinematical theory as developed earlier for nonuniformly curved systems.

Shock structures and bunching fronts in excitable reaction-diffusion systems

We report experimental results on the dynamics of excitation waves in a modified Belousov-Zhabotinsky reaction. The waves in this system obey nonmonotonic dispersion relations. This anomaly induces the stacking of excitation fronts into patterns with stable interpulse distances. The stacking process creates either a traveling shock structure or a cascade of bunching events in which metastable wave packets are formed. The direction and the speed of the shock are explained in terms of a simple geometrical analysis. We also present experimental evidence for the corresponding instabilities in two-dimensional systems. Here, wave stacking generates atypical structures in the collision of target patterns and wave bunching is accompanied by complex front deformations.