Dynamic phases of active matter systems with quenched disorder

Phase separation, edge currents, and Hall effect for active matter with Magnus dynamics

We examine run-and-tumble disks in two-dimensional systems where the particles also have a Magnus component to their dynamics. For increased activity, we find that the system forms a motility-induced phase-separated (MIPS) state with chiral edge flow around the clusters, where the direction of the current is correlated with the sign of the Magnus term. The stability of the MIPS state is non-monotonic as a function of increasing Magnus term amplitude, with the MIPS region first extending down to lower activities followed by a break up of MIPS at large Magnus amplitudes into a gel-like state. We examine the dynamics in the presence of quenched disorder and a uniform drive and find that the bulk flow exhibits a drive-dependent Hall angle. This is a result of the side jump effect produced by scattering from the pinning sites and is similar to the behavior found for skyrmions in chiral magnets with quenched disorder.

Solvable model of driven matter with pinning

We present a simple model of driven matter in a 1D medium with pinning impurities, applicable to magnetic domains walls, confined colloids, and other systems. We find rich dynamics, including hysteresis, reentrance, quasiperiodicity, and two distinct routes to chaos. In contrast to other minimal models of driven matter, the model is solvable: we derive the full phase diagram for small N, and for large N, we derive expressions for order parameters and several bifurcation curves. The model is also realistic. Its collective states match those seen in the experiments of magnetic domain walls.

Reversible, irreversible, and mixed regimes for periodically driven disks in random obstacle arrays

We examine an assembly of repulsive disks interacting with a random obstacle array under a periodic drive and find a transition from reversible to irreversible dynamics as a function of drive amplitude or disk density. At low densities and drives, the system rapidly forms a reversible state where the disks return to their exact positions at the end of each cycle. In contrast, at high amplitudes or high densities, the system enters an irreversible state where the disks exhibit normal diffusion. Between these two regimes, there can be an intermediate irreversible state where most of the system is reversible, but localized irreversible regions are present that are prevented from spreading through the system due to a screening effect from the obstacles. We also find states that we term "combinatorial reversible states" in which the disks return to their original positions after multiple driving cycles. In these states, individual disks exchange positions but form the same configurations during the subcycles of the larger reversible cycle.

Dynamic phases and combing effects for elongated particles moving over quenched disorder

We consider a two-dimensional system of elongated particles driven over a landscape containing randomly placed pinning sites. For varied pinning site density, external drive magnitude, and particle elongation, we find a wide variety of dynamic phases, including random structures, stripe or combed phases with nematic order, and clogged states. The different regimes can be identified by examining nematic ordering, cluster size, number of pinned particles, and transverse diffusion. In some regimes we find that the pinning can enhance the particle alignment, producing a nonmonotonic signature in the nematic ordering with a maximum at a particular combination of pinning density and drive. The optimal nematic occurs when a sufficient number of particles can be pinned, generating a local shear and leading to what we call a combing effect. At high drives, the combing effect is reduced when the number of pinned particles decreases. For stronger pinning, the particles form a heterogeneous clustered or clogged state that depins into a fluctuating state with high diffusion.

Transition from susceptible-infected to susceptible-infected-recovered dynamics in a susceptible-cleric-zombie-recovered active matter model

The susceptible-infected (SI) and susceptible-infected-recovered (SIR) models provide two distinct representations of epidemic evolution, distinguished by whether or not the number of susceptibles always drops to zero at long times. Here we introduce a new active matter epidemic model, the "susceptible-cleric-zombie-recovered" (SCZR) model, in which spontaneous recovery is absent but zombies can recover with probability γ via interaction with a cleric. Upon colliding with a zombie, both susceptibles and clerics enter the zombie state with probability β and α, respectively. By changing the initial fraction of clerics or their healing ability rate γ, we can tune the SCZR model between SI dynamics, in which no susceptibles or clerics remain at long times, and SIR dynamics, in which a finite number of clerics and susceptibles survive at long times. The model is relevant to certain real world diseases such as HIV where spontaneous recovery is impossible but where medical interventions by a limited number of caregivers can reduce or eliminate the spread of infection.

Pattern formation and flocking for particles near the jamming transition on resource gradient substrates

