Arne W.
Zantop
* and
Holger
Stark
*
Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstraße 36, 10623 Berlin, Germany. E-mail: a.zantop@tu-berlin.de; holger.stark@tu-berlin.de
First published on 5th July 2022
We study the interplay of steric and hydrodynamic interactions in suspensions of elongated microswimmers by simulating the full hydrodynamics of squirmer rods in the quasi two-dimensional geometry of a Hele-Shaw cell. To create pusher or puller-type squirmer rods, we concentrate the surface slip-velocity field more to the back or to the front of the rod and thereby are able to tune the rod's force-dipole strength. We study a wide range of aspect ratios and area fractions and provide corresponding state diagrams. The flow field of pusher-type squirmer rods destabilizes ordered structures and favors the disordered state at small area fractions and aspect ratios. Only when steric interactions become relevant, we observe a turbulent and dynamic cluster state, while for large aspect ratios a single swarm and jammed cluster occurs. The power spectrum of the turbulent state shows two distinct energy cascades at small and large wave numbers with power-law scaling and non-universal exponents. Pullers show a strong tendency to form swarms instead of the disordered state found for neutral and pusher rods. At large area fractions a dynamic cluster is observed and at larger aspect ratio a single swarm or jammed cluster occurs.
Active motion is always performed in non-equilibrium and, therefore, gives rise to new and interesting phenomena. For example, specially designed boundaries can rectify the random motion of active particles,19–23 or under gravity active particles develop polar order and even show inverted sedimentation profiles when they are bottom-heavy.24,25 In addition, microswimmer suspensions are subject to long-range hydrodynamic interactions with characteristic power-law decay.13,26,27 But also short-range steric interactions play an important role, in particular, for elongated particles, which align along each other.12 The combination of these interactions gives rise to numerous interesting dynamic patterns. Common examples are the formation of swarms or flocks,28–32 convection rolls and plumes,25 fluid pumps,33,34 vortices,35–38 active nematic patterns,39–41 and the emergence of the so-called active turbulence,42,43 termed in analogy to classical inertial turbulence.44,45 However, contrary to classical turbulence, where fluid flow is driven on the macroscopic scale, active turbulence is generated at the microscopic scale of the self-propelled particles and then energy is dissipated on larger scales. This mechanism causes a characteristic length scale for the formation of vortices and patterns, which is in contrast to the scale invariance of classical turbulence.46–48 In particular, the specific model parameters now determine the scaling of the velocity power spectrum, which is no longer universal.48–50 Active turbulence is found for microswimmers,51,52 active bio-filaments, which exist in the cytoplasm,39,53 and in growing tissue.54,55
On the theoretical side, the dynamics of active particles has been investigated at different levels of description. While the celebrated Vicsek model56–58 uses coarse-grained alignment rules, also models with explicit steric interactions exist.42,51,59–64 These models have also been generalized to active filaments31,65,66 and models, which implicitly include hydrodynamic interactions.67,68 Models, which directly simulate fluid flow and thereby explicitly include hydrodynamic interactions, use the method of multi-particle collision dynamics (MPCD)25,30,32,69–79 or the lattice-Boltzmann method.80–83 Continuum models combine elements of the Toner-Tu84 and Swift-Hohenberg85 equations to generate hydrodynamic equations for active suspensions.46,86,87 Microswimmers are also distinguished by their swimming mechanisms and the flow field, they generate in the surrounding fluid.11,13 While some microorganisms propel with cilia located all over their surface, common bacteria and algae propel with flagella that extent from the front or back of the cell body. In the first case, the flow field of a source dipole is realized, which decays as r−3, while bacteria and algae are termed pusher or puller-type swimmers that generate a long-range force-dipole flow field, which decays as r−2.13,27 The specific form of these hydrodynamic multipole flow fields determines the collective dynamics of the microswimmers.32,69,74,88–90
In this context, we proposed in ref. 76 the squirmer rod as a realistic microscopic model for elongated microswimmers. It consists of overlapping squirmers and thereby extends the well-known spherical squirmer model for ciliary propulsion26,91 and its implementation in MPCD.92 Restricting the surface slip-velocity to the front or the back of the squirmer rod, puller and pusher-type squirmer rods can be realized, respectively. The hydrodynamic flow field is simulated using an efficient implementation of MPCD with a reduced compressibility such that also collective dynamics in large and dense systems can be studied.93 Most recently, we presented the state diagram of neutral squirmer rods and identified with increasing area fraction and depending on the aspect ratio of the rods the disordered state, dynamic swarms, a single swarm, and ultimately a jammed cluster.94
In this article we generalize the squirmer-rod model to pushers and pullers with tunable force-dipole strength in the flow field. Compared to spheroidal squirmers it has the advantage that it better approximates real rodlike microswimmers such as E. coli. Using this model, we provide a comprehensive study of the state diagram for hydrodynamically interacting microswimmers over a wide range of aspect ratios, densities, and force-dipole strengths as the state diagrams in Fig. 3(a), 7(a) and 9(a) show. Hereby, we go well beyond previous works, which focused mainly on more dilute systems95,96 or a single aspect ratio.97 Thus, our work provides an overall view how hydrodynamics and shape determine the dynamic states of microswimmers. For example, our particle-based model nicely illustrates that the active turbulent state occurs as a compromise between the disordering hydrodynamic pusher–pusher interactions and aligning steric interactions. We also demonstrate that the two distinct energy cascades at low and large wave numbers in the power spectral density of velocity fluctuations exhibits non-universal exponents. Moreover, we find a dynamic cluster state at large densities. Besides this state all other states are also found in dry active rods.98 But, in our case, they have a specific contribution from hydrodynamics as, for example, the turbulent and swarming states show. The overall appearance of the states in our state diagrams can be summarized as follows. For pushers we observe that the swarming states of neutral squirmer rods are destabilized. Instead, for smaller aspect ratio between the disordered and dynamic cluster state, we observe the turbulent state as already mentioned. At high aspect ratios, where steric interactions become more relevant, we recover the single swarm and jammed cluster state of neutral squirmer rods. For pullers, hydrodynamic interactions stabilize the swarming state even for our smallest area fraction. Thus, compared to neutral squirmer rods, states are shifted towards lower densities. Variation of the force-dipole strength at constant aspect ratio supports all these findings.
The article is organized as follows. In Section 2 we give a brief overview of the methods used in this paper. Section 3 provides a detailed study of the different dynamic states of the squirmer rods as a function of their area fraction, aspect ratio, and force-dipole strength. We end with a summary and conclusions in Section 4.
![]() | ||
Fig. 1 (a) Sketch of the surface slip velocity of three squirmer rods with respective swimmer-type parameters χ = −1, 0, and 1 (from left to right). The arrows within the squirmer rods indicate the swimming direction ê. The parameter χ ∈ {−1, 1} can be used to smoothly vary the type of the squirmer rod from pusher to neutral to puller. (b) Force-dipole coefficient A2(χ) (blue) and active velocity v(χ)/v(0) (red) as a function of the swimmer parameter χ. (c) Schematic state diagram for neutral squirmer rods adapted from ref. 94. |
The squirmer rods propel through the axisymmetric and tangential slip velocity field at the surface of individual spherical squirmers,
vs = Bs1[(ê·![]() ![]() | (1) |
As described so far, the surface slip velocity of the squirmer rod resembles ciliated microorganisms such as Paramecium. In this realization, the profile of the slip velocity generates a flow field in the surrounding fluid, the far field of which can be described by a source dipole and an additional force quadrupole singularity, which both decay with |usd|, |ufq| ∼ r−3, as we show in ref. 76.
However, other prominent microswimmers such as E. coli bacteria or Chlamydomonas algae propel by rotating or beating flagella that extent from the back or the front of their bodies, respectively. These modes of propulsion create a pair of opposing forces that generate the more long-ranged force-dipole flow field |ufd| ∼ r−2. To generalise the squirmer rod model to these pusher and puller-type microswimmers, we concentrate the surface slip velocity either to the back or the front of the squirmer rod [see Fig. 1(a)]. This is done by multiplying the surface flow field with the envelope function
![]() | (2) |
In the present work, we consider the collective dynamics of squirmer rods confined between two parallel walls. In this so-called Hele-Shaw cell the radial decays of hydrodynamic multipoles are modified compared to the bulk fluid, such that the confined source dipole and force dipole decay as |ũsd| ∼ ρ−2 and |ũfd| ∼ ρ−3, respectively, where ρ is the polar distance76,99 (see appendix of ref. 99). As a consequence, the source dipole has the longest range in the flow field and ultimately dominates the far field. However, already in our previous work we realized that a distance of Δz = 6R between the walls, which we will use in our simulations, alters the relative strength of the source and force dipoles.76 As a consequence, the force dipole dominates the flow field at short and medium distance as we will demonstrate in Section 3.2.
The MPCD method uses a sequence of streaming and collision steps of the point-like fluid particles. In the streaming step the fluid particles move ballistically with their velocities during time step Δt. Then, the simulation box is divided into cubic cells with edge length a0 and the velocities of the fluid particles in one cell are modified randomly but keeping the mean velocity or linear momentum fixed. In this work, we use a collision rule optimized to achieve a low compressibility of the fluid.93 As in other MPCD methods, it includes angular momentum conservation and a thermostat. By choosing the MPCD fluid density n0 = 20/a03, i.e., on average 20 fluid particles per cell, and the time step , we obtain a fluid viscosity of
, which is comparable to previous work.74,76,96 Here, m0 is the mass of one fluid particle.
The immersed squirmer rods are modeled with the mass density of the fluid ρ0 = m0n0. The radius of the constituting squirmers is chosen as Rsq = 3a0 and the strength of the surface slip velocity as Bs1 = 0.1, which results in a Péclet number of Pe ≈ 350.94 Steric repulsion is realized with a Weeks–Chandler–Andersen potential103 that acts between two squirmers of different rods. To ensure that there is no significant overlap between two squirmer rods, we choose a strong force constant εWCA ≈ 104kBT.
The squirmer rods move because they acquire momentum from the surrounding fluid. In the MPCD streaming step, linear and angular momentum is transferred to the squirmer rod by collisions with the fluid particles. We achieve this by applying the so-called bounce back rule, which we also modify accordingly to implement the slip velocity in eqn (1) on the surface of the squirmer rods.100,101 Squirmer rods also contain so-called “ghost” particles, which improve the implementation of no-slip boundary conditions.101 During the MPCD collision step, they exchange momentum with the fluid, which is ascribed to the squirmer rod. Lastly, the equations of motion for the squirmer rods are integrated with a refined time step of δtMD = Δt/3 using a symplectic splitting algorithm for rigid body molecular dynamics.104
For our simulations we use two different geometries. To determine how the force dipole strength A2 depends on the swimmer-type parameter χ, we use a cubic box of linear size L = 100a0 with periodic boundary conditions along all spatial directions. In this case we use the time 104Δt to equilibrate the MPCD fluid flow fields and then average the fluctuating flow fields over additional 5 × 105Δt time steps during simulations.
To simulate the collective dynamics of the squirmer rods, we consider the quasi two-dimensional geometry of a Hele-Shaw cell of linear size L = 300a0 in the x and y direction. Along the z direction the system is confined by walls separated by a distance Δz = 6R to mimic the experimental setups using microfluidic chambers51,105 or liquid–oil interfaces.106,107 It has also been used in previous work.94,97 Additionally, this realization guarantees a strong contribution of the force-dipole interaction in the near field, as we will show in Section 3.2 and Fig. 2. In the x and y direction periodic boundary conditions are employed. In this geometry, the N squirmer rods are confined to only move in the midplane of the Hele-Shaw cell by a strong harmonic potential. The rods' initial positions are generated randomly for area fractions ϕ = NAsw/L2 < 0.6, where Asw is the two-dimensional cross section of one squirmer rod. For ϕ ≥ 0.6 squirmer rods are placed on a rectangular lattice, all with a randomly chosen orientation either parallel or anti-parallel to one of the major axis of the unit cell. We simulate for a time of 107Δt while saving snapshots every 2500Δt for further analysis. To ensure that the system is equilibrated, we omit the first 100 snapshots from the analysis.
To improve statistics in the study of the emergent turbulent patterns of pusher-type squirmer rods, we perform two additional simulation runs for all turbulent states and their neighboring points in the (α, ϕ) parameter space. Furthermore, for all these cases we also perform three simulation runs with an increased system size of L = 600a0 to investigate finite-size effects.
Indeed, we find a linear relation of the swimmer-type parameter χ and the force-dipole coefficient χ ∼ A2 [Fig. 1(b)]. This is expected since the terminal values of the envelope function f(x) in eqn (2), which determine the strength of the force dipole, are linear in χ. Furthermore, we find that due to the definition of the envelope function with 〈f(x)〉 = 1, the swimming velocity is nearly independent of χ [cf.Fig. 1(b), red curve].
Fig. 2(a) shows the flow field of a pusher-type squirmer rod (χ = −1.0, α = 3.25) swimming in Hele-Shaw geometry with wall distance Δz = 6R. The force-dipole field visibly dominates the flow field with its characteristic outwards directed streamlines along the rod and inwards directed streamlines at the side. From the simulated flow field, we extracted the radial parts of different hydrodynamic multipoles ũρ,n(ρ) following our previous work,76 and arrived at the curves shown in Fig. 2(b). The radial part of the force dipole ∼ρ−3 (red symbols) dominates the flow field up to a distance of circa 8lS, where it is exceeded by the more long-ranged field of the source dipole ∼ρ−2 (blue symbols). For smaller cell heights Δz, this crossover occurs at smaller distances. Since we are interested in exploring the effect of the hydrodynamic force dipole, we keep Δz = 6R for the rest of this work, such that the force dipole dominates the flow field close to the squirmer rod.
Most interestingly, compared to neutral squirmer rods the flow field of the pusher rod obviously suppresses the formation of dynamic swarms and impedes the single swarm state, which only occurs for large α ≥ 4.0. This is in agreement with findings of Saintillan and Shelley,88 who employ slender-body theory to show that polar and nematic order in systems of elongated pusher microswimmers is destroyed by their hydrodynamic flow fields. Likewise, explicit hydrodynamic simulations of collective dynamics of spherical squirmers show that pushers create disordered homogeneous systems,13,109 while pullers show swarming.109 An argument for this difference in the collective dynamic behavior is found in both implicit90 and explicit71 hydrodynamic simulations, which show that pushers are deflected during collision, while pullers align.
Instead of the suppressed swarming states for α ≤ 3.25, we find a very dynamic or turbulent state. A typical snapshot of a rod configuration is shown in Fig. 4(a), a video is provided in the Video 1 (ESI†). Here, the hydrodynamic dipole–dipole interactions compete with steric interactions that favor the formation of single swarms at large area fraction ϕ in the case of neutral rods.94 Towards lower ϕ the turbulent state transitions to the disordered state and towards higher aspect ratio α to cluster or swarm states.
Similar to other examples of active turbulence in theory and experiments,47,49–51,97 we find two cascades in the power spectral density of the squirmer velocities, as we will show in detail in Section 3.3.2. For the turbulent state at α = 3.25 and ϕ = 0.6 we construct a continuous velocity field v(x, y) to visualize the turbulent flow pattern and its vortices. A snapshot of the system and the resulting flow field are shown in Fig. 4(a) and (b), respectively. We also calculate the vorticity ω = (∇ × v)z, which is shown in Fig. 4(c). To easily obey the periodic boundary conditions, the vorticity was determined via a Fourier transformation. In all turbulent states, we observe that squirmer rods show a local alignment, which extends over short streaks, where squirmer rods swim side by side and head to tail. However, these streaks buckle and dissolve frequently, leaving a chaotic pattern of streaks and vortices.
For large aspect ratio α ≥ 4.0 we still find the single swarm state [Fig. 3(b) and Video 2, ESI†], where a stable cluster moves through the system, and also a jammed cluster state. Since they are also found for neutral squirmer rods,76 we deduce that they mainly arise from steric interactions. Indeed, they are also found for dry active rods.98 For aspect ratios below α < 4.75 and large area fractions clusters still form, but they deform dynamically while squirmer rods join or leave constantly [Fig. 3(c) and Video 3, ESI†]. However, the power spectral densities of this dynamic cluster state lacks the characteristic scaling behavior of the turbulent state, as we will show in Section 3.3.2. The absence of two power-law regimes (energy cascades) in the power spectra separated by a maximum, which defines a characteristic length scale, is our main criterion to distinguish the dynamic cluster state from the turbulent state. Hence, this state is intermediate between the jammed cluster and turbulent state. In the jammed cluster state observed for α = 4.75 and ϕ = 0.77 the average swimmer velocity is low, 〈vi·êi〉 ≈ 0.18v0, and rises to 〈vi·êi〉 ≈ 0.4v0 for the aspect ratio α ≤ 4 in the dynamic cluster state.
As already mentioned, the disordered state extends to higher densities compared to neutral squirmer rods. Here, velocity correlations between squirmer rods are short-ranged and the power spectral density decays without any algebraic behavior as we will show in Section 3.3.2. In other words, no patterns emerge because energy is not transported to larger scales.
All together, our results imply that the pusher-type flow fields inhibit or destabilize the steric alignment of swimmers that has been found for neutral squirmer rods. In intermediate regions of the state diagram, the competition of both effects lead to new dynamic states.
![]() | (3) |
![]() | (4) |
In Fig. 4 we compare the normalized velocity pair-correlation functions v(R) = Cv(R)/〈vi2〉 of pusher-type squirmer rods for aspect ratio α = 3.25 at different area fractions ϕ. A more complex behavior is observed for an area fraction in the range ϕ ∈ [0.4, 0.77], which allows for competing steric and hydrodynamic interactions and for which we show
v(R) in Fig. 5. The respective states are characterized as disordered (ϕ = 0.4/0.5), turbulent (ϕ = 0.6/0.7), and dynamic cluster ϕ = 0.77, as we detail also further below. Generally, for larger ϕ, the correlations extend to larger distances, which indicates that the size of ordered patterns increases with density. Distinctive for pushers we find that
v(R) is not purely positive but also shows distances where it becomes negative. In more detail, we observe the following characteristics. For distances, where steric interactions occur and dominate, R < lS,
v(R) is positive. With increasing R it decays to zero around R = lS to R = 3lS and becomes negative for area fractions ϕ up to 0.7. In the turbulent states for ϕ = 0.6 and 0.7,
v(R) returns to positive values at more than twice the first zero-crossing distance [inset of Fig. 5]. For ϕ = 0.6 we even observe a third zero crossing. In the disordered states for ϕ = 0.4 and 0.5
v(R) does not show a second zero crossing but approaches zero from the negative region. Lastly, at ϕ = 0.77 the system is in the dynamic cluster state with densely packed squirmer rods. This is in agreement with the pair-correlation function, which has the longest range and only exhibits anti-correlations for R > 6lS. However, we find that the orientational autocorrelation function 〈ê(t)·ê(t + δt)〉 (not shown) decays after the cluster moved a distance of 3–5lS, considering its mean velocity of 〈vi·êi〉 ≈ 0.4v0 for α = 4.0. This clearly indicates that the clusters are dynamic. For other aspect ratios α only small quantitative changes occur. Especially, the correlation functions for the turbulent and disordered states are representative for all α and ϕ.
The power spectral densities E(k) calculated from Cv(R) are in accordance with the findings so far. In Fig. 6(a) we compare the power spectral densities E(k) for three different aspect ratios α = 1.75, 2.5, and 3.25 in dense systems. The velocity autocorrelation functions v(R) show the same oscillating decay as already observed above for α = 3.25 and ϕ = 0.6. We called these states active turbulence since the power spectral density E(k) shows a broad peak, which separates two regions from each other with power-law decays towards low and high k. This is characteristic of active turbulence. Close to the characteristic wave number kc = 2π/lS, where energy is inserted from the active motion of the squirmer rods, E(k) decays to zero. The maxima of the power spectral density all roughly occur at 0.2kc, indicating that the pattern size scales with lS, which makes sense when the density is roughly the same.
As reported in other experimental51 and theoretical50 studies on polar active fluids, the scaling exponents for the power laws in k are not universal. Towards smaller length scales (larger k), the exponent shows only a weak dependence on the aspect ratio α with values of −1.1 and −1.3, while towards larger length scales we observe a stronger variation of the scaling exponent with aspect ratio. For short squirmer rods with α = 1.75, the power spectral density decays more rapidly with E(k) ∼ k1.9 and for a larger aspect ratio α = 3.25, the power spectral density scales with E(k) ∼ k1. This indicates that for larger aspect ratio patterns at large length scales (small k) are observed more frequently.
For disordered states the power spectra do not follow power laws [cf.Fig. 6(b)]. Instead, they quickly decay towards smaller k, which indicates again that energy is not transported to larger length scales. In contrast, for the dynamic cluster state energy is concentrated at the larger length scales (small k) as expected, while E(k) decays strongly towards smaller length scales [cf.Fig. 6(c)]. To investigate finite size effects, we compare the regular systems of size L = 300a0 to systems with the increased system size of L = 600a0 performed for α = 4 and ϕ = 0.7. In the bigger system the maximum of E(k) shifts to smaller values of k roughly proportional to the change in L which suggests that the maximum correlation length always corresponds to the system size, and no characteristic length scale emerges. At the same time, we observe a decrease of max[E(k)] indicating that the dynamic of the cluster slows down as L increases. This behavior clearly distinguishes the dynamic cluster state from active turbulence, which has an intrinsic characteristic length scale. However, we note that data for larger systems would be necessary to show that the maximum correlation length 1/argmax[E(k)] diverges for infinite system size. Though, this currently exceeds our computational capabilities. Note, the irregular behavior in the curve for L = 300a0 is due to the fact that with increasing α close to the jammed cluster state the dynamics of the clusters slows down and longer simulation times are needed to smooth the curve.
The most prominent feature of puller-type squirmer rods is that their flow fields strongly promote the formation of small swarms. Even for dilute systems with area fractions down to 0.2, the hydrodynamic interactions between the squirmer rods favor small swarms against the disordered state. Different from the purely steric alignment responsible for swarming neutral squirmer rods,94 pairs of puller-type squirmer rods form T-shaped configurations, although deformed, when they collide. Examples are indicated in the snapshot of Fig. 7(b). This leads to small groups, where several squirmer rods follow a central leader of the group, arranged in a “comet tail” and pointing inwards [Fig. 7(b)]. A similar behavior is also found in ref. 95 and for dry active rods.98 However, we observe that for the hydrodynamically interacting squirmer rods the swarming behavior extends to smaller densities and aspect ratios compared to dry self-propelled rods.61,64 Often the swarms have an asymmetric shape, which then induces curved trajectories (Video 4, ESI†). With increasing area fraction, ϕ ≥ 0.4, swarms also increase in size. However, the overall alignment in the swarms rather decreases. This is due to the low polarity of the small swarms, which then merge in an unordered fashion with increasing ϕ.
For aspect ratios α ≥ 3 we also observe the formation of a single swarm, which collects all squirmer rods in the system [Fig. 7(d) and Video 5, ESI†]. For area fractions ϕ ≥ 0.6 and sufficiently large aspect ratios the single swarm becomes jammed and forms a static cluster state [Fig. 7(e) and Video 6, ESI†]. Reducing the aspect ratio at ϕ = 0.77 to α ≤ 2.5, the jammed cluster becomes dynamic similar to the state observed in pusher-type squirmer rods [Fig. 7(c) and Video 7, ESI†].
We have already seen, that even in dilute systems with ϕ = 0.2, the puller-type flow fields locally align the squirmer rods and thereby cause swarming. This is also visible in the velocity correlation function v(R), which is more long-ranged compared to pusher-type rods [Fig. 5(a)]. Thus, squirmer rods are significantly correlated for distances up to R = 2lS. As described in Section 3.4.1, the size of swarms increases with density. For the subsequent area fraction of ϕ = 0.4 the system is already in the single swarm state at the chosen aspect ratio of 4.0. Hence, the velocity correlation function
v(R) becomes more long-ranged. At ϕ = 0.6 and ϕ = 0.77, when the system is noticeably in the jammed cluster state, we observe two different features. For ϕ = 0.6,
v(R) reaches a plateau, which indicates that the jammed cluster drifts in one direction since the squirmer rods are aligned, on average. At the highest density ϕ = 0.77 the cluster formation takes place more quickly resulting in a disordered structure without drift. Here, the velocity correlation function
v(R) decreases and does not exhibit a plateau until R/lS = 5.
In dilute systems with ϕ = 0.2, Ce(δt) decays to zero at around δt = 15lS/v0, meaning that swarms swim a distance of around 7lS until they either rotate, merge with other clusters or dissolve (Video 8, ESI†). For the single swarm state at ϕ = 0.4, we observe a more long-ranged autocorrelation. After a swarm has formed it remains stable but also goes through configurational changes so that partially straight motion is observed interrupted by rotations or turns of the swarm (Video 5, ESI†). As a result, Ce(δt) exhibits a damped oscillation, where the first minimum at δt ≈ 37lS/v0 means that the swarm has rotated by 180°. In the jammed cluster state, we again observe two different features. For ϕ = 0.6, the cluster drifts due to its overall polar order. However, also configurational changes occur, which results in a decay of Ce(δt) towards a non-zero value meaning that the cluster has reached some stable configuration within the observation window. At the higher density ϕ = 0.77, the cluster now has a more random but static arrangement of the rods. Therefore, it hardly drifts and it only rotates very slowly, which causes a very slow decrease of Ce(δt).
Comparing dense systems (ϕ = 0.77) of different aspect ratios to each other, we observe the following characteristics from high to low α [Fig. 8(c)]. For the largest aspect ratio α = 4.75 steric interactions between the rods are most dominant and we find an entirely jammed configuration with constant Ce(δt) = 1. Lowering the aspect ratio, the jammed cluster becomes more and more interrupted by occasional configurational changes, which results in the noticeable but slow decay of Ce(δt) for α = 3.25. In contrast, the dynamic cluster states at α = 1.75 and α = 2.5 show a fast exponential decay of Ce(δt) due to the dynamic rearrangements within the cluster. The decay is faster for the lower aspect ratio, where the rods reorient more easily.
For puller-type squirmer rods with χ > 0, we observe that already weak pullers with χ = 0.25 show swarming at the small area fraction of ϕ = 0.2, i. e., the same state we observed in Section 3.4.1 for the highest strength χ = 1.0. Obviously, the force-dipole field, which attracts nearby rods along the rod axis, promotes swarming. However, with increasing ϕ, swarming also extends to weak pusher rods. At higher area fractions ϕ ≥ 0.6, we find single swarm and jammed cluster states due to steric interactions between the rods, where neutral (χ = 0) and puller-type squirmer rods show qualitatively similar behavior.
For pusher-type squirmer rods, χ < 0, we observe a more diverse behavior, which we already noted for χ = 1 in Section 3.3.1. While at ϕ = 0.2 only the disordered state occurs, at ϕ = 0.4 we find a transition from swarming to turbulent to disordered state with increasing pusher strength. The analysis of the power spectral density E(k) reveals that the system with χ = −0.5 is indeed in the turbulent state. Thus, the pusher flow field destabilizes swarming clusters. As a compromise between swarming and disordered states, active turbulence occurs in conjunction with steric repulsion. Increasing density further to ϕ = 0.6, steric repulsion stabilizes swarming clusters and active turbulence is shifted to the largest pusher strength χ = −1. Finally, at ϕ = 0.77 the pusher flow field destabilizes the jammed cluster of neutral and puller-type squirmer rods, which then becomes dynamic.
In the dilute systems, all the velocity pair-correlation functions in the swarming state of puller-type swimmers (χ > 0) are nearly identical and show correlations up to circa 2lS [Fig. 9(b)]. This indicates that the mechanism behind swarm formation is fully established already at small dipole strength. Pusher-type squirmer rods (χ <0) only exhibit the disordered state and no swarms, so v(R) is more short-ranged compared to pullers. Pushers prefer to order side by side and therefore pronounced peaks at closest distance R = lS/α are observed compared to puller rods. Furthermore, we observe anti-correlations with
v(R) < 0 starting at distances R between lS and 2lS. They belong to pusher rods approaching each other head on.
At the higher area fraction ϕ = 0.6, we observe the single swarm, swarming, and turbulent states [Fig. 9(c)]. In the single swarm state observed for weak pullers (χ = 0.25), v(R) is long-ranged and shows significant correlations for all recorded distances up to R = 5lS. In the swarming state that emerges for pushers at χ = −0.25 and −0.5 the correlation length of
v(R) decreases with increasing force-dipole strength. At χ = −0.5 it is around R = 2lS. In the turbulent state at χ = −1,
v(R) becomes even more short-ranged and negative at R = lS due to the occurrence of vortices.
In Fig. 9(d) we compare the power spectral densities E(k) for three systems in the turbulent state at α = 3.25. All three spectra show a regular cascade towards large k (small scales) and a second cascade towards small k (large scales). At large k all three power spectral densities show roughly the same scaling behavior E(k) ∼ k−1.3. Furthermore, the curves for (χ, ϕ) = (−0.5, 0.4) and (χ, ϕ) = (−1.0, 0.6) coincide in the range from k = 0.2kc to kc, which includes the maximum of both spectra. For the second cascade, different scaling exponents are observed. At χ = −1 the power spectral density decays more rapidly in the system with ϕ = 0.6 following a scaling E(k) ∼ k2, while in the denser system (ϕ = 0.7) the scaling is E(k) ∼ k1. Furthermore, the maximum of E(k) shifts to a smaller wave number k, which implies a larger intrinsic length scale of the turbulent pattern. For the system with the smaller pusher strength, χ = −0.5, the scaling follows E(k) ∼ k1.4, thus, it is situated between the two cases just discussed.
For pushers with largest force-dipole strength χ = −1, we observe a new turbulent state along with other modifications in the state diagram compared to neutral squirmer rods. In dilute systems, where steric interactions are less important, the pusher flow field suppresses the formation of the swarming states so that the disordered state extends to larger area fractions and aspect ratios. When steric interactions and thereby steric alignment of the rods become more important with increasing area fraction, the turbulent state emerges at small aspect ratios. To observe it for weaker pusher strengths at constant aspect ratio, one has to decrease the area fraction, which confirms the importance of steric interactions for the turbulent state. The velocity pair-correlation function decays and exhibits negative regions after one or two swimmer lengths, which indicates a characteristic length scale. At the same time, the power spectral density of the velocity fluctuations shows two energy cascades at small and large wave numbers with power-law scaling and non-universal exponents as reported in other works on active turbulence in polar fluids.49,50,97 Increasing the area fraction further, a transition to the dynamic cluster state occurs at medium aspect ratios, which is reminiscent of the turbulent state but without the characteristic energy cascades. This is in contrast to soft active rods, where jammed states occur,42,63 which are, however, destroyed for strong enough self-propulsion.112 At larger aspect ratios, steric interactions become more important. Here, pushers resemble the dynamics of neutral squirmer rods94 and dry active rods,61,62 such that, instead of the turbulent and dynamic swarm state, a single swarm and jammed cluster are observed.
For pullers at all studied force-dipole strengths χ and in the most dilute system at ϕ = 0.2, hydrodynamic interactions already promote the swarming state instead of the disordered state observed for neutral rods. At smaller aspect ratios and the largest area fraction, the dynamic cluster state occurs, while for larger aspect ratios swarming and dynamic clustering are replaced by a single swarm and jammed cluster state.
Thus, our study clearly indicates the importance of the type of the swimmer flow field for the occurring states. Although comparable states are often found for dry self-propelled rods, for squirmer rods their occurrence in the state diagram crucially depends on the specific hydrodynamic flow field. Specifically, active turbulence is only found for pushers and puller type flow fields greatly enhance the formations of swarms. Additionally, the different scaling exponents in the turbulent state show that steric interactions, tuned by rod density, play an important role for this state. All in all, the turbulent state occurs as a compromise between disordering hydrodynamic pusher flow fields and aligning steric interactions. Therefore, with our work we contribute to the insights how various biological propulsion strategies determine the collective motion of microswimmers.
In future work, it will be interesting to investigate collective dynamics of pusher and puller-type squirmer rods also in the bulk fluid, where the hydrodynamic force-dipole has the most long-range or dominant flow field. Likewise, implementing a liquid–air interface using a slip-boundary condition at one wall of the Hele-Shaw geometry, would provide a realistic modeling of recent experiments, where bacteria move close to a fluid–air interface.113
Footnote |
† Electronic supplementary information (ESI) available: We provide eight videos of different collective dynamic states in the top view. Videos 1–3 show pusher rods in the turbulent (1), single swarm (2) and dynamic cluster state (3). Videos 4–8 show puller rods in the swarming (4 and 8), single swarm (5), jammed cluster (6), and dynamic cluster state (7). See DOI: https://doi.org/10.1039/d2sm00449f |
This journal is © The Royal Society of Chemistry 2022 |