Aggregation and segregation of confined active particles

The narrow escape problem of a chiral active particle (CAP): an optimal scheme

We report a simulation study on the narrow escape kinetics of a chiral active particle (CAP) confined to a circular domain with a narrow escape opening. The study's main objective is to optimize the CAP's escape chances as a function of the relevant parameters, such as translational and rotational speeds of the CAP, domain size, etc. We identified three regimes in the escape kinetics, namely the noise-dominated regime, the optimal regime, and the chiral activity-dominated regime. In particular, the optimal regime is characterized by an escape scheme that involves a direct passage to the domain boundary at first and then a unidirectional drift along the boundary towards the exit. Furthermore, we propose a non-dimensionalization approach to optimize the escape performance across microorganisms with varying motile characteristics. Additionally, we explore the influence of the translational and rotational noise on the CAP's escape kinetics.

Demixing in Binary Mixtures with Differential Diffusivity at High Density

Spontaneous phase separation, or demixing, is important in biological phenomena such as cell sorting. In particle-based models, an open question is whether differences in diffusivity can drive such demixing. While differential-diffusivity-induced phase separation occurs in mixtures with a packing fraction up to 0.7 [S. N. Weber et al. Binary mixtures of particles with different diffusivities demix, Phys. Rev. Lett. 116, 058301 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.058301], here we investigate whether demixing persists at even higher densities relevant for cells. For particle packing fractions between 0.7 and 1.0 the system demixes, but at packing fractions above unity the system remains mixed, exposing re-entrant behavior in the phase diagram that occurs when phase separation can no longer drive a change in entropy production at high densities. We also find that a confluent Voronoi model for tissues does not phase separate, consistent with particle-based simulations.

Nematic Torques in Scalar Active Matter: When Fluctuations Favor Polar Order and Persistence

We study the impact of nematic alignment on scalar active matter in the disordered phase. We show that nematic torques control the emergent physics of particles interacting via pairwise forces and can either induce or prevent phase separation. The underlying mechanism is a fluctuation-induced renormalization of the mass of the polar field that generically arises from nematic torques. The correlations between the fluctuations of the polar and nematic fields indeed conspire to increase the particle persistence length, contrary to what phenomenological computations predict. This effect is generic and our theory also quantitatively accounts for how nematic torques enhance particle accumulation along confining boundaries and opposes demixing in mixtures of active and passive particles.

An effective hydrodynamic description of marching locusts

A fundamental question in complex systems is how to relate interactions between individual components ('microscopic description') to the global properties of the system ('macroscopic description'). Furthermore, it is unclear whether such a macroscopic description exists and if such a description can capture large-scale properties. Here, we address the validity of a macroscopic description of a complex biological system using the collective motion of desert locusts as a canonical example. One of the world's most devastating insect plagues begins when flightless juvenile locusts form 'marching bands'. These bands display remarkable coordinated motion, moving through semiarid habitats in search of food. We investigated how well macroscopic physical models can describe the flow of locusts within a band. For this, we filmed locusts within marching bands during an outbreak in Kenya and automatically tracked all individuals passing through the camera frame. We first analyzed the spatial topology of nearest neighbors and found individuals to be isotropically distributed. Despite this apparent randomness, a local order was observed in regions of high density in the radial distribution function, akin to an ordered fluid. Furthermore, reconstructing individual locust trajectories revealed a highly aligned movement, consistent with the one-dimensional version of the Toner-Tu equations, a generalization of the Navier-Stokes equations for fluids, used to describe the equivalent macroscopic fluid properties of active particles. Using this effective Toner-Tu equation, which relates the gradient of the pressure to the acceleration, we show that the effective 'pressure' of locusts increases as a linear function of density in segments with the highest polarization (for which the one-dimensional approximation is most appropriate). Our study thus demonstrates an effective hydrodynamic description of flow dynamics in plague locust swarms.

Emergent structures in active block copolymer composites

Block copolymer melts offer unique templates to control the position and alignment of nanoparticles due to their ability to self-assemble into periodic ordered structures. Active particles are shown to coassemble with block copolymers leading to emergent organized structures. The block copolymer acts as a soft template that can control the self-propulsion of active particles, both for interface-segregated and selective nanoparticles. At moderate activities, active particles can form organized structures such as polarized trains or rotating vortices. At high activity, the contrast in the polymeric and colloidal timescales can lead to particle swarms with distorted block copolymer morphology, due to the competition between polymeric self-assembly and active Brownian self-propulsion.

Mechanochemical active ratchet

Self-propelled nanoparticles moving through liquids offer the possibility of creating advanced applications where such nanoswimmers can operate as artificial molecular-sized motors. Achieving control over the motion of nanoswimmers is a crucial aspect for their reliable functioning. While the directionality of micron-sized swimmers can be controlled with great precision, steering nano-sized active particles poses a real challenge. One of the reasons is the existence of large fluctuations of active velocity at the nanoscale. Here, we describe a mechanism that, in the presence of a ratchet potential, transforms these fluctuations into a net current of active nanoparticles. We demonstrate the effect using a generic model of self-propulsion powered by chemical reactions. The net motion along the easy direction of the ratchet potential arises from the coupling of chemical and mechanical processes and is triggered by a constant, transverse to the ratchet, force. The current magnitude sensitively depends on the amplitude and the periodicity of the ratchet potential and the strength of the transverse force. Our results highlight the importance of thermodynamically consistent modeling of chemical reactions in active matter at the nanoscale and suggest new ways of controlling dynamics in such systems.

