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

Phase separation in active binary mixtures with chemical reaction

We study motility-induced phase separation (MIPS) in active AB binary mixtures undergoing the chemical reaction A ⇌ B. Starting from the evolution equations for the density fields ρi(r⃑,t) describing MIPS, we phenomenologically incorporate the effects of the reaction through the reaction rate Γ into the equations. The steady-state domain morphologies depend on Γ and the relative activity of the species, Δ. For a sufficiently large Γ and Δ ≠ 1, the more active component of the mixture forms a droplet morphology. We characterize the morphology of domains by calculating the equal-time correlation function C(r,t) and the structure factor S(k,t), exhibiting scaling violation. The average domain size, L(t), follows a diffusive growth as L(t) ∼ t1/3 before reaching the steady state domain size, Lss. Additionally, Lss shows the scaling relation Lss∼ Γ-1/4, independent of Δ.

Kinetically arrested periodic clusters in active filament arrays

We study the dynamics and pattern formation of ordered arrays of active semi-flexible filaments, each of which is pinned at one end and free at the other. The filaments are modeled as connected chains of polar active particles with activity incorporated through local follower forces acting along their local tangent. Using Brownian dynamics simulations in two dimensions, we show that for a range of activity and filament separation, the filament array self-assembles into regularly spaced, kinetically arrested compact clusters. Activity, array geometry, filament elasticity, and grafting density are each seen to crucially influence the show size, shape, and spacing of these emergent clusters. Furthermore, cluster shapes for different grafting densities can be rescaled into self-similar forms with activity-dependent scaling exponents. We derive theoretical expressions that relate the number of filaments in a cluster and the spacing between adjacent clusters to filament activity, filament elasticity, and grafting density. Our results provide insight into the physical mechanisms involved in the initiation of clustering and suggest that steric contact forces and friction balance active forces and filament elasticity to shape and stabilize emergence clusters.

Spontaneous generation of angular momentum in chiral active crystals

We study a two-dimensional chiral active crystal composed of underdamped chiral active particles. These particles, characterized by intrinsic handedness and persistence, interact via linear forces derived from harmonic potentials. Chirality plays a pivotal role in shaping the system's behavior: it reduces displacement and velocity fluctuations while inducing cross-spatial correlations among different Cartesian components of velocity. These features distinguish chiral crystals from their non-chiral counterparts, leading to the emergence of net angular momentum, as predicted analytically. This angular momentum, driven by the torque generated by the chiral active force, exhibits a non-monotonic dependence on the degree of chirality. Additionally, it contributes to the entropy production rate, as revealed through a path-integral analysis. We investigate the dynamic properties of the crystal in both Fourier and real space. Chirality induces a non-dispersive peak in the displacement spectrum, which underlies the generation of angular momentum and oscillations in time-dependent autocorrelation functions or mean-square displacement, all of which are analytically predicted.

Non-equilibrium coexistence between a fluid and a hotter or colder crystal of granular hard disks

Non-equilibrium phase coexistence is commonly observed in both biological and artificial systems, yet understanding it remains a significant challenge. Unlike equilibrium systems, where free energy provides a unifying framework, the absence of such a quantity in non-equilibrium settings complicates their theoretical understanding. Granular materials, driven out of equilibrium by energy dissipation during collisions, serve as an ideal platform to investigate these systems, offering insights into the parallels and distinctions between equilibrium and non-equilibrium phase behavior. For example, the coexisting dense phase is typically colder than the dilute phase, a result usually attributed to greater dissipation in denser regions. In this article, we demonstrate that this is not always the case. Using a simple numerical granular model, we show that a hot solid and a cold liquid can coexist in granular systems. This counterintuitive phenomenon arises because the collision frequency can be lower in the solid phase than in the liquid phase, consistent with equilibrium results for hard-disk systems. We further demonstrate that kinetic theory can be extended to accurately predict phase temperatures even at very high packing fractions, including within the solid phase. Our results highlight the importance of collisional dynamics and energy exchange in determining phase behavior in granular materials, offering new insights into non-equilibrium phase coexistence and the complex physics underlying granular systems.

Interfacial and density fluctuations in a lattice model of motility-induced phase separation

