Hydrodynamics determines collective motion and phase behavior of active colloids in quasi-two-dimensional confinement

Tunable colloidal spinners: Active chirality and hydrodynamic interactions governed by rotating external electric fields

The rotational dynamics of microparticles in liquids have a wide range of applications, including chemical microreactors, biotechnologies, microfluidic devices, tunable heat and mass transfer, and fundamental understanding of chiral active soft matter which refers to systems composed of particles that exhibit a handedness in their rotation, breaking mirror symmetry at the microscopic level. Here, we report on the study of two effects in colloids in rotating electric fields: (i) the rotation of individual colloidal particles in rotating electric field and related to that (ii) precession of pairs of particles. We show that the mechanism responsible for the rotation of individual particles is related to the time lag between the external field applied to the particle and the particle polarization. Using numerical simulations and experiments with silica particles in a water-based solvent, we prove that the observed rotation of particle pairs and triplets is governed by the tunable rotation of individual particles and can be explained and described by the action of hydrodynamic forces. Our findings demonstrate that colloidal suspensions in rotating electric fields, under some conditions, represent a novel class of chiral soft active matter-tunable colloidal spinners. The experiments and the corresponding theoretical framework we developed open novel prospects for future studies of these systems and for their potential applications.

Collective behavior of squirmers in thin films

Bacteria in biofilms form complex structures and can collectively migrate within mobile aggregates, which is referred to as swarming. This behavior is influenced by a combination of various factors, including morphological characteristics and propulsive forces of swimmers, their volume fraction within a confined environment, and hydrodynamic and steric interactions between them. In our study, we employ the squirmer model for microswimmers and the dissipative particle dynamics method for fluid modeling to investigate the collective motion of swimmers in thin films. The film thickness permits a free orientation of non-spherical squirmers, but constraints them to form a two-layered structure at maximum. Structural and dynamic properties of squirmer suspensions confined within the slit are analyzed for different volume fractions of swimmers, motility types (e.g., pusher, neutral squirmer, puller), and the presence of a rotlet dipolar flow field, which mimics the counter-rotating flow generated by flagellated bacteria. Different states are characterized, including a gas-like phase, swarming, and motility-induced phase separation, as a function of increasing volume fraction. Our study highlights the importance of an anisotropic swimmer shape, hydrodynamic interactions between squirmers, and their interaction with the walls for the emergence of different collective behaviors. Interestingly, the formation of collective structures may not be symmetric with respect to the two walls. Furthermore, the presence of a rotlet dipole significantly mitigates differences in the collective behavior between various swimmer types. These results contribute to a better understanding of the formation of bacterial biofilms and the emergence of collective states in confined active matter.

Light-Induced Material Motion Fingerprint - A Tool Toward Selective Interfacial Sensitive Fractioning of Microparticles via Microfluidic Methods

In this article, a novel strategy is presented to selectively separate a mixture of equally sized microparticles but differences in material composition and surface properties. The principle relies on a photosensitive surfactant, which makes particles under light illumination phoretically active. The latter hovers microparticles from a planar interface and together with a superimposed fluid flow, particles experience a drift motion characteristic to its interfacial properties. The drift motion is investigated as a function of applied wavelength, demonstrating that particles composed of different material show a unique spectrally resolved light-induced motion profile. Differences in those motion profile allow a selective fractioning of a desired particle from a complex particle mixture made out of more than two equally sized different particle types. Besides that, the influence of applied wavelength is systematically studied, and discussed the origin of the spectrally resolved chemical activity of microparticles from measured photo-isomerization rates.

Collective dynamics of active dumbbells near a circular obstacle

In this article, we present the collective dynamics of active dumbbells in the presence of a static circular obstacle using Brownian dynamics simulation. The active dumbbells aggregate on the surface of a circular obstacle beyond a critical radius. The aggregation is non-uniform along the circumference, and the aggregate size increases with the activity (Pe) and the curvature radius (Ro). The dense aggregate of active dumbbells displays persistent rotational motion with a certain angular speed, which linearly increases with activity. Furthermore, we show a strong polar ordering of the active dumbbells within the aggregate. The polar ordering exhibits long-range correlation, with the correlation length corresponding to the aggregate size. Additionally, we show that the residence time of an active dumbbell on the obstacle surface increases rapidly with area fraction due to many-body interactions that lead to a slowdown of the rotational diffusion. This article further considers the dynamical behavior of a tracer particle in the solution of active dumbbells. Interestingly, the speed of the passive tracer particle displays a crossover from monotonically decreasing to increasing with the size of the tracer particle upon increasing the dumbbells' speed. Furthermore, the effective diffusion of the tracer particle displays non-monotonic behavior with the area fraction; the initial increase in diffusivity is followed by a decrease for a larger area fraction.

An Account on BiVO4 as Photocatalytic Active Matter

