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

Socioeconomic agents as active matter in nonequilibrium Sakoda-Schelling models

How robust are socioeconomic agent-based models with respect to the details of the agents' decision rule? We tackle this question by considering an occupation model in the spirit of the Sakoda-Schelling model, historically introduced to shed light on segregation dynamics among human groups. For a large class of utility functions and decision rules, we pinpoint the nonequilibrium nature of the agent dynamics, while recovering an equilibrium-like phase separation phenomenology. Within the mean-field approximation we show how the model can be mapped, to some extent, onto an active matter field description. Finally, we consider nonreciprocal interactions between two populations and show how they can lead to nonsteady macroscopic behavior. We believe our approach provides a unifying framework to further study geography-dependent agent-based models, notably paving the way for joint consideration of population and price dynamics within a field theoretic approach.

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.

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.

Growth and aging in a few phase-separating active matter systems

Via computer simulations we study evolution dynamics in systems of continuously moving active Brownian particles. The obtained results are discussed against those from the passive 2D Ising case. Following sudden quenches of random configurations to state points lying within the miscibility gaps and to the critical points, we investigate the far-from-steady-state dynamics by calculating quantities associated with structure and characteristic length scales. We also study aging for quenches into the miscibility gap and provide a quantitative picture for the scaling behavior of the two-time order-parameter correlation function. The overall structure and dynamics are consistent with expectations from the Ising model. This remains true for certain active lattice models as well, for which we present results for quenches to the critical points.

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.

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.

Signatures of irreversibility in microscopic models of flocking

Flocking in d=2 is a genuine nonequilibrium phenomenon for which irreversibility is an essential ingredient. We study a class of minimal flocking models whose only source of irreversibility is self-propulsion and use the entropy production rate (EPR) to quantify the departure from equilibrium across their phase diagrams. The EPR is maximal in the vicinity of the order-disorder transition, where reshuffling of the interaction network is fast. We show that signatures of irreversibility come in the form of asymmetries in the steady-state distribution of the flock's microstates. These asymmetries occur as consequences of the time-reversal symmetry breaking in the considered self-propelled systems, independently of the interaction details. In the case of metric pairwise forces, they reduce to local asymmetries in the distribution of pairs of particles. This study suggests a possible use of pair asymmetries both to quantify the departure from equilibrium and to learn relevant information about aligning interaction potentials from data.

Disordered Collective Motion in Dense Assemblies of Persistent Particles

We explore the emergence of nonequilibrium collective motion in disordered nonthermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, characterized by remarkable velocity correlations and irregular turbulent flows. For all persistence values, the active fluid undergoes a nonequilibrium glass transition at large density. This is accompanied by collective motion, whose nature evolves from near-equilibrium spatially heterogeneous dynamics at small persistence, to a qualitatively different intermittent dynamics when persistence is large. This latter regime involves a complex time evolution of the correlated displacement field.

Boundary fluctuation dynamics of a phase-separated domain in planar geometry

Using phase-ordering kinetics and of renormalization group theories, we derive analytically the relaxation times of the long wavelength fluctuations of a phase-separated domain boundary in the vicinity of (and below) the critical temperature in the planar Ising universality class. For a conserved order parameter, the relaxation time grows like Λ^{3} at wavelength Λ and can be expressed in terms of parameters relevant at the microscopic scale: lattice spacing, bulk diffusion coefficient of the minority phase, and temperature. These results are supported by numerical simulations of 2D Ising models, enabling in addition calculating the nonuniversal numerical prefactor. We discuss the applications of these findings to the determination of the real timescale associated with elementary Monte Carlo moves from the measurement of long wavelength relaxation times on experimental systems or molecular dynamics simulations.

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.

Analytical approach to chiral active systems: suppressed phase separation of interacting Brownian circle swimmers

We consider chirality in active systems by exemplarily studying the phase behavior of planar systems of interacting Brownian circle swimmers with a spherical shape. For this purpose, we derive a predictive field theory that is able to describe the collective dynamics of circle swimmers. The theory yields a mapping between circle swimmers and noncircling active Brownian particles and predicts that the angular propulsion of the particles leads to a suppression of their motility-induced phase separation, being in line with recent simulation results. In addition, the theory provides analytical expressions for the spinodal corresponding to the onset of motility-induced phase separation and the associated critical point as well as for their dependence on the angular propulsion of the circle swimmers. We confirm our findings by Brownian dynamics simulations. The agreement between results from theory and simulations is found to be good.

Long-Range Influence of a Pump on a Critical Fluid

A pump coupled to a conserved density generates long-range modulations, resulting from the non-equilibrium nature of the dynamics. We study how these modulations are modified at the critical point where the system exhibits intrinsic long-range correlations. To do so, we consider a pump in a diffusive fluid, which is known to generate a density profile in the form of an electric dipole potential and a current in the form of a dipolar field above the critical point. We demonstrate that while the current retains its form at the critical point, the density profile changes drastically. At criticality, in d<4 dimensions, the deviation of the density from the average is given by sgn[cos(θ)]|cos(θ)/r^{(d-1)}|^{1/δ} at large distance r from the pump and angle θ with respect to the pump's orientation. At short distances, there is a crossover to a cos(θ)/r^{d-3+η} profile. Here δ and η are Ising critical exponents. The effect of the local pump on the domain wall structure below the critical point is also considered.

