Information coding with frequency of oscillations in Belousov-Zhabotinsky encapsulated disks

Chemical Memory with Discrete Turing Patterns Appearing in the Glycolytic Reaction

Memory is an essential element in information processing devices. We investigated a network formed by just three interacting nodes representing continuously stirred tank reactors (CSTRs) in which the glycolytic reaction proceeds as a potential realization of a chemical memory unit. Our study is based on the 2-variable computational model of the reaction. The model parameters were selected such that the system has a stable limit cycle and several distinct, discrete Turing patterns characterized by stationary concentrations at the nodes. In our interpretation, oscillations represent a blank memory unit, and Turing patterns code information. The considered memory can preserve information on one of six different symbols. The time evolution of the nodes was individually controlled by the inflow of ATP. We demonstrate that information can be written with a simple and short perturbation of the inflow. The perturbation applies to only one or two nodes, and it is symbol specific. The memory can be erased with identical inflow perturbation applied to all nodes. The presented idea of pattern-coded memory applies to other reaction networks that allow for discrete Turing patterns. Moreover, it hints at the experimental realization of memory in a simple system with the glycolytic reaction.

Neuromorphic Liquids, Colloids, and Gels: A Review

Advances in flexible electronic devices and robotic software require that sensors and controllers be virtually devoid of traditional electronic components, be deformable and stretch-resistant. Liquid electronic devices that mimic biological synapses would make an ideal core component for flexible liquid circuits. This is due to their unbeatable features such as flexibility, reconfiguration, fault tolerance. To mimic synaptic functions in fluids we need to imitate dynamics and complexity similar to those that occurring in living systems. Mimicking ionic movements are considered as the simplest platform for implementation of neuromorphic in material computing systems. We overview a series of experimental laboratory prototypes where neuromorphic systems are implemented in liquids, colloids and gels.

Information Processing Using Networks of Chemical Oscillators

I believe the computing potential of systems with chemical reactions has not yet been fully explored. The most common approach to chemical computing is based on implementation of logic gates. However, it does not seem practical because the lifetime of such gates is short, and communication between gates requires precise adjustment. The maximum computational efficiency of a chemical medium is achieved if the information is processed in parallel by different parts of it. In this paper, I review the idea of computing with coupled chemical oscillators and give arguments for the efficiency of such an approach. I discuss how to input information and how to read out the result of network computation. I describe the idea of top-down optimization of computing networks. As an example, I consider a small network of three coupled chemical oscillators designed to differentiate the white from the red points of the Japanese flag. My results are based on computer simulations with the standard two-variable Oregonator model of the oscillatory Belousov−Zhabotinsky reaction. An optimized network of three interacting oscillators can recognize the color of a randomly selected point with >98% accuracy. The presented ideas can be helpful for the experimental realization of fully functional chemical computing networks.

The Concilium of Information Processing Networks of Chemical Oscillators for Determining Drug Response in Patients With Multiple Myeloma

It can be expected that medical treatments in the future will be individually tailored for each patient. Here we present a step towards personally addressed drug therapy. We consider multiple myeloma treatment with drugs: bortezomib and dexamethasone. It has been observed that these drugs are effective for some patients and do not help others. We describe a network of chemical oscillators that can help to differentiate between non-responsive and responsive patients. In our numerical simulations, we consider a network of 3 interacting oscillators described with the Oregonator model. The input information is the gene expression value for one of 15 genes measured for patients with multiple myeloma. The single-gene networks optimized on a training set containing outcomes of 239 therapies, 169 using bortezomib and 70 using dexamethasone, show up to 71% accuracy in differentiating between non-responsive and responsive patients. If the results of single-gene networks are combined into the concilium with the majority voting strategy, then the accuracy of predicting the patient’s response to the therapy increases to ∼ 85%.

Identification of the best medium for experiments on chemical computation with Belousov–Zhabotinsky reaction and ferroin-loaded Dowex beads