Photocatalytic materials are gaining popularity and research investment for developing light-driven micromotors. While most of the early work used highly stable TiO2 as a material to construct micromotors, mostly in combination with noble metals, other semiconductors offer a wider range of properties, including independence from high-energy UV light. This review focuses on our work with BiVO4 which has shown promise due to its small band gap and resulting ability to absorb blue light. Additionally, this salt's well-defined crystal structures lead to exploitable charge separation on different crystal facets, providing sufficient asymmetry to cause active propulsion. These properties have given rise to fascinating physical and chemical behaviors that show how rich and variable active matter can become. Here, we present the synthesis of different BiVO4 microparticles and their material properties that make them excellent candidates as active micromotors. A critical factor in understanding inherently asymmetric micromotors is knowledge of their flow fields. However, due to their small size and the need to use even smaller tracer particles to avoid perturbing the flow field, measuring flow fields at the microscale is a difficult task. We also present these first results, which allow us to demonstrate the correlation between chemical reactivity and the flow generated, leading to active motion. Due to the nontoxic nature of BiVO4, these visible-light-responsive microswimmers have been used to study the first steps toward applications, even in sensitive areas such as food technology. Although these initial tests are far from being realized, we have to face the fact that a single microswimmer will not be able to perform macroscale tasks. Therefore, we present the reader with the first simple studies of collective motion, hoping for many new contributions to the field. The one-step synthesis of BiVO4 clearly paves the way for studies requiring large numbers of particles. We predict that the combination of promising applications for a nontoxic material which is readily synthesized in large quantities will contribute pivotally to advance the field of active matter beyond the proof-of-concept stage.

Pair Interactions of Self-Propelled SiO2-Pt Janus Colloids: Chemically Mediated Encounters

Driven by the necessity to achieve a thorough comprehension of the bottom-up fabrication process of functional materials, this experimental study investigates the pairwise interactions or collisions between chemically active SiO2-Pt Janus colloids. These collisions are categorized based on the Janus colloids' orientations before and after they make physical contact. In addition to the hydrodynamic interactions, the Janus colloids are also known to affect each other's chemical field, resulting in chemophoretic interactions, which depend on the degree of surface anisotropy in reactivity of Janus colloid and the solute-surface interaction at play. Our study reveals that these interactions lead to a noticeable decrease in particle speed and changes in orientation that correlate with the contact duration and yield different collision types. Distinct configurations of contact during collisions were found, whose mechanisms and likelihood are found to be dependent primarily on the chemical interactions. Such estimates of collision and their characterization in dilute suspensions shall have a key impact in determining the arrangement and time scales of dynamical structures and assemblies of denser suspensions and potentially the functional materials of the future.

Persistent and responsive collective motion with adaptive time delay

It is beneficial for collective structures to simultaneously have high persistence to environmental noise and high responsivity to nontrivial external stimuli. However, without the ability to differentiate useful information from noise, there is always a tradeoff between persistence and responsivity within the collective structures. To address this, we propose adaptive time delay inspired by the adaptive behavior observed in the school of fish. This strategy is tested using particles powered by optothermal fields coupled with an optical feedback-control system. By applying the adaptive time delay with a proper threshold, we experimentally observe the responsivity of the collective structures enhanced by approximately 1.6 times without sacrificing persistence. Furthermore, we integrate adaptive time delay with long-distance transportation and obstacle-avoidance capabilities to prototype adaptive swarm microrobots. This research demonstrates the potential of adaptive time delay to address the persistence-responsivity tradeoff and lays the foundation for intelligent swarm micro/nanorobots operating in complex environments.

Controlling active turbulence by activity patterns

By patterning activity in space, one can control active turbulence. To show this, we use Doi's hydrodynamic equations of a semidilute solution of active rods. A linear stability analysis reveals the resting isotropic fluid to be unstable above an absolute pusher activity. The emergent activity-induced paranematic state displays active turbulence, which we characterize by different quantities including the energy spectrum, which shows the typical power-law decay with exponent -4. Then, we control the active turbulence by a square lattice of circular spots where activity is switched off. In the parameter space lattice constant versus surface-to-surface distance of the spots, we identify different flow states. Most interestingly, for lattice constants below the vorticity correlation length and for spot distances smaller than the nematic coherence length, we observe a multi-lane flow state, where flow lanes with alternating flow directions are separated by a street of vortices. The flow pattern displays pronounced multistability and also appears transiently at the transition to the isotropic active-turbulence state. At larger lattice constants a trapped vortex state is identified with a non-Gaussian vorticity distribution due to the low flow vorticity at the spots. It transitions to conventional active turbulence for increasing spot distance.

Isotropic active colloids: explicit vs. implicit descriptions of propulsion mechanisms

