Self-propelled particle in an external potential: existence of an effective temperature

Generalized Ornstein-Uhlenbeck model for active motion

We investigate a one-dimensional model of active motion, which takes into account the effects of persistent self-propulsion through a memory function in a dissipative-like term of the generalized Langevin equation for particle swimming velocity. The proposed model is a generalization of the active Ornstein-Uhlenbeck model introduced by G. Szamel [Phys. Rev. E 90, 012111 (2014)10.1103/PhysRevE.90.012111]. We focus on two different kinds of memory which arise in many natural systems: an exponential decay and a power law, supplemented with additive colored noise. We provide analytical expressions for the velocity autocorrelation function and the mean-squared displacement, which are in excellent agreement with numerical simulations. For both models, damped oscillatory solutions emerge due to the competition between the memory of the system and the persistence of velocity fluctuations. In particular, for a power-law model with fractional Brownian noise, we show that long-time active subdiffusion occurs with increasing long-term memory.

On the fluctuation-dissipation relation in non-equilibrium and non-Hamiltonian systems

We review generalized fluctuation-dissipation relations, which are valid under general conditions even in "nonstandard systems," e.g., out of equilibrium and/or without a Hamiltonian structure. The response functions can be expressed in terms of suitable correlation functions computed in the unperturbed dynamics. In these relations, typically, one has nontrivial contributions due to the form of the stationary probability distribution; such terms take into account the interaction among the relevant degrees of freedom in the system. We illustrate the general formalism with some examples in nonstandard cases, including driven granular media, systems with a multiscale structure, active matter, and systems showing anomalous diffusion.

Active particles in noninertial frames: How to self-propel on a carousel

Typically the motion of self-propelled active particles is described in a quiescent environment establishing an inertial frame of reference. Here we assume that friction, self-propulsion, and fluctuations occur relative to a noninertial frame and thereby the active Brownian motion model is generalized to noninertial frames. First, analytical solutions are presented for the overdamped case, both for linear swimmers and for circle swimmers. For a frame rotating with constant angular velocity ("carousel"), the resulting noise-free trajectories in the static laboratory frame are trochoids if these are circles in the rotating frame. For systems governed by inertia, such as vibrated granulates or active complex plasmas, centrifugal and Coriolis forces become relevant. For both linear and circling self-propulsion, these forces lead to out-spiraling trajectories which for long times approach a spira mirabilis. This implies that a self-propelled particle will typically leave a rotating carousel. A navigation strategy is proposed to avoid this expulsion, by adjusting the self-propulsion direction at will. For a particle, initially quiescent in the rotating frame, it is shown that this strategy only works if the initial distance to the rotation center is smaller than a critical radius R_{c} which scales with the self-propulsion velocity. Possible experiments to verify the theoretical predictions are discussed.

Mode-coupling theory for the steady-state dynamics of active Brownian particles

We present a theory for the steady-state dynamics of a two-dimensional system of spherically symmetric active Brownian particles. The derivation of the theory consists of two steps. First, we integrate out the self-propulsions and obtain a many-particle evolution equation for the probability distribution of the particles' positions. Second, we use the projection operator technique and a mode-coupling-like factorization approximation to derive an equation of motion for the density correlation function. The nonequilibrium character of the active system manifests itself through the presence of a steady-state correlation function that quantifies spatial correlations of microscopic steady-state currents of the particles. This function determines the dependence of the short-time dynamics on the activity. It also enters into the expression for the memory matrix and thus influences the long-time glassy dynamics.

Active chiral particles under confinement: surface currents and bulk accumulation phenomena

In this work, we study the stationary behavior of an assembly of independent chiral active particles under confinement by employing an extension of the active Ornstein-Uhlenbeck model. The chirality modeled by means of an effective torque term leads to a drastic reduction in the accumulation near the walls with respect to the case without handedness and to the appearance of currents parallel to the container walls accompanied by a large accumulation of particles in the inner region. In the case of two-dimensional chiral particles confined by harmonic walls, we determine the analytic form of the distribution of positions and velocities in two different situations: a rotationally invariant confining potential and an infinite channel with parabolic walls. Both these models display currents and chirality induced inner accumulation. These phenomena are further investigated by means of a more realistic description of a channel, where the wall and bulk regions are clearly separated. The corresponding current and density profiles are obtained by numerical simulations. At variance with the harmonic models, the third model shows a progressive emptying of the wall regions and the simultaneous enhancement of the bulk population. We explain such a phenomenon in terms of the combined effect of wall repulsive forces and chiral motion and provide a semiquantitative description of the current profile in terms of effective viscosity of the chiral gas.