Our study is focused on identification of the best medium for future experiments on information processing with Belousov–Zhabotinsky reaction proceeding in Dowex beads with immobilized catalyst inside. The optimum medium should be characterized by long and stable nonlinear behavior, mechanical stability and should allow for control with electric potential. We considered different types of Dowex ion-exchange resins, bead distributions and various initial concentrations of substrates: malonic acid and 1,4-cyclohexanedione. The electric potential on platinum electrodes, stabilized by a potentiostat is used to control medium evolution. A negative electric potential generates activator species HBrO2 on the working electrode according to the reaction: BrO3−  + 2e−  + 3H+  → HBrO2 + H2O, while positive electric potential attracts inhibitor species Br− to the proximity of it. We study oscillation amplitude and period stability in systems with ferroin loaded Dowex 50W-X2 and Dowex 50W-X8 beads during experiments exceeding 16 h. It has been observed, that the above mentioned resins generate a smaller number of CO2 bubbles close to the beads than Dowex 50W-X4, which makes Dowex 50W-X2 and Dowex 50W-X8 more suitable for applications in chemical computing. We report amplitude stability, oscillation frequency, merging and annihilation of travelling waves in a lattice of Dowex 50W-X8 beads (mesh size 50–100) in over 19 h long experiments with equimolar solution of malonic acid and 1,4-cyclohexanedione. This system looks as a promising candidate for chemical computing devices that can operate for a day.

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.

Towards proteinoid computers. Hypothesis paper

Proteinoids - thermal proteins - are produced by heating amino acids to their melting point and initiation of polymerisation to produce polymeric chains. Proteinoids swell in aqueous solution into hollow microspheres. The proteinoid microspheres produce endogenous burst of electrical potential spikes and change patterns of their electrical activity in response to illumination. The microspheres can interconnect by pores and tubes and form networks with a programmable growth. We speculate on how ensembles of the proteinoid microspheres can be developed into unconventional computing devices.

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.

Advanced Chemical Computing Using Discrete Turing Patterns in Arrays of Coupled Cells

We examine dynamical switching among discrete Turing patterns that enable chemical computing performed by mass-coupled reaction cells arranged as arrays with various topological configurations: three coupled cells in a cyclic array, four coupled cells in a linear array, four coupled cells in a cyclic array, and four coupled cells in a branched array. Each cell is operating as a continuous stirred tank reactor, within which the glycolytic reaction takes place, represented by a skeleton inhibitor-activator model where ADP plays the role of activator and ATP is the inhibitor. The mass coupling between cells is assumed to be operating in three possible transport regimes: (i) equal transport coefficients of the inhibitor and activator (ii) slightly faster transport of the activator, and (iii) strongly faster transport of the inhibitor. Each cellular array is characterized by two pairs of tunable parameters, the rate coefficients of the autocatalytic and inhibitory steps, and the transport coefficients of the coupling. Using stability and bifurcation analysis we identified conditions for occurrence of discrete Turing patterns associated with non-uniform stationary states. We found stable symmetric and/or asymmetric discrete Turing patterns coexisting with stable uniform periodic oscillations. To switch from one of the coexisting stable regimes to another we use carefully targeted perturbations, which allows us to build systems of logic gates specific to each topological type of the array, which in turn enables to perform advanced modes of chemical computing. By combining chemical computing techniques in the arrays with glycolytic excitable channels, we propose a cellular assemblage design for advanced chemical computing.

Experimental verification of an opto-chemical "neurocomputer"

A theoretically predicted hierarchical network of pulse coupled chemical micro-oscillators and excitable micro-cells that we call a chemical "neurocomputer" (CN) or even a chemical "brain" is tested experimentally using the Belousov-Zhabotinsky reaction. The CN consists of five functional units: (1) a central pattern generator (CPG), (2) an antenna, (3) a reader for the CPG, (4) a reader for the antenna unit, and (5) a decision making (DM) unit. A hybrid CN, in which such chemical units as readers and DM units are replaced by electronic units, is tested as well. All these variations of the CN respond intelligently to external signals, since they perform an automatic transition from a current to a new dynamic mode of the CPG, which is similar to the antenna dynamic mode that in turn is induced by external signals. In other words, we show for the first time that a network of pulse coupled chemical micro-oscillators is capable of intelligent adaptive behavior.

