Giant adsorption of microswimmers: Duality of shape asymmetry and wall curvature

Active transport of a passive colloid in a bath of run-and-tumble particles

The dispersion of a passive colloid immersed in a bath of non-interacting and non-Brownian run-and-tumble microswimmers in two dimensions is analyzed using stochastic simulations and an asymptotic theory, both based on a minimal model of swimmer-colloid collisions characterized solely by frictionless steric interactions. We estimate the effective long-time diffusivity D of the suspended colloid resulting from its interaction with the active bath, and elucidate its dependence on the level of activity (persistence length of swimmer trajectories), the mobility ratio of the colloid to a swimmer, and the number density of swimmers in the bath. We also propose a semi-analytical model for the colloid diffusivity in terms of the variance and correlation time of the net fluctuating active force on the colloid resulting from swimmer collisions. Quantitative agreement is found between numerical simulations and analytical results in the experimentally-relevant regime of low swimmer density, low mobility ratio, and high activity.

Phenotyping single-cell motility in microfluidic confinement

The movement trajectories of organisms serve as dynamic read-outs of their behaviour and physiology. For microorganisms this can be difficult to resolve due to their small size and fast movement. Here, we devise a novel droplet microfluidics assay to encapsulate single micron-sized algae inside closed arenas, enabling ultralong high-speed tracking of the same cell. Comparing two model species - Chlamydomonas reinhardtii (freshwater, 2 cilia), and Pyramimonas octopus (marine, 8 cilia), we detail their highly-stereotyped yet contrasting swimming behaviours and environmental interactions. By measuring the rates and probabilities with which cells transition between a trio of motility states (smooth-forward swimming, quiescence, tumbling or excitable backward swimming), we reconstruct the control network that underlies this gait switching dynamics. A simplified model of cell-roaming in circular confinement reproduces the observed long-term behaviours and spatial fluxes, including novel boundary circulation behaviour. Finally, we establish an assay in which pairs of droplets are fused on demand, one containing a trapped cell with another containing a chemical that perturbs cellular excitability, to reveal how aneural microorganisms adapt their locomotor patterns in real-time.

Interplay of surface interaction and magnetic torque in single-cell motion of magnetotactic bacteria in microfluidic confinement

Swimming microorganisms often experience complex environments in their natural habitat. The same is true for microswimmers in envisioned biomedical applications. The simple aqueous conditions typically studied in the lab differ strongly from those found in these environments and often exclude the effects of small volume confinement or the influence that external fields have on their motion. In this work, we investigate magnetically steerable microswimmers, specifically magnetotactic bacteria, in strong spatial confinement and under the influence of an external magnetic field. We trap single cells in micrometer-sized microfluidic chambers and track and analyze their motion, which shows a variety of different trajectories, depending on the chamber size and the strength of the magnetic field. Combining these experimental observations with simulations using a variant of an active Brownian particle model, we explain the variety of trajectories by the interplay between the wall interactions and the magnetic torque. We also analyze the pronounced cell-to-cell heterogeneity, which makes single-cell tracking essential for an understanding of the motility patterns. In this way, our work establishes a basis for the analysis and prediction of microswimmer motility in more complex environments.

Unraveling the Kinematics of Sperm Motion by Reconstructing the Flagellar Wave Motion in 3D

Sperm swim through the female reproductive tract by propagating a 3D flagellar wave that is self-regulatory in nature and driven by dynein motors. Traditional microscopy methods fail to capture the full dynamics of sperm flagellar activity as they only image and analyze sperm motility in 2D. Here, an automated platform to analyze sperm swimming behavior in 3D by using thin-lens approximation and high-speed dark field microscopy to reconstruct the flagellar waveform in 3D is presented. It is found that head-tethered mouse sperm exhibit a rolling beating behavior in 3D with the beating frequency of 6.2 Hz using spectral analysis. The flagellar waveform bends in 3D, particularly in the distal regions, but is only weakly nonplanar and ambidextrous in nature, with the local helicity along the flagellum fluctuating between clockwise and counterclockwise handedness. These findings suggest a nonpersistent flagellar helicity. This method provides new opportunities for the accurate measurement of the full motion of eukaryotic flagella and cilia which is essential for a biophysical understanding of their activation by dynein motors.

Noncentral forces mediated between two inclusions in a bath of active Brownian rods

