Motility-Induced Temperature Difference in Coexisting Phases

Spontaneous generation of angular momentum in chiral active crystals

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

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

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

Collective dynamics of intelligent active Brownian particles with visual perception and velocity alignment in 3D: spheres, rods, and worms

Many living systems, such as birds and fish, exhibit collective behaviors like flocking and swarming. Recently, an experimental system of active colloidal particles has been developed, where the motility of each particle is adjusted based on its visual detection of surrounding particles. These particles with visual-perception-dependent motility exhibit group formation and cohesion. Inspired by these behaviors, we investigate intelligent active Brownian particles (iABPs) equipped with visual perception and velocity alignment in three dimensions using computer simulations. The visual-perception-based self-steering describes the tendency of iABPs to move toward the center of mass of particles within their visual cones, while velocity alignment encourages alignment with neighboring particles. We examine how the behavior varies with the visual cone angle θ, self-propulsion speed (Péclet number Pe), and the interaction strengths of velocity alignment (Ωa) and visual-based self-steering (Ωv). Our findings show that spherical iABPs form dense clusters, worm-like clusters, milling behaviors, and dilute-gas phases, consistent with 2D studies. By reducing the simulation box size, we observe additional structures like band-like clusters and dense baitball formations. Additionally, rod-like iABPs form band-like, worm-like, radiating, and helical structures, while iABP worms exhibit band-like, streamlined, micellar-like and entangled structures. Many of these patterns resemble collective behaviors in nature, such as ant milling, fish baitballs, and worm clusters. Advances in synthetic techniques could enable nanorobots with similar capabilities, offering insights into multicellular systems through active matter.

Time-dependent propulsion of fully inertial active stochastic particles: theory and simulations

Up to now, studies on fully inertial (adding mass and moment of inertia) active Brownian particles (IABPs) have only considered a constant propulsion force. This work overcomes this by studying IABPs but with a time-dependent propulsion and analytically characterizes this system by finding its mean-square displacement and effective diffusion for any periodic time-dependent propulsion speed. To exemplify the periodic general expressions, three particular self-propulsion signals are addressed, expressly, a cosine, a square-wave, and a zig-zag propulsion force. Langevin dynamics simulations are also employed to validate the analytical findings.

Suppression of bubbles in unstable active liquids via fast evaporation

A common intuition in thermodynamics is that bubbles can spontaneously grow in unstable liquids, which will be detrimental to a variety of physical and chemical processes, such as evaporation-induced self-assembly and electrocatalysis. Here, we show that this common intuition can be significantly reversed by demonstrating a suppression of bubbles in unstable active liquids induced by fast evaporation, which is in contrast to the bubble growth in passive liquids. Such anomalous bubble suppression can be attributed to an activity-induced inversion of pressure difference between bubbles and their surrounding liquid. Moreover, this pressure flip depends on the activity as well as the thermodynamics of passive liquids, and it can generate different kinetic pathways that allow controlling the bubble dynamics in unstable liquids. Our results establish a foundation for promoting applications of unstable active liquids in various physical and chemical processes.

Harmonically trapped inertial run-and-tumble particle in one dimension

We study the nonequilibrium stationary state of a one-dimensional inertial run-and-tumble particle (IRTP) trapped in a harmonic potential. We find that the presence of inertia leads to two distinct dynamical scenarios, namely, overdamped and underdamped, characterized by the relative strength of the viscous and the trap timescales. We also find that inertial nature of the active dynamics leads to the particle being confined in specific regions of the phase plane in the overdamped and underdamped cases, which we compute analytically. Moreover, the interplay of the inertial and active timescales gives rise to several subregimes, which are characterized by very different behavior of position and velocity fluctuations of the IRTP. In particular, in the underdamped regime, both the position and velocity undergo transitions from a novel multipeaked structure in the strongly active limit to a single-peaked Gaussian-like distribution in the passive limit. On the other hand, in the overdamped scenario, the position distribution shows a transition from a U shape to a dome shape, as activity is decreased. Interestingly, the velocity distribution in the overdamped scenario shows two transitions-from a single-peaked shape with an algebraic divergence at the origin in the strongly active regime to a double-peaked one in the moderately active regime to a dome-shaped one in the passive regime.

