Active Brownian equation of state: metastability and phase coexistence

Kinetic Theory of Motility Induced Phase Separation for Active Brownian Particles

When two active Brownian particles collide, they slide along each other until they can continue their free motion. For persistence lengths much larger than the particle diameter, the directors do not change, but the collision can be modeled as producing a net displacement on the particles compared to their free motion in the absence of the encounter. With these elements, a Boltzmann-Enskog-like kinetic theory is built. A linear stability analysis of the homogeneous state predicts a density instability resulting from the effective velocity reduction of tagged particles predicted by the theory.

Phase Coexistence and Edge Currents in the Chiral Lennard-Jones Fluid

We study a model chiral fluid in two dimensions composed of Brownian disks interacting via a Lennard-Jones potential and a nonconservative transverse force, mimicking colloids spinning at a given rate. The system exhibits a phase separation between a chiral liquid and a dilute gas phase that can be characterized using a thermodynamic framework. We compute the equations of state and show that the surface tension controls interface corrections to the coexisting pressure predicted from the equal-area construction. Transverse forces increase surface tension and generate edge currents at the liquid-gas interface. The analysis of these currents shows that the rotational viscosity introduced in chiral hydrodynamics is consistent with microscopic bulk mechanical measurements. Chirality can also break the solid phase, giving rise to a dense fluid made of rotating hexatic patches. Our Letter paves the way for the development of the statistical mechanics of chiral particles assemblies.

Motility-induced coexistence of a hot liquid and a cold gas

If two phases exist at the same time, such as a gas and a liquid, they have the same temperature. This fundamental law of equilibrium physics is known to apply even to many non-equilibrium systems. However, recently, there has been much attention in the finding that inertial self-propelled particles like Janus colloids in a plasma or microflyers could self-organize into a hot gas-like phase that coexists with a colder liquid-like phase. Here, we show that a kinetic temperature difference across coexisting phases can occur even in equilibrium systems when adding generic (overdamped) self-propelled particles. In particular, we consider mixtures of overdamped active and inertial passive Brownian particles and show that when they phase separate into a dense and a dilute phase, both phases have different kinetic temperatures. Surprisingly, we find that the dense phase (liquid) cannot only be colder but also hotter than the dilute phase (gas). This effect hinges on correlated motions where active particles collectively push and heat up passive ones primarily within the dense phase. Our results answer the fundamental question if a non-equilibrium gas can be colder than a coexisting liquid and create a route to equip matter with self-organized domains of different kinetic temperatures.

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.

Kinetic temperature and pressure of an active Tonks gas

Using computer simulation and analytical theory, we study an active analog of the well-known Tonks gas, where active Brownian particles are confined to a periodic one-dimensional (1D) channel. By introducing the notion of a kinetic temperature, we derive an accurate analytical expression for the pressure and clarify the paradoxical behavior where active Brownian particles confined to 1D exhibit anomalous clustering but no motility-induced phase transition. More generally, this work provides a deeper understanding of pressure in active systems as we uncover a unique link between the kinetic temperature and swim pressure valid for active Brownian particles in higher dimensions.

Effect of speed fluctuations on the collective dynamics of active disks

Numerical simulations are performed on the collective dynamics of active disks, whose self-propulsion speed (U) varies in time, and whose orientation evolves according to rotational Brownian motion. Two protocols for the evolution of speed are considered: (i) a deterministic one involving a periodic change in U at a frequency ω; and (ii) a stochastic one in which the speeds are drawn from a power-law distribution at time-intervals governed by a Poissonian process of rate β. In the first case, an increase in ω causes the disks to go from a clustered state to a homogeneous one through an apparent phase-transition, provided that the direction of self-propulsion is allowed to reverse. Similarly, in the second case, for a fixed value of β, the extent of cluster-breakup is larger when reversals in the self-propulsion direction are permitted. Motility-induced phase separation of the disks may therefore be avoided in active matter suspensions in which the constituents are allowed to reverse their self-propulsion direction, immaterial of the precise temporal nature of the reversal (deterministic or stochastic). Equally, our results demonstrate that phase separation could occur even in the absence of a time-averaged motility of an individual active agent, provided that the rate of direction reversals is smaller than the orientational diffusion rate.

Jamming and flocking in the restricted active Potts model

We study the active Potts model with either site occupancy restriction or on-site repulsion to explore jamming and kinetic arrest in a flocking model. The incorporation of such volume exclusion features leads to a surprisingly rich variety of self-organized spatial patterns. While bands and lanes of moving particles commonly occur without or under weak volume exclusion, strong volume exclusion along with low temperature, high activity, and large particle density facilitates jams due to motility-induced phase separation. Through several phase diagrams, we identify the phase boundaries separating the jammed and free-flowing phases and study the transition between these phases which provide us with both qualitative and quantitative predictions of how jamming might be delayed or dissolved. We further formulate and analyze a hydrodynamic theory for the restricted APM which predicts various features of the microscopic model.

Inertia changes evolution of motility-induced phase separation in active matter across particle activity

