Entrainment in a Chemical Oscillator Chain with a Pacemaker

Traveling waves propagating through coupled microbeads in the Belousov-Zhabotinsky reaction

Spatio-temporal patterns, namely global oscillations (GO) and traveling waves (TW), were investigated in spherical microbeads loaded with a catalyst for the Belousov-Zhabotinsky (BZ) reaction onto the surface (2D-loaded) or the entire volume of the bead (3D-loaded). GO and TW selectively appeared in the 2D- and 3D-loaded beads, respectively, placed on a polyethylene terephthalate (PET) sheet in the catalyst-free BZ solution. We examined two types of coupling of the two beads: 2D-3D and 3D-3D couplings. In both cases, synchronization occurred when the minimum distance between the two beads, l, was shorter than the threshold. Herein, we reported not only temporal information, that is, phase difference, but also spatial information, that is, the directions of the TW propagating through the coupled BZ beads. In the synchronization for the 2D-3D coupling, TW in the 3D-loaded bead were initiated from the point near the 2D-loaded bead as a pacemaker and propagated in the opposite direction. By contrast, the directions of the TW in the 3D-loaded bead changed depending on l in the synchronization for the 3D-3D coupling. These experimental results can be quantitatively reproduced by numerical calculations based on the diffusion dynamics of an activator of the BZ reaction. Our results suggest that the features of spatio-temporal wave propagation are indicative of the configuration of the oscillators.

Switching between Two Oscillatory States Depending on the Electrical Potential

Various spatiotemporal patterns were created on the surface or in the body of cation-exchange resin beads which were loaded with the catalyst of the Belousov-Zhabotinsky (BZ) reaction. Either global oscillations (GO) or traveling waves (TW) and the switching between them were observed in the previous papers, but it was not clear how chemicals contribute to the reaction inside/around the BZ bead. In this paper, we scanned the electrical potential, E, from +1 to -1 V (negative scan) and then turned from -1 to +1 V (positive scan) to control the switching between GO and TW. We found that the electrical switching potential from TW to GO, E_TG, and from GO to TW, E_GT, depended on the scanning direction of E and the diameter of the bead, d. The present study suggests that the electrode-induced increase of the inhibitor, Br^-, and the activator, HBrO₂, around the BZ bead plays an important role in determining E_TG and E_GT.

Synchronization in reaction-diffusion systems with multiple pacemakers

Spatially extended oscillatory systems can be entrained by pacemakers, regions that oscillate with a higher frequency than the rest of the medium. Entrainment happens through waves originating at a pacemaker. Typically, biological and chemical media can contain multiple pacemaker regions, which compete with each other. In this paper, we perform a detailed numerical analysis of how wave propagation and synchronization of the medium depend on the properties of these pacemakers. We discuss the influence of the size and intrinsic frequency of pacemakers on the synchronization properties. We also study a system in which the pacemakers are embedded in a medium without any local dynamics. In this case, synchronization occurs if the coupling determined by the distance and diffusion is strong enough. The transition to synchronization is reminiscent of systems of discrete coupled oscillators.

Effects of stochastic time-delayed feedback on a dynamical system modeling a chemical oscillator

We examine how stochastic time-delayed negative feedback affects the dynamical behavior of a model oscillatory reaction. We apply constant and stochastic time-delayed negative feedbacks to a point Field-Körös-Noyes photosensitive oscillator and compare their effects. Negative feedback is applied in the form of simulated inhibitory electromagnetic radiation with an intensity proportional to the concentration of oxidized light-sensitive catalyst in the oscillator. We first characterize the system under nondelayed inhibitory feedback; then we explore and compare the effects of constant (deterministic) versus stochastic time-delayed feedback. We find that the oscillatory amplitude, frequency, and waveform are essentially preserved when low-dispersion stochastic delayed feedback is used, whereas small but measurable changes appear when a large dispersion is applied.

Chemical communication and dynamics of droplet emulsions in networks of Belousov-Zhabotinsky micro-oscillators produced by microfluidics