We numerically examine a bidisperse system of active and passive particles coupled to a resource substrate. The active particles deplete the resource at a fixed rate and move toward regions with higher resources, while all of the particles interact sterically with each other. We show that at high densities, this system exhibits a rich variety of pattern-forming phases along with directed motion or flocking as a function of the relative rates of resource absorption and consumption as well as the active to passive particle ratio. These include partial phase separation into rivers of active particles flowing through passive clusters, strongly phase separated states where the active particles induce crystallization of the passive particles, mixed jammed states, and fluctuating mixed fluid phases. For higher resource recovery rates, we demonstrate that the active particles can undergo motility-induced phase separation, while at high densities, there can be a coherent flock containing only active particles or a solid mixture of active and passive particles. The directed flocking motion typically shows a transient in which the flow switches among different directions before settling into one direction, and there is a critical density below which flocking does not occur. We map out the different phases as function of system density, resource absorption and recovery rates, and the ratio of active to passive particles.

Morphologies and dynamics of free surfaces of crystals composed of active particles

Active matter exhibits various collective motions and nonequilibrium phases, such as crystals; however, their surface properties have been poorly explored. Here, we use Brownian dynamics simulations to investigate the surface morphology and dynamics of two-dimensional active crystals during and after growth. For crystal growth on a substrate, the position and roughness of the crystal surface reach steady states at different times. In the steady state, the surface exhibits superdiffusive behaviour at the short time, and the roughness is insensitive to the roughening process and particle activity. We observe two-stage and three-stage surface roughening at different Péclet numbers. The result of dynamic scaling analysis shows that the surface is similar to anomalous roughening, which is distinct from the normal roughening typically found in conventional passive systems. Capillary wave theory for a thermal equilibrium system can describe the active surface fluctuations only in the long-wavelength regime, indicating that active particles mainly drive the surface out of equilibrium locally. These similarities and differences between the active and passive crystal surfaces are essential for understanding active crystals and interfaces.

Using active matter to introduce spatial heterogeneity to the susceptible infected recovered model of epidemic spreading

The widely used susceptible-infected-recovered (S-I-R) epidemic model assumes a uniform, well-mixed population, and incorporation of spatial heterogeneities remains a major challenge. Understanding failures of the mixing assumption is important for designing effective disease mitigation approaches. We combine a run-and-tumble self-propelled active matter system with an S-I-R model to capture the effects of spatial disorder. Working in the motility-induced phase separation regime both with and without quenched disorder, we find two epidemic regimes. For low transmissibility, quenched disorder lowers the frequency of epidemics and increases their average duration. For high transmissibility, the epidemic spreads as a front and the epidemic curves are less sensitive to quenched disorder; however, within this regime it is possible for quenched disorder to enhance the contagion by creating regions of higher particle densities. We discuss how this system could be realized using artificial swimmers with mobile optical traps operated on a feedback loop.

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

Dense polar active fluids in a disordered environment

We examine the influence of quenched disorder on the flocking transition of dense polar active matter. We consider incompressible systems of active particles with aligning interactions under the effect of either quenched random forces or random dilution. The system displays a continuous disorder-order (flocking) transition, and the associated scaling behavior is described by a new universality class which is controlled by a quenched Navier-Stokes fixed point. We determine the critical exponents through a perturbative renormalization group analysis. We show that the two forms of quenched disorder, random force and random mass (dilution), belong to the same universality class, in contrast with the situation at equilibrium.

Chemotactic smoothing of collective migration

Collective migration-the directed, coordinated motion of many self-propelled agents-is a fascinating emergent behavior exhibited by active matter with functional implications for biological systems. However, how migration can persist when a population is confronted with perturbations is poorly understood. Here, we address this gap in knowledge through studies of bacteria that migrate via directed motion, or chemotaxis, in response to a self-generated nutrient gradient. We find that bacterial populations autonomously smooth out large-scale perturbations in their overall morphology, enabling the cells to continue to migrate together. This smoothing process arises from spatial variations in the ability of cells to sense and respond to the local nutrient gradient-revealing a population-scale consequence of the manner in which individual cells transduce external signals. Altogether, our work provides insights to predict, and potentially control, the collective migration and morphology of cellular populations and diverse other forms of active matter.

The dynamics of chemically propelled dimer motor on a pinning substrate

The dynamics of self-propelled micro-motors, in a thin fluid film containing an attractive substrate, is investigated by means of a particle-based simulation. A chemically powered sphere dimer, consisting of a...

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.

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.