Applications of Information Theory Methods for Evolutionary Optimization of Chemical Computers

It is commonly believed that information processing in living organisms is based on chemical reactions. However, the human achievements in constructing chemical information processing devices demonstrate that it is difficult to design such devices using the bottom-up strategy. Here I discuss the alternative top-down design of a network of chemical oscillators that performs a selected computing task. As an example, I consider a simple network of interacting chemical oscillators that operates as a comparator of two real numbers. The information on which of the two numbers is larger is coded in the number of excitations observed on oscillators forming the network. The parameters of the network are optimized to perform this function with the maximum accuracy. I discuss how information theory methods can be applied to obtain the optimum computing structure.

Hierarchical network of pulse coupled chemical oscillators with adaptive behavior: Chemical neurocomputer

We consider theoretically a network of pulse coupled oscillators with time delays. Each oscillator is described by the Oregonator-like model for the Belousov-Zhabotinsky (BZ) reaction. Different groups of oscillators constitute five functional units: (1) a central pattern generator (CPG), (2) a "reader" unit that can identify dynamical modes of the CPG, (3) an antenna (A) unit that receives external signals and responds on them by generating different dynamical modes, (4) another reader unit for identification of the dynamical modes in the A unit, and (5) a decision making unit that switches the current dynamical mode of the CPG to the mode that is similar to the current mode in the A unit. We call this network a chemical neurocomputer, since chemical BZ reaction occurs in each micro-oscillator, while pulse connectivity of these cells is inspired by the brain.

A brief history of liquid computers

A substrate does not have to be solid to compute. It is possible to make a computer purely from a liquid. I demonstrate this using a variety of experimental prototypes where a liquid carries signals, actuates mechanical computing devices and hosts chemical reactions. We show hydraulic mathematical machines that compute functions based on mass transfer analogies. I discuss several prototypes of computing devices that employ fluid flows and jets. They are fluid mappers, where the fluid flow explores a geometrically constrained space to find an optimal way around, e.g. the shortest path in a maze, and fluid logic devices where fluid jet streams interact at the junctions of inlets and results of the computation are represented by fluid jets at selected outlets. Fluid mappers and fluidic logic devices compute continuously valued functions albeit discretized. There is also an opportunity to do discrete operation directly by representing information by droplets and liquid marbles (droplets coated by hydrophobic powder). There, computation is implemented at the sites, in time and space, where droplets collide one with another. The liquid computers mentioned above use liquid as signal carrier or actuator: the exact nature of the liquid is not that important. What is inside the liquid becomes crucial when reaction-diffusion liquid-phase computing devices come into play: there, the liquid hosts families of chemical species that interact with each other in a massive-parallel fashion. I shall illustrate a range of computational tasks, including computational geometry, implementable by excitation wave fronts in nonlinear active chemical medium. The overview will enable scientists and engineers to understand how vast is the variety of liquid computers and will inspire them to design their own experimental laboratory prototypes. This article is part of the theme issue 'Liquid brains, solid brains: How distributed cognitive architectures process information'.

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.

Towards fungal computer

We propose that fungi Basidiomycetes can be used as computing devices: information is represented by spikes of electrical activity, a computation is implemented in a mycelium network and an interface is realized via fruit bodies. In a series of scoping experiments, we demonstrate that electrical activity recorded on fruits might act as a reliable indicator of the fungi's response to thermal and chemical stimulation. A stimulation of a fruit is reflected in changes of electrical activity of other fruits of a cluster, i.e. there is distant information transfer between fungal fruit bodies. In an automaton model of a fungal computer, we show how to implement computation with fungi and demonstrate that a structure of logical functions computed is determined by mycelium geometry.

Street map analysis with excitable chemical medium