Using Brownian Dynamics simulations, we study effective interactions mediated between two identical and impermeable disks (inclusions) immersed in a bath of identical, active (self-propelled), Brownian rods in two spatial dimensions, by assuming that the self-propulsion axis of the rods may generally deviate from their longitudinal axis. When the self-propulsion is transverse (perpendicular to the rod axis), the accumulation of active rods around the inclusions is significantly enhanced, causing a more expansive steric layering (ring formation) of the rods around the inclusions, as compared with the reference case of longitudinally self-propelling rods. As a result, the transversally self-propelling rods also mediate a significantly longer ranged effective interaction between the inclusions. The bath-mediated interaction arises due to the overlaps between the active-rod rings formed around the inclusions, as they are brought into small separations. When the self-propulsion axis is tilted relative to the rod axis, we find an asymmetric imbalance of active-rod accumulation around the inclusion dimer. This leads to a noncentral interaction, featuring an anti-parallel pair of transverse force components and, hence, a bath-mediated torque on the dimer.

Emergent probability fluxes in confined microbial navigation

When the motion of a motile cell is observed closely, it appears erratic, and yet the combination of nonequilibrium forces and surfaces can produce striking examples of organization in microbial systems. While most of our current understanding is based on bulk systems or idealized geometries, it remains elusive how and at which length scale self-organization emerges in complex geometries. Here, using experiments and analytical and numerical calculations, we study the motion of motile cells under controlled microfluidic conditions and demonstrate that probability flux loops organize active motion, even at the level of a single cell exploring an isolated compartment of nontrivial geometry. By accounting for the interplay of activity and interfacial forces, we find that the boundary's curvature determines the nonequilibrium probability fluxes of the motion. We theoretically predict a universal relation between fluxes and global geometric properties that is directly confirmed by experiments. Our findings open the possibility to decipher the most probable trajectories of motile cells and may enable the design of geometries guiding their time-averaged motion.

Active particles with polar alignment in ring-shaped confinement

We study steady-state properties of active, nonchiral and chiral Brownian particles with polar alignment and steric interactions confined within a ring-shaped confinement (annulus) in two dimensions. Exploring possible interplays between polar interparticle alignment, geometric confinement and the surface curvature, being incorporated here on minimal levels, we report a surface-population reversal effect, whereby active particles migrate from the outer concave boundary of the annulus to accumulate on its inner convex boundary. This contrasts the conventional picture, implying stronger accumulation of active particles on concave boundaries relative to the convex ones. The population reversal is caused by both particle alignment and surface curvature, disappearing when either of these factors is absent. We explore the ensuing consequences for the chirality-induced current and swim pressure of active particles and analyze possible roles of system parameters, such as the mean number density of particles and particle self-propulsion, chirality, and alignment strengths.

Aggregate morphology of active Brownian particles on porous, circular walls

We study the motility-induced aggregation of active Brownian particles (ABPs) on a porous, circular wall. We observe that the morphology of aggregated dense-phase on a static wall depends on the wall porosity, particle motility, and the radius of the circular wall. Our analysis reveals two morphologically distinct, dense aggregates; a connected dense cluster that spreads uniformly on the circular wall and a localized cluster that breaks the rotational symmetry of the system. These distinct morphological states are similar to the macroscopic structures observed in aggregates on planar, porous walls. We systematically analyze the parameter regimes where the different morphological states are observed. We further extend our analysis to motile circular rings. We show that the motile ring propels almost ballistically due to the force applied by the active particles when they form a localized cluster, whereas it moves diffusively when the active particles form a continuous cluster. This property demonstrates the possibility of extracting useful work from a system of ABPs, even without artificially breaking the rotational symmetry.

Vortex formation of spherical self-propelled particles around a circular obstacle

A vortex is a common ratchet phenomenon in active systems. The spatial symmetry is usually broken by introducing asymmetric shapes or spontaneously by collective motion in the presence of hydrodynamic interactions or other alignment effects. Unexpectedly, we observe, by simulations, the formation of a vortex in the simplest model of a circular obstacle immersed in a bath of spherical self-propelled particles. No symmetry-breaking factors mentioned above are included in this model. The vortex forms only when the particle activity is high, i.e. large persistence. The obstacle size is also a key factor and the vortex only forms in a limited range of obstacle sizes. The sustainment of the vortex originates from the bias of the rotating particle cluster around the obstacle in accepting the incoming particles based on their propelling directions. Our results provide new understanding of and insights into the spontaneous symmetry-breaking and ratchet phenomena in active matter.

Fine balance of chemotactic and hydrodynamic torques: When microswimmers orbit a pillar just once