Chemical communication leading to synchronization and collective behaviour of dynamic elements, such as cell colonies, is a widespread phenomenon with biological, physical and chemical importance. Such synchronization between elements proceeds via chemical communication by emmision, interdiffusion and reception of specific messenger molecules. On a lab scale, these phenomena can be modeled by encapsulating an oscillating chemical reaction, which serves as a signal (information) sender/receiver element, inside microcompartments such as droplet emulsions, liposomes and polymersomes. Droplets can thus be regarded as single units, able to generate chemical messengers that diffuse in the environment and hence can interact with other compartments. The Belousov-Zhabotinsky (BZ) reaction is a well-known chemical oscillator largely used as a model for complex nonlinear phenomena, including chemical, physical and biological examples. When the BZ-reaction is encapsulated inside microcompartments, its chemical intermediates can serve as messengers by diffusing among different microcompartments, to trigger specific reactions leading to a collective behavior between the elements. The geometry and constitution of the diffusion pathways play an important role in governing the collective behaviour of the system. In this context, microfluidics is not only a versatile tool for mastering the encapsulation process of the BZ-reaction in monodisperse microcompartments, but also for creating geometries and networks with well defined boundaries. The individual compartments can be engineered with selected properties using different surfactants in the case of simple emulsions, or with specific membrane properties in the case of liposomes. Furthermore, it enables the arrangement of these microcompartments in various geometric configurations, where the diffusive coupling pathways between individual compartments are both spatially and chemically well-defined. In this tutorial paper, we review a number of articles reporting various approaches to generate networks of compartmentalized Belousov-Zhabotinsky (BZ) chemical oscillators using microfluidics. In contrast to biological cellular networks, the dynamical characteristics of the BZ-reaction is well-known and, when confined in microcompartments arranged in different configurations with a pure interdiffusive coupling, these communicative microreactors can serve to mimic various types of bio-physical networks, aiding to comprehend the concept of chemical communication.

Synchronization of Coupled Oscillators on a Two-Dimensional Plane

The effect of the transfer rate of signal molecules on coupled chemical oscillators arranged on a two-dimensional plane was systematically investigated in this paper. A microreactor equipped with a surface acoustic wave (SAW) mixer was applied to adjust the transfer rate of the signal molecules in the microreactor. The SAW mixer with adjustable input powers provided a simple means to generate different mixing rates in the microreactor. A robust synchronization of the oscillators was found at an input radio frequency power of 20 dBm, with which the chemical waves were initiated at a fixed site of the oscillator system. With increasing input power, the frequency of the chemical waves was increased, which agreed well with the prediction given by the time-delayed phase oscillator model. Results from the finite element simulation agreed well with the experimental results.

Experimental study of synchronization of coupled electrical self-oscillators and comparison to the Sakaguchi-Kuramoto model

We explore the collective phase dynamics of Wien-bridge oscillators coupled resistively. We carefully analyze the behavior of two coupled oscillators, obtaining a transformation from voltage to effective phase. From the phase dynamics we show that the coupling can be quantitatively described by Sakaguchi's modification to the Kuramoto model. We also examine an ensemble of oscillators whose frequencies are taken from a flat distribution within a fixed frequency interval. We characterize in detail the synchronized cluster, its initial formation, as well as its effect on unsynchronized oscillators, all as a function of a global coupling strength.

Phase-lag synchronization in networks of coupled chemical oscillators

Chemical oscillators with a broad frequency distribution are photochemically coupled in network topologies. Experiments and simulations show that the network synchronization occurs by phase-lag synchronization of clusters of oscillators with zero- or nearly zero-lag synchronization. Symmetry also plays a role in the synchronization, the extent of which is explored as a function of coupling strength, frequency distribution, and the highest frequency oscillator location. The phase-lag synchronization occurs through connected synchronized clusters, with the highest frequency node or nodes setting the frequency of the entire network. The synchronized clusters successively "fire," with a constant phase difference between them. For low heterogeneity and high coupling strength, the synchronized clusters are made up of one or more clusters of nodes with the same permutation symmetries. As heterogeneity is increased or coupling strength decreased, the phase-lag synchronization occurs partially through clusters of nodes sharing the same permutation symmetries. As heterogeneity is further increased or coupling strength decreased, partial synchronization and, finally, independent unsynchronized oscillations are observed. The relationships between these classes of behavior are explored with numerical simulations, which agree well with the experimentally observed behavior.