Activity induced delocalization and freezing in self-propelled systems

We study a system of interacting active particles, propelled by colored noises, characterized by an activity time τ, and confined by a single-well anharmonic potential. We assume pair-wise repulsive forces among particles, modelling the steric interactions among microswimmers. This system has been experimentally studied in the case of a dilute suspension of Janus particles confined through acoustic traps. We observe that already in the dilute regime - when inter-particle interactions are negligible - increasing the persistent time, τ, pushes the particles away from the potential minimum, until a saturation distance is reached. We compute the phase diagram (activity versus interaction length), showing that the interaction does not suppress this delocalization phenomenon but induces a liquid- or solid-like structure in the densest regions. Interestingly a reentrant behavior is observed: a first increase of τ from small values acts as an effective warming, favouring fluidization; at higher values, when the delocalization occurs, a further increase of τ induces freezing inside the densest regions. An approximate analytical scheme gives fair predictions for the density profiles in the weakly interacting case. The analysis of non-equilibrium heat fluxes reveals that in the region of largest particle concentration equilibrium is restored in several aspects.

Approximating microswimmer dynamics by active Brownian motion: Energetics and efficiency

We consider the dynamics of a microswimmer and show that they can be approximated by active Brownian motion. The swimmer is modeled by coupled overdamped Langevin equations with periodic driving. We compare the energy dissipation of the real swimmer to that of the active Brownian motion model, finding that the latter can massively underestimate the complete dissipation. This discrepancy is related to the inability to infer the full dissipation from partial observation of the complete system. We introduce an efficiency that measures how much of the dissipated energy is spent on forward propulsion.

Stationary superstatistics distributions of trapped run-and-tumble particles

We present an analysis of the stationary distributions of run-and-tumble particles trapped in external potentials in terms of a thermophoretic potential that emerges when trapped active motion is mapped to trapped passive Brownian motion in a fictitious inhomogeneous thermal bath. We elaborate on the meaning of the non-Boltzmann-Gibbs stationary distributions that emerge as a consequence of the persistent motion of active particles. These stationary distributions are interpreted as a class of distributions in nonequilibrium statistical mechanics known as superstatistics. Our analysis provides an original insight on the link between the intrinsic nonequilibrium nature of active motion and the well-known concept of local equilibrium used in nonequilibrium statistical mechanics and contributes to the understanding of the validity of the concept of effective temperature. Particular cases of interest, regarding specific trapping potentials used in other theoretical or experimental studies, are discussed. We point out as an unprecedented effect, the emergence of new modes of the stationary distribution as a consequence of the coupling of persistent motion in a trapping potential that varies highly enough with position.

Active particles under confinement and effective force generation among surfaces

We consider the effect of geometric confinement on the steady-state properties of a one-dimensional active suspension subject to thermal noise. The random active force is modeled by an Ornstein-Uhlenbeck process and the system is studied both numerically, by integrating the Langevin governing equations, and analytically by solving the associated Fokker-Planck equation under suitable approximations. The comparison between the two approaches displays a fairly good agreement and in particular, we show that the Fokker-Planck approach can predict the structure of the system both in the wall region and in the bulk-like region where the surface forces are negligible. The simultaneous presence of thermal noise and active forces determines the formation of a layer, extending from the walls towards the bulk, where the system exhibits polar order. We relate the presence of such ordering to the mechanical pressure exerted on the container's walls and show how it depends on the separation of the boundaries and determines a Casimir-like attractive force mediated by the active suspension.

Active dumbbells: Dynamics and morphology in the coexisting region

With the help of molecular dynamics simulations we study an ensemble of active dumbbells in purely repulsive interaction. We derive the phase diagram in the density-activity plane and we characterise the various phases with liquid, hexatic and solid character. The analysis of the structural and dynamical properties, such as enstrophy, mean-square displacement, polarisation, and correlation functions, shows the continuous character of liquid and hexatic phases in the coexisting region when the activity is increased starting from the passive limit.

Hidden entropy production and work fluctuations in an ideal active gas

