Next Article in Journal
Experimental Study on Critical Parameters Degradation of Nano PDSOI MOSFET under TDDB Stress
Next Article in Special Issue
Anisotropic Diffusion of Elongated Particles in Active Coherent Flows
Previous Article in Journal
Analysis and Optimization of Dynamic and Static Characteristics of the Compliant-Amplifying Mechanisms
Previous Article in Special Issue
Dielectrophoresis-Based Selective Droplet Extraction Microfluidic Device for Single-Cell Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy by Neighbor Distance as a New Measure for Characterizing Spatiotemporal Orders in Microscopic Collective Systems

1
University of Michigan—Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, Shanghai 200240, China
2
Department of Mechanical Engineering, Carnegie Mellon University, 5000 Forbes Ave, Pittsburgh, PA 15213, USA
3
College of Engineering, University of Michigan, Ann Arbor, MI 48109, USA
4
The Academy for Engineering and Technology, Fudan University, Shanghai 200433, China
5
Physical Intelligence Department, Max Planck Institute for Intelligent Systems, 70569 Stuttgart, Germany
6
Department of Physics, University of Stuttgart, 70569 Stuttgart, Germany
7
Department of Physics, Banaras Hindu University, Varanasi 221005, India
8
Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Forschungszentrum Jülich, 52425 Jülich, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Micromachines 2023, 14(8), 1503; https://doi.org/10.3390/mi14081503
Submission received: 15 June 2023 / Revised: 21 July 2023 / Accepted: 24 July 2023 / Published: 26 July 2023
(This article belongs to the Special Issue Active Colloidal and Micro Systems for Propulsion and Manipulation)

Abstract

:
Collective systems self-organize to form globally ordered spatiotemporal patterns. Finding appropriate measures to characterize the order in these patterns will contribute to our understanding of the principles of self-organization in all collective systems. Here we examine a new measure based on the entropy of the neighbor distance distributions in the characterization of collective patterns. We study three types of systems: a simulated self-propelled boid system, two active colloidal systems, and one centimeter-scale robotic swarm system. In all these systems, the new measure proves sensitive in revealing active phase transitions and in distinguishing steady states. We envision that the entropy by neighbor distance could be useful for characterizing biological swarms such as bird flocks and for designing robotic swarms.

1. Introduction

Active matter systems consisting of self-propelled or driven individuals can self-organize and form large-scale ordered spatiotemporal structures [1,2]. They serve as model systems for non-equilibrium mechanics as their individual units break the time-reversal symmetry [3] and local detailed balance [4]. They are ubiquitous across many scales, from bird flocks [5], fish schools [6], and ant colonies [7] to bacteria colonies [8,9], artificial micro-nano machines [10,11,12,13,14,15,16,17], and molecular motors [18,19]. Though these collective systems differ widely in the sizes and levels of intelligence of their individual units, they share the common characteristics of self-organization, i.e., local rules and behaviors determining global properties and functions. Therefore, one of the fundamental challenges in the study of collective systems is constructing appropriate measures based on local parameters that reflect global properties or functions. Addressing this challenge will contribute to our understanding of the principles of self-organization in all collective systems and is meaningful for both understanding the living systems and designing artificial systems.
Collective systems are often analyzed from the perspective of the statistical mechanics of particles. In this approach, the change in the spatiotemporal patterns of a collective system can be characterized by an order parameter, akin to the characterization of phase transition in equilibrium systems. In the seminal work of Vicsek, the average velocity was used as the order parameter to study the phase transition of collectives of self-propelled particles [20]. This average velocity is the vectorial mean of individual self-propelled particles and thus reflects the degree of alignment of the whole collective. This average velocity, however, does not reflect the degree of order based on the positions of individual units. Positional orders of static structures are mostly analyzed according to the symmetries of point groups, a typical example being crystal structures. Analyzing dynamic spatiotemporal patterns of collective systems poses challenges for this approach because the relative positions between individual units vary over time.
An alternative approach is to analyze the spatiotemporal patterns from the perspectives of information. This approach has a long history but is under resurgence recently. The connection between structure and information can be traced back to the idea of an “aperiodic crystal” in living systems by Erwin Schrödinger, which foretold the discovery of DNA [21]. Following the ideas of Kolmogorov–Chaitin (KC) complexity, proposed in the 1960s to quantify the amount of information that one computer program carries, Martiniani et al. recently proposed computable information density (CID) to measure the information contents contained in the structures of non-equilibrium systems [22]. Following the spirit of KC complexity, Mackay and Cartwright noted that structures also have repetition and redundancy, just the same as computer programs, and they proposed assembly complexity to measure the amount of information a structure carries [23]. Other researchers seek to apply the formalism of Shannon entropy to characterize structures. For example, Frenkel et al. proposed Voronoi entropy (VE) to quantify the symmetry and orderliness of two-dimensional microdroplet clusters and compared it with continuous symmetry measure (CSM) [24]. As another example, the Shannon entropy could be used to quantify the degree of order in crystallography and the complexity of the inorganic crystal [25,26]. In addition, information theory has been widely applied to crystallography and has given rise to chaotic crystallography [27]. Building on these ideas, we recently proposed the entropy by neighbor distance H N D i s t as a measure of the order of the collective pattern generated by micron-scale spinning rafts driven by magnetic fields [28]. This measure can distinguish various patterns of hundreds of units at different rotating frequencies of magnetic fields and is even sensitive to different local symmetries of individual rafts. As a natural extension, we want to ask the following question: How well does this new measure generalize to other collective systems?
Here, we extend the new measure entropy by neighbor distance H NDist to three types of systems: a simulated self-propelled boid system, two active colloidal systems, and one centimeter-scale robotic swarm system. In the simulated boid system, we compare H NDist with the average velocity v a in characterizing the change in spatiotemporal patterns. We compare the situation where neighbors are defined by metric pairwise distances with the situation where neighbors are defined by topology and find that H NDist is more sensitive for the latter case. In active colloidal systems, we demonstrate that H NDist correlates well with the system’s response to external stimuli and can also distinguish different steady-state patterns. We analyze the impact of bin size in the calculation of H NDist . In the centimeter-scale robotic swarm system, we show that H NDist can track the overall progress of self-organization and distinguish different steady-state patterns. Moreover, the analysis of circular agents assembled into the configuration of the 26 letters in the alphabet shows the shape dependence of H NDist . We expect that these results will further stimulate efforts to apply H NDist in other areas, particularly in biological systems and microrobotic systems.

