Scattering of biflagellate microswimmers from surfaces

Trapping and scattering of a multiflagellated bacterium by a hard surface

Thiovulum majus, which is one of the fastest known bacteria, swims using hundreds of flagella. Unlike typical pusher cells, which swim in circular paths over hard surfaces, T. majus localize near hard boundaries by turning their flagella to exert a net force normal to the surface. To probe the torques that stabilize this hydrodynamically bound state, the trajectories of several thousand collisions between a T. majus cell and a wall of a quasi-two-dimensional microfluidic chamber are analyzed. Measuring the fraction of cells escaping the wall either to the left or to the right of the point of contact-and how this probability varies with incident angle and time spent in contact with the surface-maps the scattering dynamics onto a first passage problem. These measurements are compared to the prediction of a Fokker-Planck equation to fit the angular velocity of a cell in contact with a hard surface. This analysis reveals a bound state with a narrow basin of attraction in which cells orient their flagella normal to the surface. The escape angle predicted by matching these near field dynamics with the far-field hydrodynamics is consistent with observation. We discuss the significance of these results for the ecology of T. majus and their self-organization into active chiral crystals.

Dispersion of run-and-tumble microswimmers through disordered media

Understanding the transport properties of microorganisms and self-propelled particles in porous media has important implications for human health as well as microbial ecology. In free space, most microswimmers perform diffusive random walks as a result of the interplay of self-propulsion and orientation decorrelation mechanisms such as run-and-tumble dynamics or rotational diffusion. In an unstructured porous medium, collisions with the microstructure result in a decrease in the effective spatial diffusivity of the particles from its free-space value. Here, we analyze this problem for a simple model system consisting of noninteracting point particles performing run-and-tumble dynamics through a two-dimensional disordered medium composed of a random distribution of circular obstacles, in the absence of Brownian diffusion or hydrodynamic interactions. The particles are assumed to collide with the obstacles as hard spheres and subsequently slide on the obstacle surface with no frictional resistance while maintaining their orientation, until they either escape or tumble. We show that the variations in the long-time diffusivity can be described by a universal dimensionless hindrance function f(ϕ,Pe) of the obstacle area fraction ϕ and Péclet number Pe, or ratio of the swimmer run length to the obstacle size. We analytically derive an asymptotic expression for the hindrance function valid for dilute media (Peϕ≪1), and its extension to denser media is obtained using stochastic simulations. As we explain, the model is also easily generalized to describe dispersion in three dimensions.

The three-dimensional coarse-graining formulation of interacting elastohydrodynamic filaments and multi-body microhydrodynamics

Elastic filaments are vital to biological, physical and engineering systems, from cilia driving fluid in the lungs to artificial swimmers and micro-robotics. Simulating slender structures requires intricate balance of elastic, body, active and hydrodynamic moments, all in three dimensions. Here, we present a generalized three-dimensional (3D) coarse-graining formulation that is efficient, simple-to-implement, readily extendable and usable for a wide array of applications. Our method allows for simulation of collections of 3D elastic filaments, capable of full flexural and torsional deformations, coupled non-locally via hydrodynamic interactions, and including multi-body microhydrodynamics of structures with arbitrary geometry. The method exploits the exponential mapping of quaternions for tracking 3D rotations of each interacting element in the system, allowing for computation times up to 150 times faster than a direct quaternion implementation. Spheres are used as a 'building block' of both filaments and solid microstructures for straightforward and intuitive construction of arbitrary three-dimensional geometries present in the environment. We highlight the strengths of the method in a series of non-trivial applications including bi-flagellated swimming, sperm-egg scattering and particle transport by cilia arrays. Applications to lab-on-a-chip devices, multi-filaments, mono-to-multi flagellated microorganisms, Brownian polymers, and micro-robotics are straightforward. A Matlab code is provided for further customization and generalizations.

Ameboid cell migration through regular arrays of micropillars under confinement

