The Tris(2,2'-Bipyridyl)Ruthenium-Catalysed Belousov–Zhabotinsky Reaction

Investigating the nonlinear dynamics of photosensitive Belousov-Zhabotinsky gels via bifurcation analyses

Controlling the dynamics of active stimuli-responsive smart materials is essential to replicate the biomimetic functionalities at different length scales for a variety of biological systems-based applications. Photosensitive Belousov-Zhabotinsky (BZ) gels, powered by a nonlinear chemical oscillator, called a BZ reaction are one of the stimuli-responsive smart materials in demand due to their ability to continuously transduce chemical oscillations into mechanical deformations. The chemical oscillations in a BZ reaction and subsequent mechanical oscillations in photosensitive BZ gels occur due to the redox cycle of photosensitive ruthenium complex-based catalysts. In this work, our objective is to identify how the behavior of photosensitive BZ gels can be tuned and used for biomimetic applications by investigating its dynamical characteristics using bifurcation analyses. Specifically, we use the normal form approach and perform linear and nonlinear stability analyses to identify high-order bifurcations by computing higher-order Lyapunov and frequency coefficients. We revealed the existence of domains that encompass coexisting stable and unstable limit cycles (LCs), which merge to form a semi-stable LC at the limit point of cycle (LPC). Their existence shows how a slight variation in the BZ gel recipe can significantly alter its dynamics. Subsequently, we quantify the amplitude and frequency of oscillations in different domains under the effect of variation of BZ reaction formulations. We believe that the outcomes of our work serve as an efficient template for the design and control of BZ gel-based applications. The usage of a normal form and a systematic representation of nonlinear dynamics allow our framework to be extended for other nonlinear dynamical systems.

Experimental studies of spiral wave teleportation in a light sensitive Belousov-Zhabotinsky system

Cardiac arrythmias are a form of heart disease that contributes toward making heart disease a significant cause of death globally. Irregular rhythms associated with cardiac arrythmias are thought to arise due to singularities in the heart tissue that generate reentrant waves in the underlying excitable medium. A normal approach to removing such singularities is to apply a high voltage electric shock, which effectively resets the phase of the cardiac cells. A concern with the use of this defibrillation technique is that the high-energy shock can cause lasting damage to the heart tissue. Various theoretical works have investigated lower-energy alternatives to defibrillation. In this work, we demonstrate the effectiveness of a low-energy defibrillation method in an experimental 2D Belousov-Zhabotinsky (BZ) system. When implemented as a 2D spatial reaction, the BZ reaction serves as an effective analog of general excitable media and supports regular and reentrant wave activity. The defibrillation technique employed involves targeted low-energy perturbations that can be used to "teleport" and/or annihilate singularities present in the excitable BZ medium.

Co-Complexes-Based Self-Oscillating Gels Driven by the Belousov-Zhabotinsky Reaction

We report the synthesis of novel cobalt complexes-based catalysts designed for the oscillatory Belousov-Zhabotinsky (BZ) reaction. For the first time, we introduce cobalt complex-based self-oscillating gels that demonstrate autonomous color oscillations within a BZ reagent solution, functioning without the need for any external stimuli. We created acrylamide-based self-oscillating gels containing immobilized tris(2,2'-bipyridine)cobalt(II) or tris(1,10-phenanthroline)cobalt(II) complexes and gels containing covalently bound (5-acrylamido-1,10-phenanthroline)bis(2,2'-bipyridine)cobalt(II), (5-acrylamido-1,10-phenanthroline)bis(1,10-phenanthroline) cobalt(II), or tris(5-acrylamido-1,10-phenanthroline)cobalt(II) complexes. When the BZ reaction takes place within the gels, it results in the observation of moving chemical waves and reversible color changes. We believe that Co-complexes-based self-oscillating gels have potential applications in the design of soft actuators and chemical devices for signal processing.

Fabrication of submillimeter-sized spherical self-oscillating gels and control of their isotropic volumetric oscillatory behaviors

In this study, we established a fabrication method and analyzed the volumetric self-oscillatory behaviors of submillimeter-sized spherical self-oscillating gels. We validated that the manufactured submillimeter-sized spherical self-oscillating gels exhibited isotropic volumetric oscillations during the Belousov-Zhabotinsky (BZ) reaction. In addition, we experimentally elucidated that the volumetric self-oscillatory behaviors (i.e., period and amplitude) and the oscillatory profiles depended on the following parameters: (1) the molar composition of N-(3-aminopropyl)methacrylamide hydrochloride (NAPMAm) in the gels and (2) the concentration of Ru(bpy)₃-NHS solution containing an active ester group on conjugation. These clarified relationships imply that controlling the amount of Ru(bpy)₃ in the gel network could influence the gel volumetric oscillation during the BZ reaction. These results of submillimeter-sized and spherical self-oscillating gels bridge knowledge gaps in the current field because the gels with corresponding sizes and shapes have not been systematically explored yet. Therefore, our study could be a cornerstone for diverse applications of (self-powered) gels in various scales and shapes, including soft actuators exhibiting life-like functions.

Plasticity in networks of active chemical cells with pulse coupling

A method for controlling the coupling strength is proposed for pulsed coupled active chemical micro-cells. The method is consistent with Hebb's rules. The effect of various system parameters on this "spike-timing-dependent plasticity" is studied. In addition to networks of two and three coupled active cells, the effect of this "plasticity" on the dynamic modes of a network of four pulse-coupled chemical micro-cells unidirectionally coupled in a circle is studied. It is shown that the proposed adjustment of the coupling strengths leads to spontaneous switching between network eigenmodes.

Chirality observed in a driven ruthenium-catalyzed Belousov-Zhabotinsky reaction model

Chirality is commonly associated with the spatial geometry of the atoms composing molecules, the biochemistry of living organisms, and spin properties. In sharp contrast, here we report chirality found in numerically computed stability diagrams of a chemical reaction governed by purely classical (that is, not quantum) equations, namely in a photochemically periodically perturbed ruthenium-catalyzed Belousov-Zhabotinsky reaction model. This novel chirality offers opportunities to explore hitherto unsuspected properties of purely classical chemical oscillators.

Chemical micro-oscillators based on the Belousov–Zhabotinsky reaction

The results of studies on the development of micro-oscillators (MOs) based on the Belousov –Zhabotinsky (BZ) oscillatory chemical reaction are integrated and systematized. The mechanisms of the BZ reaction and the methods of immobilization of the catalyst of the BZ reaction in micro-volumes are briefly discussed. Methods for creating BZ MOs based on water microdroplets in the oil phase and organic and inorganic polymer microspheres are considered. Methods of control and management of the dynamics of BZ MO networks are described, including methods of MO synchronization. The prospects for the design of neural networks of MOs with intelligent-like behaviour are outlined. Such networks present a new area of nonlinear chemistry, including, in particular, the creation of a chemical ‘computer’.

The bibliography includes 250 references.

Novel modes of synchronization in star networks of coupled chemical oscillators

Photochemically coupled micro-oscillators are studied experimentally and computationally in star networks to investigate the modes and mechanisms of synchronization. The micro-oscillators are catalyst-loaded beads that are placed in catalyst-free Belousov-Zhabotinsky (BZ) solutions. The properties of the photochemical coupling between the oscillators are determined by the composition of the BZ reaction mixtures, and both excitatory coupling and inhibitory coupling are studied. Synchronization of peripheral oscillators coupled through a hub oscillator is exhibited at coupling strengths leading to novel modes of synchronization of the hub with the peripheral oscillators. A theoretical analysis provides insights into the mechanism of the synchronization. The heterogeneous peripheral oscillators have different phase velocities that give rise to a phase divergence; however, the perturbation from the hub acts to realign the phases by delaying the faster oscillators more than the slower oscillators.

Quint points lattice in a driven Belousov-Zhabotinsky reaction model

We report the discovery of a regular lattice of exceptional quint points in a periodically driven oscillator, namely, in the frequency-amplitude control parameter space of a photochemically periodically perturbed ruthenium-catalyzed Belousov-Zhabotinsky reaction model. Quint points are singular boundary points where five distinct stable oscillatory phases coalesce. While spikes of the activator show a smooth and continuous variation, the spikes of the inhibitor show an intricate but regular branching into a myriad of stable phases that have fivefold contact points. Such boundary points form a wide parameter lattice as a function of the frequency and amplitude of light absorption. These findings revise current knowledge about the topology of the control parameter space of a celebrated prototypical example of an oscillating chemical reaction.

Covalent Post-Assembly Modification: A Synthetic Multipurpose Tool in Supramolecular Chemistry

The use of covalent post-assembly modification (PAM) in supramolecular chemistry has grown significantly in recent years, to the point where PAM is now a versatile synthesis tool for tuning, modulating, and expanding the functionality of self-assembled complexes and materials. PAM underpins supramolecular template-synthesis strategies, enables modular derivatization of supramolecular assemblies, permits the covalent 'locking' of unstable structures, and can trigger controlled structural transformations between different assembled morphologies. This Review discusses key examples of PAM spanning a range of material classes, including discrete supramolecular complexes, self-assembled soft nanostructures and hierarchically ordered polymeric and framework materials. As such, we hope to highlight how PAM has continued to evolve as a creative and functional addition to the synthetic chemist's toolbox for constructing bespoke self-assembled complexes and materials.

Transition from spiral wave chimeras to phase cluster states

Photochemically coupled Belousov-Zhabotinsky micro-oscillators are studied in experiments and simulations. Generally good agreement between the experimental and simulated dynamical behavior is found, with spiral wave chimeras exhibited at small values of the time delay in the coupling between the oscillators, spiral wave core splitting at higher values, and phase cluster states replacing the spiral wave dynamics at the highest values of the time delay. Spiral wave chimera dynamics is exhibited experimentally for much of the time delay range, while spiral wave phase cluster states are exhibited more in the model simulations. In addition to comparing the experimental and simulation behavior, we explore the novel spiral wave phase cluster states and develop a mechanism for this new and unusual dynamical behavior.

The Belousov-Zhabotinsky Reaction in Thermoresponsive Core-Shell Hydrogel Microspheres with a Tris(2,2'-bipyridyl)ruthenium Catalyst in the Core

The Belousov-Zhabotinsky (BZ) reaction shows temporal or spatiotemporal structures such as redox oscillation of the catalyst, [ruthenium(II)tris(2,2'-bipyridine)][PF₆]₂ ([Ru(bpy)₃][PF₆]₂). In this study, autonomously oscillating hydrogel microspheres (microgels) were investigated, which show swelling/deswelling oscillation induced by the redox oscillation of the BZ reaction inside the gel. Despite the periodically and autonomously induced oscillation that does not require an external stimulus, it has not been possible to perform any manipulation of the oscillatory behavior over time. The results of the present study show that it is possible to reversibly switch the microgel oscillations from an "on" active state of the BZ reaction to an "off" inactive state by changing the temperature in combination with thermoresponsive microgels. To realize on-demand switching, the construction of double-shell structures is crucial; the thermoresponsive first shell allows the microgels to modulate the diffusion of the substrates or intermediates in the BZ reaction, while the second shell maintains colloidal stability under high temperatures and high ion concentrations. The functionalized double-shell microgels were prepared via multistep seeded precipitation polymerization. The oscillatory switching behavior of the BZ reaction was observed directly and evaluated by ultraviolet-visible (UV-vis) spectroscopy. The central concept of this study, i.e., "on-off switching" can be expected to benefit the development of advanced bioinspired materials.

Measurements and Simulations of the Acidity Dependence of the Kinetics of the Iron-Catalyzed Belousov-Zhabotinsky Reaction: Proton-Catalysis in the Electron Transfer Reaction Involving the [Fe(phen)3]3+ Species

The acidity dependence of the iron-catalyzed bromate-malonic acid Belousov-Zhabotinsky reaction was studied in the range 0.36 M < [H₂SO₄]₀ < 1.20 M, and the temporal evolutions of the oscillation patterns were analyzed. The experimental results show that the period times PT _i decrease exponentially with increasing acidity and that the period times parallel the decrease of the reduction times RT with increasing acidity. Simulations using the reactions of the commonly accepted core reaction mechanism failed to match the measurements even in a qualitative fashion. However, we found that compelling agreement between the experiments and the simulations can be achieved over the entire range with the inclusion of second-order proton-catalysis of the oxidation of bromomalonic acid (BrMA) by the [Fe(phen)₃]³⁺ species in the reaction identified in this paper as reaction 9 (R9), and this [H⁺] dependence is informative about the species involved in the outer sphere electron transfer reaction. The trication [Fe(phen)₃]³⁺ species is stabilized by ion pairing and solvation, and one may anticipate the presence of [Fe(phen)₃(HSO₄) _n(H₂O) _m]^{(3- n)+} species ( n = 0-3). Our results suggest that the removal of aggregating HSO₄^- ions by protonation creates a better oxidant and facilitates the approach of the reductant BrMA, and the second-order [H⁺] dependence further suggests that BrMA is primarily oxidized by a doubly charged [Fe(phen)₃(HSO₄)₁(L) _k]²⁺ species. Considering the complexity of the BZ system and the uncertainties in the many reaction rate constants, we were somewhat surprised to find this high level of agreement by (just) the replacement of R9 by R9'. In fact, the near-quantitative agreement presents a powerful corroboration of the core reaction mechanism of the BrMA-rich BZ reaction, and the replacement of R9 by R9' extends the validity of this core reaction mechanism to acidities above and below the typical acidity of BZ reactions ([H⁺] ≈ 1 M).

Autonomous cycling between excitatory and inhibitory coupling in photosensitive chemical oscillators

Photochemically coupled Belousov-Zhabotinsky micro-oscillators are studied in experiments and simulations. The photosensitive oscillators exhibit excitatory or inhibitory coupling depending on the surrounding reaction mixture composition, which can be systematically varied. In-phase or out-of-phase synchronization is observed with predominantly excitatory or inhibitory coupling, respectively, and complex frequency cycling between excitatory and inhibitory coupling is found between these extremes. The dynamical behavior is characterized in terms of the corresponding phase response curves, and a map representation of the dynamics is presented.

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.

Desynchronization of stochastically synchronized chemical oscillators

Experimental and theoretical studies are presented on the design of perturbations that enhance desynchronization in populations of oscillators that are synchronized by periodic entrainment. A phase reduction approach is used to determine optimal perturbation timing based upon experimentally measured phase response curves. The effectiveness of the perturbation waveforms is tested experimentally in populations of periodically and stochastically synchronized chemical oscillators. The relevance of the approach to therapeutic methods for disrupting phase coherence in groups of stochastically synchronized neuronal oscillators is discussed.

Creation and perturbation of planar networks of chemical oscillators

Methods for creating custom planar networks of diffusively coupled chemical oscillators and perturbing individual oscillators within the network are presented. The oscillators consist of the Belousov-Zhabotinsky (BZ) reaction contained in an emulsion. Networks of drops of the BZ reaction are created with either Dirichlet (constant-concentration) or Neumann (no-flux) boundary conditions in a custom planar configuration using programmable illumination for the perturbations. The differences between the observed network dynamics for each boundary condition are described. Using light, we demonstrate the ability to control the initial conditions of the network and to cause individual oscillators within the network to undergo sustained period elongation or a one-time phase delay.

Insights into collective cell behaviour from populations of coupled chemical oscillators

Biological systems such as yeast show coordinated activity driven by chemical communication between cells. Experiments with coupled chemical oscillators can provide insights into the collective behaviour.

Link weight evolution in a network of coupled chemical oscillators

Link weight evolution is studied in a network of coupled chemical oscillators. Oscillators are perturbed by adjustments in imposed light intensity based on excitatory or inhibitory links to other oscillators undergoing excitation. Experimental and modeling studies demonstrate that the network is capable of producing sustained coordinated activity. The individual nodes of the network exhibit incoherent firing events; however, a dominant frequency can be discerned within the collective signal by Fourier analysis. The introduction of spike-timing-dependent plasticity yields a network that evolves to a stable unimodal link weight distribution.

Disproportionation of bromous acid HOBrO by direct O-transfer and via anhydrides O(BrO)2 and BrO-BrO2. An ab initio study of the mechanism of a key step of the Belousov-Zhabotinsky oscillating reaction

The results are reported of an ab initio study of the thermochemistry and of the kinetics of the HOBrO disproportionation reaction 2HOBrO (2) ⇄ HOBr (1) + HBrO(3) (3), reaction ( R4' ), in gas phase (MP2(full)/6-311G*) and aqueous solution (SMD(MP2(full)/6-311G*)). The reaction energy of bromous acid disproportionation is discussed in the context of the coupled reaction system R2-R4 of the FKN mechanism of the Belousov-Zhabotinsky reaction and considering the acidities of HBr and HOBrO(2). The structures were determined of ten dimeric aggregates 4 of bromous acid, (HOBrO)(2), of eight mixed aggregates 5 formed between the products of disproportionation, (HOBr)(HOBrO(2)), and of four transition states structures 6 for disproportionation by direct O-transfer. It was found that the condensation of two HOBrO molecules provides facile access to bromous acid anhydride 7, O(BrO)(2). A discussion of the potential energy surface of Br(2)O(3) shows that O(BrO)(2) is prone to isomerization to the mixed anhydride 8, BrO-BrO(2), and to dissociation to 9, BrO, and 10, BrO(2), and their radical pair 11. Hence, three possible paths from O(BrO)(2) to the products of disproportionation, HOBr and HOBrO(2), are discussed: (1) hydrolysis of O(BrO)(2) along a path that differs from its formation, (2) isomerization of O(BrO)(2) to BrO-BrO(2) followed by hydrolysis, and (3) O(BrO)(2) dissociation to BrO and BrO(2) and their reactions with water. The results of the potential energy surface analysis show that the rate-limiting step in the disproportionation of HOBrO consists of the formation of the hydrate 12a of bromous acid anhydride 7 via transition state structure 14a. The computed activation free enthalpy ΔG(act)(SMD) = 13.6 kcal/mol for the process 2·2a → [14a](‡) → 12a corresponds to the reaction rate constant k(4) = 667.5 M(-1) s(-1) and is in very good agreement with experimental measurements. The potential energy surface analysis further shows that anhydride 7 is kinetically and thermodynamically unstable with regard to hydrolysis to HOBr and HOBrO(2) via transition state structure 14b. The transition state structure 14b is much more stable than 14a, and, hence, the formation of the "symmetrical anhydride" from bromous acid becomes an irreversible reaction for all practical purposes because 7 will instead be hydrolyzed as a "mixed anhydride" to afford HOBr and HOBrO(2). The mixed anhydride 8, BrO-BrO(2), does not play a significant role in bromous acid disproportionation.