Diffusion, subdiffusion, and localization of active colloids in random post lattices

Behavioral transition of a fish school in a crowded environment

In open water, social fish gather to form schools, in which fish generally align with each other. In this work, we study how this social behavior evolves when perturbed by artificial obstacles. We measure the behavior of a group of zebrafish in the presence of a periodic array of pillars. When the pillar density is low, the fish regroup with a typical interdistance and a well-polarized state with parallel orientations, similarly to their behavior in open-water conditions. Above a critical density of pillars, their social interactions, which are mostly based on vision, are screened and the fish spread randomly through the aquarium, orienting themselves along the free axes of the pillar lattice. The abrupt transition from natural to artificial orientation happens when the pillar interdistance is comparable to the social distance of the fish, i.e., their most probable interdistance. We develop a stochastic model of the relative orientation between fish pairs, taking into account alignment, antialignment, and tumbling, from a distribution biased by the environment. This model provides a good description of the experimental probability distribution of the relative orientation between the fish and captures the behavioral transition. Using the model to fit the experimental data provides qualitative information on the evolution of cognitive parameters, such as the alignment or the tumbling rates, as the pillar density increases. At high pillar density, we find that the artificial environment imposes its geometrical constraints to the fish school, drastically increasing the tumbling rate.

Microbes in porous environments: from active interactions to emergent feedback

Microbes thrive in diverse porous environments-from soil and riverbeds to human lungs and cancer tissues-spanning multiple scales and conditions. Short- to long-term fluctuations in local factors induce spatio-temporal heterogeneities, often leading to physiologically stressful settings. How microbes respond and adapt to such biophysical constraints is an active field of research where considerable insight has been gained over the last decades. With a focus on bacteria, here we review recent advances in self-organization and dispersal in inorganic and organic porous settings, highlighting the role of active interactions and feedback that mediates microbial survival and fitness. We discuss open questions and opportunities for using integrative approaches to advance our understanding of the biophysical strategies which microbes employ at various scales to make porous settings habitable.

Chiral active particles are sensitive reporters to environmental geometry

Chiral active particles (CAPs) are self-propelling particles that break time-reversal symmetry by orbiting or spinning, leading to intriguing behaviors. Here, we examined the dynamics of CAPs moving in 2D lattices of disk obstacles through active Brownian dynamics simulations and granular experiments with grass seeds. We find that the effective diffusivity of the CAPs is sensitive to the structure of the obstacle lattice, a feature absent in achiral active particles. We further studied the transport of CAPs in obstacle arrays under an external field and found a reentrant directional locking effect, which can be used to sort CAPs with different activities. Finally, we demonstrated that parallelogram lattices of obstacles without mirror symmetry can separate clockwise and counter-clockwise CAPs. The mechanisms of the above three novel phenomena are qualitatively explained. As such, our work provides a basis for designing chirality-based tools for single-cell diagnosis and separation, and active particle-based environmental sensors.

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.

Active Brownian particles in random and porous environments

The transport of active particles may occur in complex environments, in which it emerges from the interplay between the mobility of the active components and the quenched disorder of the environment. Here, we explore the structural and dynamical properties of active Brownian particles (ABPs) in random environments composed of fixed obstacles in three dimensions. We consider different arrangements of the obstacles. In particular, we consider two particular situations corresponding to experimentally realizable settings. First, we model pinning particles in (non-overlapping) random positions and, second, in a percolating gel structure and provide an extensive characterization of the structure and dynamics of ABPs in these complex environments. We find that the confinement increases the heterogeneity of the dynamics, with new populations of absorbed and localized particles appearing close to the obstacles. This heterogeneity has a profound impact on the motility induced phase separation exhibited by the particles at high activity, ranging from nucleation and growth in random disorder to a complex phase separation in porous environments.

Boundary design regulates the diffusion of active matter in heterogeneous environments