Modeling the couplings between active particles often neglects the possible many-body effects that control the propulsion mechanism. Accounting for such effects requires the explicit modeling of the molecular details at the origin of activity. Here, we take advantage of a recent two-dimensional model of isotropic active particles whose propulsion originates from the interactions between solute particles in the bath. The colloid catalyzes a chemical reaction in its vicinity, which results in a local phase separation of solute particles, and the density fluctuations of solute particles cause the enhanced diffusion of the colloid. In this paper, we investigate an assembly of such active particles, using (i) an explicit model, where the microscopic dynamics of the solute particles is accounted for; and (ii) an implicit model, whose parameters are inferred from the explicit model at infinite dilution. In the explicit solute model, the long-time diffusion coefficient of the active colloids strongly decreases with density, an effect which is not captured by the derived implicit model. This suggests that classical models, which usually decouple pair interactions from activity, fail to describe collective dynamics in active colloidal systems driven by solute-solute interactions.

Hydrodynamic Anisotropy of Depletion in Nonequilibrium

Active colloids in a bath of inert particles of smaller size cause anisotropic depletion. The active hydrodynamics of this nonequilibrium phenomenon, which is fundamentally different from its equilibrium counterpart and passive particles in an active bath, remains scarcely understood. Here we combine mesoscale hydrodynamic simulation as well as theoretical analysis to examine the physical origin for the active depletion around a self-propelled noninteractive colloid. Our results elucidate that the variable hydrodynamic effect critically governs the microstructure of the depletion zone. Three characteristic states of anisotropic depletion are identified, depending on the strength and stress of activity. This yields a state diagram of depletion in the two-parameter space, captured by developing a theoretical model with the continuum kinetic theory and leading to a mechanistic interpretation of the hydrodynamic anisotropy of depletion. Furthermore, we demonstrate that such depletion in nonequilibrium results in various clusters with ordered organization of squirmers, which follows a distinct principle contrary to that of the entropy scenario of depletion in equilibrium. The findings might be of immediate interest to tune the hydrodynamics-mediated anisotropic interactions and active nonequilibrium organizations in the self-propulsion systems.

Optimal navigation of a smart active particle: directional and distance sensing

We employ Q learning, a variant of reinforcement learning, so that an active particle learns by itself to navigate on the fastest path toward a target while experiencing external forces and flow fields. As state variables, we use the distance and direction toward the target, and as action variables the active particle can choose a new orientation along which it moves with constant velocity. We explicitly investigate optimal navigation in a potential barrier/well and a uniform/ Poiseuille/swirling flow field. We show that Q learning is able to identify the fastest path and discuss the results. We also demonstrate that Q learning and applying the learned policy works when the particle orientation experiences thermal noise. However, the successful outcome strongly depends on the specific problem and the strength of noise.

Hydrodynamic interactions between squirmers near walls: far-field dynamics and near-field cluster stability

Confinement increases contacts between microswimmers in dilute suspensions and affects their interactions. In particular, boundaries have been shown experimentally to lead to the formation of clusters that would not occur in bulk fluids. To what extent does hydrodynamics govern these boundary-driven encounters between microswimmers? We consider theoretically the symmetric boundary-mediated encounters of model microswimmers under gravity through far-field interaction of a pair of weak squirmers, as well as the lubrication interactions occurring after contact between two or more squirmers. In the far field, the orientation of microswimmers is controlled by the wall and the squirming parameter. The presence of a second swimmer influences the orientation of the original squirmer, but for weak squirmers, most of the interaction occurs after contact. We thus analyse next the near-field reorientation of circular groups of squirmers. We show that a large number of swimmers and the presence of gravity can stabilize clusters of pullers, while the opposite is true for pushers; to be stable, clusters of pushers thus need to be governed by other interactions (e.g. phoretic). This simplified approach to the phenomenon of active clustering enables us to highlight the hydrodynamic contribution, which can be hard to isolate in experimental realizations.

Particle-based mesoscopic model for phase separation in a binary fluid mixture

A mesoscopic simulation model to study the phase separation in a binary fluid mixture in three dimensions (3D) is presented here by augmenting the existing particle-based multiparticle collision dynamics (MPCD) algorithm. The approach describes the nonideal equation of the fluid state by incorporating the excluded-volume interaction between the two components within the framework of stochastic collision, which depends on the local fluid composition and velocity. Calculating the nonideal contribution to the pressure both from simulation and analytics shows the model to be thermodynamically consistent. A phase diagram to explore the range of parameters that give rise to phase separation in the model is investigated. The interfacial width and phase growth obtained from the model agree with the literature for a wide range of temperatures and parameters.

Field-driven cluster formation in two-dimensional colloidal binary mixtures

We study size- and charge-asymmetric oppositely charged colloids driven by an external electric field. The large particles are connected by harmonic springs, forming a hexagonal-lattice network, while the small particles are free of bonds and exhibit fluidlike motion. We show that this model exhibits a cluster formation pattern when the external driving force exceeds a critical value. The clustering is accompanied with stable wave packets in vibrational motions of the large particles.

The interplay between chemo-phoretic interactions and crowding in active colloids

