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

Tunable assembly of confined Janus microswimmers in sub-kHz AC electric fields under gravity

Active systems comprising micron-sized self-propelling units, also termed microswimmers, are promising candidates for the bottom-up assembly of small structures and reconfigurable materials. Here we leverage field-driven colloidal assembly to induce structural transformations in dense layers of microswimmers driven by an alternating current (AC) electric field and confined in a microfabricated trap under the influence of gravity. By varying the electric field frequency, we realize significant structural transformations, from a gas-like state at high frequencies to dynamically rearranging dense crystalline clusters at lower frequencies, characterized by vorticity in their dynamics. We demonstrate the ability to switch between these states on-demand, showing that the clustering mechanism differs from motility-induced phase separation. Our results offer a valuable framework for designing high-density active matter systems with controllable structural properties, envisioned to advance the development of artificial materials with self-healing and reconfiguration capabilities.

Confined active particles with spatially dependent Lorentz force: An odd twist to the "best Fokker-Planck approximation"

We derive a version of the so-called "best Fokker-Planck approximation" (BFPA) to describe the spatial properties of interacting active Ornstein-Uhlenbeck particles in arbitrary spatial dimensions. In doing so, we also take into account the odd-diffusive contribution of the Lorentz force acting on a charged particle in a spatially dependent magnetic field, sticking to the overdamped limit. While the BFPA itself does not turn out to be widely useful, our general approach allows us to deduce an appropriate generalization of the Fox approximation, which we use to characterize the stationary behavior of a single active particle in an external potential by deriving analytic expressions for configurational probability distributions (or effective potentials). In agreement with computer simulations, our theory predicts that the Lorentz force reduces the effective attraction and thus the probability to find an active particle in the vicinity of a repulsive wall. Even for an inhomogeneous magnetic field, our theoretical findings provide useful qualitative insights, specifically regarding the location of accumulation regions.

Proposal of a quantum version of active particles via a nonunitary quantum walk

The main aim of the present paper is to define an active particle in a quantum framework as a minimal model of quantum active matter and investigate the differences and similarities of quantum and classical active matter. Although the field of active matter has been expanding, most research has been conducted on classical systems. Here, we propose a truly deterministic quantum active-particle model with a nonunitary quantum walk as the minimal model of quantum active matter. We aim to reproduce results obtained previously with classical active Brownian particles; that is, a Brownian particle, with finite energy take-up, becomes active and climbs up a potential wall. We realize such a system with nonunitary quantum walks. We introduce new internal states, the ground state | G ⟩ and the excited state | E ⟩ , and a new nonunitary operator N ( g ) N(g) for an asymmetric transition between | G ⟩ and | E ⟩ . The non-Hermiticity parameter ( g ) promotes the transition to the excited state; hence, the particle takes up energy from the environment. For our quantum active particle, we successfully observe that the movement of the quantum walker becomes more active in a nontrivial manner as we increase the non-Hermiticity parameter ( g ) , which is similar to the classical active Brownian particle. We also observe three unique features of quantum walks, namely, ballistic propagation of peaks in one dimension, the walker staying on the constant energy plane in two dimensions, and oscillations originating from the resonant transition between the ground state | G ⟩ and the excited state | E ⟩ both in one and two dimensions.

Inertial Active Matter with Coulomb Friction

Friction is central to the motion of active (self-propelled) objects such as bacteria, animals, and robots. While in a viscous fluid friction is described by Stokes's law, objects in contact with other solid bodies are often governed by more complex empirical friction laws. Here, we study active particles subject to Coulomb friction using a combination of active granular experiments and simulations, supported by theoretical predictions. The interplay of friction and activity forces induces a rich behavior resulting in three distinct dynamical regimes. While for low activity Brownian motion is recovered, for large activity we observe a dynamical stop and go regime that continuously switches from diffusion and accelerated motion. For greater activity, we observe a supermobile dynamical regime characterized by a fully accelerated motion which is described by an anomalous scaling of the diffusion coefficient with the activity. These findings cannot be observed with Stokes viscous forces typical of active swimmers but are central in dry active objects.

Dynamical crossovers and correlations in a harmonic chain of active particles