Physical boundaries play a key role in governing the overall transport properties of nearby self-propelled particles. In this work, we develop dispersion theories and conduct Brownian dynamics simulations to predict the coupling between surface accumulation and effective diffusivity of active particles in boundary-rich media. We focus on three models that are well-understood for passive systems: particle transport in (i) an array of fixed volume-excluding obstacles; (ii) a pore with spatially heterogeneous width; and (iii) a tortuous path with kinks and corners. While the impact of these entropic barriers on passive particle transport is well established, we find that these classical models of porous media flows break down due to the unique interplay between activity and the microstructure of the internal geometry. We study the activity-induced slowdown of effective diffusivity by formulating a Smoluchowski description of long-time self diffusivity which contains contributions from the density and fluctuation fields of the active particles. Particle-based and finite element simulations corroborate this perspective and reveal important nonequilibrium considerations of active transport.

Motor-driven advection competes with crowding to drive spatiotemporally heterogeneous transport in cytoskeleton composites

The cytoskeleton-a composite network of biopolymers, molecular motors, and associated binding proteins-is a paradigmatic example of active matter. Particle transport through the cytoskeleton can range from anomalous and heterogeneous subdiffusion to superdiffusion and advection. Yet, recapitulating and understanding these properties-ubiquitous to the cytoskeleton and other out-of-equilibrium soft matter systems-remains challenging. Here, we combine light sheet microscopy with differential dynamic microscopy and single-particle tracking to elucidate anomalous and advective transport in actomyosin-microtubule composites. We show that particles exhibit multi-mode transport that transitions from pronounced subdiffusion to superdiffusion at tunable crossover timescales. Surprisingly, while higher actomyosin content increases the range of timescales over which transport is superdiffusive, it also markedly increases the degree of subdiffusion at short timescales and generally slows transport. Corresponding displacement distributions display unique combinations of non-Gaussianity, asymmetry, and non-zero modes, indicative of directed advection coupled with caged diffusion and hopping. At larger spatiotemporal scales, particles in active composites exhibit superdiffusive dynamics with scaling exponents that are robust to changing actomyosin fractions, in contrast to normal, yet faster, diffusion in networks without actomyosin. Our specific results shed important new light on the interplay between non-equilibrium processes, crowding and heterogeneity in active cytoskeletal systems. More generally, our approach is broadly applicable to active matter systems to elucidate transport and dynamics across scales.

Chromatographic separation of active polymer-like worm mixtures by contour length and activity

The convective transport rate of polymers through confined geometries depends on their size, allowing for size-based separation of polymer mixtures (chromatography). Here, we investigate whether mixtures of active polymers can be separated in a similar manner based on their activity. We use thin, living Tubifex tubifex worms as a model system for active polymers and study the transport of these worms by an imposed flow through a channel filled with a hexagonal pillar array. The transport rate through the channel depends strongly on the degree of activity, an effect that we assign to the different distribution of conformations sampled by the worms depending on their activity. Our results demonstrate a unique way to sort mixtures of active polymers based on their activity and provide a versatile and convenient experimental system to investigate the hydrodynamics of active polymers.

Repulsive torques alone trigger crystallization of constant speed active particles

We investigate the possibility for self-propelled particles to crystallize without reducing their intrinsic speed. We illuminate how, in the absence of any force, the competition between self-propulsion and repulsive torques determines the macroscopic phases of constant-speed active particles. This minimal model expands upon existing approaches for an improved understanding of crystallization of active matter.

Susceptibility of Orientationally Ordered Active Matter to Chirality Disorder

We investigate the susceptibility of long-range ordered phases of two-dimensional dry aligning active matter to population disorder, taken in the form of a distribution of intrinsic individual chiralities. Using a combination of particle-level models and hydrodynamic theories derived from them, we show that while in finite systems all ordered phases resist a finite amount of such chirality disorder, the homogeneous ones (polar flocks and active nematics) are unstable to any amount of disorder in the infinite-size limit. On the other hand, we find that the inhomogeneous solutions of the coexistence phase (bands) may resist a finite amount of chirality disorder even asymptotically.