Many motile microorganisms communicate with each other and their environments via chemical signaling which leads to long-range interactions mediated by self-generated chemical gradients. However, consequences of the interplay between crowding and chemotactic interactions on their collective behavior remain poorly understood. In this work, we use Brownian dynamics simulations to investigate the effect of packing fraction on the formation of non-equilibrium structures in a monolayer of diffusiophoretic self-propelled colloids as a model for chemically active particles. Focusing on the case when a chemical field induces attractive positional and repulsive orientational interactions, we explore dynamical steady-states of active colloids of varying packing fractions and degrees of motility. In addition to collapsed, active gas, and dynamical clustering steady-states reported earlier for low packing fractions, a new phase-separated state emerges. The phase separation results from a competition between long-range diffusiophoretic interactions and motility and is observed at moderate activities and a wide range of packing fractions. Our analysis suggests that the fraction of particles in the largest cluster is a suitable order parameter for capturing the transition from an active gas and dynamical clustering states to a phase-separated state.

Squirmer hydrodynamics near a periodic surface topography

The behaviour of microscopic swimmers has previously been explored near large-scale confining geometries and in the presence of very small-scale surface roughness. Here, we consider an intermediate case of how a simple microswimmer, the tangential spherical squirmer, behaves adjacent to singly and doubly periodic sinusoidal surface topographies that spatially oscillate with an amplitude that is an order of magnitude less than the swimmer size and wavelengths that are also within an order of magnitude of this scale. The nearest neighbour regularised Stokeslet method is used for numerical explorations after validating its accuracy for a spherical tangential squirmer that swims stably near a flat surface. The same squirmer is then introduced to different surface topographies. The key governing factor in the resulting swimming behaviour is the size of the squirmer relative to the surface topography wavelength. For instance, directional guidance is not observed when the squirmer is much larger, or much smaller, than the surface topography wavelength. In contrast, once the squirmer size is on the scale of the topography wavelength, limited guidance is possible, often with local capture in the topography troughs. However, complex dynamics can also emerge, especially when the initial configuration is not close to alignment along topography troughs or above topography crests. In contrast to sensitivity in alignment and topography wavelength, reductions in the amplitude of the surface topography or variations in the shape of the periodic surface topography do not have extensive impacts on the squirmer behaviour. Our findings more generally highlight that the numerical framework provides an essential basis to elucidate how swimmers may be guided by surface topography.

Hydrodynamic effects on the liquid-hexatic transition of active colloids

We study numerically the role of hydrodynamics in the liquid-hexatic transition of active colloids at intermediate activity, where motility induced phase separation (MIPS) does not occur. We show that in the case of active Brownian particles (ABP), the critical density of the transition decreases upon increasing the particle's mass, enhancing ordering, while self-propulsion has the opposite effect in the activity regime considered. Active hydrodynamic particles (AHP), instead, undergo the liquid-hexatic transition at higher values of packing fraction [Formula: see text] than the corresponding ABP, suggesting that hydrodynamics have the net effect of disordering the system. At increasing densities, close to the hexatic-liquid transition, we found in the case of AHP the appearance of self-sustained organized motion with clusters of particles moving coherently.

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.

Bacterial active matter

Bacteria are among the oldest and most abundant species on Earth. Bacteria successfully colonize diverse habitats and play a significant role in the oxygen, carbon, and nitrogen cycles. They also form human and animal microbiota and may become sources of pathogens and a cause of many infectious diseases. Suspensions of motile bacteria constitute one of the most studied examples of active matter: a broad class of non-equilibrium systems converting energy from the environment (e.g., chemical energy of the nutrient) into mechanical motion. Concentrated bacterial suspensions, often termed active fluids, exhibit complex collective behavior, such as large-scale turbulent-like motion (so-called bacterial turbulence) and swarming. The activity of bacteria also affects the effective viscosity and diffusivity of the suspension. This work reports on the progress in bacterial active matter from the physics viewpoint. It covers the key experimental results, provides a critical assessment of major theoretical approaches, and addresses the effects of visco-elasticity, liquid crystallinity, and external confinement on collective behavior in bacterial suspensions.

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.

Collective motion of cells modeled as ring polymers

In this article, we use a coarse-grained model of disjoint semi-flexible ring polymers to investigate computationally the spatiotemporal collective behavior of cell colonies. A ring polymer in this model is self-propelled by a motility force along the cell's polarity, which depends on its historical kinetics. Despite the repulsive interaction between the cells, a collective behavior sets in as a result of cells pushing against each other. This cooperative motion emerges as the amplitude of the motility force is increased and/or their areal density is increased. The degree of collectivity, characterized by the average cluster size, the velocity field order parameter, and the polarity field nematic order parameter, is found to increase with increasing the amplitude of the motility force and area coverage of the cells. Furthermore, the degree of alignment exhibited by the cell velocity field within a cluster is found to be stronger than that exhibited by the cell polarity. Comparison between the collective behavior of elongated cells and that of circular cells, at the same area coverage and motility force, shows that elongated cells exhibit a stronger collective behavior than circular cells, in agreement with earlier studies of self-propelled anisotropic particles. An investigation of two-cell collisions shows that while two clustered cells move in tandem, their polarities are misaligned. As such the cells push against each other while moving coherently.