Scaling of the entropy production rate in a φ^{4} model of active matter

In active φ^{4} field theories the nonequilibrium terms play an important role in describing active phase separation; however, they are irrelevant, in the renormalization group sense, at the critical point. Their irrelevance makes the critical exponents the same as those of the Ising universality class. Despite their irrelevance, they contribute to a nontrivial scaling of the entropy production rate at criticality. We consider the nonequilibrium dynamics of a nonconserved scalar field φ (Model A) driven out-of-equilibrium by a persistent noise that is correlated on a finite timescale τ, as in the case of active baths. We perform the computation of the density of entropy production rate σ and we study its scaling near the critical point. We find that similar to the case of active Model A, and although the nonlinearities responsible for nonvanishing entropy production rates in the two models are quite different, the irrelevant parameter τ makes the critical dynamics irreversible.

Dynamics of active particles with space-dependent swim velocity

We study the dynamical properties of an active particle subject to a swimming speed explicitly depending on the particle position. The oscillating spatial profile of the swim velocity considered in this paper takes inspiration from experimental studies based on Janus particles whose speed can be modulated by an external source of light. We suggest and apply an appropriate model of an active Ornstein Uhlenbeck particle (AOUP) to the present case. This allows us to predict the stationary properties, by finding the exact solution of the steady-state probability distribution of particle position and velocity. From this, we obtain the spatial density profile and show that its form is consistent with the one found in the framework of other popular models. The reduced velocity distribution highlights the emergence of non-Gaussianity in our generalized AOUP model which becomes more evident as the spatial dependence of the velocity profile becomes more pronounced. Then, we focus on the time-dependent properties of the system. Velocity autocorrelation functions are studied in the steady-state combining numerical and analytical methods derived under suitable approximations. We observe a non-monotonic decay in the temporal shape of the velocity autocorrelation function which depends on the ratio between the persistence length and the spatial period of the swim velocity. In addition, we numerically and analytically study the mean square displacement and the long-time diffusion coefficient. The ballistic regime, observed in the small-time region, is deeply affected by the properties of the swim velocity landscape which induces also a crossover to a sub-ballistic but superdiffusive regime for intermediate times. Finally, the long-time diffusion coefficient decreases as the amplitude of the swim velocity oscillations increases because the diffusion is mainly determined by those regions where the particles are slow.

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.

Correlated escape of active particles across a potential barrier

We study the dynamics of one-dimensional active particles confined in a double-well potential, focusing on the escape properties of the system, such as the mean escape time from a well. We first consider a single-particle both in near and far-from-equilibrium regimes by varying the persistence time of the active force and the swim velocity. A non-monotonic behavior of the mean escape time is observed with the persistence time of the activity, revealing the existence of an optimal choice of the parameters favoring the escape process. For small persistence times, a Kramers-like formula with an effective potential obtained within the unified colored noise approximation is shown to hold. Instead, for large persistence times, we developed a simple theoretical argument based on the first passage theory, which explains the linear dependence of the escape time with the persistence of the active force. In the second part of the work, we consider the escape on two active particles mutually repelling. Interestingly, the subtle interplay of active and repulsive forces may lead to a correlation between particles, favoring the simultaneous jump across the barrier. This mechanism cannot be observed in the escape process of two passive particles. Finally, we find that in the small persistence regime, the repulsion favors the escape, such as in passive systems, in agreement with our theoretical predictions, while for large persistence times, the repulsive and active forces produce an effective attraction, which hinders the barrier crossing.

Doi-Peliti field theory of free active Ornstein-Uhlenbeck particles

We derive a Doi-Peliti field theory for free active Ornstein-Uhlenbeck particles, or, equivalently, free inertial Brownian particles, and present a way to diagonalize the quadratic part of the action and calculate the propagator. Unlike previous coarse-grained approaches this formulation correctly tracks particle identity and yet can easily be expanded to include potentials and arbitrary reactions.

Collective effects in confined active Brownian particles

We investigate a two-dimensional system of active particles confined to a narrow annular domain. Despite the absence of explicit interactions among the velocities or the active forces of different particles, the system displays a transition from a disordered and stuck state to an ordered state of global collective motion where the particles rotate persistently clockwise or anticlockwise. We describe this behavior by introducing a suitable order parameter, the velocity polarization, measuring the global alignment of the particles' velocities along the tangential direction of the ring. We also measure the spatial velocity correlation function and its correlation length to characterize the two states. In the rotating phase, the velocity correlation displays an algebraic decay that is analytically predicted together with its correlation length, while in the stuck regime, the velocity correlation decays exponentially with a correlation length that increases with the persistence time. In the first case, the correlation (and, in particular, its correlation length) does not depend on the active force but the system size only. The global collective motion, an effect caused by the interplay between finite-size, periodicity, and persistent active forces, disappears as the size of the ring becomes infinite, suggesting that this phenomenon does not correspond to a phase transition in the usual thermodynamic sense.

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.