From quasiperiodicity to chaos in the Belousov–Zhabotinskii reaction. I. Experiment

Glycolytic Wave Patterns in a Simple Reaction-diffusion System with Inhomogeneous Influx: Dynamic Transitions

An inhomogeneous profile of chemostatted species generates a rich variety of patterns in glycolytic waves depicted in a Selkov reaction-diffusion framework here. A key role played by diffusion amplitude and symmetry in the chemostatted species profile in dictating the fate of local spatial dynamics involving periodic, quasiperiodic, and chaotic patterns and transitions among them are investigated systematically. More importantly, various dynamic transitions, including wave propagation direction changes, are illustrated in interesting situations. Besides numerical results, our analytical formulation of the amplitude equation connecting complex Ginzburg-Landau and Lambda-omega representation shed light on the phase dynamics of the system. This systematic study of the glycolytic reaction-diffusion wave is in line with previous experimental results in open spatial reactor and will provide a knowledge about the dynamics that shape and control biological information processing and related phenomena.

Synchronization of thermoacoustic quasiperiodic oscillation by periodic external force

Quasiperiodic oscillations can occur in nonequilibrium systems where two or more frequency components are generated simultaneously. Many studies have explored the synchronization of periodic and chaotic oscillations; however, the synchronization of quasiperiodic oscillations has not received much attention. This study experimentally documents forced synchronization of the quasiperiodic state and the internally locked state of a thermoacoustic oscillator system. This system consists of a gas-filled resonance tube with a nonuniform cross-sectional area. The thermoacoustic oscillator was designed and built in such a way that nonlinear interactions between the fundamental acoustic oscillation mode and the third mode of the gas column are controlled by a temperature difference that is locally created in the resonance tube. Bifurcation diagrams were mapped out by changing the forcing strength and frequency. Separated Arnold tongues were found and both modes were entrained to the external force through complete synchronization. A saddle-node bifurcation was observed in the route from partial to complete synchronization when the forcing strength was relatively weak. However, a Hopf (torus-death) bifurcation was observed when the forcing was relatively strong. In the internally locked state, the bifurcation occurred after the internal locking was broken down by the external force.

Dynamic Environments as a Tool to Preserve Desired Output in a Chemical Reaction Network

Current efforts to design functional molecular systems have overlooked the importance of coupling out-of-equilibrium behaviour with changes in the environment. Here, the authors use an oscillating reaction network and demonstrate that the application of environmental forcing, in the form of periodic changes in temperature and in the inflow of the concentration of one of the network components, removes the dependency of the periodicity of this network on temperature or flow rates and enforces a stable periodicity across a wide range of conditions. Coupling a system to a dynamic environment can thus be used as a simple tool to regulate the output of a network. In addition, the authors show that coupling can also induce an increase in behavioural complexity to include quasi-periodic oscillations.

Nonlinear behavior and fluctuation-induced dynamics in the photosensitive Belousov-Zhabotinsky reaction

The photosensitive Belousov-Zhabotinsky (pBZ) reaction has been used extensively to study the properties of chemical oscillators. In particular, recent experiments revealed the existence of complex spatiotemporal dynamics for systems consisting of coupled micelles (V < 10^-21 L) or droplets (V ≈ [10^-8-10^-11] L) in which the pBZ reaction takes place. These results have been mostly understood in terms of reaction-diffusion models. However, in view of the small size of the droplets and micelles, large fluctuations of concentrations are to be expected. In this work, we investigate the role of fluctuations on the dynamics of a single droplet with stochastic simulations of an extension of the Field-Körös-Noyes (FKN) model taking into account the photosensitivity. The birhythmicity and chaotic behaviors predicted by the FKN model in the absence of fluctuations become transient or intermittent regimes whose lifetime decreases with the size of the droplet. Simple oscillations are more robust and can be observed even in small systems (V > 10^-12 L), which justifies the use of deterministic models in microfluidic systems of coupled oscillators. The simulations also reveal that fluctuations strongly affect the efficiency of inhibition by light, which is often used to control the kinetics of these systems: oscillations are found for parameter values for which they are supposed to be quenched according to deterministic predictions.

Dynamic transitions in a model of the hypothalamic-pituitary-adrenal axis

Dynamic properties of a nonlinear five-dimensional stoichiometric model of the hypothalamic-pituitary-adrenal (HPA) axis were systematically investigated. Conditions under which qualitative transitions between dynamic states occur are determined by independently varying the rate constants of all reactions that constitute the model. Bifurcation types were further characterized using continuation algorithms and scale factor methods. Regions of bistability and transitions through supercritical Andronov-Hopf and saddle loop bifurcations were identified. Dynamic state analysis predicts that the HPA axis operates under basal (healthy) physiological conditions close to an Andronov-Hopf bifurcation. Dynamic properties of the stress-control axis have not been characterized experimentally, but modelling suggests that the proximity to a supercritical Andronov-Hopf bifurcation can give the HPA axis both, flexibility to respond to external stimuli and adjust to new conditions and stability, i.e., the capacity to return to the original dynamic state afterwards, which is essential for maintaining homeostasis. The analysis presented here reflects the properties of a low-dimensional model that succinctly describes neurochemical transformations underlying the HPA axis. However, the model accounts correctly for a number of experimentally observed properties of the stress-response axis. We therefore regard that the presented analysis is meaningful, showing how in silico investigations can be used to guide the experimentalists in understanding how the HPA axis activity changes under chronic disease and/or specific pharmacological manipulations.