Migrating cells often encounter a wide variety of topographic features-including the presence of obstacles-when navigating through crowded biological environments. Unraveling the impact of topography and crowding on the dynamics of cells is key to better understand many essential physiological processes such as the immune response. We study the impact of geometrical cues on ameboid migration of HL-60 cells differentiated into neutrophils. A microfluidic device is designed to track the cells in confining geometries between two parallel plates with distance h, in which identical micropillars are arranged in regular pillar forests with pillar spacing e. We observe that the cells are temporarily captured near pillars, with a mean contact time that is independent of h and e. By decreasing the vertical confinement h, we find that the cell velocity is not affected, while the persistence reduces; thus, cells are able to preserve their velocity when highly squeezed but lose the ability to control their direction of motion. At a given h, we show that by decreasing the pillar spacing e in the weak lateral confinement regime, the mean escape time of cells from effective local traps between neighboring pillars grows. This effect, together with the increase of cell-pillar contact frequency, leads to the reduction of diffusion constant D. By disentangling the contributions of these two effects on D in numerical simulations, we verify that the impact of cell-pillar contacts on cell diffusivity is more pronounced at smaller pillar spacing.

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.

Optimal motion control of three-sphere based low-Reynolds number swimming microrobot

Microrobots with their promising applications are attracting a lot of attention currently. A microrobot with a triangular mechanism was previously proposed by scientists to overcome the motion limitations in a low-Reynolds number flow; however, the control of this swimmer for performing desired manoeuvres has not been studied yet. Here, we have proposed some strategies for controlling its position. Considering the constraints on arm lengths, we proposed an optimal controller based on quadratic programming. The simulation results demonstrate that the proposed optimal controller can steer the microrobot along the desired trajectory as well as minimize fluctuations of the actuators length.

Microswimmers near corrugated, periodic surfaces

We explore hydrodynamic interactions between microswimmers and corrugated, or rough, surfaces, as found often in biological systems and microfluidic devices. Using the Lorentz reciprocal theorem for viscous flows we derive exact expressions for the roughness-induced velocities up to first order in the surface-height fluctuations and provide solutions for the translational and angular velocities valid for arbitrary surface shapes. We apply our theoretical predictions to elucidate the impact of a periodic, wavy surface on the velocities of microswimmers modeled in terms of a superposition of Stokes singularities. Our findings, valid in the framework of a far-field analysis, show that the roughness-induced velocities vary non-monotonically with the wavelength of the surface. For wavelengths comparable to the swimmer-surface distance a pusher can experience a repulsive contribution due to the reflection of flow fields at the edge of a surface cavity, which decreases the overall attraction to the wall.

Steric scattering of rod-like swimmers in low Reynolds number environments

Microbes form integral components of all natural ecosystems. In most cases, the surrounding micro-environment has physical variations that affect the movements of micro-swimmers, including solid objects of varying size, shape and density. As swimmers move through viscous environments, a combination of hydrodynamic and steric forces are known to significantly alter their trajectories in a way that depends on surface curvature. In this work, our goal was to clarify the role of steric forces when rod-like swimmers interact with solid objects comparable to cell size. We imaged hundreds-of-thousands of scattering interactions between swimming bacteria and micro-fabricated pillars with radii from ∼1 to ∼10 cell lengths. Scattering interactions were parameterized by the angle of the cell upon contact with the pillar, and primarily produced forward-scattering events that fell into distinct chiral distributions for scattering angle - no hydrodynamic trapping was observed. The chirality of a scattering event was a stochastic variable whose probability smoothly and symmetrically depended on the contact angle. Neglecting hydrodynamics, we developed a model that only considers contact forces and torques for a rear-pushed thin-rod scattering from a cylinder - the model predictions were in good agreement with measured data. Our results suggest that alteration of bacterial trajectories is subject to distinct mechanisms when interacting with objects of different size; primarily steric for objects below ∼10 cell lengths and requiring incorporation of hydrodynamics at larger scales. These results contribute to a mechanistic framework in which to examine (and potentially engineer) microbial movements through natural and synthetic environments that present complex steric structure.

Control of Helical Navigation by Three-Dimensional Flagellar Beating

Helical swimming is a ubiquitous strategy for motile cells to generate self-gradients for environmental sensing. The model biflagellate Chlamydomonas reinhardtii rotates at a constant 1-2 Hz as it swims, but the mechanism is unclear. Here, we show unequivocally that the rolling motion derives from a persistent, nonplanar flagellar beat pattern. This is revealed by high-speed imaging and micromanipulation of live cells. We construct a fully 3D model to relate flagellar beating directly to the free-swimming trajectories. For realistic geometries, the model reproduces both the sense and magnitude of the axial rotation of live cells. We show that helical swimming requires further symmetry breaking between the two flagella. These functional differences underlie all tactic responses, particularly phototaxis. We propose a control strategy by which cells steer toward or away from light by modulating the sign of biflagellar dominance.