Effect of confinement and topology: 2-TIPS vs. MIPS

2-TIPS (two temperature induced phase separation) refers to the phase separation phenomenon observed in mixtures of active and passive particles which are modelled using scalar activity. The active particles are connected to a thermostat at high temperature while the passive particles are connected to the thermostat at low temperature and the relative temperature difference between "hot" and "cold" particles is taken as the measure of the activity χ of the non-equilibrium system. The study of such binary mixtures of hot and cold particles under various kinds of confinement is an important problem in many physical and biological processes. The nature and extent of phase separation are heavily influenced by the geometry of confinement, activity, and density of the non-equilibrium binary mixture. Investigating such 3D binary mixtures confined by parallel walls, we observe that the active and passive particles phase separate, but the extent of phase separation is reduced compared to bulk phase separation at high densities and enhanced at low densities. However, when the binary mixture of active and passive particles is confined inside a spherical cavity, the phase separation is radial for small radii of the confining sphere and the extent of phase separation is higher compared to their bulk counterparts. Confinement leads to interesting properties in the passive (cold) region like enhanced layering and high compression in the direction parallel to the confining wall. In 2D, both the bulk and confined systems of the binary mixture show a significant decrement in the extent of phase separation at higher densities. This observation is attributed to the trapping of active particles inside the passive cluster, which increases with density. Thus the 2D systems show structures more akin to dense-dilute phase co-existence, which is observed in motility induced phase separation in 2D active systems. The binary mixture constrained on the spherical surface also shows similar phase co-existence. Our analyses reveal that the coexistent densities observed in 2-TIPS on the spherical surface agree with the findings of previous studies on MIPS in active systems on a sphere.

Motility induced phase separation of deformable cells

Using a multi-phase field model, we examine how particle deformability, which is a proxy for cell stiffness, affects motility induced phase separation (MIPS). We show that purely repulsive deformable, i.e., squishy, cells phase separate more effectively than their rigid counterparts. This can be understood as due to the fact that deformability increases the effective duration of collisions. In addition, the dense regions become increasingly disordered as deformability increases. Our results contextualize the applicability of MIPS to biological systems and have implications for how cells in biological systems may self-organize.

Effective temperature for an intermittent bistable potential

Thermodynamics of far-from-equilibrium systems often require measurement of effective parameters such as temperature. Whether such approach is valid for the general case of resetting protocols, active systems, or of confined systems under time-varying fields is still under investigation. We report on the effect of switching ON-OFF of an asymmetric bistable potential to the mean first passage time (MFPT) of a probed particle to go from one potential minima to the other. Experimental results coupled with numerical simulations shows the potential becoming more symmetric at slow switching. Moreover, the MFPT deviates from equilibrium condition with an effective temperature, Teff < T, at slow switching but approaches room temperature, T, at fast switching. For each switching rate, we quantify how far the system is from equilibrium by measuring deviation from a detailed balance like relation and the net circulation of flux present in phase-space. Both analysis suggest equilibrium condition are met at high switching.

Mixtures of self-propelled particles interacting with asymmetric obstacles

In the presence of an obstacle, active particles condensate into a surface "wetting" layer due to persistent motion. If the obstacle is asymmetric, a rectification current arises in addition to wetting. Asymmetric geometries are therefore commonly used to concentrate microorganisms like bacteria and sperms. However, most studies neglect the fact that biological active matter is diverse, composed of individuals with distinct self-propulsions. Using simulations, we study a mixture of "fast" and "slow" active Brownian disks in two dimensions interacting with large half-disk obstacles. With this prototypical obstacle geometry, we analyze how the stationary collective behavior depends on the degree of self-propulsion "diversity," defined as proportional to the difference between the self-propulsion speeds, while keeping the average self-propulsion speed fixed. A wetting layer rich in fast particles arises. The rectification current is amplified by speed diversity due to a superlinear dependence of rectification on self-propulsion speed, which arises from cooperative effects. Thus, the total rectification current cannot be obtained from an effective one-component active fluid with the same average self-propulsion speed, highlighting the importance of considering diversity in active matter.

Steady-state moments under resetting to a distribution

The nonequilibrium steady state emerging from stochastic resetting to a distribution is studied. We show that for a range of processes, the steady-state moments can be expressed as a linear combination of the moments of the distribution of resetting positions. The coefficients of this series are universal in the sense that they do not depend on the resetting distribution, only the underlying dynamics. We consider the case of a Brownian particle and a run-and-tumble particle confined in a harmonic potential, where we derive explicit closed-form expressions for all moments for any resetting distribution. Numerical simulations are used to verify the results, showing excellent agreement.

Active Brownian particles in external force fields: Field-theoretical models, generalized barometric law, and programmable density patterns

We investigate the influence of external forces on the collective dynamics of interacting active Brownian particles in two as well as three spatial dimensions. Via explicit coarse graining, we derive predictive models, i.e., models that give a direct relation between the models' coefficients and the bare parameters of the system, that are applicable for space- and time-dependent external force fields. We study these models for the cases of gravity and harmonic traps. In particular, we derive a generalized barometric formula for interacting active Brownian particles under gravity that is valid for low to high concentrations and activities of the particles. Furthermore, we show that one can use an external harmonic trap to induce motility-induced phase separation in systems that, without external fields, remain in a homogeneous state. This finding makes it possible to realize programmable density patterns in systems of active Brownian particles. Our analytic predictions are found to be in very good agreement with Brownian dynamics simulations.