Interpolating the radial distribution function in a two-dimensional fluid across a wide temperature range

Calculations of pair correlations in fluids usually require resource-intensive simulations or integral equations, while existing simple approximations lack accuracy. Here, we show that the pair correlation function for monolayer fluid-like systems can be decomposed into correlation peaks defined using Voronoi cells. Being properly normalized, these peaks exhibit a universal form, weak temperature dependence, and resemble those of an ideal gas, except for the first peak. As a result, we propose a simple and accurate approach to interpolate the pair correlation functions, suitable for molecular, colloids, and cellular fluids.

Density and inertia effects on two-dimensional active semiflexible filament suspensions

We examine the influence of density on the transition between chain and spiral structures in planar assemblies of active semiflexible filaments, utilizing detailed numerical simulations. We focus on how increased density, and higher Péclet numbers, affect the activity-induced transition spiral state in a semiflexible, self-avoiding active chain. Our findings show that increasing the density causes the spiral state to break up, reverting to a motile chain-like shape. This results in a density-dependent reentrant phase transition from spirals back to open chains. We attribute this phenomenon to an inertial effect observed at the single polymer level, where increased persistence length due to inertia has been shown in recent three-dimensional studies to cause polymers to open up. Our two-dimensional simulations further reveal that a reduction in the damping coefficient leads to partial unwinding of the spirals, forming longer arms. In suspension, interactions among these extended arms can trigger a complete unwinding of the spirals, driven by the combined effects of density and inertia.

Constrained motion of self-propelling eccentric disks linked by a spring

It has been supposed that the interplay of elasticity and activity plays a key role in triggering the non-equilibrium behaviors in biological systems. However, the experimental model system is missing to investigate the spatiotemporally dynamical phenomena. Here, a model system of an active chain, where active eccentric-disks are linked by a spring, is designed to study the interplay of activity, elasticity, and friction. Individual active chain exhibits longitudinal and transverse motions; however, it starts to self-rotate when pinning one end and self-beat when clamping one end. In addition, our eccentric-disk model can qualitatively reproduce such behaviors and explain the unusual self-rotation of the first disk around its geometric center. Furthermore, the structure and dynamics of long chains were studied via simulations without steric interactions. It was found that a hairpin conformation emerges in free motion, while in the constrained motions, the rotational and beating frequencies scale with the flexure number (the ratio of self-propelling force to bending rigidity), χ, as ∼(χ)4/3. Scaling analysis suggests that it results from the balance between activity and energy dissipation. Our findings show that topological constraints play a vital role in non-equilibrium synergy behaviors.

Inertia and activity: spiral transitions in semi-flexible, self-avoiding polymers

We consider a two-dimensional, tangentially active, semi-flexible, self-avoiding polymer to find a dynamical re-entrant transition between motile open chains and spinning achiral spirals with increasing activity. Utilizing probability distributions of the turning number, we ascertain the comparative stability of the spiral structure and present a detailed phase diagram within the activity inertia plane. The onset of spiral formation at low activity levels is governed by a torque balance and is independent of inertia. At higher activities, however, inertial effects lead to spiral destabilization, an effect absent in the overdamped limit. We further delineate alterations in size and shape by analyzing the end-to-end distance distribution and the radius of gyration tensor. The Kullback-Leibler divergence from equilibrium distributions exhibits a non-monotonic relationship with activity, reaching a peak at the most compact spirals characterized by the most persistent spinning. As inertia increases, this divergence from equilibrium diminishes.

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

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

Fluctuation Theorems for Heat Exchanges between Passive and Active Baths

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

Inertial active harmonic particle with memory induced spreading by viscoelastic suspension