Polar swimmers induce several phases in active nematics

Swimming bacteria in passive nematics in the form of lyotropic liquid crystals are defined as a new class of active matter known as living liquid crystals in recent studies. It has also been shown that liquid crystal solutions are promising candidates for trapping and detecting bacteria. We ask the question, can a similar class of matter be designed for background nematics which are also active? Hence, we developed a minimal model for the mixture of polar particles in active nematics. It is found that the active nematics in such a mixture are highly sensitive to the presence of polar particles and show the formation of large scale higher order structures for a relatively low polar particle density. Upon increasing the density of polar particles, different phases of active nematics are found and it is observed that the system shows two phase transitions. The first phase transition is a first order transition from quasi-long-ranged ordered active nematics to disordered active nematics with larger scale structures. On further increasing density of polar particles, the system transitions to a third phase, where polar particles form large, mutually aligned clusters. These clusters sweep the whole system and enforce local order in the nematics. The current study can be helpful for detecting the presence of very low densities of polar swimmers in active nematics and can be used to design and control different structures in active nematics.

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.

Chemotactic migration of bacteria in porous media

Chemotactic migration of bacteria-their ability to direct multicellular motion along chemical gradients-is central to processes in agriculture, the environment, and medicine. However, current understanding of migration is based on studies performed in bulk liquid, despite the fact that many bacteria inhabit tight porous media such as soils, sediments, and biological gels. Here, we directly visualize the chemotactic migration of Escherichia coli populations in well-defined 3D porous media in the absence of any other imposed external forcing (e.g., flow). We find that pore-scale confinement is a strong regulator of migration. Strikingly, cells use a different primary mechanism to direct their motion in confinement than in bulk liquid. Furthermore, confinement markedly alters the dynamics and morphology of the migrating population-features that can be described by a continuum model, but only when standard motility parameters are substantially altered from their bulk liquid values to reflect the influence of pore-scale confinement. Our work thus provides a framework to predict and control the migration of bacteria, and active matter in general, in complex environments.

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.

Active Brownian particles moving through disordered landscapes

Disordered media are ubiquitous in systems where self-propelled particles are present, ranging from biological settings to synthetic systems, like in active microfluidic devices. Here we investigate the behavior of active Brownian particles that have an internal energy depot and move through a landscape with a quenched frictional disorder. We consider the cases of very fast internal relaxation processes and the limit of strong disorder. Analytical calculations of the mean-square displacement in the fast-relaxation approximation is shown to agree well with numerically integrated energy depot dynamics and predict normal dispersion for a bounded drag coefficient and anomalous dispersion for power-law dependence of the drag on spatial coordinates. Furthermore, we show that in the strongly disordered limit the self-propulsion speed can, for practical purposes, be considered a fluctuating quantity. Distributions of self-propulsion speeds are investigated numerically for different parameter choices.

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.

Flocking with a q-fold discrete symmetry: Band-to-lane transition in the active Potts model

We study the q-state active Potts model (APM) on a two-dimensional lattice in which self-propelled particles have q internal states corresponding to the q directions of motion. A local alignment rule inspired by the ferromagnetic q-state Potts model and self-propulsion via biased diffusion according to the internal particle states elicits collective motion at high densities and low noise. We formulate a coarse-grained hydrodynamic theory with which we compute the phase diagrams of the APM for q=4 and q=6 and analyze the flocking dynamics in the coexistence region, where the high-density (polar liquid) phase forms a fluctuating stripe of coherently moving particles on the background of the low-density (gas) phase. A reorientation transition of the phase-separated profiles from transversal band motion to longitudinal lane formation is found, which is absent in the Vicsek model and the active Ising model. The origin of this reorientation transition is revealed by a stability analysis: for large velocities the transverse diffusivity approaches zero and stabilizes lanes. Computer simulations corroborate the analytical predictions of the flocking and reorientation transitions and validate the phase diagrams of the APM.

Speed inhomogeneity accelerates information transfer in polar flock

