High-Frequency Oscillations in the Belousov−Zhabotinsky Reaction

Cold-induced degradation of core clock proteins implements temperature compensation in the Arabidopsis circadian clock

The period of circadian clocks is maintained at close to 24 hours over a broad range of physiological temperatures due to temperature compensation of period length. Here, we show that the quantitative control of the core clock proteins TIMING OF CAB EXPRESSION 1 [TOC1; also known as PSEUDO-RESPONSE REGULATOR 1 (PRR1)] and PRR5 is crucial for temperature compensation in Arabidopsis thaliana. The prr5 toc1 double mutant has a shortened period at higher temperatures, resulting in weak temperature compensation. Low ambient temperature reduces amounts of PRR5 and TOC1. In low-temperature conditions, PRR5 and TOC1 interact with LOV KELCH PROTEIN 2 (LKP2), a component of the E3 ubiquitin ligase Skp, Cullin, F-box (SCF) complex. The lkp2 mutations attenuate low temperature-induced decrease of PRR5 and TOC1, and the mutants display longer period only at lower temperatures. Our findings reveal that the circadian clock maintains its period length despite ambient temperature fluctuations through temperature- and LKP2-dependent control of PRR5 and TOC1 abundance.

Artificial temperature-compensated biological clock using temperature-sensitive Belousov-Zhabotinsky gels

The circadian rhythm is a fundamental physiological function for a wide range of organisms. The molecular machinery for generating rhythms has been elucidated over the last few decades. Nevertheless, the mechanism for temperature compensation of the oscillation period, which is a prominent property of the circadian rhythm, is still controversial. In this study, we propose a new mechanism through a chemically synthetic approach (i.e., we realized temperature compensation by the Belousov-Zhabotinsky (BZ) gels). The BZ gels are prepared by embedding a metal catalyst of the BZ reaction into the gel polymer. We made the body of BZ gels using a temperature-sensitive polymer gel, which enabled temperature compensation of the oscillation by using temperature dependence of volume. Moreover, we constructed a simple mathematical model for the BZ oscillation in temperature-sensitive gels. The model can reproduce temperature compensation of BZ gels, even though all reactions are temperature sensitive according to the Arrhenius rule. Our finding hints that a soft body coupling may be underlying temperature-compensated biological functions, including circadian rhythms.

Pattern formation in a four-ring reaction-diffusion network with heterogeneity

In networks of nonlinear oscillators, symmetries place hard constraints on the system that can be exploited to predict universal dynamical features and steady states, providing a rare generic organizing principle for far-from-equilibrium systems. However, the robustness of this class of theories to symmetry-disrupting imperfections is untested in free-running (i.e., non-computer-controlled) systems. Here, we develop a model experimental reaction-diffusion network of chemical oscillators to test applications of the theory of dynamical systems with symmeries in the context of self-organizing systems relevant to biology and soft robotics. The network is a ring of four microreactors containing the oscillatory Belousov-Zhabotinsky reaction coupled to nearest neighbors via diffusion. Assuming homogeneity across the oscillators, theory predicts four categories of stable spatiotemporal phase-locked periodic states and four categories of invariant manifolds that guide and structure transitions between phase-locked states. In our experiments, we observed that three of the four phase-locked states were displaced from their idealized positions and, in the ensemble of measurements, appeared as clusters of different shapes and sizes, and that one of the predicted states was absent. We also observed the predicted symmetry-derived synchronous clustered transients that occur when the dynamical trajectories coincide with invariant manifolds. Quantitative agreement between experiment and numerical simulations is found by accounting for the small amount of experimentally determined heterogeneity in intrinsic frequency. We further elucidate how different patterns of heterogeneity impact each attractor differently through a bifurcation analysis. We show that examining bifurcations along invariant manifolds provides a general framework for developing intuition about how chemical-specific dynamics interact with topology in the presence of heterogeneity that can be applied to other oscillators in other topologies.

Fast-Moving Self-Propelled Droplets of a Nanocatalyzed Belousov-Zhabotinsky Reaction

Self-sustained locomotion by virtue of an internalized chemical reaction is a characteristic feature of living systems and has inspired researchers to develop such self-moving biomimetic systems. Here, we harness a self-oscillating Belousov-Zhabotinsky (BZ) reaction, a well-known chemical oscillator, with enhanced kinetics by virtue of our graphene-based catalytic mats, to elucidate the spontaneous locomotion of BZ reaction droplets. Specifically, our nanocatalysts comprise ruthenium nanoparticle decorations on graphene oxide, reduced graphene oxide, and graphene nanosheets, thereby creating 0D-2D heterostructures. We demonstrate that when these nanocatalyzed droplets of the BZ reaction are placed in an oil-surfactant medium, they exhibit a macroscopic translatory motion at the velocities of few millimeters per second. This motion is brought about by the combination of enhanced kinetics of the BZ reaction and the Marangoni effect. Our investigations reveal that the velocity of locomotion increases with the electrical conductivity of our nanocomposites. Moreover, we also show that the positive feedback generated by the reaction-diffusion phenomena results in an asymmetric distribution of surface tension that, in turn, facilitates the self-propelled motion of the BZ droplets. Finally, we explore a system of multiple nanocatalyzed BZ droplets and reveal a variety of motions caused by their mutual interactions. Our findings suggest that through the use of 0D-2D hybrid nanomaterials, it is possible to design fast-moving self-propelled synthetic objects for a variety of biomimetic applications.