2. Results and Discussion

2.1. Simulated Boid System

We adopt a discrete-agent model with point-like boids whose dynamics are governed by three rules: cohesion, separation, and alignment [29]. Cohesion and separation cause the attractive and repulsive interactions between neighboring boids, respectively, and the alignment rule causes a boid to align its heading direction with the average of its neighbors. The sum of these interactions is scaled by a coefficient S, the steering factor, and then added to the velocity of the boid in order to get the new velocity of the boid at the next time step. The steering factor S adjusts the relative weight between the influence of the neighbors and the inertial of the boid. We use two methods to define neighbors in our simulation. The first one is through metric distance: If two boids are within a certain threshold distance, they are neighbors. The second one is through topology: We construct Voronoi tessellations to determine neighboring pairs. The initial position and velocity direction of all boids are random (assuming uniform distribution throughout the entire interval), and the amplitude of the velocity is kept constant. At each step of the update on velocity, we introduce noise in the direction of velocity to investigate the effect on the collective patterns. A detailed description of the model of the simulation is included in the method section. Two representative videos are included in the Supplementary Video S1.
First, we use the result of a representative single simulation run to demonstrate the temporal evolution of the collective patterns (Figure 1). This run is performed with the condition of zero noise, and the neighbors are defined by topology. Figure 1a shows the patterns of boids at several key time steps. In the beginning, the positions and the directions of velocities of boids are random. As the simulation progresses, the directions of the boids are gradually aligned. Figure 1b,c show the quantification of these patterns of boids. Figure 1c shows that v a starts at zero in the initial state and continues to increase until reaching near the maximum value of 1 after ~40,000 steps, while the topological structures revealed by H NDist do not appear to change much before 40,000 steps, and the value of H NDist only starts to decrease after 40,000 steps, and then reaches a plateau after 80,000 steps (Figure 1b). This change in H NDist is reflected by the increase in uniformity of the distribution of boids in the patterns. The comparison of Figure 1b,c shows that the velocity alignment is faster than the change in topology structures in the patterns of the boids model. It demonstrates that the order of the velocity direction and the topological structure are not equivalent, so these two measures characterize different aspects of the collective patterns.
Next, we explore the effect of noise in boid velocity updates and the effect of two different methods of defining neighbors on the characterization of collective patterns by H NDist and v a (Figure 2). Noise is added to the angle component of the velocity update at each time step, v i t + Δ t = V v ^ + θ i + Δ θ θ ^ , where V is the constant speed of all boids, and a random variable Δ θ is added to the angular component of the velocity. Then, Δ θ is drawn from a uniform distribution η , η , and η is the noise factor that is varied from 0 to 7°. We run the simulations with different noise factors and take the data of the last 1000 steps when the simulations reach steady states.
In the situation where the neighbors are defined by topology (Figure 2a), after the final steady-states are reached, the magnitude of the averaged velocity v a decreases as the noise factor η increases, whereas entropy by neighbor distance H NDist   increases as η increases. The onset of the change in v a and H NDist also differ: H NDist is sensitive to the initial small noise, whereas v a is robust to the small noise and only decreases appreciably when η reaches 3 ° ~ 4 ° . Intuitively, noise influences the velocity direction directly, but the small noise does not greatly change the velocity alignment, so v a is relatively stable with the small noise. However, because H NDist is based on the topological structure, a small fluctuation may cause the change of the topology. For example, the number of neighbors in a hexagonal lattice could change from 6 to 5 or 7 because of a small shift in the position of one vertex. This sensitivity explains why H NDist increase dramatically at the beginning. Therefore, these two measures are again complementary: one for the topological order and the other for the alignment order.
In situations where neighbors are defined by metric distances (Figure 2b), v a still exhibits a significant decrease as the noise factor η increases, whereas H NDist shows little changes as η increases. We suspect that the reason for the negligible change in H NDist is that boids may interact with a lot more neighbors when the neighbors are defined by metric distances than when the neighbors are defined by topological structures, and hence they tend to form big clusters in the final steady state. This tendency of cluster formation is relatively robust to the disturbance of noise. Therefore, H NDist shows negligible changes.

2.2. Active Colloids