Free and enclosed inertial active gas

In this work, the free expansion of an inertial active gas in three dimensions made of spherical non-interactive active Brownian particles with both translational and rotational inertia (IABPs) is studied. After elucidating the active particles' orientational correlation in three dimensions by employing a Fokker-Planck formalism, the diffusion, mean-square speed, persistence length, reorientation time, Swim and Reynolds pressures and total pressure of this system, are obtained theoretically and corroborated by performing Langevin dynamics simulations. Afterwards, a numerical study on particles' distribution and the mechanical pressure exerted by the active gas enclosed in a cubic box and its dependence on inertia is also carried out. This experiment highlights two important observations: first, as inertia in the system grows while fixing activity, a more uniform particle distribution within the box is achieved. In other words, the classical accumulation of active particles at the walls is seen to be suppressed by inertia. Second, an active gas with translational and rotational inertiae and made of spherical particles still has a state equation which is offered here. This is supported by the fact that both the mechanical pressure definition and the bulk pressure definition as the trace of the swim and Reynolds stress tensors, coincide in the thermodynamic limit.

Stationary states of an active Brownian particle in a harmonic trap

We study the stationary states of an overdamped active Brownian particle (ABP) in a harmonic trap in two dimensions via mathematical calculations and numerical simulations. In addition to translational diffusion, the ABP self-propels with a certain velocity, whose magnitude is constant, but its direction is subject to Brownian rotation. In the limit where translational diffusion is negligible, the stationary distribution of the particle's position shows a transition between two different shapes, one with maximum and the other with minimum density at the center, as the trap stiffness is increased. We show that this nonintuitive behavior is captured by the relevant Fokker-Planck equation, which, under minimal assumptions, predicts a continuous phase transition-like change between the two different shapes. As the translational diffusion coefficient is increased, both these distributions converge into the equilibrium, Boltzmann form. Our simulations support the analytical predictions and also show that the probability distribution of the orientation angle of the self-propulsion velocity undergoes a transition from unimodal to bimodal forms in this limit. We also extend our simulations to a three-dimensional trap and find similar behavior.

Phase behavior and dynamics in a colloid-polymer mixture under spherical confinement

From studies via molecular dynamics simulations, we report results on structure and dynamics in mixtures of active colloids and passive polymers that are confined inside a spherical container with a repulsive boundary. All interactions in the fully passive limit are chosen in such a way that in equilibrium coexistence between colloid-rich and polymer-rich phases occurs. For most part of the studies the chosen compositions give rise to Janus-like structure: nearly one side of the sphere is occupied by the colloids and the rest by the polymers. This partially wet situation mimics approximately a neutral wall in the fully passive scenario. Following the introduction of a velocity-aligning activity to the colloids, the shape of the polymer-rich domain changes to that of an ellipsoid, around the long axis of which the colloid-rich domain attains a macroscopic angular momentum. In the steady state, the orientation of this axis evolves via diffusion, magnitude of which depends upon the strength of activity, but only weakly.

Collective behavior of chiral active particles with anisotropic interactions in a confined space

Extensive studies so far have indicated that chirality, anisotropic interactions and spatial confinement play important roles in collective dynamics in active matter systems. However, how the overall interplay of these crucial factors affects the novel phases and macroscopic properties remains less explored. Here, using Langevin dynamics simulations, we investigate the self-organization of a chiral active system composed of amphiphilic Janus particles, where the embedded anisotropic interaction orientation is assumed to be either the same or just opposite to the direction of active force. A wealth of dynamic phases are observed including formation of phase separation, clustering state, homogeneous state, spiral vortex flow, swarm and spatiotemporal oscillation. By tuning self-propelled angular speed and anisotropic interaction strength, we identify the non-equilibrium phase diagrams, and reveal the very non-trivial modulation of both vortex and swarm patterns. Intriguingly, we find that strong chirality-alignment-confinement coupling yields a self-driven spatial and temporal organization periodically oscillating between a counterclockwise vortex and a clockwise one. Our work provides a new understanding of the novel self-assembly arising in such a confined system and enables new strategies for achieving ordered dynamic structures.

Binary Mixtures in Linear Convection Arrays

We numerically investigated the dynamics of a mixture of finite-size active and passive disks in a linear array of two-dimensional convection rolls. The interplay of advection and steric interactions produces a number of interesting effects, like the stirring of a passive colloidal fluid by a small fraction of slow active particles, or the separation of the mixture active and passive colloidal fractions by increasing the motility of the active one, which eventually clusters in stagnation areas along the array walls. These mechanisms are quantitatively characterized by studying the dependence of the diffusion constants of the active and passive particles on the parameters of the active mixture fraction.

Directed Autonomous Motion and Chiral Separation of Self-Propelled Janus Particles in Convection Roll Arrays

Self-propelled Janus particles exhibit autonomous motion thanks to engines of their own. However, due to the randomly changing direction of such motion they are of little use for emerging nanotechnological and biomedical applications. Here, we numerically show that the motion of chiral active Janus particles can be directed, subjecting them to a linear array of convection rolls. The rectification power of self-propulsion motion here can be made to be more than 60%, which is much larger than earlier reports. We show that rectification of a chiral Janus particle's motion leads to conspicuous segregation of dextrogyre and levogyre active particles from a racemic binary mixture. Further, we demonstrate how efficiently the rectification effect can be exploited to separate dextrogyre and levogyre particles when their intrinsic torques are distributed with Gaussian statistics.

