A lattice Boltzmann model for squirmers

Emergent collective dynamics of pusher and puller squirmer rods: swarming, clustering, and turbulence

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.

Alignment and propulsion of squirmer pusher-puller dumbbells

The properties of microswimmer dumbbells composed of pusher-puller pairs are investigated by mesoscale hydrodynamic simulations employing the multiparticle collision dynamics approach for the fluid. An individual microswimmer is represented by a squirmer, and various active-stress combinations in a dumbbell are considered. The squirmers are connected by a bond, which does not impose any geometrical restriction on the individual rotational motion. Our simulations reveal a strong influence of the squirmers' flow fields on the orientation of their propulsion directions, their fluctuations, and the swimming behavior of a dumbbell. The properties of pusher-puller pairs with equal magnitude of the active stresses dependent only weakly on the stress magnitude. This is similar to dumbbells of microswimmers without hydrodynamic interactions. However, for non-equal stress magnitudes, the active stress implies strong orientational correlations of the swimmers' propulsion directions with respect to each other as well as the bond vector. The orientational coupling is most pronounced for pairs with large differences of the active stress magnitude. The alignment of the squirmer propulsion directions with respect to each other is preferentially orthogonal in dumbbells with a strong pusher and weak puller, and antiparallel in the opposite case when the puller dominates. These strong correlations affect the active motion of dumbbells which is faster for strong pushers and slower for strong pullers.

Gyrotactic cluster formation of bottom-heavy squirmers

Squirmers that are bottom-heavy experience a torque that aligns them along the vertical so that they swim upwards. In a suspension of many squirmers, they also interact hydrodynamically via flow fields that are initiated by their swimming motion and by gravity. Swimming under the combined action of flow field vorticity and gravitational torque is called gyrotaxis. Using the method of multi-particle collision dynamics, we perform hydrodynamic simulations of a many-squirmer system floating above the bottom surface. Due to gyrotaxis they exhibit pronounced cluster formation with increasing gravitational torque. The clusters are more volatile at low values but compactify to smaller clusters at larger torques. The mean distance between clusters is mainly controlled by the gravitational torque and not the global density. Furthermore, we observe that neutral squirmers form clusters more easily, whereas pullers require larger gravitational torques due to their additional force-dipole flow fields. We do not observe clustering for pusher squirmers. Adding a rotlet dipole to the squirmer flow field induces swirling clusters. At high gravitational strengths, the hydrodynamic interactions with the no-slip boundary create an additional vertical alignment for neutral squirmers, which also supports cluster formation.

Interactions in active colloids

The past two decades have seen a remarkable progress in the development of synthetic colloidal agents which are capable of creating directed motion in an unbiased environment at the microscale. These self-propelling particles are often praised for their enormous potential to self-organize into dynamic nonequilibrium structures such as living clusters, synchronized super-rotor structures or self-propelling molecules featuring a complexity which is rarely found outside of the living world. However, the precise mechanisms underlying the formation and dynamics of many of these structures are still barely understood, which is likely to hinge on the gaps in our understanding of how active colloids interact. In particular, besides showing comparatively short-ranged interactions which are well known from passive colloids (Van der Waals, electrostatic etc), active colloids show novel hydrodynamic interactions as well as phoretic and substrate-mediated 'osmotic' cross-interactions which hinge on the action of the phoretic field gradients which are induced by the colloids on other colloids in the system. The present article discusses the complexity and the intriguing properties of these interactions which in general are long-ranged, non-instantaneous, non-pairwise and non-reciprocal and which may serve as key ingredients for the design of future nonequilibrium colloidal materials. Besides providing a brief overview on the state of the art of our understanding of these interactions a key aim of this review is to emphasize open key questions and corresponding open challenges.

Multi-particle collision dynamics with a non-ideal equation of state. II. Collective dynamics of elongated squirmer rods