We investigate the self-propulsion of an inertial active particle confined in a two-dimensional harmonic trap. The particle is suspended in a non-Newtonian or viscoelastic suspension with a friction kernel that decays exponentially with a time constant characterizing the memory timescale or transient elasticity of the medium. By solving the associated non-Markovian dynamics, we identify two regimes in parameter space distinguishing the oscillatory and non-oscillatory behavior of the particle motion. By simulating the particle trajectories and exactly calculating the steady-state probability distribution functions and mean square displacement; interestingly, we observe that with an increase in the memory time scale, the effective temperature of the environment increases. As a consequence, the particle becomes energetic and spread away from the center, covering larger space inside the confinement. On the other hand, with an increase in the duration of the activity, the particle becomes trapped by the harmonic confinement.

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.

Influence of anomalous agents on the dynamics of an active system

Swarming behavior in systems of self-propelled particles, whether biological or artificial, has received increased attention in recent years. Here, we show that even a small number of particles with anomalous behavior can change dramatically collective dynamics of the swarming system and can impose unusual behavior and transitions between dynamic states. Our results pave the way to practical approaches and concepts of multiagent dynamics in groups of flocking animals: birds, insects, and fish, i.e., active and living soft matter.

Force renormalization for probes immersed in an active bath

Langevin equations or generalized Langevin equations (GLEs) are popular models for describing the motion of a particle in a fluid medium in an effective manner. Here we examine particles immersed in an inherently nonequilibrium fluid, i.e., an active bath, which are subject to an external force. Specifically, we consider two types of forces that are highly relevant for microrheological studies: A harmonic, trapping force and a constant, "drag" force. We study such systems by molecular simulations and use the simulation data to extract an effective GLE description. We find that within this description, in an active bath, the external force in the GLE is not equal to the physical external force, but rather a renormalized external force, which can be significantly smaller. The effect cannot be attributed to the mere temperature renormalization, which is also observed.

When an active bath behaves as an equilibrium one

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

Structure of the active Fokker-Planck equation

This paper solves in one and two dimensions the steady noninteractive active Fokker-Planck (FP) equation and finds that its velocity distribution admits, under limiting cases, a dual behavior. Briefly, when the inertial relaxation time is smaller than the orientation time, the active FP equation admits a bimodal shape, whereas the inverse condition is seen to admit a Gaussian one. Once the velocity distribution functions are available, they are used to find their effect on the system's transport properties, such as its mean-square speed. In the process, a useful mathematical identity for the first kind Bessel function as a sum of bimodal exponential functions is spotted.

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.

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

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

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.

Dynamics of active particles with translational and rotational inertia

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

Flocking without Alignment Interactions in Attractive Active Brownian Particles

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

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

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

From a microscopic inertial active matter model to the Schrödinger equation

Active field theories, such as the paradigmatic model known as 'active model B+', are simple yet very powerful tools for describing phenomena such as motility-induced phase separation. No comparable theory has been derived yet for the underdamped case. In this work, we introduce active model I+, an extension of active model B+ to particles with inertia. The governing equations of active model I+ are systematically derived from the microscopic Langevin equations. We show that, for underdamped active particles, thermodynamic and mechanical definitions of the velocity field no longer coincide and that the density-dependent swimming speed plays the role of an effective viscosity. Moreover, active model I+ contains an analog of the Schrödinger equation in Madelung form as a limiting case, allowing one to find analoga of the quantum-mechanical tunnel effect and of fuzzy dark matter in active fluids. We investigate the active tunnel effect analytically and via numerical continuation.

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.

Cooperation and competition in the collective drive by motor proteins: mean active force, fluctuations, and self-load

We consider the dynamics of a bio-filament under the collective drive of motor proteins. They are attached irreversibly to a substrate and undergo stochastic attachment-detachment with the filament to produce a directed force on it. We establish the dependence of the mean directed force and force correlations on the parameters describing the individual motor proteins using analytical theory and direct numerical simulations. The effective Langevin description for the filament motion gives mean-squared displacement, asymptotic diffusion constant, and mobility leading to an effective temperature. Finally, we show how competition between motor protein extensions generates a self-load, describable in terms of the effective temperature, affecting the filament motion.

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.

Two-temperature activity induces liquid-crystal phases inaccessible in equilibrium

