Abstract
The exponential growth of the power of modern digital computers is based upon the miniaturization of vast nanoscale arrays of electronic switches, but this will be eventually constrained by fabrication limits and power dissipation. Chemical processes have the potential to scale beyond these limits by performing computations through chemical reactions, yet the lack of well-defined programmability limits their scalability and performance. Here, we present a hybrid digitally programmable chemical array as a probabilistic computational machine that uses chemical oscillators using Belousov-Zhabotinsky reaction partitioned in interconnected cells as a computational substrate. This hybrid architecture performs efficient computation by distributing information between chemical and digital domains together with inbuilt error correction logic. The efficiency is gained by combining digital logic with probabilistic chemical logic based on nearest neighbour interactions and hysteresis effects. We demonstrated the computational capabilities of our hybrid processor by implementing one- and two-dimensional Chemical Cellular Automata demonstrating emergent dynamics of life-like entities called Chemits. Additionally, we demonstrate hybrid probabilistic logic as a viable logic for solving combinatorial optimization problems.
Similar content being viewed by others
Introduction
The exponential increase in computing power has been driven by the vast growth of transistors on silicon chips1. This growth was made possible by developments in fabrication technology reducing the feature sizes of the transistors2, but this paradigm is currently approaching the limits imposed by physics as quantum effects become more pronounced3. To overcome these limitations, various novel computational architectures based on physical and chemical processes have been proposed4,5,6,7,8. Quantum computers show the potential to solve problems that are intractable on classic computing machines but currently suffer from scalability issues due to error correction9. Alternative computing substrates and unconventional computation paradigms are being developed based on mapping computational logic to various physical phenomena10. These tend to emulate transistor-based logic gates and other circuit components into the physical domain using architectures based on Boolean circuits11 or discovered using artificial intelligence12. Other classic computational architectures which utilize the true nature of physical phenomena include reaction-diffusion13 and neuromorphic computers14. These architectures and their algorithms have been designed to solve a specific set of abstract mathematical problems15,16,17,18. The challenge is to develop new platforms that can take advantage of chemical substrates but can be easily programmable.
The key concept of hybrid electronic–chemical computation is inspired by the previously developed concepts of heterotic19 as well as physical computation20 where an amalgamation of different computational systems such as digital, chemical, optical, etc., are used together for efficient computation. Herein, we present a probabilistic information processor based on a hybrid electronic–chemical computational architecture that is digitally addressable which utilizes both analogue and their digital-equivalent states for information processing. The efficiency is gained by distributing the computational task into different computational substrates which are efficient in performing specific operations using their natural propensity. This device is built from a programmable chemical array-based21 around a new type of hybrid chemical state machine for information processing. The hybrid computational array operates using the Belousov–Zhabotinsky (BZ) reaction, which is an oscillating chemical reaction22,23. The cells are arranged in a rectangular array and are individually programmable with a central oscillator driver and have four tuneable gates that allow the oscillations to be coupled between adjacent cells in the grid, which defines the basic computational element. By programming cell and interfacial controls, the emerging chemical oscillations of the ith cell (together with neighbouring cells defined with index j) are defined by the analogue chemical state (\(\underline{{\rm {C{S}}}_{i}^{t}}\)) that can be monitored and mapped to a digitally representable chemical state (\({{\rm {C{S}}}_{i}^{t}}\)) for each cell at time step t. The chemical state evolution in the analogue domain occurs via physical processes such as chemical oscillations, while the digital state is based on a finite state logic. Both the analogue and digital representations of chemical states are synchronized using a clock signal created by the chemical oscillations. Thus, this approach uses a hybrid electronic–chemical logic where the digital domain applies deterministic logic using finite state machine implemented on digital states (\({{\rm {C{S}}}_{i}^{t}}\)) and the chemical domain performs analogue computation using interacting chemical oscillations. A simplified representation of the time evolution of the hybrid electronic–chemical logic in the analogue domain using chemical states (\(\underline{{\rm {C{S}}}_{i}^{t}}\)) can be defined by
where F represents the digital operation based on the readout of the chemical states and G represents interactions in the physical system. F can be defined as a digital finite state machine (FSM) that reads the current analogue signal from the interacting physical system between nearest neighbours and performs operations back into the analogue domain and \(\underline{{\rm {C{S}}}_{j}^{t}}\) represents chemical states of neighbouring cells of the primary cell with chemical state \(\underline{{\rm {C{S}}}_{i}^{t}}\). The physical system evolves in a high dimensional space controlled by the complex myriad interactions between the chemical oscillations and their hysteresis effects together with the effect of localized digital operations and leads to the probabilistic outcomes of the emerging chemical states. Hence, the computation is performed by iterating single-step operations involving analogue and digital states with finite state logic defined for a specific problem.
In our experiments, using a digitally programmable input–output (I/O) system, we showed that it is possible to have an error correction system whereby the oscillatory chemical state can be reinforced, and the effects of phase shift can be eliminated so that the hybrid electronic control system amplifies and stabilizes the chemical states (see Fig. 1). As a proof of programmability and computation, we showed that the system could embody a hybrid electronic–chemical form of cellular automaton (CA) and increase the connectivity in the configuration space way beyond the digital feature space even for simple rules of the elementary CA24. The probabilistic nature of the hybrid computational architecture was demonstrated by implementing a two-dimensional probabilistic chemical cellular automata (CCA) which shows emergent dynamics similar to that seen in Conway’s Game of Life25. Additionally, we also implemented solutions to combinatorial optimization problems such as number partitioning26, Boolean satisfiability27, and the travelling salesman problem28 using a hybrid computational approach.
Results
The chemical computing platform
The chemical array exploits the excitability of the non-linear chemical oscillations of the Belousov –Zhabotinsky (BZ) reaction via localized spatial control22,23. The BZ reaction is highly excitable, can be maintained far from equilibrium, and it can be both spatially and temporally addressable. To exploit these features, we designed an experimental setup where we can programme the addressable medium using an electronically controlled input system (see Fig. 1a, b). The experimental architecture consists of a 3D-printed 1D and 2D grid of interconnected reactors supported on an array of motors equipped with magnetic heads. At the centre of each reactor and the interface of neighbouring cells, magnetic stirrers are placed to match the position of the motor shaft. Each motor is individually addressable, and its speed is controlled using a pulse width modulation (PWM) signal generated using a microcontroller. The schematic diagram and the physical implementation of the two-dimensional experimental setup are shown in Fig. 2, Supplementary Information Section 1 and Supplementary Video 1. At the start of each experiment, the reagents required to initiate the BZ reaction (solutions of malonic acid, potassium bromate, an iron-based redox catalyst and sulfuric acid) are added to the reactors using an automated liquid handling system interface (see the “Methods” section and Supplementary Information Section 1). The role of the central cell stirrer is to initiate and then maintain the chemical oscillations and their amplitude in a programmable way by varying the stirring speed. The mass transfer due to the hydrodynamic coupling between the neighbouring cells leads to interactions between the chemical oscillations of the BZ reactions (see Fig. 1c, d) and Supplementary Video 2. By tuning the speed of the interfacial stirrers (PWM levels), we can control and programme the strength of intercellular couplings between two interacting cells and limit them to their nearest neighbours.
The BZ oscillations induced in a single cell are extremely sensitive to the composition of local redox species and the time of the actuation of the stirrer. This extreme sensitivity to the initial state, localized fluctuations, and time of actuation causes the phase of the chemical oscillations in individual cells to show significant drift with time. We observed that these unfavourable phase shifts between individual cells potentially limited the programmability of the system, and an error correction process to prevent decoherence was needed. To develop a hybrid computational logic, we investigated the phenomenological behaviour of oscillations in single and coupled neighbouring cells. The BZ system combined with the stirrer acts like a forced-damped oscillator. Chemical oscillations appear when the cell stirrer is active (forced oscillator) and once the stirrer is deactivated, the amplitude of oscillations starts falling (damped oscillator) and disappears completely after a few cycles. In the case of coupled neighbouring cells, no interaction occurs between them when the interfacial stirrer is off. On activating the interfacial stirrer, bidirectional interaction occurs and weakly coupled neighbouring cell oscillations come in phase independent of their initial phase differences. Additionally, the interactions in both one- and two-dimensional geometry are confined to nearest neighbours only (see Supplementary Information Section 2). In the digital domain, for the finite state logic, two chemical states were defined as \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=0\) or \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=1\). To eradicate the errors in defining digital states due to the strong hysteresis effects in the chemical domain, the chemical state \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=0\) was defined by low-amplitude chemical oscillations instead of no oscillations. The low oscillation modes were generated by stirrers in pulsing mode. The high amplitude chemical oscillation was represented in the digital domain as \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=1\) and was created in the chemical domain by continuous mode at a high stirring rate. Additionally, to address the potential for errors resulting from state-decoherence (phase shift) in the oscillations among the cells, we found that the introduction of a global ‘clock’ signal (SYNC) was possible by creating weak coupling between neighbouring cells via interfacial stirrers. The weak coupling keeps chemical oscillations over all the cells synchronized (see Supplementary Information Section 2). These global weak oscillations are used for defining clocking signals which allows us to prevent unwanted dephasing (see Fig. 3).
As the BZ reaction proceeds and the malonic acid ‘fuel’ is consumed, the amplitudes of the oscillations decrease with time (see Fig. 3a, b). The deviation between the oscillations in the interconnected network of cells is shown in Fig. 3c. To consistently define discrete chemical states based on the observed colour amplitudes, we trained a convolutional neural network (CNN) on a dataset of time-dependent images labelled with discrete states. The three discrete states based on time-dependent colour classification are red (R), light blue (LB) and blue (B) which are also referred to as CNN states. In the presence of only weak oscillations, two distinct CNN states (Red and Light Blue) are recognized, while in the presence of strong oscillations, three distinct states (Red, Light Blue and Blue) emerge (see Fig. 3d), and we use them to describe the pattern of oscillations in a finite state machine called as rFSM (recognition finite state machine) to determine the chemical states. An interpreted simplified example of clock states corresponding to Fig. 3d is shown in Fig. 3e (see Supplementary Information Section 4 for actual implementation). In our current design, the rFSM logic consists of two different digital states: \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=0\) and \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=1\). The emergence of state \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=0\) occurs when a weak wave pattern is observed (R → LB → R) and the emergence of state \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=1\) occurs when a stronger oscillation is observed (R → LB → B → LB → R). From the chemical state readout, a feedback loop can be completed by implementing a deterministic digital state machine (D) that reads out the chemical states of all or a subset of cells as inputs and returns the PWM levels of the stirrers. These new PWM levels were then applied to the cell and interfacial stirrers (see Supplementary Information Sections 3 and 4 for the implementation of rFSM).
Probabilistic logic and one-dimensional chemical cellular automata (1D-CCA)
A single information processing loop can be represented by a combination of four state machines (C,T,D,P) acting on digital (\({{{\rm {CS}}}}_{i}^{t}\)) and an analogue representation of chemical states (\(\underline{{{{{{{\rm{CS}}}}}}}_{i}^{t}}\)) (see Fig. 4a and b). The state machine C is probabilistic and represents the evolution of \(\underline{{{{{{{\rm{CS}}}}}}}_{i}^{t}}\) while D represents a deterministic digital state machine that reads \({{{\rm {CS}}}}_{i}^{t}\) and updates stirrer states (DS). The state machines T and P represent analogue to digital chemical state conversion (\({{{{{{\rm{CS}}}}}}}_{i}^{t}=T(\underline{{{{{{{\rm{CS}}}}}}}_{i}^{t}})\)) and the physical effects of stirrers represent digital to analogue chemical states, respectively. The time evolution of the digital and analogue representation of chemical states in hybrid probabilistic computation can be represented by
where K represents a hybrid state machine comprising of four state machines C, T, D, P and i, j represents the central and neighbouring cells and the state machine T is included in C for simplicity. IC defines the initial conditions. The emergence of the new chemical state has both implicit and explicit dependence on the previous chemical states. The implicit dependence comes from the hysteresis effect of the oscillations and the explicit dependence comes from the physical interaction of the stirrer state (P) on the oscillations which depends on the previous chemical states via finite state logic \(D({{{{{{\rm{CS}}}}}}}_{i}^{t},\, {{{{{{\rm{CS}}}}}}}_{j}^{t})\).
To demonstrate the working principle of the closed-loop hybrid electronic–chemical logic and the programmability with clocking logic, we implemented the elementary cellular automata (CA) rules (Rule 30, 110 and 250, see Supplementary Information Section 4) in a fully deterministic way, see Rule 30 as an example Fig. 4c and Supplementary Video 3. In the deterministic mode as there is a one-to-one mapping between stirrer (DS) and chemical states (\({{{\rm {CS}}}}_{i}^{t}\)), information loops through digital and chemical domains where the chemical state machine mirrors the digital state machine precisely (see Fig. 4a). Next, we introduced the probabilistic computational logic by introducing hybrid automaton rules (see Fig. 4b). Consequently, the new chemical states in the analogue domain emerge probabilistically where CA rule 30 is modified by creating asymmetric actuation on interfacial stirrers as an example (see Fig. 4c) (and Supplementary Information Section 4). The probabilistic outcomes of the chemical states can be seen in a single cell with the effect of stirrer speed on the emerging peak amplitudes (see Fig. 4d). By introducing new automaton rules that exploit the probabilistic computation mode enabled by the high dimensional space associated with the chemical state machine (C), larger configuration spaces can be explored.
Inspired by the elementary CA rules, we also developed a family of one-dimensional chemical cellular automata (1D-CCA) rules based on a closed-loop state machine utilizing PWM and chemical states. A CCA rule can be defined on an experimental platform based on the observed chemical states and selected PWM levels of interfacial and cell stirrers. In the absence of interfacial stirrers, we can recreate elementary CA rules within a closed feedback loop due to one-to-one mapping between the PWM states of cell stirrers and observed chemical states. However, when the interfacial stirrers are active, due to hydrodynamic coupling between the neighbouring cells and coupled hysteresis effects, one-to-one mapping between PWM and chemical states does not exist and novel patterns can emerge out. As an example, we describe the 1D-CCA state machines combined with a simple phenomenological model which was then used to simulate the emergence of patterns in a one-dimensional geometry (see Supplementary Information Section 5). 1D-CCA rules which define the state machine to act on stirrer based on chemical states comprised of two different values, {Crule}–{Irule}. Crule updates the central cell stirrer PWM state based on the chemical states of nearest-neighbouring cells \(({C}_{i}^{t},\, {C}_{i-1}^{t},\, {C}_{i+1}^{t})\) similar to the elementary CA rule table. Irule updates the PWM states of the two interfacial stirrers based on the chemical states of the two connecting cells \(({C}_{i-1}^{t},\, {C}_{i}^{t})\) and \(({C}_{i}^{t},\, {C}_{i+1}^{t})\). See Supplementary Information Section 5 for quantification of input and chemical states, probabilistic model, and examples of 1D-CCA using a phenomenological model.
Two-dimensional chemical cellular automata (2D-CCA)
By extending the CCA into two dimensions it is possible to explore the dynamics and emergence of complex patterns based on local rules defined by the hybrid electronic–chemical logic. Inspired by Conway’s Game of Life and using a similar concept described in 1D-CCA (based on chemical states and PWM levels of interfacial and cell stirrers), as a discrete model of computation, we developed a two-dimensional probabilistic cellular automaton with extended multi-cellular life-like entities. In the 2D-CCA, we defined the basic emergent life-like units as chemical entities (Chemits), which are comprised of a combination of five nearest-neighbouring cells in the von Neumann neighbourhood (see Fig. 5a, b) for experimental and pictorial representation of the Chemit. The objective is to observe emergent dynamics of Chemits such as replication-propagation in the presence of weak fluctuating oscillations over the whole cell array. The positions and dynamics of Chemits are defined by the combination of digital and analogue representation of chemical states as described previously. We implemented a 2D-CCA digital state machine that takes the chemical states (\({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=\{{{{{\mathrm{0,1}}}}}\}\)) emerging from the probabilistic outcomes of the chemical information processing and outputs 4-state PWM logic with states {\({{S}}_{{0}},\, {{S}}_{{1}},\, {{S}}_{{2}},\, {{S}}_{{3}}\)} which defines the stirrer operations. S0 corresponds to the inactive stirrer state, S1 introduces the random fluctuations (weak chemical oscillations on randomly selected cells), S2 creates the Chemit core and S3 introduces the extended interacting body of the Chemit which leads to the interactions with the surroundings. The chemical oscillations created by PWM state S1 in the absence of Chemits only leads to the \({{{{{{\rm{CS}}}}}}}_{{i}}^{{t}}=0\) chemical state and hence cannot create a new Chemit. The Chemit core defined by a high PWM value (S2) creates high \({{\rm {CS}}}_{{i}}^{{t}}=1\) states. The pictorial representation of the complete closed-loop probabilistic logic of two-dimensional CCA representing the emergence of Chemits is shown in Fig. 5c.
In the experiments, based on the local rules the Chemits show emergent dynamics such as propagation, replication, and competition analogous to life-like species (see Supplementary Video 4). The emergence of a high chemical state at a specific cell occurs probabilistically due to the interaction of Chemit cells with local neighbours in the analogue chemical domain and hysteresis effects. These high chemical states at nearest neighbours or next-nearest neighbours in different configurations lead to propagation, replication, and competition events (see Fig. 5b). Propagation and replication events occur when a high chemical state occurs at nearest and next-nearest neighbours. When two Chemits interact with each other, a competition event occurs, and for one Chemit, has a 50% survival chance. However, in some cases, multiple strong oscillations causing high chemical states to occur at nearest and next-nearest neighbours, which leads to random selection among multiple events, where a propagation, a replication, or a competition event among all the neighbours is selected randomly. Figure 5a shows a series of snapshots from the experiments showing the propagation and replication dynamics of Chemit in a 7 × 7 two-dimensional array with periodic boundary conditions. The complete description of the 2D-CCA and pseudo-code is described in Supplementary Information Section 6. Over a 7 × 7 experimental array together with periodic boundary conditions, we observed the emergence of peaks in the population of the Chemits due to sudden localized replication events which later fall due to competition events within the constrained space or resource (see Fig. 5d). The total number of digitally representable high chemical states at a given step leading to Chemits dynamics are shown in Fig. 5e.
The population and the propagation dynamics of the Chemits with given control parameters are governed by the initial number of Chemits, available spatial resources, and random fluctuations and it is interesting that the system shows highly complex emergent behaviour as observed in Conway’s Game of Life25, but instantiated in a physical device in a hybrid manner. To investigate the emergent behaviour over a larger spatial scale, we developed a simple chemical probabilistic state machine based on the observed phenomenological model to simulate the dynamics of Chemits behaviour up to 150 × 150 cell array (see Supplementary Information Section 6 and Supplementary Video 5). Similar to the experimental observations, the simulations show the sudden formation of local population clusters due to fast replication and as well as annihilation due to competition events at local clusters of large populations. We observed unstable population dynamics of Chemits on an array with a smaller number of cells, however, with an increase in the number of cells the Chemits population stabilizes at different levels which depends on the available spatial resource (see Fig. 5f). We further investigated the population dynamics by varying the initial population of Chemits over a 100 × 100 array and observed convergence in the population at the steady state independent of the initial population (see Fig. 5g) and Supplementary Information Section 6 for further characterization. These simulations demonstrate a strong correlation between the global population dynamics and the local probabilistic rules emerging from digital and chemical state machines. The emergent population dynamics and steady-state kinetics of Chemits show similarities to population behaviours observed in evolutionary biology and can be further extended towards computation operations for information processing.
Solving combinatorial optimization problems using hybrid computational machine
Building on the dynamic feedback loop between electronic and chemical states, we implemented a hybrid electronic–chemical information processing algorithm to solve quadratic combinatorial optimization problems. Here, the idea is to utilize the probabilistic logic from the chemical analogue medium in the hybrid processor to reach the problem solution more efficiently than that of purely deterministic logic on an electronic processor. In the context of quantum adiabatic optimization, various combinatorial optimization problems such as number partitioning, satisfiability (SAT) and Hamilton cycles can be formulated as energy/cost minimization problems on an Ising lattice29,30. Inspired by the Ising or equivalent quadratic unconstrained binary optimization (QUBO) formulations of these problems, we implemented a hybrid electronic–chemical state machine capable of performing energy minimization using chemical states or PWM states equivalent to Ising spin variables. The generalized Hamiltonian up to a quadratic coupling is given by
where Qi could represent the chemical state of the cell, \({h}^{(0)}\) is an offset energy term, \({h}_{i}^{(1)}\) defines the self-interaction term of the spin equivalent state \({Q}_{i}\) and \({h}_{{ij}}^{(2)}\) defines the coupling between states \({Q}_{i}\) and \({Q}_{j}\). The Ising formulation of the optimization problem can be represented by a connected graph with self-interactions and pairwise couplings for Hamiltonian formulation and mapping to the chemical array (see Supplementary Information Section 7). The sign of the coupling coefficients describes the positive (ferromagnetic type) and negative (anti-ferromagnetic type) couplings and the magnitude describes the coupling strength. In the ideal case of implementation of a computation algorithm in a hybrid system such as the proposed hybrid processor, the Ising spin equivalent chemical states should flip according to the interactions defined by local couplings. It should let the problem Hamiltonian reach the global minimum energy configuration and the solutions can be interpreted from the corresponding chemical and digital states.
As the simplest example to explore this computational problem, we have implemented an electronic–chemical hybrid state machine where the two chemical states map directly to Ising spins (−1, +1) and demonstrate the “proof-of-principle computation” through the information loop without any neighbouring interactions (see Supplementary Information Section 7 for pseudo codes and implementation details). Based on the emerging chemical states, the configurational energy is estimated from the Hamiltonian in-silico and compared with the lowest energy configuration so far and iterated until the minimal energy configuration is achieved. We show that the hybrid information processor is able to reach the energy minima successfully looping information between analogue and digital domains with higher processing occurring in the digital domain.
To increase the proportion of information processing in the chemical domain, we devised an extended approach where we utilized a combination of the digital state machine coupled with a probabilistic chemical state machine where the PWM states of the cell stirrers map directly to the Ising spin variables. The flow chart of the complete computational scheme is shown in Fig. 6a. For the pairwise neighbouring interactions generating probabilistic outcomes, a lookup table of chemical states as shown in Fig. 6b was created and the problem was mapped on the platform (see Fig. 6c). As a result, the emergence of the new chemical states not only depends on PWM states of stirrers but also the interactions within chemical states defined by the hybrid state machine. At each step, a comparison was made between ideal states from the lookup table and the emerging chemical states and was utilized for the acceptance of the step towards energy minimization. The probabilistic outcome of the new chemical states arises from the combination of analogue and digital processing. This in turn leads to the lowest energy state due to higher connectivity in the configuration space. This improvement in efficiency will become more evident with the scaling and complexity of the problem with large local minimum configurations. By distributing the algorithm between the digital and chemical logic, we demonstrate large-scale combinatorial problems can be solved efficiently31.
To achieve this advantage, we map the variables of the problem to the cells such that all the coupling interactions can be introduced between neighbours, see Fig. 6c which shows the mapping of a fully connected four-number partitioning problem (primary spin cell shown in blue) employing multiple instantiations of the same spins (auxiliary cells shown in red) to accomplish pairwise coupling between all the variables of the Hamiltonian. At each step, after flipping the PWM states of the cell randomly, pairwise operations to estimate the energy change were performed in a parallel approach. In the chemical decision-making step, if the emergence of the new chemical state is consistent with the lookup table, the energy change was recorded as such else was recorded with a negative sign.
The overall energy change was estimated via summation over all the pairwise spins, and the flipping was accepted or rejected similarly to the first hybrid logic scheme according to the principle of energy minimization. Once the energy associated with the current PWM states estimated from the Hamiltonian reaches the global minimum, the solution to the optimization problem is interpreted from the PWM states, see Fig. 6d and e for the solution to the four-number partitioning problem for the number set \(S=\{1,\, 3,\, 4,\, 8\}\) using the two-hybrid algorithms. The spin state configurations of the primary spin cells for the 4-number partitioning problem using the second hybrid algorithm with the corresponding energies are shown in Fig. 6f (see Supplementary Video 6). Other examples including number partitioning, Boolean Satisfiability and Travelling Salesman Problems are demonstrated in Supplementary Information Section 7.
To investigate the influence of probabilistic logic from the chemical state machine in the hybrid approach, we quantified the algorithmic performance by defining hybrid logic as a combination of the digital algorithm coupled with analogue chemical processing. The probabilistic decision-making occurring in the chemical state machine is tuneable by selecting the PWM levels (see Fig. 4d). We formulated an 8-number partitioning problem on the number set \(S=\{1,\, 3,\, 4,\, 9,\, 3,\, 5,\, 3,\, 6\}\) and estimated the success probability with all possible starting (28 = 256) configurations. By calculating the stationary distribution of the possible configurations, we could estimate the probability of finding the configuration with global minimum energy at different values of the deterministic index (see Supplementary Information Section 7). We observe by the reduction in the deterministic index taking advantage of the probabilistic chemical state machine; the success probability distribution shrinks leading to higher chances of finding the solution independent of the initial configuration (see Fig. 6g and h). The distribution shows for a given 8-number partition problem, the hybrid approach finds a solution with all starting configurations however even though the deterministic algorithm shows a high probability of success, many configurations get trapped in local minima and give no solution.
Discussion
Various unconventional computational approaches such as DNA computation, reaction-diffusion computation, optical computation, and reservoir computation utilize features of the physical systems towards computation which can have different computational principles. DNA computation performs computation using DNA sequences and requires a problem-specific design in general. Reaction-diffusion computers in the absence of a digital substrate are extremely difficult to programme and instantiating generalized problems cannot be achieved easily. Similarly, Reservoir computation in the absence of well-defined architecture could lead to strong variation in local rules within the physical substrate and therefore cannot be generalized for any computation. On the other hand, optical and quantum computers are extremely fast and show a promising future in computation, however, currently, they suffer from various scalability problems such as quantum error correction, are extremely expensive, bulky in size and manufacturing is complex. As an intermediate solution to various problems, we demonstrated a hybrid approach toward information processing and computation where an electronically programmable chemical medium is used to process information with a tuneable probabilistic logic utilizing the natural behaviour of physicochemical processes (see Supplementary Information Section 7, Qualification of hybrid electronic–chemical computation). This hybrid approach with the distribution of information between the digital and analogue chemical state machines is generic and implementable on any digitally programmable physicochemical process. Additionally, this could be further miniaturized to utilize stronger probabilistic effects at the molecular scale together with digital programmability and instantiation. The modern computation algorithms such as machine learning (such as deep neural networks) require non-linear behaviour which is inherently present in chemical systems. This makes hybrid-computation architecture highly resource-efficient for a specific set of problems where non-linearities and probabilistic behaviour are crucial. A highly scalable architecture using hybrid computation principles as described in the manuscript could be utilized for deep learning applications (see Supplementary Information Section 8 for further discussion). The proposed architecture using BZ reaction as a programmable analogue system exhibits hybrid cellular automata rules using both digital and analogue states. We also discovered large-scale emergent dynamics of life-like entities (Chemits) resulting from the local probabilistic interactions between digital and chemical state machines. As a further extension, we also demonstrated hybrid approaches toward solving combinatorial optimization problems by minimizing Ising Hamiltonian formulation utilizing chemical states and nearest-neighbouring couplings. Importantly, we showed the distribution of information processing between the digital and the chemical domains, demonstrating that chemistry is taking an active part in the computation. Using the hybrid computational approach, larger and more complex computational problems can be mapped on the substrate. When the connectivity between the variable is larger than the possible nearest neighbour connections, multiple instantiations of the same spin can also be created to solve the problem. The use of digital components to input information into the chemical domain can be made more scalable by employing optical inputs using spatial light modulators (SLMs) or via electrochemical actuation on high-density CMOS electrode arrays. These approaches are highly scalable and fully programmable using state-of-the-art digital electronics. Similarly, the readout protocols can be substituted by a high-density CMOS sensor array to read the open circuit potential (OCP) and weak continuous flow system to stabilize cell oscillations over a long period. In that case, the use of CNN can be completely avoided reducing the role of digital computation. These input and output strategies can make the computational substrate energetically computationally efficient. The computational power in the hybrid system can be further enhanced by increasing the role of analogue chemical logic while designing new computational algorithms (see Supplementary Sections 7 and 8 for further discussions).
Methods
Experimental platform
The overall hybrid electronic–chemical computational platform consists of three main control domains: (a) chemical domain, (b) experimental setup, and (c) digital domain. (a) Chemical domain which consists of stock solutions required for the BZ reaction that was pumped using syringe pumps sequentially in the right proportion into the mixing chamber. The mixing chamber contains a magnetic stirrer bar that rotates at 140 RPM constantly to ensure the stock solutions are well mixed. Using another pair of syringe pumps, the reaction mixture in the mixing chamber was then transferred to the 3D-printed experimental arena with stirrers in the (b) experimental setup. In this experimental setup, the rotation of stirrers is controlled by DC motors equipped with Neodymium-based permanent magnets located at the bottom of the arena. Each motor’s speed and direction can be individually addressed by the supported electronics control. The BZ chemical oscillations occurring in the experimental arena on the response of stirrer actuation were then observed and recorded by a camera. These temporal oscillatory patterns were then passed into the (c) digital domain where further information processing occurs. The oscillatory patterns were classified into three different states using a convolutional neural network (CNN). There are three different classification states, RED, LIGHT BLUE, and BLUE. These classified states were used for creating a global chemical clock over all the cells as well as the programmable chemical states (CS) for computation. The chemical clocking logic is used over all the experiments as a sync signal for a single feedback loop step. These patterns in the given time frame were then converted into the observed chemical states using a finite state machine (FSM). This FSM reads the temporal CNN states over a given time and returns the digital CA state based on the observed oscillatory behaviour and resets. The hybrid electronic–chemical computational logic is then implemented on these chemical states using various problem-dependent state machines. Once an experiment was finished, the remaining solution was drained into the waste container using a pair of syringe pumps and the experimental arena underwent a series of rinsing and cleaning cycles to get ready for the next experiment.
Preparation of chemical solutions
The stock solutions for the automated platforms were prepared as follows. Ferroin (0.1 M) solution was prepared by dissolving 2.78 g of ferrous sulfate heptahydrate and 5.40 g of 1,10-phenanthroline in 10 mL of deionised water. The solution was then further diluted to 0.001 M for the experiment. Sulfuric acid (1.0 M) solution was prepared by diluting 56 mL of concentrated H2SO4 to 1 L of deionised water. Potassium bromate (0.5 M) KBrO3 solution was prepared by dissolving 83.5 g of KBrO3 in 1 L of 1 M H2SO4. Malonic acid (1.0 M) solution was prepared by dissolving 104 g of CH2(COOH)2 in 1 L of deionised water. Deionised water PURELAB® Option-S/R 7/15 was used as the source of all water used in the experiments and preparation of stock solutions.
State recognition using convolutional neural network
We utilized a supervised learning strategy where are large dataset was created on images of individual cells which were categorised as “red”, “light blue”, and “blue”. Two databases were created for one- and two-dimensional experimental setups. The first database for a one-dimensional platform contains >13,000 images, and the second database from a two-dimensional platform contains >7000 images. We trained a convolutional neural network (CNN) using TensorFlow 1.X (Conv2D) based on the post-processed dataset for the classification of three distinct colour states.
Simulations
All the simulations and further analysis were performed using Wolfram Mathematica and all the source codes are available on the repository. All the simulations regarding the dynamics of Chemits were performed 25 times at each parameter value to obtain statistically significant data.
Data availability
Most of the source data which includes experimental and simulated data used for generating the figures in the main text and supplementary information are available on Zenodo32 https://doi.org/10.5281/zenodo.10732131. Additional data due to large size are available upon request to the corresponding author at Lee.Cronin@glasgow.ac.uk.
Code availability
The code used to operate the platform as various implemented simulation models are available at https://github.com/croningp/BZComputation33 (https://doi.org/10.5281/zenodo.10723443).
Change history
18 March 2024
A Correction to this paper has been published: https://doi.org/10.1038/s41467-024-46859-8
References
Waldrop, M. M. The chips are down for Moore’s law. Nature 530, 144–147 (2016).
Yuan Taur et al. CMOS scaling into the nanometer regime. Proc. IEEE 85, 486–504 (1997).
Hu, C. Future CMOS scaling and reliability. Proc. IEEE 81, 682–689 (1993).
Arute, F. et al. Quantum supremacy using a programmable superconducting processor. Nature 574, 505–510 (2019).
Gorecki, J. et al. Chemical computing with reaction-diffusion processes. Philos. Trans. R. Soc. A. 373, 20140219 (2015).
Rueckes, T. et al. Carbon nanotube-based nonvolatile random access memory for molecular computing. Science 289, 94–97 (2000).
Schrauwen, B., Verstraeten, D. & Campenhout, J. An overview of reservoir computing: theory, applications and implementations. In Proc. 15th European Symposium on Artificial Neural Networks 482 (2007).
Adleman, L. M. Molecular computation of solutions to combinatorial problems. Science 266, 1021–1024 (1994).
Bennett, C. H. & DiVincenzo, D. P. Quantum information and computation. Nature 404, 247–255 (2000).
Fang, Y., Yashin, V. V., Levitan, S. P. & Balazs, A. C. Pattern recognition with “materials that compute”. Sci. Adv. 2, e1601114 (2016).
Katsikis, G., Cybulski, J. S. & Prakash, M. Synchronous universal droplet logic and control. Nat. Phys. 11, 588–596 (2015).
Lin, X. et al. All-optical machine learning using diffractive deep neural networks. Science 361, 1004–1008 (2018).
Adamatzky, A., De Lacy Costello, B. & Asai, T. Reaction-diffusion Computers (Elsevier, 2005).
Torrejon, J. et al. Neuromorphic computing with nanoscale spintronic oscillators. Nature 547, 428–431 (2017).
Kuhnert, L., Agladze, K. I. & Krinsky, V. I. Image processing using light-sensitive chemical waves. Nature 337, 244–247 (1989).
Parrilla-Gutierrez, J. M. et al. A programmable chemical computer with memory and pattern recognition. Nat. Commun. 11, 1442 (2020).
Steinbock, O., Tóth, Á. & Showalter, K. Navigating complex labyrinths: optimal paths from chemical waves. Science 267, 868–871 (1995).
Tsompanas, M.-A., Fullarton, C. & Adamatzky, A. Belousov–Zhabotinsky liquid marbles in robot control. Sens. Actuators B: Chem. 295, 194–203 (2019).
Kendon, V., Sebald, A. & Stepney, S. Heterotic computing: past, present and future. Philos. Trans. R. Soc. A. 373, 20140225 (2015).
Horsman, D., Stepney, S., Wagner, R. C. & Kendon, V. When does a physical system compute? Proc. R. Soc. A. 470, 20140182 (2014).
Tompkins, N. et al. Testing Turing’s theory of morphogenesis in chemical cells. Proc. Natl Acad. Sci. USA 111, 4397–4402 (2014).
Horvath, V., Kutner, D. J., Chavis Iii, J. T. & Epstein, I. R. Pulse-coupled BZ oscillators with unequal coupling strengths. Phys. Chem. Chem. Phys. 17, 4664–4676 (2015).
Petrov, V., Gáspár, V., Masere, J. & Showalter, K. Controlling chaos in the Belousov—Zhabotinsky reaction. Nature 361, 240–243 (1993).
Wolfram, S. Computation theory of cellular automata. Commun. Math. Phys. 96, 15–57 (1984).
Gardner, M. Mathematical games. Sci. Am. 197, 126–141 (1957).
Korf, R. E. A complete anytime algorithm for number partitioning. Artif. Intell. 106, 181–203 (1998).
Hansen, P. & Jaumard, B. Algorithms for the maximum satisfiability problem. Computing 44, 279–303 (1990).
Applegate, D. L., Bixby, R. E., Chvatal, V. & Cook, W. J. The Traveling Salesman Problem: a Computational Study (Princeton University Press, 2011).
Lucas, A. Ising formulations of many NP problems. Front. Phys. 2, 5 (2014).
Guo, S. Y. et al. A molecular computing approach to solving optimization problems via programmable microdroplet arrays. Matter 4, 1107–1124 (2021).
Pierangeli, D., Marcucci, G. & Conti, C. Large-scale photonic ising machine by spatial light modulation. Phys. Rev. Lett. 122, 213902 (2019).
Sharma, A. et al. A Programmable Hybrid Digital Chemical Information Processor based on the Belousov–Zhabotinsky Reaction https://doi.org/10.5281/zenodo.10732131 (Zenodo, 2024).
Sharma, A. et al. A programmable hybrid digital chemical information processor based on the Belousov–Zhabotinsky reaction https://doi.org/10.5281/zenodo.10723443 (Zenodo, 2024).
Acknowledgements
We would like to thank Hessam Mehr and Liam Wilbraham of the University of Glasgow for their discussions. The authors gratefully acknowledge financial support from the EPSRC (Grant Nos. EP/H024107/1, EP/I033459/1, EP/J00135X/1, EP/J015156/1, EP/K021966/1, EP/K023004/1, EP/K038885/1, EP/L015668/1, and EP/L023652/1), the ERC (project 670467 SMART-POM), and the DARPA molecular informatics project.
Author information
Authors and Affiliations
Contributions
L.C. conceived the original idea and together with A.S. designed the project and the research plan. L.C. designed the reactor array, A.S. and M.T.-K.N. designed and built the robotic platform with help from J.M.P.G. J.M.P.G. implemented the computer vision and chemical clocking algorithms, and M.T.-K.N. created the training dataset. M.T.-K.N. and A.S. implemented the 1D-CCA. A.S., Y.J. and M.T.-K.N. implemented the 2D-CCA. A.S. and Y.J. implemented the hybrid computation and M.T.-K.N., Y.J. and A.S. performed the experiments. A.S. did the data analysis, and developed models and ran the simulations with help from Y.J. A.S. helped benchmark the system with L.C. Finally, A.S. and L.C. wrote the paper with help from the rest of the authors.
Corresponding author
Ethics declarations
Competing interests
The work presented here has been filed as a patent PCT/EP2023/057243.
Peer review
Peer review information
Nature Communications thanks Konrad Szacilowski, and the other, anonymous, reviewer for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sharma, A., Ng, M.TK., Parrilla Gutierrez, J.M. et al. A programmable hybrid digital chemical information processor based on the Belousov-Zhabotinsky reaction. Nat Commun 15, 1984 (2024). https://doi.org/10.1038/s41467-024-45896-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-45896-7
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.