Cell segregation via differential collision modes between heterotypic cell populations

In tissue development and regeneration, the establishment of sharp boundaries between heterotypic cells is essential for the differentiation of tissue functions. During the dynamic rearrangements of constituent cells that result from cell division and collective migration, the segregation boundary encounters various challenges. Several studies have suggested that cortical actomyosin structures play a crucial role in the maintenance of the boundary interface of segregated cell populations, implicating actin-mediated stresses. Examining physical cellular properties such as motility, traction, and intercellular stress, we investigated the formation and maintenance of the stable segregation between epithelial and mesenchymal cell populations devoid of heterotypic adhesions. At the contact boundary, the homotypic adhesion-mediated epithelial aggregates exerted collision-mediated compression against the surrounding mesenchymal cells. Our results demonstrated that heterotypic cell populations established a robust interfacial boundary by accumulating stress from active collisions and repulsions between two dissimilar cell types. Furthermore, the moment of the heterotypic collisions was identified by the existence of a sharp rise in maximum shear stress within the cell cluster.

Motility-induced phase separation of self-propelled soft inertial disks

The phase diagram of the phenomenon of motility-induced phase separation (MIPS) for a collection of self-propelled interacting disks over a large inertial range is explored using active Langevin dynamics simulation with particular emphasis on disk softness and effective size. It is shown that the parabola-like phase boundary between the homogeneous and MIPS states in the semi-log space of disk softness and effective size moves towards the hard disk limit with increase in inertia, before complete disappearance in the limit of large inertia. With increase in effective size of the disks, re-entrant phase separation, that is the system phase-separating from a homogeneous phase and eventually re-entering the homogeneous phase, is reported. The structural and the dynamical properties of the different phases are investigated in the considered inertial range. The particular shape of the phase boundary and the re-entrant behavior is explained based on several qualitative and quantitative results. Unlike most of the earlier studies on MIPS, which consider hard particle limits, our findings may be directly applicable to soft active matter for a range of physical and biological systems.

Passive objects in confined active fluids: A localization transition

We study how walls confining active fluids interact with asymmetric passive objects placed in their bulk. We show that the objects experience nonconservative long-ranged forces mediated by the active bath. To leading order, these forces can be computed using a generalized image theorem. The walls repel asymmetric objects, irrespective of their microscopic properties or their orientations. For circular cavities, we demonstrate how this may lead to the localization of asymmetric objects in the center of the cavity, something impossible for symmetric ones.

Activity Induced Rigidity of Liquid Droplets

Here we show that encapsulating active Janus particles within a drop renders it more resistant to deformation. This drop is deformed under the action of an extensional flow. Such deformation is primarily resisted by the drop interfacial tension. When the particles are active under the action of laser illumination, the deformation decreases signaling an increase in effective tension or Laplace pressure. This increase is attributed to the activity of the particles. Our results using numerous drop sizes, particle number densities, and active velocities show that the obtained increase agrees surprisingly well, over an extended range, with a standard expression for the pressure engendered by an ensemble of active particles, proposed years ago but not tested yet in three dimensions.

Diffusion coefficient and power spectrum of active particles with a microscopically reversible mechanism of self-propelling

Catalytically active macromolecules are envisioned as key building blocks in development of artificial nanomotors. However, theory and experiments report conflicting findings regarding their dynamics. The lack of consensus is mostly caused by a limited understanding of specifics of self-propulsion mechanisms at the nanoscale. Here, we study a generic model of a self-propelled nanoparticle that does not rely on a particular mechanism. Instead, its main assumption is the fundamental symmetry of microscopic dynamics of chemical reactions: the principle of microscopic reversibility. Significant consequences of this assumption arise if we subject the particle to an action of an external time-periodic force. The particle diffusion coefficient is then enhanced compared to the unbiased dynamics. The enhancement can be controlled by the force amplitude and frequency. We also derive the power spectrum of particle trajectories. Among new effects stemming from the microscopic reversibility are the enhancement of the spectrum at all frequencies and sigmoid-shaped transitions and a peak at characteristic frequencies of rotational diffusion and external forcing. The microscopic reversibility is a generic property of a broad class of chemical reactions, therefore we expect that the presented results will motivate new experimental studies aimed at testing of our predictions. This could provide new insights into dynamics of catalytic macromolecules.

Mechanical feedback controls the emergence of dynamical memory in growing tissue monolayers

The growth of a tissue, which depends on cell–cell interactions and biologically relevant processes such as cell division and apoptosis, is regulated by a mechanical feedback mechanism. We account for these effects in a minimal two-dimensional model in order to investigate the consequences of mechanical feedback, which is controlled by a critical pressure, p c. A cell can only grow and divide if its pressure, due to interaction with its neighbors, is less than p c. Because temperature is not a relevant variable, the cell dynamics is driven by self-generated active forces (SGAFs) that arise due to cell division. We show that even in the absence of intercellular interactions, cells undergo diffusive behavior. The SGAF-driven diffusion is indistinguishable from the well-known dynamics of a free Brownian particle at a fixed finite temperature. When intercellular interactions are taken into account, we find persistent temporal correlations in the force–force autocorrelation function (FAF) that extends over a timescale of several cell division times. The time-dependence of the FAF reveals memory effects, which increases as p c increases. The observed non-Markovian effects emerge due to the interplay of cell division and mechanical feedback and are inherently a non-equilibrium phenomenon.

Coexisting orbits and chaotic dynamics of a confined self-propelled particle