A collection of self-propelled particles (SPPs) shows coherent motion and exhibits a true long-range-ordered state in two dimensions. Various studies show that the presence of spatial inhomogeneities can destroy the usual long-range ordering in the system. However, the effects of inhomogeneity due to the intrinsic properties of the particles are barely addressed. In this paper we consider a collection of polar SPPs moving at inhomogeneous speed (IS) on a two-dimensional substrate, which can arise due to varying physical strengths of the individual particles. To our surprise, the IS not only preserves the usual long-range ordering present in homogeneous speed models but also induces faster ordering in the system. Furthermore, the response of the flock to an external perturbation is also faster, compared to the Vicsek-like model systems, due to the frequent update of neighbors of each SPP in the presence of the IS. Therefore, our study shows that an IS can promote information transfer in a moving flock.

Motile dissenters disrupt the flocking of active granular matter

We report flocking in the dry active granular matter of millimeter-sized two-step-tapered rods without an intervening medium. The system undergoes the flocking phase transition at a threshold area fraction of ∼0.12 having high orientational correlations between the particles. However, the one-step-tapered rods do not flock and are used as the motile dissenters in the flock-forming granular matter. At the critical fraction of dissenters of ∼0.3, the flocking order of the system gets completely destroyed. The variance of the system's order parameter shows a maximum near the dissenter fraction f∼0.05, suggesting a finite-size crossover between the ordered and disordered phases. Our experiments bring out the disruption of the cooperative behavior in heterogeneous active systems with possible implications in real-life examples.

Detecting depinning and nonequilibrium transitions with unsupervised machine learning

Using numerical simulations of a model disk system, we demonstrate that a machine learning generated order-parameter-like measure can detect depinning transitions and different dynamic flow phases in systems driven far from equilibrium. We specifically consider monodisperse passive disks with short range interactions undergoing a depinning phase transition when driven over quenched disorder. The machine learning derived order-parameter-like measure identifies the depinning transition as well as different dynamical regimes, such as the transition from a flowing liquid to a phase separated liquid-solid state that is not readily distinguished with traditional measures such as velocity-force curves or Voronoi tessellation. The order-parameter-like measure also shows markedly distinct behavior in the limit of high density where jamming effects occur. Our results should be general to the broad class of particle-based systems that exhibit depinning transitions and nonequilibrium phase transitions.

Enhanced dynamics of active Brownian particles in periodic obstacle arrays and corrugated channels

We study the motion of an active Brownian particle (ABP) using the overdamped Langevin dynamics on a two-dimensional substrate with periodic array of obstacles and in a quasi-one-dimensional corrugated channel comprised of periodically arrayed obstacles. The periodic arrangement of the obstacles enhances the persistent motion of the ABP in comparison to its motion in the free space. Persistent motion increases with the activity of the ABP. We note that the periodic arrangement induces directionality in ABP motion at late time, and it increases with the size of the obstacles. We also note that the ABP exhibits a super-diffusive dynamics in the corrugated channel. The transport property is independent of the shape of the channel; rather it depends on the packing fraction of the obstacles in the system. However, the ABP shows the usual diffusive dynamics in the quasi-one-dimensional channel with flat boundary.

Flow and clogging of particles in shaking random obstacles

Transport of three types of particles (passive particles, active particles, and polar particles) is investigated in a random obstacle array in the presence of a dc drift force. The obstacles are static or synchronously shake along the given direction. When the obstacles are static, the average velocity is a peaked function of the dc drift force (negative differential mobility) for low particle density, while the average velocity monotonically increases with the dc drift force (positive differential mobility) for high particle density. Under the same conditions, passive particles are most likely to pass through the obstacles, while polar particles are easily trapped by the obstacles. The polar alignment can strongly reduce the overall mobility of particles. When the obstacles shake along the given direction, the optimal shaking frequency or amplitude can maximize the average velocity. It is more effective to reduce clogging for the transverse shaking than that for the longitudinal shaking.

Dynamic phases, stratification, laning, and pattern formation for driven bidisperse disk systems in the presence of quenched disorder

Using numerical simulations, we examine the dynamics of driven two-dimensional bidisperse disks flowing over quenched disorder. The system exhibits a series of distinct dynamical phases as a function of applied driving force and packing fraction, including a phase-separated state as well as a smectic state with liquid-like or polycrystalline features. At low driving forces, we find a clogged phase with an isotropic density distribution, while at intermediate driving forces the disks separate into bands of high and low densities with either liquid-like or polycrystalline structure in the high-density bands. In addition to the density phase separation, we find that in some cases there is a fractionation of the disk species, particularly when the disk size ratio is large. The species phase-separated regimes form a variety of patterns such as large disks separated by chains of smaller disks. Our results show that the formation of laning states can be enhanced by tuning the ratio of disk radius of the two species, due to the clumping of small disks in the interstitial regions between the large disks. This system could be experimentally realized using sterically interacting colloidal particles suspended in a viscous fluid driven over random pinning arrays or granular matter suspended in fluid moving over a random landscape.