Combined excitatory and inhibitory coupling in a 1-D array of Belousov–Zhabotinsky droplets

Experimental and theoretical studies of the coupling between Belousov–Zhabotinsky droplets in oil as a function of malonic acid, drop size, drop spacing, and time.

Self-arrangement of cellular circadian rhythms through phase-resetting in plant roots

We discovered a striped pattern of gene expression with circadian rhythms in growing plant roots using bioluminescent imaging of gene expression. Our experimental analysis revealed that the stripe wave in the bioluminescent image originated at the root tip and was caused by a continuous phase resetting of circadian oscillations. Some complex stripe waves containing arrhythmic regions were also observed. We succeeded in describing the formation mechanisms of these patterns using a growing phase oscillator network with a phase-resetting boundary condition.

Chemomechanical synchronization in heterogeneous self-oscillating gels

Using computational modeling, we introduce patches of self-oscillating gels undergoing the Belousov-Zhabotinsky (BZ) reaction into a nonreactive polymer network and thereby demonstrate how these BZ gels can be harnessed to impart remarkable functionality to the entire system. By first focusing on two adjacent patches of BZ gels, we show that the patches' oscillations can become synchronized in phase or out of phase, with the oscillation frequency depending on the synchronization mode and the spatial separation between these domains. We then apply these results to an array of five adjacent BZ patches and by varying the distance between these pieces, we dramatically alter the dynamical behavior of the patterned gel. For example, the sample can be made to exhibit a unidirectional traveling wave or display a concerted expansion and contraction, properties that are valuable for creating gel-based devices, such as micropumps and microactuators. The findings point to a "modular" design approach, which can impart different functionality simply by arranging identical pieces of BZ gels into distinct spatial arrangements within a polymer matrix.

Resonance clustering in globally coupled electrochemical oscillators with external forcing

Experiments are carried out with a globally coupled, externally forced population of limit-cycle electrochemical oscillators with an approximately unimodal distribution of heterogeneities. Global coupling induces mutually entrained (at frequency omega1) states; periodic forcing produces forced-entrained (omegaF) states. As a result of the interaction of mutual and forced entrainment, resonant cluster states occur with equal spacing of frequencies that have discretized frequencies following a resonance rule omegan congruent with nomega1-(n-1)omegaF. Resonance clustering requires an optimal, intermediate global coupling strength; at weak coupling the clusters have smaller sizes and do not strictly follow the resonance rule, while at strong coupling the population behaves similar to a single, giant oscillator.

Synchronization of plant circadian oscillators with a phase delay effect of the vein network

Synchronization phenomena in coupled circadian oscillators of plant leaves were investigated experimentally using bioluminescence technology for a clock gene. Analyzing the phase of circadian oscillation, the phase-wave propagations and the phase delay caused by the vein network were observed. We describe these phase dynamics using a two-layer model with coupled Stuart-Landau equations. Global synchronization of circadian oscillators in the leaf is also investigated.

Loss of coherence in a population of diffusively coupled oscillators

The authors investigate the relationship between the natural frequency distribution of diffusively coupled chemical oscillators and their entrainment by pacemakers. The system consists of micrometer-sized catalyst beads which are coupled to their neighbors by diffusion of the activator/inhibitor species through the catalyst-free Belousov-Zhabotinsky (BZ) reaction solution. The frequency distribution is measured as a function of the beads' number of neighbors. With the maximum number of neighbors, either target waves or disordered patterns are observed in the reaction domain and there is a shift to higher frequencies than those observed in the natural frequency distribution. The loss of coherence between neighbor oscillators is quantified by a decrease in the phase synchronization index. The experimental results are reproduced in simulations which demonstrate that the decrease in the degree of synchronization is correlated with the appearance of a small fraction of permanently excited beads in BZ populations of high mean frequency and/or large width.