We explore the dynamics of a tracer in an active particle harmonic chain, investigating the influence of interactions. Our analysis involves calculating mean-squared displacements (MSDs) and space-time correlations through Green's function techniques and numerical simulations. Depending on chain characteristics, i.e., different time scales determined by interaction stiffness and persistence of activity, tagged-particle MSDs exhibit ballistic, diffusive, and single-file diffusion (SFD) scaling over time, with crossovers explained by our analytic expressions. Our results reveal transitions in bulk particle displacement distributions from an early-time bimodal to late-time Gaussian, passing through regimes of unimodal distributions with finite support and negative excess kurtosis and longer-tailed distributions with positive excess kurtosis. The distributions exhibit data collapse, aligning with ballistic, diffusive, and SFD scaling in the appropriate time regimes. However, at much longer times, the distributions become Gaussian. Finally, we derive analytic expressions for steady-state static and dynamic two-point displacement correlations. We verify these from simulations and highlight the differences from the equilibrium results.

The fluctuation-dissipation relation holds for a macroscopic tracer in an active bath

The fluctuation-dissipation relation (FDR) links thermal fluctuations and dissipation at thermal equilibrium through temperature. Extending it beyond equilibrium conditions in pursuit of broadening thermodynamics is often feasible, albeit with system-dependent specific conditions. We demonstrate experimentally that a generalized FDR holds for a harmonically trapped tracer colliding with self-propelled walkers. The generalized FDR remains valid across a large spectrum of active fluctuation frequencies, extending from underdamped to critically damped dynamics, which we attribute to a single primary channel for energy input and dissipation in our system.

Active Brownian information engine: Self-propulsion induced colossal performance

The information engine is a feedback mechanism that extorts work from a single heat bath using the mutual information earned during the measurement. We consider an overdamped active Ornstein-Uhlenbeck particle trapped in a 1D harmonic oscillator. The particle experiences fluctuations from an inherent thermal bath with a diffusion coefficient (D) and an active reservoir, with characteristic correlation time (τa) and strength (Da). We design a feedback-driven active Brownian information engine (ABIE) and analyze its best performance criteria. The optimal functioning criteria, the information gained during measurement, and the excess output work are reliant on the dispersion of the steady-state distribution of the particle's position. The extent of enhanced performance of such ABIE depends on the relative values of two underlying time scales of the process, namely, thermal relaxation time (τr) and the characteristic correlation time (τa). In the limit of τa/τr → 0, one can achieve the upper bound on colossal work extraction as ∼0.202γ(D+Da) (γ is the friction coefficient). The excess amount of extracted work reduces and converges to its passive counterpart (∼0.202γD) in the limit of τa/τr → high. Interestingly, when τa/τr = 1, half the upper bound of excess work is achieved irrespective of the strength of either reservoirs, thermal or active. Finally, we look into the average displacement of active Brownian particles in each feedback cycle, which surpasses its thermal analog due to the broader marginal probability distribution.

Counterion Blockade in a Heterogeneously Charged Single-File Water Channel

The Possion-Nernst-Planck theories fail to describe the ionic transport in Angstrom channels, where conduction deviates from Ohm's law, which is attributed to the dehydration/self-energy barrier and dissociation of Bjerrum ion pairs in previous work. Here, we find that the cations can be strongly bound to the surface charge, which blocks the ionic transport in a single-file water channel, causing nonlinear current-voltage curves. The presence of free ions significantly increases the probability of bound ions being released, resulting in an ionic current. We find that ionic conduction gradually becomes Ohmic as the surface charge density increases, but the conduction amplitude decreases due to the increased friction from the bound ions. We rationalize the ionic transport using 1D Kramers' escape theory framework, which describes nonlinear ionic current and the impact of surface charge density on the I-V curves. Our results show that the strong Coulomb interaction between the counterion and surface charge may cause ionic blockade in Angstrom channels.

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.

Motility driven glassy dynamics in confluent epithelial monolayers

As wounds heal, embryos develop, cancer spreads, or asthma progresses, the cellular monolayer undergoes a glass transition between solid-like jammed and fluid-like flowing states. During some of these processes, the cells undergo an epithelial-to-mesenchymal transition (EMT): they acquire in-plane polarity and become motile. Thus, how motility drives the glassy dynamics in epithelial systems is critical for the EMT process. However, no analytical framework that is indispensable for deeper insights exists. Here, we develop such a theory inspired by a well-known glass theory. One crucial result of this work is that the confluency affects the effective persistence time-scale of active force, described by its rotational diffusivity, Deffr. Deffr differs from the bare rotational diffusivity, Dr, of the motile force due to cell shape dynamics, which acts to rectify the force dynamics: Deffr is equal to Dr when Dr is small and saturates when Dr is large. We test the theoretical prediction of Deffr and how it affects the relaxation dynamics in our simulations of the active Vertex model. This novel effect of Deffr is crucial to understanding the new and previously published simulation data of active glassy dynamics in epithelial monolayers.