Collections of self-propelled particles that move persistently by continuously consuming free energy are a paradigmatic example of active matter. In these systems, unlike Brownian "hot colloids," the breakdown of detailed balance yields a continuous production of entropy at steady state, even for an ideal active gas. We quantify the irreversibility for a noninteracting active particle in two dimensions by treating both conjugated and time-reversed dynamics. By starting with underdamped dynamics, we identify a hidden rate of entropy production required to maintain persistence and prevent the rapidly relaxing momenta from thermalizing, even in the limit of very large friction. Additionally, comparing two popular models of self-propulsion with identical dissipation on average, we find that the fluctuations and large deviations in work done are markedly different, providing thermodynamic insight into the varying extents to which macroscopically similar active matter systems may depart from equilibrium.

Complex structures generated by competing interactions in harmonically confined colloidal suspensions

We investigate the structural properties of colloidal particle systems interacting via an isotropic pair potential and confined by a three-dimensional harmonic potential. The interaction potential has a repulsive-attractive-repulsive profile that varies with the interparticle distance (also known as a 'mermaid' potential). We performed Langevin dynamics simulations to find the equilibrium configurations of the system. We show that particles can self-assemble in complex structural patterns, such as compact disks, fringed disks, rods, spherical clusters with superficial entrances among others. Also, for particular values of the parameters of the interaction potential, we could identify that some configurations were formed by quasi two-dimensional (2D) structures which are stable for 2D systems.

Active processes in one dimension

We consider the thermal and athermal overdamped motion of particles in one-dimensional geometries where discrete internal degrees of freedom (spin) are coupled with the translational motion. Adding a driving velocity that depends on the time-dependent spin constitutes the simplest model of active particles (run-and-tumble processes) where the violation of the equipartition principle and of the Sutherland-Einstein relation can be studied in detail even when there is generalized reversibility. We give an example (with four spin values) where the irreversibility of the translational motion manifests itself only in higher-order (than two) time correlations. We derive a generalized telegraph equation as the Smoluchowski equation for the spatial density for an arbitrary number of spin values. We also investigate the Arrhenius exponential law for run-and-tumble particles; due to their activity the slope of the potential becomes important in contrast to the passive diffusion case and activity enhances the escape from a potential well (if that slope is high enough). Finally, in the absence of a driving velocity, the presence of internal currents such as in the chemistry of molecular motors may be transmitted to the translational motion and the internal activity is crucial for the direction of the emerging spatial current.

Activated barrier crossing dynamics of a Janus particle carrying cargo

We numerically study the escape kinetics of a self-propelled Janus particle, carrying a cargo, from a meta-stable state.

Effective equilibrium states in mixtures of active particles driven by colored noise

We consider the steady-state behavior of pairs of active particles having different persistence times and diffusivities. To this purpose we employ the active Ornstein-Uhlenbeck model, where the particles are driven by colored noises with exponential correlation functions whose intensities and correlation times vary from species to species. By extending Fox's theory to many components, we derive by functional calculus an approximate Fokker-Planck equation for the configurational distribution function of the system. After illustrating the predicted distribution in the solvable case of two particles interacting via a harmonic potential, we consider systems of particles repelling through inverse power-law potentials. We compare the analytic predictions to computer simulations for such soft-repulsive interactions in one dimension and show that at linear order in the persistence times the theory is satisfactory. This work provides the toolbox to qualitatively describe many-body phenomena, such as demixing and depletion, by means of effective pair potentials.

Memory effects in active particles with exponentially correlated propulsion

The Ornstein-Uhlenbeck particle (OUP) model imagines a microscopic swimmer propelled by an active force which is correlated with itself on a finite time scale. Here we investigate the influence of external potentials on an ideal suspension of OUPs, in both one and two spatial dimensions, with particular attention paid to the pressure exerted on "confining walls." We employ a mathematical connection between the local density of OUPs and the statistics of their propulsion force to demonstrate the existence of an equation of state in one dimension. In higher dimensions we show that active particles generate a nonconservative force field in the surrounding medium. A simplified far-from-equilibrium model is proposed to account for OUP behavior in the vicinity of potentials. Building on this, we interpret simulations of OUPs in more complicated situations involving asymmetrical and spatially curved potentials, and characterize the resulting inhomogeneous stresses in terms of competing active length scales.

Entropy Production and Fluctuation Theorems for Active Matter

