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.

Synchronization scenarios induced by delayed communication in arrays of diffusively coupled autonomous chemical oscillators

We study the impact of delayed feedbacks in the collective synchronization of ensembles of identical and autonomous micro-oscillators. To this aim, we consider linear arrays of Belousov-Zhabotinsky (BZ) oscillators confined in micro-compartmentalised systems, where the delayed feedback mimics natural lags that can arise due to the confinement properties and mechanisms driving the inter-oscillator communication. The micro-oscillator array is modeled as a set of Oregonator-like kinetics coupled via mass exchange of the chemical messengers. Changes in the synchronization patterns are explored by varying the delayed feedback introduced in the messenger species Br2. A direct transition from anti-phase to in-phase synchronization and back to the initial anti-phase scheme is observed by progressively increasing the time delay from zero to the value T₀, which is the oscillation period characterising the system without any delayed coupling. The route from anti- to in-phase oscillations (and back) consists of regimes where windows of in-phase oscillations are periodically broken by anti-phase beats. Similarities between these phase transition dynamics and synchronization scenarios characterising the coordination of oscillatory limb movements are finally discussed.

Distance dependent types of coupling of chemical micro-oscillators immersed in a water-in-oil microemulsion

A system of micro-spheres immersed in a water-in-oil microemulsion (ME) is studied both theoretically and experimentally.

Impact of PDMS-Based Microfluidics on Belousov-Zhabotinsky Chemical Oscillators

Sub-nanoliter volumes of the Belousov-Zhabotinsky (BZ) reaction are sealed in microfluidic devices made from polydimethylsiloxane (PDMS). Bromine, which is a BZ reaction intermediate that participates in the inhibitory pathway of the reaction, is known to permeate into PDMS, and it has been suggested that PDMS and bromine can react ( J. Phys. Chem. A. 108, 2004, 1325-1332). We characterize the extent to which PDMS affects BZ oscillations by varying the volume of the PDMS surrounding the BZ reactors. We measure how the oscillation period varies with PDMS volume and compare with a theoretical reaction-diffusion model, concluding that bromine reacts with PDMS. We demonstrate that minimizing the amount of PDMS by making the samples as thin as possible maximizes the number of oscillations before the BZ reaction reaches equilibrium and ceases to oscillate. We also demonstrate that the deleterious effects of the PDMS-BZ interactions are somewhat mitigated by imposing constant chemical boundary conditions through using a light-sensitive catalyst, ruthenium, in combination with patterned illumination. Furthermore, we show that light can modulate the frequency and phase of the BZ oscillators contained in a PDMS matrix by 20-30%.

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.

Membrane Structure Drives Synchronization Patterns in Arrays of Diffusively Coupled Self-Oscillating Droplets

Networks of diffusively coupled inorganic oscillators, confined in nano- and microcompartments, are effective for predicting and understanding the global dynamics of those systems where the diffusion of activatory or inhibitory signals regulates the communication among different individuals. By taking advantage of a microfluidic device, we study the dynamics of arrays of diffusively coupled Belousov-Zhabotinsky (BZ) oscillators encapsulated in water-in-oil single emulsions. New synchronization patterns are induced and controlled by modulating the structural and chemical properties of the phospholipid-based biomimetic membranes via the introduction of specific dopants. Doping molecules do not alter the membrane basic backbone, but modify the lamellarity (and, in turn, the permeability) or interact chemically with the reaction intermediates. A transition from two-period clusters showing 1:2 period-locking to one-period antiphase synchronization is observed by decreasing the membrane lamellarity. An unsynchronized scenario is found when the dopant is able to interfere with chemical communication by reacting with the chemical messengers.

Microfluidic compartmentalization of diffusively coupled oscillators in multisomes induces a novel synchronization scenario

Multisome compartments encapsulating the Belousov–Zhabotinsky reaction produced by microfluidics arranged in 1D arrays showed a novel type of global synchronization.