In equilibrium hard-rod fluids, and in effective hard-rod descriptions of anisotropic soft-particle systems, the transition from the isotropic (I) phase to the nematic phase (N) is observed above the rod aspect ratio L/D=3.70 as predicted by Onsager. We examine the fate of this criterion in a molecular dynamics study of a system of soft repulsive spherocylinders rendered active by coupling half the particles to a heat bath at a higher temperature than that imposed on the other half. We show that the system phase-separates and self-organizes into various liquid-crystalline phases that are not observed in equilibrium for the respective aspect ratios. In particular, we find a nematic phase for L/D=3 and a smectic phase for L/D=2 above a critical activity.

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.

Collective behavior of soft self-propelled disks with rotational inertia

We investigate collective properties of a large system of soft self-propelled inertial disks with active Langevin dynamics simulation in two dimensions. Rotational inertia of the disks is found to favor motility induced phase separation (MIPS), due to increased effective persistence of the disks. The MIPS phase diagram in the parameter space of rotational inertia and disk softness is reported over a range of values of translation inertia and self-propulsion strength of the disks. Our analytical prediction of the phase boundary between the homogeneous (no-MIPS) and MIPS state in the limit of small and large rotational inertia is found to agree with the numerical data over a large range of translational inertia. Shape of the high density MIPS phase is found to change from circular to rectangular one as the system moves away from the phase boundary. Structural and dynamical properties of the system, measured by several physical quantities, are found to be invariant in the central region of the high density MIPS phase, whereas they are found to vary gradually near the peripheral region of the high density phase. Importantly, the width of the peripheral region near the phase boundary is much larger compared to the narrow peripheral region far away from the phase boundary. Rich dynamics of the disks inside the high density MIPS phase is addressed. Spatial correlation of velocity of the disks is found to increase with rotational inertia and disk hardness. However, temporal correlation of the disks' velocity is found to be a function of rotational inertia, while it is independent of disk softness.

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.

Motility-induced phase separation of self-propelled soft inertial disks

The phase diagram of the phenomenon of motility-induced phase separation (MIPS) for a collection of self-propelled interacting disks over a large inertial range is explored using active Langevin dynamics simulation with particular emphasis on disk softness and effective size. It is shown that the parabola-like phase boundary between the homogeneous and MIPS states in the semi-log space of disk softness and effective size moves towards the hard disk limit with increase in inertia, before complete disappearance in the limit of large inertia. With increase in effective size of the disks, re-entrant phase separation, that is the system phase-separating from a homogeneous phase and eventually re-entering the homogeneous phase, is reported. The structural and the dynamical properties of the different phases are investigated in the considered inertial range. The particular shape of the phase boundary and the re-entrant behavior is explained based on several qualitative and quantitative results. Unlike most of the earlier studies on MIPS, which consider hard particle limits, our findings may be directly applicable to soft active matter for a range of physical and biological systems.

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.

Passive probe particle in an active bath: can we tell it is out of equilibrium?

We study a passive probe immersed in a fluid of active particles. Despite the system's non-equilibrium nature, the trajectory of the probe does not exhibit non-equilibrium signatures: its velocity distribution remains Gaussian, the second fluctuation dissipation theorem is not fundamentally violated, and the motion does not indicate breaking of time reversal symmetry. To tell that the probe is out of equilibrium requires examination of its behavior in tandem with that of the active fluid: the kinetic temperature of the probe does not equilibrate to that of the surrounding active particles. As a strategy to diagnose non-equilibrium from probe trajectories alone, we propose to examine their response to a small perturbation which reveals a non-equilibrium signature through a violation of the first fluctuation dissipation theorem.

Percolation transition in phase-separating active fluid

The motility-induced phase separation exhibited by active particles with repulsive interactions is well known. We show that the interaction softness of active particles destabilizes the highly ordered dense phase, leading to the formation of a porous cluster which spans the system. This soft limit can also be achieved if the particle motility is increased beyond a critical value, at which the system clearly exhibits all the characteristics of a standard percolation transition. We also show that in the athermal limit, active particles exhibit similar transitions even at low motility. With these additional new phases, the phase diagram of repulsive active particles is revealed to be richer than what was previously conceived.