The effects of inertia in active matter and motility-induced phase separation (MIPS) have attracted growing interest but still remain poorly studied. We studied MIPS behavior in the Langevin dynamics across a broad range of particle activity and damping rate values with molecular dynamic simulations. Here we show that the MIPS stability region across particle activity values consists of several domains separated by discontinuous or sharp changes in susceptibility of mean kinetic energy. These domain boundaries have fingerprints in the system's kinetic energy fluctuations and characteristics of gas, liquid, and solid subphases, such as the number of particles, densities, or the power of energy release due to activity. The observed domain cascade is most stable at intermediate damping rates but loses its distinctness in the Brownian limit or vanishes along with phase separation at lower damping values.

Classical Nucleation Theory for Active Fluid Phase Separation

Classical Nucleation Theory (CNT), linking rare nucleation events to the free-energy landscape of a growing nucleus, is central to understanding phase-change kinetics in passive fluids. Nucleation in nonequilibrium systems is much harder to describe because there is no free energy, but instead a dynamics-dependent quasipotential that typically must be found numerically. Here we extend CNT to a class of active phase-separating systems governed by a minimal field-theoretic model (Active Model B+). In the small noise and supersaturation limits that CNT assumes, we compute analytically the quasipotential, and hence, nucleation barrier, for liquid-vapor phase separation. Crucial to our results, detailed balance, although broken microscopically by activity, is restored along the instanton trajectory, which in CNT involves the nuclear radius as the sole reaction coordinate.

Collective behavior of soft self-propelled disks with rotational inertia

We investigate collective properties of a large system of soft self-propelled inertial disks with active Langevin dynamics simulation in two dimensions. Rotational inertia of the disks is found to favor motility induced phase separation (MIPS), due to increased effective persistence of the disks. The MIPS phase diagram in the parameter space of rotational inertia and disk softness is reported over a range of values of translation inertia and self-propulsion strength of the disks. Our analytical prediction of the phase boundary between the homogeneous (no-MIPS) and MIPS state in the limit of small and large rotational inertia is found to agree with the numerical data over a large range of translational inertia. Shape of the high density MIPS phase is found to change from circular to rectangular one as the system moves away from the phase boundary. Structural and dynamical properties of the system, measured by several physical quantities, are found to be invariant in the central region of the high density MIPS phase, whereas they are found to vary gradually near the peripheral region of the high density phase. Importantly, the width of the peripheral region near the phase boundary is much larger compared to the narrow peripheral region far away from the phase boundary. Rich dynamics of the disks inside the high density MIPS phase is addressed. Spatial correlation of velocity of the disks is found to increase with rotational inertia and disk hardness. However, temporal correlation of the disks' velocity is found to be a function of rotational inertia, while it is independent of disk softness.

Active Refrigerators Powered by Inertia

We present the operational principle for a refrigerator that uses inertial effects in active Brownian particles to locally reduce their (kinetic) temperature by 2 orders of magnitude below the environmental temperature. This principle exploits the peculiar but so-far unknown shape of the phase diagram of inertial active Brownian particles to initiate motility-induced phase separation in the targeted cooling regime only. Remarkably, active refrigerators operate without requiring isolating walls opening the route toward using them to systematically absorb and trap, e.g., toxic substances from the environment.

Percolation transition in phase-separating active fluid

The motility-induced phase separation exhibited by active particles with repulsive interactions is well known. We show that the interaction softness of active particles destabilizes the highly ordered dense phase, leading to the formation of a porous cluster which spans the system. This soft limit can also be achieved if the particle motility is increased beyond a critical value, at which the system clearly exhibits all the characteristics of a standard percolation transition. We also show that in the athermal limit, active particles exhibit similar transitions even at low motility. With these additional new phases, the phase diagram of repulsive active particles is revealed to be richer than what was previously conceived.

Unjamming and emergent nonreciprocity in active ploughing through a compressible viscoelastic fluid

A dilute suspension of active Brownian particles in a dense compressible viscoelastic fluid, forms a natural setting to study the emergence of nonreciprocity during a dynamical phase transition. At these densities, the transport of active particles is strongly influenced by the passive medium and shows a dynamical jamming transition as a function of activity and medium density. In the process, the compressible medium is actively churned up - for low activity, the active particle gets self-trapped in a cavity of its own making, while for large activity, the active particle ploughs through the medium, either accompanied by a moving anisotropic wake, or leaving a porous trail. A hydrodynamic approach makes it evident that the active particle generates a long-range density wake which breaks fore-aft symmetry, consistent with the simulations. Accounting for the back-reaction of the compressible medium leads to (i) dynamical jamming of the active particle, and (ii) a dynamical non-reciprocal attraction between two active particles moving along the same direction, with the trailing particle catching up with the leading one in finite time. We emphasize that these nonreciprocal effects appear only when the active particles are moving and so manifest in the vicinity of the jamming-unjamming transition.

Critical behavior of active Brownian particles: Connection to field theories