Active particles in explicit solvent: Dynamics of clustering for alignment interaction

We study the dynamics of clustering in systems containing active particles that are immersed in an explicit solvent. For this, we have adopted a hybrid simulation method, consisting of molecular dynamics and multiparticle collision dynamics. In our model, the overlap-avoiding passive interaction of an active particle with another active particle or a solvent particle has been taken care of via variants of the Lennard-Jones potential. Dynamic interactions among the active particles have been incorporated via a Vicsek-like alignment rule in self-propulsion that facilitates clustering. We quantify the effects of activity and importance of hydrodynamics on the dynamics of clustering via variations of relevant system parameters. Results are obtained for low overall density of active particles, for which the state point is close to the vapor branch of the coexistence curve, and thus the morphology consists of disconnected clusters. In such a situation, the mechanism of growth switches among particle diffusion, diffusive coalescence, and ballistic aggregation, depending upon the presence or absence of active and hydrodynamic interactions providing different kinds of mobilities to the clusters. Corresponding growth laws have been quantified and discussed in the context of appropriate theoretical pictures. Our results suggest that multiparticle collision dynamics is an effective method for the investigation of hydrodynamic phenomena in phase-separating active matter systems.

Migration of active filaments under Poiseuille flow in a microcapillary tube

We present a comprehensive study of active filaments confined in a cylindrical channel under Poiseuille flow. The activity drives the filament towards the channel boundary, whereas external fluid flow migrates the filament away from the boundary. This migration further shifts towards the centre for higher flow strength. The migration behaviour of the filaments is presented in terms of the alignment order parameter that shows the alignment grows with shear and activity. Further, we have also addressed the role of length of filament on the migration behaviour, which suggests higher migration for larger filaments. Moreover, we discuss the polar ordering of filaments as a function of distance from the centre of channel that displays upstream motion near the boundary and downstream motion at the centre of the tube.

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.

Can playing Spirograph lead to an ordered structure in self-propelled particles?

Understanding the local dynamics of microorganisms infecting a cell could help us develop efficient strategies to counter their aggregation. In the present study we have introduced a simple model of self-propelled particles (SPPs) with constant linear velocity, both in 2 and 3 dimensions, which captures the essential features of a microorganism's aggregation as well the dynamics around an attractive point (AP). The static behavior shows the presence of an icosahedral structure for a finite number of SPPs, and a hexagonal closed packed structure for an infinite number of SPPs, which was confirmed using Steinhardt bond order parameters for a 3-dimensional model. For a single SPP the dynamic behaviour involves the formation of orbits around the AP, which can be categorised into three dynamical regions based on the strength of coupling between the AP and SPP. For weak coupling we observe a rosette-like trajectory reminiscent of the pattern formed by the Spirograph toy. For intermediate coupling, circular trajectories were observed, and for very strong coupling the SPP was static and was always aligned with the AP. The radial distance from the AP to SPP was determined by the angular velocities of the SPP for the rosette-like region whereas for the circular and static regions, it was determined by the coupling constant. Even for a finite number of SPPs we observed the same behavior as long as the SPPs could rotate around the AP without colliding with each other.

Spontaneous propulsion of an isotropic colloid in a phase-separating environment

The motion of active colloids is generally achieved through their anisotropy, as exemplified by Janus colloids. Recently, there was a growing interest in the propulsion of isotropic colloids, which requires some local symmetry breaking. Although several mechanisms for such propulsion were proposed, little is known about the role played by the interactions within the environment of the colloid, which can have a dramatic effect on its propulsion. Here, we propose a minimal model of an isotropic colloid in a bath of solute particles that interact with each other. These interactions lead to a spontaneous phase transition close to the colloid, to directed motion of the colloid over very long timescales and to significantly enhanced diffusion, in spite of the crowding induced by solute particles. We determine the range of parameters where this effect is observable in the model, and we propose an effective Langevin equation that accounts for it and allows one to determine the different contributions at stake in self-propulsion and enhanced diffusion.

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.

Cyclic force driven colloidal self-assembly near a solid surface

HYPOTHESIS

Self-assembled colloidal mobility out of a non-equilibrium system can depend on many external and interparticle forces including hydrodynamic forces. While the driving forces guiding colloidal suspension, translation and self-assembly are different and unique, hydrodynamic forces are always present and can significantly influence particle motion. Unfortunately, these interparticle hydrodynamic interactions are typically overlooked.

EXPERIMENTS

Here, we studied the collective behavior of colloidal particles (4.0 µm PMMA), located near the solid surface in a fluid medium confined in a cylindrical cell (3.0 mm diameter, 0.25 mm height) which was rotated vertically at a low rotational speed (20 rpm). The observed colloidal behavior was then validated through a Stokesian dynamics simulation where the concept of hydrodynamic contact force or lubrication interactions are avoided which is not physically intuitive and mathematically cumbersome. Rather, we adopted hard-sphere like colloidal collision or mobility model, while adopting other useful simplification and approximations.

FINDINGS