Mode-coupling theory for tagged-particle motion of active Brownian particles

We derive a mode-coupling theory (MCT) to describe the dynamics of a tracer particle that is embedded in a dense system of active Brownian particles (ABPs) in two spatial dimensions. The ABP undergo translational and rotational Brownian motion and are equipped with a fixed self-propulsion speed along their orientational vector that describes their active motility. The resulting equations of motion for the tagged-particle density-correlation functions describe the various cases of tracer dynamics close to the glass transition: that of a single active particle in a glass-forming passive host suspensions, that of a passive colloidal particle in a suspension of ABP, and that of active tracers in a bath of active particles. Numerical results are presented for these cases assuming hard-sphere interactions among the particles. The qualitative and quantitative accuracy of the theory is tested against event-driven Brownian dynamics (ED-BD) simulations of active and passive hard disks. Simulation and theory are found in quantitative agreement, provided one adjusts the overall density (as known from the passive description of glassy dynamics), and allows for a rescaling of self-propulsion velocities in the active host system. These adjustments account for the fact that ABP-MCT generally overestimates the tendency for kinetic arrest. We confirm in the simulations a peculiar feature of the transient and stationary dynamical density-correlation functions regarding their lack of symmetry under time reversal, demonstrating the nonequilibrium nature of the system and how it manifests itself in the theory.

Active matter shepherding and clustering in inhomogeneous environments

We consider a mixture of active and passive run-and-tumble disks in an inhomogeneous environment where only half of the sample contains quenched disorder or pinning. The disks are initialized in a fully mixed state of uniform density. We identify several distinct dynamical phases as a function of motor force and pinning density. At high pinning densities and high motor forces, there is a two-step process initiated by a rapid accumulation of both active and passive disks in the pinned region, which produces a large density gradient in the system. This is followed by a slower species phase separation process where the inactive disks are shepherded by the active disks into the pin-free region, forming a nonclustered fluid and producing a more uniform density with species phase separation. For higher pinning densities and low motor forces, the dynamics becomes very slow and the system maintains a strong density gradient. For weaker pinning and large motor forces, a floating clustered state appears, and the time-averaged density of the system is uniform. We illustrate the appearance of these phases in a dynamic phase diagram.

Interactions of E. coli with cylindrical micro-pillars of different geometric modifications

Understanding the behavior of bacteria at the proximity of different surfaces is of great importance and interest. Despite recent exciting progress in geometric control of bacterial behavior around surfaces, a detailed comparison on the interaction of bacteria with cylindrical surfaces of different geometric modifications is still missing. Here, we investigated how bacteria interacted with cylindrical micro-pillars and modified cylindrical micro-pillars with sprocket, gear, and flower-like wall surface features. Using phase-contrast microscopy, we examined the motion of bacteria around the micro-pillars, and observed different responses of bacteria to each geometric modification. In addition, we extracted the trajectories of the bacteria and characterized several parameters (instantaneous velocity v, change of direction δ, approaching angle ϕ) to quantitatively compare the effects of the geometric modifications on the micro-pillars. We found that sharp spikes showed the largest effect, compared to smooth surface, convex and concave ripples. Lastly, we carried out numerical simulations, which explained the experimental observations and showed that the observed effects were due to the geometric modifications.

High-dimensional percolation criticality and hints of mean-field-like caging of the random Lorentz gas

The random Lorentz gas (RLG) is a minimal model for transport in disordered media. Despite the broad relevance of the model, theoretical grasp over its properties remains weak. For instance, the scaling with dimension d of its localization transition at the void percolation threshold is not well controlled analytically nor computationally. A recent study [Biroli et al., Phys. Rev. E 103, L030104 (2021)2470-004510.1103/PhysRevE.103.L030104] of the caging behavior of the RLG motivated by the mean-field theory of glasses has uncovered physical inconsistencies in that scaling that heighten the need for guidance. Here we first extend analytical expectations for asymptotic high-d bounds on the void percolation threshold and then computationally evaluate both the threshold and its criticality in various d. In high-d systems, we observe that the standard percolation physics is complemented by a dynamical slowdown of the tracer dynamics reminiscent of mean-field caging. A simple modification of the RLG is found to bring the interplay between percolation and mean-field-like caging down to d=3.

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

We examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.

Directional locking and the influence of obstacle density on skyrmion dynamics in triangular and honeycomb arrays

We numerically examine the dynamics of a single skyrmion driven over triangular and honeycomb obstacle arrays at zero temperature. The skyrmion Hall angleθ_sk, defined as the angle between the applied external drive and the direction of the skyrmion motion, increases in quantized steps or continuously as a function of the applied drive. For the obstacle arrays studied in this work, the skyrmion exhibits two main directional locking angles ofθ_sk= -30° and -60°. We show that these directions are privileged due to the obstacle landscape symmetry, and coincide with channels along which the skyrmion may move with few or no obstacle collisions. Here we investigate how changes in the obstacle density can modify the skyrmion Hall angles and cause some dynamic phases to appear or grow while other phases vanish. This interesting behavior can be used to guide skyrmions along designated trajectories via regions with different obstacle densities. For fixed obstacle densities, we investigate the evolution of the lockedθ_sk= -30° and -60° phases as a function of the Magnus force, and discuss possibilities for switching between these phases using topological selection.

Breakdown of Ergodicity and Self-Averaging in Polar Flocks with Quenched Disorder

We show that spatial quenched disorder affects polar active matter in ways more complex and far reaching than heretofore believed. Using simulations of the 2D Vicsek model subjected to random couplings or a disordered scattering field, we find in particular that ergodicity is lost in the ordered phase, the nature of which we show to depend qualitatively on the type of quenched disorder: for random couplings, it remains long-range ordered, but qualitatively different from the pure (disorderless) case. For random scatterers, polar order varies with system size but we find strong non-self-averaging, with sample-to-sample fluctuations dominating asymptotically, which prevents us from elucidating the asymptotic status of order.

Time-dependent inertia of self-propelled particles: The Langevin rocket

Many self-propelled objects are large enough to exhibit inertial effects but still suffer from environmental fluctuations. The corresponding basic equations of motion are governed by active Langevin dynamics, which involve inertia, friction, and stochastic noise for both the translational and orientational degrees of freedom coupled via the self-propulsion along the particle orientation. In this paper, we generalize the active Langevin model to time-dependent parameters and explicitly discuss the effect of time-dependent inertia for achiral and chiral particles. Realizations of this situation are manifold, ranging from minirockets (which are self-propelled by burning their own mass), to dust particles in plasma (which lose mass by evaporating material), to walkers with expiring activity. Here we present analytical solutions for several dynamical correlation functions, such as mean-square displacement and orientational and velocity autocorrelation functions. If the parameters exhibit a slow power law in time, we obtain anomalous superdiffusion with a nontrivial dynamical exponent. Finally, we constitute the "Langevin rocket" model by including orientational fluctuations in the traditional Tsiolkovsky rocket equation. We calculate the mean reach of the Langevin rocket and discuss different mass ejection strategies to maximize it. Our results can be tested in experiments on macroscopic robotic or living particles or in self-propelled mesoscopic objects moving in media of low viscosity, such as complex plasma.

Emergence of dynamic vortex glasses in disordered polar active fluids

In equilibrium, disorder conspires with topological defects to redefine the ordered states of matter in systems as diverse as crystals, superconductors, and liquid crystals. Far from equilibrium, however, the consequences of quenched disorder on active condensed matter remain virtually uncharted. Here, we reveal a state of strongly disordered active matter with no counterparts in equilibrium: a dynamical vortex glass. Combining microfluidic experiments and theory, we show how colloidal flocks collectively cruise through disordered environments without relaxing the topological singularities of their flows. The resulting state is highly dynamical but the flow patterns, shaped by a finite density of frozen vortices, are stationary and exponentially degenerated. Quenched isotropic disorder acts as a random gauge field turning active liquids into dynamical vortex glasses. We argue that this robust mechanism should shape the collective dynamics of a broad class of disordered active matter, from synthetic active nematics to collections of living cells exploring heterogeneous media.