Active biological systems reside far from equilibrium, dissipating heat even in their steady state, thus requiring an extension of conventional equilibrium thermodynamics and statistical mechanics. In this Letter, we have extended the emerging framework of stochastic thermodynamics to active matter. In particular, for the active Ornstein-Uhlenbeck model, we have provided consistent definitions of thermodynamic quantities such as work, energy, heat, entropy, and entropy production at the level of single, stochastic trajectories and derived related fluctuation relations. We have developed a generalization of the Clausius inequality, which is valid even in the presence of the non-Hamiltonian dynamics underlying active matter systems. We have illustrated our results with explicit numerical studies.

Mode-coupling theory for active Brownian particles

We present a mode-coupling theory (MCT) for the slow dynamics of two-dimensional spherical active Brownian particles (ABPs). The ABPs are characterized by a self-propulsion velocity v_{0} and by their translational and rotational diffusion coefficients D_{t} and D_{r}, respectively. Based on the integration-through-transients formalism, the theory requires as input only the equilibrium static structure factors of the passive system (where v_{0}=0). It predicts a nontrivial idealized-glass-transition diagram in the three-dimensional parameter space of density, self-propulsion velocity, and rotational diffusivity that arise because at high densities, the persistence length of active swimming ℓ_{p}=v_{0}/D_{r} interferes with the interaction length ℓ_{c} set by the caging of particles. While the low-density dynamics of ABPs is characterized by a single Péclet number Pe=v_{0}^{2}/D_{r}D_{t}, close to the glass transition the dynamics is found to depend on Pe and ℓ_{p} separately. At fixed density, increasing the self-propulsion velocity causes structural relaxation to speed up, while decreasing the persistence length slows down the relaxation. The active-MCT glass is a nonergodic state that is qualitatively different from the passive glass. In it, correlations of initial density fluctuations never fully decay, but also an infinite memory of initial orientational fluctuations is retained in the positions.

Active Brownian equation of state: metastability and phase coexistence

As a result of the competition between self-propulsion and excluded volume interactions, purely repulsive self-propelled spherical particles undergo a motility-induced phase separation (MIPS). We carry out a systematic computational study, considering several interaction potentials, systems confined by hard walls or with periodic boundary conditions, and different initial conditions. This approach allows us to identify that, despite its non-equilibrium nature, the equations of state of Active Brownian Particles (ABP) across MIPS verify the characteristic properties of first-order liquid-gas phase transitions, meaning, equality of pressure of the coexisting phases once a nucleation barrier has been overcome and, in the opposite case, hysteresis around the transition as long as the system remains in the metastable region. Our results show that the equations of state of ABPs account for their phase behaviour, providing a firm basis to describe MIPS as an equilibrium-like phase transition.

Pressure and flow of exponentially self-correlated active particles

Microscopic swimming particles, which dissipate energy to execute persistent directed motion, are a classic example of a nonequilibrium system. We investigate the noninteracting Ornstein-Uhlenbeck Particle (OUP), which is propelled through a viscous medium by a force which is correlated over a finite time. We obtain an exact expression for the steady-state phase-space density of a single OUP confined by a quadratic potential, and use the result to explore more complex geometries, both through analytical approximations and numerical simulations. In a "Casimir"-style setup involving two narrowly spaced walls, we describe a particle-trapping phenomenon, which leads to a repulsive effective interaction between the walls, while in a two-dimensional annulus geometry, we observe net stresses which resemble the Laplace pressure.

Nonequilibrium mode-coupling theory for dense active systems of self-propelled particles

The physics of active systems of self-propelled particles, in the regime of a dense liquid state, is an open puzzle of great current interest, both for statistical physics and because such systems appear in many biological contexts. We develop a nonequilibrium mode-coupling theory (MCT) for such systems, where activity is included as a colored noise with the particles having a self-propulsion force f₀ and a persistence time τ_p. Using the extended MCT and a generalized fluctuation-dissipation theorem, we calculate the effective temperature T_eff of the active fluid. The nonequilibrium nature of the systems is manifested through a time-dependent T_eff that approaches a constant in the long-time limit, which depends on the activity parameters f₀ and τ_p. We find, phenomenologically, that this long-time limit is captured by the potential energy of a single, trapped active particle (STAP). Through a scaling analysis close to the MCT glass transition point, we show that τ_α, the α-relaxation time, behaves as τ_α ∼ f₀^-2γ, where γ = 1.74 is the MCT exponent for the passive system. τ_α may increase or decrease as a function of τ_p depending on the type of active force correlations, but the behavior is always governed by the same value of the exponent γ. Comparison with the numerical solution of the nonequilibrium MCT and simulation results give excellent agreement with scaling analysis.