Green algae scatter off sharp viscosity gradients

We study the behaviour of the green alga Chlamydomonas reinhardtii (CR) in the presence of neighbouring regions of different viscosity. We show that the velocity and angular diffusion of the algae decreases when the viscosity of the surrounding medium is increased. We report on a phenomenon occurring when the algae try to cross from a region of low viscosity to a highly viscous one, which causes CR to re-orient and scatter away from the interface if it is approached at a sufficiently small angle. We highlight that the effect does not occur for CR crossing from high to low viscosity regions. Lastly we show that algae do not concentrate in the region of high viscosity despite them swimming slower there. On the contrary, they concentrate in the region of low viscosity or maintain a uniform concentration profile, depending on the viscosity ratio between the two regions.

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.

Motile Bacteria at Oil-Water Interfaces: Pseudomonas aeruginosa

Bacteria are important examples of active or self-propelled colloids. Because of their directed motion, they accumulate near interfaces. There, they can become trapped and swim adjacent to the interface via hydrodynamic interactions, or they can adsorb directly and swim in an adhered state with complex trajectories that differ from those in bulk in both form and spatiotemporal implications. We have adopted the monotrichous bacterium Pseudomonas aeruginosa PA01 as a model species and have studied its motion at oil-aqueous interfaces. We have identified conditions in which bacteria swim persistently without restructuring the interface, allowing detailed and prolonged study of their motion. In addition to characterizing the ensemble behavior of the bacteria, we have observed a gallery of distinct trajectories of individual swimmers on and near fluid interfaces. We attribute these diverse swimming behaviors to differing trapped states for the bacteria in the fluid interface. These trajectory types include Brownian diffusive paths for passive adsorbed bacteria, curvilinear trajectories including curly paths with radii of curvature larger than the cell body length, and rapid pirouette motions with radii of curvature comparable to the cell body length. Finally, we see interfacial visitors that come and go from the interfacial plane. We characterize these individual swimmer motions. This work may impact nutrient cycles for bacteria on or near interfaces in nature. This work will also have implications in microrobotics, as active colloids in general and bacteria in particular are used to carry cargo in this burgeoning field. Finally, these results have implications in engineering of active surfaces that exploit interfacially trapped self-propelled colloids.

Medical micro/nanorobots in complex media

Medical micro/nanorobots have received tremendous attention over the past decades owing to their potential to be navigated into hard-to-reach tissues for a number of biomedical applications ranging from targeted drug/gene delivery, bio-isolation, detoxification, to nanosurgery. Despite the great promise, the majority of the past demonstrations are primarily under benchtop or in vitro conditions. Many developed micro/nanoscale propulsion mechanisms are based on the assumption of a homogeneous, Newtonian environment, while realistic biological environments are substantially more complex. Moving toward practical medical use, the field of micro/nanorobotics must overcome several major challenges including propulsion through complex media (such as blood, mucus, and vitreous) as well as deep tissue imaging and control in vivo. In this review article, we summarize the recent research efforts on investigating how various complexities in biological environments impact the propulsion of micro/nanoswimmers. We also highlight the emerging technological approaches to enhance the locomotion of micro/nanorobots in complex environments. The recent demonstrations of in vivo imaging, control and therapeutic medical applications of such micro/nanorobots are introduced. We envision that continuing materials and technological innovations through interdisciplinary collaborative efforts can bring us steps closer to the fantasy of "swallowing a surgeon".

Confinement and activity regulate bacterial motion in porous media

Understanding how bacteria move in porous media is critical to applications in healthcare, agriculture, environmental remediation, and chemical sensing. Recent work has demonstrated that E. coli, which moves by run-and-tumble dynamics in a homogeneous medium, exhibits a new form of motility when confined in a disordered porous medium: hopping-and-trapping motility, in which cells perform rapid, directed hops punctuated by intervals of slow, undirected trapping. Here, we use direct visualization to shed light on how these processes depend on pore-scale confinement and cellular activity. We find that hopping is determined by pore-scale confinement, and is independent of cellular activity; by contrast, trapping is determined by the competition between pore-scale confinement and cellular activity, as predicted by an entropic trapping model. These results thus help to elucidate the factors that regulate bacterial motion in porous media, and could help aid the development of new models of motility in heterogeneous environments.

Relating Rheotaxis and Hydrodynamic Actuation using Asymmetric Gold-Platinum Phoretic Rods