Statistical mechanics of passive Brownian particles in a fluctuating harmonic trap

We consider passive Brownian particles trapped in an "imperfect" harmonic trap. The trap is imperfect because it is randomly turned off and on, and as a result particles fail to equilibrate. Another way to think about this is to say that a harmonic trap is time dependent on account of its strength evolving stochastically in time. Particles in such a system are passive and activity arises through external control of a trapping potential, thus, no internal energy is used to power particle motion. A stationary Fokker-Planck equation of this system can be represented as a third-order differential equation, and its solution, a stationary distribution, can be represented as a superposition of Gaussian distributions for different strengths of a harmonic trap. This permits us to interpret a stationary system as a system in equilibrium with quenched disorder.

Extremely persistent dense active fluids

We study the dynamics of dense three-dimensional systems of active particles for large persistence times τp at constant average self-propulsion force f. These systems are fluid counterparts of previously investigated extremely persistent systems, which in the large persistence time limit relax only on the time scale of τp. We find that many dynamic properties of the systems we study, such as the mean-squared velocity, the self-intermediate scattering function, and the shear-stress correlation function, become τp-independent in the large persistence time limit. In addition, the large τp limits of many dynamic properties, such as the mean-square velocity and the relaxation times of the scattering function, and the shear-stress correlation function, depend on f as power laws with non-trivial exponents. We conjecture that these systems constitute a new class of extremely persistent active systems.

Recovering the activity parameters of an active fluid confined in a sphere

The properties of an active fluid, for example, a bacterial bath or a collection of microtubules and molecular motors, can be accessed through the dynamics of passive particle probes. Here, in the perspective of analyzing experimental situations of confinement in droplets, we consider the kinematics of a negatively buoyant probe particle in an active fluid, both confined within a spherical domain. The active bath generates a fluctuating flow that pushes the particle with a velocity that is modeled as a colored stochastic noise, characterized by two parameters, the intensity and memory time of the active flow. When the particle departs a little from the bottom of the spherical domain, the configuration is well approximated by a particle in a two-dimensional harmonic trap subjected to the colored noise, in which case an analytical solution exists, which is the base for quantitative analysis. We numerically simulate the dynamics of the particle and use the planar, two-dimensional mean square displacement to recover the activity parameters of the bath. This approach yields satisfactory results as long as the particle remains relatively confined; that is, as long as the intensity of the colored noise remains low.

Deep learning probability flows and entropy production rates in active matter

Active matter systems, from self-propelled colloids to motile bacteria, are characterized by the conversion of free energy into useful work at the microscopic scale. They involve physics beyond the reach of equilibrium statistical mechanics, and a persistent challenge has been to understand the nature of their nonequilibrium states. The entropy production rate and the probability current provide quantitative ways to do so by measuring the breakdown of time-reversal symmetry. Yet, their efficient computation has remained elusive, as they depend on the system's unknown and high-dimensional probability density. Here, building upon recent advances in generative modeling, we develop a deep learning framework to estimate the score of this density. We show that the score, together with the microscopic equations of motion, gives access to the entropy production rate, the probability current, and their decomposition into local contributions from individual particles. To represent the score, we introduce a spatially local transformer network architecture that learns high-order interactions between particles while respecting their underlying permutation symmetry. We demonstrate the broad utility and scalability of the method by applying it to several high-dimensional systems of active particles undergoing motility-induced phase separation (MIPS). We show that a single network trained on a system of 4,096 particles at one packing fraction can generalize to other regions of the phase diagram, including to systems with as many as 32,768 particles. We use this observation to quantify the spatial structure of the departure from equilibrium in MIPS as a function of the number of particles and the packing fraction.

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.

Time-dependent properties of run-and-tumble particles: Density relaxation

We characterize collective diffusion of hardcore run-and-tumble particles (RTPs) by explicitly calculating the bulk-diffusion coefficient D(ρ,γ) for arbitrary density ρ and tumbling rate γ, in systems on a d-dimensional periodic lattice. We study two minimal models of RTPs: Model I is the standard version of hardcore RTPs introduced in [Phys. Rev. E 89, 012706 (2014)10.1103/PhysRevE.89.012706], whereas model II is a long-ranged lattice gas (LLG) with hardcore exclusion, an analytically tractable variant of model I. We calculate the bulk-diffusion coefficient analytically for model II and numerically for model I through an efficient Monte Carlo algorithm; notably, both models have qualitatively similar features. In the strong-persistence limit γ→0 (i.e., dimensionless ratio r_{0}γ/v→0), with v and r_{0} being the self-propulsion speed and particle diameter, respectively, the fascinating interplay between persistence and interaction is quantified in terms of two length scales: (i) persistence length l_{p}=v/γ and (ii) a "mean free path," being a measure of the average empty stretch or gap size in the hopping direction. We find that the bulk-diffusion coefficient varies as a power law in a wide range of density: D∝ρ^{-α}, with exponent α gradually crossing over from α=2 at high densities to α=0 at low densities. As a result, the density relaxation is governed by a nonlinear diffusion equation with anomalous spatiotemporal scaling. In the thermodynamic limit, we show that the bulk-diffusion coefficient-for ρ,γ→0 with ρ/γ fixed-has a scaling form D(ρ,γ)=D^{(0)}F(ρav/γ), where a∼r_{0}^{d-1} is particle cross section and D^{(0)} is proportional to the diffusion coefficient of noninteracting particles; the scaling function F(ψ) is calculated analytically for model II (LLG) and numerically for model I. Our arguments are independent of dimensions and microscopic details.

Searching for a partially absorbing target by a run-and-tumble particle in a confined space

A random search of a partially absorbing target by a run-and-tumble particle in a confined one-dimensional space is investigated. We analytically obtain the mean searching time, which shows a nonmonotonic behavior as a function of the self-propulsion speed of the active particle, indicating the existence of an optimal speed, when the absorption strength of the target is finite. In the limit of large and small absorption strengths, respectively, asymptotes of the mean searching time and the optimal speed are found. We also demonstrate that the first-passage problem of a diffusive run-and-tumble particle in high dimensions can be mapped into a one-dimensional problem with a partially absorbing target. Finally, as a practical application exploiting the existence of the optimal speed, we propose a filtering device to extract active particles with a desired speed and evaluate how the resolution of the filtering device depends on the absorption strength.

Time-delayed Duffing oscillator in an active bath

During recent decades active particles have attracted an incipient attention as they have been observed in a broad class of scenarios, ranging from bacterial suspension in living systems to artificial swimmers in nonequilibirum systems. The main feature of these particles is that they are able to gain kinetic energy from the environment, which is widely modeled by a stochastic process due to both (Gaussian) white and Ornstein-Uhlenbeck noises. In the present work, we study the nonlinear dynamics of the forced, time-delayed Duffing oscillator subject to these noises, paying special attention to their impact upon the maximum oscillations amplitude and characteristic frequency of the steady state for different values of the time delay and the driving force. Overall, our results indicate that the role of the time delay is substantially modified with respect to the situation without noise. For instance, we show that the oscillations amplitude grows with increasing noise strength when the time delay acts as a damping term in absence of noise, whereas the trajectories eventually become aperiodic when the oscillations are sustained by the time delay. In short, the interplay among the noises, forcing, and time delay gives rise to a rich dynamics: a regular and periodic motion is destroyed or restored owing to the competition between the noise and the driving force depending on time delay values, whereas an erratic motion insensitive to the driving force emerges when the time delay makes the motion aperiodic. Interestingly, we also show that, for a sufficient noise strength and forcing amplitude, an approximately periodic interwell motion is promoted by means of stochastic resonance.

Non-Markovian Modeling of Nonequilibrium Fluctuations and Dissipation in Active Viscoelastic Biomatter

Based on a Hamiltonian that incorporates the elastic coupling between a tracer particle and the embedding active viscoelastic biomatter, we derive a generalized non-Markovian Langevin model for the nonequilibrium mechanical tracer response. Our analytical expressions for the frequency-dependent tracer response function and the tracer positional autocorrelation function agree quantitatively with experimental data for red blood cells and actomyosin networks with and without adenosine triphosphate over the entire frequency range and in particular reproduce the low-frequency violation of the fluctuation-dissipation theorem. The viscoelastic power laws, the elastic constants and effective friction coefficients extracted from the experimental data allow straightforward physical interpretation.

Non-reciprocity across scales in active mixtures

In active matter, particles typically experience mediated interactions, which are not constrained by Newton's third law and are therefore generically non-reciprocal. Non-reciprocity leads to a rich set of emerging behaviors that are hard to account for starting from the microscopic scale, due to the absence of a generic theoretical framework out of equilibrium. Here we consider bacterial mixtures that interact via mediated, non-reciprocal interactions (NRI) like quorum-sensing and chemotaxis. By explicitly relating microscopic and macroscopic dynamics, we show that, under conditions that we derive explicitly, non-reciprocity may fade upon coarse-graining, leading to large-scale equilibrium descriptions. In turn, this allows us to account quantitatively, and without fitting parameters, for the rich behaviors observed in microscopic simulations including phase separation, demixing, and multi-phase coexistence. We also derive the condition under which non-reciprocity survives coarse-graining, leading to a wealth of dynamical patterns. Again, our analytical approach allows us to predict the phase diagram of the system starting from its microscopic description. All in all, our work demonstrates that the fate of non-reciprocity across scales is a subtle and important question.

Enhanced diffusion of tracer particles in nonreciprocal mixtures

We study the diffusivity of a tagged particle in a binary mixture of Brownian particles with nonreciprocal interactions. Numerical simulations reveal that, for a broad class of interaction potentials, nonreciprocity can significantly increase the long-time diffusion coefficient of tracer particles and that this diffusion enhancement is associated with a breakdown of the Einstein relation. These observations are quantified and confirmed via two different and complementary analytical approaches: (i) a linearized stochastic density field theory, which is particularly accurate in the limit of soft interactions, and (ii) a reduced two-body description, which is exact at leading order in the density of particles. The latter reveals that diffusion enhancement can be attributed to the formation of transiently propelled dimers of particles, whose cohesion and speed are controlled by the nonreciprocal interactions.

Survival probabilities and first-passage distributions of self-propelled particles in spherical cavities

A model of self-propelled motion in a closed compartment containing simple or complex fluids is formulated in this paper in terms of the dynamics of a point particle moving in a spherical cavity under the action of random thermal forces and exponentially correlated noise. The particle's time evolution is governed by a generalized Langevin equation (GLE) in which the memory function, connected to the thermal forces by a fluctuation-dissipation relation, is described by Jeffrey's model of viscoelasticity (which reduces to a model of ordinary viscous dynamics in a suitable limit). The GLE is transformed exactly to a Fokker-Planck equation that in spherical polar coordinates is in turn found to admit of an exact solution for the particle's probability density function under absorbing boundary conditions at the surface of the sphere. The solution is used to derive an expression (that is also exact) for the survival probability of the particle in the sphere, starting from its center, which is then used to calculate the distribution of the particle's first-passage times to the boundary. The behavior of these quantities is investigated as a function of the Péclet number and the persistence time of the athermal forces, providing insight into the effects of nonequilibrium fluctuations on confined particle motion in three dimensions.

Work Fluctuations for a Harmonically Confined Active Ornstein-Uhlenbeck Particle

We study the active work fluctuations of an active Ornstein-Uhlenbeck particle in the presence of a confining harmonic potential. We tackle the problem analytically both for stationary and generic uncorrelated initial states. Our results show that harmonic confinement can induce singularities in the active work rate function, with linear stretches at large positive and negative active work, at sufficiently large active and harmonic force constants. These singularities originate from big jumps in the displacement and in the active force, occurring at the initial or ending points of trajectories and marking the relevance of boundary terms in this problem.

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.

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.

Information Engine in a Nonequilibrium Bath

Information engines can convert thermal fluctuations of a bath at temperature T into work at rates of order k_{B}T per relaxation time of the system. We show experimentally that such engines, when in contact with a bath that is out of equilibrium, can extract much more work. We place a heavy, micron-scale bead in a harmonic potential that ratchets up to capture favorable fluctuations. Adding a fluctuating electric field increases work extraction up to ten times, limited only by the strength of the applied field. Our results connect Maxwell's demon with energy harvesting and demonstrate that information engines in nonequilibrium baths can greatly outperform conventional engines.

Effects of the self-propulsion parity on the efficiency of a fuel-consuming active heat engine

We propose a thermodynamically consistent, analytically tractable model of steady-state active heat engines driven by both temperature difference and a constant chemical driving. While the engine follows the dynamics of the active Ornstein-Uhlenbeck particle, its self-propulsion stems from the mechanochemical coupling with the fuel consumption dynamics, allowing for both even- and odd-parity self-propulsion forces. Using the standard methods of stochastic thermodynamics, we show that the entropy production of the engine satisfies the conventional Clausius relation, based on which we define the efficiency of the model that is bounded from above by the second law of thermodynamics. Using this framework, we obtain exact expressions for the efficiency at maximum power. The results show that the engine performance has a nonmonotonic dependence on the magnitude of the chemical driving and that the even-parity (odd-parity) engines perform better when the size of the engine is smaller (larger) than the persistence length of the active particle. We also discuss the existence of a tighter upper bound on the efficiency of the odd-parity engines stemming from the detailed structure of the entropy production.

Dynamic phase transition induced by active molecules in a supercooled liquid

The purpose of this work is to use active particles to investigate the effect of facilitation on supercooled liquids. To this end we examine the behavior of a model supercooled liquid that is doped with a mixture of active particles and slowed particles. To simulate the facilitation mechanism, the activated particles are subjected to a force that follows the mobility of their most mobile neighboring molecule, while the slowed particles experience a friction force. Upon activation, we observe a fluidization of the entire medium along with a significant increase in dynamic heterogeneity. This effect is reminiscent of the fluidization observed experimentally when introducing molecular motors into soft materials. Interestingly, when the characteristic time τ_{μ}, used to define the mobility in the facilitation mechanism, matches the physical time t^{*} that characterizes the spontaneous cooperativity of the material, we observe a phase transition accompanied by structural aggregation of the active molecules. This transition is characterized by a sharp increase in fluidization and dynamic heterogeneity.

Stationary states of an active Brownian particle in a harmonic trap

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

Single active particle in a harmonic potential: Question about the existence of the Jarzynski relation

The interest in active matter stimulates the need to generalize thermodynamic description and relations to active matter systems, which are intrinsically out of equilibrium. One important example is the Jarzynski relation, which links the exponential average of work done in an arbitrary process connecting two equilibrium states with the difference of the free energies of these states. Using a simple model system, a single thermal active Ornstein-Uhlenbeck particle in a harmonic potential, we show that if the standard stochastic thermodynamics definition of work is used, the Jarzynski relation is not generally valid for processes connecting stationary states of active matter systems.

Most probable paths for active Ornstein-Uhlenbeck particles

Fluctuations play an important role in the dynamics of stochastic systems. In particular, for small systems, the most probable thermodynamic quantities differ from their averages because of the fluctuations. Using the Onsager Machlup variational formalism we analyze the most probable paths for nonequilibrium systems, in particular, active Ornstein-Uhlenbeck particles, and investigate how the entropy production along these paths differs from the average entropy production. We investigate how much information about their nonequilibrium nature can be obtained from their extremum paths and how these paths depend on the persistence time and their swim velocities. We also look at how the entropy production along the most probable paths varies with the active noise and how it differs from the average entropy production. This study would be useful to design artificial active systems with certain target trajectories.

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.

Boundary design regulates the diffusion of active matter in heterogeneous environments

Physical boundaries play a key role in governing the overall transport properties of nearby self-propelled particles. In this work, we develop dispersion theories and conduct Brownian dynamics simulations to predict the coupling between surface accumulation and effective diffusivity of active particles in boundary-rich media. We focus on three models that are well-understood for passive systems: particle transport in (i) an array of fixed volume-excluding obstacles; (ii) a pore with spatially heterogeneous width; and (iii) a tortuous path with kinks and corners. While the impact of these entropic barriers on passive particle transport is well established, we find that these classical models of porous media flows break down due to the unique interplay between activity and the microstructure of the internal geometry. We study the activity-induced slowdown of effective diffusivity by formulating a Smoluchowski description of long-time self diffusivity which contains contributions from the density and fluctuation fields of the active particles. Particle-based and finite element simulations corroborate this perspective and reveal important nonequilibrium considerations of active transport.

Most probable path of active Ornstein-Uhlenbeck particles

Using the path integral representation of the nonequilibrium dynamics, we compute the most probable path between arbitrary starting and final points that is followed by an active particle driven by persistent noise. We focus our attention on the case of active particles immersed in harmonic potentials, where the trajectory can be computed analytically. Once we consider the extended Markovian dynamics where the self-propulsive drive evolves according to an Ornstein-Uhlenbeck process, we can compute the trajectory analytically with arbitrary conditions on position and self-propulsion velocity. We test the analytical predictions against numerical simulations and we compare the analytical results with those obtained within approximated equilibriumlike dynamics.

Tuning nonequilibrium phase transitions with inertia

In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation-dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. This density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.

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.

Configurational temperature in active matter. I. Lines of invariant physics in the phase diagram of the Ornstein-Uhlenbeck model

This paper shows that the configurational temperature of liquid-state theory, T_{conf}, defines an energy scale, which can be used for adjusting model parameters of active Ornstein-Uhlenbeck particle (AOUP) models in order to achieve approximately invariant structure and dynamics upon a density change. The required parameter changes are calculated from the variation of a single configuration's T_{conf} for a uniform scaling of all particle coordinates. The resulting equations are justified theoretically for models involving a potential-energy function with hidden scale invariance. The validity of the procedure is illustrated by computer simulations of the Kob-Andersen binary Lennard-Jones AOUP model, showing the existence of lines of approximate invariance of the reduced-unit radial distribution function and time-dependent mean-square displacement.

Mode-coupling theory for the dynamics of dense underdamped active Brownian particle system

We present a theory to study the inertial effect on glassy dynamics of the underdamped active Brownian particle (UABP) system. Using the assumption of the nonequilibrium steady-state, we obtain an effective Fokker-Planck equation for the probability distribution function (PDF) as a function of positions and momentums. With this equation, we achieve the evolution equation of the intermediate scattering function through the Zwanzig-Mori projection operator method and the mode-coupling theory (MCT). Theoretical analysis shows that the inertia of the particle affects the memory function and corresponding glass transition by influencing the structure factor and a velocity correlation function. The theory provides theoretical support and guidance for subsequent simulation work.

Entropy production of active particles formulated for underdamped dynamics

The present work investigates the effect of inertia on the entropy production rate Π for all canonical models of active particles for different dimensions and the type of confinement. To calculate Π, the link between the entropy production and dissipation of heat rate is explored, resulting in a simple and intuitive expression. By analyzing the Kramers equation, alternative formulations of Π are obtained and the virial theorem for active particles is derived. Exact results are obtained for particles in an unconfined environment and in a harmonic trap. In both cases, Π is independent of temperature. For the case of a harmonic trap, Π attains a maximal value for τ=ω^{-1}, where τ is the persistence time and ω is the natural frequency of an oscillator. For active particles in one-dimensional box, or other nonharmonic potentials, thermal fluctuations are found to reduce Π.

On the origin of universal cell shape variability in confluent epithelial monolayers

Cell shape is fundamental in biology. The average cell shape can influence crucial biological functions, such as cell fate and division orientation. But cell-to-cell shape variability is often regarded as noise. In contrast, recent works reveal that shape variability in diverse epithelial monolayers follows a nearly universal distribution. However, the origin and implications of this universality remain unclear. Here, assuming contractility and adhesion are crucial for cell shape, characterized via aspect ratio (r), we develop a mean-field analytical theory for shape variability. We find that all the system-specific details combine into a single parameter α that governs the probability distribution function (PDF) of r; this leads to a universal relation between the standard deviation and the average of r. The PDF for the scaled r is not strictly but nearly universal. In addition, we obtain the scaled area distribution, described by the parameter μ. Information of α and μ together can distinguish the effects of changing physical conditions, such as maturation, on different system properties. We have verified the theory via simulations of two distinct models of epithelial monolayers and with existing experiments on diverse systems. We demonstrate that in a confluent monolayer, average shape determines both the shape variability and dynamics. Our results imply that cell shape distribution is inevitable, where a single parameter describes both statics and dynamics and provides a framework to analyze and compare diverse epithelial systems. In contrast to existing theories, our work shows that the universal properties are consequences of a mathematical property and should be valid in general, even in the fluid regime.

Harvesting information to control nonequilibrium states of active matter

We propose to use a correlated noise bath to drive an optically trapped Brownian particle that mimics active biological matter. Due to the flexibility and precision of our setup, we are able to control the different parameters that drive the stochastic motion of the particle with unprecedented accuracy, thus reaching strongly correlated regimes that are not easily accessible with real active matter. In particular, by using the correlation time (i.e., the "color") of the noise as a control parameter, we can trigger transitions between two nonequilibrium steady states with no expended work, but only a calorific cost. Remarkably, the measured heat production is directly proportional to the spectral entropy of the correlated noise, in a fashion that is reminiscent of Landauer's principle. Our procedure can be viewed as a method for harvesting information from the active fluctuations.

Active Refrigerators Powered by Inertia

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

Anomalous Cooling and Overcooling of Active Colloids

The phenomenon that a system at a hot temperature cools faster than at a warm temperature, referred to as the Mpemba effect, has recently been realized for trapped colloids. Here, we investigate the cooling and heating process of a self-propelled active colloid using numerical simulations and theoretical calculations with a model that can be directly tested in experiments. Upon cooling, activity induces a Mpemba effect and the active particle transiently escapes an effective temperature description. At the end of the cooling process the notion of temperature is recovered and the system can exhibit even smaller temperatures than its final temperature, a surprising phenomenon which we refer to as activity-induced overcooling.

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.

Inertial active Ornstein-Uhlenbeck particle in the presence of a magnetic field

We consider an inertial active Ornstein-Uhlenbeck particle in an athermal bath. The particle is charged, constrained to move in a two-dimensional harmonic trap, and a magnetic field is applied perpendicular to the plane of motion. The steady-state correlations and the mean-square displacement are studied when the particle is confined as well as when it is set free from the trap. With the help of both numerical simulation and analytical calculations, we observe that inertia plays a crucial role in the dynamics in the presence of a magnetic field. In a highly viscous medium where the inertial effects are negligible, the magnetic field has no influence on the correlated behavior of position as well as velocity. In the time asymptotic limit, the overall displacement of the confined harmonic particle gets enhanced by the presence of a magnetic field and saturates for a stronger magnetic field. On the other hand, when the particle is set free, the overall displacement gets suppressed and approaches zero when the strength of the field is very high. Interestingly, it is seen that in the time asymptotic limit, the confined harmonic particle behaves like a passive particle and becomes independent of the activity, especially in the presence of a very strong magnetic field. Similarly, for a free particle the mean-square displacement in the long time limit becomes independent of activity even for a longer persistence of noise cor- relation in the dynamics.

From predicting to learning dissipation from pair correlations of active liquids

Active systems, which are driven out of equilibrium by local non-conservative forces, can adopt unique behaviors and configurations. Towards designing such materials, an important challenge is to precisely connect the static structure of active systems to the dissipation of energy induced by the local driving. Here, we use tools from liquid-state theories and machine learning to take on these challenges. We first demonstrate analytically for an isotropic active matter system that dissipation and pair correlations are closely related when driving forces behave like an active temperature. We then extend a nonequilibrium mean-field framework for predicting these pair correlations which, unlike most existing approaches, is applicable even for strongly interacting particles and far from equilibrium, to predict dissipation in these systems. Based on this theory, we reveal analytically a robust relation between dissipation and structure which holds even as the system approaches a nonequilibrium phase transition. Finally, we construct a neural network which maps static configurations of particles to their dissipation rate without any prior knowledge of the underlying dynamics. Our results open novel perspectives on the interplay between dissipation and organization out-of-equilibrium.

Mean-field theory for the structure of strongly interacting active liquids

Active systems, which are driven out of equilibrium by local non-conservative forces, exhibit unique behaviors and structures with potential utility for the design of novel materials. An important and difficult challenge along the path towards such a goal is to precisely predict how the structure of active systems is modified as their driving forces push them out of equilibrium. Here, we use tools from liquid-state theories to approach this challenge for a classic minimal isotropic active matter model. First, we construct a nonequilibrium mean-field framework which can predict the structure of systems of weakly interacting particles. Second, motivated by equilibrium solvation theories, we modify this theory to extend it with surprisingly high accuracy to strongly interacting particles, distinguishing it from most existing similarly tractable approaches. Our results provide insight into spatial organization in strongly interacting out-of-equilibrium systems and strategies to control them.

Motion of a self-propelled particle with rotational inertia

Overdamped active Brownian motion of self-propelled particles in a liquid has been fairly well studied. However, there are a variety of situations in which the overdamped approximation is not justified, for instance, when self-propelled particles move in a low-viscosity medium or when their rotational diffusivity is enhanced by internal active processes or external control. Examples of various origins include biofilaments driven by molecular motors, living and artificial microflyers and interfacial surfers, field-controlled and superfluid microswimmers, vibration-driven granular particles and autonomous mini-robots with sensorial delays, etc. All of them extend active Brownian motion to the underdamped case, i.e., to active Langevin motion, which takes into account inertia. Despite a rich experimental background, there is a gap in the theory in the field where rotational inertia significantly affects the random walk of active particles on all time scales. In particular, although the well-known models of active Brownian and Ornstein-Uhlenbeck particles include a memory effect of the direction of motion, they are not applicable in the underdamped case, because the rotational inertia, which they do not account for, can partially prevent "memory loss" with increasing rotational diffusion. We describe the two-dimensional motion of a self-propelled particle with both translational and rotational inertia and velocity fluctuations. The proposed generalized analytical equations for the mean kinetic energy, mean-square displacement and noise-averaged trajectory of the self-propelled particle are confirmed by numerical simulations in a wide range of self-propulsion velocities, moments of inertia, rotational diffusivities, medium viscosities and observation times.