We analyze motility-induced phase separation and bubbly phase separation in a two-dimensional lattice model of self-propelled particles. We compare systems where the dense (liquid) phase has slab and droplet geometries. We find that interfacial fluctuations of the slab are well-described by capillary wave theory, despite the existence of bubbles in the dense phase. We attribute this to a separation of time scales between bubble expulsion and interfacial relaxation. We also characterize the dependence of liquid and vapor densities on the curvature of the liquid droplet, as well as the density fluctuations inside the phases. The vapor phase behaves similarly to an equilibrium system, displaying a Laplace pressure effect that shifts its density, and Gaussian density fluctuations. The liquid phase has large non-Gaussian fluctuations, but this is not accompanied by a large density shift, contrary to the equilibrium case. Nevertheless, the shift of the vapor density can be used to infer an effective surface tension that appears to also quantify capillary wave fluctuations.

Motility-induced collapse of an active Brownian particle polymer chain

The ability of particles to transform absorbed energy into translational movements brings peculiar order into nonequilibrium matter. Connected together into a chain, these particles collectively behave completely differently from well-known equilibrium polymers. Examples of such systems vary from nanoscale to macroscopic objects. Herein, we demonstrate the ability of self-propelled monomer units to cause chain collapse under conditions where there is no explicit attraction between particles. The resulting conformation is heterogeneous, characterized by a dense condensed core surrounded by elongated loops and stabilized by an effective temperature gradient. Polymer segments in the condensed phase display crumpled packing. Using molecular dynamics analysis, we show that this effect is similar to motility-induced phase separation in active media.

Ratchet effect and jamming in dense mixtures of active and passive colloids in narrow pores

Using the framework of generalized exclusion processes, we study mixtures of passive and active particles interacting by steric repulsion. The particles move in a pore with a periodically modulated aperture, which is modeled by a quasi-one-dimensional channel with a periodic tooth-shaped profile. Internal driving of the active particles induces a ratchet current of these particles. In the current-density diagram, we observe three main regimes: of free flow, of thermally activated processes, and of spinodal decomposition. When the density of particles is increased, we observe a transition to a jammed state, where the ratchet current is substantially reduced. In time evolution, the transition to a jammed state is seen as a sudden drop of current at a certain time. The probability distribution of these jamming times follows an exponential law. The average jamming time depends exponentially on the density of active particles. The coefficient in this exponential is nearly independent of the switching rate of the active particles as well as the presence or absence of passive particles. Due to the interaction, the current of active particles imposes a drag on the passive particles. In the limit of both large systems and long times, the current of passive particles always has the same direction as the ratchet current of active particles. However, during the evolution of the system, we observe a very slow (logarithmic in time) approach to the asymptotic value, sometimes accompanied by current reversal, i.e., the current of active and passive particles may go in opposite directions.

The influence of active agent motility on SIRS epidemiological dynamics

Active Brownian disks moving in two dimensions that exchange information about their internal state stochastically are chosen to model epidemic spread in a self-propelled population of agents under the susceptible-infected-recovered-susceptible (SIRS) framework. The state of infection of an agent, or disk, governs its self-propulsion speed; consequently, the activity of the agents in the system varies in time. Two different protocols (one-to-one and one-to-many) are considered for the transmission of disease from the infected to susceptible populations. The effectiveness of the two protocols are practically identical at high values of the infection transmission rate. The one-to-many protocol, however, outperforms the one-to-one protocol at lower values of the infection transmission rate. Salient features of the macroscopic SIRS model are revisited, and compared to predictions from the agent-based model. Lastly, the motility induced phase separation in a population of such agents with a fluctuating fraction of active disks is found to be well-described by theories governing phase separation in a mixture of active and passive particles with a constant fraction of passive disks.

Percolation transitions in a binary mixture of active Brownian particles with different softness

Homogeneous active Brownian particle (ABP) systems with purely repulsive interactions are considered to exhibit a simple phase behavior, but various physical attributes of active entities can lead to variation in the collective dynamics. Recent studies have shown that even homogeneous ABPs exhibit complex behavior due to an interplay between particle softness and motility. However, the heterogeneity in the composition of ABPs has not been explored yet. In this paper, we study the structural properties of a binary mixture of ABPs with different particle softness by varying the relative softness and composition. We found that upon varying the motility parameter, the system underwent a motility-induced phase separation (MIPS) followed by a percolation transition similar to the homogeneous systems. However, we observed a novel feature: the formation of a space-filling structure made of particles with higher stiffness, within the dense cluster of MIPS containing both types of particles. Our systematic analysis shows that this structure formation occurs only if the difference in softness of both types of particles is sufficiently large. Furthermore, the presence of a non-linear scaling for different compositions of binary ABPs suggests that there is a complex relationship between the composition and the structural properties. Our study demonstrates that the composition heterogeneity of ABPs can lead to complex phase behavior.

Aggregate morphing of self-aligning soft active disks in semi-confined geometry