We investigate theoretically the dynamics of a confined active swimmer with velocity and orientation axis coupled to each other via a self-alignment torque. For an isotropic harmonic potential, this system is known to exhibit two distinct dynamical phases: a climbing one, where the particle is oriented radially and undergoes angular Brownian motion, and a circularly orbiting phase. Here we show that for nonradially symmetric confinement an assortment of complex phenomena emerge. For an elliptic harmonic potential the orbiting phase splits into several periodic orbits with a diversity of shapes: ovals, lemniscates, and generalized lemniscates with multiple lobes. These orbits can coexist in the parameter space and decay into one another induced by noise. For anharmonic confining potentials, we report transitions from periodic to chaotic dynamics, as one changes the intensity of the self-alignment torque and noise-induced complex orbits. These results demonstrate that the combination of the shape of the trapping potential and self-alignment torque can induce a rich variety of nontrivial dynamical states of a confined active particle.

Pattern formation and phase transition in the collective dynamics of a binary mixture of polar self-propelled particles

The collective behavior of a binary mixture of polar self-propelled particles (SPPs) with different motile properties is studied. The binary mixture consists of slow-moving SPPs (sSPPs) of fixed velocity v_{s} and fast-moving SPPs (fSPPs) of fixed velocity v_{f}. These SPPs interact via a short-range interaction irrespective of their types. They move following certain position and velocity update rules similar to the Vicsek model (VM) under the influence of an external noise η. The system is studied at different values of v_{f} keeping v_{s}=0.01 constant for a fixed density ρ=0.5. Different phase-separated collective patterns that appear in the system over a wide range of noise η are characterized. The fSPPs and the sSPPs are found to be orientationally phase synchronized at the steady state. We studied an orientational order-disorder transition varying the angular noise η and identified the critical noise η_{c} for different v_{f}. Interestingly, both the species exhibit continuous transition for v_{f}<100v_{s} and discontinuous transition for v_{f}>100v_{s}. A new set of critical exponents is determined for the continuous transitions. However, the binary model is found to be nonuniversal as the values of the critical exponents depend on the velocity. The effect of interaction radius on the system behavior is also studied.

Disordered boundaries destroy bulk phase separation in scalar active matter

We show that disordered boundaries destroy bulk phase separation in scalar active systems in dimension d<d_{c}=3. This is in strong contrast with the equilibrium case where boundaries have no impact on the bulk of phase-separated systems. The underlying mechanism is revealed by considering a localized deformation of an otherwise flat wall, from which the case of a disordered boundary can be inferred. We find long-ranged correlations of the density field as well as a cascade of eddies which we show prevent bulk phase separation in low enough dimensions. The results are derived for dilute systems as well as in the presence of interactions, under the sole condition that the density field is the unique hydrodynamic mode. Our theoretical calculations are validated by numerical simulations of microscopic active systems.

Searching for structural predictors of plasticity in dense active packings

In amorphous solids subject to shear or thermal excitation, so-called structural indicators have been developed that predict locations of future plasticity or particle rearrangements. An open question is whether similar tools can be used in dense active materials, but a challenge is that under most circumstances, active systems do not possess well-defined solid reference configurations. We develop a computational model for a dense active crowd attracted to a point of interest, which does permit a mechanically stable reference state in the limit of infinitely persistent motion. Previous work on a similar system suggested that the collective motion of crowds could be predicted by inverting a matrix of time-averaged two-particle correlation functions. Seeking a first-principles understanding of this result, we demonstrate that this active matter system maps directly onto a granular packing in the presence of an external potential, and extend an existing structural indicator based on linear response to predict plasticity in the presence of noisy dynamics. We find that the strong pressure gradient necessitated by the directed activity, as well as a self-generated free boundary, strongly impact the linear response of the system. In low-pressure regions the linear-response-based indicator is predictive, but it does not work well in the high-pressure interior of our active packings. Our findings motivate and inform future work that could better formulate structure-dynamics predictions in systems with strong pressure gradients.

Active matter dynamics in confined microfluidic environments

The field of active matter is a nascent area of research in soft condensed matter physics, which is drawing on the expertise of researchers from diverse disciplines. Small scale active particles-both inorganic and biological-display non-trivial emergent dynamics and interactions that could help us understand complex biological processes and phenomena. Recently, using microfluidic technologies, several research groups have performed important experimental and theoretical studies to understand the behavior of self-propelled particles and molecular active matter within confined environments-to glean a fundamental understanding of the cellular processes occurring under ultra-low Reynolds number conditions. In this chapter, we would like to review applications of microfluidics in active matter research, highlighting a few important theoretical and experimental investigations. We will conclude the discussion with a note on the future of this field mentioning a few open questions that are at the forefront of our minds.

Escape kinetics of self-propelled particles from a circular cavity

We numerically investigate the mean exit time of an inertial active Brownian particle from a circular cavity with single or multiple exit windows. Our simulation results witness distinct escape mechanisms depending on the relative amplitudes of the thermal length and self-propulsion length compared to the cavity and pore sizes. For exceedingly large self-propulsion lengths, overdamped active particles diffuse on the cavity surface, and rotational dynamics solely governs the exit process. On the other hand, the escape kinetics of a very weakly damped active particle is largely dictated by bouncing effects on the cavity walls irrespective of the amplitude of self-propulsion persistence lengths. We show that the exit rate can be maximized for an optimal self-propulsion persistence length, which depends on the damping strength, self-propulsion velocity, and cavity size. However, the optimal persistence length is insensitive to the opening windows' size, number, and arrangement. Numerical results have been interpreted analytically based on qualitative arguments. The present analysis aims at understanding the transport controlling mechanism of active matter in confined structures.

