Alignment-Induced Self-Organization of Autonomously Steering Microswimmers: Turbulence, Vortices, and Jets
Abstract
Systems of motile microorganisms exhibit a multitude of collective phenomena, including motility-induced phase separation and turbulence. Sensing of the environment and adaptation of movement plays an essential role in the emergent behavior. We study the collective motion of wet self-steering polar microswimmers, which align their propulsion direction hydrodynamically with that of their neighbors, by mesoscale hydrodynamics simulations. The simulations of the employed squirmer model reveal a distinct dependence on the swimmer flow field, i.e., pullers versus pushers. The collective motion of pushers is characterized by active turbulence, with nearly homogeneous density and a Gaussian velocity distribution. Pullers exhibit a strong tendency for clustering and display velocity and vorticity distributions with fat exponential tails; their dynamics is chaotic, with a temporal appearance of vortex rings and fluid jets. Our results show that the collective behavior of intelligent microswimmers is very diverse and still offers many surprises to be discovered.
I Introduction
Emergence of dynamic structures and patterns is an essential feature of biological active motile systems. Examples include microbial swarms [1, 2] on the cellular level, to schools of fish [3, 4], flocks of birds [5, 6], and collective motion in human crowds [7] on a macroscopic level. Also, in artificial active systems consisting of synthetic self-propelled particles and microrobots, the collective dynamics of constituent objects is of fundamental importance for their application in engineering and medicine to achieve a large spectrum of functionalities [8, 9, 10, 11].
A fundamental aspect in such systems is the active and autonomous motion of the constituting particles. While activity and self-propulsion can give rise to several novel types of collective behaviors, such as motility-induced phase separation (MIPS) [12, 13] and active turbulence [14, 15, 16], the fact that biological microswimmers are not only motile, but also gather information about their environment and adapt their motion by self-steering, remains largely unexplored and yet to be elucidated [17, 18, 19].
Many living organisms are immersed in a fluid medium, and their collective behavior is strongly affected or even dominated by hydrodynamics [20, 21]. The hydrodynamic environment is not just the background medium in which aquatic microorganisms are based, but it is rather essential for locomotion on the individual level as well as inter-organism interactions [1]. On the mesoscale, life has adapted to low-Reynolds hydrodynamics, e.g., in the emergence of bacterial turbulence [22, 23, 16] and in coordinated cell migration during embryogenesis [24, 25]; similar physical laws govern the swarming microrobots [8, 26]. However, studying large-scale hydrodynamic systems is challenging as the fluid adds a large number of degrees of freedom, while the consideration of large-scale systems is unavoidable to accurately capture long-range hydrodynamics interactions and to minimize potential finite-size effects. So far, only limited studies have been conducted in this direction, particularly in three dimensions.
The goal of our current endeavor is to unravel the emergent collective behavior of systems, which combine two essential components of living and artificial active systems, self-steering and hydrodynamics. In the context of active matter, the Vicsek model [27, 28] is among the earliest and simplest models for the collective directional motion of self-propulsion and self-steering particles due to mutual alignment. While the role of alignment interactions has been extensively investigated for dry active systems [28], hydrodynamic interactions themselves have rather been viewed as a physical alignment mechanism in wet systems [29, 30]. However, hydrodynamic propulsion can in fact also destabilize polar order [31, 32, 33], which implies the necessity for an additional stabilization mechanism for alignment. To achieve stable alignment, we consider a hydrodynamic extension of the Vicsek model. Our active agents, modeled as squirmers [34, 35, 36, 37, 38], sense the propulsion direction of neighboring agents and adapt their propulsion direction accordingly by hydrodynamic self-steering [39, 40], with a “slow” temporal response due to limited maneuverability. We perform large-scale simulations in three-dimensional systems, capturing the fluid environment by the multiparticle collision dynamics (MPC) technique [41, 42, 43], a particle-resolved mesoscale hydrodynamic simulation approach.
We observe and characterize the emergence of swarming dynamics in polar active fluids, consisting either of pusher or puller microswimmers. Our results show that hydrodynamic interactions destabilize polar order, giving rise to rich collective spatio-temporal behavior beyond the simple symmetry breaking of the dry Vicsek model. Pusher systems feature active turbulence with non-universal scaling exponents in the kinetic energy spectrum, revealing a route toward active turbulence via self-steering. Systems of pullers are accompanied by formation of dense, swarming clusters driven by hydrodynamic interactions. In particular, formation of toroidal structures are observed in the vorticity field, which are characterized by enhanced spatial vorticity-velocity cross correlations at short distances. This demonstrates that the formation of vortex rings is a direct consequence of strong active jets caused by propulsion and alignment.
II Result
II.1 Hydrodynamic Vicsek model and Polar Order
We consider a system of spherical active microswimmers with radius and instantaneous orientation , . For self-propulsion and self-steering, the squirmer model is employed, where the self-propulsion is achieved via an axisymmetric surface-slip boundary condition, characterized by the speed and the active stress [see Materials and Methods, Eq. (4)]. Self-steering of microswimmers is modeled via an adaptive non-axisymmetric surface-slip boundary condition [see Materials and Methods, Eqs. (4) and (5)), and Fig. 1(b),(c)], which enables the squirmer to rotate and consequently reorient toward a desired direction with the limited angular velocity [44]
(1) |
where characterizes the strength of adaptation. Accordingly, we introduce two dimensionless parameters, the Péclet number and the maneuverability in the form
(2) |
where and are the diameter and the (thermal) rotational diffusion coefficient of a squirmer, respectively.
In Eq. (1), sensed information of each particle about the orientation of other neighboring particles is represented by a vector . In the interest of the investigation and understanding of the generic collective behavior of an active polar fluid, we focus here on a “minimal” model, where only the swimming and steering mechanism via the adaptive surface flow field is explicitly considered, whereas sensing and information processing is taken into account implicitly by the sensed-information vector in Eq. (1), see Ref. [40] for a more detailed discussion. As a representative example of information exchanges between intelligent microswimmers, we consider a Vicsek-type alignment interaction, where each microswimmer aims at adapting its orientation and propulsion direction to the average orientation of neighboring particles, see Fig. 1(a) for illustration. Specifically, we employ a non-additive rule of orientation adaptation, which results in non-reciprocal interactions between microswimmers, where [45]
(3) |
Here, is the set of neighbors of the -th particle in its alignment range , and is the number of neighbors. As apparent from the definition in Eq. (3), , which serves as an input signal triggering adaptive surface flows according to Eqs. (4) and (5), is typically not a unit vector. Therefore, the magnitude of the adaptation force depends on the strength of orientational order of the neighboring microswimmers. We emphasize that in our approach the steering is achieved solely via the modification of the surface flow fields, mimicking the autonomous behaviors of microorganisms, in contrast to external driving forces, see, e.g., Ref. [46]. This difference leads to very different behaviors in wet and dry systems.
For the MPC fluid dynamics simulations, we consider a MPC variant with angular momentum conservation [47], see Supplementary Materials (SM), Sec. S-I for more details. Our highly parallelized, GPU-accelerated implementation employed for the simulations is based on the framework proposed in Ref. [48]. For an accurate characterization of emergent behaviors, we consider large system sizes up to , where is the length of the cubic simulation box and is the side length of a MPC collision cell, and up to squirmers. For the squirmers, we choose a sensing range , and a strength of the active stress and for pushers and pullers, respectively. For most simulations, we consider the packing fraction (based on squirmer radius), or (based on sensing range), if not explicitly stated otherwise. For comparison, we also perform simulations of a dry system of aligning self-steering “intelligent” active Brownian particles (iABPs) of the same packing fraction.
Results for the global polar order parameter are shown in Fig. 1(d). While the systems are disordered for small maneuverability in both dry and wet cases, a global polar order emerges only in dry systems for large , indicating that hydrodynamic interactions destabilize the polar order. Instead of simple polar ordering, our squirmer systems show swarming dynamics with and without density modulation depending on the swimming mechanism (pusher or puller).
II.2 Pushers: Active Turbulence via Self-Steering
Systems of rearly propelled squirmers (pushers) exhibit active turbulence, see Fig. 2(a)-I for a snapshot (also Movie S1), and their collective motion features vortical flows. Figure 2(a)-II displays the fluid velocity field reflecting vortical structures and fluctuations in the magnitude of the vorticity. Accordingly, equal-time spatial velocity correlations become negative at large distances (less than ) for large , see Fig. 2(b). However, no pronounced density fluctuations are visible, see Fig. 2(a)-I. This is confirmed by the analysis of the local density distribution via from Voronoi tessellation [49]. As shown in Fig. 2(c), the distribution exhibits a peak near the global squirmer density independent of .
To characterize the dynamics, we examine the kinetic energy spectrum as a representative indicator of active turbulence [15]. We determine the energy spectra for both the squirmer and the fluid. For the squirmers, we first calculate spatial velocity correlation functions and then perform a Fourier transform to obtain the energy spectrum. For the fluid, we calculate the energy spectrum directly from the (Eulerian) velocity field.
Before proceeding to a more detailed analysis of systems of self-steering squirmers, we briefly discuss as reference case the dynamics of pushers without self-steering, i.e., squirmers with . Interestingly, we do not observe any significant collective behavior without self-steering for as well as , with and , in our simulations (see SM, Fig. S1). It seems that thermal fluctuations and consequently rotational diffusion of the microswimmers suppress any structure formation at the considered small packing fraction, compare also Ref. [50] for a related lattice Boltzmann simulation study of pushers. Moreover, the energy spectrum of squirmers strongly deviates from the fluid energy spectrum for large and intermediate (Fig. S1), indicating that the fluid flows only mildly affect the squirmer dynamics. As also pointed out in Ref. [23], presumably even higher densities of force-dipoles are necessary to promote the emergence of collective motion engaging both the fluid and microswimmers of pushers without self-steering, solely based on hydrodynamic alignment. Indeed, the number density explored in Refs. [51, 52] is typically one or two orders of magnitude larger than in our case due to different geometry of active particles and their use of force dipoles.
In sharp contrast, systems of self-steering squirmers display pronounced self-organization, as demonstrated in Fig. 2. With increasing maneuverability , a scaling regime emerges in the fluid power spectrum for wave numbers , see Fig. 2(d). Corresponding energy spectra of the squirmer motion exhibit the same scaling behavior with the same exponents – with scaling extending even to lower wave numbers for – confirming the emergence of active turbulence, where the dynamics of squirmers and fluid are strongly correlated on larger length scales. The exponents of the energy spectrum, are found to be non-universal, depending on and increasing with , roughly in the range for , see Fig. S2. Moreover, with increase , not only increases, but also the scaling regimes extends to large length scales (smaller , as more squirmers participate in the self-organized vortex structures, and the peak height in the energy spectrum increases (up to for ). This indicates that squirmers attain much higher velocities, which is quantitatively confirmed by the mean-square displacement (MSD), see Fig. 3. In the ballistic regime, the MSD , hence, the increasing amplitude in Fig. 3 is a direct measure of the increasing speed as increases for , i.e., in the regime where a broad scaling region can be identified in the energy spectrum (Fig. 2(d)). The Péclet number merely affects turbulent dynamics on large length scales, but the ratio determines the scaling exponent as shown in Fig. S2. We also probe density effects varying the number of squirmers. We first notice that the fluid energy spectra seem to converge for high densities, as shown in Fig. 2(e), in accordance with the previous report [52]. As the packing fraction decreases from to corresponding to and , downward shifts in the fluid energy spectra are observed, indicating that fluid stirring by squirmers are not strong at low densities. Consequently, the scaling regime in the energy spectra of squirmers shrinks as the density decreases, e.g., to a narrow range of at .
For , or length scales smaller than about , differences between the fluid and squirmer energy spectra appear. In this regime, near-field hydrodynamic flows play a significant role. Moreover, the properties of the energy spectra depend on , indicating that noise effects are significant at these small length scales, see Fig. S3 for fluid thermal energy spectra [53]. Therefore, in wet systems and small densities, a strong self-steering of active particles is crucial for the emergence of a large-scale coherent collective motion. Otherwise disorder and hydrodynamic instabilities on small length-scales may prevail [31]. Such an observation should also apply to systems where the alignment is mediated via steric repulsion of elongated body shapes [15] or strong hydrodynamic force dipoles [52]. In any cases, in active turbulence of wet systems, collective fluid flows induced by microswimmers build up the collective behavior of microswimmers on large length-scales with fast dynamics. In sharp contrast, active turbulence in dry systems [54, 55] requires densely packed active particles, because a speedup mechanism as in wet systems is lacking as steric repulsions may only result in a slow-down of active particles. Instead, in dry active turbulence, the scaling regime develops in large , or equivalently, small regimes due to chaotic interparticle collisions on small length scales.
II.3 Pullers: Swarming Dynamics via Self-Steering
In a systems of self-steering pullers, a rich swarming dynamics develops, as shown in Fig. 4. The self-organization is characterized by the formation of morphologically complex clusters of microswimmers, which, on larger length scales, exhibit visually chaotic movements and exchange constitutive squirmers with each others, see Fig. 4(a)-I (also Movie S2). Still, the puller system exhibits a velocity field with vortical structures (see Fig. 4(a)-II), surprisingly similar to pushers.
Density Modulation.
The local density distribution is calculated by employing a Voronoi tessellation [49]. In sharp contrast to pusher systems (Fig. 2(c)), puller systems exhibit significant density modulations, see Fig. 4(b). For large , , the density distribution is broad, while for , a clear tendency of segregation is observed (Movie S3) with a low-density peak at , and a high-density peak at . For small in the range , non-mobile clusters appear (see Fig. S4). We focus here on the dynamic clusters, as formation of a dense static cluster may be affected by a depletion of MPC fluid particles inside the cluster – which is related to the (weak) compressibility of the MPC fluid, and artificially enhance cluster stability [56]. For , unimodal distributions are recovered, but with a peak at a density smaller than the global density of , and with fatter tails than those in pusher systems, which reflects a clustering tendency reported for pullers [56, 57].
Decoupling of Puller Orientation and Velocity.
The emergence of a high-density regime for needs to be distinguished from motility-induced phase separation (MIPS) in dry systems. First, squirmers in a (small) cluster exhibit local polar order, giving rise to a coherent directional motion of the cluster. Second, the pullers actually swim faster on average than their bare self-propulsion speed , as shown in Fig. 4(c) – as in the system of pushers (see Fig. 2(e)). In particular, decoupling of self-propulsion and velocity of squirmers is far more significant than in MIPS. Even situations, where pullers are driven backward occur frequently (Movie S4). Indeed, as shown in Fig. 4(d), the probability that the velocity of a squirmer is anti-parallel to its orientation, i.e., , is even higher than for the parallel case, i.e., . Therefore, we conclude that hydrodynamic interactions between self-steering pullers dominate over the self-propulsion forces. The attractive hydrodynamic interactions between aligned pullers in a head-to-tail configuration promotes the formation of dense clusters. Yet, as we will demonstrate below, dense clusters are not static, but exhibit a highly dynamical morphology.
Decoupling of Puller and Fluid Dynamics.
The kinetic energy spectra for squirmer motion and the fluid in puller systems are presented in Fig. 3. While for the swarming of pullers, exhibits a power-law decay in the intermediate wave-number regime () for the fluid, as in pusher systems, several features significantly deviate from those of pushers, which again demonstrates the uniqueness of puller swarming. Most of all, a pronounced mismatch between squirmer and fluid energy spectra is observed in the intermediate regime of (see Fig. S5(a) and Fig. S7(a) for the corresponding correlation functions), which indicates that pullers are not simply driven by the fluid flow. The exponents of the three scaling regimes in Fig. 4(e) show universal behavior for the investigated range of . Notably, in the case of the squirmer energy spectra for , which correspond the length scales of the high-density intra-cluster regime where steric repulsion and near-field hydrodynamic interactions dominate, the obtained value of the exponent is in accordance with the Kolmogorov exponent for classical hydrodynamic turbulence and with active turbulence of dense quasi-two-dimensional pusher systems [15]. For larger length scales, an exponent is obtained, which persists under variation of global densities, see Fig. S6. A similar value has been found previously in Lattice Boltzmann simulations of extended force dipoles [52], although at higher densities and for pushers.
Formation of vortex ring.
An even more striking feature is observed in configurations of the vorticity field. Visual observation of the time evolution of the system indicates a typical dynamical behavior in the morphology of clusters, which involves the pulsatile transformation of aggregates from a spherical shape into jellyfish-like arrangements. In terms of fluid mechanics, the jellyfish-like morphology suggests formation of a vortex ring, which is indeed confirmed by emergence of toroidal structures in the vorticity fields extracted from the simulations, see Fig. 5(a). As shown in Fig. 5(b), the vorticity field exhibits whirling patterns within regions where the magnitude of the vorticity field is large.
For a detailed illustration, we consider the small system size , where only a single cluster emerges (Movie S5). As shown in Fig. 5(c)-I, the dynamics initiates with formation of a cluster. Then, due to alignment, squirmers rotate, and a polar order emerges within the cluster, see Fig. 5(c)-II, and also Fig. S5(b) for the corresponding order parameter. Such an ordered structure gives rise to a strong collective fluid flow, which generates a pronounced jet in front of the cluster, see the yellow surface in Fig. 5(e)-I. Notably, the jet flow is self-generated via active stirring of microswimmers in this case, instead of external perturbation as in passive hydrodynamic fluids. Subsequently, a spread-out motion of squirmers is initiated (Fig. 5(c)-III). Simultaneously, a vortex ring is formed around the cluster (blue ring in Fig. 5(e)-I), while squirmers are moving forward. As shown in Fig. 5, the velocity field is indeed wrapping around the vortex ring, in accordance with Fig. 5(b). Then, the pullers continue to spread out, rolling about the region where the vortex ring forms, see Fig. 5(c)-IV. While swimmers at the cluster center swim forward, they are dragged backward at the periphery, as shown in Fig. 5(e)-II, which contributes to the anomalous behavior in the distribution of (see Fig. 4(d)). Eventually, the cluster dissolves, ending its life cycle.
II.4 Velocity-Vorticity Coupling
The sequential time evolution described so far indicates a strong coupling between the fluid velocity and the vorticity field in a puller cluster. As shown in Fig. 6(a), the rotation of the vorticity field for pullers is indeed centered at regions with a strong velocity field.
For a more quantitative characterization, we examine the distribution of the Cartesian components of the velocity and vorticity fields. In pusher systems, both the velocity and vorticity fields (Fig. S7(b) and Fig. 6(b), respectively) exhibit a Gaussian distribution, which is an indicator of active turbulence [15]. For pullers, both deviate from a Gaussian, demonstrating that the swarming dynamics of pullers is not active turbulence. Specifically, the velocity distribution of pullers exhibits “fat” exponential tails, as shown in Fig. S7(b), in line with the emergence of stronger fluid flows than “expected”, i.e., the occurrence of jet plumes induced by aligned pullers. Also the vorticity field for pullers shows a broader distribution than that of pushers as displayed in Fig. 6(b). Furthermore, the vorticity distribution for pullers in Fig. 6(b) shows a rather sharp peak at , the average of three Cartesian components () of , indicating a weak separation between regions with strong and weak vorticity.
A more fundamental difference in velocity-vorticity coupling is revealed by a cross correlations between the magnitude of vorticity and that of the velocity field ( and ), as defined in Eq. (15). In Fig. 6(c), the various cross correlations for pushers exhibit a pronounced peak at , , and for , , and , respectively, while the cross correlation at is not strong. Hence, for pushers, the velocity field of a vortex is weak at the center but strong at intermediate regions between the center and periphery. In sharp contrast, for pullers, the cross correlation between vorticity and velocity fields is already strong at small distances, which demonstrates that a strong velocity field generates a strong vorticity in the immediate vicinity of the jet flow. Moreover, the cross correlation decays faster than for pushers, assuming negative values before approaching zero.
III Discussion and Conclusions
We have studied self-organization and dynamics in three-dimensional wet systems of self-steering squirmers, which aim for alignment of their orientation with that of their neighbors. We demonstrate that alignment via hydrodynamic self-steering gives rise to a rich collective behavior in such polar active fluids, depending on the type of active stresses, i.e., whether the microswimmers are pushers or pullers. In both cases, an essential role of hydrodynamics is the breaking of long-range polar order, which causes the emergence of chaotic behavior.
For pushers, the particle distribution is quite homogeneous, the distribution of the Cartesian velocity components is Gaussian, and the kinetic energy spectrum displays a peak and subsequent power-law decay with increasing wave vector, which indicates active-turbulent behavior. An intriguing feature is that strong self-steering enhances the coherent movement of microswimmers, which leads to collective particle motion with speeds much faster then the individual swim speed – as can be seen in the increasing magnitude of the peak in the energy spectrum, combined with an extension of the scaling regime toward large length scales. This implies that large-scale flows are induced by the collective motion, which drag the microswimmers along and supersede their individual motion. Thus, the polarity field and the fluid flow field are strongly coupled.
For pullers, another type of self-organization emerges, which is strictly distinguished from motility-induced phase separation (MIPS) of dry ABP systems as well as active turbulence of pushers. The particle density is now found to be very inhomogeneous, as the pullers tend to form clusters. However, these clusters are not static, but appear as quite unstable. The particle alignment inside the cluster generates a strong fluid jet and a vortex ring, which pulls apart the cluster and leads to its disintegration. These strong flows imply that the probability of fast fluid flows is enhanced, which is reflected in the emergence of fat tails in the velocity distribution.
In contrast, wet systems of self-propelled squirmers without self-steering display no interesting collective behavior in three dimensions, not even at high squirmer volume fractions.
Our numerical observations for ensembles of self-steering pullers challenges current theoretical views on collective behaviors in wet active systems. So far, it has been typically assumed that the polarity and velocity fields of active fluids are essentially identical, based on the assumption of a nearly homogeneous distribution of active particles. Heterogeneous densities have been observed recently in models of compressible polar active fluids for bacterial suspensions [58, 59]; however, the mechanism is entirely different in this case, as hydrodynamic interactions are not considered, and clustering is driven by a strong dependence of self-propulsion speed on the local density.
Our results demonstrate that in wet systems of self-steering microswimmers in three dimensions, the interplay of the particle density and polarity, and of the fluid velocity field can give rise to a surprisingly rich variety of emergent behaviors – already for a highly simplified model systems with only a single particle type.
IV Materials and Methods
IV.1 Mesoscale Fluid Model: Multi-Particle Collision Dynamics
We adopt the multiparticle collision dynamics (MPC) method [42, 43], a particle-based mesoscale simulation approach, as model for the fluid. Specifically, we employ the stochastic rotation variant of MPC with angular momentum conservation (MPC-SRD+a) [47] and the cell-level Maxwell-Boltzmann scaling thermostat [60]. The algorithm consists of alternating streaming and collision steps. In the streaming step, the MPC point particles of mass propagate ballistically over the collision time interval , denoted as collision time. In the collision step, fluid particles are sorted into the cells of a cubic lattice of lattice constant defining the collision environment. Then, their relative velocities, with respect to the center-of-mass velocity of the collision cell, are rotated around a randomly oriented axes by a fixed angle . The algorithm conserves mass, linear, and angular momentum on the collision-cell level, while thermal fluctuations are automatically incorporated. More details are described in Refs. [38, 15] and the SM, which refers to Refs. [61, 62, 63, 64] additionally. Our GPU-based, highly parallelized implementation employed for the simulations is proposed in Ref. [48].
IV.2 Hydrodynamic Self-Steering
Hydrodynamic self-steering of squirmers is achieved via adaptive surface flow fields, given as
(4) | ||||
(5) |
where and are the polar and azimuthal angles in a body-fixed reference frame. The parameter characterizes the swim speed and the strength of the force dipole, where for pullers and for pushers [21, 66]. and control the magnitudes of non-axisymmetric surface-flow components, leading to rotational motion of the body [39, 40]. Specifically, the non-axisymmetric flow fields enable the adaptive motion of Eq. (1) via
(6) | ||||
(7) |
where and are the unit vectors along the axis of the body-fixed reference frame. For the self-propulsion, we use and , which correspond to and , and and , respectively. The values of the self-steering strength are varied from to yielding .
IV.3 Steric squirmer interaction
Steric repulsion between two squirmers is described by the separation-shifted Lennard-Jones potential
(8) |
for and zero otherwise, where indicates the surface-to-surface distance between the two squirmer. To avoid loss of hydrodynamic interactions when two squirmers contact with each other, we also include a virtual safety distance [67, 15], which leads to the effective distance , where denotes the center-to-center distance and is the squirmer diameter. We choose . Numerically, the equations of motion for the rigid-body dynamics of the squirmers are solved by the velocity-Verlet algorithm, see SM for more details.
IV.4 Spatial correlation and Energy Spectrum
Squirmers: Particle-based approach.
The spatial velocity correlation function of the squirmers is defined as
(9) |
where . Here, we use the velocity averaged over a short time interval , for which we consider . The energy spectrum can be calculated via Fourier transformation. Here we consider the Fourier sine transform [68]
(10) |
Fluid: Field-based approach.
We first extract the fluid velocity field from the simulation data by introducing a grid dividing the whole system into cells. The velocities of all MPC particles, averaged over a short time interval as for squirmers, are additionally averaged over each cell to obtain with for and . Then the discrete Fourier transform
(11) |
is performed. The energy spectrum is calculated straightforwardly via [69]
(12) |
which is then averaged over all directions of to obtain . Then, the spatial velocity correlation function is obtained via the Fourier transformation
(13) |
from which we calculate by averaging over all directions of . To reduce noise effects at small length scales, velocity fields are averaged over boxes with side lengths , , or .
The vorticity field is then defined as
(14) |
Numerically, the vorticity field is calculated by the five-point stencil method from the velocity field. The vorticity spatial correlation are also obtained from Eqs. (11)-(13). Moreover, a gliding time average is performed for the velocity with a time window , corresponding to .
Cross correlation.
We again utilize Fourier transformation to calculate cross correlations between velocity and vorticity fields. Specifically, we first calculate the magnitudes of velocity and vorticity fields, which are then shifted by their average values, i.e., and . Then, from the Fourier transforms of the fields and , the cross correlation is obtained via
(15) |
where superscript indicates the complex conjugate. is obtained by averaging over all directions.
References
- Pedley and Kessler [1992] T. J. Pedley and J. O. Kessler, Hydrodynamic phenomena in suspensions of swimming microorganisms, Annu. Rev. Fluid Mech. 24, 313 (1992).
- Be’er and Ariel [2019] A. Be’er and G. Ariel, A statistical physics view of swarming bacteria, Mov. Ecol. 7, 9 (2019).
- Katz et al. [2011] Y. Katz, K. Tunstrøm, C. C. Ioannou, C. Huepe, and I. D. Couzin, Inferring the structure and dynamics of interactions in schooling fish, Proc. Natl. Acad. Sci. USA 108, 18720 (2011).
- Gautrais et al. [2012] J. Gautrais, F. Ginelli, R. Fournier, S. Blanco, M. Soria, H. Chaté, and G. Theraulaz, Deciphering interactions in moving animal groups, PLOS Comput. Biol. 8, 1 (2012).
- Ballerini et al. [2008] M. Ballerini, N. Cabibbo, R. Candelier, A. Cavagna, E. Cisbani, I. Giardina, A. Orlandi, G. Parisi, A. Procaccini, M. Viale, and V. Zdravkovic, Empirical investigation of starling flocks: a benchmark study in collective animal behaviour, Anim. Behav. 76, 201 (2008).
- Bialek et al. [2012] W. Bialek, A. Cavagna, I. Giardina, T. Mora, E. Silvestri, M. Viale, and A. M. Walczak, Statistical mechanics for natural flocks of birds, Proc. Natl. Acad. Sci. USA 109, 4786 (2012).
- Helbing et al. [2000] D. Helbing, I. Farkas, and T. Vicsek, Simulating dynamical features of escape panic, Nature , 487 (2000).
- Xie et al. [2019] H. Xie, M. Sun, X. Fan, Z. Lin, W. Chen, L. Wang, L. Dong, and Q. He, Reconfigurable magnetic microrobot swarm: Multimode transformation, locomotion, and manipulation, Sci. Robot. 4, eaav8006 (2019).
- Wang et al. [2021] Q. Wang, K. F. Chan, K. Schweizer, X. Du, D. Jin, S. C. H. Yu, B. J. Nelson, and L. Zhang, Ultrasound doppler-guided real-time navigation of a magnetic microswarm for active endovascular delivery, Sci. Adv. 7, eabe5914 (2021).
- Chen and Bechinger [2022] C.-J. Chen and C. Bechinger, Collective response of microrobotic swarms to external threats, N. J. Phys. 24, 033001 (2022).
- Wang et al. [2023] X. Wang, P.-C. Chen, K. Kroy, V. Holubec, and F. Cichos, Spontaneous vortex formation by microswimmers with retarded attractions, Nat. Commun. 14, 56 (2023).
- Fily and Marchetti [2012] Y. Fily and M. C. Marchetti, Athermal phase separation of self-propelled particles with no alignment, Phys. Rev. Lett. 108, 235702 (2012).
- Buttinoni et al. [2013] I. Buttinoni, J. Bialké, F. Kümmel, H. Löwen, C. Bechinger, and T. Speck, Dynamical clustering and phase separation in suspensions of self-propelled colloidal particles, Phys. Rev. Lett. 110, 238301 (2013).
- Wensink et al. [2012] H. H. Wensink, J. Dunkel, S. Heidenreich, K. Drescher, R. E. Goldstein, H. Löwen, and J. M. Yeomans, Meso-scale turbulence in living fluids, Proc. Natl. Acad. Sci. USA 109, 14308 (2012).
- Qi et al. [2022] K. Qi, E. Westphal, G. Gompper, and R. G. Winkler, Emergence of active turbulence in microswimmer suspensions due to active hydrodynamic stress and volume exclusion, Commun. Phys. 5, 49 (2022).
- Aranson [2022] I. S. Aranson, Bacterial active matter, Rep. Prog. Phys. 85, 076601 (2022).
- Ziepke et al. [2022] A. Ziepke, I. Maryshev, I. S. Aranson, and E. Frey, Multi-scale organization in communicating active matter, Nat. Commun. 13, 6727 (2022).
- Sawicki et al. [2023] J. Sawicki, R. Berner, S. A. M. Loos, M. Anvari, R. Bader, W. Barfuss, N. Botta, N. Brede, I. Franović, D. J. Gauthier, S. Goldt, A. Hajizadeh, P. Hövel, O. Karin, P. Lorenz-Spreen, C. Miehl, J. Mölter, S. Olmi, E. Schöll, A. Seif, P. A. Tass, G. Volpe, S. Yanchuk, and J. Kurths, Perspectives on adaptive dynamical systems, Chaos 33, 071501 (2023).
- Negi et al. [2024] R. S. Negi, R. G. Winkler, and G. Gompper, Collective behavior of self-steering active particles with velocity alignment and visual perception, Phys. Rev. Res. 6, 013118 (2024).
- Koch and Subramanian [2011] D. L. Koch and G. Subramanian, Collective hydrodynamics of swimming microorganisms: Living fluids, Annu. Rev. Fluid Mech. 43, 637 (2011).
- Elgeti et al. [2015] J. Elgeti, R. G. Winkler, and G. Gompper, Physics of microswimmers—single particle motion and collective behavior: a review, Rep. Prog. Phys. 78, 056601 (2015).
- Dombrowski et al. [2004] C. Dombrowski, L. Cisneros, S. Chatkaew, R. E. Goldstein, and J. O. Kessler, Self-concentration and large-scale coherence in bacterial dynamics, Phys. Rev. Lett. 93, 098103 (2004).
- Drescher et al. [2011] K. Drescher, J. Dunkel, L. H. Cisneros, S. Ganguly, and R. E. Goldstein, Fluid dynamics and noise in bacterial cell-cell and cell-surface scattering, Proc. Natl. Acad. Sci. USA 108, 10940 (2011).
- Deneke et al. [2019] V. E. Deneke, A. Puliafito, D. Krueger, A. V. Narla, A. De Simone, L. Primo, M. Vergassola, S. De Renzis, and S. Di Talia, Self-organized nuclear positioning synchronizes the cell cycle in drosophila embryos, Cell 177, 925 (2019).
- Hernández-López et al. [2023] C. Hernández-López, A. Puliafito, Y. Xu, Z. Lu, S. D. Talia, and M. Vergassola, Two-fluid dynamics and micron-thin boundary layers shape cytoplasmic flows in early Drosophila embryos, Proc. Natl. Acad. Sci. USA 120, e2302879120 (2023).
- Ceron et al. [2023] S. Ceron, G. Gardi, K. Petersen, and M. Sitti, Programmable self-organization of heterogeneous microrobot collectives, Proc. Natl. Acad. Sci. USA 120, e2221913120 (2023).
- Vicsek et al. [1995] T. Vicsek, A. Czirók, E. Ben-Jacob, I. Cohen, and O. Shochet, Novel type of phase transition in a system of self-driven particles, Phys. Rev. Lett. 75, 1226 (1995).
- Chaté [2020] H. Chaté, Dry aligning dilute active matter, Annu. Rev. Condens. Matter Phys. 11, 189 (2020).
- Reinken et al. [2018] H. Reinken, S. H. L. Klapp, M. Bär, and S. Heidenreich, Derivation of a hydrodynamic theory for mesoscale dynamics in microswimmer suspensions, Phys. Rev. E 97, 022613 (2018).
- Liu et al. [2021a] Z. T. Liu, Y. Shi, Y. Zhao, H. Chaté, X. qing Shi, and T. H. Zhang, Activity waves and freestanding vortices in populations of subcritical quincke rollers, Proc. Natl. Acad. Sci. USA 118, e2104724118 (2021a).
- Aditi Simha and Ramaswamy [2002] R. Aditi Simha and S. Ramaswamy, Hydrodynamic fluctuations and instabilities in ordered suspensions of self-propelled particles, Phys. Rev. Lett. 89, 058101 (2002).
- Llopis and Pagonabarraga [2010] I. Llopis and I. Pagonabarraga, Hydrodynamic interactions in squirmer motion: Swimming with a neighbour and close to a wall, J. Non-Newtonian Fluid Mech. 165, 946 (2010).
- Clopés et al. [2020] J. Clopés, G. Gompper, and R. G. Winkler, Hydrodynamic interactions in squirmer dumbbells: active stress-induced alignment and locomotion, Soft Matter 16, 10676 (2020).
- Lighthill [1952] M. J. Lighthill, On the squirming motion of nearly spherical deformable bodies through liquids at very small Reynolds numbers, Comm. Pure Appl. Math. 5, 109 (1952).
- Blake [1971] J. R. Blake, A spherical envelope approach to ciliary propulsion, J. Fluid Mech. 46, 199 (1971).
- Ishikawa et al. [2006] T. Ishikawa, M. P. Simmonds, and T. J. Pedley, Hydrodynamic interaction of two swimming model micro-organisms, J. Fluid Mech. 568, 119 (2006).
- Götze and Gompper [2010] I. O. Götze and G. Gompper, Mesoscale simulations of hydrodynamic squirmer interactions, Phys. Rev. E 82, 041921 (2010).
- Theers et al. [2016a] M. Theers, E. Westphal, G. Gompper, and R. G. Winkler, Modeling a spheroidal microswimmer and cooperative swimming in a narrow slit, Soft Matter 12, 7372 (2016a).
- Pak and Lauga [2014] O. S. Pak and E. Lauga, Generalized squirming motion of a sphere, J. Eng. Math. 88, 1 (2014).
- Goh et al. [2023] S. Goh, R. G. Winkler, and G. Gompper, Hydrodynamic pursuit by cognitive self-steering microswimmers, Commun. Phys. 6, 310 (2023).
- Malevanets and Kapral [1999] A. Malevanets and R. Kapral, Mesoscopic model for solvent dynamics, J. Chem. Phys. 110, 8605 (1999).
- Kapral [2008] R. Kapral, Multiparticle collision dynamics: Simulations of complex systems on mesoscale, Adv. Chem. Phys. 140, 89 (2008).
- Gompper et al. [2009] G. Gompper, T. Ihle, D. M. Kroll, and R. G. Winkler, Multi-particle collision dynamics: A particle-based mesoscale simulation approach to the hydrodynamics of complex fluids, Adv. Polym. Sci. 221, 1 (2009).
- Goh et al. [2022] S. Goh, R. G. Winkler, and G. Gompper, Noisy pursuit and pattern formation of self-steering active particles, New J. Phys. 24, 093039 (2022).
- Chepizhko et al. [2021] O. Chepizhko, D. Saintillan, and F. Peruani, Revisiting the emergence of order in active matter, Soft Matter 17, 3113 (2021).
- Bera et al. [2022] A. Bera, S. Sahoo, S. Thakur, and S. K. Das, Active particles in explicit solvent: Dynamics of clustering for alignment interaction, Phys. Rev. E 105, 014606 (2022).
- Noguchi and Gompper [2008] H. Noguchi and G. Gompper, Transport coefficients of off-lattice mesoscale-hydrodynamics simulation techniques, Phys. Rev. E 78, 016706 (2008).
- Westphal et al. [2024] E. Westphal, S. Goh, R. G. Winkler, and G. Gompper, Htmpc: A heavily templated C++ library for large scale particle-based mesoscale hydrodynamics simulations using multiparticle collision dynamics, arXiv preprint arXiv:2406.15236 10.48550/arXiv.2406.15236 (2024).
- Rycroft [2009] C. H. Rycroft, VORO++: A three-dimensional Voronoi cell library in C++, Chaos 19, 041111 (2009).
- Gascó et al. [2023] A. Gascó, I. Pagonabarraga, and A. Scagliarini, Three-dimensional active turbulence in microswimmer suspensions: simulations and modelling, arXiv preprint arXiv:2304.03662 10.48550/arXiv.2304.03662 (2023).
- Stenhammar et al. [2017] J. Stenhammar, C. Nardini, R. W. Nash, D. Marenduzzo, and A. Morozov, Role of correlations in the collective behavior of microswimmer suspensions, Phys. Rev. Lett. 119, 028005 (2017).
- Bárdfalvy et al. [2019] D. Bárdfalvy, H. Nordanger, C. Nardini, A. Morozov, and J. Stenhammar, Particle-resolved lattice Boltzmann simulations of 3-dimensional active turbulence, Soft Matter 15, 7747 (2019).
- Huang et al. [2012] C.-C. Huang, G. Gompper, and R. G. Winkler, Hydrodynamic correlations in multiparticle collision dynamics fluids, Phys. Rev. E 86, 056711 (2012).
- Wensink and Löwen [2012] H. H. Wensink and H. Löwen, Emergent states in dense systems of active rods: from swarming to turbulence, J. Phys.: Condens. Matter 24, 460130 (2012).
- Keta et al. [2024] Y.-E. Keta, J. U. Klamser, R. L. Jack, and L. Berthier, Emerging mesoscale flows and chaotic advection in dense active matter, Phys. Rev. Lett. 132, 218301 (2024).
- Theers et al. [2018] M. Theers, E. Westphal, K. Qi, R. G. Winkler, and G. Gompper, Clustering of microswimmers: interplay of shape and hydrodynamics, Soft Matter 14, 8590 (2018).
- Zantop and Stark [2022] A. W. Zantop and H. Stark, Emergent collective dynamics of pusher and puller squirmer rods: swarming, clustering, and turbulence, Soft Matter 18, 6179 (2022).
- Worlitzer et al. [2021a] V. M. Worlitzer, G. Ariel, A. Be’er, H. Stark, M. Bär, and S. Heidenreich, Motility-induced clustering and meso-scale turbulence in active polar fluids, New J. Phys. 23, 033012 (2021a).
- Worlitzer et al. [2021b] V. M. Worlitzer, G. Ariel, A. Be’er, H. Stark, M. Bär, and S. Heidenreich, Turbulence-induced clustering in compressible active fluids, Soft Matter 17, 10447 (2021b).
- Huang et al. [2010] C.-C. Huang, A. Chatterji, G. Sutmann, G. Gompper, and R. G. Winkler, Cell-level canonical sampling by velocity scaling for multiparticle collision dynamics simulations, J. Comput. Phys. 229, 168 (2010).
- Omelyan [1998] I. P. Omelyan, Algorithm for numerical integration of the rigid-body equations of motion, Phys. Rev. E 58, 1169 (1998).
- Lamura et al. [2001] A. Lamura, G. Gompper, T. Ihle, and D. M. Kroll, Multiparticle collision dynamics: Flow around a circular and a square cylinder, Europhys. Lett. 56, 319 (2001).
- Padding et al. [2005] J. T. Padding, A. Wysocki, H. Löwen, and A. A. Louis, Stick boundary conditions and rotational velocity auto-correlation functions for colloidal particles in a coarse-grained representation of the solvent, J. Phys.: Condens. Matter 17, S3393 (2005).
- Theers and Winkler [2014] M. Theers and R. G. Winkler, Effects of thermal fluctuations and fluid compressibility on hydrodynamic synchronization of microrotors at finite oscillatory Reynolds number: A multiparticle collision dynamics simulation study, Soft Matter 10, 5894 (2014).
- Theers and Winkler [2015] M. Theers and R. G. Winkler, Bulk viscosity of multiparticle collision dynamics fluids, Phys. Rev. E 91, 033309 (2015).
- Shaebani et al. [2020] M. R. Shaebani, A. Wysocki, R. G. Winkler, G. Gompper, and H. Rieger, Computational models for active matter, Nat. Rev. Phys. 2, 181 (2020).
- Theers et al. [2016b] M. Theers, E. Westphal, G. Gompper, and R. G. Winkler, From local to hydrodynamic friction in Brownian motion: A multiparticle collision dynamics simulation study, Phys. Rev. E 93, 032604 (2016b).
- Batchelor [1959] G. K. Batchelor, The theory of homogeneous turbulence (University Press, Cambridge, 1959).
- Liu et al. [2021b] Z. Liu, W. Zeng, X. Ma, and X. Cheng, Density fluctuations and energy spectra of 3D bacterial suspensions, Soft Matter 17, 10806 (2021b).
Acknowledgements
Funding: The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. (www.gauss-centre.eu) for funding this project by providing computing time through the John von Neumann Institute for Computing (NIC) on the GCS Supercomputer JUWELS at Jülich Supercomputing Centre (JSC). Author contributions: G.G. and R.G.W. designed the research. E.W. wrote the simulation code. S.G. performed the simulations and analyzed the data. S.G., R.G.W., and G.G. discussed the results and wrote the manuscript together. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are available from the corresponding author upon reasonable request.