Active colloidal systems include self-propelled microparticles [30,31,32] and external field-driven microparticles [33,34,35]. They not only serve as model systems for research on non-equilibrium physics [36,37,38] but also are used as microrobots for biomedical and environmental applications [39,40,41]. One of the striking phenomena in active self-propelled particles is that they could experience so-called motility-induced phase separation (MIPS) with purely repulsive interactions, which is impossible for passive colloidal particles without attraction. This phenomenon has been observed in both simulations and experiments [7,42,43].
One of the representative models for MIPS is the light-activated colloidal particles. For example, Palacci et al. used a suspension of synthetic photoactivated colloidal particles [44], and they observed that when the blue light is on, homogeneously distributed particles began to assemble into clusters with an average size of 35 particles (~10 μm). They claim that the osmotically driven motion and steric hindrances (collisions) are necessary for the formation of “living crystals”. We tracked the particles in a video of living crystals and performed analysis using H NDist (Figure 3a and Supplementary Video S2). H NDist captures the three stages of the whole process: clustering (light on), dispersed (light off), and again clustering (light on). Turning on the light causes the particles to move, which induces MIPS, making the structure more ordered and clustered, and thus causing H NDist to decrease, and the reverse is true for the dispersed phase.
Further, we explore the influence of bin sizes on calculating H NDist . Figure 3b shows that the trends differ when the bin sizes are smaller or larger than the radius R of the particles (3 px). For example, at 23 s (light off) the H NDist increases or decreases for the bin sizes smaller or larger than R, respectively. In addition, the downward trend that appeared in 34 s (light on) is not well reflected in the three curves whose bin sizes are larger than R. The analysis of H NDist , norm also shows that the trend in the three curves whose bin sizes are smaller than R is similar (Figure 3c).
To explain the effect of the bin sizes on H NDist , we check the histograms. When the bin size is too large, small changes in separation distances are not accounted for, resulting in H NDist being not sensitive enough, such as the comparison before and after 23 s (light off). Figure 3d shows that after 23 s, the distribution near the first peak increases when the bin size is 1 px, which contributes to the increase of H NDist , whereas it is difficult to capture these small changes in distributions when the bin size is 10 px. Furthermore, by comparing the histograms at the beginning and at the end, we find that the initial neighbor distance distribution is concentrated in the range below 100 px. If the bin size is too large, the corresponding H NDist of the initial dispersed state will be smaller than the final aggregated state, which is inconsistent with our intuition. Therefore, we think that bin sizes smaller than the radius of particles are good choices for the analysis of collective microscopic patterns.
Figure 3. The analysis of light-activated colloidal surfers system using entropies by neighbor distances. (a) Temporal evolution of H NDist (bin size = 1 px). (b) H NDist of different bin sizes (px). (c) Temporal evolution of H NDist , norm with different bin sizes (px). (d) Histogram of neighbor distance distribution at four critical moments with bin size 1 and 10 px. The original video is from the Ref. [44].
Figure 3. The analysis of light-activated colloidal surfers system using entropies by neighbor distances. (a) Temporal evolution of H NDist (bin size = 1 px). (b) H NDist of different bin sizes (px). (c) Temporal evolution of H NDist , norm with different bin sizes (px). (d) Histogram of neighbor distance distribution at four critical moments with bin size 1 and 10 px. The original video is from the Ref. [44].
Micromachines 14 01503 g003
Another typical type of experimental model of active colloids is the electric-field-activated colloidal particles. Snezhko et al. used pear-shaped dielectric particles confined inside a cylindrical cell to introduce chiral rollers with activity-controlled curvatures of their trajectories and spontaneous handedness of their motion [45]. By controlling activity through variations of the static (DC) electric field applied perpendicular to the bottom surface of the cell, they showed emergent dynamic phases, ranging from a gas of spinners to aster-like vortices. We tracked the particles in two of their published videos and demonstrate that H NDist could distinguish a gas-like steady state and an aggregated steady state easily (Figure 4a and Supplementary Video S3).
Furthermore, the two states remain distinguishable as long as the bin sizes are smaller than the radius of the particles (R = 10 px), which reflects the robustness of H NDist , norm in distinguishing patterns with large differences in their orders (Figure 4b). In examining the non-normalized H NDist , we find that H NDist decreases with increasing bin sizes (Figure 4c). Different bin sizes cause the overlapping of curves of the two states, and the fluctuation of non-normalized H NDist of the gas state is greater than that of the aggregated state. Figure 4d shows that the standard deviation (std) of the gas state is nearly three times that of aggregated state. In addition, the values of standard deviations are stable under different bin sizes, which provide additional support for choosing the bin size smaller than R, as the previous case suggests.
Figure 4. The analysis of active chiral fluids with two different states using entropies by neighbor distances. (a) Temporal evolution of H NDist , norm . (b) H NDist , norm of different bin sizes (px): solid lines represent the gas state and dotted lines represent the aggregated state. (c) Non-normalized H NDist of different bin sizes (px): solid lines represent the gas state and dotted lines represent the aggregated state. (d) The standard deviations (std) of H NDist of the time series. The original videos are from the reference [45].
Figure 4. The analysis of active chiral fluids with two different states using entropies by neighbor distances. (a) Temporal evolution of H NDist , norm . (b) H NDist , norm of different bin sizes (px): solid lines represent the gas state and dotted lines represent the aggregated state. (c) Non-normalized H NDist of different bin sizes (px): solid lines represent the gas state and dotted lines represent the aggregated state. (d) The standard deviations (std) of H NDist of the time series. The original videos are from the reference [45].
Micromachines 14 01503 g004