Simulations of flow fields around microscopic objects typically require methods that both solve the Navier-Stokes equations and also include thermal fluctuations. One such method popular in the field of soft-matter physics is the particle-based simulation method of multi-particle collision dynamics (MPCD). However, in contrast to the typically incompressible real fluid, the fluid of the traditional MPCD methods obeys the ideal-gas equation of state. This can be problematic because most fluid properties strongly depend on the fluid density. In a recent article, we proposed an extended MPCD algorithm and derived its non-ideal equation of state and an expression for the viscosity. In the present work, we demonstrate its accuracy and efficiency for the simulations of the flow fields of single squirmers and of the collective dynamics of squirmer rods. We use two exemplary squirmer-rod systems for which we compare the outcome of the extended MPCD method to the well-established MPCD version with an Andersen thermostat. First, we explicitly demonstrate the reduced compressibility of the MPCD fluid in a cluster of squirmer rods. Second, for shorter rods, we show the interesting result that in simulations with the extended MPCD method, dynamic swarms are more pronounced and have a higher polar order. Finally, we present a thorough study of the state diagram of squirmer rods moving in the center plane of a Hele-Shaw geometry. From a small to large aspect ratio and density, we observe a disordered state, dynamic swarms, a single swarm, and a jammed cluster, which we characterize accordingly.

An effective and efficient model of the near-field hydrodynamic interactions for active suspensions of bacteria

Active suspensions of microswimmers demonstrate novel emergent behaviors (self-organizations, active turbulence, etc.) on macroscopic length scales. For such systems with, minimally, thousands of microswimmers, direct numerical simulations of the hydrodynamic interactions are computationally infeasible, and reduced models are needed. We demonstrated that existing models are not satisfactory in describing the hydrodynamic interactions for microswimmers in close proximity with even qualitatively erroneous predictions, indicating a pressing need for an adequate model. We propose a model that is both physically effective and computationally efficient in describing such hydrodynamics. The main novelty of our model is the description of hydrodynamic interactions through a resistance tensor, as opposed to an effective steric interaction in existing models.

Near-field hydrodynamic interactions in active fluids are essential to determine many important emergent behaviors observed, but have not been successfully modeled so far. In this work, we propose an effective model capturing the essence of the near-field hydrodynamic interactions through a tensorial coefficient of resistance, validated numerically by a pedagogic model system consisting of an Escherichia coli bacterium and a passive sphere. In a critical test case that studies the scattering angle of the bacterium–sphere pair dynamics, we prove that the near-field hydrodynamics can make a qualitative difference even for this simple two-body system: Calculations based on the proposed model reveal a region in parameter space where the bacterium is trapped by the passive sphere, a phenomenon that is regularly observed in experiments but cannot be explained by any existing model. In the end, we demonstrate that our model also leads to efficient simulation of active fluids with tens of thousands of bacteria, sufficiently large for investigations of many emergent behaviors.

High-order simulation scheme for active particles driven by stress boundary conditions

We study the dynamics and interactions of elliptic active particles in a two dimensional solvent. The particles are self-propelled through prescribing a fluid stress at one half of the fluid-particle boundary. The fluid is treated explicitly solving the Stokes equation through a discontinuous Galerkin scheme, which allows to simulate strictly incompressible fluids. We present numerical results for a single particle and give an outlook on how to treat suspensions of interacting active particles.

Hydrodynamic interactions in squirmer dumbbells: active stress-induced alignment and locomotion

Hydrodynamic interactions are fundamental for the dynamics of swimming self-propelled particles. Specifically, bonds between microswimmers enforce permanent spatial proximity and, thus, enhance emergent correlations by microswimmer-specific flow fields. We employ the squirmer model to study the swimming behavior of microswimmer dumbbells by mesoscale hydrodynamic simulations, where the squirmers' rotational motion is geometrically unrestricted. An important aspect of the applied particle-based simulation approach-the multiparticle collision dynamics method-is the intrinsic account for thermal fluctuations. We find a strong effect of active stress on the motility of dumbbells. In particular, pairs of strong pullers exhibit orders of magnitude smaller swimming efficiency than pairs of pushers. This is a consequence of the inherent thermal fluctuations in combination with the strong coupling of the squirmers' rotational motion, which implies non-exponentially decaying auto- and cross-correlation functions of the propulsion directions, and active stress-dependent characteristic decay times. As a consequence, specific stationary-state relative alignments of the squirmer propulsion directions emerge, where pullers are preferentially aligned in an antiparallel manner along the bond vector, whereas pushers are preferentially aligned normal to the bond vector with a relative angle of approximately 60° at weak active stress, and one of the propulsion directions is aligned with the bond at strong active stress. The distinct differences between dumbbells comprised of pusher or pullers suggest means to control microswimmer assemblies for future microbot applications.

Towards an analytical description of active microswimmers in clean and in surfactant-covered drops

