Anne-Déborah C.
Nguindjel†
,
Pieter J.
de Visser†
,
Mitch
Winkens†
and
Peter A.
Korevaar
*
Institute for Molecules and Materials, Radboud University, Nijmegen, The Netherlands. E-mail: p.korevaar@science.ru.nl
First published on 29th September 2022
Living organisms employ chemical self-organization to build structures, and inspire new strategies to design synthetic systems that spontaneously take a particular form, via a combination of integrated chemical reactions, assembly pathways and physicochemical processes. However, spatial programmability that is required to direct such self-organization is a challenge to control. Thermodynamic equilibrium typically brings about a homogeneous solution, or equilibrium structures such as supramolecular complexes and crystals. This perspective addresses out-of-equilibrium gradients that can be driven by coupling chemical reaction, diffusion and hydrodynamics, and provide spatial differentiation in the self-organization of molecular, ionic or colloidal building blocks in solution. These physicochemical gradients are required to (1) direct the organization from the starting conditions (e.g. a homogeneous solution), and (2) sustain the organization, to prevent it from decaying towards thermodynamic equilibrium. We highlight four different concepts that can be used as a design principle to establish such self-organization, using chemical reactions as a driving force to sustain the gradient and, ultimately, program the characteristics of the gradient: (1) reaction–diffusion coupling; (2) reaction–convection; (3) the Marangoni effect and (4) diffusiophoresis. Furthermore, we outline their potential as attractive pathways to translate chemical reactions and molecular/colloidal assembly into organization of patterns in solution, (dynamic) self-assembled architectures and collectively moving swarms at the micro-, meso- and macroscale, exemplified by recent demonstrations in the literature.
Non-covalent, reversible associations embedded in the design of the building blocks provide a powerful tool to program assembly processes. Under thermodynamic control, the assemblies will spread out into a homogeneous solution, or form equilibrium structures, such as supramolecular assemblies or crystals (Fig. 1). Developing ‘life-like’ systems however involves a paradigm where self-organization occurs out-of-equilibrium.11–15 For example, the construction of the cells’ cytoskeleton shows how coupling of assembly steps to chemical reactions drives self-organization in time and space.16 First, the tubulin building blocks are activated by a chemical fuel, guanosine triphosphate (GTP), in order to assemble into microtubules. Over time, hydrolysis of GTP occurs, and tubulin becomes inactive again. Hence, the self-organization of the microtubules relies on the supply of fresh GTP, which enables the cell to control the build-up and disassembly of the cytoskeleton. The transient nature of microtubule assemblies has inspired the development of synthetic systems17,18 – varying from 1D nanofibers19–23 to colloidal assemblies24 and macroscale gels25 – that form and fall apart over time, as programmed by the kinetics of the underlying chemical reactions.
This Perspective addresses the challenge to arrive at spatial programmability in the out-of-equilibrium organization of matter that is of interest to systems chemistry, ranging from supramolecular assemblies to colloids, hydrogels, crystalline structures and flow patterns in chemical reaction networks. In general, such organization requires the spontaneous emergence of concentration gradients that provide a spatial differentiation in the system. To avoid decay back to thermodynamic equilibrium (e.g. a homogeneous solution, emulsion, or non-organized dispersion of colloids), these gradients need to be sustained under out-of-equilibrium conditions. Such sustainment requires at first a gradient source; a chemical reaction that takes place in solution, a source element that releases a chemical compound, or a sink element that depletes the compound. Even though patterns can emerge around gradient sources that are a priori localized at a fixed position, in an entirely self-organizing system, one needs a mechanism by which the production, release or depletion of the compound gets localized such, that the gradient of this compound is sustained. In this Perspective, we outline four physicochemical principles that enable design strategies to arrive at sustained gradients in solution-based systems: reaction–diffusion coupling (Fig. 1b), reaction–convection (Fig. 1c), the Marangoni effect (Fig. 1d) and diffusiophoresis (Fig. 1e). We selected these principles as they provide unique pathways to direct the self-organization of systems that either comprise concentration patterns in solution (Fig. 1g), self-assembled structures (Fig. 1h) or dynamic swarms of colloids or (micro)droplets (Fig. 1i), at the micro-, meso- or macroscale.
The first experimental demonstration of a Turing pattern in a synthetic system was provided by Castets et al. in the chlorite-iodide-malonic acid (CIMA) reaction, where starch, used as indicator in the gel, provides a selective complexation and thereby generates the required distinction in diffusion rates between the activator (iodide) and the inhibitor (chlorite).36–39 Furthermore, Turing patterns have been reported in the ferrocyanide–iodate–sulphite reaction40–42 and, under specific conditions, in the Belousov–Zhabotinsky (BZ) reaction.43–45 In this reaction that naturally displays reaction–diffusion patterns, oscillation results from the autocatalytic oxidation of a redox catalyst by bromate ions (activator), a process that is inhibited by bromide ions. A reducing agent (e.g. malonic acid) restarts the cycle, by reducing the catalyst again. Classically, the BZ reaction provides waves at the cm-scale that spiral through the reaction medium in an unstirred batch configuration. However, Epstein et al. showed how the classical BZ reaction can be further tailored in water-in-oil micro-emulsions.46,47 Here, differentiation in diffusivity is obtained as polar species are transferred via the aqueous micro-droplets, whereas non-polar intermediates are transferred more rapidly as they diffuse through the oil-based continuous phase, resulting into a large diversity of spatial organizations, including Turing patterns in the 100 micrometer-range.48
Even though most of these pattern-forming reactions are found by chance,49 it is clear that the emergence of reaction–diffusion patterns relies on the kinetic reaction rate constants and diffusion constants involved,41,50,51 and various models have been developed to simulate reaction–diffusion patterns.11,52–55 Furthermore, generic design principles have been proposed, such as a stepwise method through which an oscillatory reaction occurring in a homogeneous, well-mixed solution can be translated into reaction–diffusion patterns in open one-side-fed reactors.56,57 To this end, a complexing agent is introduced which selectively binds the activator to decrease its diffusion rate (a strategy that was suggested earlier by Lengyel and Epstein38), as well as an antagonist reaction to control its accumulation. All in all, examples of complexity emerging from reaction–diffusion systems have been demonstrated which range from nano- and microscale structures to patterns that span centimeters, and from “static” Turing patterns to dynamic patterns that continuously change over time.49,58–63 Here, we highlight some of the more recent strategies that have been developed to couple reaction and diffusion and generate self-organizing structures and patterns.
Fig. 3 Complex static structures emerge from reaction–diffusion. (a) DMAc functions as the inhibitor in a Turing-like reaction–diffusion system, resulting in patterned PBI-Cu membranes (b). Reprinted (adapted) with permission from ref. 64. Copyright 2021 American Chemical Society. (c) Various polycrystalline biomorphic structures grow from reaction–diffusion triggered BaCO3–SiO2 co-precipitation. Due to the narrow pH range in which passivating SiO2 precipitation takes place, crystal formation is affected by the acidic growth front. This causes the structures to ‘blossom’ towards the bulk solution if it has high pH (top) and ‘curl’ inwards if bulk pH is lower (bottom). (d) BaCO3–SiO2 co-precipitate structure that is formed in the blossoming regime. (e) Ion-substituted structures, exemplifying the shapes obtained in either the blossoming (top row) or curling (bottom row) regime. Reproduced from ref. 72 with permission from the Wiley-VCH publishing group, Copyright 2020. |
Furthermore, inorganic structures with biomorphic shapes such as corals, vases and helices can be generated upon coupling the precipitation reactions of barium carbonate (BaCO3) nanocrystals and amorphous silica (SiO2). García-Ruiz showed how polycrystalline biomorphic structures grow when CO2 diffuses into an aqueous solution of sodium metasilicate (Na2SiO3) and barium chloride (BaCl2).65,66 These chemicals participate in the following (overall) precipitation reactions:
Ba2+ + CO2 + H2O → BaCO3 + 2H+ | (1) |
SiO32− + 2H+ → SiO2 + H2O | (2) |
For example, Yoshida et al. have exploited the BZ reaction to direct biomimetic hydrogels that are able to move without any external driving force.76 An oscillating hydrogel was produced by including a Ru(bpy)3 catalyst – which can participate in the BZ reaction – into a poly(NIPAAm-co-Ru(bpy)3-co-AMPS) polymer network. Upon immersion of the gel in the catalyst-free BZ reaction medium, the BZ reaction starts and the catalyst inside the hydrogel oscillates between Ru(III) and Ru(II) oxidation states. These oxidation states, in turn, influence the hydrophilicity of the gel, and thereby drive its periodical swelling (Ru(III)) and deswelling (Ru(II)), which results into peristaltic motion of the gel (Fig. 4).77 Alternatively, walking motion can be achieved using gels that swell asymmetrically due to an internal gradient in Ru(bpy)3 and AMPS concentration, and self-oscillating polymer brushes – polymer chains grafted onto a glass substrate – can be used to transport nanoscale objects.76,78 Using a similar poly(NiPAAm-co-Ru(bpy)3-co-AMPS) hydrogel immersed in catalyst-free BZ-solution, Lu et al.79 exploit the catalyst's photosensitivity. When weakly illuminated, the frequency of oscillations in the gel increases with light intensity, while above a certain intensity threshold, stronger illumination reduces the oscillation frequency instead. When placed inside a capillary of which the left and right halves are illuminated with different intensities, the authors were able to cause both photophobic and phototropic motion of the gel, as the imbalance between oscillation frequencies in the weakly and brightly illuminated parts of the gel dictated its movement.
Fig. 4 Functional oscillating systems based on the Belousov–Zhabotinsky (BZ) reaction. By incorporating the BZ reaction into a hydrogel, the hydrophilicity change resulting from Ru2+–Ru3+ switching can be exploited to cause periodic swelling and deswelling, giving rise to peristaltic motion. Reproduced from ref. 77 with permission from the Wiley-VCH publishing group, Copyright 2008. |
Additionally, BZ reactions can be used as a basis for chemical communication amongst compartmentalized systems such as micro-droplets or -particles.80–82 For example, Litschel et al. used this concept to create a synthetic ‘central pattern generator’, a motif found in the nervous system that produces pulses which alternatively travel down the left and right sides of its spinal column, inducing cyclic muscle contractions.83 A microfluidic device was created that comprises a linear array of nanoliter wells, connected via microchannels and loaded with the BZ reagent solution. Coupling of the BZ reactions in neighboring wells occurs either by diffusion of the hydrophobic inhibitor Br2 through the PDMS walls (inhibitory coupling), or by water-soluble species, which can propagate oscillations along the microchannels (excitatory coupling). A photosensitive catalyst was included to suppress spontaneous oscillation by initially exposing the wells to light. As a result, BZ oscillations selectively start when the light is turned off at the beginning of the array (the ‘pacemaker’), which excites the subsequent wells, such that traveling waves of BZ oscillations progress down the array. Even if the two parallel arrays are excited simultaneously, their inhibitory coupling upon exchange of Br2 through the PDMS walls causes a phase-shift, such that they end up oscillating out of phase. One can envision such a system being coupled to previously mentioned examples of contracting BZ-gels, creating synthetic analogues of muscles that constantly contract and relax upon actuation by neural signals.
For example, Chen et al. reported a system where rapidly diffusing activator (input) strands displace fluorescent signal strands from reporter hybrids.84 Subsequently, the signal strands can be sequestered by slowly diffusing inhibitor hybrids. Various fluorescent patterns can then be generated from gels with a homogeneous loading of these reporter hybrids. For example, single and concentric rings spread outward from a cavity in the gel that contains both activator and inhibitor strands, while destructive interference patterns appear when one cavity is filled with activator and another with inhibitor, and constructive interference is observed when orthogonal activators are added into separate cavities. Zenk et al. reported a setup that generates sustained DNA concentration gradients and patterns.85 In their design, single and double DNA strands diffuse from opposite liquid reservoirs into a hydrogel, and release an output DNA strand via strand-displacement. Throughout the system, ‘recapturer complexes’ are present at constant concentration in order to sequester the output strands. Together, this constant production and depletion of the output strands maintains a stable, yet out-of-equilibrium concentration profile in the hydrogel. This approach provides the foundation for the experimental realization of (programmed) DNA-based pattern generating systems which thus far have only been demonstrated in silico, such as the pattern-processing modules used by Scalise et al. in their simulations.86
In addition to these DNA strand displacement reactions, Zadorin et al. have demonstrated how positional information can be established in DNA-based systems with chemical reactions that produce and deplete DNA strands.87 In the example highlighted in Fig. 5a, autocatalysis of activator DNA strand A1 is established via template strand TA1. Upon binding of the activator, TA1 templates the extension of A1 into a strand that comprises two sequences of A1 – a reaction that is catalyzed by DNA polymerase. Next, a nicking enzyme cleaves the extended strand such that two A1 strands are released from this autocatalytic system. Finally, A1 is converted into waste when DNA polymerase acts upon it while it is bound to repressor strand R1. By combining autocatalysis and inhibition, a weak gradient in repressor strand can be transferred into a strong, almost stepwise increase in A1 concentration – in analogy to Wolpert's scenario – with the repressor R1 functioning as the morphogen. By coupling two orthogonal autocatalytic motifs via a morphogen that contains both repressor strands, a ‘French flag pattern’ with three chemically distinct zones can be generated based on a single morphogen gradient (Fig. 5b). This gradient is then exploited to selectively bind colloidal beads in a particular area via DNA hybridization.
Fig. 5 DNA-based reaction–diffusion systems. (a) DNA reaction networks for autocatalysis, repression and degradation. (b) A bifunctional repressor strand functions as the morphogen in a demonstration of DNA-based positional information. The repressor gradient is amplified by the autocatalytic reactions, resulting in three zones with distinct chemical compositions, i.e. a ‘French flag pattern’.87 (c) A DNA-based density-dependent autocatalytic reaction is triggered as a network of beads, microtubules and kinesin motors contracts. (d) As a result, a fluorescent front is initiated from the sites where a critical bead density is reached.88 |
Self-organization processes in life, such as embryonic development, often involve a coupling of chemical as well as mechanical systems. Such coupling can pave the way towards new control mechanisms in artificial self-shaping materials. In work by Senoussi et al., DNA reaction networks were implemented in active matter-based gels that, upon contraction, initiate propagating reaction fronts.88 To this end, two different types of microbeads are functionalized at their surface, each with one type of single DNA strand, where bead type B templates the formation of the activator for bead type D, and vice versa (Fig. 5c). These beads are loaded into a solution with an activator strand-degrading exonuclease, as well as microtubules and kinesin motors, which together form an active gel that concentrates the microbeads upon contraction. When the microbeads are brought close together, they easily exchange their activator strands in cross-catalysis. Next, activator species B is coupled to a motif (TBA) that enables autocatalysis of activator A. As a result, an autocatalytic front of A emerges, which is triggered from the location where the active gel has started to contract (Fig. 5d).
Finally, the spatiotemporal patterns that emerge from the Min protein system provide an intriguing platform to study in vitro pattern formation. The Min protein system is a biological reaction–diffusion system in Escherichia coli that comprises the MinD and MinE proteins, and forms patterns that localize the middle plane to assist in bacteria cell division. MinD switches from a solubilized monomeric state in solution to a membrane-bound dimeric state upon binding ATP (Fig. 6a). MinD at the membrane sequesters MinE from the solution, which elevates the ATPase activity of MinD. As result ATP is hydrolyzed to ADP, and MinD dissociates from the membrane again. In the reaction–diffusion patterns that emerge from this system, MinD serves as the activator, and MinE as the inhibitor. To introduce external control, comparable to the photocatalytic BZ reactions, Glock et al. introduced an azobenzene-bound peptide derived from MinE. This peptide binds to MinD when the azobenzene moiety is in the cis state, activating its ATPase activity and thereby triggering its dissociation from the membrane (Fig. 6b).89 Importantly, the binding and unbinding of the peptide-azobenzene derivative to MinD can be controlled by promoting the cis or trans state upon exposure to UV and blue light, respectively. Notably, entrainment of the spatiotemporal patterns of the MinDE system was obtained after periodic illumination by UV, which remained stable for more than 20 illumination cycles (Fig. 6c). When the periodic illumination was stopped, the spatiotemporal patterns returned to their ground state after a few oscillations. The relative simplicity of this network and the manipulation technique offers possibilities for control of nonlinear dynamics of natural reaction–diffusion systems. This opens pathways to establish patterns in synthetic biology contexts, such as the development of a synthetic cell.
Fig. 6 Entrainment of Min-protein patterns based on a photoswitchable MinE peptide. (a) Monomeric MinD in solution forms a dimer that attaches to the membrane upon binding ATP. (b) Upon illumination with UV, the MinE peptide associates with membrane-bound MinD, causing hydrolysis of ATP followed by release of MinD from the membrane. (c) Periodic illumination of the sample synchronizes the patterns with the light pulses. Reproduced from ref. 89 with permission from the Wiley-VCH publishing group, Copyright 2018. |
Solutal buoyancy convection is described by the solutal Rayleigh number, the ratio between the rate of convective transport due to solutal buoyancy and the rate of diffusive transport:
(1) |
Importantly, if Ra ≫ Rac, strong convection causes a rapid equilibration of the system, resulting in a limited lifetime of the convection patterns. Hence, experimental systems with intermediate Ra values are required to establish solutal buoyancy convection patterns.
Experiments in our lab with the pH indicator bromothymol blue showed similar convection patterns (Fig. 8). We have activated a photoacid generator to decrease the pH in the top layer of the solution upon exposure to UV. Starting from pH ≈ 8, the acidification resulted in a change of the pH indicator bromothymol blue from green to yellow. Convection patterns emerged over the course of 5 to 10 minutes, varying in shape from stripes (Fig. 8a) to branches and dots (Fig. 8b). Intriguingly, we observed that the patterns can be varied by varying the depth of the solution, i.e. Lc.
Surprisingly, despite the abundance of chemical reactions that potentially give rise to density differences, the implementation of chemo-hydrodynamic effects in chemical systems as a design principle, rather than a phenomenon as such, is quite rare. However, recent theoretical studies emphasize the untapped potential of reaction–convection: self-generated convection can accelerate chemical reactions,104 and enhance nonlinear kinetics such as chemical oscillations.105
To engineer self-organizing convection patterns in solution, careful design of the chemical reactions is important. Enzymatic reactions with a constant production rate – in order to generate a steady density gradient – can relatively easily be tuned via the enzyme concentration. Furthermore, it is relatively easy to spatially restrict the region of chemical activity by immobilizing the enzyme(s), e.g. at the top or at the bottom of a reaction container. For example, Zhang et al. use an enzymatic reaction network comprised of two enzymes, glucose oxidase (GOx) and horseradish peroxidase (HRP).106 GOx utilizes atmospheric oxygen at the top of the reaction container to oxidize glucose into gluconic acid and hydrogen peroxide. Due to gluconic acid being ∼0.2% more dense than glucose, a vertical density gradient emerges.107 The patterns are visualized via hydrogen peroxide, a side product of the glucose oxidation. HRP uses hydrogen peroxide to oxidize 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) into ABTS˙+. This green product visualizes the solutal buoyancy convection patterns. The patterns that appear were demonstrated to depend strongly on the liquid depth and evolve over time into different configurations.
Apart from self-organization into convection cells of various geometrical shapes, solutal buoyancy convection can be tailored to self-organize flow patterns into so-called chemical pumps.108–111 By immobilizing enzymes at fixed locations in a container, substrates and products can be ‘pumped’ over centimetre-long distances if the density differences between substrates and products are large enough. The pump velocity is governed by the substrate concentrations and the reaction rate. Intricate flow patterns can be designed by combining sequential enzymatic reactions. For example, Maiti et al. engineered a bifurcating flow pattern using three different enzymes: glucose oxidase (GOx), catalase (CAT) and acid phosphatase (AP).110 These three enzymes were separately immobilized in pieces of agarose gel and spaced apart 6–7 mm in a triangular configuration inside a reaction container. GOx and CAT take part in a reaction cascade while AP catalyses a separate reaction:
(2) |
(3) |
Marangoni flows are used in a wide variety of systems, and the examples selected here outline the diversity of approaches to create complex spatial organization. We focus on flows induced by the soluto-capillary effect, as – contrary to thermally-induced flows – they can be controlled by chemical reactions, and thereby offer a greater ability to program Marangoni flow in self-organizing systems.
A supramolecular approach was used by Cheng et al. in their attempt to extend the lifetime of a Marangoni flow induced by a gradient of sodium dodecyl sulfate (SDS) at an air–water interface.134 Similar to the boat experiments presented above, SDS is gradually released in solution from a hydrogel source, giving rise to a Marangoni flow sustained for approximatively 45 s, before the surface tension is evened out. The system is altered to prevent interfacial saturation with the introduction of cyclodextrin as a complexation agent in the medium. By sequestering the surfactant in a supramolecular complex that readily dissolves in the bulk solution, the system is sustained far from equilibrium. The simultaneous release and capture of SDS dramatically extends the time for which the Marangoni flow is sustained by a factor 50.
To further control the way the flow is established and the period for which it is sustained (and so for which ΔC > 0), chemical reactions that impact the surface tension are particularly useful. Here, instead of relying on deposition, a gradient is established from an initially homogeneous interface by introducing reactions that locally produce or consume surface active compounds. The chemistry of fatty acid compounds is mainly used in such examples, because their acid/base properties are directly coupled with their surface activity. Indeed, the difference in surface tension between the protonated (high γ) vs. the deprotonated (low γ) forms of these compounds is significant enough to generate instabilities, meaning that the presence of a pH gradient can eventually lead to Marangoni flow. Lagzi et al.135 and later Suzuno et al.136 take advantage of the properties of 2-hexyldecanoic acid to solve mazes at a liquid-air interface. In both examples, a pH gradient is established in a solution containing the fatty acid by placing a hydrogel pre-soaked in concentrated acid at the exit of the maze. In that setting, protons diffuse away from the hydrogel source and the subsequent protonation of the fatty acid occurs unevenly at the interface, giving rise to a Marangoni flow that drags objects floating at velocities up to 10 mm s−1 directly towards the maze exit within seconds.
In a work recently published by our research group, the need for an external source of protons is eliminated upon using a photochemical reaction.137 Starting from a homogeneous aqueous solution containing a photoacid generator and the surface active sodium oleate, a pH gradient is established upon local irradiation of the interface (Fig. 10a). The rapid acidification of the medium leads to the protonation of oleate into oleic acid. Apart from displaying a higher surface tension than its deprotonated counterpart, oleic acid tends to phase-separate into oil droplets, which enhances the gradient. We have demonstrated that the Marangoni flow, by continuously bringing unreacted photoacid and oleate towards the irradiated area could sustain itself over extended periods of time of up to 30 minutes. Moreover, careful investigation into the ratio of sodium oleate over photoacid introduced in the system demonstrated that the flow velocity and duration could be accurately controlled.
Fig. 10 (a) Self-sustained Marangoni flow at a 2D air–water interface, triggered by local UV exposure. Photoacid generator 1-(2-nitrophenyl)ethyl sulfate (NPE sulfate) produces an acid, which protonates the surfactant oleate. Thereby, the surface tension increases at the site of UV exposure, where phase separated oleic acid droplets accumulate. Adapted from ref. 137 with permission from the Wiley-VCH publishing group, Copyright 2021. (b) Propulsion of octanol droplets in an aqueous solution (i.e. 3D system): an asymmetric occupation of lipid (of which the formation is catalyzed by the octanol-water interface) at the droplet interface results into an internal Marangoni flow in the droplet, which drives their propulsion. Adapted from ref. 147. |
Considerable advantages of chemistry-induced gradients can be anticipated with the way chemical reactions can be coupled with one another into networks. The reversibility of acid/base chemistry can also be an asset in efforts to increase the complexity of the system and possibly generating feedback loops.
Over the past decade, photoswitches have gained traction for their ability to remotely modify surface properties, because of their light sensitivity and their high reversibility.138,139 In a variety of studies, amphiphilic molecules containing a photoswitchable azobenzene moiety have been synthetized to control interfacial mass transfer. These systems, exemplified by the work of Diguet et al., rely on the drastic change in surface activity displayed by the trans (low γ) and cis (high γ) isomers (in some cases +10 mN m−1)140 to generate stimuli-responsive Marangoni flows.141 In a work by Lv et al., the trajectory of polystyrene microparticles present at an air–water interface is manipulated by using such a photoswitch.142 Upon light activation (laser irradiation at 442 nm) of an interface covered with the trans surfactant, local cis isomerisation leads to the transport and trapping of the particles at the beam location. The use of multiple laser beams in close proximity reveals the possibility to manipulate adjacent surface tension gradients to control the positioning of multiple particles.
Self-propulsion of droplets immersed in a micellar continuous phase can also occur when the micelles present in the system are able to partially solubilize the droplet phase. This phenomenon is exemplified by Izri et al. where water droplets are immersed in solutions of squalane and monoolein.145 Monoolein, which acts as a surfactant, ensures the stability of the emulsion by covering the droplets interface. When its concentration exceeds the critical micelle concentration (CMC), surfactant molecules present at the water-squalane interface migrate to the bulk to join the newly formed micelles, which in turn begin to solubilize water in the squalene running phase, thereby displaying the above mentioned symmetry breaking. A surface tension gradient emerges at the droplets interface because they are no longer homogeneously covered with surfactant, and self-propulsed motion is observed with velocities ranging between 10 and 50 μm s−1.
In a similar system, Thutupalli et al. demonstrate that chemical reactions involving monoolein at water droplets interface can be at the origin of self-propulsion and collective behaviour.146 The aqueous droplets, immersed in squalane and monoolein, contain bromine, or reagents leading to bromine formation. At the droplet interface, bromine readily reacts with the monoolein double bond, to form a less active brominated surfactant. The subsequent surface tension imbalance results in active droplets, of which the internal and external velocity fields can be computed by taking into account the relation between surface tension and the bromination density. Collective behaviour, signified here by the polar alignment of the velocities of neighbouring droplets, is observed experimentally and explored theoretically, revealing that the interaction of self-propulsed droplets in a high density environment may lead to the formation of ordered structures.
In a work performed by Babu et al., another type of chemical system is used to trigger the motion of organic droplets.147,148 The authors show that the lipid-forming reaction between 1-hexanethiol and 2-methacryloyloxyethyl phosphorylcholine can be exploited to induce the propulsion of octanol droplets immersed in an aqueous medium (Fig. 10b). The lipid formation is catalysed by the octanol–water interface, leading to a homogeneous lipid coverage of the octanol droplets. Over time, the lipid concentration increases in the medium and they begin to self-assemble into micelles that partially solubilize the octanol in water. The recruitment of lipids into micelles results in a non-homogeneous coverage of the droplets and thus in symmetry breaking which leads to internal Marangoni flow. This work demonstrates that for a specific droplet/surfactant couple, a critical propulsion concentration, beyond the CMC value, needs to be reached to observe motion in the system.
Lastly, a recent study by Cheon et al. demonstrates that the self-propulsion of oil droplets in micellar solutions can be drastically enhanced by the addition of silica particles.149 They hypothesize that the adsorption of these particles at the rear of the droplet interface helps to further polarize the droplet and thereby enhance the Marangoni induced-propulsion. While the precise mechanism is not clarified yet, this effect is observed for a wide range of surfactants and oil and so, opens new perspective in the domain of self-propelled droplets.
In a source-sink configuration, the surface tension gradient is sustained by the interaction between an object that generates a Marangoni flow (source) and another one that absorbs surface active compounds (sink). Our research group has used such dynamics to obtain self-organized patterns between source and drain droplets deposited at an 2D air–water interface.150 Here, source droplets composed of the surfactant tetraethyleneglycol monododecyl ether (C12E4) are used not only to generate a surface tension gradient by releasing the surfactant, but also for their ability to grow centimetre-long, multi-lamellar wire-like structures called myelin figures. We showed that to form patterns, the positioning of these sources can be mediated by the presence of liquid organic droplets that act as sinks by absorbing C12E4 from the air–water interface (Fig. 11a). The resulting Marangoni flow accelerates the growth of the myelin figures from the source droplet, while simultaneously directing them towards the sink droplets. Upon adsorption at the sink, the myelin figures impose a force on the sink droplet, drawing it towards the source until the attractive force is balanced with the repelling Marangoni flow (Fig. 11b and c).151 Along with a kinetic model, we demonstrated experimentally that the composition of the sink determines its capacity for surfactant uptake, and thereby on the self-organization of the myelins-droplet networks.
Fig. 11 Self-organized source-sink system driven by Marangoni flows. (a) Optical microscopy recordings of an air–water interface where a newly deposited oleic acid sink droplet (orange) positions itself between two C12E4 source droplets (red). The interaction between the three floating droplets is mediated by the attractive and repulsive forces originating from the local surface tension gradients at the interface. The source droplets grow myelins and release surfactant molecules in solution, inducing a surface tension gradient. The repelling Marangoni flow that they generate is balanced by the sink droplet that absorbs the surfactant and tethers to the myelins, leading to self-organization into a linear pattern. The scale bar is 2 mm. (b) Simulation of the dynamic surface tension based on C12E4 release rate from the source (Φsource), its absorption rate by the sink (Φdrain) and its desorption from the interface (Φwater). The blue curve represents a situation where the deposition of the sink induces a local surface tension gradient because it absorbs the surfactant faster than it desorbs from the interface. In opposition, the red curve shows a situation where the desorption is faster than the absorption. (c) The schemes present the impact of the dynamic surface tension on the droplets’ positioning. The top, absorption-dominant situation leads to attraction between source and sink. The bottom, desorption-dominant situation leads to repulsion of the sink away from the source, followed by tethering of the myelins to the sink, which keeps the sink in position. Adapted from ref. 150. |
A different type of source-sink configuration is developed by Meredith et al. where organic droplets immersed in an aqueous Triton X-100 micellar solution display predator-prey behaviour (Fig. 12a).152 In this 3D system, the difference in micellar solubilisation rate between droplets of 1-bromooctane (BOct) and ethoxynonafluorobutane (EFB), respectively used as source and sink, is exploited to give rise to complex chase dynamics. In a first step, the BOct predator droplets deplete their surroundings from Triton X-100 by recruiting the surfactant to form BOct-filled micelles (Fig. 12b). This produces a gradient of filled micelles that tends to repel other predator droplets, as the depletion of surfactant generates an asymmetric surfactant occupation at the droplet interface, which drives its motion away from the other BOct droplet. On the contrary, EFB prey droplets are able to have non-reciprocal interactions with their predators because they are resistant to micellar solubilisation and thus induce no repulsion. When EFB preys are in proximity to their BOct predators, they primarily interact with the filled micelle gradient by taking up the BOct that they carry, emptying the previously filled micelles. The adsorption of these newly-formed empty micelles at the front of BOct droplets triggers the internal Marangoni propulsion that is at the origin of the chase, as the prey droplet flees from the predator towards regions with an even lower concentration of solubilized oil. This dynamic behaviour is sustained by the gradient in BOct-filled micelles and can therefore be modulated by the Triton X-100 concentration in the aqueous phase. Increasing the density of predators and prey in a system also leads to self-assembled clusters displaying remarkably complex interactions based on the same non-reciprocal exchanges.
Fig. 12 Self-organized motion in two-body systems. (a) Chase dynamics between BOct (red) and EFB (blue) droplets immersed in a micellar Triton X-100 solution. The BOct predator droplet chases after the EFB prey droplet. Once the predator has caught the prey they keep moving as a pair in the solution. (b) The scheme details the non-reciprocal mechanism behind the chase: (1) BOct is solubilized in the Triton micelles, creating a repelling gradient of filled micelles around it (top scheme). (2) BOct-filled micelles diffuse in the solution and the transported oil is taken up by the EFB droplets, releasing new free micelles at the back of the prey. (3) The free micelles are adsorbed at the front of the predator droplets, triggering the Marangoni propulsion of the predator towards the prey (bottom scheme). (4) The preys also escape the predator by using Marangoni propulsion, moving towards areas rich in empty micelles.152 (c) Self-driven droplet vehicle actuated by the Marangoni effect. The vehicle is made of two droplets (head in yellow, tail in purple) adsorbed on a solid substrate. The surface tension difference between head (high γ) and tail (low γ) generates a convective flow that enables the vehicle to move into a confined space until the tail merges with the head. (d) Table recapitulating the type of vehicle interaction (fusion, delay fusion, ejection, chase) that can be obtained when the surface tension of the tail droplet is decreased while the head droplet is kept constant. (c and d) Reproduced from ref. 154 with permission from the Wiley-VCH publishing group, Copyright 2021. |
Self-driven droplet vehicles developed by Cira et al.,153 and further explored by Li et al.154 are fascinating examples of the versatility present in two-body systems. These vehicles consist of two liquid droplets with distinct surface tensions at their air–liquid interface, that are deposited on a solid substrate, one serving as head and the other as tail (Fig. 12c). When in contact, the convective flow that is generated because of the surface tension difference between the two liquids constitutes the driving force of the vehicles (Fig. 12d), and Li et al.154 demonstrate that this surface tension difference must exceed 5 mN m−1 while the size ratio of the droplets (Dhead/Dtail) must be above 1.6 to display propulsion. Interestingly, they revealed that increasing Δγ did not systematically induce an increase in propulsion velocity and that the latter was in fact determined by the surface tension coefficient /Δγ, where represents the average of the initial γhead and γtail, and Δγ is the difference between the initial γhead and γtail. Finally, they show that by carefully choosing the liquids, the vehicles can efficiently deposit colloidal particles at solid surfaces into patterns at the microscale.
In a work published by Liu et al., decanol droplets are driven through complex pathways at a liquid interface by surface tension gradients induced with organic solvent vapours.155 The decanol droplets are deposited on aqueous ethanol solutions while carbon tetrachloride (CCl4) vapour is released from droplets suspended on a glass substrate placed a few millimetres above the interface (Fig. 13a). The CCl4 vapour locally increases the surface tension and therefore generates Marangoni flows that attract the decanol droplets, towards the vapour-rich area. Intriguingly, the use of multiple droplets at the interface reveals the possibility to generate dynamic droplet patterns. Here, CCl4 is placed on the glass substrate to form the letter T, H or U, and when placed over the interface, the decanol droplets organize according to that shape while maintaining a distance between one another. This out-of-equilibrium process is controlled by the attractive Marangoni flow, driven by the CCl4 vapour, and the repulsive interaction between the decanol droplets, which is driven by the outward Marangoni flow when the droplets are in close contact. Together, this shows promise for the emergence of collective effects in multi-droplets systems.
Fig. 13 Collective behaviour in multi-droplets systems. (a) Decanol droplets at air–water interface self-organize into various letter-shaped patterns under the influence of CCl4 vapours. The latter are used to increase the local surface tension of the interface and subsequently attract the decanol droplets. The patterns are then sustained by the simultaneous action of the Marangoni effect and the repulsion exerted by each droplet. Reprinted (adapted) with permission from ref. 155. Copyright 2019 American Chemical Society. (b) Dynamic self-organized patterns are observed when ethyl salicylate droplets are deposited over an aqueous SDS solution. The initially random arrangement (I) evolves to successively form a ring (II), a chain-like structure (III) and eventually a compact cluster (IV). The organic droplets partially solubilize in the surfactant-rich medium over time, each generating its own Marangoni flow. The transformation of the patterns is therefore understood to be the result of the complex interaction between the local surface tension gradients produced by each individual droplet, and capillary interaction. The scale bar is 10 mm. Reprinted from Colloids and Surfaces A, 566, J. Čejková, K. Schwarzenberger, K. Eckert, S. Tanaka, Dancing performance of organic droplets in aqueous surfactant solutions, 141–147, Copyright 2019, with permission from Elsevier. |
In experiments involving TiO2 particles at a silicon oil–water interface, Wittmann et al. observed Marangoni-induced collective motion upon UV irradiation of the system.156 The particles are adsorbed onto the interface because of the presence of the surfactant Triton X and they catalyse the photo-degradation of H2O2 in the aqueous solution. They report that the spherical particles, which individually display no motility due to their symmetry, are able to collectively move towards an area exposed to UV when they are present as a monolayer, indicating that a collective effect is at the origin of the phenomenon. Throughout their investigation, they rule out possible thermal or phoretic effects and argue that the inhomogeneous repartitioning of the particles at the oil–water interface along with the interface response to the photochemical activity is responsible for symmetry-breaking events leading to a Marangoni flow. The system appears to be maintained in an out-of-equilibrium state by continuous UV exposure, and therefore continuous H2O2 decomposition, as the particles’ motion is interrupted as soon as the irradiation is stopped. Additionally, they demonstrate the existence of a particle density dependence on flow behaviours. In low density settings, the TiO2 particles gather in the irradiated area within seconds and then immediately bounce away from it while in high density settings, they gather and accumulate in the irradiated area to form a stable cluster.
Intriguing collective effects are observed by Tanaka et al., and further explored by Čejková et al. when droplets of decanol or ethyl salicylate (ES) are respectively deposited at the interface of sodium decanoate or SDS aqueous solutions.157,158 The organic droplets first display individual random motion, but then evolve into self-organized dynamic clusters. Three consecutive cluster configurations are showcased in both systems (Fig. 13b). First, a geometrical configuration is obtained where the droplets arrange themselves into a semi-static pattern, circular for ES droplets and squared for the decanol droplets. Over time, these patterns become less stable and eventually collapse, followed by the emergence of a second configuration. Here, the cluster appears to be oscillating because droplets form chain-like structures where they continuously permutate with one another. In the third configuration, the droplets merge via capillary attraction into compact clusters that remain static for hours. The complex interaction between each droplet and the surfactant solution seems to be at the origin of the unique configurations displayed in both systems. Indeed, a study of the decanol droplets in decanoate solution reveals that partial dissolution of the droplets in the medium decreases the local surface tension and thus generates a Marangoni flow away from them. In this context, it is hypothesized that the varying intensity of the Marangoni flow generated from each droplet mediates the interaction, repulsive or attractive, between the members of the cluster.
Fig. 15 Self-organisation driven by diffusiophoresis. (a) Motion pattern of negatively charged tracer particles along a sustained pH gradient between benzoic acid and Ca(OH)2 sources.162 (b) Electroosmotic flow drives the motion of tracer particles amongst CaCO3 microparticles on a negatively charged substrate. Reprinted (adapted) with permission from ref. 163. Copyright 2012 American Chemical Society. (c) Oscillatory clustering and dispersion of Ag3PO4 microparticles in an aqueous H2O2 solution, where reaction (1) and (2) alternate and drive the electroosmotic flow.170 (d) Diffusiophoretic colloidal swarms based on ZnO nanorods and sulfonated polystyrene microbeads, that are coupled via an ion exchange reaction that provide positive and negative feedback to ZnO/polystyrene clustering. (e) Collective swarm behaviour of ZnO/polystyrene particles (orange) that accumulate in the edge of the geometry (grey) with the highest curvature.171 |
Furthermore, McDermott et al.163 reported self-generated flow patterns that emerge from calcium carbonate micropumps, where CaCO3 microparticles slowly dissolve and generate an ionic gradient of Ca2+, HCO3− and OH− ions in the surrounding solution. Due to the high diffusion constant of OH−, an electric field of 1–10 V cm−1 is generated, which is oriented away from the CaCO3 microparticles. This electric field attracts negatively charged sulfate-functionalized polystyrene tracer particles towards the CaCO3 microparticles. Interestingly, as the tracer particles approach the surface underneath the CaCO3 microparticles, they are rapidly ejected radially outward in solution, as shown in Fig. 15b. The flow pattern of the tracer particles exemplifies the electroosmotic flow, which originates from the electric field that is generated around the dissolving CaCO3 microparticles, and exerts its influence on the counterions of the negatively charged substrate.
In addition to particles that move in external concentration gradients, self-organization of dynamic colloidal swarms can emerge from a homogeneous dispersion of particles via autophoresis.160 Here, the colloids self-generate a gradient of ions, which provides the electric field that is required for their propulsion and leads to organized collective behaviour. This concept has been exploited in self-assembly,164,165 dynamic clustering,166 schooling167,168 and predator-prey behaviour of colloidal systems.169 Recently, it has been shown that autophoretic systems can be classified as being either self-diffusiophoretic, when a net ion production occurs, or self-electrophoretic, when the net ion production is zero (i.e., ions are released at one side of the particle, and consumed at the other side).160 Self-diffusiophoretic colloids tend to move slower in a denser population, and form clusters that grow over time, whereas the motility of self-electrophoretic colloid clusters is independent of their density, and results in dynamic clusters with limited sizes.160
Collective, oscillatory behaviour of colloidal swarms has been reported by Altemose et al.170 Silver orthophosphate microparticles (Ag3PO4) are dispersed in an aqueous hydrogen peroxide (H2O2) solution that is exposed to continuous UV illumination. Two opposing reactions occur: (1) reduction of silver ions (Ag+) to silver metal (Ag) under UV illumination, and (2) oxidation of Ag to Ag+ by H2O2 (Fig. 15c). As OH−, which is formed upon the photo-reduction, has a higher diffusion coefficient than the positive ions that are produced, an electric field is generated which is directed away from the Ag3PO4 microparticles. This electric field acts on the positive counterions of the underlying glass substrate, resulting in an electroosmotic flow away from individual Ag3PO4 microparticles, which drives their mutual repulsion in the swarm. However, the re-oxidation (reaction (2)) consumes OH−, and thereby reverses the electroosmotic flow. This causes the Ag3PO4 microparticles to attract each other and form compacted clusters. At the onset of the experiment, microparticles form stable clusters in the presence of H2O2 and UV, and after approx. 5–10 min, instabilities emerge that cause the clusters to disperse. Subsequently, the combination of both reactions causes the microparticles to periodically cluster and disperse. Intriguingly, these oscillations can be synchronized among neighbouring clusters through the diffusion of ions, and oscillation frequencies can be controlled via the H2O2 concentration.
Recently, Wu et al.171 reported diffusiophoretic colloidal swarms inspired by a microbial symbiosis system, where the waste product of one chemical species is the nutrient that drives the motility of the other, and vice versa (Fig. 15d). To this end, ZnO nanorods and sulfonated polystyrene microbeads are coupled via an ion exchange reaction. ZnO nanorods release Zn2+ and OH− in solution, resulting in self-diffusiophoresis that propels the nanorods with velocities up to 2 μm s−1. As the polystyrene microbeads are negatively charged, they are attracted via electrophoresis along the Zn2+ gradient towards the ZnO nanorods. At the same time, the porous polystyrene microbeads take up Zn2+ and release H+via an ion exchange reaction. The H+ gradient that is generated around the negatively charged PS beads causes them to attract each other via an electroosmotic flow at the substrate. Furthermore, the high H+ concentration around the PS bead clusters lowers the local pH, which enhances Zn2+ release as the polystyrene beads approach the ZnO nanorods. Together, the ion exchange between ZnO nanorods and polystyrene microbeads generates a positive feedback that enhances the motility of both. However, close to the glass substrate, underneath the ZnO nanorods/polystyrene microbeads conglomerates, an outward electroosmotic flow is generated, which disassembles the aggregates. Together, the dynamic coupling of the ZnO nanorod/polystyrene microbead system through interconnected ion-exchange reactions allows for hierarchical and dynamic swarm behaviour (Fig. 15e). The clustering of the aggregates is enhanced when high concentrations of Zn2+ and H+ accumulate at the edge of the sample chamber. As a result, the ZnO and polystyrene particles accumulate at all edges and corners present in the geometry of the system. Over the course of the accumulation, the highest curvature confines the ions most successfully. Once a threshold accumulation is reached, the mixture acts as a collective swarm and all colloids migrate towards the sharpest corner – in analogy to the quorum sensing routines displayed by microorganisms.
We anticipate that molecular systems which translate chemical reactions into gradients can strengthen the link between systems chemistry – i.e. the study of multicomponent molecular networks and their emergent properties – and advanced materials that operate at the micro-, meso- or macroscale.173–175 Whereas classical responsive systems such as hydrogels typically respond to a single parameter (e.g. swelling/contraction upon pH change), sustained physicochemical gradients allow integration of multiple chemical input signals into a specific self-organizing response.176,177 This might be exploited to direct the bottom-up construction of morphogenetic materials, for example in interactive biomaterials that adapt their form or texture in concert with bio-molecular cues from the surrounding tissue. Alternatively, self-organization can generate read-out patterns in diagnostic systems that distinguish and report signal inputs derived from complex samples. Furthermore, in chemical systems that solve computational problems, such as maze-solving,178 logic operations179 or pattern recognition,180 the information is typically transferred via diffusion among the ‘computing cores’ where chemical reactions take place that represent the computational steps. Whereas spread of molecular signals via diffusion by itself is uncontrolled, self-organizing gradients can autonomously determine how chemical input signals are guided, for example via reaction–diffusion waves, traveling droplets, growing filaments or organizing colloidal swarms, amongst the computing cores, and ultimately be transferred into material-based responses. Taken together, we foresee that self-organization emerging from chemical systems presents new potential to process complex molecular information into matter-based functional behaviour.
Footnote |
† These authors contributed equally. |
This journal is © the Owner Societies 2022 |