Upon particles settling in a circular orbit, they hydrodynamically interact with each other and evolve in different structures depending on the pattern of gravity forces. Their agglomeration is a function of the applied rotation scheme, either forming colloidal clusters or lanes. While evolving into dynamic structures, colloids also laterally migrate away from the surface.

Crystallisation and polymorph selection in active Brownian particles

We explore crystallisation and polymorph selection in active Brownian particles with numerical simulation. In agreement with previous work (Wysocki et al. in Europhys Lett 105:48004, 2014), we find that crystallisation is suppressed by activity and occurs at higher densities with increasing Péclet number ([Formula: see text]). While the nucleation rate decreases with increasing activity, the crystal growth rate increases due to the accelerated dynamics in the melt. As a result of this competition, we observe the transition from a nucleation and growth regime at high [Formula: see text] to "spinodal nucleation" at low [Formula: see text]. Unlike the case of passive hard spheres, where preference for FCC over HCP polymorphs is weak, activity causes the annealing of HCP stacking faults, thus strongly favouring the FCC symmetry at high [Formula: see text]. When freezing occurs more slowly, in the nucleation and growth regime, this tendency is much reduced and we see a trend towards the passive case of little preference for either polymorph.

Hydrodynamics of immiscible binary fluids with viscosity contrast: a multiparticle collision dynamics approach

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

Collective Entrainment and Confinement Amplify Transport by Schooling Microswimmers

Microswimmers can serve as cargo carriers that move deep inside complex flow networks. When a school collectively entrains the surrounding fluid, their transport capacity can be enhanced. This effect is quantified with good agreement between experiments with self-propelled droplets and a confined Brinkman squirmer model. The volume of liquid entrained can be much larger than the droplet itself, amplifying the effective cargo capacity over an order of magnitude, even for dilute schools. Hence, biological and engineered swimmers can efficiently transport materials into confined environments.

A Review on the Some Issues of Multiphase Flow with Self-Driven Particles

Multiphase flow with self-driven particles is ubiquitous and complex. Exploring the flow properties has both important academic meaning and engineering value. This review emphasizes some recent studies on multiphase flow with self-driven particles: the hydrodynamic interactions between self-propelled/self-rotary particles and passive particles; the aggregation, phase separation and sedimentation of squirmers; the influence of rheological properties on its motion; and the kinematic characteristics of axisymmetric squirmers. Finally, some open problems, challenges, and future directions are highlighted.

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.

Polarization-density patterns of active particles in motility gradients

The colocalization of density modulations and particle polarization is a characteristic emergent feature of motile active matter in activity gradients. We employ the active-Brownian-particle model to derive precise analytical expressions for the density and polarization profiles of a single Janus-type swimmer in the vicinity of an abrupt activity step. Our analysis allows for an optional (but not necessary) orientation-dependent propulsion speed, as often employed in force-free particle steering. The results agree well with measurement data for a thermophoretic microswimmer presented in the companion paper [Söker et al., Phys. Rev. Lett. 126, 228001 (2021)10.1103/PhysRevLett.126.228001], and they can serve as a template for more complex applications, e.g., to motility-induced phase separation or studies of physical boundaries. The essential physics behind our formal results is robustly captured and elucidated by a schematic two-species "run-and-tumble" model.

Squirming in a viscous fluid enclosed by a Brinkman medium

Cell motility plays important roles in a range of biological processes, such as reproduction and infections. Studies have hypothesized that the ulcer-causing bacterium Helicobacter pylori invades the gastric mucus layer lining the stomach by locally turning nearby gel into sol, thereby enhancing its locomotion through the biological barrier. In this work, we present a minimal theoretical model to investigate how heterogeneity created by a swimmer affects its own locomotion. As a generic locomotion model, we consider the swimming of a spherical squirmer in a purely viscous fluid pocket (representing the liquified or degelled region) surrounded by a Brinkman porous medium (representing the mucus gel). The use of the squirmer model enables an exact, analytical solution to this hydrodynamic problem. We obtain analytical expressions for the swimming speed, flow field, and power dissipation of the swimmer. Depending on the details of surface velocities and fluid properties, our results reveal the existence of a minimum threshold size of mucus gel that a swimmer needs to liquify in order to gain any enhancement in swimming speed. The threshold size can be as much as approximately 30% of the swimmer size. We contrast these predictions with results from previous models and highlight the significant role played by the details of surface actuations. In addition to their biological implications, these results could also inform the design of artificial microswimmers that can penetrate into biological gels for more effective drug delivery.

Smoothed profile method for direct numerical simulations of hydrodynamically interacting particles