Dynamical attributes of nanocatalyzed self-oscillating reactions via bifurcation analyses

Self-oscillating chemical reactions that undergo reaction-diffusion (RD) phenomena have shown great potential for designing stimuli-responsive materials. Belousov-Zhabotinsky (BZ) reactions are one such class of reactions that exhibit nonlinear chemical oscillations due to redox cycles of the metal-ion catalyst by virtue of Hopf bifurcation. Using bifurcation analyses, here we investigate the BZ reactions, catalyzed by 0D-2D catalytic nanomats and bare nanosheets, which are known to exhibit enhanced dynamic response due to catalysts' heterogeneity. Specifically, we incorporate the nanocatalysts' activity in the kinetic model of the BZ reactions and, subsequently, use catalysts' activity as the bifurcation parameter for analyses. By computing higher-order Lyapunov and frequency coefficients, we have revealed new oscillatory regimes in the bifurcation diagram, including re-entrant regions where sustained oscillations are unexpectedly suppressed, even with high catalytic activity. In addition, we also calculate the amplitude and frequency of BZ oscillations in each of these regions as a function of nanocatalysts' activity. We believe that our current findings can be used to harness the nonlinearity of RD-based dynamical systems to provide unique functionalities to active stimuli-response systems.

Queiroz JP, Faria RB. The Chlorate-Nitrous Acid-Iodine-Iodide Oscillating Reaction

This work presents a new oscillating reaction based on chlorate and observed in a CSTR at room temperature. This can be the first member of a new family of oscillating reactions. In addition, it is also the first oscillating reaction to use nitrous acid as a reactant. Four different behaviors were observed: simple oscillations, mixed mode oscillations, bursts, and quasiperiodicity. The period of oscillations is very short, which is around 1 s. Together with the fact that it also shows fast bursts, it opens the possibility that it can be used to simulate fast biological events, like the neuron's communications signals.

Effect of Reaction Parameters on the Wavelength of Pulse Waves in the Belousov-Zhabotinsky Reaction-Diffusion System

The wavelength of Belousov-Zhabotinsky (BZ) traveling waves is the key factor that limits the scale of BZ self-oscillating gel motors. To achieve control of the wavelength, it is necessary to evaluate the wavelength dependence on species concentrations and temperature. In this work, the effect of reaction parameters on the wavelength of BZ pulse waves was studied. The most effective way to reduce the wavelength of pulse waves is to increase the concentration of organic species and/or the temperature. Decreasing the concentration of bromate, hydrogen ion, or metal catalyst also reduces the wavelength of pulse waves. This work provides a convenient and direct method to produce sub-millimeter BZ waves, which could be applied to designing BZ wave-driven small-scale gel motors as well as providing insight into other emergent behaviors of self-oscillating gels.

Resonant Behavior in a Periodically Forced Nonisothermal Oregonator

Nonisothermal chemical oscillators are poorly studied systems because chemical oscillations are conventionally studied under isothermal conditions. Coupling chemical reactions with heat generation and removal in a nonisothermal oscillatory system can lead to a highly nontrivial nonlinear dynamic behavior. For the current study, we considered the three-variable Oregonator model with the temperature incorporated as a variable (not a parameter), thus adding an energy balance to the set of equations. The effect of temperature on reaction rates is included through the temperature-dependent reaction rate coefficients (Arrhenius law). To model a continuous operation in a laboratory environment, the system was subjected to external forcing through the coolant temperature and infrared irradiation. By conducting numerical simulations and parametric studies, we found that the system is capable of a resonant behavior exhibiting induced oscillations. Our findings indicate that an external source of heat (e.g., via an infrared light emitting diode) can be used to induce a Hopf bifurcation under resonant conditions in an experimental Belousov-Zhabotinsky reactor.

Multivariate statistical analysis of chemical and electrochemical oscillators for an accurate frequency selection