Complex and chaotic oscillations in a model for the catalytic hydrogen peroxide decomposition under open reactor conditions

Numerous periodic and aperiodic dynamic states obtained in a model for hydrogen peroxide decomposition in the presence of iodate and hydrogen ions (the Bray-Liebhafsky reaction) realized in an open reactor (CSTR), where the flow rate was the control parameter, have been investigated numerically. Between two Hopf bifurcation points, different simple and complex oscillations and different routes to chaos were observed. In the region of the mixed-mode evolution of the system, the transitions between two successive mixed-mode simple states are realized by period-doubling of the initial state leading to a chaotic window in which the next dynamic state emerges mixed with the initial one. It appears in increasing proportions in concatenated patterns until total domination. Thus, with increasing flow rate the period-doubling route to chaos was obtained, whereas with decreasing flow rate the peak-adding route to chaos was obtained. Moreover, in very narrow regions of flow rates, chaotic mixtures of mixed-mode patterns were observed. This evolution of patterns repeats until the end of the mixed-mode region at high flow rates that corresponds to chaotic mixtures of one large and many small amplitude oscillations. Starting from the reverse Hopf bifurcation point and decreasing the flow rate, simple small amplitude sinusoidal oscillations were encountered and then the period-doubling route to chaos. With a further decreasing flow rate, the mixed-mode oscillations emerge inside the chaotic window.

Line-defects-mediated complex-oscillatory spiral waves in a chemical system

In this paper, we summarize our experimental observations on complex-oscillatory spiral waves that arise in a Belousov-Zhabotinsky (BZ) reaction-diffusion system. The observed wave structures generically bear line defects across which the phase of local oscillation changes by a multiple of 2 pi. The local oscillation at every spatial point along a line defect of period-2 (P-2) oscillatory media is period-1 (P-1) oscillatory. For the homogeneous BZ reaction can be excitable, simply periodic, complex periodic, or chaotic as the control parameters are tuned, a number of different complex wave states are revealed. A two-dimensional phase diagram, which includes domains of P-2 oscillatory spirals, intermittently breathing spirals, period-3 (P-3) oscillatory spirals, two different types of mixed-mode periodic spirals, and line-defect-mediated turbulence, is constructed. Several different transitions among different dynamic states are described systematically. In all cases, line defects are found to play an important role.

A new route to chaos: sequences of topological torus bifurcations

We consider a sequence of topological torus bifurcations (TTBs) in a nonlinear, quasiperiodic Mathieu equation. The sequence of TTBs and an ensuing transition to chaos are observed by computing the principal Lyapunov exponent over a range of the bifurcation parameter. We also consider the effect of the sequence on the power spectrum before and after the transition to chaos. We then describe the topology of the set of knotted tori that are present before the transition to chaos. Following the transition, solutions evolve on strange attractors that have the topology of fractal braids in Poincare sections. We examine the topology of fractal braids and the dynamics of solutions that evolve on them. We end with a brief discussion of the number of TTBs in the cascade that leads to chaos.

Dynamical behavior of a complex fluid near an out-of-equilibrium transition: approaching simple rheological chaos

We report here an extensive study of sustained oscillations of the viscosity of a complex fluid near an out-of-equilibrium transition. Using well defined protocols, we perform rheological measurements of the onion texture near a layering transition in a Couette flow. This complex fluid exhibits sustained oscillations of the viscosity, on a large time scale (500 s) at controlled stress. These oscillations are directly correlated to an oscillating microstructural change of the texture of the fluid. We observe a great diversity of dynamical behavior and show that there is a coupling with spatial effects in the inverted Delta v direction. This is in agreement with a careful analysis of the temporal series of the viscosity with the dynamical system theory. This analysis indicates that the observed dynamical responses do not strictly correspond to three-dimensional chaotic states, probably because some spatiotemporal effects are present and are likely to play an important role.

Oscillating viscosity in a lyotropic lamellar phase under shear flow

We report some time-dependent behavior of lyotropic lamellar phase under shear flow. At fixed stress, near a layering instability, the system presents an oscillating shear rate. We build up a new stress versus shear rate diagram that includes temporal behavior. This diagram is made of two distinct branches of stationary states which correspond, respectively, to disordered and ordered multilamellar vesicle phases. When increasing the shear stress, prior to the transition to the ordered structural state, sustained oscillations of the viscosity are recorded. They correspond to periodic structural change of the entire sample between a disordered and a ordered state of multilamellar vesicles.