A general method is presented for computing the motions of hydrodynamically interacting particles in various kinds of host fluids for arbitrary Reynolds numbers. The method follows the standard procedure for performing direct numerical simulations (DNS) of particulate systems, where the Navier-Stokes equation must be solved consistently with the motion of the rigid particles, which defines the temporal boundary conditions to be satisfied by the Navier-Stokes equation. The smoothed profile (SP) method provides an efficient numerical scheme for coupling the continuum fluid mechanics with the dispersed moving particles, which are allowed to have arbitrary shapes. In this method, the sharp boundaries between solid particles and the host fluid are replaced with a smeared out thin shell (interfacial) region, which can be accurately resolved on a fixed Cartesian grid utilizing a SP function with a finite thickness. The accuracy of the SP method is illustrated by comparison with known exact results. In the present paper, the high degree of versatility of the SP method is demonstrated by considering several types of active and passive particle suspensions.

Hydrodynamic collision between a microswimmer and a passive particle in a micro-channel

Microswimmers interacting with passive particles in confinement are common in many systems, e.g., spermatozoa encountering other cells or debris in the female reproductive tract or active particles interacting with polymers and tracers in microfluidic channels. The behaviour of such systems is driven by simultaneous, three way hydrodynamic interactions between the microswimmer, the passive particle and the microchannel walls. Therefore, in this work we investigate the hydrodynamic collision between a model microswimmer and a passive particle using three different methods: (i) the point particle approach, (ii) analytical calculations based on method of reflections, and (iii) lattice Boltzmann numerical simulations. We show that the hydrodynamic collision is essentially an asymmetric process - the trajectory of the microswimmer is altered only in an intermediate stage while the passive particle undergoes a three stage displacement with a net displacement towards or away from the microchannel walls. The path of the passive particle is a simple consequence of the velocity field generated by the swimmer: an open triangle in bulk fluid and a loop-like trajectory in confinement. We demonstrate the generality of our findings and conclude that the net displacement of the passive particle due to collision may be capitalised in order to develop applications such as size separation of colloidal particles and deposition of particles in the microchannel interiors.

Phase separation of an active colloidal suspension via quorum-sensing

We present the Brownian dynamics simulation of an active colloidal suspension in two dimensions, where the self-propulsion speed of a colloid is regulated according to the local density sensed by it. The role of concentration-dependent motility in the phase-separation of colloids and their dynamics is investigated in detail. Interestingly, the system phase separates at a very low packing fraction (Φ≈ 0.125) at higher self-propulsion speeds (Pe), into a dense phase coexisting with a homogeneous phase and attains a long-range crystalline order beyond the transition point. The transition point is quantified here from the local density profiles and local and global-bond order parameters. We have shown that the characteristics of the phase diagram are qualitatively akin to the active Brownian particle (ABP) model. Moreover, our investigation reveals that the density-dependent motility amplifies the slow-down of the directed speed, which facilitates phase-separation even at low packing fractions. The effective diffusivity shows a crossover from quadratic rise to a power-law behavior of exponent 3/2 with Pe in the phase-separated regime. Furthermore, we have shown that the effective diffusion decreases exponentially with packing fraction in the phase-separated regime, while it shows a linear decrease in the single phase regime.

Controlled locomotion of a droplet propelled by an encapsulated squirmer

We work out the propulsion of a viscous drop which is driven by two mechanisms: the active velocity of an encapsulated squirmer and an externally applied force acting on the squirmer. Of particular interest is the existence of a stable comoving state of drop and squirmer, allowing for controlled manipulation of the viscous drop by external forcing. The velocities of droplet and squirmer, as well as the conditions for a stable comoving state are worked out analytically for the axisymmetric configuration with a general displacement of the squirmer from the center of the droplet.

FastTrack: An open-source software for tracking varying numbers of deformable objects

Analyzing the dynamical properties of mobile objects requires to extract trajectories from recordings, which is often done by tracking movies. We compiled a database of two-dimensional movies for very different biological and physical systems spanning a wide range of length scales and developed a general-purpose, optimized, open-source, cross-platform, easy to install and use, self-updating software called FastTrack. It can handle a changing number of deformable objects in a region of interest, and is particularly suitable for animal and cell tracking in two-dimensions. Furthermore, we introduce the probability of incursions as a new measure of a movie’s trackability that doesn’t require the knowledge of ground truth trajectories, since it is resilient to small amounts of errors and can be computed on the basis of an ad hoc tracking. We also leveraged the versatility and speed of FastTrack to implement an iterative algorithm determining a set of nearly-optimized tracking parameters—yet further reducing the amount of human intervention—and demonstrate that FastTrack can be used to explore the space of tracking parameters to optimize the number of swaps for a batch of similar movies. A benchmark shows that FastTrack is orders of magnitude faster than state-of-the-art tracking algorithms, with a comparable tracking accuracy. The source code is available under the GNU GPLv3 at https://github.com/FastTrackOrg/FastTrack and pre-compiled binaries for Windows, Mac and Linux are available at http://www.fasttrack.sh.