Active matter commensuration and frustration effects on periodic substrates

We show that self-driven particles coupled to a periodic obstacle array exhibit active matter commensuration effects that are absent in the Brownian limit. As the obstacle size is varied for sufficiently large activity, a series of commensuration effects appear in which the motility induced phase separation produces commensurate crystalline states, while for other obstacle sizes we find frustrated or amorphous states. The commensuration effects are associated with peaks in the amount of sixfold ordering and the maximum cluster size. When a drift force is added to the system, the mobility contains peaks and dips similar to those found in transport studies for commensuration effects in superconducting vortices and colloidal particles.

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.

Active Brownian and inertial particles in disordered environments: Short-time expansion of the mean-square displacement

We consider an active Brownian particle moving in a disordered two-dimensional energy or motility landscape. The averaged mean-square displacement (MSD) of the particle is calculated analytically within a systematic short-time expansion. As a result, for overdamped particles, both an external random force field and disorder in the self-propulsion speed induce ballistic behavior adding to the ballistic regime of an active particle with sharp self-propulsion speed. Spatial correlations in the force and motility landscape contribute only to the cubic and higher-order powers in time for the MSD. Finally, for inertial particles two superballistic regimes are found where the scaling exponent of the MSD with time is α=3 and α=4. We confirm our theoretical predictions by computer simulations. Moreover, they are verifiable in experiments on self-propelled colloids in random environments.

Directional locking effects for active matter particles coupled to a periodic substrate

Directional locking occurs when a particle moving over a periodic substrate becomes constrained to travel along certain substrate symmetry directions. Such locking effects arise for colloids and superconducting vortices moving over ordered substrates when the direction of the external drive is varied. Here we study the directional locking of run-and-tumble active matter particles interacting with a periodic array of obstacles. In the absence of an external biasing force, we find that the active particle motion locks to various symmetry directions of the substrate when the run time between tumbles is large. The number of possible locking directions depends on the array density and on the relative sizes of the particles and the obstacles. For a square array of large obstacles, the active particle only locks to the x, y, and 45^{∘} directions, while for smaller obstacles, the number of locking angles increases. Each locking angle satisfies θ=arctan(p/q), where p and q are integers, and the angle of motion can be measured using the ratio of the velocities or the velocity distributions in the x and y directions. When a biasing driving force is applied, the directional locking behavior is affected by the ratio of the self-propulsion force to the biasing force. For large biasing, the behavior resembles that found for directional locking in passive systems. For large obstacles under biased driving, a trapping behavior occurs that is nonmonotonic as a function of increasing run length or increasing self-propulsion force, and the trapping diminishes when the run length is sufficiently large.

Crowding-Enhanced Diffusion: An Exact Theory for Highly Entangled Self-Propelled Stiff Filaments

We study a strongly interacting crowded system of self-propelled stiff filaments by event-driven Brownian dynamics simulations and an analytical theory to elucidate the intricate interplay of crowding and self-propulsion. We find a remarkable increase of the effective diffusivity upon increasing the filament number density by more than one order of magnitude. This counterintuitive "crowded is faster" behavior can be rationalized by extending the concept of a confining tube pioneered by Doi and Edwards for highly entangled, crowded, passive to active systems. We predict a scaling theory for the effective diffusivity as a function of the Péclet number and the filament number density. Subsequently, we show that an exact expression derived for a single self-propelled filament with motility parameters as input can predict the nontrivial spatiotemporal dynamics over the entire range of length and timescales. In particular, our theory captures short-time diffusion, directed swimming motion at intermediate times, and the transition to complete orientational relaxation at long times.