Belousov-Zhabotinsky (BZ) thin layer solution is a fruitful substrate for designing unconventional computing devices. A range of logical circuits, wet electronic devices, and neuromorphic prototypes have been constructed. Information processing in BZ computing devices is based on interaction of oxidation (excitation) wave fronts. Dynamics of the wave fronts propagation is programed by geometrical constraints and interaction of colliding wave fronts is tuned by illumination. We apply the principles of BZ computing to explore a geometry of street networks. We use two-variable Oregonator equations, the most widely accepted and verified in laboratory experiments BZ models, to study propagation of excitation wave fronts for a range of excitability parameters, with gradual transition from excitable to subexcitable to nonexcitable. We demonstrate a pruning strategy adopted by the medium with decreasing excitability when wider and ballistically appropriate streets are selected. We explain mechanics of streets selection and pruning. The results of the paper will be used in future studies of studying dynamics of cities and characterizing geometry of street networks.

A 'reader' unit of the chemical computer

We suggest the main principals and functional units of the parallel chemical computer, namely, (i) a generator (which is a network of coupled oscillators) of oscillatory dynamic modes, (ii) a unit which is able to recognize these modes (a 'reader') and (iii) a decision-making unit, which analyses the current mode, compares it with the external signal and sends a command to the mode generator to switch it to the other dynamical regime. Three main methods of the functioning of the reader unit are suggested and tested computationally: (a) the polychronization method, which explores the differences between the phases of the generator oscillators; (b) the amplitude method which detects clusters of the generator and (c) the resonance method which is based on the resonances between the frequencies of the generator modes and the internal frequencies of the damped oscillations of the reader cells. Pro and contra of these methods have been analysed.

Chemical memory with states coded in light controlled oscillations of interacting Belousov–Zhabotinsky droplets

The information storing potential of droplets, in which an oscillatory, photosensitive Belousov–Zhabotinsky (BZ) reaction proceeds, is investigated experimentally.

Cancer classification with a network of chemical oscillators

We discuss chemical information processing considering dataset classifiers formed with a network of interacting droplets.

On Emulation of Flueric Devices in Excitable Chemical Medium

Flueric devices are fluidic devices without moving parts. Fluidic devices use fluid as a medium for information transfer and computation. A Belousov-Zhabotinsky (BZ) medium is a thin-layer spatially extended excitable chemical medium which exhibits travelling excitation wave-fronts. The excitation wave-fronts transfer information. Flueric devices compute via jets interaction. BZ devices compute via excitation wave-fronts interaction. In numerical model of BZ medium we show that functions of key flueric devices are implemented in the excitable chemical system: signal generator, and, xor, not and nor Boolean gates, delay elements, diodes and sensors. Flueric devices have been widely used in industry since late 1960s and are still employed in automotive and aircraft technologies. Implementation of analog of the flueric devices in the excitable chemical systems opens doors to further applications of excitation wave-based unconventional computing in soft robotics, embedded organic electronics and living technologies.

Microfluidic platform for reproducible self-assembly of chemically communicating droplet networks with predesigned number and type of the communicating compartments

We report a microfluidic system for individually tailored generation and incubation of core-shell liquid structures with multiple cores that chemically communicate with each other via lipid membranes. We encapsulate an oscillating reaction-diffusion Belousov-Zhabotinsky (BZ) medium inside the aqueous droplets and study the propagation of chemical wave-fronts through the membranes. We further encapsulate the sets of interconnected BZ-droplets inside oil-lipid shells in order to i) chemically isolate the structures and ii) confine them via tunable capillary forces which leads to self-assembly of predesigned topologies. We observe that doublets (pairs) of droplets encapsulated in the shell exhibit oscillation patterns that evolve in time. We collect statistical data from tens of doublets all created under precisely controlled, almost identical conditions from which we conclude that the different types of transitions between the patterns depend on the relative volumes of the droplets within a chemically coupled pair. With this we show that the volume of the compartment is an important control parameter in designing chemical networks, a feature previously appreciated only by theory. Our system not only allows for new insights into the dynamics of geometrically complex and interacting chemical systems but is also suitable for generating autonomous chemically interconnected microstructures with possible future use, e.g., as smart biosensors or drug-release capsules.

Binary full adder, made of fusion gates, in a subexcitable Belousov-Zhabotinsky system