Mode coupling theory for nonequilibrium glassy dynamics of thermal self-propelled particles

We present a mode coupling theory study for the relaxation and glassy dynamics of a system of strongly interacting self-propelled particles, wherein the self-propulsion force is described by Ornstein-Uhlenbeck colored noise and thermal noises are included. Our starting point is an effective Smoluchowski equation governing the distribution function of particle positions, from which we derive a memory function equation for the time dependence of density fluctuations in nonequilibrium steady states. With the basic assumption of the absence of macroscopic currents and standard mode coupling approximation, we can obtain expressions for the irreducible memory function and other relevant dynamic terms, wherein the nonequilibrium character of the active system is manifested through an averaged diffusion coefficient D[combining macron] and a nontrivial structural function S₂(q) with q being the magnitude of wave vector q. D[combining macron] and S₂(q) enter the frequency term and the vertex term for the memory function, and thus influence both the short time and the long time dynamics of the system. With these equations obtained, we study the glassy dynamics of this thermal self-propelled particle system by investigating the Debye-Waller factor f_q and relaxation time τ_α as functions of the persistence time τ_p of self-propulsion, the single particle effective temperature T_eff as well as the number density ρ. Consequently, we find the critical density ρ_c for given τ_p shifts to larger values with increasing magnitude of propulsion force or effective temperature, in good accordance with previously reported simulation work. In addition, the theory facilitates us to study the critical effective temperature T for fixed ρ as well as its dependence on τ_p. We find that T increases with τ_p and in the limit τ_p → 0, it approaches the value for a simple passive Brownian system as expected. Our theory also well recovers the results for passive systems and can be easily extended to more complex systems such as active-passive mixtures.

Directional transport of colloids inside a bath of self-propelling walkers

We present a setup in which passive colloids inside a solvent are moved to the boundaries of the container. The directional transport is facilitated by self-propelling microparticles ("walkers") with an activity gradient, which reduces their propulsion in the vicinity of bounding walls. An attractive interaction leads to the adsorption of walkers onto the colloid-surfaces in regions of low walker activity. It is shown that the activity gradient generates a free energy gradient which in turn acts as a driving force on the passive colloids. We carry out molecular dynamics simulations and present approaches to a theoretical description of the involved processes. Although the simulation data are not reproduced on a fully quantitative level, their qualitative features are covered by the model. The effect described here may be applied to facilitate a directional transport of drugs or to eliminate pollutants.

Driven Brownian particle as a paradigm for a nonequilibrium heat bath: Effective temperature and cyclic work extraction

We apply the concept of a frequency-dependent effective temperature based on the fluctuation-dissipation ratio to a driven Brownian particle in a nonequilibrium steady state. Using this system as a thermostat for a weakly coupled harmonic oscillator, the oscillator thermalizes according to a canonical distribution at the respective effective temperature across the entire frequency spectrum. By turning the oscillator from a passive thermometer into a heat engine, we realize the cyclic extraction of work from a single thermal reservoir, which is feasible only due to its nonequilibrium nature.

Study of dynamic heterogeneity of an active particle system

We have studied spatial and temporal dynamic heterogeneity (DH) in a system of hard-sphere particles, subjected to active forces with constant amplitude and random direction determined by rotational diffusion with correlation time τ. We have used a variety of observables to characterize the DH behavior, including the deviation from standard Stokes-Einstein (SE) relation, a non-Gaussian parameter α_{2}(Δt) for the distribution of particle displacement within a certain time interval Δt, a four-point susceptibility χ_{4}(Δt,ΔL) for the correlation in dynamics between any two points in space separated by distance ΔL within some time window Δt, and a vector spatial-temporal correlation function S_{vec}(R,Δt) for vector displacements within time interval Δt of particle pairs originally separated by R. By mapping the particle motion into a continuous-time random walk with constant jump length, we can obtain the average waiting time 〈t_{x}〉∝D_{s}^{-1} and persistence time 〈t_{p}〉∝η, with D_{s} the self-diffusion coefficient and η the shear viscosity, such that the observable λ=〈t_{p}〉/〈t_{x}〉∝D_{s}η can be calculated as a function of the control parameter τ to show how it deviates from its SE value λ_{0}. Interestingly, we find λ/λ_{0} shows a nonmonotonic behavior for large volume fraction φ_{a}, wherein λ/λ_{0} undergoes a minimum at a certain intermediate value of τ, indicating that both small and large particle activity may lead to strong DH. Such a reentrance phenomenon is further demonstrated in terms of the non-Gaussian parameters α_{2}, four-point susceptibility χ_{4}, and vector spatiotemporal correlation functions S_{vec}, respectively. Detail analysis shows that it is the competition between the dual roles of particle activity, namely, activity-induced higher effective temperature and activity-induced clustering, that leads to such nontrivial nonmonotonic behaviors. In addition, we find that DH may also show a maximum level at an intermediate value of φ_{a} if τ is large enough, implying that a more crowded system may be less heterogeneous than a less crowded one for a system with high particle activity.