Dynamics of Reaction-Diffusion Oscillators in Star and other Networks with Cyclic Symmetries Exhibiting Multiple Clusters

We experimentally and theoretically investigate the dynamics of inhibitory coupled self-driven oscillators on a star network in which a single central hub node is connected to k peripheral arm nodes. The system consists of water-in-oil Belousov-Zhabotinsky ∼100  μm emulsion drops contained in storage wells etched in silicon wafers. We observed three dynamical attractors by varying the number of arms in the star graph and the coupling strength: (i) unlocked, uncorrelated phase shifts between all oscillators; (ii) locked, arm hubs synchronized in phase with a k-dependent phase shift between the arm and central hub; and (iii) center silent, a central hub stopped oscillating and the arm hubs oscillated without synchrony. We compare experiment to theory. For case (ii), we identified a logarithmic dependence of the phase shift on star degree, and were able to discriminate between contributions to the phase shift arising from star topology and oscillator chemistry.

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.

Towards Functional Droplet Architectures: a Belousov-Zhabotinsky Medium for Networks

The confluence of droplet-compartmentalised chemical systems and architectures composed of interacting droplets points towards a novel technology mimicking core features of the cellular architecture that dominates biology. A key challenge to achieve such a droplet technology is long-term stability in conjunction with interdroplet communication. Here, we probed the parameter space of the Belousov-Zhabotinsky (BZ) medium, an extensively studied model for non-equilibrium chemical reactions, pipetted as 2.5 mm droplets in hexadecane oil. The presence of asolectin lipids enabled the formation of arrays of contacted BZ droplets, of which the wave patterns were characterised over time. We utilised laser-cut acrylic templates with over 40 linear oil-filled slots in which arrays are formed by pipetting droplets of the desired BZ composition, enabling parallel experiments and automated image analysis. Using variations of conventional malonic acid BZ medium, wave propagation over droplet-droplet interfaces was not observed. However, a BZ medium containing both malonic acid and 1,4-cyclohexanedione was found to enable inter-droplet wave propagation. We anticipate that the chemical excitation properties of this mixed-substrate BZ medium, in combination with the droplet stability of the networks demonstrated here for nearly 400 droplets in a template-defined topology, will facilitate the development of scalable functional droplet networks.

Engineering reaction-diffusion networks with properties of neural tissue

We present an experimental system of networks of coupled non-linear chemical reactors, which we theoretically model within a reaction-diffusion framework. The networks consist of patterned arrays of diffusively coupled nanoliter-scale reactors containing the Belousov-Zhabotinsky (BZ) reaction. Microfluidic fabrication techniques are developed that provide the ability to vary the network topology and the reactor coupling strength and offer the freedom to choose whether an arbitrary reactor is inhibitory or excitatory coupled to its neighbor. This versatile experimental and theoretical framework can be used to create a wide variety of chemical networks. Here we design, construct and characterize chemical networks that achieve the complexity of central pattern generators (CPGs), which are found in the autonomic nervous system of a variety of organisms.

Dynamics of a 1D array of inhibitory coupled chemical oscillators in microdroplets with global negative feedback

We have investigated the effect of global negative feedback (GNF) on the dynamics of a 1D array of water microdroplets (MDs) filled with the reagents of the photosensitive oscillatory Belousov–Zhabotinsky (BZ) reaction.

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.

Optical Communication among Oscillatory Reactions and Photo-Excitable Systems: UV and Visible Radiation Can Synchronize Artificial Neuron Models

Neuromorphic engineering promises to have a revolutionary impact in our societies. A strategy to develop artificial neurons (ANs) is to use oscillatory and excitable chemical systems. Herein, we use UV and visible radiation as both excitatory and inhibitory signals for the communication among oscillatory reactions, such as the Belousov-Zhabotinsky and the chemiluminescent Orban transformations, and photo-excitable photochromic and fluorescent species. We present the experimental results and the simulations regarding pairs of ANs communicating by either one or two optical signals, and triads of ANs arranged in both feed-forward and recurrent networks. We find that the ANs, powered chemically and/or by the energy of electromagnetic radiation, can give rise to the emergent properties of in-phase, out-of-phase, anti-phase synchronizations and phase-locking, dynamically mimicking the communication among real neurons.

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.

