Pattern Formation in the Bromate-Sulfite-Ferrocyanide Reaction

Recent advances in the temporal and spatiotemporal dynamics induced by bromate–sulfite-based pH-oscillators

The bromate–sulfite reaction-based pH-oscillators represent one of the most useful subgroup among the chemical oscillators. They provide strong H+-pulses which can generate temporal oscillations in other systems coupled to them and they show wide variety of spatiotemporal dynamics when they are carried out in different gel reactors. Some examples are discussed. When pH-dependent chemical and physical processes are linked to a bromate–sulfite-based oscillator, rhythmic changes can appear in the concentration of some cations and anions, in the distribution of the species in a pH-sensitive stepwise complex formation, in the oxidation number of the central cation in a chelate complex, in the volume or the desorption-adsorption ability of a piece of gel. These reactions are quite suitable for generating spatiotemporal patterns in open reactors. Many reaction–diffusion phenomena, moving and stationary patterns, have been recently observed experimentally using different reactor configurations, which allow exploring the effect of different initial and boundary conditions. Here, we summarize the most relevant aspects of these experimental and numerical studies on bromate–sulfite reaction-based reaction–diffusion systems.

Insights from chemical systems into Turing-type morphogenesis

In 1952, Alan Turing proposed a theory showing how morphogenesis could occur from a simple two morphogen reaction-diffusion system [Turing, A. M. (1952) Phil. Trans. R. Soc. Lond. A 237, 37-72. (doi:10.1098/rstb.1952.0012)]. While the model is simple, it has found diverse applications in fields such as biology, ecology, behavioural science, mathematics and chemistry. Chemistry in particular has made significant contributions to the study of Turing-type morphogenesis, providing multiple reproducible experimental methods to both predict and study new behaviours and dynamics generated in reaction-diffusion systems. In this review, we highlight the historical role chemistry has played in the study of the Turing mechanism, summarize the numerous insights chemical systems have yielded into both the dynamics and the morphological behaviour of Turing patterns, and suggest future directions for chemical studies into Turing-type morphogenesis. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

A simple hydrogel device with flow-through channels to maintain dissipative non-equilibrium phenomena

The development of autonomous chemical systems that could imitate the properties of living matter, is a challenging problem at the meeting point of materials science and nonequilibrium chemistry. Here we design a multi-channel gel reactor in which out-of-equilibrium conditions are maintained by antagonistic chemical gradients. Our device is a rectangular hydrogel with two or more channels for the flows of separated reactants, which diffuse into the gel to react. The relative position of the channels acts as geometric control parameters, while the concentrations of the chemicals in the channels and the variable composition of the hydrogel, which affects the diffusivity of the chemicals, can be used as chemical control parameters. This flexibility allows finding easily the optimal conditions for the development of nonequilibrium phenomena. We demonstrate this straightforward operation by generating diverse spatiotemporal patterns in different chemical reactions. The use of additional channels can create interacting reaction zones.

Designing Stationary Reaction-Diffusion Patterns in pH Self-Activated Systems

Since Alan Turing's 1952 pioneering work, reaction-diffusion (RD) processes are regarded as prototype mechanisms for pattern formation in living systems. Though suspected in many aspects of morphogenetic development, pure RD patterns have not yet been demonstrated in living organisms. The first observations of an autonomous development of stationary chemical patterns were made in the early 1990s. In this Account, we discuss the recent developments for producing stationary pH RD patterns in open spatial reactors. The theoretical analysis of the early experiments anticipated the possibility of finding Turing patterns in a wide range of oscillatory reactions if one could control the kinetic and diffusional rate of some key species. However, no experimentally effective method to produce stationary Turing patterns was attained before 2009, and the number of systems stagnated at two until then. The two precursor reaction systems benefited from unplanned favorable chemical properties of the RD media. Theoretical studies point out that appropriate diffusion rate differences are necessary to produce stationary patterns since a competition between an effective short distance self-activation and a long distance inhibitory process is required. This differential diffusion would naturally lead to differential exchange rates between the RD system and its feed environment, an aspect somewhat overlooked in theoretical and in primal experimental approaches. Our pattern design method takes this aspect into account. A slower diffusion of a self-activated species (here, protons), produced in the RD part of the spatial reactor, generates the accumulation of this species compared to the other species. This accumulation has to be at least partly compensated by an independent scavenging reaction. The above requirement naturally brought us to focus on two-substrate pH oscillatory reactions. Stationary RD patterns are now well documented in six pH driven reaction systems. Furthermore, the coupling with a pH dependent metal ion complexing agent led to stationary patterns in calcium ion concentration. Our effective semiempirical design method does not require a detailed knowledge of the reaction kinetics; thus it is applicable to a broad spectrum of reactions and even to synthetic biological systems. It is based on simple dynamic arguments and on general topological characteristics of a nonequilibrium phase diagram. We first illustrate our method with numerical simulations, based on a realistic but idealized general model of the two-substrate pH-oscillator reaction family, and provide a refined view of the topology of the resulting phase diagrams. Then, we exemplify its effectiveness by observations made in distinct pH self-activated systems. Analogies and differences between experiments and the model calculations are pointed out. Besides standard hexagonal arrays of spots and parallel stripes, hitherto undocumented dynamic phenomena, such as randomly blinking areas and complex dynamic and stationary filamentous structures, were observed. The main challenge, to find low-mobility complexing agents that would selectively and reversibly bind a species controlling the self-activatory kinetic path of the reaction, was readily overcome in multiple ways by anions of weak acids: not only by polymeric substances but in some cases by a pH color indicator or even smaller molecules, depending on their proton binding affinity.