Crystallization and flow in active patch systems

Based upon recent experiments in which Janus particles are made into active swimmers by illuminating them with laser light, we explore the effect of applying a light pattern on the sample, thereby creating activity inducing zones or active patches. We simulate a system of interacting Brownian diffusers that become active swimmers when moving inside an active patch and analyze the structure and dynamics of the ensuing stationary state. We find that, in some respects, the effect of spatially inhomogeneous activity is qualitatively similar to a temperature gradient. For asymmetric patches, however, this analogy breaks down because the ensuing stationary state is specific to partial active motion.

Brownian motion of a circle swimmer in a harmonic trap

We study the dynamics of a Brownian circle swimmer with a time-dependent self-propulsion velocity in an external temporally varying harmonic potential. For several situations, the noise-free swimming paths, the noise-averaged mean trajectories, and the mean-square displacements are calculated analytically or by computer simulation. Based on our results, we discuss optimal swimming strategies in order to explore a maximum spatial range around the trap center. In particular, we find a resonance situation for the maximum escape distance as a function of the various frequencies in the system. Moreover, the influence of the Brownian noise is analyzed by comparing noise-free trajectories at zero temperature with the corresponding noise-averaged trajectories at finite temperature. The latter reveal various complex self-similar spiral or rosette-like patterns. Our predictions can be tested in experiments on artificial and biological microswimmers under dynamical external confinement.

Escape rate of active particles in the effective equilibrium approach

The escape rate of a Brownian particle over a potential barrier is accurately described by the Kramers theory. A quantitative theory explicitly taking the activity of Brownian particles into account has been lacking due to the inherently out-of-equilibrium nature of these particles. Using an effective equilibrium approach [Farage et al., Phys. Rev. E 91, 042310 (2015)PLEEE81539-375510.1103/PhysRevE.91.042310] we study the escape rate of active particles over a potential barrier and compare our analytical results with data from direct numerical simulation of the colored noise Langevin equation. The effective equilibrium approach generates an effective potential that, when used as input to Kramers rate theory, provides results in excellent agreement with the simulation data.

Critical phenomena in active matter

We investigate the effect of self-propulsion on a mean-field order-disorder transition. Starting from a φ^{4} scalar field theory subject to an exponentially correlated noise, we exploit the unified colored-noise approximation to map the nonequilibrium active dynamics onto an effective equilibrium one. This allows us to follow the evolution of the second-order critical point as a function of the noise parameters: the correlation time τ and the noise strength D. Our results suggest that the universality class of the model remains unchanged. We also estimate the effect of Gaussian fluctuations on the mean-field approximation finding an Ornstein-Zernike-like expression for the static structure factor at long wavelengths. Finally, to assess the validity of our predictions, we compare the mean-field theoretical results with numerical simulations of active Lennard-Jones particles in two and three dimensions, finding good qualitative agreement at small τ values.

The nonequilibrium glassy dynamics of self-propelled particles