Geometric confinements are frequently encountered in the biological world and strongly affect the stability, topology, and transport properties of active suspensions in viscous flow. Based on a far-field analytical model, the low-Reynolds-number locomotion of a self-propelled microswimmer moving inside a clean viscous drop or a drop covered with a homogeneously distributed surfactant, is theoretically examined. The interfacial viscous stresses induced by the surfactant are described by the well-established Boussinesq-Scriven constitutive rheological model. Moreover, the active agent is represented by a force dipole and the resulting fluid-mediated hydrodynamic couplings between the swimmer and the confining drop are investigated. We find that the presence of the surfactant significantly alters the dynamics of the encapsulated swimmer by enhancing its reorientation. Exact solutions for the velocity images for the Stokeslet and dipolar flow singularities inside the drop are introduced and expressed in terms of infinite series of harmonic components. Our results offer useful insights into guiding principles for the control of confined active matter systems and support the objective of utilizing synthetic microswimmers to drive drops for targeted drug delivery applications.

Squirmer rods as elongated microswimmers: flow fields and confinement

Microswimmers or active elements, such as bacteria and active filaments, have an elongated shape, which determines their individual and collective dynamics. There is still a need to identify what role long-range hydrodynamic interactions play in their fascinating dynamic structure formation. We construct rods of different aspect ratios using several spherical squirmer model swimmers. With the help of the mesoscale simulation method of multi-particle collision dynamics we analyze the flow fields of these squirmer rods both in a bulk fluid and in Hele-Shaw geometries of different slab widths. Based on the hydrodynamic multipole expansion either for bulk or confinement between two parallel plates, we categorize the different multipole contributions of neutral as well as pusher-type squirmer rods. We demonstrate how confinement alters the radial decay of the flow fields for a given force or source multipole moment compared to the bulk fluid.

Frequency-dependent higher-order Stokes singularities near a planar elastic boundary: Implications for the hydrodynamics of an active microswimmer near an elastic interface

The emerging field of self-driven active particles in fluid environments has recently created significant interest in the biophysics and bioengineering communities owing to their promising future for biomedical and technological applications. These microswimmers move autonomously through aqueous media, where under realistic situations they encounter a plethora of external stimuli and confining surfaces with peculiar elastic properties. Based on a far-field hydrodynamic model, we present an analytical theory to describe the physical interaction and hydrodynamic couplings between a self-propelled active microswimmer and an elastic interface that features resistance toward shear and bending. We model the active agent as a superposition of higher-order Stokes singularities and elucidate the associated translational and rotational velocities induced by the nearby elastic boundary. Our results show that the velocities can be decomposed in shear and bending related contributions which approach the velocities of active agents close to a no-slip rigid wall in the steady limit. The transient dynamics predict that contributions to the velocities of the microswimmer due to bending resistance are generally more pronounced than those due to shear resistance. Bending can enhance (suppress) the velocities resulting from higher-order singularities whereas the shear related contribution decreases (increases) the velocities. Most prominently, we find that near an elastic interface of only energetic resistance toward shear deformation, such as that of an elastic capsule designed for drug delivery, a swimming bacterium undergoes rotation of the same sense as observed near a no-slip wall. In contrast to that, near an interface of only energetic resistance toward bending, such as that of a fluid vesicle or liposome, we find a reversed sense of rotation. Our results provide insight into the control and guidance of artificial and synthetic self-propelling active microswimmers near elastic confinements.

Hydrodynamic mobility reversal of squirmers near flat and curved surfaces

Self-propelled particles have been experimentally shown to orbit spherical obstacles and move along surfaces. Here, we theoretically and numerically investigate this behavior for a hydrodynamic squirmer interacting with spherical objects and flat walls using three different methods of approximately solving the Stokes equations: The method of reflections, which is accurate in the far field; lubrication theory, which describes the close-to-contact behavior; and a lattice Boltzmann solver that accurately accounts for near-field flows. The method of reflections predicts three distinct behaviors: orbiting/sliding, scattering, and hovering, with orbiting being favored for lower curvature as in the literature. Surprisingly, it also shows backward orbiting/sliding for sufficiently strong pushers, caused by fluid recirculation in the gap between the squirmer and the obstacle leading to strong forces opposing forward motion. Lubrication theory instead suggests that only hovering is a stable point for the dynamics. We therefore employ lattice Boltzmann to resolve this discrepancy and we qualitatively reproduce the richer far-field predictions. Our results thus provide insight into a possible mechanism of mobility reversal mediated solely through hydrodynamic interactions with a surface.