Many researchers and engineers face the challenge of tracking objects from very different systems across several fields of research. We observed that despite this diversity the core of the tracking task is very general and can be formalized. We thus introduce the notion of incursions—i.e. to what extent an object can enter a neighbor’s space—which can be defined on a statistical basis and captures the interplay between the acquisition rate, the objects’ dynamics and the geometrical characteristics of the scene, including density. To validate this approach, we compiled a dataset from various fields of Physics, Biology and human activities to serve as a benchmark for general-purpose tracking softwares. This dataset is open and accepts new submissions. We also developped a software called FastTrack that is able to track most of the movies in the dataset by proposing standard image processing tools and state-of-the-art implementation of the matching algorithm, which is at the core of the tracking task. Besides, it is open-source, simple to install and use and has an ergonomic interface to obtain fast and reliable results. FastTrack is particularly convenient for small-scale research projects, typically when the development of a dedicated software is overkill.

Disorder-Induced Long-Ranged Correlations in Scalar Active Matter

We study the impact of quenched random potentials and torques on scalar active matter. Microscopic simulations reveal that motility-induced phase separation is replaced in two dimensions by an asymptotically homogeneous phase with anomalous long-ranged correlations and nonvanishing steady-state currents. Using a combination of phenomenological models and a field-theoretical treatment, we show the existence of a lower-critical dimension d_{c}=4, below which phase separation is only observed for systems smaller than an Imry-Ma length scale. We identify a weak-disorder regime in which the structure factor scales as S(q)∼1/q^{2}, which accounts for our numerics. In d=2, we predict that, at larger scales, the behavior should cross over to a strong-disorder regime. In d>2, these two regimes exist separately, depending on the strength of the potential.

Multi-particle collision dynamics with a non-ideal equation of state. I

The method of multi-particle collision dynamics (MPCD) and its different implementations are commonly used in the field of soft matter physics to simulate fluid flow at the micron scale. Typically, the coarse-grained fluid particles are described by the equation of state of an ideal gas, and the fluid is rather compressible. This is in contrast to conventional fluids, which are incompressible for velocities much below the speed of sound, and can cause inhomogeneities in density. We propose an algorithm for MPCD with a modified collision rule that results in a non-ideal equation of state and a significantly decreased compressibility. It allows simulations at less computational costs compared to conventional MPCD algorithms. We derive analytic expressions for the equation of state and the corresponding compressibility as well as shear viscosity. They show overall very good agreement with simulations, where we determine the pressure by simulating a quiet bulk fluid and the shear viscosity by simulating a linear shear flow and a Poiseuille flow.

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.

Bacterial biomechanics-From individual behaviors to biofilm and the gut flora

Bacteria inhabit a variety of locations and play important roles in the environment and health. Our understanding of bacterial biomechanics has improved markedly in the last decade and has revealed that biomechanics play a significant role in microbial biology. The obtained knowledge has enabled investigation of complex phenomena, such as biofilm formation and the dynamics of the gut flora. A bottom-up strategy, i.e., from the cellular to the macroscale, facilitates understanding of macroscopic bacterial phenomena. In this Review, we first cover the biomechanics of individual bacteria in the bulk liquid and on surfaces as the base of complex phenomena. The collective behaviors of bacteria in simple environments are next introduced. We then introduce recent advances in biofilm biomechanics, in which adhesion force and the flow environment play crucial roles. We also review transport phenomena in the intestine and the dynamics of the gut flora, focusing on that in zebrafish. Finally, we provide an overview of the future prospects for the field.

Unravelling the role of phoretic and hydrodynamic interactions in active colloidal suspensions

Active fluids comprise a variety of systems composed of elements immersed in a fluid environment which can convert some form of energy into directed motion; as such they are intrinsically out-of-equilibrium in the absence of any external force. A fundamental problem in the physics of active matter concerns the understanding of how the characteristics of autonomous propulsion and agent-agent interactions determine the collective dynamics of the system. We study numerically the suspensions of self-propelled diffusiophoretic colloids, in (quasi)-2d configurations, accounting for both dynamically resolved solute-mediated phoretic interactions and solvent-mediated hydrodynamic interactions. Our results show that the system displays different scenarios at changing the colloid-solute affinity and it develops a cluster phase in the chemoattractive case. We study the statistics of cluster sizes and cluster morphologies for different magnitudes of colloidal activity. Finally, we provide evidences that hydrodynamics plays a relevant role in the aggregation kinetics and cluster morphology, significantly hindering cluster growth.

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.

Competing active and passive interactions drive amoebalike crystallites and ordered bands in active colloids

Swimmers and self-propelled particles are physical models for the collective behavior and motility of a wide variety of living systems, such as bacteria colonies, bird flocks, and fish schools. Such artificial active materials are amenable to physical models which reveal the microscopic mechanisms underlying the collective behavior. Here we study colloids in a dc electric field. Our quasi-two-dimensional system of electrically driven particles exhibits a rich and exotic phase behavior exhibiting passive crystallites, motile crystallites, an active gas, and banding. Amongst these are two mesophases, reminiscent of systems with competing interactions. At low field strengths activity suppresses demixing, leading to motile crystallites. Meanwhile, at high field strengths, activity drives partial demixing to traveling bands. We parametrize a particulate simulation model which reproduces the experimentally observed phases.