We study the dependence of alignment and confinement on the aggregate morphology of self-aligning soft disks(particles) in a planer box (two dimensional) geometry confined along y direction using Langevin dynamics simulations. We show that when the box width decreases, the aggregate wall accumulation becomes non-uniform and displays non-monotonic behaviour in terms of phase behavior and height of these aggregates with an increase in alignment strength. Additionally, we identify two distinct categories of wall aggregates: layered and non-layered structures each exhibiting distinct local structural properties. For non-layered structures, local speed of the particles stay nearly constant as we move away from the boundary, while for layered structures, it increases with distance from the boundary. Our analysis shows that active pressure difference is a useful indicator for different aggregate morphologies and the peaks in the pressure curve are indicative of the average and minimum height of the structure.

Statistical properties of microphase and bubbly phase-separated active fluids

In phase-separated active fluids, the Ostwald process can go into reverse, leading to either microphase separation or bubbly phase separation. We show that the latter is formed of two macroscopic regions that are occupied by the homogeneous fluid and by the microphase separated one. Within the microphase-separated fluid, the relative rate of the Ostwald process, coalescence, and nucleation determines whether the size distribution of mesoscopic domains is narrowly peaked or displays a broad range of sizes before attaining a cutoff independent of system size. Our results are obtained via large-scale simulations of a minimal field theory for active phase separation and reproduced by an effective model in which the degrees of freedom are the locations and sizes of the microphase-separated domains.

Active string fluids and gels formed by dipolar active Brownian particles in 3D

Self-propelled particles possessing permanent magnetic dipole moments occur naturally in magnetotactic bacteria and can be built into man-made systems such as active colloids or micro-robots. Yet, the interplay between self-propulsion and anisotropic dipole-dipole interactions on dynamic self-assembly in three dimensions (3D) remains poorly understood. We conduct Brownian dynamics simulations of active dipolar particles in 3D, focusing on the low-density regime, where dipolar hard spheres tend to form chain-like aggregates and percolated networks with increasing dipolar coupling strength. We find that strong active forces override dipolar attractions, effectively inhibiting chain-like aggregation and network formation. Conversely, activating particles with low to moderate forces results in a fluid composed of active chains and rings. At strong dipolar coupling strengths, this active fluid transitions into an active gel, consisting of a percolated network of active chains. Although the overall structure of the active gel remains interconnected, the network experiences more frequent configurational rearrangements due to the reduced bond lifetime of active dipolar particles. Consequently, particles exhibit enhanced translational and rotational diffusion within the active fluid of strings and active gels compared to their passive counterparts. We quantify the influence of activity on aggregate topology as they transition from branched structures to unconnected chains and rings. Our findings are summarized in a state diagram, delineating the impact of dipolar coupling strength and active force magnitude on the system.

Noise-Induced Phase Separation and Time Reversal Symmetry Breaking in Active Field Theories Driven by Persistent Noise

Within the Landau-Ginzburg picture of phase transitions, scalar field theories develop phase separation because of a spontaneous symmetry-breaking mechanism. This picture works in thermodynamics but also in the dynamics of phase separation. Here we show that scalar nonequilibrium field theories undergo phase separation just because of nonequilibrium fluctuations driven by a persistent noise. The mechanism is similar to what happens in motility-induced phase separation where persistent motion introduces an effective attractive force. We observe that noise-induced phase separation occurs in a region of the phase diagram where disordered field configurations would otherwise be stable at equilibrium. Measuring the local entropy production rate to quantify the time-reversal symmetry breaking, we find that such breaking is concentrated on the boundary between the two phases.

Phase separation kinetics and cluster dynamics in two-dimensional active dumbbell systems

Molecular dynamics simulations were employed to investigate the phase separation process of a two-dimensional active Brownian dumbbell model. We evaluated the time dependence of the typical size of the dense component using the scaling properties of the structure factor, along with the averaged number of clusters and their radii of gyration. The growth observed is faster than in active disk models, and this effect is further enhanced under stronger activity. Next, we focused on studying the hexatic order of the clusters. The length associated with the orientational order increases algebraically with time and faster than for spherical active Brownian particles. Under weak active forces, most clusters exhibit a uniform internal orientational order. However, under strong forces, large clusters consist of domains with different orientational orders. We demonstrated that the latter configurations are not stable, and given sufficient time to evolve, they eventually achieve homogeneous configurations as well. No gas bubbles are formed within the clusters, even when there are patches of different hexatic order. Finally, attention was directed towards the geometry and motion of the clusters themselves. By employing a tracking algorithm, we showed that clusters smaller than the typical size at the observation time exhibit regular shapes, while larger ones display fractal characteristics. In between collisions or break-ups, the clusters behave as solid bodies. Their centers of mass undergo circular motion, with radii increasing with the cluster size. The angular velocity of the center of mass equals that of the constituents with respect to their center of mass. These observations were rationalised with a simple mechanical model.