Confined diffusion in a random Lorentz gas environment

We study the diffusive behavior of biased Brownian particles in a two dimensional confined geometry filled with the freezing obstacles. The transport properties of these particles are investigated for various values of the obstacle density η and the scaling parameter f, which is the ratio of work done to the particles to available thermal energy. We show that, when the thermal fluctuations dominate over the external force, i.e., small f regime, particles get trapped in the given environment when the system percolates at the critical obstacle density η_{c}≈1.2. However, as f increases, we observe that particle trapping occurs prior to η_{c}. In particular, we find a relation between η and f which provides an estimate of the minimum η up to a critical scaling parameter f_{c} beyond which the Fick-Jacobs description is invalid. Prominent transport features like nonmonotonic behavior of the nonlinear mobility, anomalous diffusion, and greatly enhanced effective diffusion coefficient are explained for various strengths of f and η. Also, it is interesting to observe that particles exhibit different kinds of diffusive behaviors, i.e., subdiffusion, normal diffusion, and superdiffusion. These findings, which are genuine to the confined and random Lorentz gas environment, can be useful to understand the transport of small particles or molecules in systems such as molecular sieves and porous media, which have a complex heterogeneous environment of the freezing obstacles.

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".

Overdamped Brownian dynamics in piecewise-defined energy landscapes

We study the overdamped Brownian dynamics of particles moving in piecewise-defined potential energy landscapes U(x), where the height Q of each section is obtained from the exponential distribution p(Q)=aβexp(-aβQ), where β is the reciprocal thermal energy, and a>0. The averaged effective diffusion coefficient 〈D_{eff}〉 is introduced to characterize the diffusive motion: 〈x^{2}〉=2〈D_{eff}〉t. A general expression for 〈D_{eff}〉 in terms of U(x) and p(Q) is derived and then applied to three types of energy landscape: flat sections, smooth maxima, and sharp maxima. All three cases display a transition between subdiffusive and diffusive behavior at a=1, and a reduction to free diffusion as a→∞. The behavior of 〈D_{eff}〉 around the transition is investigated and found to depend heavily upon the shape of the maxima: Energy landscapes made up of flat sections or smooth maxima display power-law behavior, while for landscapes with sharp maxima, strongly divergent behavior is observed. Two aspects of the subdiffusive regime are studied: the growth of the mean squared displacement with time and the distribution of mean first-passage times. For the former, agreement between Brownian dynamics simulations and a coarse-grained equivalent was observed, but the results deviated from the random barrier model's predictions. The discrepancy could be a finite-time effect. For the latter, agreement between the characteristic exponent calculated numerically and that predicted by the random barrier model is observed in the large-amplitude limit.

A practical guide to active colloids: choosing synthetic model systems for soft matter physics research

Synthetic active colloids that harvest energy stored in the environment and swim autonomously are a popular model system for active matter. This emerging field of research sits at the intersection of materials chemistry, soft matter physics, and engineering, and thus cross-talk among researchers from different backgrounds becomes critical yet difficult. To facilitate this interdisciplinary communication, and to help soft matter physicists with choosing the best model system for their research, we here present a tutorial review article that describes, in appropriate detail, six experimental systems of active colloids commonly found in the physics literature. For each type, we introduce their background, material synthesis and operating mechanisms and notable studies from the soft matter community, and comment on their respective advantages and limitations. In addition, the main features of each type of active colloid are summarized into two useful tables. As materials chemists and engineers, we intend for this article to serve as a practical guide, so those who are not familiar with the experimental aspects of active colloids can make more informed decisions and maximize their creativity.

Numerical simulations of self-diffusiophoretic colloids at fluid interfaces

