Pawan
Kumar
a,
Dezső
Horváth
b and
Ágota
Tóth
*a
aDepartment of Physical Chemistry and Materials Science, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary. E-mail: atoth@chem.u-szeged.hu
bDepartment of Applied and Environmental Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary
First published on 23rd May 2023
Self-assembly functionalizes active constituents to perform rhythmic activities. Here, our results show that the capillary-Marangoni interaction of irregularly moving gel beads develops complex patterns at the air–liquid interface. The collective behavior of the self-assembled structures exhibits breathing dynamics, polygonal oscillating rings, and cluster synchrony of chains. Interestingly, the trapping of soft particles generates relay synchronization of a rotor. Swarming of clusters is found to form rhythmic shrinking and expanding multiple-ring patterns. The development of self-organized spatiotemporal patterns of our active gel system provides a new way of creating collective oscillations.
Self-assembly is a widely studied phenomenon where the communication between particle–particle or particle–solvent can result in large-scale, ordered architectures. The nature of forces plays a crucial role in determining the final structure and the assemblies. Capillary interaction at the air–liquid interface aggregates the floating or pinned objects found in everyday examples, such as cereal clumping in a milk bowl17 and bubbles clustering in the sink.18 Micro- and nanoparticles can aggregate and create ordered arrays or patterns of particles at interfaces.19
Nowadays, the emergent dynamical self-assembly, where the organized patterns self-propel, has been a focus issue of active matter. Attractive-repulsive forces of a pair of oil droplets can lead to atomic-like motion.20,21 In the context of complex patterns, previous studies of active droplets have identified spontaneous rotating clusters22 and dynamical shuffling of droplets.23 Recent findings have also reported programmable self-assemblies23–25 and dissipative structures to develop temporal aggregation,26 where particles disappear upon switching off the external stimuli. Moreover, synchronization of self-moving patterns and shuffling of packed aggregation have been established, mainly in the presence of stimuli such as light,27 ultrasound,28 and electromagnetic field.29 The simultaneity of self-assembly in conjunction with the periodic oscillations of each constituent, forming collective synchrony of patterns, is hard to achieve in active particles and not much known in the absence of forcing.
In this paper, we present how the collective spatiotemporal periodic oscillations of self-organized structures can synchronize active gels. Furthermore, a simple model is developed to show the collective oscillations of breathing patterns.
The chitosan–gadolinium salt sol [0.1 M GdCl3 in the 0.75 or 1.5 w/v% of CS] in 0.2 M CH3COOH and sodium hydroxide 3.0 M were prepared with deionized water (Purite 300). For active droplets (CGE, chitosan–gadolinium–ethanol), 1.6 mL ethanol was added in the 10 mL of chitosan–gadolinium salt sol with 0.008 w/v% methylene blue to aid visualization. The mixed CGE solution was stirred for 2 hours and stored overnight in room condition. The passive droplets (CG) were prepared by adding 1.6 mL of deionized water in the chitosan–gadolinium salt sol to maintain the same concentration of chitosan and gadolinium compared with the CGE solutions. The density and viscosity of the solutions, summarized in Table S1 in the ESI,† were measured with a density meter (Anton Paar DMA-550) and with a rotational viscometer (Anton Paar ViscoQC-300), respectively. Surface tension measurements were carried out using Krüss tensiometer.
Experiments were performed at temperature 24 ± 1 °C using circular polyvinyl Petri dishes of inner diameter 8.7 cm or 13.7 cm filled with 60 mL or 150 mL of 3.0 M alkaline solution, respectively. CG or CGE droplets were dripped from a height of ≈1 cm on the alkaline surface using a peristaltic pump (Ismatic Reglo) through the Tygon tube with a needle of inner diameter 0.6 mm (or 1.65 mm). Similar behavior was obtained on decreasing and increasing alkaline solution height to 0.5 cm and 1.5 cm, respectively. Those droplets that adhered to the boundaries or fell onto others during dripping were manually removed. During an experiment we have focused our attention on clusters that are sufficiently far form the rest of the beads and the dish wall, i.e., they can be considered isolated. The characterization of these self-assemblies was based on performing 15–20 experiments for each population size of 4–10.
An Osmo action camera was positioned above the container to record the experiments as videos (30 or 60 frames per second with 1080p resolution), illuminated with a diffusive LED and/or white ceiling light source. The size of raw images (1024 × 768 pixels) were captured using a CCD camera equipped with Vivitar lens. Data were analyzed using MatLab scripts and Fiji ImageJ software with MultiTracker or Mtrack2 Plugins.
![]() | (1) |
![]() | (2) |
![]() | ||
Fig. 2 (A) Side view of the CG bead at (i) t = 2 s and (ii) t = 19 s. Field of view: 0.59 × 0.20 cm2. (B) The temporal evolution of bead height h (blue), and diameter d (red). (C) Relative velocity vrel of the CG beads as a function of the ratio of distance from their centroid r and bead radius d/2. The inset images show beads (i) far from or (ii) close to each other, Δt = 1 s. The concentration of CG(sol) is 0.75 w/v% and di = 1.65 mm. The measured data are shown with gray filled circles, the binned data with black hollow circles and their error bars. The red-dashed curve corresponds to the fit of eqn (1) and (2), the blue solid curve to the fit of eqn (1) assuming G(2r/d) = 1. |
We have previously shown that self-propulsion of chitosan hydrogel beads could be realized by adding surface-active material ethanol, as a fuel to the chitosan solution, and we observed continuous and intermittent (stop-and-run) motion of the bead before reaching the final stationary state.36 Therefore to make active droplets, we infused ethanol in the CG solutions. When we place the CGE droplets on the alkaline solution, the ethanol diffuses out of the droplet decreasing the surface tension of the liquid around it. In addition, small asymmetry in the shape of the droplet introduces a spatial gradient in the surface tension, and the arising Marangoni force propel the bead along the surface. Initially, the CGE beads continuously move in an unorganized manner independently of one another at high speed, as depicted in Fig. 3A. Note that the sharp falls in speed occur because of the collisions with the Petri dish wall and the beads, verified by comparing the x–y trajectories with the videos. The depletion of ethanol from the droplet gradually decreases the speed of the beads and produces irregular spikes of run-and-stop motion because the weaker flux of ethanol out of the bead cannot maintain a continuous propulsion, similar to our study with a single chitosan bead.36 At this stage, the beads still remain unsynchronized, illustrated by the speed-time series in the inset Fig. 3A. The cessation of individual run phases indicates that geometric asymmetry in the bead cannot create sufficient surface tension gradient any more. A new source in the gradient of surface tension arises, however, from the crowdedness of beads. The ethanol concentration around the beads will be greater in areas adjacent to many beads, leading to smaller surface tension. This then propels the beads away from each other with speeds that are now comparable to that resultant from the capillary attraction. The interplay of these antagonist forces may result in a self-assembly of the nearest beads.
For the case of two CGE beads, two situations can occur: the beads either move away due to the Marangoni effect or form a pair and self-propel together (see Fig. S2 in the ESI†). For 3–5 beads more complex behavior, synchronization is observed as shown in Fig. 3B for population size of 4 beads, and Fig. S3 (ESI†) for 3 and 5 beads. The inset images in Fig. 3B display the breathing cycle (see Video S1, ESI†), where all beads approach (i and iii) and separate (ii) with their speed is being minimal at the closest and the farthest positions. The beads far from the assembly die in a few seconds, but the organized pattern remains active for a longer duration. Fig. 3C shows that the nearest-neighbor distances of the beads vary irregularly, but the phase remains in sync. The lower threshold of ∼0.63 cm indicates when collisions occur, and the average difference between the maximal and minimal distances of 0.11 ± 0.03 cm provides the repulsion range of the four beads. Our experimental characterization supports that in the assembly state, the repulsion domain is shorter than the range of capillary attraction.
In order to construct a simple model which can recapture the essential dynamics of our experimental system, we have to consider the following. The experimental capillary length lc is 4.5 mm, which yields an Eötvös number E0 = d2/lc2 = 1.2 considering the final size of the beads. This means that surface tension effects are comparable to that of gravitational at the interface around a bead. The transient self-organized oscillation of beads only emerges after the intensive self-propulsion of independent active beads has ceased. By this stage a background ethanol concentration is built in the aqueous phase around the beads acting as sources. Although the drift velocity of the beads is small, the Marangoni number (Ma = ud/D) falls in the range of 2000–8000, indicating that diffusion is negligible on the time scale of significant Marangoni flow with D = 2 × 10−5 cm2 s−1. Hence the temporal evolution of the concentration field can be neglected on the time scale of oscillatory motion, and the beads can be approximated as point sources with constant concentration. The constant course diffusion is characterized by a concentration field with the complementary error function,37 hence, the asymptotic concentration gradient takes the form of exp(−r2), which provides the functional form for the gradient of surface tension and for that of the repelling Marangoni force. The dimensionless governing equation can be expressed as
![]() | (3) |
![]() | (4) |
Besides breathing behavior, more complex dynamics can be observed as factors, such as spatial positions, collision and adhesion/pairing of the beads, influence the spatiotemporal dynamics. For example, the inset image of Fig. 4A shows that beads 2, 4, and 5 form a chain, but 1 and 3 remain free (see Video S3, ESI†). The speed evolution in Fig. 4A illustrates that beads 1–3 oscillate with a higher speed than beads 4 and 5. When beads aggregate, chemical gradients push them away from the central core and capillary interaction organize them into a ring pattern (see Fig. S5 in the ESI†). The self-organized ring pattern can open and close at different positions. The ring shown in Fig. 4B(i) closes and Fig. 4B(ii) opens at the same position resembling inhale and exhale activities (see Video S4, ESI†). Fig. 4B(iii) illustrates that beads 2–7 oscillate with irregular amplitude in synchrony, but bead 1 remains in a nearly steady state or translates with the assembly (x–y trajectories in Fig. S5 in the ESI†). We have also found that the oscillators far from the ring opening mode do not sustain the oscillations for long. Such partial synchrony, where some beads exhibit periodic dynamics in synchrony and the rest stay in a steady state, is also obtained for the self-assembly of three to five beads.
Doubling the CS(sol) concentration in the CGE droplets produces similar dynamic behavior. This concentration increase, however, increases the sol viscosity, resulting in longer bead life because of the stronger cohesive force in the polymer network. An increase in the sol density and a decrease in gel diameter produces stronger capillary forces, hence, greater approaching speed of beads (see Fig. S6 in the ESI†). The longer active propulsion of beads and the arising capillary force together result in various dynamic ordered structures during the transition from initial unordered to the final clustered state. Moreover, the collective oscillations are sustained for a longer duration. Two novel type of patterns form: Fig. 4C(i) represents the chain shape, where beads 3 and 4 adhere, while 1 with 2 and 5 with 6 align at the opposite side. Collectively, they exhibit cluster synchrony, where the outer groups in Fig. 4C(ii) 1,2 and Fig. 4C(iii) 5,6 oscillate with different frequencies, while 3,4 move with the cluster (see Video S5, ESI†). Interestingly, beads oscillate in period-2 dynamics as the outward motion is slow due to slowly increasing and decreasing speed, followed by sharp spikes when the beads move inwards.
Beads' shuffling, collision of a bead to the assembly, or a ring opening–closing process (see Fig. S7 in the ESI†) can induce the trapping of a bead inside a ring of five beads, as depicted in Fig. 4D(i). This self-organized assembly is found to rotate clockwise or anti-clockwise, depending on the initial propel direction of the beads. The speed dynamics in Fig. 4D(ii–iii) show that the beads exhibit oscillatory motion (see Video S6, ESI†). Relay synchronization is also observed, where the speed of bead 2 is in anti-phase with bead 3 but in-phase with bead 4. Note that the anti- and in-phase motion of beads have mild phase drifting, but overall, they maintain the synchronized behavior. When we consider beads 1 and 2, the intermittency of the in-phase and anti-phase region is noticed, and for beads 1 and 5, near anti-phase synchrony is maintained during the oscillatory state of beads (see Fig. S7 in the ESI†). The depletion of ethanol with time decreases the speed of the rotor, which distorts the phases of oscillations in the later stage. We have observed rotating oscillatory pattern only in the case of six beads. Trapping is not sustained for a population size of less than six, while for larger population (N > 6), beads in the ring are crowded and do not have sufficient space to generate the oscillatory rotation. In the case of multiple beads, various clusters form and oscillate synchronously, as shown in Fig. S8 (ESI†). The size of the larger clusters (see Fig. S9, ESI†) increases with time due to the capillary attraction among the beads and smaller clusters.
We have also explored the communication between two groups of gel beads. The interaction of inactive CG particles shows that both groups remain intact in the static assembly after collision (see Fig. S10 in the ESI†). On the other hand, the elongated, elliptical-like configurations of the active groups with population ratios 18:
15 and 39
:
90—shown by red:blue dots in Fig. 5A(i) and (ii), respectively—transform into circular shape (see Fig. 5A(iii and iv)), indicated by the black dotted lines. In the two groups of beads (red and blue), pairwise interactions dominate as the beads are distributed along chain-like formations. Furthermore, the swarm radii of the particles—defined as
, where N is the number of beads,
i their position vectors, and
m the position of the collective center of mass—decreases due to clustering of the two assemblies (see Fig. 5B). The disk-like shape of the collective group with multiple ring patterns, as shown in the inset images of Fig. 5B(i–iii), exhibits oscillatory behavior. The decrease in the swarm radii indicates disk shrinking, and the expansion of the disk occurs on increasing Rs (see Video S7, ESI†). Similar behavior is obtained for an assembly of 33 gel beads (see Fig. S11 in the ESI†). During these activities, the active motion of the particles can shuffle the ring connections. When the motion ceases and the system reaches the near equilibrium state, the final configuration of the disk develops multiple islands devoid of beads. Moreover, the 2-D Voronoi diagram using the centroid position of the beads (see Fig. 5C), indicates that polygonal cells with 4–8 edges spread all over the space. Fig. 5D shows that pentagonal and hexagonal cells dominate. We have performed different sets of experiments for population sizes of 99–163 particles, and in all cases a similar probability distribution is observed as shown in Fig. S12 (ESI†). The total Voronoi entropy
is 1.59, where P(s) is the probability of polygons with s edges. For a fully regular structure Svor = 0, and for random distribution of points Svor = 1.71 has been reported.38 In our case, the high Svor characterizes a disordered pattern.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3sm00461a |
This journal is © The Royal Society of Chemistry 2023 |