Alignment-mediated segregation in an active-passive mixture

We report segregation between the athermal active and passive particles mediated by the local alignment interaction in a confined space. The competition between the alignment interaction and self-propulsion force results in a transition between disordered and ordered phases. We show that as the coordination between the particles increases, they form an ordered mill, which helps the particles to aggregate into isotropic clusters. As a result, particles segregate into active core and passive shells. This segregation phenomenon is adversely affected by the packing fraction and the size dispersion between active and passive particles. We show that this adverse effect can be overcome by incorporating higher coordination in the system. We report that the monodispersed system is more desirable for segregation in a binary mixture than a bidispersed system, as the latter favors the mixed state.

Dependency of active pressure and equation of state on stiffness of wall

Autonomous motion and motility are hallmarks of active matter. Active agents, such as biological cells and synthetic colloidal particles, consume internal energy or extract energy from the environment to generate self-propulsion and locomotion. These systems are persistently out of equilibrium due to continuous energy consumption. It is known that pressure is not always a state function for generic active matter. Torque interaction between active constituents and confinement renders the pressure of the system a boundary-dependent property. The mechanical pressure of anisotropic active particles depends on their microscopic interactions with a solid wall. Using self-propelled dumbbells confined by solid walls as a model system, we perform numerical simulations to explore how variations in the wall stiffness influence the mechanical pressure of dry active matter. In contrast to previous findings, we find that mechanical pressure can be independent of the interaction of anisotropic active particles with walls, even in the presence of intrinsic torque interaction. Particularly, the dependency of pressure on the wall stiffness vanishes when the stiffness is above a critical level. In such a limit, the dynamics of dumbbells near the walls are randomized due to the large torque experienced by the dumbbells, leading to the recovery of pressure as a state variable of density.

Morphologies and dynamics of the interfaces between active and passive phases

Active matters exhibit interesting collective behaviors and novel phases, which provide an important platform for the study of nonequilibrium physics. Mixtures of active and passive particles have been intensively studied in motility-induced phase separation, but the morphology of the active-passive interface has been poorly explored. In this work, we investigate the interface morphology in two-dimensional mixtures of active and passive particles using Brownian dynamics simulations. By systematically changing the Péclet number (Pe) and area fraction (ρ), we obtain the phase diagram of the active-passive interface, including rough sharp, rough invasive and flat interdiffusive interfaces. For a sharp interface, dynamic scaling analysis in the propagation stage shows that the roughness exponent α, the growth exponent β, the time exponent κ, and the dynamic exponent z satisfy z = α/(β - κ). Such anomalous scaling indicates that the roughening behavior does not belong to the conventional universality classes with Family-Vicsek scaling for the growth of passive interfaces. On the other hand, the interface in the middle-wavelength regime during the morphology relaxation stage can be described by capillary wave theory. The mean interface position propagates with time as t^1/2, which is robust at different ρ and Pe values in the propagation stage and exhibits superdiffusion in the morphology relaxation stage. These similarities and differences between the active-inactive interfaces and passive interfaces cast light on the interfacial growth of active matter.

Emergent spiral vortex of confined biased active particles

Confinement is known to have profound effects on the collective dynamics of many active systems. Here, we investigate a modeled active system in circular confinement consisting of biased active particles, where the direction of active force deviates a biased angle from the principle orientation of the anisotropic interaction. We find that such particles can spontaneously form a spiral vortex with two concentric and counter-rotating regions near the boundary. The emerged vortex can be measured by the vortex order parameter which shows nonmonotonic dependencies on both the biased angle and the strength of the anisotropic interaction. Our work can provide an understanding of such dynamic behaviors and enable different strategies for designing ordered collective behaviors.

Absorption of self-propelled particles into a dense porous medium

We study the absorption of self-propelled particles into a finite-size dense porous medium, which is mimicked by an obstacle array. We find that, depending on the competition of the propelling strength versus the repulsive barrier formed by obstacles and the contrast between the characteristic time scales of permeation and propelling persistence, the absorption process exhibits three distinct types of behavior. In Type I and II behavior, the propelling strength is not large enough to surmount the barrier, and hence particles transport in the medium by barrier-hopping dynamics. The initial permeation of particles toward the medium center is phenomenologically similar to a normal slow diffusion process. But, surprisingly, after the initial permeation process, a concentrated nucleus of particle aggregates forms and grows at the medium center in Type I, due to the long propelling persistence. Such an abnormal "nucleation" phenomenon does not appear in Type II, in which the propelling persistence is low. When the propelling strength is very high (Type III), particles transport smoothly in the medium, hence the initial slow diffusion process disappears and small particle clusters form and merge randomly in the medium. Our results provide a foundation for applications of active objects in a complex environment and also suggest the possible usage of a porous medium, for example, in the selection or sorting of active matter.

Cluster and conquer: the morphodynamics of invasion of a compliant substrate by active rods