The dynamics of active colloids is very sensitive to the presence of boundaries and interfaces which therefore can be used to control their motion. Here we analyze the dynamics of active colloids adsorbed at a fluid-fluid interface. By using a mesoscopic numerical approach which relies on an approximated numerical solution of the Navier-Stokes equation, we show that when adsorbed at a fluid interface, an active colloid experiences a net torque even in the absence of a viscosity contrast between the two adjacent fluids. In particular, we study the dependence of this torque on the contact angle of the colloid with the fluid-fluid interface and on its surface properties. We rationalize our results via an approximate approach which accounts for the appearance of a local friction coefficient. By providing insight into the dynamics of active colloids adsorbed at fluid interfaces, our results are relevant for two-dimensional self assembly and emulsion stabilization by means of active colloids.

Topotaxis of active Brownian particles

Recent experimental studies have demonstrated that cellular motion can be directed by topographical gradients, such as those resulting from spatial variations in the features of a micropatterned substrate. This phenomenon, known as topotaxis, has been extensively studied for topographical gradients at the subcellular scale, but can also occur in the presence of a spatially varying density of cell-sized features. Such a large-scale topotaxis has recently been observed in highly motile cells that persistently crawl within an array of obstacles with smoothly varying lattice spacing. We introduce a toy model of large-scale topotaxis, based on active Brownian particles. Using numerical simulations and analytical arguments, we demonstrate that topographical gradients introduce a spatial modulation of the particles' persistence, leading to directed motion toward regions of higher persistence. Our results demonstrate that persistent motion alone is sufficient to drive large-scale topotaxis and could serve as a starting point for more detailed studies on self-propelled particles and cells.

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.

Tuning the Random Walk of Active Colloids: From Individual Run-and-Tumble to Dynamic Clustering

Active particles such as swimming bacteria or self-propelled colloids spontaneously self-organize into large-scale dynamic structures. The emergence of these collective states from the motility pattern of the individual particles, typically a random walk, is yet to be probed in a well-defined synthetic system. Here, we report the experimental realization of tunable colloidal motion that reproduces run-and-tumble and Lévy trajectories. We utilize the Quincke effect to achieve controlled sequences of repeated particle runs and random reorientations. We find that a population of these random walkers exhibit behaviors reminiscent of bacterial suspensions such as dynamic clusters and mesoscale turbulentlike flows.

Enhanced propagation of motile bacteria on surfaces due to forward scattering

How motile bacteria move near a surface is a problem of fundamental biophysical interest and is key to the emergence of several phenomena of biological, ecological and medical relevance, including biofilm formation. Solid boundaries can strongly influence a cell's propulsion mechanism, thus leading many flagellated bacteria to describe long circular trajectories stably entrapped by the surface. Experimental studies on near-surface bacterial motility have, however, neglected the fact that real environments have typical microstructures varying on the scale of the cells' motion. Here, we show that micro-obstacles influence the propagation of peritrichously flagellated bacteria on a flat surface in a non-monotonic way. Instead of hindering it, an optimal, relatively low obstacle density can significantly enhance cells' propagation on surfaces due to individual forward-scattering events. This finding provides insight on the emerging dynamics of chiral active matter in complex environments and inspires possible routes to control microbial ecology in natural habitats.

Generalized Ornstein-Uhlenbeck model for active motion

We investigate a one-dimensional model of active motion, which takes into account the effects of persistent self-propulsion through a memory function in a dissipative-like term of the generalized Langevin equation for particle swimming velocity. The proposed model is a generalization of the active Ornstein-Uhlenbeck model introduced by G. Szamel [Phys. Rev. E 90, 012111 (2014)10.1103/PhysRevE.90.012111]. We focus on two different kinds of memory which arise in many natural systems: an exponential decay and a power law, supplemented with additive colored noise. We provide analytical expressions for the velocity autocorrelation function and the mean-squared displacement, which are in excellent agreement with numerical simulations. For both models, damped oscillatory solutions emerge due to the competition between the memory of the system and the persistence of velocity fluctuations. In particular, for a power-law model with fractional Brownian noise, we show that long-time active subdiffusion occurs with increasing long-term memory.