Fluctuation Theorems for Heat Exchanges between Passive and Active Baths

In addition to providing general constraints on probability distributions, fluctuation theorems allow us to infer essential information on the role played by temperature in heat exchange phenomena. In this numerical study, we measure the temperature of an out-of-equilibrium active bath using a fluctuation theorem that relates the fluctuations in the heat exchanged between two baths to their temperatures. Our setup consists of a single particle moving between two wells of a quartic potential accommodating two different baths. The heat exchanged between the two baths is monitored according to two definitions: as the kinetic energy carried by the particle whenever it jumps from one well to the other and as the work performed by the particle on one of the two baths when immersed in it. First, we consider two equilibrium baths at two different temperatures and verify that a fluctuation theorem featuring the baths temperatures holds for both heat definitions. Then, we introduce an additional Gaussian coloured noise in one of the baths, so as to make it effectively an active (out-of-equilibrium) bath. We find that a fluctuation theorem is still satisfied with both heat definitions. Interestingly, in this case the temperature obtained through the fluctuation theorem for the active bath corresponds to the kinetic temperature when considering the first heat definition, while it is larger with the second one. We interpret these results by looking at the particle jump phenomenology.

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.

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.

When an active bath behaves as an equilibrium one

Active scalar baths consisting of active Brownian particles are characterized by a non-Gaussian velocity distribution, a kinetic temperature, and a diffusion coefficient that scale with the square of the active velocity v_{0}. While these results hold in overdamped active systems, inertial effects lead to normal velocity distributions, with kinetic temperature and diffusion coefficient increasing as ∼v_{0}^{α} with 1<α<2. Remarkably, the late-time diffusivity and mobility decrease with mass. Moreover, we show that the equilibrium Einstein relation is asymptotically recovered with inertia. In summary, the inertial mass restores an equilibriumlike behavior.

Phase behaviour and dynamics of three-dimensional active dumbbell systems

We present a comprehensive numerical study of the phase behavior and dynamics of a three-dimensional active dumbbell system with attractive interactions. We demonstrate that attraction is essential for the system to exhibit nontrivial phases. We construct a detailed phase diagram by exploring the effects of the system's activity, density, and attraction strength. We identify several distinct phases, including a disordered, a gel, and a completely phase-separated phase. Additionally, we discover a novel dynamical phase, that we name percolating network, which is characterized by the presence of a spanning network of connected dumbbells. In the phase-separated phase we characterize numerically and describe analytically the helical motion of the dense cluster.

Active-parameter polydispersity in the 2d ABP Yukawa model

In experiments and simulations of passive as well as active matter the most commonly studied kind of parameter polydispersity is that of varying particles size. This paper investigates by simulations the effects of introducing polydispersity in other parameters for two-dimensional active Brownian particles with Yukawa pair interactions. Polydispersity is studied separately in the translational and rotational diffusion coefficients, as well as in the swim velocityv0. Uniform and binary parameter distributions are considered in the homogeneous and the motility-induced phase-separation (MIPS) phases. We find only minute changes in structure and dynamics upon the introduction of parameter polydispersity, even for situations involving 50% polydispersity. The reason for this is not clear. An exception is the case ofv0polydispersity for which the average radial distribution function with changing polydispersity shows significant variations in the MIPS phase. Even in this case, however, the dynamics is only modestly affected. As a possible application of our findings, we suggest that a temporary introduction of polydispersity into a single-component active-matter model characterized by a very long equilibration time, i.e. a glass-forming active system, may be used to equilibrate the system efficiently by particle swaps.

Characterizing different motility-induced regimes in active matter with machine learning and noise

We examine motility-induced phase separation (MIPS) in two-dimensional run-and-tumble disk systems using both machine learning and noise fluctuation analysis. Our measures suggest that within the MIPS state there are several distinct regimes as a function of density and run time, so that systems with MIPS transitions exhibit an active fluid, an active crystal, and a critical regime. The different regimes can be detected by combining an order parameter extracted from principal component analysis with a cluster stability measurement. The principal component-derived order parameter is maximized in the critical regime, remains low in the active fluid, and has an intermediate value in the active crystal regime. We demonstrate that machine learning can better capture dynamical properties of the MIPS regimes compared to more standard structural measures such as the maximum cluster size. The different regimes can also be characterized via changes in the noise power of the fluctuations in the average speed. In the critical regime, the noise power passes through a maximum and has a broad spectrum with a 1/f^{1.6} signature, similar to the noise observed near depinning transitions or for solids undergoing plastic deformation.

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.