The colonisation of a soft passive material by motile cells such as bacteria is common in biology. The resulting colonies of the invading cells are often observed to exhibit intricate patterns whose morphology and dynamics can depend on a number of factors, particularly the mechanical properties of the substrate and the motility of the individual cells. We use simulations of a minimal 2D model of self-propelled rods moving through a passive compliant medium consisting of particles that offer elastic resistance before being plastically displaced from their equilibrium positions. It is observed that the clustering of active (self-propelled) particles is crucial for understanding the morphodynamics of colonisation. Clustering enables motile colonies to spread faster than they would have as isolated particles. The colonisation rate depends non-monotonically on substrate stiffness with a distinct maximum at a non-zero value of substrate stiffness. This is observed to be due to a change in the morphology of clusters. Furrow networks created by the active particles have a fractal-like structure whose dimension varies systematically with substrate stiffness but is less sensitive to particle activity. The power-law growth exponent of the furrowed area is smaller than unity, suggesting that, to sustain such extensive furrow networks, colonies must regulate their overall growth rate.

Rotation reversal of a ratchet gear powered by active particles

Rotation of a gear powered by active particles is numerically investigated in a circular chamber. Due to the nonequilibrium properties of active particles, net gear rotation is achieved in a bath composed of self-propelling particles. Our setup can convert the random motion of active particles into the directional rotation of the ratchet gear. The direction of rotation is determined by the asymmetry of the gear and the persistence length (the ratio of the self-propulsion speed to the rotation diffusion coefficient) of active particles. Remarkably, the direction of rotation for large persistence length is opposite to the direction of rotation for small persistence length. Therefore, for a given asymmetric gear, we can observe the rotation reversal when tuning the system parameters (e.g., the self-propulsion speed, the rotation diffusion coefficient, and the packing fraction of active particles). Our findings are relevant to the experimental pursuit of rectifying random motion to directional motion in active matter.

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.

Statistical mechanics of active Ornstein-Uhlenbeck particles

We study the statistical properties of active Ornstein-Uhlenbeck particles (AOUPs). In this simplest of models, the Gaussian white noise of overdamped Brownian colloids is replaced by a Gaussian colored noise. This suffices to grant this system the hallmark properties of active matter, while still allowing for analytical progress. We study in detail the steady-state distribution of AOUPs in the small persistence time limit and for spatially varying activity. At the collective level, we show AOUPs to experience motility-induced phase separation both in the presence of pairwise forces or due to quorum-sensing interactions. We characterize both the instability mechanism leading to phase separation and the resulting phase coexistence. We probe how, in the stationary state, AOUPs depart from their thermal equilibrium limit by investigating the emergence of ratchet currents and entropy production. In the small persistence time limit, we show how fluctuation-dissipation relations are recovered. Finally, we discuss how the emerging properties of AOUPs can be characterized from the dynamics of their collective modes.

Tunable depletion force in active and crowded environments

We adopt two-dimensional Langevin dynamics simulations to study the effective interactions between two passive colloids in a bath crowded with active particles. We mainly pay attention to the significant effects of active particle size, crowding-activity coupling, and chirality. First, a transition of depletion force from repulsion to attraction is revealed by varying particle size. Moreover, larger active crowders with sufficient activity can generate strong attractive force, which is in contrast to the cage effect in passive media. It is interesting that the attraction induced by large active crowders follows a linear scaling with the persistence length of active particles. Second, the effective force also experiences a transition from repulsion to attraction as volume fraction increases, as a consequence of the competition between the two contrastive factors of activity and crowding. As bath volume fraction is relatively small, activity generates a dominant repulsion force, while as the bath becomes concentrated, crowding-induced attraction becomes overwhelming. Lastly, in a chiral bath, we observe a very surprising oscillation phenomenon of active depletion force, showing an evident quasiperiodic variation with increasing chirality. Aggregation of active particles in the vicinity of the colloids is carefully examined, which serves as a reasonable picture for our observations. Our findings provide an inspiring strategy for the tunable active depletion force by crowding, activity, and chirality.

Collective motion in large deviations of active particles

We analyze collective motion that occurs during rare (large deviation) events in systems of active particles, both numerically and analytically. We discuss the associated dynamical phase transition to collective motion, which occurs when the active work is biased towards larger values, and is associated with alignment of particles' orientations. A finite biasing field is needed to induce spontaneous symmetry breaking, even in large systems. Particle alignment is computed exactly for a system of two particles. For many-particle systems, we analyze the symmetry breaking by an optimal-control representation of the biased dynamics, and we propose a fluctuating hydrodynamic theory that captures the emergence of polar order in the biased state.

Theory for the Casimir effect and the partitioning of active matter

Active Brownian particles (ABPs) distribute non-homogeneously near surfaces, and understanding how this depends on system properties-size, shape, activity level, etc.-is essential for predicting and exploiting the behavior of active matter systems. Active particles accumulate at no-flux surfaces owing to their persistent swimming, which depends on their intrinsic swim speed and reorientation time, and are subject to confinement effects when their run or persistence length is comparable to the characteristic size of the confining geometry. It has been observed in simulations that two parallel plates experience a "Casimir effect" and attract each other when placed in a dilute bath of ABPs. In this work, we provide a theoretical model based on the Smoluchowski equation and a macroscopic mechanical momentum balance to analytically predict this attractive force. We extend this method to describe the concentration partitioning of active particles between a confining channel and a reservoir, showing that the ratio of the concentration in the channel to that in the bulk increases as either run length increases or channel height decreases. The theoretical results agree well with Brownian dynamics simulations and finite element calculations.

Spatiotemporal dynamics of a self-propelled system with opposing alignment and repulsive forces