Quincke rotor dynamics in confinement: rolling and hovering

The Quincke effect is an electrohydrodynamic instability which gives rise to a torque on a dielectric particle in a uniform DC electric field. Previous studies reported that a sphere initially resting on the electrode rolls with steady velocity. We experimentally find that in strong fields the rolling becomes unsteady, with time-periodic velocity. Furthermore, we find another regime, where the rotating sphere levitates in the space between the electrodes. Our experimental results show that the onset of Quincke rotation strongly depends on particle confinement and the threshold for rolling is higher compared to rotation in the hovering state.

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.

Ideal circle microswimmers in crowded media

Microswimmers are exposed in nature to crowded environments and their transport properties depend in a subtle way on the interaction with obstacles. Here, we investigate a model for a single ideal circle swimmer exploring a two-dimensional disordered array of impenetrable obstacles. The microswimmer moves on circular orbits in the freely accessible space and follows the surface of an obstacle for a certain time upon collision. Depending on the obstacle density and the radius of the circular orbits, the microswimmer displays either long-range transport or is localized in a finite region. We show that there are transitions from two localized states to a diffusive state each driven by an underlying static percolation transition. We determine the non-equilibrium state diagram and calculate the mean-square displacements and diffusivities by computer simulations. Close to the transition lines transport becomes subdiffusive which is rationalized as a dynamic critical phenomenon.

Anomalous mobility of a driven active particle in a steady laminar flow

We study, via extensive numerical simulations, the force-velocity curve of an active particle advected by a steady laminar flow, in the nonlinear response regime. Our model for an active particle relies on a colored noise term that mimics its persistent motion over a time scale [Formula: see text]. We find that the active particle dynamics shows non-trivial effects, such as negative differential and absolute mobility (NDM and ANM, respectively). We explore the space of the model parameters and compare the observed behaviors with those obtained for a passive particle ([Formula: see text]) advected by the same laminar flow. Our results show that the phenomena of NDM and ANM are quite robust with respect to the details of the considered noise: in particular for finite [Formula: see text] a more complex force-velocity relation can be observed.

Negative differential mobility in interacting particle systems

Driven particles in the presence of crowded environment, obstacles, or kinetic constraints often exhibit negative differential mobility (NDM) due to their decreased dynamical activity. Based on the empirical studies of conserved lattice gas model, two species exclusion model and other interacting particle systems we propose a new mechanism for complex many-particle systems where slowing down of certain non-driven degrees of freedom by the external field can give rise to NDM. To prove that the slowing down of the non-driven degrees is indeed the underlying cause, we consider several driven diffusive systems including two species exclusion models, misanthrope process, and show from the exact steady state results that NDM indeed appears when some non-driven modes are slowed down deliberately. For clarity, we also provide a simple pedagogical example of two interacting random walkers on a ring which conforms to the proposed scenario.

Clogging and depinning of ballistic active matter systems in disordered media

We numerically examine ballistic active disks driven through a random obstacle array. Formation of a pinned or clogged state occurs at much lower obstacle densities for the active disks than for passive disks. As a function of obstacle density, we identify several distinct phases including a depinned fluctuating cluster state, a pinned single-cluster or jammed state, a pinned multicluster state, a pinned gel state, and a pinned disordered state. At lower active disk densities, a drifting uniform liquid forms in the absence of obstacles, but when even a small number of obstacles are introduced, the disks organize into a pinned phase-separated cluster state in which clusters nucleate around the obstacles, similar to a wetting phenomenon. We examine how the depinning threshold changes as a function of disk or obstacle density and find a crossover from a collectively pinned cluster state to a disordered plastic depinning transition as a function of increasing obstacle density. We compare this to the behavior of nonballistic active particles and show that as we vary the activity from completely passive to completely ballistic, a clogged phase-separated state appears in both the active and passive limits, while for intermediate activity, a readily flowing liquid state appears and there is an optimal activity level that maximizes the flux through the sample.