Dynamical Pattern Formation without Self-Attraction in Quorum-Sensing Active Matter: The Interplay between Nonreciprocity and Motility

We study a minimal model involving two species of particles interacting via quorum-sensing rules. Combining simulations of the microscopic model and linear stability analysis of the associated coarse-grained field theory, we identify a mechanism for dynamical pattern formation that does not rely on the standard route of intraspecies effective attractive interactions. Instead, our results reveal a highly dynamical phase of chasing bands induced only by the combined effects of self-propulsion and nonreciprocity in the interspecies couplings. Turning on self-attraction, we find that the system may phase separate into a macroscopic domain of such chaotic chasing bands coexisting with a dilute gas. We show that the chaotic dynamics of bands at the interfaces of this phase-separated phase results in anomalously slow coarsening.

Entropy production and collective excitations of crystals out of equilibrium: The concept of entropons

We study the collective vibrational excitations of crystals under out-of-equilibrium steady conditions that give rise to entropy production. Their excitation spectrum comprises equilibriumlike phonons of thermal origin and additional collective excitations called entropons because each of them represents a mode of spectral entropy production. Entropons coexist with phonons and dominate them when the system is far from equilibrium while they are negligible in near-equilibrium regimes. The concept of entropons has been recently introduced and verified in a special case of crystals formed by self-propelled particles. Here we show that entropons exist in a broader class of active crystals that are intrinsically out of equilibrium and characterized by the lack of detailed balance. After a general derivation, several explicit examples are discussed, including crystals consisting of particles with alignment interactions and frictional contact forces.

Control of Active Brownian Particles: An Exact Solution

Control of stochastic systems is a challenging open problem in statistical physics, with a wealth of potential applications from biology to granulates. Unlike most cases investigated so far, we aim here at controlling a genuinely out-of-equilibrium system, the two dimensional active Brownian particles model in a harmonic potential, a paradigm for the study of self-propelled bacteria. We search for protocols for the driving parameters (stiffness of the potential and activity of the particles) bringing the system from an initial passivelike steady state to a final activelike one, within a chosen time interval. The exact analytical results found for this prototypical model of self-propelled particles brings control techniques to a wider class of out-of-equilibrium systems.

Chiral active matter in external potentials

We investigate the interplay between chirality and confinement induced by the presence of an external potential. For potentials having radial symmetry, the circular character of the trajectories induced by the chiral motion reduces the spatial fluctuations of the particle, thus providing an extra effective confining mechanism, that can be interpreted as a lowering of the effective temperature. In the case of non-radial potentials, for instance, with an elliptic shape, chirality displays a richer scenario. Indeed, the chirality can break the parity symmetry of the potential that is always fulfilled in the non-chiral system. The probability distribution displays a strong non-Maxwell-Boltzmann shape that emerges in cross-correlations between the two Cartesian components of the position, that vanishes in the absence of chirality or when radial symmetry of the potential is restored. These results are obtained by considering two popular models in active matter, i.e. chiral Active Brownian particles and chiral active Ornstein-Uhlenbeck particles.

Dynamics of Motility-Induced Clusters: Coarsening beyond Ostwald Ripening

We study the dynamics of clusters of active Brownian disks generated by motility-induced phase separation, by applying an algorithm that we devised to track cluster trajectories. We identify an aggregation mechanism that goes beyond Ostwald ripening but also yields a dynamic exponent characterizing the cluster growth z=3, in the timescales explored numerically. Clusters of mass M self-propel with enhanced diffusivity D∼Pe^{2}/sqrt[M]. Their fast motion drives aggregation into large fractal structures, which are patchworks of diverse hexatic orders, and coexist with regular, orientationally uniform, smaller ones. To bring out the impact of activity, we perform a comparative study of a passive system that evidences major differences with the active case.

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.

Collective Flows Drive Cavitation in Spinner Monolayers

Hydrodynamic interactions can give rise to a collective motion of rotating particles. This, in turn, can lead to coherent fluid flows. Using large scale hydrodynamic simulations, we study the coupling between these two in spinner monolayers at weakly inertial regime. We observe an instability, where the initially uniform particle layer separates into particle void and particle rich areas. The particle void region corresponds to a fluid vortex, and it is driven by a surrounding spinner edge current. We show that the instability originates from a hydrodynamic lift force between the particle and fluid flows. The cavitation can be tuned by the strength of the collective flows. It is suppressed when the spinners are confined by a no-slip surface, and multiple cavity and oscillating cavity states are observed when the particle concentration is reduced.