Emergent collective behavior of active Brownian particles with visual perception

Systems comprised of self-steering active Brownian particles are studied via simulations for a minimal cognitive flocking model. The dynamics of the active Brownian particles is extended by an orientational response with limited maneuverability to an instantaneous visual input of the positions of neighbors within a vision cone and a cut-off radius. The system exhibits large-scale self-organized structures, which depend on selected parameter values, and, in particular, the presence of excluded-volume interactions. The emergent structures in two dimensions, such as worms, worm-aggregate coexistence, and hexagonally close-packed structures, are analysed and phase diagrams are constructed. The analysis of the particle's mean-square displacement shows ABP-like dynamics for dilute systems and the worm phase. In the limit of densely packed structures, the active diffusion coefficient is significantly smaller and depends on the number of particles in the cluster. Our analysis of the cluster-growth dynamics shows distinct differences to processes in systems of short-range attractive colloids in equilibrium. Specifically, the characteristic time for the growth and decay of clusters of a particular size is longer than that of isotropically attractive colloids, which we attribute to the non-reciprocal nature of the directed visual perception. Our simulations reveal a strong interplay between ABP-characteristic interactions, such as volume exclusion and rotational diffusion, and cognitive-based interactions and navigation.

Self-sustained non-equilibrium co-existence of fluid and solid states in a strongly coupled complex plasma system

A complex (dusty) plasma system is well known as a paradigmatic model for studying the kinetics of solid-liquid phase transitions in inactive condensed matter. At the same time, under certain conditions a complex plasma system can also display characteristics of an active medium with the micron-sized particles converting energy of the ambient environment into motility and thereby becoming active. We present a detailed analysis of the experimental complex plasmas system that shows evidence of a non-equilibrium stationary coexistence between a cold crystalline and a hot fluid state in the structure due to the conversion of plasma energy into the motion energy of microparticles in the central region of the system. The plasma mediated non-reciprocal interaction between the dust particles is the underlying mechanism for the enormous heating of the central subsystem, and it acts as a micro-scale energy source that keeps the central subsystem in the molten state. Accurate multiscale simulations of the system based on combined molecular dynamics and particle-in-cell approaches show that strong structural nonuniformity of the system under the action of electostatic trap makes development of instabilities a local process. We present both experimental tests conducted with a complex plasmas system in a DC glow discharge plasma and a detailed theoretical analysis.

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.

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

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

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.

Heating leads to liquid-crystal and crystalline order in a two-temperature active fluid of rods

We report phase separation and liquid-crystal ordering induced by scalar activity in a system of soft repulsive spherocylinders (SRSs) of shape anisotropy L/D=5 using molecular dynamics (MD) simulations. Activity is introduced by increasing the temperature of half of the SRSs (labeled hot) while maintaining the temperature of the other half constant at a lower value (labeled cold). The difference between the two temperatures scaled by the lower temperature provides a measure of the activity. Starting from different equilibrium initial phases, we find that activity leads to segregation of the hot and cold particles. Activity also drives the cold particles through a phase transition to a more ordered state and the hot particles to a state of less order compared to the initial equilibrium state. The cold components of a homogeneous isotropic structure acquire nematic and, at higher activity, crystalline order. Similarly, the cold zone of a nematic initial state undergoes smectic and crystal ordering above a critical value of activity while the hot component turns isotropic. We find that the hot particles occupy a larger volume and exert an extra kinetic pressure, confining, compressing, and provoking an ordering transition of the cold-particle domains.

Anomalous transport phenomenon of a charged Brownian particle under a thermal gradient and a magnetic field

There is a growing interest in the stochastic processes of nonequilibrium systems subject to nonconserved forces, such as the magnetic forces acting on charged particles and the chiral self-propelled force acting on active particles. In this paper, we consider the stationary transport of noninteracting Brownian particles under a constant magnetic field in a position-dependent temperature background. We demonstrate the existence of the Nernst-like stationary density current perpendicular to both the temperature gradient and magnetic field, induced by the intricate coupling between the nonconserved force and the multiplicative noises due to the position-dependent temperature.

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.