Effect of concurrent alignment and repulsion is studied in the purview of a confined active matter system using a modified force-based Vicsek model. On alteration of the alignment and the repulsive force parameters, a low alignment random phase, a midrange alignment milling phase, and a high alignment oscillatory phase are identified. Based on the particle aggregations, the milling phase is further classified into three subphases, two of which are spatial patterns: one consisting of compact ring-shaped mills and the other incorporating both rings and clusters. A correlation function based on the inner product of spatial velocity fluctuations of the particles shows a high correlation length for the ringed milling and the rings-clusters hybrid milling state. On analyzing temporal velocity fluctuations of particles through chaos detection techniques, low alignment and high alignment states are indicative of chaos, while the middle order alignment is symbolic of periodicity. The extent of synchronization of the particles' motion is analyzed through a Hilbert transform-based mean frequency approach, leading to the detection of a weak chimera state in the case of the spatial structures. The ringed milling state shows a unique category of weak chimera consisting of multiple oscillator groups showcasing different synchronization frequencies coexisting with desynchronized oscillators.

Modeling multicellular dynamics regulated by extracellular-matrix-mediated mechanical communication via active particles with polarized effective attraction

Collective cell migration is crucial to many physiological and pathological processes such as embryo development, wound healing, and cancer invasion. Recent experimental studies have indicated that the active traction forces generated by migrating cells in a fibrous extracellular matrix (ECM) can mechanically remodel the ECM, giving rise to bundlelike mesostructures bridging individual cells. Such fiber bundles also enable long-range propagation of cellular forces, leading to correlated migration dynamics regulated by the mechanical communication among the cells. Motivated by these experimental discoveries, we develop an active-particle model with polarized effective attractions (APPA) to investigate emergent multicellular migration dynamics resulting from ECM-mediated mechanical communications. In particular, the APPA model generalizes the classic active-Brownian-particle (ABP) model by imposing a pairwise polarized attractive force between the particles, which depends on the instantaneous dynamic states of the particles and mimics the effective mutual pulling between the cells via the fiber bundle bridge. The APPA system exhibits enhanced aggregation behaviors compared to the classic ABP system, and the contrast is more apparent at lower particle densities and higher rotational diffusivities. Importantly, in contrast to the classic ABP system where the particle velocities are not correlated for all particle densities, the high-density phase of the APPA system exhibits strong dynamic correlations, which are characterized by the slowly decaying velocity correlation functions with a correlation length comparable to the linear size of the high-density phase domain (i.e., the cluster of particles). The strongly correlated multicellular dynamics predicted by the APPA model is subsequently verified in in vitro experiments using MCF-10A cells. Our studies indicate the importance of incorporating ECM-mediated mechanical coupling among the migrating cells for appropriately modeling emergent multicellular dynamics in complex microenvironments.

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.

Entropy production of an active particle in a box

A run-and-tumble particle in a one-dimensional box (infinite potential well) is studied. The steady state is analytically solved and analyzed, revealing the emergent length scale of the boundary layer where particles accumulate near the walls. The mesoscopic steady state entropy production rate of the system is derived from coupled Fokker-Planck equations with a linear reaction term, resulting in an exact analytic expression. The entropy production density is shown to peak at the walls. Additionally, the derivative of the entropy production rate peaks at a system size proportional to the length scale of the accumulation boundary layer, suggesting that the behavior of the entropy production rate and its derivatives as a function of the control parameter may signify a qualitative behavior change in the physics of active systems, such as phase transitions.

Aggregate morphology of active Brownian particles on porous, circular walls

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

Morphological transitions of active Brownian particle aggregates on porous walls

Motility-induced wall aggregation of Active Brownian Particles (ABPs) is a well-studied phenomenon. Here, we study the aggregation of ABPs on porous walls, which allows the particles to penetrate through at large motility. We show that the active aggregates undergo a morphological transition from a connected dense-phase to disconnected droplets with an increase in wall porosity and the particle self-motility, similar to wetting-dewetting transitions in equilibrium fluids. We show that both morphologically distinct states are stable, and independent of initial conditions at least in some parameter regions. Our analysis reveals that changes in wall porosity affect the intrinsic properties of the aggregates and changes the effective wall-aggregate interfacial tension, consistent with the appearance of the morphological transition. Accordingly, a close analysis of the density, as well as orientational distribution, indicates that the underlying reason for such morphological transitions is not necessarily specific to the systems with porous walls, and it can be possible to observe in a larger class of confined, active systems by tuning the properties of confining walls.

Testing a thermodynamic approach to collective animal behavior in laboratory fish schools

Collective behaviors displayed by groups of social animals are observed frequently in nature. Understanding and predicting the behavior of complex biological systems is dependent on developing effective descriptions and models. While collective animal systems are characteristically nonequilibrium, we can employ concepts from equilibrium statistical mechanics to motivate the measurement of material-like properties in laboratory animal aggregates. Here, we present results from a new set of experiments that utilize high speed footage of two-dimensional schooling events, particle tracking, and projected static and dynamic light fields to observe and control the behavior of negatively phototaxic fish schools (Hemigrammus bleheri). First, we use static light fields consisting of dark circular regions to produce visual stimuli that confine the schools to a range of areas. We find that schools have a maximum density which is independent of group size, and that a swim pressurelike quantity, Π increases linearly with number density, suggesting that unperturbed schools exist on an isotherm. Next, we use dynamic light fields where the radius of the dark region shrinks linearly with time to compress the schools. We find that an effective temperature parameter depends on the compression time and our results are thus consistent with the school having a constant heat flux. These findings further evidence the utility of effective thermodynamic descriptions of nonequilibrium systems in collective animal behavior.