The effect of experimental parameters on the frequency of chemical oscillators has been systematically studied since the first observations of clock reactions. The approach is mainly based on univariate changes in one specific parameter while others are kept constant. The frequency is then monitored and the effect of each parameter is discussed separately. This type of analysis, however, does not take into account the multiple interactions among the controllable parameters and the synergic responses on the oscillation frequency. We have carried out a multivariate statistical analysis of chemical (BZ-ferroin catalyzed reaction) and electrochemical (Cu/Cu₂O cathodic deposition) oscillators and identified the contributions of the experimental parameters on frequency variations. The BZ reaction presented a strong dependence on the initial concentration of sodium bromate and temperature, resulting in a frequency increase. The concentration of malonic acid, the organic substrate, affects the system but with lower intensity compared with the combination of sodium bromate and temperature. On the other hand, the Cu/Cu₂O electrochemical oscillator was shown to be less sensitive to changes in the temperature. The applied current density and pH were the two parameters which most perturbed the system. Interestingly, the frequency behaved nonmonotonically with a quadratic dependence. The multivariate analysis of both oscillators exhibited significant differences - while the homogenous oscillator displayed a linear dependence with the factors, the heterogeneous one revealed a more complex dependence with quadratic terms. Our results may contribute, for instance, in the synthesis of self-organized materials in which an accurate frequency selection is required and, depending on its value, different physicochemical properties are obtained.

Thermal switch of oscillation frequency in Belousov-Zhabotinsky liquid marbles

External control of oscillation dynamics in the Belousov-Zhabotinsky (BZ) reaction is important for many applications including encoding computing schemes. When considering the BZ reaction, there are limited studies dealing with thermal cycling, particularly cooling, for external control. Recently, liquid marbles (LMs) have been demonstrated as a means of confining the BZ reaction in a system containing a solid-liquid interface. BZ LMs were prepared by rolling 50 μl droplets in polyethylene (PE) powder. Oscillations of electrical potential differences within the marble were recorded by inserting a pair of electrodes through the LM powder coating into the BZ solution core. Electrical potential differences of up to 100 mV were observed with an average period of oscillation ca 44 s. BZ LMs were subsequently frozen to -1°C to observe changes in the frequency of electrical potential oscillations. The frequency of oscillations reduced upon freezing to 11 mHz cf. 23 mHz at ambient temperature. The oscillation frequency of the frozen BZ LM returned to 23 mHz upon warming to ambient temperature. Several cycles of frequency fluctuations were able to be achieved.

Identification of two aliphatic position isomers between α- and β-ketoglutaric acid by using a Briggs-Rauscher oscillating system

This paper reports a novel method for identification of two aliphatic position isomers between α-ketoglutaric acid (α-KA) and β-ketoglutaric acid (β-KA) by their different perturbation effects on a Briggs-Rauscher oscillating system, in which tetraaza-macrocyclic complex [NiL](ClO4)2 is used as the catalyst. The ligand L in the complex is 5,7,7,12,14,14-hexamethyl-1,4,8,11-tetraazacyclotetradeca-4,11-diene. The experimental results have shown that addition of α-KA into the system does not affect the oscillating patterns, while the presence of β-KA in a dynamic system influences the oscillatory amplitude. A more interesting feature is that, in the presence of a higher concentration of β-KA, there are damped oscillations after the initial spike, followed by quenching (more exactly: very small oscillations) of the oscillations before the subsequent regeneration of oscillations. A qualitative approach was thus established by employing a Briggs-Rauscher system for identification of these two isomers. The concentrations of these two isomers that can be distinguished lie over the range between 5.0 × 10(-6) and 2.5 × 10(-3) mol/L. A reaction mechanism based on the FCA model has been proposed. An explanation is that β-KA reacts with HOO(•) radicals to form acetone, whereas the α-KA does not.

Noise induced oscillations and coherence resonance in a generic model of the nonisothermal chemical oscillator

Oscillating chemical reactions are common in biological systems and they also occur in artificial non-biological systems. Generally, these reactions are subject to random fluctuations in environmental conditions which translate into fluctuations in the values of physical variables, for example, temperature. We formulate a mathematical model for a nonisothermal minimal chemical oscillator containing a single negative feedback loop and study numerically the effects of stochastic fluctuations in temperature in the absence of any deterministic limit cycle or periodic forcing. We show that noise in temperature can induce sustained limit cycle oscillations with a relatively narrow frequency distribution and some characteristic frequency. These properties differ significantly depending on the noise correlation. Here, we have explored white and colored (correlated) noise. A plot of the characteristic frequency of the noise induced oscillations as a function of the correlation exponent shows a maximum, therefore indicating the existence of autonomous stochastic resonance, i.e. coherence resonance.

Effect of temperature on precision of chaotic oscillations in nickel electrodissolution

We investigate the effects of temperature on complexity features of chaotic electrochemical oscillations using the anodic electrodissolution of nickel in sulfuric acid. The precision of the "period" of chaotic oscillation is characterized by phase diffusion coefficient (D). It is shown that reduced phase diffusion coefficient (D/frequency) exhibits Arrhenius-type dependency on temperature with apparent activation energy of 108 kJ/mol. The reduced Lyapunov exponent of the attractor exhibits no considerable dependency on temperature. These results suggest that the precision of electrochemical oscillations deteriorates with increase in temperature and the variation of phase diffusion coefficient does not necessarily correlate with that of Lyapunov exponent. Modeling studies qualitatively simulate the behavior observed in the experiments: the precision of oscillations in the chaotic Ni dissolution model can be tuned by changes of a time scale parameter of an essential variable, which is responsible for the development of chaotic behavior.