In an excitable thin-layer Belousov-Zhabotinsky (BZ) medium a localized perturbation leads to the formation of omnidirectional target or spiral waves of excitation. A subexcitable BZ medium responds to asymmetric local perturbation by producing traveling localized excitation wave-fragments, distant relatives of dissipative solitons. The size and life span of an excitation wave-fragment depend on the illumination level of the medium. Under the right conditions the wave-fragments conserve their shape and velocity vectors for extended time periods. I interpret the wave-fragments as values of Boolean variables. When two or more wave-fragments collide they annihilate or merge into a new wave-fragment. States of the logic variables, represented by the wave-fragments, are changed in the result of the collision between the wave-fragments. Thus, a logical gate is implemented. Several theoretical designs and experimental laboratory implementations of Boolean logic gates have been proposed in the past but little has been done cascading the gates into binary arithmetical circuits. I propose a unique design of a binary one-bit full adder based on a fusion gate. A fusion gate is a two-input three-output logical device which calculates the conjunction of the input variables and the conjunction of one input variable with the negation of another input variable. The gate is made of three channels: two channels cross each other at an angle, a third channel starts at the junction. The channels contain a BZ medium. When two excitation wave-fragments, traveling towards each other along input channels, collide at the junction they merge into a single wave-front traveling along the third channel. If there is just one wave-front in the input channel, the front continues its propagation undisturbed. I make a one-bit full adder by cascading two fusion gates. I show how to cascade the adder blocks into a many-bit full adder. I evaluate the feasibility of my designs by simulating the evolution of excitation in the gates and adders using the numerical integration of Oregonator equations.

Understanding Networks of Computing Chemical Droplet Neurons Based on Information Flow

In this paper, we present general methods that can be used to explore the information processing potential of a medium composed of oscillating (self-exciting) droplets. Networks of Belousov–Zhabotinsky (BZ) droplets seem especially interesting as chemical reaction-diffusion computers because their time evolution is qualitatively similar to neural network activity. Moreover, such networks can be self-generated in microfluidic reactors. However, it is hard to track and to understand the function performed by a medium composed of droplets due to its complex dynamics. Corresponding to recurrent neural networks, the flow of excitations in a network of droplets is not limited to a single direction and spreads throughout the whole medium. In this work, we analyze the operation performed by droplet systems by monitoring the information flow. This is achieved by measuring mutual information and time delayed mutual information of the discretized time evolution of individual droplets. To link the model with reality, we use experimental results to estimate the parameters of droplet interactions. We exemplarily investigate an evolutionary generated droplet structure that operates as a NOR gate. The presented methods can be applied to networks composed of at least hundreds of droplets.

Chemical computing with reaction-diffusion processes

Chemical reactions are responsible for information processing in living organisms. It is believed that the basic features of biological computing activity are reflected by a reaction-diffusion medium. We illustrate the ideas of chemical information processing considering the Belousov-Zhabotinsky (BZ) reaction and its photosensitive variant. The computational universality of information processing is demonstrated. For different methods of information coding constructions of the simplest signal processing devices are described. The function performed by a particular device is determined by the geometrical structure of oscillatory (or of excitable) and non-excitable regions of the medium. In a living organism, the brain is created as a self-grown structure of interacting nonlinear elements and reaches its functionality as the result of learning. We discuss whether such a strategy can be adopted for generation of chemical information processing devices. Recent studies have shown that lipid-covered droplets containing solution of reagents of BZ reaction can be transported by a flowing oil. Therefore, structures of droplets can be spontaneously formed at specific non-equilibrium conditions, for example forced by flows in a microfluidic reactor. We describe how to introduce information to a droplet structure, track the information flow inside it and optimize medium evolution to achieve the maximum reliability. Applications of droplet structures for classification tasks are discussed.

Chemical computing based on Turing patterns in two coupled cells with equal transport coefficients

Two diffusively coupled reaction cells with a nonlinear reaction are used to perform chemical computing based on targeted perturbations switching between two Turing patterns defining two states of a logical device.