The coherent motions of thermal active Brownian particles

Active matter exhibits many intriguing non-equilibrium characteristics, for instance, without any attractive and aligned interactions, the active Brownian particle (ABP) system undergoing motility-induced phase separation forms a high-density phase with both structural ordering and dynamical coherence. Recently, the velocity correlation among the particles in this high-density phase was found in non-thermal overdamped ABP systems. However, it seemed to disappear if thermal noises were included, bringing some confusion about the generality of the consistency between structures and dynamics in ABPs. Here, we demonstrate that the thermal noises imposing a large random term on the instantaneous velocity of ABPs hinder the observation of the inherent correlation in the motions of ABPs. By averaging the instantaneous velocity (or equivalently, calculating the displacement), we show that the inherent motions of thermal-fluctuated ABPs are highly coherent. Whether there is thermal noise or not, the inherent collective motions of ABPs do exist, and the collective motion domains are consistent spatially with the ordered clusters of ABPs in the high-density phase. At the boundary of these ordered clusters, the active forces of the particles tend to point inward and compress to sustain these clusters, thus the particles in the clusters move coherently to form some vortex-like or aligned velocity domains.

Dynamics of active particles with translational and rotational inertia

Inertial effects affecting both the translational and rotational dynamics are inherent to a broad range of active systems at the macroscopic scale. Thus, there is a pivotal need for proper models in the framework of active matter to correctly reproduce experimental results, hopefully achieving theoretical insights. For this purpose, we propose an inertial version of the active Ornstein-Uhlenbeck particle (AOUP) model accounting for particle mass (translational inertia) as well as its moment of inertia (rotational inertia) and derive the full expression for its steady-state properties. The inertial AOUP dynamics introduced in this paper is designed to capture the basic features of the well-established inertial active Brownian particle model, i.e. the persistence time of the active motion and the long-time diffusion coefficient. For a small or moderate rotational inertia, these two models predict similar dynamics at all timescales and, in general, our inertial AOUP model consistently yields the same trend upon changing the moment of inertia for various dynamical correlation functions.

Flocking without Alignment Interactions in Attractive Active Brownian Particles

Within a simple model of attractive active Brownian particles, we predict flocking behavior and challenge the widespread idea that alignment interactions are necessary to observe this collective phenomenon. Here, we show that even nonaligning attractive interactions can lead to a flocking state. Monitoring the velocity polarization as the order parameter, we reveal the onset of a first-order transition from a disordered phase, characterized by several small clusters, to a flocking phase, where a single flocking cluster is emerging. The scenario is confirmed by studying the spatial connected correlation function of particle velocities, which reveals scale-free behavior in flocking states and exponential-like decay for nonflocking configurations. Our predictions can be tested in microscopic and macroscopic experiments showing flocking, such as animals, migrating cells, and active colloids.

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.

The interplay between chemo-phoretic interactions and crowding in active colloids

Many motile microorganisms communicate with each other and their environments via chemical signaling which leads to long-range interactions mediated by self-generated chemical gradients. However, consequences of the interplay between crowding and chemotactic interactions on their collective behavior remain poorly understood. In this work, we use Brownian dynamics simulations to investigate the effect of packing fraction on the formation of non-equilibrium structures in a monolayer of diffusiophoretic self-propelled colloids as a model for chemically active particles. Focusing on the case when a chemical field induces attractive positional and repulsive orientational interactions, we explore dynamical steady-states of active colloids of varying packing fractions and degrees of motility. In addition to collapsed, active gas, and dynamical clustering steady-states reported earlier for low packing fractions, a new phase-separated state emerges. The phase separation results from a competition between long-range diffusiophoretic interactions and motility and is observed at moderate activities and a wide range of packing fractions. Our analysis suggests that the fraction of particles in the largest cluster is a suitable order parameter for capturing the transition from an active gas and dynamical clustering states to a phase-separated state.

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.

Glassy Dynamics in Chiral Fluids

Chiral active matter is enjoying a rapid increase of interest, spurred by the rich variety of asymmetries that can be attained in, e.g., the shape or self-propulsion mechanism of active particles. Though this has already led to the observance of so-called chiral crystals, active chiral glasses remain largely unexplored. A possible reason for this could be the naive expectation that interactions dominate the glassy dynamics and the details of the active motion become increasingly less relevant. Here, we show that quite the opposite is true by studying the glassy dynamics of interacting chiral active Brownian particles. We demonstrate that when our chiral fluid is pushed to glassy conditions, it exhibits highly nontrivial dynamics, especially compared to a standard linear active fluid such as common active Brownian particles. Despite the added complexity, we are still able to present a full rationalization for all identified dynamical regimes. Most notably, we introduce a new "hammering" mechanism, unique to rapidly spinning particles in high-density conditions, that can fluidize a chiral active solid.