Design of localized spatiotemporal pH patterns by means of antagonistic chemical gradients

Spatially localized moving and stationary pH patterns are generated in two-side-fed reaction-diffusion systems. The patterns are sandwiched between two quiescent zones and positioned by the antagonistic gradients of the reactants of the self-activatory process. Spatial bistability, spatiotemporal oscillations, and formation of stationary Turing patterns have been predicted by numerical simulations and observed in experiments performed by using different hydrogen ion autocatalytic chemical systems. The formation of stationary patterns due to long-range inhibition is promoted by a large molecular weight hydrogen ion binding polymer.

Spatially localized moving and stationary pH patterns are generated in two-side-fed reaction-diffusion systems.

Determination of the diffusion coefficient of hydrogen ion in hydrogels

The role of diffusion in chemical pattern formation has been widely studied due to the great diversity of patterns emerging in reaction-diffusion systems, particularly in H⁺-autocatalytic reactions where hydrogels are applied to avoid convection. A custom-made conductometric cell is designed to measure the effective diffusion coefficient of a pair of strong electrolytes containing sodium ions or hydrogen ions with a common anion. This together with the individual diffusion coefficient for sodium ions, obtained from PFGSE-NMR spectroscopy, allows the determination of the diffusion coefficient of hydrogen ions in hydrogels. Numerical calculations are also performed to study the behavior of a diffusion-migration model describing ionic diffusion in our system. The method we present for one particular case may be extended for various hydrogels and diffusing ions (such as hydroxide) which are relevant e.g. for the development of pH-regulated self-healing mechanisms and hydrogels used for drug delivery.

Modeling Studies of Inhomogeneity Effects during Laser Flash Photolysis Experiments: A Reaction-Diffusion Approach

This work presents a rigorous mathematical study of the effect of unavoidable inhomogeneities in laser flash photolysis experiments. There are two different kinds of inhomegenities: the first arises from diffusion, whereas the second one has geometric origins (the shapes of the excitation and detection light beams). Both of these are taken into account in our reported model, which gives rise to a set of reaction-diffusion type partial differential equations. These equations are solved by a specially developed finite volume method. As an example, the aqueous reaction between the sulfate ion radical and iodide ion is used, for which sufficiently detailed experimental data are available from an earlier publication. The results showed that diffusion itself is in general too slow to influence the kinetic curves on the usual time scales of laser flash photolysis experiments. However, the use of the absorbances measured (e.g., to calculate the molar absorption coefficients of transient species) requires very detailed mathematical consideration and full knowledge of the geometrical shapes of the excitation laser beam and the separate detection light beam. It is also noted that the usual pseudo-first-order approach to evaluating the kinetic traces can be used successfully even if the usual large excess condition is not rigorously met in the reaction cell locally.

Kinetic and Diffusion-Driven Instabilities in the Bromate-Sulfite-Ferrocyanide System

The spatiotemporal dynamics of the bromate-sulfite-ferrocyanide (BSF) reaction-diffusion system in a open one-side-fed reactor (OSFR) is investigated by numerical simulations. The results of the simulations are compared with experiments performed in an annular shape OSFR. Both kinetic and diffusion-driven instabilities are identified in the model. There are two hydrogen ion consuming pathways in the mechanism: the partial oxidation of sulfite to dithionate and the oxidation of ferrocyanide by bromate ions. Their dynamical effects are similar, as they support the same negative feedback loop via sulfite ion. However, the time scale of the oxidation of ferrocyanide by bromate ions can be conveniently controlled by the input feed concentrations, thus it provides a more flexible way to find spatiotemporal oscillations. Long-range activation due to the relative fast diffusion of hydrogen ions compared to the other reactants can also result in oscillations in this mechanism. We show that the spatial extent of the reaction-diffusion medium along the direction of the diffusive feed (the thickness) acts as a general control parameter of the dynamics. Oscillations, either originated in kinetic or in diffusive instabilities, can only develop in a narrow range of the thickness. This property explains the experimentally often observed spatial localization of the oscillations. A reciprocal relationship is found between two main control parameters of the dynamics, which are the thickness and the hydrogen ion input feed concentration.