2.3. Robotic Swarm

Swarm robotics is the study of a large number of robots that can spontaneously determine their next move by exchanging information, such as their position and speed, with neighboring individuals, and self-organize into the desired configuration without sending and receiving direct control commands from external supervisors such as a base station or a human overseer [46]. This bionic intelligent control system takes inspiration from the self-organized behaviors of social animals and has important applications in many fields ranging from farming and food management, to military defense and space systems [47].
Nagpal et al. demonstrate a thousand-robot swarm (dubbed kilobots) capable of large-scale, flexible self-assembly of two-dimensional shapes entirely through programmable local interactions and local sensing [48]. Specifically, they developed a collective algorithm that enables 1024 robots, each with limited capabilities, to cooperatively assemble into some pre-designed shapes. We picked two of their published videos: the assembly of the letter K-shape and starfish shape, and we tracked the robots during the self-assembly process (Figure 5a,b and Supplementary Video S4). H NDist show a continuous gradual increase in both self-assembly processes and thus can be a sensitive measure for the progress of the programmable self-assembly. Moreover, the final values of H NDist for the two shapes differ, which suggests H NDist may be used as an indicator for reconfigurable kilobot self-assemblies as well.
To further demonstrate that H NDist is shape-dependent, we simulate robot swarms assembled into configurations of 26 letters in the alphabet. We first create a densely packed lattice composed of circles of equal diameter, and then obtain the configurations of letters by masking the lattices with corresponding shapes (Figure 5c). These letters show different entropy values (Figure 5d). These differences are due to the circles along the edges: Even though the configurations are all based on the square lattice, their edge profiles are quite different, and the units on edges can show significantly different neighbor distances, resulting in different neighbor distance distributions. This finding may be useful for the design of control algorithms for robotic swarms.
Figure 5. The analysis of the self-assembly process in the kilobot swarm using H NDist . (a) K-shape. (b) Starfish shape. The original videos are from the Ref. [48] (c) simulated robot swarm assembled into shapes of letters. (d) H NDist of 26 letters in the alphabet.
Figure 5. The analysis of the self-assembly process in the kilobot swarm using H NDist . (a) K-shape. (b) Starfish shape. The original videos are from the Ref. [48] (c) simulated robot swarm assembled into shapes of letters. (d) H NDist of 26 letters in the alphabet.
Micromachines 14 01503 g005

3. Conclusions and Discussion

In summary, we extend H NDist into three types of active systems and prove its usefulness in characterizing the spatiotemporal patterns of these systems. In the active colloid systems and robotic swarm systems, we have analyzed published videos of experiments and showed H NDist as a useful indicator of pattern changes and a measure to distinguish different steady states. In the simulated boid system, we showed that H NDist is a sensitive measure of phase transitions and can reveal the existence of noise. We further found that H NDist is more sensitive to situations where neighbors are defined by the topological structures instead of metric distances. Field studies of starling flocks show that the interactions of birds depend on topological distance rather than metric distance, and each bird only interacts with six to seven neighbors on average [49]. Accordingly, we surmise that H NDist could be a sensitive measure for biological systems, including bird flocks.
We think that the reason why H NDist could accurately reveal the changes of patterns in a variety of systems relates to the dominance of pairwise contribution to the configuration entropy in a system. The configuration entropy of a system contains contributions from multi-body correlation functions, and Duane Wallace et al. have shown that in simple liquids the two-body excess entropy, which only involves the pair correlation function, accounts for more than 90% of the configuration entropy [50,51,52]. We are currently investigating the fundamental relation between H NDist and two-body excess entropy, and the results will be published in the future.
We also anticipate useful applications of H NDist in swarm robotics because a local measure such as H NDist should in principle be very sensitive to the self-organization of robots who only interact locally. For example, in recent years researchers have been developing decentralized algorithms for robotic control so that they could self-organize into certain patterns through local perception and communication between individuals to reduce the data transmission and communication cost between individuals and a central control station [46,53,54]. H NDist could be used as an important measure for characterizing the overall configuration and for fine-tuning local interaction to achieve specific functions.

4. Methods

4.1. Boids Simulation

We adopted a simplified discrete-agent model with point-like boids whose dynamics are governed by the three rules: cohesion, separation, and alignment. Specifically, cohesion means each boid will tend to gather together and move to the center of its neighbors. Separation means if two boids get too close, they will tend to stay away from each other. The alignment refers to boids will move along the same direction as their neighbors. The size of the simulation box is 1000 × 1000 pixels. We employ a periodic boundary condition.
Algorithm
The interaction between boids includes separation, alignment and cohesion force:
F i , sep = j φ r i r j r i j 2 F i , alg = j φ v j n i v i F i , coh = j φ r j n i r i
where φ is the set containing all neighbors of i th boid.
To eliminate the influence of the absolute value of three kinds of interaction, we normalized these three parameters to obtain the total force:
F i = F i , sep F i , sep + F i , alg F i , alg + F i , coh F i , coh
Thus, we define the direction of movement at the moment t + Δ t as:
θ i t + Δ t = arctan v i t + S F i v i t + S F i
where S is the steering factor, which controls the relative weight between the influence from neighbors and the influence from the velocity of the boid. We tune S in the range from 0.001 to 0.1.
The velocity at moment t + Δ t has contributions from v ^ and θ ^ ,which represent the radial and angular components separately:
v i t + Δ t = V v ^ + θ i + Δ θ θ ^   Δ θ η , η
where V is the constant speed, Δ θ is a random variable drawn from a uniform distribution, and η is the noise factor in degree.
Then we could get the exact position of i th agent for the next moment by:
x i t + Δ t = x i t + v i t + Δ t Δ t
Definition of neighbors who has mutual interaction
We defined neighbors who have influence on each other in two different ways. For metric distance standard in two-dimension, the set of neighbors φ includes all units within the circle whose radius is R center at particle i . The second way is defined by topology structure. Particularly, we perform Voronoi tessellation, and the neighbors are defined by a convex polygon.