We study the detention statistics of self-propelling droplet microswimmers attaching to microfluidic pillars. These droplets show negative autochemotaxis: they shed a persistent repulsive trail of spent fuel that biases them to detach from pillars in a specific size range after orbiting them just once. We have designed a microfluidic assay recording microswimmers in pillar arrays of varying diameter, derived detention statistics via digital image analysis, and interpreted these statistics via the Langevin dynamics of an active Brownian particle model. By comparing data from orbits with and without residual chemical field, we can independently estimate quantities such as hydrodynamic and chemorepulsive torques, chemical coupling constants and diffusion coefficients, as well as their dependence on environmental factors such as wall curvature. This type of analysis is generalizable to many kinds of microswimmers.

Tuning the self-organization of confined active particles by the steepness of the trap

We consider the collective dynamics of self-propelling particles in two dimensions. They can align themselves according to the direction of propulsion of their neighbours, together with small rotational fluctuations. They also interact with each other via soft, isotropic, repulsive potentials. The particles are confined in a circular trap. The steepness of the trap is tuneable. The average packing fraction of the particles is low. When the trap is steep, particles flock along its boundary. They form a polar cluster that spreads over the boundary. The cluster is not spatially ordered. We show that when the steepness is decreased beyond a threshold value, the cluster becomes round and compact and eventually spatial order (hexagonal) emerges in addition to the pre-established polar order. We investigate the kinetics of such ordering. We find that while rotating around the centre of the trap along its circular boundary, the cluster needs to roll around its centre of mass to be spatially ordered. We have studied the stability of the order when the trap is suddenly switched off. We find that for the particles with velocity alignment interaction, the decay of the spatial order is much slower than the particles without the alignment interaction.

Maximum in density heterogeneities of active swimmers

Suspensions of unicellular microswimmers such as flagellated bacteria or motile algae can exhibit spontaneous density heterogeneities at large enough concentrations. We introduce a novel model for biological microswimmers that creates the flow field of the corresponding microswimmers, and takes into account the shape anisotropy of the swimmer's body and stroke-averaged flagella. By employing multiparticle collision dynamics, we directly couple the swimmer's dynamics to the fluid's. We characterize the nonequilibrium phase diagram, as the filling fraction and Péclet number are varied, and find density heterogeneities in the distribution of both pullers and pushers, due to hydrodynamic instabilities. We find a maximum degree of clustering at intermediate filling fractions and at large Péclet numbers resulting from a competition of hydrodynamic and steric interactions between the swimmers. We develop an analytical theory that supports these results. This maximum might represent an optimum for the microorganisms' colonization of their environment.

Active fluids at circular boundaries: swim pressure and anomalous droplet ripening

We investigate the swim pressure exerted by non-chiral and chiral active particles on convex or concave circular boundaries. Active particles are modeled as non-interacting and non-aligning self-propelled Brownian particles. The convex and concave circular boundaries are used to model a fixed inclusion immersed in an active bath and a cavity (or container) enclosing the active particles, respectively. We first present a detailed analysis of the role of convex versus concave boundary curvature and of the chirality of active particles in their spatial distribution, chirality-induced currents, and the swim pressure they exert on the bounding surfaces. The results will then be used to predict the mechanical equilibria of suspended fluid enclosures (generically referred to as 'droplets') in a bulk with active particles being present either inside the bulk fluid or within the suspended droplets. We show that, while droplets containing active particles behave in accordance with standard capillary paradigms when suspended in a normal bulk, those containing a normal fluid exhibit anomalous behaviors when suspended in an active bulk. In the latter case, the excess swim pressure results in non-monotonic dependence of the inside droplet pressure on the droplet radius; hence, revealing an anomalous regime of behavior beyond a threshold radius, in which the inside droplet pressure increases upon increasing the droplet size. Furthermore, for two interconnected droplets, mechanical equilibrium can occur also when the droplets have different sizes. We thus identify a regime of anomalous droplet ripening, where two unequal-sized droplets can reach a final state of equal size upon interconnection, in stark contrast with the standard Ostwald ripening phenomenon, implying shrinkage of the smaller droplet in favor of the larger one.

Curvature-Guided Motility of Microalgae in Geometric Confinement

Microorganisms, such as bacteria and microalgae, often live in habitats consisting of a liquid phase and a plethora of interfaces. The precise ways in which these motile microbes behave in their confined environment remain unclear. Using experiments and Brownian dynamics simulations, we study the motility of a single Chlamydomonas microalga in an isolated microhabitat with controlled geometric properties. We demonstrate how the geometry of the habitat controls the cell's navigation in confinement. The probability of finding the cell swimming near the boundary increases with the wall curvature, as seen for both circular and elliptical chambers. The theory, utilizing an asymmetric dumbbell model of the cell and steric wall interactions, captures this curvature-guided navigation quantitatively with no free parameters.