We study the glassy dynamics taking place in dense assemblies of athermal active particles that are driven solely by a nonequilibrium self-propulsion mechanism. Active forces are modeled as an Ornstein-Uhlenbeck stochastic process, characterized by a persistence time and an effective temperature, and particles interact via a Lennard-Jones potential that yields well-studied glassy behavior in the Brownian limit, which is obtained as the persistence time vanishes. By increasing the persistence time, the system departs more strongly from thermal equilibrium and we provide a comprehensive numerical analysis of the structure and dynamics of the resulting active fluid. Finite persistence times profoundly affect the static structure of the fluid and give rise to nonequilibrium velocity correlations that are absent in thermal systems. Despite these nonequilibrium features, for any value of the persistence time we observe a nonequilibrium glass transition as the effective temperature is decreased. Surprisingly, increasing departure from thermal equilibrium is found to promote (rather than suppress) the glassy dynamics. Overall, our results suggest that with increasing persistence time, microscopic properties of the active fluid change quantitatively, but the general features of the nonequilibrium glassy dynamics observed with decreasing the effective temperature remain qualitatively similar to those of thermal glass-formers.

Linear response to leadership, effective temperature, and decision making in flocks

Large collections of autonomously moving agents, such as animals or micro-organisms, are able to flock coherently in space even in the absence of a central control mechanism. While the direction of the flock resulting from this critical behavior is random, this can be controlled by a small subset of informed individuals acting as leaders of the group. In this article we use the Vicsek model to investigate how flocks respond to leadership and make decisions. Using a combination of numerical simulations and continuous modeling we demonstrate that flocks display a linear response to leadership that can be cast in the framework of the fluctuation-dissipation theorem, identifying an effective temperature reflecting how promptly the flock reacts to the initiative of the leaders. The linear response to leadership also holds in the presence of two groups of informed individuals with competing interests, indicating that the flock's behavioral decision is determined by both the number of leaders and their degree of influence.

How Far from Equilibrium Is Active Matter?

Active matter systems are driven out of thermal equilibrium by a lack of generalized Stokes-Einstein relation between injection and dissipation of energy at the microscopic scale. We consider such a system of interacting particles, propelled by persistent noises, and show that, at small but finite persistence time, their dynamics still satisfy a time-reversal symmetry. To do so, we compute perturbatively their steady-state measure and show that, for short persistent times, the entropy production rate vanishes. This endows such systems with an effective fluctuation-dissipation theorem akin to that of thermal equilibrium systems. Last, we show how interacting particle systems with viscous drags and correlated noises can be seen as in equilibrium with a viscoelastic bath but driven out of equilibrium by nonconservative forces, hence providing energetic insight into the departure of active systems from equilibrium.

Configurational entropy and effective temperature in systems of active Brownian particles

We propose a method to determine the effective density of states and configurational entropy in systems of active Brownian particles by measuring the probability distribution function of potential energy at varying temperatures. Assuming that the entropy is a continuous and monotonically increasing function of energy, we provide support that two-dimensional systems of purely repulsive active Brownian spheres can be mapped onto an equilibrium system with a Boltzmann-like distribution and an effective temperature. We find that the effective temperature depends even for a large number of particles on system size, suggesting that active systems are non-extensive. In addition, the effective Helmholtz free energy can be derived from the configurational entropy. We verify our results regarding the configurational entropy by using thermodynamic integration of the effective Helmholtz free energy with respect to temperature.

Theory for the dynamics of dense systems of athermal self-propelled particles

We present a derivation of a recently proposed theory for the time dependence of density fluctuations in stationary states of strongly interacting, athermal, self-propelled particles. The derivation consists of two steps. First, we start from the equation of motion for the joint distribution of particles' positions and self-propulsions and we integrate out the self-propulsions. In this way we derive an approximate, many-particle evolution equation for the probability distribution of the particles' positions. Second, we use this evolution equation to describe the time dependence of steady-state density correlations. We derive a memory function representation of the density correlation function and then we use a factorization approximation to obtain an approximate expression for the memory function. In the final equation of motion for the density correlation function the nonequilibrium character of the active system manifests itself through the presence of a new steady-state correlation function that quantifies spatial correlations of the velocities of the particles. This correlation function enters into the frequency term, and thus it describes the dependence of the short-time dynamics on the properties of the self-propulsions. More importantly, the correlation function of particles' velocities enters into the vertex of the memory function and through the vertex it modifies the long-time glassy dynamics.

Effect of self-propulsion on equilibrium clustering