Controlled Fluidization, Mobility, and Clogging in Obstacle Arrays Using Periodic Perturbations

We show that the clogging susceptibility and flow of particles moving through a random obstacle array can be controlled with a transverse or longitudinal ac drive. The flow rate can vary over several orders of magnitude, and we find both an optimal frequency and an optimal amplitude of driving that maximizes the flow. For dense arrays, at low ac frequencies, a heterogeneous creeping clogged phase appears in which rearrangements between different clogged configurations occur. At intermediate frequencies, a high-mobility fluidized state forms, and, at high frequencies, the system reenters a heterogeneous frozen clogged state. These results provide a technique for optimizing flow through heterogeneous media that could also serve as the basis for a particle separation method.

Laning and clustering transitions in driven binary active matter systems

It is well known that a binary system of nonactive disks that experience driving in opposite directions exhibits jammed, phase separated, disordered, and laning states. In active matter systems, such as a crowd of pedestrians, driving in opposite directions is common and relevant, especially in conditions which are characterized by high pedestrian density and emergency. In such cases, the transition from laning to disordered states may be associated with the onset of a panic state. We simulate a laning system containing active disks that obey run-and-tumble dynamics, and we measure the drift mobility and structure as a function of run length, disk density, and drift force. The activity of each disk can be quantified based on the correlation timescale of the velocity vector. We find that in some cases, increasing the activity can increase the system mobility by breaking up jammed configurations; however, an activity level that is too high can reduce the mobility by increasing the probability of disk-disk collisions. In the laning state, the increase of activity induces a sharp transition to a disordered strongly fluctuating state with reduced mobility. We identify a novel drive-induced clustered laning state that remains stable even at densities below the activity-induced clustering transition of the undriven system. We map out the dynamic phase diagrams highlighting transitions between the different phases as a function of activity, drive, and density.

Optimized Diffusion of Run-and-Tumble Particles in Crowded Environments

We study the transport of self-propelled particles in dynamic complex environments. To obtain exact results, we introduce a model of run-and-tumble particles (RTPs) moving in discrete time on a d-dimensional cubic lattice in the presence of diffusing hard-core obstacles. We derive an explicit expression for the diffusivity of the RTP, which is exact in the limit of low density of fixed obstacles. To do so, we introduce a generalization of Kac's theorem on the mean return times of Markov processes, which we expect to be relevant for a large class of lattice gas problems. Our results show the diffusivity of RTPs to be nonmonotonic in the tumbling probability for low enough obstacle mobility. These results prove the potential for the optimization of the transport of RTPs in crowded and disordered environments with applications to motile artificial and biological systems.

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.

Statistical analysis of phase formation in 2D colloidal systems

Colloidal systems offer unique opportunities for the study of phase formation and structure since their characteristic length scales are accessible to visible light. As a model system the two-dimensional assembly of colloidal magnetic and non-magnetic particles dispersed in a ferrofluid (FF) matrix is studied by transmission optical microscopy. We present a method to statistically evaluate images with thousands of particles and map phases by extraction of local variables. Different lattice structures and long-range connected branching chains are observed, when tuning the effective magnetic interaction and varying particle ratios.

Negative differential mobility and trapping in active matter systems

Using simulations, we examine the average velocity as a function of applied drift force for active matter particles moving through a random obstacle array. We find that for low drift force, there is an initial flow regime where the mobility increases linearly with drive, while for higher drift forces a regime of negative differential mobility appears in which the velocity decreases with increasing drive due to the trapping of active particles behind obstacles. A fully clogged regime exists at very high drift forces when all the particles are permanently trapped behind obstacles. We find for increasing activity that the overall mobility is nonmonotonic, with an enhancement of the mobility for small levels of activity and a decrease in mobility for large activity levels. We show how these effects evolve as a function of disk and obstacle density, active run length, drift force, and motor force.