We explore the behavior of micron-scale autophoretic Janus (Au/Pt) rods, having various Au/Pt length ratios, swimming near a wall in an imposed background flow. We find that their ability to robustly orient and move upstream, i.e., to rheotax, depends strongly on the Au/Pt ratio, which is easily tunable in synthesis. Numerical simulations of swimming rods actuated by a surface slip show a similar rheotactic tunability when varying the location of the surface slip versus surface drag. The slip location determines whether swimmers are pushers (rear actuated), pullers (front actuated), or in between. Our simulations and modeling show that pullers rheotax most robustly due to their larger tilt angle to the wall, which makes them responsive to flow gradients. Thus, rheotactic response infers the nature of difficult to measure flow fields of an active particle, establishes its dependence on swimmer type, and shows how Janus rods can be tuned for flow responsiveness.

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

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

Bacterial hopping and trapping in porous media

Diverse processes-e.g. bioremediation, biofertilization, and microbial drug delivery-rely on bacterial migration in disordered, three-dimensional (3D) porous media. However, how pore-scale confinement alters bacterial motility is unknown due to the opacity of typical 3D media. As a result, models of migration are limited and often employ ad hoc assumptions. Here we reveal that the paradigm of run-and-tumble motility is dramatically altered in a porous medium. By directly visualizing individual Escherichia coli, we find that the cells are intermittently and transiently trapped as they navigate the pore space, exhibiting diffusive behavior at long time scales. The trapping durations and the lengths of "hops" between traps are broadly distributed, reminiscent of transport in diverse other disordered systems; nevertheless, we show that these quantities can together predict the long-time bacterial translational diffusivity. Our work thus provides a revised picture of bacterial motility in complex media and yields principles for predicting cellular migration.

Diffusion of active particles in a complex environment: Role of surface scattering

Experiments have shown that self-propelled particles can slide along the surface of a circular obstacle without becoming trapped over long times. Using simulations and theory, we study the impact of boundary conditions on the diffusive transport of active particles in an obstacle lattice. We find that particle dynamics with sliding boundary conditions result in large diffusivities even at high obstacle density, unlike classical specular reflection. These dynamics are very well described by a model based on run-and-tumble particles with microscopically derived reorientation functions arising from obstacle-induced tumbles. This model, however, fails to describe fine structure in the diffusivity at high obstacle density predicted by simulations for pusherlike collisions. Using a simple deterministic model, we show that this structure results from particles being guided by the lattice. Our results thus show how nonclassical surface scattering introduces a dependence on the lattice geometry at high densities. We discuss implications for the study of bacteria in complex environments.

Flow characteristics of Chlamydomonas result in purely hydrodynamic scattering

It has long been believed that eukaryotic flagellated swimming cells feel solid boundaries through direct ciliary contact. Specifically, based on observations of behavior of green alga Chlamydomonas reinhardtii it has been reported that it is their "flagella [that] prevent the cell body from touching the surface" [Kantsler et al., Proc. Natl. Acad. Sci. USA 110, 1187 (2013)PNASA60027-842410.1073/pnas.1210548110]. Here, via investigation of a model swimmer whose flow field closely resembles that of C. reinhardtii, we show that the scattering from a wall can be purely hydrodynamic and that no mechanical or flagellar force is needed for sensing and escaping the boundary.

State diagram of a three-sphere microswimmer in a channel

Geometric confinements are frequently encountered in soft matter systems and in particular significantly alter the dynamics of swimming microorganisms in viscous media. Surface-related effects on the motility of microswimmers can lead to important consequences in a large number of biological systems, such as biofilm formation, bacterial adhesion and microbial activity. On the basis of low-Reynolds-number hydrodynamics, we explore the state diagram of a three-sphere microswimmer under channel confinement in a slit geometry and fully characterize the swimming behavior and trajectories for neutral swimmers, puller- and pusher-type swimmers. While pushers always end up trapped at the channel walls, neutral swimmers and pullers may further perform a gliding motion and maintain a stable navigation along the channel. We find that the resulting dynamical system exhibits a supercritical pitchfork bifurcation in which swimming in the mid-plane becomes unstable beyond a transition channel height while two new stable limit cycles or fixed points that are symmetrically disposed with respect to the channel mid-height emerge. Additionally, we show that an accurate description of the averaged swimming velocity and rotation rate in a channel can be captured analytically using the method of hydrodynamic images, provided that the swimmer size is much smaller than the channel height.

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.

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.