Improving estimation of entropy production rate for run-and-tumble particle systems by high-order thermodynamic uncertainty relation

Entropy production plays an important role in the regulation and stability of active matter systems, and its rate quantifies the nonequilibrium nature of these systems. However, entropy production is hard to experimentally estimate even in some simple active systems like molecular motors or bacteria, which may be modeled by the run-and-tumble particle (RTP), a representative model in the study of active matters. Here we resolve this problem for an asymmetric RTP in one dimension, first constructing a finite-time thermodynamic uncertainty relation (TUR) for a RTP, which works well in the short observation time regime for entropy production estimation. Nevertheless, when the activity dominates, i.e., the RTP is far from equilibrium, the lower bound for entropy production from TUR turns out to be trivial. We address this issue by introducing a recently proposed high-order thermodynamic uncertainty relation (HTUR), in which the cumulant generating function of current serves as a key ingredient. To exploit the HTUR, we adopt a method to analytically obtain the cumulant generating function of the current we study, with no need to explicitly know the time-dependent probability distribution. The HTUR is demonstrated to be able to estimate the steady state energy dissipation rate accurately because the cumulant generating function covers higher-order statistics of the current, including rare and large fluctuations besides its variance. Compared to the conventional TUR, the HTUR could give significantly improved estimation of energy dissipation, which can work well even in the far from equilibrium regime. We also provide a strategy based on the improved bound to estimate the entropy production from a moderate amount of trajectory data for experimental feasibility.

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.

Critical behavior of quorum-sensing active particles

It is still a debated issue whether all critical active particles belong to the same universality class. Here we numerically study the critical behavior of quorum sensing active particles that represents the archetypal model for interpreting motility-induced phase separation. Mean-field theory predicts that this model should undergo a full phase separation if particles slow-down enough when sensing the presence of their neighbors and that the coexistence line terminates in a critical point. By performing large-scale numerical simulations, we confirm this scenario, locate the critical point and use finite-size scaling analysis to show that the static and dynamic critical exponents of this active system substantially agree with those of the Ising universality class.

Reentrant phase separation of a sparse collection of nonreciprocally aligning self-propelled disks

We study a model of aligning self-propelled disks that nonreciprocally reorient the self-propulsion directions along the interparticle separation and towards the other disks. In the limit of small inertia and large softness, where conventional motility-induced phase separation is absent, we demonstrate that the homogeneous system at a small area fraction phase-separates into clusters and a low-density phase that, eventually, reenters the homogeneous phase with a monotonic increase in alignment strength. The disks inside the clusters move with a finite space-dependent speed, constantly shuttling between clusters through the surrounding low-density homogeneous phase while maintaining the hexatic structure properties within the clusters. The area fraction gradually increases from the periphery towards the center of the clusters with a negligible correlation of the velocity and propulsion direction inside the clusters. The novel collective behavior of reentrant phase separation is found to follow from both the limits of hard disks and extremely small inertia, tending towards the overdamped limit. However, important differences in the structural and dynamical properties are shown in the limit of hard disks and extremely small inertia, as compared to that for soft disks at finite inertia. We show that the cluster phase is associated with an effective temperature for a wide range of values of alignment strength, whereas an effective temperature is associated with the specific range of alignment in the low-density phase. We believe that the reentrant phase behavior in the limit of small area fraction and the remarkable properties of the clusters should be useful in understanding a wide range of physics issues, ranging from clogging and unclogging to information exchange and transport, in biological and synthetic self-propelled systems.

Hydrodynamic effects on the liquid-hexatic transition of active colloids

We study numerically the role of hydrodynamics in the liquid-hexatic transition of active colloids at intermediate activity, where motility induced phase separation (MIPS) does not occur. We show that in the case of active Brownian particles (ABP), the critical density of the transition decreases upon increasing the particle's mass, enhancing ordering, while self-propulsion has the opposite effect in the activity regime considered. Active hydrodynamic particles (AHP), instead, undergo the liquid-hexatic transition at higher values of packing fraction [Formula: see text] than the corresponding ABP, suggesting that hydrodynamics have the net effect of disordering the system. At increasing densities, close to the hexatic-liquid transition, we found in the case of AHP the appearance of self-sustained organized motion with clusters of particles moving coherently.

Intermittent attractive interactions lead to microphase separation in nonmotile active matter