We explore the relation between active Brownian particles, a minimal particle-based model for active matter, and scalar field theories. Both show a liquid-gas-like phase transition toward stable coexistence of a dense liquid with a dilute active gas that terminates in a critical point. However, a comprehensive mapping between the particle-based model parameters and the effective coefficients governing the field theories has not been established yet. We discuss conflicting recent numerical results for the critical exponents of active Brownian particles in two dimensions. Starting from the intermediate effective hydrodynamic equations, we then present a construction for a scalar order parameter for active Brownian particles that yields the active model B+. We argue that a crucial ingredient is the coupling between density and polarization in the particle current. The renormalization flow close to two dimensions exhibits a pair of perturbative fixed points that limit the attractive basin of the Wilson-Fisher fixed point, with the perspective that the critical behavior of active Brownian particles in two dimensions is governed by a strong-coupling fixed point different from Wilson-Fisher and not necessarily corresponding to Ising universality.

Dynamical anomalies and structural features of Active Brownian Particles characterised by two repulsive length scales

In this work we study a two-dimensional system composed by Active Brownian Particles (ABPs) interacting via a repulsive potential with two-length-scales, a soft shell and a hard-core. Depending on the ratio between the strength of the soft shell barrier and the activity, we find two regimes: If this ratio is much larger or smaller than 1, the observed behavior is comparable with ABPs interacting via a single length-scale potential. If this ratio is similar to 1, the two length-scales are relevant for both structure and dynamical properties. On the structural side, when the system exhibits a motility induced phase separation, the dense phase is characterised by new and more complex structures compared with the hexatic phase observed in single length-scale systems.On the dynamical side, as far as we are aware, this is the first representation of an anomalous dynamics in active particles.

A broad perspective to particle-laden fluid interfaces systems: from chemically homogeneous particles to active colloids

Particles adsorbed to fluid interfaces are ubiquitous in industry, nature or life. The wide range of properties arising from the assembly of particles at fluid interface has stimulated an intense research activity on shed light to the most fundamental physico-chemical aspects of these systems. These include the mechanisms driving the equilibration of the interfacial layers, trapping energy, specific inter-particle interactions and the response of the particle-laden interface to mechanical perturbations and flows. The understanding of the physico-chemistry of particle-laden interfaces becomes essential for taking advantage of the particle capacity to stabilize interfaces for the preparation of different dispersed systems (emulsions, foams or colloidosomes) and the fabrication of new reconfigurable interface-dominated devices. This review presents a detailed overview of the physico-chemical aspects that determine the behavior of particles trapped at fluid interfaces. This has been combined with some examples of real and potential applications of these systems in technological and industrial fields. It is expected that this information can provide a general perspective of the topic that can be exploited for researchers and technologist non-specialized in the study of particle-laden interfaces, or for experienced researcher seeking new questions to solve.

Intrinsic structure perspective for MIPS interfaces in two-dimensional systems of active Brownian particles

Suspensions of active Brownian particles (ABPs) undergo motility-induced phase separation (MIPS) over a wide range of mean density and activity strength, which implies the spontaneous aggregation of particles due to the persistence of their direction of motion, even in the absence of an explicit attraction. Both similarities and qualitative differences have been obtained when the MIPS is analysed in the same terms as a liquid-gas phase coexistence in an equilibrium attractive system. Negative values of the mechanical surface tension have been reported, from the total forces across the interface, while stable fluctuations of the interfacial line could be interpreted as a positive capillary surface tension; in equilibrium liquid surfaces, these two magnitudes are equal. We present here the analysis of 2D-ABP interfaces in terms of the intrinsic density and force profiles, calculated with the particle distance to the instantaneous interfacial line. Our results provide new insight into the origin of MIPS from the local rectification of the random active force on the particles near the interface. As has been reported, this effect acts as an external potential that produces a pressure gradient across the interface, such that the mechanical surface tension of the MIPS cannot be described as that of equilibrium coexisting phases; however, our analysis shows that most of that effect comes from the tightly caged particles at the dense (inner) side of the MIPS interface, rather than from the free moving particles at the outer side that collide with the dense cluster. Moreover, a clear correlation appears between the decay of the hexatic order parameter at the dense slab and the end of the MIPS as the strength of the active force is lowered. We show that, using the strong active forces required for MIPS, the interfacial structure and properties are very similar for ABPs with purely repulsive interactions (the Weeks-Chandler-Andersen-Lennard-Jones (WCA-LJ) model truncated at its minimum) and when the interaction includes a range of the LJ attractive force.

Unified analysis of topological defects in 2D systems of active and passive disks

We provide a comprehensive quantitative analysis of localized and extended topological defects in the steady state of 2D passive and active repulsive Brownian disk systems. We show that, both in and out-of-equilibrium, the passage from the solid to the hexatic is driven by the unbinding of dislocations, in quantitative agreement with the KTHNY singularity. Instead, extended clusters of defects largely dominate below the solid-hexatic critical line. The latter percolate in the liquid phase very close to the hexatic-liquid transition, both for continuous and discontinuous transitions, in the homogeneous liquid regime. At critical percolation the clusters of defects are fractal with statistical and geometric properties that are independent of the activity and compatible with the universality class of uncorrelated critical percolation. We also characterize the spatial organization of point-like defects and we show that the disclinations are not free, but rather always very near more complex defect structures. At high activity, the bulk of the dense phase generated by Motility-Induced Phase Separation is characterized by a density of point-like defects, and statistics and morphology of defect clusters, set by the amount of activity and not the packing fraction. Hexatic domains within the dense phase are separated by grain-boundaries along which a finite network of topological defects resides, interrupted by gas bubbles in cavitation. This structure is dynamic in the sense that the defect network allows for an unzipping mechanism that leaves free space for gas bubbles to appear, close, and even be released into the dilute phase.

Phase separation of self-propelled disks with ferromagnetic and nematic alignment

We present a comprehensive study of a model system of repulsive self-propelled disks in two dimensions with ferromagnetic and nematic velocity alignment interactions. We characterize the phase behavior of the system as a function of the alignment and self-propulsion strength, featuring orientational order for strong alignment and motility-induced phase separation (MIPS) at moderate alignment but high enough self-propulsion. We derive a microscopic theory for these systems yielding a closed set of hydrodynamic equations from which we perform a linear stability analysis of the homogenous disordered state. This analysis predicts MIPS in the presence of aligning torques. The nature of the continuum theory allows for an explicit quantitative comparison with particle-based simulations, which consistently shows that ferromagnetic alignment fosters phase separation, while nematic alignment does not alter either the nature or the location of the instability responsible for it. In the ferromagnetic case, such behavior is due to an increase of the imbalance of the number of particle collisions along different orientations, giving rise to the self-trapping of particles along their self-propulsion direction. On the contrary, the anisotropy of the pair correlation function, which encodes this self-trapping effect, is not significantly affected by nematic torques. Our work shows the predictive power of such microscopic theories to describe complex active matter systems with different interaction symmetries and sheds light on the impact of velocity-alignment interactions in motility-induced phase separation.

Wetting Transition of Active Brownian Particles on a Thin Membrane

We study nonequilibrium analogues of surface phase transitions in a minimal model of active particles in contact with a purely repulsive potential barrier that mimics a thin porous membrane. Under conditions of bulk motility-induced phase separation, the interaction strength ϵ_{w} of the barrier controls the affinity of the dense phase for the barrier region. We uncover clear signatures of a wetting phase transition as ϵ_{w} is varied. In common with its equilibrium counterpart, the character of this transition depends on the system dimensionality: a continuous transition with large density fluctuations and gas bubbles is uncovered in 2D while 3D systems exhibit a sharp transition absent of large correlations.

Turbulence-induced clustering in compressible active fluids

We study a novel phase of active polar fluids, which is characterized by the continuous creation and destruction of dense clusters due to self-sustained turbulence. This state arises due to the interplay between self-advection of the aligned swimmers and their defect topology. The typical cluster size is determined by the characteristic vortex size. Our results are obtained by investigating a continuum model of compressible polar active fluids, which incorporates typical experimental observations in bacterial suspensions, in particular a non-monotone dependence of speed on density.

Tracer dynamics in crowded active-particle suspensions

We discuss the dynamics of active Brownian particles (ABPs) in crowded environments through the mean-squared displacement (MSD) of active and passive tracer particles in both active and passive host systems. Exact equations for the MSD are derived using a projection operator technique, extending to dense systems the known solution for a single ABP. The interaction of the tracer particle with the host particles gives rise to strong memory effects. Evaluating these approximately in the framework of a recently developed mode-coupling theory for active Brownian particles (ABP-MCT), we discuss the various dynamical regimes that emerge: While self-propelled motion gives rise to super-diffusive MSD, at high densities, this competes with an interaction-induced sub-diffusive regime. The predictions of the theory are shown to be in good agreement with results obtained from an event-driven Brownian dynamics (ED-BD) simulation scheme for the dynamics of two-dimensional active Brownian hard disks.

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

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

Dynamic overlap concentration scale of active colloids

By introducing the notion of a dynamic overlap concentration scale, we identify additional universal features of the mechanical properties of active colloids. We codify these features by recognizing that the characteristic length scale of an active particle's trajectory, the run length, introduces a concentration scale ϕ^{*}. Large-scale simulations of repulsive active Brownian particles (ABPs) confirm that this run-length dependent concentration, the trajectory-space analog of the overlap concentration in polymer solutions, delineates distinct concentration regimes in which interparticle collisions alter particle trajectories. Using ϕ^{*} and concentration scales associated with colloidal jamming, the mechanical equation of state for ABPs collapses onto a set of principal curves that contain several overlooked features. The inclusion of these features qualitatively alters previous predictions of the behavior for active colloids, as we demonstrate by computing the spinodal for a suspension of purely repulsive ABPs. Our findings suggest that dynamic overlap concentration scales should help unravel the behavior of active and driven systems.

Shear-induced orientational ordering in an active glass former

Dense assemblies of self-propelled particles that can form solid-like states also known as active or living glasses are abundant around us, covering a broad range of length scales and timescales: from the cytoplasm to tissues, from bacterial biofilms to vehicular traffic jams, and from Janus colloids to animal herds. Being structurally disordered as well as strongly out of equilibrium, these systems show fascinating dynamical and mechanical properties. Using extensive molecular dynamics simulation and a number of distinct dynamical and mechanical order parameters, we differentiate three dynamical steady states in a sheared model active glassy system: 1) a disordered state, 2) a propulsion-induced ordered state, and 3) a shear-induced ordered state. We supplement these observations with an analytical theory based on an effective single-particle Fokker-Planck description to rationalize the existence of the shear-induced orientational ordering behavior in an active glassy system without explicit aligning interactions of, for example, Vicsek type. This ordering phenomenon occurs in the large persistence time limit and is made possible only by the applied steady shear. Using a Fokker-Planck description with parameters that can be measured independently, we make testable predictions for the joint distribution of single-particle position and orientation. These predictions match well with the joint distribution measured from direct numerical simulation. Our results are of relevance for experiments exploring the rheological response of dense active colloids and jammed active granular matter systems.

Interaction from structure using machine learning: in and out of equilibrium

Prediction of pair potential given a typical configuration of an interacting classical system is a difficult inverse problem. There exists no exact result that can predict the potential given the structural information. We demonstrate that using machine learning (ML) one can get a quick but accurate answer to the question: "which pair potential lead to the given structure (represented by pair correlation function)?" We use artificial neural network (NN) to address this question and show that this ML technique is capable of providing very accurate prediction of pair potential irrespective of whether the system is in a crystalline, liquid or gas phase. We show that the trained network works well for sample system configurations taken from both equilibrium and out of equilibrium simulations (active matter systems) when the later is mapped to an effective equilibrium system with a modified potential. We show that the ML prediction about the effective interaction for the active system is not only useful to make prediction about the MIPS (motility induced phase separation) phase but also identifies the transition towards this state.

Machine learning for phase behavior in active matter systems

We demonstrate that deep learning techniques can be used to predict motility-induced phase separation (MIPS) in suspensions of active Brownian particles (ABPs) by creating a notion of phase at the particle level. Using a fully connected network in conjunction with a graph neural network we use individual particle features to predict to which phase a particle belongs. From this, we are able to compute the fraction of dilute particles to determine if the system is in the homogeneous dilute, dense, or coexistence region. Our predictions are compared against the MIPS binodal computed from simulation. The strong agreement between the two suggests that machine learning provides an effective way to determine the phase behavior of ABPs and could prove useful for determining more complex phase diagrams.

Phase behavior and surface tension of soft active Brownian particles

We study quasi two-dimensional, monodisperse systems of active Brownian particles (ABPs) for a range of activities, stiffnesses, and densities. We develop a microscopic, analytical method for predicting the dense phase structure formed after motility-induced phase separation (MIPS) has occurred, including the dense cluster's area fraction, interparticle pressure, and radius. Our predictions are in good agreement with our Brownian dynamics simulations. We, then, derive a continuum model to investigate the relationship between the predicted interparticle pressure, the swim pressure, and the macroscopic pressure in the momentum equation. We find that formulating the point-wise macroscopic pressure as the interparticle pressure and modeling the particle activity through a spatially variant body force - as opposed to a volume-averaged swim pressure - results in consistent predictions of pressure in both the continuum model and the microscopic theory. This formulation of pressure also results in nearly zero surface tension for the phase separated domains, irrespective of activity, stiffness, and area fraction. Furthermore, using Brownian dynamics simulations and our continuum model, we showed that both the interface width and surface tension, are intrinsic characteristics of the system. On the other hand, if we were to exclude the body force induced by activity, we find that the resulting surface tension values are linearly dependent on the size of the simulation, in contrast to the statistical mechanical definition of surface tension.

Influence of interaction softness on phase separation of active particles

Using a minimal model of active Brownian particles, we study the effect of a crucial parameter, namely the softness of the interparticle repulsion, on motility-induced phase separation. We show that an increase in particle softness reduces the ability of the system to phase separate and the system exhibits a delayed transition. After phase separation, the system state properties can be explained by a single relevant length scale, the effective interparticle distance. We estimate this length scale analytically and use it to rescale the state properties at dense phase for systems with different interaction softness. Using this length scale, we provide a scaling relation for the time taken to phase separate which shows a high sensitivity to the interaction softness.

Phase Diagram of Active Brownian Spheres: Crystallization and the Metastability of Motility-Induced Phase Separation

Motility-induced phase separation (MIPS), the phenomenon in which purely repulsive active particles undergo a liquid-gas phase separation, is among the simplest and most widely studied examples of a nonequilibrium phase transition. Here, we show that states of MIPS coexistence are in fact only metastable for three-dimensional active Brownian particles over a very broad range of conditions, decaying at long times through an ordering transition we call active crystallization. At an activity just above the MIPS critical point, the liquid-gas binodal is superseded by the crystal-fluid coexistence curve, with solid, liquid, and gas all coexisting at the triple point where the two curves intersect. Nucleating an active crystal from a disordered fluid, however, requires a rare fluctuation that exhibits the nearly close-packed density of the solid phase. The corresponding barrier to crystallization is surmountable on a feasible timescale only at high activity, and only at fluid densities near maximal packing. The glassiness expected for such dense liquids at equilibrium is strongly mitigated by active forces, so that the lifetime of liquid-gas coexistence declines steadily with increasing activity, manifesting in simulations as a facile spontaneous crystallization at extremely high activity.

Characterization of MIPS in a suspension of repulsive active Brownian particles through dynamical features

We study a two-dimensional system composed by Active Brownian Particles (ABPs), focusing on the onset of Motility Induced Phase Separation (MIPS), by means of molecular dynamics simulations. For a pure hard-disk system with no translational diffusion, the phase diagram would be completely determined by their density and Péclet number. In our model, two additional effects are present: translational noise and the overlap of particles; we study the effects of both in the phase space. As we show, the second effect can be mitigated if we use, instead of the standard Weeks-Chandler-Andersen potential, a stiffer potential: the pseudo-hard sphere potential. Moreover, in determining the boundary of our phase space, we explore different approaches to detect MIPS and conclude that observing dynamical features, via the non-Gaussian parameter, is more efficient than observing structural ones, such as through the local density distribution function. We also demonstrate that the Vogel-Fulcher equation successfully reproduces the decay of the diffusion as a function of density, with the exception of very high densities. Thus, in this regard, the ABP system behaves similar to a fragile glass.

How to study a persistent active glassy system

We explore glassy dynamics of dense assemblies of soft particles that are self-propelled by active forces. These forces have a fixed amplitude and a propulsion direction that varies on a timescaleτ_p, the persistence timescale. Numerical simulations of such active glasses are computationally challenging when the dynamics is governed by large persistence times. We describe in detail a recently proposed scheme that allows one to study directly the dynamics in the large persistence time limit, on timescales around and well above the persistence time. We discuss the idea behind the proposed scheme, which we call 'activity-driven dynamics', as well as its numerical implementation. We establish that our prescription faithfully reproduces all dynamical quantities in the appropriate limitτ_p→ ∞. We deploy the approach to explore in detail the statistics of Eshelby-like plastic events in the steady state dynamics of a dense and intermittent active glass.

Universality class of the motility-induced critical point in large scale off-lattice simulations of active particles

We perform large-scale computer simulations of an off-lattice two-dimensional model of active particles undergoing a motility-induced phase separation (MIPS) to investigate the system's critical behaviour close to the critical point of the MIPS curve. By sampling steady-state configurations for large system sizes and performing finite size scaling analysis we provide exhaustive evidence that the critical behaviour of this active system belongs to the Ising universality class. In addition to the scaling observables that are also typical of passive systems, we study the critical behaviour of the kinetic temperature difference between the two active phases. This quantity, which is always zero in equilibrium, displays instead a critical behavior in the active system which is well described by the same exponent of the order parameter in agreement with mean-field theory.

Phase Separation and Multibody Effects in Three-Dimensional Active Brownian Particles

Simulation studies of the phase diagram of repulsive active Brownian particles in three dimensions reveal that the region of motility-induced phase separation between a high and low density phase is enclosed by a region of gas-crystal phase separation. Near-critical loci and structural crossovers can additionally be identified in analogy with simple fluids. Motivated by the striking similarity to the behavior of equilibrium fluids with short-ranged pairwise attractions, we show that a direct mapping to pair potentials in the dilute limit implies interactions that are insufficiently attractive to engender phase separation. Instead, this is driven by the emergence of multibody effects associated with particle caging that occurs at sufficiently high number density. We quantify these effects via information-theoretical measures of n-body effective interactions extracted from the configurational structure.

Motility-Induced Microphase and Macrophase Separation in a Two-Dimensional Active Brownian Particle System

As a result of nonequilibrium forces, purely repulsive self-propelled particles undergo macrophase separation between a dense and a dilute phase. We present a thorough study of the ordering kinetics of such motility-induced phase separation (MIPS) in active Brownian particles in two dimensions, and we show that it is generically accompanied by microphase separation. The growth of the dense phase follows a law akin to the one of liquid-gas phase separation. However, it is made of a mosaic of hexatic microdomains whose size does not coarsen indefinitely, leaving behind a network of extended topological defects from which microscopic dilute bubbles arise. The characteristic length of these finite-size structures increases with activity, independently of the choice of initial conditions.

Frictional active Brownian particles

Frictional forces affect the rheology of hard-sphere colloids, at high shear rate. Here we demonstrate, via numerical simulations, that they also affect the dynamics of active Brownian particles and their motility-induced phase separation. Frictional forces increase the angular diffusivity of the particles, in the dilute phase, and prevent colliding particles from resolving their collision by sliding one past to the other. This leads to qualitatively changes of motility-induced phase diagram in the volume-fraction motility plane. While frictionless systems become unstable towards phase separation as the motility increases only if their volume fraction overcomes a threshold, frictional systems become unstable regardless of their volume fraction. These results suggest the possibility of controlling the motility-induced phase diagram by tuning the roughness of the particles.

Rod-assisted heterogeneous nucleation in active suspensions

Motility induced phase separation as well as the nucleation process in active particle systems has gained extensive research attention very recently. Most studies so far have considered homogeneous cases without the influence of foreign seeds or impurities; however, the heterogeneous nucleation process, widely studied in passive systems, has not been systematically investigated yet. Here we study the heterogeneous nucleation process and phase behaviors of a suspension of active Brownian particles by introducing a rod-like passive seed. We found that such a seed can exponentially accelerate the nucleation rate and thus readily induce phase separation of a dilute active system, while a homogeneous one with the same volume fraction still maintains a single phase. It is observed that the seed would automatically detach from the dense phase after the completion of phase separation instead of staying inside as an impurity. Interestingly, we found that the phase behavior is re-entrant with the activity: single-phase states exist at both high and low activities, with phase separated states in between. Our results demonstrate that heterogeneous nucleation in an active system can show novel behaviors with respect to its passive counterpart, and pave the way for more future studies in relevant fields.

Different-shaped micro-objects driven by active particle aggregations

The dynamics of passive micro-objects in an active bath has been receiving much attention. However, the influence of the shapes of micro-objects remains unclear. Here, we use 2D simulation to investigate the interaction between active Brownian particles and different-shaped passive micro-objects. We show that active particles accumulate around micro-objects and self-assemble into living aggregations at a high active velocity and high volume fraction. The shapes of micro-objects affect the distributions of the aggregations. In turn, the different distribution of aggregations influences the motion of micro-objects and induces abnormal diffusive behaviors. We further demonstrate that polar distributed aggregations at a high active velocity and the inhibition of the active bath at a low active velocity induce the counterintuitive anisotropic enhanced diffusion of rods, and the steric interaction between active particles induces the reverse translation-rotation coupled diffusion of chevrons.

Clustering and phase separation in mixtures of dipolar and active particles

The self-assembly of colloidal particles in dynamic environments has become an important field of study because of potential applications in fabricating out-of-equilibrium materials. We investigate the phase behavior of mixtures of passive dipolar colloids and active soft spheres using Brownian dynamics simulations in two dimensions. The phase behaviors exhibited include dipolar percolated network, dipolar string-fluid, isotropic fluid, and a phase-separated state. We find that the clustering of dipolar colloids is enhanced in the presence of slow-moving active particles compared to the clustering of dipolar particles mixed with passive particles. When the active particle motility is high, the chains of dipolar particles are either broken into short chains or pushed into dense clusters. Motility-induced phase separation into dense and dilute phases is also present. The area fraction of particles in the dilute phase increases as the fraction of active particles in the system decreases, while the area fraction of particles in the dense phase remains constant. Our findings are relevant to the development of reconfigurable self-assembled materials.

Collective forces in scalar active matter

Large-scale collective behavior in suspensions of active particles can be understood from the balance of statistical forces emerging beyond the direct microscopic particle interactions. Here we review some aspects of the collective forces that can arise in suspensions of self-propelled active Brownian particles: wall forces under confinement, interfacial forces, and forces on immersed bodies mediated by the suspension. Even for non-aligning active particles, these forces are intimately related to a non-uniform polarization of particle orientations induced by walls and bodies, or inhomogeneous density profiles. We conclude by pointing out future directions and promising areas for the application of collective forces in synthetic active matter, as well as their role in living active matter.

Quorum-sensing active particles with discontinuous motility

We develop a dynamic mean-field theory for polar active particles that interact through a self-generated field, in particular one generated through emitting a chemical signal. While being a form of chemotactic response, it is different from conventional chemotaxis in that particles discontinuously change their motility when the local concentration surpasses a threshold. The resulting coupled equations for density and polarization are linear and can be solved analytically for simple geometries, yielding inhomogeneous density profiles. Specifically, here we consider a planar and circular interface. Our theory thus explains the observed coexistence of dense aggregates with an active gas. There are, however, differences from the more conventional picture of liquid-gas coexistence based on a free energy, most notably the absence of a critical point. We corroborate our analytical predictions by numerical simulations of active particles under confinement and interacting through volume exclusion. Excellent quantitative agreement is reached through an effective translational diffusion coefficient. We finally show that an additional response to the chemical gradient direction is sufficient to induce vortex clusters. Our results pave the way to engineer motility responses in order to achieve aggregation and collective behavior even at unfavorable conditions.

Active Hard Spheres in Infinitely Many Dimensions

Few equilibrium-even less so nonequilibrium-statistical-mechanical models with continuous degrees of freedom can be solved exactly. Classical hard spheres in infinitely many space dimensions are a notable exception. We show that, even without resorting to a Boltzmann distribution, dimensionality is a powerful organizing device for exploring the stationary properties of active hard spheres evolving far from equilibrium. In infinite dimensions, we exactly compute the stationary state properties that govern and characterize the collective behavior of active hard spheres: the structure factor and the equation of state for the pressure. In turn, this allows us to account for motility-induced phase separation. Finally, we determine the crowding density at which the effective propulsion of a particle vanishes.

Linear Response Theory and Green-Kubo Relations for Active Matter

We address the question of how interacting active systems in a nonequilibrium steady state respond to an external perturbation. We establish an extended fluctuation-dissipation theorem for active Brownian particles (ABP), which highlights the role played by the local violation of detailed balance due to activity. By making use of a Markovian approximation we derive closed Green-Kubo expressions for the diffusivity and mobility of ABP and quantify the deviations from the Stokes-Einstein relation. We compute the linear response function to an external force using unperturbed simulations of ABP and compare the results with the analytical predictions of the transport coefficients. Our results show the importance of the interplay between activity and interactions in the departure from equilibrium linear response.

Motility-Induced Temperature Difference in Coexisting Phases

Unlike in thermodynamic equilibrium where coexisting phases always have the same temperature, here we show that systems comprising "active" self-propelled particles can self-organize into two coexisting phases at different kinetic temperatures, which are separated from each other by a sharp and persistent temperature gradient. Contrasting previous studies that have focused on overdamped descriptions of active particles, we show that a "hot-cold coexistence" occurs if and only if accounting for inertia, which is significant, e.g., in activated dusty plasmas, microflyers, whirling fruits, or beetles at interfaces. Our results exemplify a route to use active particles to create a self-sustained temperature gradient across coexisting phases. This phenomenon is fundamentally beyond equilibrium physics and is accompanied by a slow coarsening law with an exponent significantly smaller than the universal 1/3 exponent seen in both equilibrium systems and overdamped active Brownian particles.

Which interactions dominate in active colloids?

Despite mounting evidence that the same gradients, which active colloids use for swimming, induce important cross-interactions (phoretic interactions), they are still ignored in most many-body descriptions, perhaps to avoid complexity and a zoo of unknown parameters. Here we derive a simple model, which reduces phoretic far-field interactions to a pair-interaction whose strength is mainly controlled by one genuine parameter (swimming speed). The model suggests that phoretic interactions are generically important for autophoretic colloids (unless effective screening of the phoretic fields is strong) and should dominate over hydrodynamic interactions for the typical case of half-coating and moderately nonuniform surface mobilities. Unlike standard minimal models, but in accordance with canonical experiments, our model generically predicts dynamic clustering in active colloids at a low density. This suggests that dynamic clustering can emerge from the interplay of screened phoretic attractions and active diffusion.

Synthetic Chemotaxis and Collective Behavior in Active Matter

The ability to navigate in chemical gradients, called chemotaxis, is crucial for the survival of microorganisms. It allows them to find food and to escape from toxins. Many microorganisms can produce the chemicals to which they respond themselves and use chemotaxis for signaling, which can be seen as a basic form of communication, allowing ensembles of microorganisms to coordinate their behavior, for example, during embryogenesis, biofilm formation, or cellular aggregation. For example, Dictyostelium cells use signaling as a survival strategy: when starving, they produce certain chemicals toward which other cells show taxis. This leads to aggregation of the cells resulting in a multicellular aggregate that can sustain long starvation periods. Remarkably, the past decade has led to the development of synthetic microswimmers, which can self-propel through a solvent, analogously to bacteria and other microorganisms. The mechanism underlying the self-propulsion of synthetic microswimmers like camphor boats, droplet swimmers, and in particular autophoretic Janus colloids involves the production of certain chemicals. As we will discuss in this Account, the same chemicals (phoretic fields) involved in the self-propulsion of a (Janus) microswimmer also act on other ones and bias their swimming direction toward (or away from) the producing microswimmer. Synthetic microswimmers therefore provide a synthetic analogue to motile microorganisms interacting by taxis toward (or away from) self-produced chemical fields. In this Account, we review recent progress in the theoretical description of synthetic chemotaxis mainly based on simulations and field theoretical descriptions. We will begin with single motile particles leaving chemical trails behind with which they interact themselves, leading to effects like self-trapping or self-avoidance. Besides these self-interactions, in ensembles of synthetic motile particles each particle also responds to the chemicals produced by other particles, inducing chemical (or phoretic) cross-interactions. When these interactions are attractive, they commonly lead to clusters, even at low particle density. These clusters may either proceed toward macrophase separation, resembling Dictyostelium aggregation, or, as shown very recently, lead to dynamic clusters of self-limited size (dynamic clustering) as seen in experiments in autophoretic Janus colloids. Besides the classical case where chemical interactions are attractive, this Account discusses, as its main focus, repulsive chemical interactions, which can create a new and less known avenue to pattern formation in active systems leading to a variety of pattern, including clusters which are surrounded by shells of chemicals, traveling waves and more complex continuously reshaping patterns. In all these cases "synthetic signalling" can crucially determine the collective behavior of synthetic microswimmer ensembles and can be used as a design principle to create patterns in motile active particles.