Chemotactic droplet swimmers in complex geometries

Chemotaxis¹ and auto-chemotaxis are key mechanisms in the dynamics of micro-organisms, e.g. in the acquisition of nutrients and in the communication between individuals, influencing the collective behaviour. However, chemical signalling and the natural environment of biological swimmers are generally complex, making them hard to access analytically. We present a well-controlled, tunable artificial model to study chemotaxis and autochemotaxis in complex geometries, using microfluidic assays of self-propelling oil droplets in an aqueous surfactant solution (Herminghaus et al 2014 Soft Matter 10 7008-22; Krüger et al 2016 Phys. Rev. Lett. 117). Droplets propel via interfacial Marangoni stresses powered by micellar solubilisation. Moreover, filled micelles act as a chemical repellent by diffusive phoretic gradient forces. We have studied these chemotactic effects in a series of microfluidic geometries, as published in Jin et al (2017 Proc. Natl Acad. Sci. 114 5089-94): first, droplets are guided along the shortest path through a maze by surfactant diffusing into the maze from the exit. Second, we let auto-chemotactic droplet swimmers pass through bifurcating microfluidic channels and record anticorrelations between the branch choices of consecutive droplets. We present an analytical Langevin model matching the experimental data. In a previously unpublished experiment, pillar arrays of variable sizes and shapes provide a convex wall interacting with the swimmer and, in the case of attachment, bending its trajectory and forcing it to revert to its own trail. We observe different behaviours based on the interplay of wall curvature and negative autochemotaxis, i.e. no attachment for highly curved interfaces, stable trapping at large pillars, and a narrow transition region where negative autochemotaxis makes the swimmers detach after a single orbit.

Self-propulsion against a moving membrane: Enhanced accumulation and drag force

Self-propulsion (SP) is a main feature of active particles (AP), such as bacteria or biological micromotors, distinguishing them from passive colloids. A renowned consequence of SP is accumulation at static interfaces, even in the absence of hydrodynamic interactions. Here we address the role of SP in the interaction between AP and a moving semipermeable membrane. In particular, we implement a model of noninteracting AP in a channel crossed by a partially penetrable wall, moving at a constant velocity c. With respect to both the cases of passive colloids with c>0 and AP with c=0, the AP with finite c show enhancement of accumulation in front of the obstacle and experience a largely increased drag force. This effect is understood in terms of an effective potential localised at the interface between particles and membrane, of height proportional to cτ/ξ, where τ is the AP's reorientation time and ξ the width characterizing the surface's smoothness (ξ→0 for hard core obstacles). An approximate analytical scheme is able to reproduce the observed density profiles and the measured drag force, in very good agreement with numerical simulations. The effects discussed here can be exploited for automatic selection and filtering of AP with desired parameters.

Three-body correlations and conditional forces in suspensions of active hard disks

Self-propelled Brownian particles show rich out-of-equilibrium physics, for instance, the motility-induced phase separation (MIPS). While decades of studying the structure of liquids have established a deep understanding of passive systems, not much is known about correlations in active suspensions. In this work we derive an approximate analytic theory for three-body correlations and forces in systems of active Brownian disks starting from the many-body Smoluchowski equation. We use our theory to predict the conditional forces that act on a tagged particle and their dependence on the propulsion speed of self-propelled disks. We identify preferred directions of these forces in relation to the direction of propulsion and the positions of the surrounding particles. We further relate our theory to the effective swimming speed of the active disks, which is relevant for the physics of MIPS. To test and validate our theory, we additionally run particle-resolved computer simulations, for which we explicitly calculate the three-body forces. In this context, we discuss the modeling of active Brownian swimmers with nearly hard interaction potentials. We find very good agreement between our simulations and numerical solutions of our theory, especially for the nonequilibrium pair-distribution function. For our analytical results, we carefully discuss their range of validity in the context of the different levels of approximation we applied. This discussion allows us to study the individual contribution of particles to three-body forces and to the emerging structure. Thus, our work sheds light on the collective behavior, provides the basis for further studies of correlations in active suspensions, and makes a step towards an emerging liquid state theory.

Complex self-propelled rings: a minimal model for cell motility

Collective behavior of active matter is observed for self-propelled particles, such as vibrated disks and active Brownian particles, as well as for cytoskeletal filaments in motile cells. Here, a system of quasi two-dimensional penetrable self-propelled rods inside rigid rings is used to construct a complex self-propelled particle. The rods interact sterically with each other and with a stationary or mobile ring via a separation-shifted Lennard-Jones potential. They either have a sliding attachment to the inside of the ring at one of their ends, or can move freely within the ring confinement. We study the inner structure and dynamics of the mobile self-propelled rings. We find that these complex particles cannot only be characterized as active Brownian particles, but can also exhibit cell-like motility: random walks, persistent motion, circling, and run-and-circle motion.

Wetting Transitions Displayed by Persistent Active Particles

A lattice model for active matter is studied numerically, showing that it displays wetting transitions between three distinctive phases when in contact with an impenetrable wall. The particles in the model move persistently, tumbling with a small rate α, and interact via exclusion volume only. When increasing the tumbling rates α, the system transits from total wetting to partial wetting and unwetting phases. In the first phase, a wetting film covers the wall, with increasing heights when α is reduced. The second phase is characterized by wetting droplets on the wall with a periodic spacing between them. Finally, the wall dries with few particles in contact with it. These phases present nonequilibrium transitions. The first transition, from partial to total wetting, is continuous and the fraction of dry sites vanishes continuously when decreasing the tumbling rate α. For the second transition, from partial wetting to dry, the mean droplet distance diverges logarithmically when approaching the critical tumbling rate, with saturation due to finite-size effects.

Scattering of biflagellate microswimmers from surfaces

We use a three-bead-spring model to investigate the dynamics of biflagellate microswimmers near a surface. While the primary dynamics and scattering are governed by geometric-dependent direct contact, the fluid flows generated by the swimmer locomotion are important in orienting it toward or away from the surface. Flagellar noise and in particular cell spinning about the main axis help a surface-trapped swimmer escape, whereas the time a swimmer spends at the surface depends on the incident angle. The dynamics results from a nuanced interplay of direct collisions, hydrodynamics, noise, and the swimmer geometry. We show that to correctly capture the dynamics of a biflagellate swimmer, minimal models need to resolve the shape asymmetry.

Geometry-driven collective ordering of bacterial vortices

Controlling the phases of matter is a challenge that spans from condensed materials to biological systems. Here, by imposing a geometric boundary condition, we study the controlled collective motion of Escherichia coli bacteria. A circular microwell isolates a rectified vortex from disordered vortices masked in the bulk. For a doublet of microwells, two vortices emerge but their spinning directions show transition from parallel to anti-parallel. A Vicsek-like model for confined self-propelled particles gives the point where the two spinning patterns occur in equal probability and one geometric quantity governs the transition as seen in experiments. This mechanism shapes rich patterns including chiral configurations in a quadruplet of microwells, thus revealing a design principle of active vortices.

Ideal bulk pressure of active Brownian particles

The extent to which active matter might be described by effective equilibrium concepts like temperature and pressure is currently being discussed intensely. Here, we study the simplest model, an ideal gas of noninteracting active Brownian particles. While the mechanical pressure exerted onto confining walls has been linked to correlations between particles' positions and their orientations, we show that these correlations are entirely controlled by boundary effects. We also consider a definition of local pressure, which describes interparticle forces in terms of momentum exchange between different regions of the system. We present three pieces of analytical evidence which indicate that such a local pressure exists, and we show that its bulk value differs from the mechanical pressure exerted on the walls of the system. We attribute this difference to the fact that the local pressure in the bulk does not depend on boundary effects, contrary to the mechanical pressure. We carefully examine these boundary effects using a channel geometry, and we show a virial formula for the pressure correctly predicts the mechanical pressure even in finite channels. However, this result no longer holds in more complex geometries, as exemplified for a channel that includes circular obstacles.

Swim pressure on walls with curves and corners

The concept of swim pressure quantifies the average force exerted by microswimmers on confining walls in nonequilibrium. Here we explore how the swim pressure depends on the wall curvature and on the presence of sharp corners in the wall. For active Brownian particles at high dilution, we present a coherent framework which describes the force and torque on passive particles of arbitrary shape, in the limit of large particles compared to the persistence length of the swimmer trajectories. The resulting forces can be used to derive, for example, the activity-induced depletion interaction between two disks, as well as to optimize the shape of a tracer particle for high swimming velocity. Our predictions are verifiable in experiments on passive obstacles exposed to a bath of bacteria or artificial microswimmers.