4.2. Data Analysis

Calculation of average velocity  v a
The average velocity is the sum of the velocity of all boids:
v a = 1 N V i S v i
where S is the set of all boids.
Calculation of Entropy by neighbor distance  H N D i s t
In calculating H N D i s t , the neighbor always refers to topological neighbors. After determining neighbors by Voronoi tessellation, we obtain the distance between each agent and its neighboring agents. Then we count the neighbor distances of all agents to obtain the neighbor distance distribution of the overall configuration. We apply the formular of the Shannon entropy to calculate the entropy by neighbor distance
H N D i s t = i p i log 2 p i
where p i = X i / X is the probability of a neighbor distance that falls within a distance interval (a bin) labeled by index i , X is the total count of all neighbor distances of all agents, and X i is the count of the neighbor distances in bin i.
For the convenience of comparison, the normalized H N D i s t was introduced as:
H N D i s t , n o r m = i p i log 2 p i w i l o g 2 i w i
where w i is the bin size.
Particle tracking
The image sequences were analyzed to investigate the structure and dynamics of the self-propelled active colloidal particles and robotic swarm. Particle tracking was performed by Python package Deeptrack (https://github.com/softmatterlab/DeepTrack2 (accessed on 23 July 2023)) [55] and Trackpy (https://soft-matter.github.io/trackpy/dev/index.html (accessed on 23 July 2023)) to identify the 2D location of each agent. The xy coordinates of the agents were used to calculate H N D i s t and characterize the structural order in the systems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/mi14081503/s1, Video S1. Boid simulation. Video S2. H N D i s t analysis of a video of “living crystals”. Video S3. H N D i s t analysis of a video of active chiral spinners. Video S4. H N D i s t analysis of two videos of kilobots self-organization.

Author Contributions

Conceptualization, W.W., G.G., V.K. and P.M.; simulation, Z.W. and S.Z.; video processing, Y.F. and J.Y.; data analysis and summary Y.F., writing—original draft preparation, Y.F.; writing—review and editing, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Commission of Shanghai Municipality, project number 23ZR1433700, and the UM-SJTU JI start-up fund. This research was partially funded by SERB India via the grant SRG/2022/001088 and IoE BHU.

Data Availability Statement

We have uploaded our code to Github (https://github.com/wowoshu/boids (accessed on 23 July 2023)).

Acknowledgments

Y.F. thank Yuan He, Yi Zhang, and Zhenkun Cai for the inspiring discussion on video processing. G.G. thanks Max Planck society for support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bowick, M.J.; Fakhri, N.; Marchetti, M.C.; Ramaswamy, S. Symmetry, Thermodynamics, and Topology in Active Matter. Phys. Rev. X 2022, 12, 010501. [Google Scholar] [CrossRef]
  2. Marchetti, M.C.; Joanny, J.F.; Ramaswamy, S.; Liverpool, T.B.; Prost, J.; Rao, M.; Simha, R.A. Hydrodynamics of Soft Active Matter. Rev. Mod. Phys. 2013, 85, 1143–1189. [Google Scholar] [CrossRef] [Green Version]
  3. Lauga, E. Life around the Scallop Theorem. Soft Matter 2011, 7, 3060–3065. [Google Scholar] [CrossRef] [Green Version]
  4. Battle, C.; Broedersz, C.P.; Fakhri, N.; Geyer, V.F.; Howard, J.; Schmidt, C.F.; MacKintosh, F.C. Broken Detailed Balance at Mesoscopic Scales in Active Biological Systems. Science 2016, 352, 604–607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Vicsek, T.; Zafeiris, A. Collective Motion. Phys. Rep. 2012, 517, 71–140. [Google Scholar] [CrossRef] [Green Version]
  6. Makris, N.C.; Ratilal, P.; Symonds, D.T.; Jagannathan, S.; Lee, S.; Nero, R.W. Fish Population and Behavior Revealed by Instantaneous Continental Shelf-Scale Imaging. Science 2006, 311, 660–663. [Google Scholar] [CrossRef]
  7. Anderson, C.; Fernandez-Nieves, A. Social Interactions Lead to Motility-Induced Phase Separation in Fire Ants. Nat. Commun. 2022, 13, 6710. [Google Scholar] [CrossRef]
  8. Worlitzer, V.M.; Jose, A.; Grinberg, I.; Bär, M.; Heidenreich, S.; Eldar, A.; Ariel, G.; Be’er, A. Biophysical Aspects Underlying the Swarm to Biofilm Transition. Sci. Adv. 2022, 8, eabn8152. [Google Scholar] [CrossRef]
  9. Zhang, H.P.; Be’er, A.; Florin, E.-L.; Swinney, H.L. Collective Motion and Density Fluctuations in Bacterial Colonies. Proc. Natl. Acad. Sci. USA 2010, 107, 13626–13630. [Google Scholar] [CrossRef]
  10. Yu, J.; Wang, B.; Du, X.; Wang, Q.; Zhang, L. Ultra-Extensible Ribbon-like Magnetic Microswarm. Nat. Commun. 2018, 9, 3260. [Google Scholar] [CrossRef] [Green Version]
  11. Zheng, J.; Chen, J.; Jin, Y.; Wen, Y.; Mu, Y.; Wu, C.; Wang, Y.; Tong, P.; Li, Z.; Hou, X.; et al. Photochromism from Wavelength-Selective Colloidal Phase Segregation. Nature 2023, 617, 499–506. [Google Scholar] [CrossRef]
  12. Bricard, A.; Caussin, J.-B.; Desreumaux, N.; Dauchot, O.; Bartolo, D. Emergence of Macroscopic Directed Motion in Populations of Motile Colloids. Nature 2013, 503, 95–98. [Google Scholar] [CrossRef] [Green Version]
  13. Meredith, C.H.; Castonguay, A.C.; Chiu, Y.-J.; Brooks, A.M.; Moerman, P.G.; Torab, P.; Wong, P.K.; Sen, A.; Velegol, D.; Zarzar, L.D. Chemical Design of Self-Propelled Janus Droplets. Matter 2022, 5, 616–633. [Google Scholar] [CrossRef]
  14. Che, S.; Zhang, J.; Mou, F.; Guo, X.; Kauffman, J.E.; Sen, A.; Guan, J. Light-Programmable Assemblies of Isotropic Micromotors. Research 2022, 2022, 9816562. [Google Scholar] [CrossRef]
  15. Li, L.; Yu, Z.; Liu, J.; Yang, M.; Shi, G.; Feng, Z.; Luo, W.; Ma, H.; Guan, J.; Mou, F. Swarming Responsive Photonic Nanorobots for Motile-Targeting Microenvironmental Mapping and Mapping-Guided Photothermal Treatment. Nano-Micro Lett. 2023, 15, 141. [Google Scholar] [CrossRef]
  16. Mou, F.; Zhang, J.; Wu, Z.; Du, S.; Zhang, Z.; Xu, L.; Guan, J. Phototactic Flocking of Photochemical Micromotors. iScience 2019, 19, 415–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Mou, F.; Li, X.; Xie, Q.; Zhang, J.; Xiong, K.; Xu, L.; Guan, J. Active Micromotor Systems Built from Passive Particles with Biomimetic Predator–Prey Interactions. ACS Nano 2020, 14, 406–414. [Google Scholar] [CrossRef] [PubMed]
  18. Harada, Y.; Noguchi, A.; Kishino, A.; Yanagida, T. Sliding Movement of Single Actin Filaments on One-Headed Myosin Filaments. Nature 1987, 326, 805–808. [Google Scholar] [CrossRef]
  19. Badoual, M.; Jülicher, F.; Prost, J. Bidirectional Cooperative Motion of Molecular Motors. Proc. Natl. Acad. Sci. USA 2002, 99, 6696–6701. [Google Scholar] [CrossRef] [PubMed]
  20. Vicsek, T.; Czirók, A.; Ben-Jacob, E.; Cohen, I.; Shochet, O. Novel Type of Phase Transition in a System of Self-Driven Particles. Phys. Rev. Lett. 1995, 75, 1226–1229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Schrödinger, E. What Is Life? Cambridge University Press: Cambridge, UK, 1944; ISBN 978-1-107-60466-7. [Google Scholar]
  22. Martiniani, S.; Chaikin, P.M.; Levine, D. Quantifying Hidden Order out of Equilibrium. Phys. Rev. X 2019, 9, 011031. [Google Scholar] [CrossRef] [Green Version]
  23. Cartwright, J.H.E.; Mackay, A.L. Beyond Crystals: The Dialectic of Materials and Information. Phil. Trans. R. Soc. A 2012, 370, 2807–2822. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Frenkel, M.; Fedorets, A.A.; Dombrovsky, L.A.; Nosonovsky, M.; Legchenkova, I.; Bormashenko, E. Continuous Symmetry Measure vs Voronoi Entropy of Droplet Clusters. J. Phys. Chem. C 2021, 125, 2431–2436. [Google Scholar] [CrossRef]
  25. Krivovichev, S. Topological Complexity of Crystal Structures: Quantitative Approach. Acta Cryst A 2012, 68, 393–398. [Google Scholar] [CrossRef] [PubMed]
  26. Krivovichev, S.V. Which Inorganic Structures Are the Most Complex? Angew. Chem. Int. Ed. 2014, 53, 654–661. [Google Scholar] [CrossRef] [PubMed]
  27. Varn, D.P.; Crutchfield, J.P. Chaotic Crystallography: How the Physics of Information Reveals Structural Order in Materials. Curr. Opin. Chem. Eng. 2015, 7, 47–56. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, W.; Gardi, G.; Malgaretti, P.; Kishore, V.; Koens, L.; Son, D.; Gilbert, H.; Wu, Z.; Harwani, P.; Lauga, E.; et al. Order and Information in the Patterns of Spinning Magnetic Micro-Disks at the Air-Water Interface. Sci. Adv. 2022, 8, eabk0685. [Google Scholar] [CrossRef] [PubMed]
  29. Reynolds, C.W. Flocks, Herds and Schools: A Distributed Behavioral Model. SIGGRAPH Comput. Graph. 1987, 21, 25–34. [Google Scholar] [CrossRef] [Green Version]
  30. Bär, M.; Großmann, R.; Heidenreich, S.; Peruani, F. Self-Propelled Rods: Insights and Perspectives for Active Matter. Annu. Rev. Condens. Matter Phys. 2020, 11, 441–466. [Google Scholar] [CrossRef] [Green Version]
  31. Pan, Q.; He, Y. Recent Advances in Self-Propelled Particles. Sci. China Chem. 2017, 60, 1293–1304. [Google Scholar] [CrossRef]
  32. Wang, W.; Duan, W.; Ahmed, S.; Sen, A.; Mallouk, T.E. From One to Many: Dynamic Assembly and Collective Behavior of Self-Propelled Colloidal Motors. Acc. Chem. Res. 2015, 48, 1938–1946. [Google Scholar] [CrossRef]
  33. Chen, Y.; Chen, D.; Liang, S.; Dai, Y.; Bai, X.; Song, B.; Zhang, D.; Chen, H.; Feng, L. Recent Advances in Field-Controlled Micro–Nano Manipulations and Micro–Nano Robots. Adv. Intell. Syst. 2022, 4, 2100116. [Google Scholar] [CrossRef]
  34. Shields, C.W.; Velev, O.D. The Evolution of Active Particles: Toward Externally Powered Self-Propelling and Self-Reconfiguring Particle Systems. J. Chem. Prac. 2017, 3, 539–559. [Google Scholar] [CrossRef] [Green Version]
  35. Wang, Q.; Zhang, L. External Power-Driven Microrobotic Swarm: From Fundamental Understanding to Imaging-Guided Delivery. ACS Nano 2021, 15, 149–174. [Google Scholar] [CrossRef] [PubMed]
  36. Cates, M.E.; Tailleur, J. Motility-Induced Phase Separation. Annu. Rev. Condens. Matter Phys. 2015, 6, 219–244. [Google Scholar] [CrossRef] [Green Version]
  37. Fu, Y.; Yu, H.; Zhang, X.; Malgaretti, P.; Kishore, V.; Wang, W. Microscopic Swarms: From Active Matter Physics to Biomedical and Environmental Applications. Micromachines 2022, 13, 295. [Google Scholar] [CrossRef] [PubMed]
  38. Theurkauff, I.; Cottin-Bizonne, C.; Palacci, J.; Ybert, C.; Bocquet, L. Dynamic Clustering in Active Colloidal Suspensions with Chemical Signaling. Phys. Rev. Lett. 2012, 108, 268303. [Google Scholar] [CrossRef] [PubMed]
  39. Aziz, A.; Pane, S.; Iacovacci, V.; Koukourakis, N.; Czarske, J.; Menciassi, A.; Medina-Sánchez, M.; Schmidt, O.G. Medical Imaging of Microrobots: Toward In Vivo Applications. ACS Nano 2020, 14, 10865–10893. [Google Scholar] [CrossRef]
  40. Chen, H.; Zhang, H.; Xu, T.; Yu, J. An Overview of Micronanoswarms for Biomedical Applications. ACS Nano 2021, 15, 15625–15644. [Google Scholar] [CrossRef]
  41. Nelson, B.J.; Kaliakatsos, I.K.; Abbott, J.J. Microrobots for Minimally Invasive Medicine. Annu. Rev. Biomed. Eng. 2010, 12, 55–85. [Google Scholar] [CrossRef] [Green Version]
  42. Buttinoni, I.; Bialké, J.; Kümmel, F.; Löwen, H.; Bechinger, C.; Speck, T. Dynamical Clustering and Phase Separation in Suspensions of Self-Propelled Colloidal Particles. Phys. Rev. Lett. 2013, 110, 238301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wittkowski, R.; Tiribocchi, A.; Stenhammar, J.; Allen, R.J.; Marenduzzo, D.; Cates, M.E. Scalar Φ4 Field Theory for Active-Particle Phase Separation. Nat. Commun. 2014, 5, 4351. [Google Scholar] [CrossRef] [Green Version]
  44. Palacci, J.; Sacanna, S.; Steinberg, A.P.; Pine, D.J.; Chaikin, P.M. Living Crystals of Light-Activated Colloidal Surfers. Science 2013, 339, 936–940. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, B.; Sokolov, A.; Snezhko, A. Reconfigurable Emergent Patterns in Active Chiral Fluids. Nat. Commun. 2020, 11, 4401. [Google Scholar] [CrossRef] [PubMed]
  46. Vásárhelyi, G.; Virágh, C.; Somorjai, G.; Nepusz, T.; Eiben, A.E.; Vicsek, T. Optimized Flocking of Autonomous Drones in Confined Environments. Sci. Robot. 2018, 3, eaat3536. [Google Scholar] [CrossRef] [PubMed]
  47. Arnold, R.; Carey, K.; Abruzzo, B.; Korpela, C. What Is A Robot Swarm: A Definition for Swarming Robotics. In Proceedings of the 2019 IEEE 10th Annual Ubiquitous Computing, Electronics & Mobile Communication Conference (UEMCON), New York, NY, USA, 10–12 October 2019; pp. 0074–0081. [Google Scholar]
  48. Rubenstein, M.; Cornejo, A.; Nagpal, R. Programmable Self-Assembly in a Thousand-Robot Swarm. Science 2014, 345, 795–799. [Google Scholar] [CrossRef] [PubMed]
  49. Ballerini, M.; Cabibbo, N.; Candelier, R.; Cavagna, A.; Cisbani, E.; Giardina, I.; Lecomte, V.; Orlandi, A.; Parisi, G.; Procaccini, A.; et al. Interaction Ruling Animal Collective Behavior Depends on Topological Rather than Metric Distance: Evidence from a Field Study. Proc. Natl. Acad. Sci. USA 2008, 105, 1232–1237. [Google Scholar] [CrossRef]
  50. Wallace, D.C. Statistical Mechanical Theory of Liquid Entropy. Int. J. Quantum Chem. 1994, 52, 425–435. [Google Scholar] [CrossRef] [Green Version]
  51. Piaggi, P.M.; Parrinello, M. Entropy Based Fingerprint for Local Crystalline Order. J. Chem. Phys. 2017, 147, 114112. [Google Scholar] [CrossRef]
  52. Wallace, D.C. On the Role of Density Fluctuations in the Entropy of a Fluid. J. Chem. Phys. 1987, 87, 2282–2284. [Google Scholar] [CrossRef]
  53. Floreano, D.; Wood, R.J. Science, Technology and the Future of Small Autonomous Drones. Nature 2015, 521, 460–466. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Brambilla, M.; Ferrante, E.; Birattari, M.; Dorigo, M. Swarm Robotics: A Review from the Swarm Engineering Perspective. Swarm Intell. 2013, 7, 1–41. [Google Scholar] [CrossRef] [Green Version]
  55. Midtvedt, B.; Helgadottir, S.; Argun, A.; Pineda, J.; Midtvedt, D.; Volpe, G. Quantitative Digital Microscopy with Deep Learning. Appl. Phys. Rev. 2021, 8, 011310. [Google Scholar] [CrossRef]
Figure 1. The analysis of one representative simulation run of the boid system (S = 0.1, topological neighbors), using entropy by neighbor distance H NDist and magnitude of the averaged velocity v a (a) representative patterns at several key timesteps: arrows represent the velocity directions of boids. (b) the evolution of v a , (c) the evolution of H NDist .
Figure 1. The analysis of one representative simulation run of the boid system (S = 0.1, topological neighbors), using entropy by neighbor distance H NDist and magnitude of the averaged velocity v a (a) representative patterns at several key timesteps: arrows represent the velocity directions of boids. (b) the evolution of v a , (c) the evolution of H NDist .
Micromachines 14 01503 g001
Figure 2. The entropy by neighbor distances H NDist and the magnitude of the averaged velocities v a of final steady-state patterns in the simulated boids systems as a function of the noise factor η . The steering factor S is 0.001 in all cases. The error bars indicate the standard deviations obtained in at least three independent runs. (a) Neighbors are defined by topology. (b) Neighbors are defined by metric distance.
Figure 2. The entropy by neighbor distances H NDist and the magnitude of the averaged velocities v a of final steady-state patterns in the simulated boids systems as a function of the noise factor η . The steering factor S is 0.001 in all cases. The error bars indicate the standard deviations obtained in at least three independent runs. (a) Neighbors are defined by topology. (b) Neighbors are defined by metric distance.
Micromachines 14 01503 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fu, Y.; Wu, Z.; Zhan, S.; Yang, J.; Gardi, G.; Kishore, V.; Malgaretti, P.; Wang, W. Entropy by Neighbor Distance as a New Measure for Characterizing Spatiotemporal Orders in Microscopic Collective Systems. Micromachines 2023, 14, 1503. https://doi.org/10.3390/mi14081503

AMA Style

Fu Y, Wu Z, Zhan S, Yang J, Gardi G, Kishore V, Malgaretti P, Wang W. Entropy by Neighbor Distance as a New Measure for Characterizing Spatiotemporal Orders in Microscopic Collective Systems. Micromachines. 2023; 14(8):1503. https://doi.org/10.3390/mi14081503

Chicago/Turabian Style

Fu, Yulei, Zongyuan Wu, Sirui Zhan, Jiacheng Yang, Gaurav Gardi, Vimal Kishore, Paolo Malgaretti, and Wendong Wang. 2023. "Entropy by Neighbor Distance as a New Measure for Characterizing Spatiotemporal Orders in Microscopic Collective Systems" Micromachines 14, no. 8: 1503. https://doi.org/10.3390/mi14081503

APA Style

Fu, Y., Wu, Z., Zhan, S., Yang, J., Gardi, G., Kishore, V., Malgaretti, P., & Wang, W. (2023). Entropy by Neighbor Distance as a New Measure for Characterizing Spatiotemporal Orders in Microscopic Collective Systems. Micromachines, 14(8), 1503. https://doi.org/10.3390/mi14081503

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop