Spatial velocity correlations in inertial systems of active Brownian particles

Emergent mesoscale correlations in active solids with noisy chiral dynamics

We present the linear response theory for an elastic solid composed of active Brownian particles with intrinsic individual chirality, deriving both a normal mode formulation and a continuum elastic formulation. Using this theory, we compute analytically the velocity correlations and energy spectra under different conditions, showing an excellent agreement with simulations. We generate the corresponding phase diagram, identifying chiral and achiral disordered regimes (for high chirality or noise levels), as well as chiral and achiral states with mesoscopic-range order (for low chirality and noise). The chiral ordered states display mesoscopic spatial correlations and oscillating time correlations, but no wave propagation. In the high chirality regime, we find a peak in the elastic energy spectrum that leads to a non-monotonic behavior with increasing noise strength that is consistent with the emergence of the 'hammering state' recently identified in chiral glasses. Finally, we show numerically that our theory, despite its linear response nature, can be applied beyond the idealized homogeneous solid assumed in our derivations. Indeed, by increasing the level of activity, we show that it remains a good approximation of the system dynamics until just below the melting transition. In addition, we show that there is still an excellent agreement between our analytical results and simulations when we extend our results to heterogeneous solids composed of mixtures of active particles with different intrinsic chirality and noise levels. The derived linear response theory is therefore robust and applicable to a broad range of real-world active systems. Our work provides a thorough analytical and numerical description of the emergent states in a densely packed system of chiral self-propelled Brownian disks, thus allowing a detailed understanding of the phases and dynamics identified in a minimal chiral active system.

Microscopic theory for nonequilibrium correlation functions in dense active fluids

One of the key hallmarks of dense active matter in the liquid, supercooled, and solid phases is the so-called equal-time velocity correlations. Crucially, these correlations can emerge spontaneously, i.e., they require no explicit alignment interactions, and therefore represent a generic feature of dense active matter. This indicates that for a meaningful comparison or possible mapping between active and passive liquids one not only needs to understand their structural properties, but also the impact of these velocity correlations. This has already prompted several simulation and theoretical studies, though they are mostly focused on athermal systems and thus overlook the effect of translational diffusion. Here, we present a fully microscopic method to calculate nonequilibrium correlations in two-dimensional systems of thermal active Brownian particles (ABPs). We use the integration through transients formalism together with (active) mode-coupling theory and analytically calculate qualitatively consistent static structure factors and active velocity correlations. We complement our theoretical results with simulations of both thermal and athermal ABPs which exemplify the disruptive role that thermal noise has on velocity correlations.

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

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

Time-dependent properties of run-and-tumble particles. II. Current fluctuations

We investigate steady-state current fluctuations in two models of hardcore run-and-tumble particles (RTPs) on a periodic one-dimensional lattice of L sites, for arbitrary tumbling rate γ=τ_{p}^{-1} and density ρ; model I consists of standard hardcore RTPs, while model II is an analytically tractable variant of model I, called a long-ranged lattice gas (LLG). We show that, in the limit of L large, the fluctuation of cumulative current Q_{i}(T,L) across the ith bond in a time interval T≫1/D grows first subdiffusively and then diffusively (linearly) with T: 〈Q_{i}^{2}〉∼T^{α} with α=1/2 for 1/D≪T≪L^{2}/D and α=1 for T≫L^{2}/D, where D(ρ,γ) is the collective- or bulk-diffusion coefficient; at small times T≪1/D, exponent α depends on the details. Remarkably, regardless of the model details, the scaled bond-current fluctuations D〈Q_{i}^{2}(T,L)〉/2χL≡W(y) as a function of scaled variable y=DT/L^{2} collapse onto a universal scaling curve W(y), where χ(ρ,γ) is the collective particle mobility. In the limit of small density and tumbling rate, ρ,γ→0, with ψ=ρ/γ fixed, there exists a scaling law: The scaled mobility γ^{a}χ(ρ,γ)/χ^{(0)}≡H(ψ) as a function of ψ collapses onto a scaling curve H(ψ), where a=1 and 2 in models I and II, respectively, and χ^{(0)} is the mobility in the limiting case of a symmetric simple exclusion process; notably, the scaling function H(ψ) is model dependent. For model II (LLG), we calculate exactly, within a truncation scheme, both the scaling functions, W(y) and H(ψ). We also calculate spatial correlation functions for the current and compare our theory with simulation results of model I; for both models, the correlation functions decay exponentially, with correlation length ξ∼τ_{p}^{1/2} diverging with persistence time τ_{p}≫1. Overall, our theory is in excellent agreement with simulations and complements the prior findings [T. Chakraborty and P. Pradhan, Phys. Rev. E 109, 024124 (2024)1539-375510.1103/PhysRevE.109.024124].

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.

Stiffening and dynamics of a two-dimensional active elastic solid

This work deals with the mechanical properties and dynamics of an active elastic solid defined as a two-dimensional network of active stochastic particles interacting by nonlinear hard springs. By proposing a discrete model, it is numerically found that when activity in the system is turned on, the active solid stiffens as a function of propulsion forces, thus deviating from equilibrium mechanics. To understand this effect a minimal stochastic model is offered, and a physical explanation based on spatial symmetry-breaking is put forward. In addition, the dynamics of the active solid in the absence of an external stress is also studied. From this, three main features are observed to emerge, namely, a collective behavior within the active solid, a time-density fluctuation, and oscillating dynamics of the internal stresses towards a steady state.

Free and enclosed inertial active gas

In this work, the free expansion of an inertial active gas in three dimensions made of spherical non-interactive active Brownian particles with both translational and rotational inertia (IABPs) is studied. After elucidating the active particles' orientational correlation in three dimensions by employing a Fokker-Planck formalism, the diffusion, mean-square speed, persistence length, reorientation time, Swim and Reynolds pressures and total pressure of this system, are obtained theoretically and corroborated by performing Langevin dynamics simulations. Afterwards, a numerical study on particles' distribution and the mechanical pressure exerted by the active gas enclosed in a cubic box and its dependence on inertia is also carried out. This experiment highlights two important observations: first, as inertia in the system grows while fixing activity, a more uniform particle distribution within the box is achieved. In other words, the classical accumulation of active particles at the walls is seen to be suppressed by inertia. Second, an active gas with translational and rotational inertiae and made of spherical particles still has a state equation which is offered here. This is supported by the fact that both the mechanical pressure definition and the bulk pressure definition as the trace of the swim and Reynolds stress tensors, coincide in the thermodynamic limit.

Entropons as collective excitations in active solids

The vibrational dynamics of solids is described by phonons constituting basic collective excitations in equilibrium crystals. Here, we consider a non-equilibrium active solid, formed by self-propelled particles, which bring the system into a non-equilibrium steady-state. We identify novel vibrational collective excitations of non-equilibrium (active) origin, which coexist with phonons and dominate over them when the system is far from equilibrium. These vibrational excitations are interpreted in the framework of non-equilibrium physics, in particular, stochastic thermodynamics. We call them "entropons" because they are the modes of spectral entropy production (at a given frequency and wave vector). The existence of entropons could be verified in future experiments on dense self-propelled colloidal Janus particles and granular active matter, as well as in living systems, such as dense cell monolayers.

The coherent motions of thermal active Brownian particles

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

Dynamics of active particles with translational and rotational inertia

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

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

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

Phase behavior of active and passive dumbbells

We report phase separation in a mixture of "hot" and "cold" three-dimensional dumbbells which interact by Lennard-Jones potential. We also have studied the effect of asymmetry of dumbbells and the variation of ratio of "hot" and "cold" dumbbells on their phase separation. The ratio of the temperature difference between hot and cold dumbbells to the temperature of cold dumbbells is a measure of the activity χ of the system. From constant density simulations of symmetric dumbbells, we observe that the "hot" and "cold" dumbbells phase separate at higher activity ratio (χ>5.80) compared to that of a mixture of hot and cold Lennard-Jones monomers (χ>3.44). We find that, in the phase-separated system, the hot dumbbells have high effective volume and hence high entropy which is calculated by two-phase thermodynamic method. The high kinetic pressure of hot dumbbells forces the cold dumbbells to form dense clusters such that at the interface the high kinetic pressure of hot dumbbells is balanced by the virial pressure of cold dumbbells. We find that phase separation pushes the cluster of cold dumbbells to have solidlike ordering. Bond orientation order parameters reveal that the cold dumbbells form solidlike ordering consisting of predominantly face-centered cubic and hexagonal-close packing packing, but the individual dumbbells have random orientations. The simulation of the nonequilibrium system of symmetric dumbbells at different ratios of number of hot dumbbells to cold dumbbells reveals that the critical activity of phase separation decreases with increase in fraction of hot dumbbells. The simulation of equal mixture of hot and cold asymmetric dumbbells revealed that the critical activity of phase separation was independent of the asymmetry of dumbbells. We also observed that the clusters of cold asymmetric dumbbells showed both crystalline and noncrystalline order depending on the asymmetry of dumbbells.

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.

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

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

Hydrodynamic effects on the liquid-hexatic transition of active colloids

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

Transition in collective motion decision making

Collective decision making in a biological motion group requires fast and robust transmission of information. Typically, directional switching information propagation across the whole group obeys a linear dispersion law. However, conventional dynamic collective motion models, such as the Vicsek model and the Couzin model did not take into account ultrafast directional synchronous motions. In the present paper, a multiparticle model is proposed based on inertial spin self-propel action, which can provide adequate description of such group motion. By considering both spin mechanism and collision avoidance, the proposed self-propelled particle spin model can nicely describe collective motion with fast directional switching. By analyzing the order parameter of the group-velocity synchronization, a mechanism of group decision making is revealed, which is based on the difference between two clusters of divergent leaders, showing a transition from the compromising phase (i.e., following the group average) to the preferred phase (i.e., aligning to a leader cluster). The finding provides new insight to the decision-making process of followers when they face with divergent leaders in group motion.

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.

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

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

Inertial particle under active fluctuations: Diffusion and work distributions

We study the underdamped motion of a passive particle in an active environment. Using the phase space path integral method we find the probability distribution function of position and velocity for a free and a harmonically bound particle. The environment is characterized by an active noise which is described as the Ornstein-Uhlenbeck process (OUP). Taking two similar, yet slightly different OUP models, it is shown how inertia along with other relevant parameters affect the dynamics of the particle. Further we investigate the work fluctuations of a harmonically trapped particle by considering the trap center being pulled at a constant speed. Finally, the fluctuation theorem of work is validated with an effective temperature in the steady-state limit.

Dynamics of active particles with space-dependent swim velocity

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

Collective motion of run-and-tumble repulsive and attractive particles in one-dimensional systems

Active matter deals with systems whose particles consume energy at the individual level in order to move. To unravel features such as the emergence of collective structures, several models have been suggested, such as the on-lattice model of run-and-tumble particles implemented via the persistent exclusion process (PEP). In our work, we study a one-dimensional system of run-and-tumble repulsive or attractive particles, both on-lattice and off-lattice. Additionally, we implement cluster motility dynamics in the on-lattice case (since in the off-lattice case, cluster motility arises from the individual particle dynamics). While we observe important differences between discrete and continuous dynamics, few common features are of particular importance. Increasing particle density drives aggregation across all different systems explored. For non-attractive particles, the effects of particle activity on aggregation are largely independent of the details of the dynamics. In contrast, once attractive interactions are introduced, the steady-state, which is completely determined by the interplay between these and the particles' activity, becomes highly dependent on the details of the dynamics.

Alignment destabilizes crystal order in active systems

We combine numerical and analytical methods to study two-dimensional active crystals formed by permanently linked swimmers and with two distinct alignment interactions. The system admits a stationary phase with quasi-long-range translational order, as well as a moving phase with quasi-long-range active force director and velocity order. The translational order in the moving phase is significantly influenced by alignment interaction. For Vicsek-like alignment, the translational order is short ranged, whereas the bond-orientational order is quasi-long ranged, implying a moving hexatic phase. For elasticity-based alignment, the translational order is quasi-long ranged parallel to the motion and short ranged in the perpendicular direction, whereas the bond orientational order is long ranged. We also generalize these results to higher dimensions.

Alignment-mediated segregation in an active-passive mixture

We report segregation between the athermal active and passive particles mediated by the local alignment interaction in a confined space. The competition between the alignment interaction and self-propulsion force results in a transition between disordered and ordered phases. We show that as the coordination between the particles increases, they form an ordered mill, which helps the particles to aggregate into isotropic clusters. As a result, particles segregate into active core and passive shells. This segregation phenomenon is adversely affected by the packing fraction and the size dispersion between active and passive particles. We show that this adverse effect can be overcome by incorporating higher coordination in the system. We report that the monodispersed system is more desirable for segregation in a binary mixture than a bidispersed system, as the latter favors the mixed state.

Active Ornstein-Uhlenbeck model for self-propelled particles with inertia

Self-propelled particles, which convert energy into mechanical motion, exhibit inertia if they have a macroscopic size or move inside a gaseous medium, in contrast to micron-sized overdamped particles immersed in a viscous fluid. Here we study an extension of the active Ornstein-Uhlenbeck model, in which self-propulsion is described by colored noise, to access these inertial effects. We summarize and discuss analytical solutions of the particle's mean-squared displacement and velocity autocorrelation function for several settings ranging from a free particle to various external influences, like a linear or harmonic potential and coupling to another particle via a harmonic spring. Taking into account the particular role of the initial particle velocity in a nonstationary setup, we observe all dynamical exponents between zero and four. After the typical inertial time, determined by the particle's mass, the results inherently revert to the behavior of an overdamped particle with the exception of the harmonically confined systems, in which the overall displacement is enhanced by inertia. We further consider an underdamped model for an active particle with a time-dependent mass, which critically affects the displacement in the intermediate time-regime. Most strikingly, for a sufficiently large rate of mass accumulation, the particle's motion is completely governed by inertial effects as it remains superdiffusive for all times.

Orbital magnetism of an active particle in viscoelastic suspension

We consider an active (self-propelling) particle in a viscoelastic fluid. The particle is charged and constrained to move in a two-dimensional harmonic trap. Its dynamics is coupled to a constant magnetic field applied perpendicular to its plane of motion via Lorentz force. Due to the finite activity, the generalized fluctuation-dissipation relation (GFDR) breaks down, driving the system away from equilibrium. While breaking GFDR, we have shown that the system can have finite classical orbital magnetism only when the dynamics of the system contains finite inertia. The orbital magnetic moment has been calculated exactly. Remarkably, we find that when the elastic dissipation timescale of the medium is larger (smaller) than the persistence timescale of the self-propelling particle, it is diamagnetic (paramagnetic). Therefore, for a given strength of the magnetic field, the system undergoes a transition from diamagnetic to paramagnetic state (and vice versa) simply by tuning the timescales of underlying physical processes, such as active fluctuations and viscoelastic dissipation. Interestingly, we also find that the magnetic moment, which vanishes at equilibrium, behaves nonmonotonically with respect to increasing persistence of self-propulsion, which drives the system out of equilibrium.

Kinetic Monte Carlo Algorithms for Active Matter Systems

We study kinetic Monte Carlo (KMC) descriptions of active particles. We show that, when they rely on purely persistent, active steps, their continuous-time limit is ill-defined, leading to the vanishing of trademark behaviors of active matter such as the motility-induced phase separation, ratchet effects, as well as to a diverging mechanical pressure. We then show how, under an appropriate scaling, mixing passive steps with active ones leads to a well-defined continuous-time limit that however differs from standard active dynamics. Finally, we propose new KMC algorithms whose continuous-time limits lead to the dynamics of active Ornstein-Uhlenbeck, active Brownian, and run-and-tumble particles.

Phase behavior and surface tension of soft active Brownian particles

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

Collective effects in confined active Brownian particles

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

Time-dependent inertia of self-propelled particles: The Langevin rocket

Many self-propelled objects are large enough to exhibit inertial effects but still suffer from environmental fluctuations. The corresponding basic equations of motion are governed by active Langevin dynamics, which involve inertia, friction, and stochastic noise for both the translational and orientational degrees of freedom coupled via the self-propulsion along the particle orientation. In this paper, we generalize the active Langevin model to time-dependent parameters and explicitly discuss the effect of time-dependent inertia for achiral and chiral particles. Realizations of this situation are manifold, ranging from minirockets (which are self-propelled by burning their own mass), to dust particles in plasma (which lose mass by evaporating material), to walkers with expiring activity. Here we present analytical solutions for several dynamical correlation functions, such as mean-square displacement and orientational and velocity autocorrelation functions. If the parameters exhibit a slow power law in time, we obtain anomalous superdiffusion with a nontrivial dynamical exponent. Finally, we constitute the "Langevin rocket" model by including orientational fluctuations in the traditional Tsiolkovsky rocket equation. We calculate the mean reach of the Langevin rocket and discuss different mass ejection strategies to maximize it. Our results can be tested in experiments on macroscopic robotic or living particles or in self-propelled mesoscopic objects moving in media of low viscosity, such as complex plasma.