In equilibrium, colloidal suspensions governed by short-range attractive and long-range repulsive interactions form thermodynamically stable clusters. Using Brownian dynamics computer simulations, we investigate how this equilibrium clustering is affected when such particles are self-propelled. We find that the clustering process is stable under self-propulsion. For the range of interaction parameters studied and at low particle density, the cluster size increases with the speed of self-propulsion (activity) and for higher activity the cluster size decreases, showing a nonmonotonic variation of cluster size with activity. This clustering behavior is distinct from the pure kinetic (or motility-induced) clustering of self-propelling particles which is observed at significantly higher activities and densities. We present an equilibrium model incorporating the effect of activity as activity-induced attraction and repulsion by imposing that the strength of these interactions depend on activity superlinearly. The model explains the cluster size dependence of activity obtained from simulations semiquantitatively. Our predictions are verifiable in experiments on interacting synthetic colloidal microswimmers.

Glassy dynamics of athermal self-propelled particles: Computer simulations and a nonequilibrium microscopic theory

We combine computer simulations and analytical theory to investigate the glassy dynamics in dense assemblies of athermal particles evolving under the sole influence of self-propulsion. Our simulations reveal that when the persistence time of the self-propulsion is increased, the local structure becomes more pronounced, whereas the long-time dynamics first accelerates and then slows down. We explain these surprising findings by constructing a nonequilibrium microscopic theory that yields nontrivial predictions for the glassy dynamics. These predictions are in qualitative agreement with the simulations and reveal the importance of steady-state correlations of the local velocities to the nonequilibrium dynamics of dense self-propelled particles.

Rotational and translational diffusion in an interacting active dumbbell system

We study the dynamical properties of a two-dimensional ensemble of self-propelled dumbbells with only repulsive interactions. This model undergoes a phase transition between a homogeneous and a segregated phase and we focus on the former. We analyze the translational and rotational mean-square displacements in terms of the Péclet number, describing the relative role of active forces and thermal fluctuations, and of particle density. We find that the four distinct regimes of the translational mean-square displacement of the single active dumbbell survive at finite density for parameters that lead to a separation of time scales. We establish the Péclet number and density dependence of the diffusion constant in the last diffusive regime. We prove that the ratio between the diffusion constant and its value for the single dumbbell depends on temperature and active force only through the Péclet number at all densities explored. We also study the rotational mean-square displacement proving the existence of a rich behavior with intermediate regimes only appearing at finite density. The ratio of the rotational late-time diffusion constant and its vanishing density limit depends on the Péclet number and density only. At low Péclet number it is a monotonically decreasing function of density. At high Péclet number it first increases to reach a maximum and then decreases as a function of density. We interpret the latter result advocating the presence of large-scale fluctuations close to the transition, at large-enough density, that favor coherent rotation inhibiting, however, rotational motion for even larger packing fractions.

Gaussian memory in kinematic matrix theory for self-propellers

We extend the kinematic matrix ("kinematrix") formalism [Phys. Rev. E 89, 062304 (2014)], which via simple matrix algebra accesses ensemble properties of self-propellers influenced by uncorrelated noise, to treat Gaussian correlated noises. This extension brings into reach many real-world biological and biomimetic self-propellers for which inertia is significant. Applying the formalism, we analyze in detail ensemble behaviors of a 2D self-propeller with velocity fluctuations and orientation evolution driven by an Ornstein-Uhlenbeck process. On the basis of exact results, a variety of dynamical regimes determined by the inertial, speed-fluctuation, orientational diffusion, and emergent disorientation time scales are delineated and discussed.

Dynamics of a homogeneous active dumbbell system

We analyze the dynamics of a two-dimensional system of interacting active dumbbells. We characterize the mean-square displacement, linear response function, and deviation from the equilibrium fluctuation-dissipation theorem as a function of activity strength, packing fraction, and temperature for parameters such that the system is in its homogeneous phase. While the diffusion constant in the last diffusive regime naturally increases with activity and decreases with packing fraction, we exhibit an intriguing nonmonotonic dependence on the activity of the ratio between the finite-density and the single-particle diffusion constants. At fixed packing fraction, the time-integrated linear response function depends nonmonotonically on activity strength. The effective temperature extracted from the ratio between the integrated linear response and the mean-square displacement in the last diffusive regime is always higher than the ambient temperature, increases with increasing activity, and, for small active force, monotonically increases with density while for sufficiently high activity it first increases and next decreases with the packing fraction. We ascribe this peculiar effect to the existence of finite-size clusters for sufficiently high activity and density at the fixed (low) temperatures at which we worked. The crossover occurs at lower activity or density the lower the external temperature. The finite-density effective temperature is higher (lower) than the single dumbbell one below (above) a crossover value of the Péclet number.