Configurable NOR gate arrays from Belousov-Zhabotinsky micro-droplets

We investigate the Belousov-Zhabotinsky (BZ) reaction in an attempt to establish a basis for computation using chemical oscillators coupled via inhibition. The system consists of BZ droplets suspended in oil. Interdrop coupling is governed by the non-polar communicator of inhibition, Br2. We consider a linear arrangement of three droplets to be a NOR gate, where the center droplet is the output and the other two are inputs. Oxidation spikes in the inputs, which we define to be TRUE, cause a delay in the next spike of the output, which we read to be FALSE. Conversely, when the inputs do not spike (FALSE) there is no delay in the output (TRUE), thus producing the behavior of a NOR gate. We are able to reliably produce NOR gates with this behavior in microfluidic experiment.

Extended source model for diffusive coupling

Motivated by the prevailing approach to diffusion coupling phenomena which considers point-like diffusing sources, we derived an analogous expression for the concentration rate of change of diffusively coupled extended containers. The proposed equation, together with expressions based on solutions to the diffusion equation, is intended to be applied to the numerical solution of systems exclusively composed of ordinary differential equations, however is able to account for effects due the finite size of the coupled sources.

Chimeras in networks with purely local coupling

Chimera states in spatially extended networks of oscillators have some oscillators synchronized while the remainder are asynchronous. These states have primarily been studied in networks with nonlocal coupling, and more recently in networks with global coupling. Here, we present three networks with only local coupling (diffusive, to nearest neighbors) which are numerically found to support chimera states. One of the networks is analyzed using a self-consistency argument in the continuum limit, and this is used to find the boundaries of existence of a chimera state in parameter space.

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.

Tunable diffusive lateral inhibition in chemical cells

The Belousov-Zhabotinsky (BZ) reaction has become the prototype of nonlinear chemical dynamics. Microfluidic techniques provide a convenient method for emulsifying BZ solutions into monodispersed drops with diameters of tens to hundreds of microns, providing a unique system in which reaction-diffusion theory can be quantitatively tested. In this work, we investigate monolayers of microfluidically generated BZ drops confined in close-packed two-dimensional (2D) arrays through experiments and finite element simulations. We describe the transition from oscillatory to stationary chemical states with increasing coupling strength, controlled by independently varying the reaction chemistry within a drop and diffusive flux between drops. For stationary drops, we studied how the ratio of stationary oxidized to stationary reduced drops varies with coupling strength. In addition, using simulation, we quantified the chemical heterogeneity sufficient to induce mixed stationary and oscillatory patterns.

Light-driven dynamic Archimedes spirals and periodic oscillatory patterns of topological solitons in anisotropic soft matter

Oscillatory and excitable systems commonly exhibit formation of dynamic non-equilibrium patterns. For example, rotating spiral patterns are observed in biological, chemical, and physical systems ranging from organization of slime mold cells to Belousov-Zhabotinsky reactions, and to crystal growth from nuclei with screw dislocations. Here we describe spontaneous formation of spiral waves and a large variety of other dynamic patterns in anisotropic soft matter driven by low-intensity light. The unstructured ambient or microscope light illumination of thin liquid crystal films in contact with a self-assembled azobenzene monolayer causes spontaneous formation, rich spatial organization, and dynamics of twisted domains and topological solitons accompanied by the dynamic patterning of azobenzene group orientations within the monolayer. Linearly polarized incident light interacts with the twisted liquid crystalline domains, mimicking their dynamics and yielding patterns in the polarization state of transmitted light, which can be transformed to similar dynamic patterns in its intensity and interference color. This shows that the delicate light-soft-matter interaction can yield complex self-patterning of both. We uncover underpinning physical mechanisms and discuss potential uses.