Nonmotile active matter exhibits a wide range of nonequilibrium collective phenomena yet examples are crucially lacking in the literature. We present a microscopic model inspired by the bacteria Neisseria meningitidis in which diffusive agents feel intermittent attractive forces. Through a formal coarse-graining procedure, we show that this truly scalar model of active matter exhibits the time-reversal-symmetry breaking terms defining the Active Model B+ class. In particular, we confirm the presence of microphase separation by solving the kinetic equations numerically. We show that the switching rate controlling the interactions provides a regulation mechanism tuning the typical cluster size, e.g., in populations of bacteria interacting via type IV pili.

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.

Active depletion torque between two passive rods

The active depletion torque experienced by two anisotropic objects in an active bath is a conceptional generalization of the equilibrium entropic torque. Using Brownian dynamics simulations, we compute the active depletion torque suffered by two passive rods immersed in an ensemble of active Brownian particles. Our results demonstrate that the active depletion torque is qualitatively different from its passive counterpart. Interestingly, we find that the active depletion torque can be greatly affected by the external constraint applied on the rotational degree of freedom of the rods, and even the direction may be changed with the orientational constraint, which is in contrast to the equilibrium depletion torque. The main reason for the remarkable features of the active depletion torque is that the active particles can significantly accumulate in the vicinity of the rods due to persistent self-propulsion, which is sensitively dependent on the constraint strength and the rod configurations. Our findings could be relevant for understanding the self-assembly and dynamics of anisotropic macromolecules in living environments.

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.

Biased-angle effect on diffusion dynamics and phase separation in anisotropic active particle system

A deep understanding for collective behavior in an active matter system with complex interactions has far-reaching impact in biology. In the present work, we adopt Langevin dynamics simulations to investigate diffusion dynamics and phase separation in an anisotropic active particle system with a tunable biased angle α defined as the deviation between the active force direction and anisotropic orientation. Our results demonstrate that the biased angle can induce super-rotational diffusion dynamics characterized by a power-law relationship between the mean square angle displacement (MSAD) and the time interval Δ t in the form of MSAD ∼ Δ t β with β > 1 and also result in non-trivial phase separation kinetics. As activity is dominant, nucleation time shows a non-monotonic dependence on the biased angle. Moreover, there arises a distinct transition of phase separation, from spinodal decomposition without apparent nucleation time to binodal decomposition with prominent nucleation delay. A significant inhibition effect occurs at right and obtuse angles, where the remarkable super-rotational diffusion prevents particle aggregation, leading to a slow nucleation process. As active force is competitive to anisotropic interactions, the system is almost homogeneous, while, intriguingly, we observe a re-entrant phase separation as a small acute angle is introduced. The prominent super-rotational diffusion under small angles provides an optimum condition for particle adsorption and cluster growth and, thus, accounts for the re-entrance of phase separation. A consistent scenario for the physical mechanism of our observations is achieved by properly considering the modulation of the biased angle on the interplay between activity and anisotropic interactions.

Kick effect of enzymes causes filament compression

We investigate the influence of enzymes on the structure and dynamics of a filament by dissipative particle dynamics simulations. Enzyme exerts a kick force on the filament monomer. We pay particular attention to two factors: the magnitude of kick force and enzyme concentration. Large kick force as well as high enzyme concentration prefers a remarkable compression of the filament reminiscent of the effective depletion interaction owing to an effective increase in enzyme size and the reduction of solvent quality. Additionally, the kick effect gives rise to an increase of enzyme density from the center-of-mass of the filament to its periphery. Moreover, the increase of enzyme concentration and kick force also causes a decrease in relaxation time. Our finding is helpful to understand the role of catalytic force in chemo-mechano-biological function and the filament behavior under chemical reaction via kick-induced change of solvent quality.

Competing effects of rotational diffusivity and activity on finite-sized clusters

Colloidal particles interacting via short-range attraction and long-range repulsion are known to stabilize finite-sized clusters under equilibrium conditions. In this work, the effect of self-propulsion speed (activity) and rotational diffusivity (D_r) on the phase behavior of such particles is investigated using Brownian dynamics simulations. The system exhibits rich phase behavior consisting of clusters of different kinds. The cluster size varies non-monotonically with activity: increasing first and decreasing at higher activity, thus driving cluster-to-fluid phase transition. Rotational diffusivity also facilitates the formation of clusters. Larger clusters could be stabilized at lowD_rvalues while at highD_rvalues, clusters are stable even at higher activities. The analysis of the static structure factor of the system confirms that rotational diffusivity delays the cluster-to-fluid transition driven by activity.

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.