New chaotic behavior and its effective control in Belousov–Zhabotinsky reaction

A new "spiral-in-to-spiral-out" mixed-mode type chaos of the Belousov–Zhabotinsky (BZ) reaction in a continuous flow stirred tank reactor (CSTR) is observed. We discussed its evolution and characterization in detail. We controlled the new chaos into steady periodic behavior by two control methods: (i) a self-adaptive delayed feedback controlling scheme; and (ii) a linear self-interacting feedback controlling scheme.Key words: BZ-CSTR reaction, chaos, control.

Integration of non-linear cellular mechanisms regulating microvascular perfusion

It is becoming increasingly evident that interactions between the different cell types present in the vessel wall and the physical forces that result from blood flow are highly complex. This short article will review evidence that irregular fluctuations in vascular resistance are generated by non-linearity in the control mechanisms intrinsic to the smooth muscle cell and can be classified as chaotic. Non-linear systems theory has provided insights into the mechanisms involved at the cellular level by allowing the identification of dominant control variables and the construction of one-dimensional iterative maps to model vascular dynamics. Experiments with novel peptide inhibitors of gap junctions have shown that the coordination of aggregate responses depends on direct intercellular communication. The sensitivity of chaotic trajectories to perturbation may nevertheless generate a high degree of variability in the response to pharmacological interventions and altered perfusion conditions.

Simplification of the quasiperiodic route to chaos in agonist-induced vasomotion by iterative circle maps

We have shown that the patterns of vasomotion induced by histamine in isolated rabbit ear resistance arteries can be described in terms of iterative circle maps that model the dynamics of coupled nonlinear oscillators. Cyclopiazonic acid (CPA), an inhibitor of the sarcoplasmic reticulum Ca(2+)-adenosinetriphosphatase pump, consistently transformed chaotic behavior into characteristic periodic oscillations known as mixed-mode responses, which consist of mixtures of large- and small-amplitude excursions and represent frequency-locked states. Quasiperiodicity, which reflects the interaction of oscillators with incommensurate frequencies, was also observed, although in a smaller number of experiments. The patterns of mixed-mode complexes found at different CPA concentrations allowed the derivation of firing numbers, i.e., number of large oscillations/sum of number of small and large oscillations, and the sequences in which they emerged conformed to Farey arithmetic. Two-dimensional return maps derived by Poincaré section of phase space representations of the dynamics were used to compute the mean number of rotations per iteration on the circle, i.e., the winding number. Plots of winding number against firing number revealed a devil's staircase-type structure. Experiments with verapamil, a voltage-operated L-type Ca(2+)-channel antagonist, confirmed that influx of extracellular Ca2+ was essential to sustain chaos, quasiperiodicity, and mixed-mode responses. Nonlinear coupling between cytosolic and membrane events in rabbit ear arteries thus results in a self-organized dynamics that collapses to that predicted by the theory of simple circle maps.

Infinite period and Hopf bifurcations for the pH-regulated oscillations in a semibatch reactor (H(2)O(2)-Cu(2+)-S(2)O(2-) (3)-NaOH system)

Dynamic behavior of the pH-regulated oscillations has been studied for the hydrogen peroxide oxidation of thiosulfate ions in the presence of trace amounts of copper(II) ions in a semibatch reactor. A solution of 0.08 M Na(2)S(2)O(3) and 0.112 M NaOH was flowed at 0.160 mL/min into 300 mL of solution containing the H(2)O(2) and Cu(2+) in a vessel. There exists a critical value of the H(2)O(2) or Cu(2+) concentrations below which the system does not oscillate. The oscillations appear due to an infinite period bifurcation at low initial concentrations of the H(2)O(2). The initial concentration of Cu(2+) may be considered as a bifurcation parameter in this case. Increase of the initial hydrogen peroxide concentration causes the pH-regulated oscillations through a nondegenerate supercritical Hopf bifurcation. The classification of bifurcations is based on the analysis of the behavior of oscillation amplitude and period at different initial concentrations of the H(2)O(2) and Cu(2+). Our results show a possibility to distinguish different scenarios for the appearance of transient oscillations in semibatch experiments. (c) 1996 American Institute of Physics.

Deterministic chaos in the Belousov-Zhabotinsky reaction: Experiments and simulations

An account of the experimental discovery of complex dynamical behavior in the continuous-flow, stirred tank reactor (CSTR) Belousov-Zhabotinsky (BZ) reaction, as well as numerical simulations based on the BZ chemistry are given. The most recent four- and three-variable models that are deduced from the well-accepted, updated chemical mechanism of the BZ reaction and which exhibit robust chaotic states are summarized. Chaos has been observed in experiments and simulations embedded in the regions of complexities at both low and high flow rates. The deterministic nature of the observed aperiodicities at low flow rates is unequivocally established. However, controversy still remains in the interpretation of certain aperiodicities observed at high flow rates.