Ideal bulk pressure of active Brownian particles

Motility-Induced Clustering of Active Particles under Soft Confinement

We investigate the structural and dynamic properties of active Brownian particles (APs) confined within a soft annulus-shaped channel. Depending on the strength of the confinement and the Péclet number, we observe a novel reentrant behavior that is not present in unconfined systems. Our findings are substantiated by numerical simulations and analytical considerations, revealing that this behavior arises from the strong coupling between the Péclet number and the effective confining dimensionality of the APs. Our work highlights the peculiarities of soft boundaries for APs and how clogging can be avoided under such conditions.

Partial and complete wetting of droplets of active Brownian particles

We study wetting droplets formed of active Brownian particles in contact with a repulsive potential barrier, in a wedge geometry. Our numerical results demonstrate a transition between partially wet and completely wet states, as a function of the barrier height, analogous to the corresponding surface phase transition in passive fluids. We analyse partially wet configurations characterised by a nonzero contact angle θ between the droplet surface and the barrier including the average density profile and its fluctuations. These findings are compared with two equilibrium systems: a Lennard-Jones fluid and a simple contour model for a liquid-vapour interface. We locate the wetting transition where cos(θ) = 1, and the neutral state where cos(θ) = 0. We discuss the implications of these results for possible definitions of surface tensions in active fluids.

Kinetic temperature and pressure of an active Tonks gas

Using computer simulation and analytical theory, we study an active analog of the well-known Tonks gas, where active Brownian particles are confined to a periodic one-dimensional (1D) channel. By introducing the notion of a kinetic temperature, we derive an accurate analytical expression for the pressure and clarify the paradoxical behavior where active Brownian particles confined to 1D exhibit anomalous clustering but no motility-induced phase transition. More generally, this work provides a deeper understanding of pressure in active systems as we uncover a unique link between the kinetic temperature and swim pressure valid for active Brownian particles in higher dimensions.

Effect of confinement and topology: 2-TIPS vs. MIPS

2-TIPS (two temperature induced phase separation) refers to the phase separation phenomenon observed in mixtures of active and passive particles which are modelled using scalar activity. The active particles are connected to a thermostat at high temperature while the passive particles are connected to the thermostat at low temperature and the relative temperature difference between "hot" and "cold" particles is taken as the measure of the activity χ of the non-equilibrium system. The study of such binary mixtures of hot and cold particles under various kinds of confinement is an important problem in many physical and biological processes. The nature and extent of phase separation are heavily influenced by the geometry of confinement, activity, and density of the non-equilibrium binary mixture. Investigating such 3D binary mixtures confined by parallel walls, we observe that the active and passive particles phase separate, but the extent of phase separation is reduced compared to bulk phase separation at high densities and enhanced at low densities. However, when the binary mixture of active and passive particles is confined inside a spherical cavity, the phase separation is radial for small radii of the confining sphere and the extent of phase separation is higher compared to their bulk counterparts. Confinement leads to interesting properties in the passive (cold) region like enhanced layering and high compression in the direction parallel to the confining wall. In 2D, both the bulk and confined systems of the binary mixture show a significant decrement in the extent of phase separation at higher densities. This observation is attributed to the trapping of active particles inside the passive cluster, which increases with density. Thus the 2D systems show structures more akin to dense-dilute phase co-existence, which is observed in motility induced phase separation in 2D active systems. The binary mixture constrained on the spherical surface also shows similar phase co-existence. Our analyses reveal that the coexistent densities observed in 2-TIPS on the spherical surface agree with the findings of previous studies on MIPS in active systems on a sphere.

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.

Intrinsic structure perspective for MIPS interfaces in two-dimensional systems of active Brownian particles

Suspensions of active Brownian particles (ABPs) undergo motility-induced phase separation (MIPS) over a wide range of mean density and activity strength, which implies the spontaneous aggregation of particles due to the persistence of their direction of motion, even in the absence of an explicit attraction. Both similarities and qualitative differences have been obtained when the MIPS is analysed in the same terms as a liquid-gas phase coexistence in an equilibrium attractive system. Negative values of the mechanical surface tension have been reported, from the total forces across the interface, while stable fluctuations of the interfacial line could be interpreted as a positive capillary surface tension; in equilibrium liquid surfaces, these two magnitudes are equal. We present here the analysis of 2D-ABP interfaces in terms of the intrinsic density and force profiles, calculated with the particle distance to the instantaneous interfacial line. Our results provide new insight into the origin of MIPS from the local rectification of the random active force on the particles near the interface. As has been reported, this effect acts as an external potential that produces a pressure gradient across the interface, such that the mechanical surface tension of the MIPS cannot be described as that of equilibrium coexisting phases; however, our analysis shows that most of that effect comes from the tightly caged particles at the dense (inner) side of the MIPS interface, rather than from the free moving particles at the outer side that collide with the dense cluster. Moreover, a clear correlation appears between the decay of the hexatic order parameter at the dense slab and the end of the MIPS as the strength of the active force is lowered. We show that, using the strong active forces required for MIPS, the interfacial structure and properties are very similar for ABPs with purely repulsive interactions (the Weeks-Chandler-Andersen-Lennard-Jones (WCA-LJ) model truncated at its minimum) and when the interaction includes a range of the LJ attractive force.

Dynamic overlap concentration scale of active colloids

By introducing the notion of a dynamic overlap concentration scale, we identify additional universal features of the mechanical properties of active colloids. We codify these features by recognizing that the characteristic length scale of an active particle's trajectory, the run length, introduces a concentration scale ϕ^{*}. Large-scale simulations of repulsive active Brownian particles (ABPs) confirm that this run-length dependent concentration, the trajectory-space analog of the overlap concentration in polymer solutions, delineates distinct concentration regimes in which interparticle collisions alter particle trajectories. Using ϕ^{*} and concentration scales associated with colloidal jamming, the mechanical equation of state for ABPs collapses onto a set of principal curves that contain several overlooked features. The inclusion of these features qualitatively alters previous predictions of the behavior for active colloids, as we demonstrate by computing the spinodal for a suspension of purely repulsive ABPs. Our findings suggest that dynamic overlap concentration scales should help unravel the behavior of active and driven systems.

Dependency of active pressure and equation of state on stiffness of wall

Autonomous motion and motility are hallmarks of active matter. Active agents, such as biological cells and synthetic colloidal particles, consume internal energy or extract energy from the environment to generate self-propulsion and locomotion. These systems are persistently out of equilibrium due to continuous energy consumption. It is known that pressure is not always a state function for generic active matter. Torque interaction between active constituents and confinement renders the pressure of the system a boundary-dependent property. The mechanical pressure of anisotropic active particles depends on their microscopic interactions with a solid wall. Using self-propelled dumbbells confined by solid walls as a model system, we perform numerical simulations to explore how variations in the wall stiffness influence the mechanical pressure of dry active matter. In contrast to previous findings, we find that mechanical pressure can be independent of the interaction of anisotropic active particles with walls, even in the presence of intrinsic torque interaction. Particularly, the dependency of pressure on the wall stiffness vanishes when the stiffness is above a critical level. In such a limit, the dynamics of dumbbells near the walls are randomized due to the large torque experienced by the dumbbells, leading to the recovery of pressure as a state variable of density.

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.

Polarization-density patterns of active particles in motility gradients

The colocalization of density modulations and particle polarization is a characteristic emergent feature of motile active matter in activity gradients. We employ the active-Brownian-particle model to derive precise analytical expressions for the density and polarization profiles of a single Janus-type swimmer in the vicinity of an abrupt activity step. Our analysis allows for an optional (but not necessary) orientation-dependent propulsion speed, as often employed in force-free particle steering. The results agree well with measurement data for a thermophoretic microswimmer presented in the companion paper [Söker et al., Phys. Rev. Lett. 126, 228001 (2021)10.1103/PhysRevLett.126.228001], and they can serve as a template for more complex applications, e.g., to motility-induced phase separation or studies of physical boundaries. The essential physics behind our formal results is robustly captured and elucidated by a schematic two-species "run-and-tumble" model.

Vorticity Determines the Force on Bodies Immersed in Active Fluids

When immersed into a fluid of active Brownian particles, passive bodies might start to undergo linear or angular directed motion depending on their shape. Here we exploit the divergence theorem to relate the forces responsible for this motion to the density and current induced by-but far away from-the body. In general, the force is composed of two contributions: due to the strength of the dipolar field component and due to particles leaving the boundary, generating a nonvanishing vorticity of the polarization. We derive and numerically corroborate results for periodic systems, which are fundamentally different from unbounded systems with forces that scale with the area of the system. We demonstrate that vorticity is localized close to the body and to points at which the local curvature changes, enabling the rational design of particle shapes with desired propulsion properties.

Active particles with polar alignment in ring-shaped confinement

We study steady-state properties of active, nonchiral and chiral Brownian particles with polar alignment and steric interactions confined within a ring-shaped confinement (annulus) in two dimensions. Exploring possible interplays between polar interparticle alignment, geometric confinement and the surface curvature, being incorporated here on minimal levels, we report a surface-population reversal effect, whereby active particles migrate from the outer concave boundary of the annulus to accumulate on its inner convex boundary. This contrasts the conventional picture, implying stronger accumulation of active particles on concave boundaries relative to the convex ones. The population reversal is caused by both particle alignment and surface curvature, disappearing when either of these factors is absent. We explore the ensuing consequences for the chirality-induced current and swim pressure of active particles and analyze possible roles of system parameters, such as the mean number density of particles and particle self-propulsion, chirality, and alignment strengths.

Coexistence of active Brownian disks: van der Waals theory and analytical results

At thermal equilibrium, intensive quantities like temperature and pressure have to be uniform throughout the system, restricting inhomogeneous systems composed of different phases. The paradigmatic example is the coexistence of vapor and liquid, a state that can also be observed for active Brownian particles steadily driven away from equilibrium. Recently, a strategy has been proposed that allows to predict phase equilibria of active particles [Solon et al., Phys. Rev. E 97, 020602(R) (2018)2470-004510.1103/PhysRevE.97.020602]. Here we elaborate on this strategy and formulate it in the framework of a van der Waals theory for active disks. For a given equation of state, we derive the effective free energy analytically and show that it yields coexisting densities in very good agreement with numerical results. We discuss the interfacial tension and the relation to Cahn-Hilliard models.

Nonequilibrium grand-canonical ensemble built from a physical particle reservoir

We introduce a nonequilibrium grand-canonical ensemble defined by considering the stationary state of a driven system of particles put in contact with a particle reservoir. When an additivity assumption holds for the large deviation function of density, a chemical potential of the reservoir can be defined. The grand-canonical distribution then takes a form similar to the equilibrium one. At variance with equilibrium, though, the probability weight is "renormalized" by a contribution coming from the contact, with respect to the canonical probability weight of the isolated system. A formal grand-canonical potential can be introduced in terms of a scaled cumulant generating function, defined as the Legendre-Fenchel transform of the large deviation function of density. The role of the formal Legendre parameter can be played, physically, by the chemical potential of the reservoir when the latter can be defined, or by a potential energy difference applied between the system and the reservoir. Static fluctuation-response relations naturally follow from the large deviation structure. Some of the results are illustrated on two different explicit examples, a gas of noninteracting active particles and a lattice model of interacting particles.

Chemically Propelled Nano and Micromotors in the Body: Quo Vadis?

The active delivery of drugs to disease sites in response to specific biomarkers is a holy grail in theranostics. If successful, it would greatly diminish the therapeutic dosage and reduce collateral cytotoxicity. In this context, the development of nano and micromotors that are able to harvest local energy to move directionally is an important breakthrough. However, serious hurdles remain before such active systems can be employed in vivo in therapeutic applications. Such motors and their energy sources must be safe and biocompatible, they should be able to move through complex body fluids, and have the ability to reach specific cellular targets. Given the complexity in the design and deployment of nano and micromotors, it is also critically important to show that they are significantly superior to inactive "smart" nanoparticles in theranostics. Furthermore, receiving regulatory approval requires the ability to scale-up the production of nano and micromotors with uniformity in structure, function, and activity. In this essay, the limitations of the current nano and micromotors and the issues that need to be resolved before such motors are likely to find theranostic applications are discussed.

Aggregate morphology of active Brownian particles on porous, circular walls

We study the motility-induced aggregation of active Brownian particles (ABPs) on a porous, circular wall. We observe that the morphology of aggregated dense-phase on a static wall depends on the wall porosity, particle motility, and the radius of the circular wall. Our analysis reveals two morphologically distinct, dense aggregates; a connected dense cluster that spreads uniformly on the circular wall and a localized cluster that breaks the rotational symmetry of the system. These distinct morphological states are similar to the macroscopic structures observed in aggregates on planar, porous walls. We systematically analyze the parameter regimes where the different morphological states are observed. We further extend our analysis to motile circular rings. We show that the motile ring propels almost ballistically due to the force applied by the active particles when they form a localized cluster, whereas it moves diffusively when the active particles form a continuous cluster. This property demonstrates the possibility of extracting useful work from a system of ABPs, even without artificially breaking the rotational symmetry.

Morphological transitions of active Brownian particle aggregates on porous walls

Motility-induced wall aggregation of Active Brownian Particles (ABPs) is a well-studied phenomenon. Here, we study the aggregation of ABPs on porous walls, which allows the particles to penetrate through at large motility. We show that the active aggregates undergo a morphological transition from a connected dense-phase to disconnected droplets with an increase in wall porosity and the particle self-motility, similar to wetting-dewetting transitions in equilibrium fluids. We show that both morphologically distinct states are stable, and independent of initial conditions at least in some parameter regions. Our analysis reveals that changes in wall porosity affect the intrinsic properties of the aggregates and changes the effective wall-aggregate interfacial tension, consistent with the appearance of the morphological transition. Accordingly, a close analysis of the density, as well as orientational distribution, indicates that the underlying reason for such morphological transitions is not necessarily specific to the systems with porous walls, and it can be possible to observe in a larger class of confined, active systems by tuning the properties of confining walls.

Surface Tensions between Active Fluids and Solid Interfaces: Bare vs Dressed

We analyze the surface tension exerted at the interface between an active fluid and a solid boundary in terms of tangential forces. Focusing on active systems known to possess an equation of state for the pressure, we show that interfacial forces are of a more complex nature. Using a number of macroscopic setups, we show that the surface tension is a combination of an equation-of-state abiding part and of setup-dependent contributions. The latter arise from generic setup-dependent steady currents which "dress" the measurement of the "bare" surface tension. The former shares interesting properties with its equilibrium counterpart, and can be used to generalize the Young-Laplace law to active systems. Finally, we show how a suitably designed probe can directly access this bare surface tension, which can also be computed using a generalized virial formula.

Active matter therapeutics

Nanotherapies based on micelles, liposomes, polymersomes, nanocapsules, magnetic nanoparticles, and noble metal nanoparticles have been at the forefront of drug delivery in the past few decades. Some of these nanopharmaceuticals have been commercially applied to treat a wide range of diseases, from dry eye syndrome to cancer. However, the majority involve particles that are passive, meaning that they do not change shape, and they lack motility; the static features can limit their therapeutic efficacy. In this review, we take a critical look at an emerging field that seeks to utilize active matter for therapeutics. In this context, active matter can be broadly referred to as micro or nanosized constructs that energetically react with their environment or external fields and translate, rotate, vibrate or change shape. Essentially, the recent literature suggests that such particles could significantly augment present-day drug delivery, by enhancing transport and increasing permeability across anatomical barriers by transporting drugs within solid tumor microenvironments or disrupting cardiovascular plaque. We discuss examples of such particles and link the transport and permeability properties of active matter to potential therapeutic applications in the context of two major diseases, namely cancer and heart disease. We also discuss potential challenges, opportunities, and translational hurdles.

Collective forces in scalar active matter

Large-scale collective behavior in suspensions of active particles can be understood from the balance of statistical forces emerging beyond the direct microscopic particle interactions. Here we review some aspects of the collective forces that can arise in suspensions of self-propelled active Brownian particles: wall forces under confinement, interfacial forces, and forces on immersed bodies mediated by the suspension. Even for non-aligning active particles, these forces are intimately related to a non-uniform polarization of particle orientations induced by walls and bodies, or inhomogeneous density profiles. We conclude by pointing out future directions and promising areas for the application of collective forces in synthetic active matter, as well as their role in living active matter.

Pressure and diffusion of active matter with inertia

It has been discovered that active matter generates novel physical quantities such as the swim pressure. This quantity arises from the exchange of extra momentum between active particles and the boundaries of the system. Given its origin, this quantity can exist at different scales; hence microorganisms and larger organisms like fish or birds generate their own swim pressure. For larger organisms or for high swimming speeds, inertia cannot necessarily be neglected; hence in this paper, we start by calculating analytically the effect of finite translational and rotational particles' inertia on the diffusion of a system of noninteracting spherical active Brownian particles. From this analysis, an enhanced diffusion coefficient due to rotational inertia is obtained, and an alternative effective persistence length and an alternative reorientation time, both sensitive to rotational inertia, are also identified. Afterwards, and to see the implications of finite inertia on bulk properties, the pressure of this system is elucidated by calculating its respective swim and Reynolds pressures. It is found that their sum becomes asymptotically sensitive to the square root of its rotational inertia. To validate our analytical results, Langevin dynamics simulations are also performed showing an excellent agreement between our theoretical predictions and the numerical results.

Microscopic origins of the swim pressure and the anomalous surface tension of active matter

The unique pressure exerted by active particles-the "swim" pressure-has proven to be a useful quantity in explaining many of the seemingly confounding behaviors of active particles. However, its use has also resulted in some puzzling findings including an extremely negative surface tension between phase separated active particles. Here, we demonstrate that this contradiction stems from the fact that the swim pressure is not a true pressure. At a boundary or interface, the reduction in particle swimming generates a net active force density-an entirely self-generated body force. The pressure at the boundary, which was previously identified as the swim pressure, is in fact an elevated (relative to the bulk) value of the traditional particle pressure that is generated by this interfacial force density. Recognizing this unique mechanism for stress generation allows us to define a much more physically plausible surface tension. We clarify the utility of the swim pressure as an "equivalent pressure" (analogous to those defined from electrostatic and gravitational body forces) and the conditions in which this concept can be appropriately applied.

Local stress and pressure in an inhomogeneous system of spherical active Brownian particles

The stress of a fluid on a confining wall is given by the mechanical wall forces, independent of the nature of the fluid being passive or active. At thermal equilibrium, an equation of state exists and stress is likewise obtained from intrinsic bulk properties; even more, stress can be calculated locally. Comparable local descriptions for active systems require a particular consideration of active forces. Here, we derive expressions for the stress exerted on a local volume of a systems of spherical active Brownian particles (ABPs). Using the virial theorem, we obtain two identical stress expressions, a stress due to momentum flux across a hypothetical plane, and a bulk stress inside of the local volume. In the first case, we obtain an active contribution to momentum transport in analogy to momentum transport in an underdamped passive system, and we introduce an active momentum. In the second case, a generally valid expression for the swim stress is derived. By simulations, we demonstrate that the local bulk stress is identical to the wall stress of a confined system for both, non-interacting ABPs as well as ABPs with excluded-volume interactions. This underlines the existence of an equation of state for a system of spherical ABPs. Most importantly, our calculations demonstrated that active stress is not a wall (boundary) effect, but is caused by momentum transport. We demonstrate that the derived stress expression permits the calculation of the local stress in inhomogeneous systems of ABPs.

Microscopic analysis of thermo-orientation in systems of off-centre Lennard-Jones particles

When fluids of anisotropic molecules are placed in temperature gradients, the molecules may align themselves along the gradient: this is called thermo-orientation. We discuss the theory of this effect in a fluid of particles that interact by a spherically symmetric potential, where the particles' centres of mass do not coincide with their interaction centres. Starting from the equations of motion of the molecules, we show how a simple assumption of local equipartition of energy can be used to predict the thermo-orientation effect, recovering the result of Wirnsberger et al. [Phys. Rev. Lett. 120, 226001 (2018)]. Within this approach, we show that for particles with a single interaction centre, the thermal centre of the molecule must coincide with the interaction centre. The theory also explains the coupling between orientation and kinetic energy that is associated with this non-Boltzmann distribution. We discuss deviations from this local equipartition assumption, showing that these can occur in linear response to a temperature gradient. We also present numerical simulations showing significant deviations from the local equipartition predictions, which increase as the centre of mass of the molecule is displaced further from its interaction centre.

Hydrodynamic interaction induced breakdown of the state properties of active fluids

The mechanical pressure of active fluids in which swimmers are modeled by soft run-and-tumble spheres is investigated by dissipative particle dynamics simulations. The incremental pressure (Π) with respect to the system pressure with inactive swimmers comprises the direct contribution of the swimmers (π) and the indirect contribution of fluids associated with hydrodynamic interactions (HIs). The pressure can be determined from the bulk and confining wall and the former is always less than the latter. The π of dilute active dispersions is proportional to their active diffusivity while Π grows generally with propulsive force and run time. However, Π is always substantially less than π because of negative contributions to pressure by HIs. The wall pressure depends on the swimmer-wall interactions, verifying that pressure is not a state function for active spheres due to the HIs. Owing to the distinct flow patterns, Π varies with the swim-type (pusher and puller) subject to the same run-and-tumble parameters at high concentrations.

Active fluids at circular boundaries: swim pressure and anomalous droplet ripening

We investigate the swim pressure exerted by non-chiral and chiral active particles on convex or concave circular boundaries. Active particles are modeled as non-interacting and non-aligning self-propelled Brownian particles. The convex and concave circular boundaries are used to model a fixed inclusion immersed in an active bath and a cavity (or container) enclosing the active particles, respectively. We first present a detailed analysis of the role of convex versus concave boundary curvature and of the chirality of active particles in their spatial distribution, chirality-induced currents, and the swim pressure they exert on the bounding surfaces. The results will then be used to predict the mechanical equilibria of suspended fluid enclosures (generically referred to as 'droplets') in a bulk with active particles being present either inside the bulk fluid or within the suspended droplets. We show that, while droplets containing active particles behave in accordance with standard capillary paradigms when suspended in a normal bulk, those containing a normal fluid exhibit anomalous behaviors when suspended in an active bulk. In the latter case, the excess swim pressure results in non-monotonic dependence of the inside droplet pressure on the droplet radius; hence, revealing an anomalous regime of behavior beyond a threshold radius, in which the inside droplet pressure increases upon increasing the droplet size. Furthermore, for two interconnected droplets, mechanical equilibrium can occur also when the droplets have different sizes. We thus identify a regime of anomalous droplet ripening, where two unequal-sized droplets can reach a final state of equal size upon interconnection, in stark contrast with the standard Ostwald ripening phenomenon, implying shrinkage of the smaller droplet in favor of the larger one.

Mechanical pressure, surface excess, and polar order of a dilute rod-like nanoswimmer suspension: role of swimmer-wall interactions

The mechanical pressure, surface excess, and polar order of a dilute rod-like nanoswimmer suspension confined by two parallel plates are explored by dissipative particle dynamics. The accumulation and preferred orientation of swimmers near the walls are distinctly shown through the density and polar order distributions for various active force, Fa, values and rod lengths. As Fa is increased, it is interesting to observe that there exists a maximum of the polar order, revealing that the dominant mechanism of the swimmer behavior can be altered by the coupling between the active force and the rod-wall interaction. As a result, the influences of the active force on the swim pressure Π(w)a contributed by the swimmers directly and the surface excess Γ* can be classified into two scaling regimes, natural rotation (weak propulsion) and forced rotation (strong propulsion). Π(w)a and Γ* are proportional to Fa2 in the former regime but become proportional to Fa in the latter regime. For all rod-wall repulsions, the swim pressure of active rods in confined systems Π(w)a always differs from that in unbounded systems Π(b)a which is simply proportional to Fa2 associated with the active diffusivity. That is, unlike thermal equilibrium systems, Π(w)a is not a state function because of the presence of the wall-torque.

Active ideal sedimentation: exact two-dimensional steady states

We consider an ideal gas of active Brownian particles that undergo self-propelled motion and both translational and rotational diffusion under the influence of gravity. We solve analytically the corresponding Smoluchowski equation in two space dimensions for steady states. The resulting one-body density is given as a series, where each term is a product of an orientation-dependent Mathieu function and a height-dependent exponential. A lower hard wall is implemented as a no-flux boundary condition. Numerical evaluation of the suitably truncated analytical solution shows the formation of two different spatial regimes upon increasing Peclet number. These regimes differ in their mean particle orientation and in their variation of the orientation-averaged density with height.

Generic Long-Range Interactions Between Passive Bodies in an Active Fluid

A single nonspherical body placed in an active fluid generates currents via breaking of time-reversal symmetry. We show that, when two or more passive bodies are placed in an active fluid, these currents lead to long-range interactions. Using a multipole expansion, we characterize their leading-order behaviors in terms of single-body properties and show that they decay as a power law with the distance between the bodies, are anisotropic, and do not obey an action-reaction principle. The interactions lead to rich dynamics of the bodies, illustrated by the spontaneous synchronized rotation of pinned nonchiral bodies and the formation of traveling bound pairs. The occurrence of these phenomena depends on tunable properties of the bodies, thus opening new possibilities for self-assembly mediated by active fluids.

Three-body correlations and conditional forces in suspensions of active hard disks

Self-propelled Brownian particles show rich out-of-equilibrium physics, for instance, the motility-induced phase separation (MIPS). While decades of studying the structure of liquids have established a deep understanding of passive systems, not much is known about correlations in active suspensions. In this work we derive an approximate analytic theory for three-body correlations and forces in systems of active Brownian disks starting from the many-body Smoluchowski equation. We use our theory to predict the conditional forces that act on a tagged particle and their dependence on the propulsion speed of self-propelled disks. We identify preferred directions of these forces in relation to the direction of propulsion and the positions of the surrounding particles. We further relate our theory to the effective swimming speed of the active disks, which is relevant for the physics of MIPS. To test and validate our theory, we additionally run particle-resolved computer simulations, for which we explicitly calculate the three-body forces. In this context, we discuss the modeling of active Brownian swimmers with nearly hard interaction potentials. We find very good agreement between our simulations and numerical solutions of our theory, especially for the nonequilibrium pair-distribution function. For our analytical results, we carefully discuss their range of validity in the context of the different levels of approximation we applied. This discussion allows us to study the individual contribution of particles to three-body forces and to the emerging structure. Thus, our work sheds light on the collective behavior, provides the basis for further studies of correlations in active suspensions, and makes a step towards an emerging liquid state theory.

Memory effects in active particles with exponentially correlated propulsion

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

Entropy Production and Fluctuation Theorems for Active Matter

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

Van't Hoff's law for active suspensions: the role of the solvent chemical potential

We extend Van't Hoff's law for the osmotic pressure to a suspension of active Brownian particles. The propelled particles exert a net reaction force on the solvent, and thereby either drive a measurable solvent flow from the connecting solvent reservoir through the semipermeable membrane, or increase the osmotic pressure and cause the suspension to rise to heights as large as micrometers for experimentally realized microswimmers described in the literature. The increase in osmotic pressure is caused by the background solvent being, in contrast to passive suspensions, no longer at the chemical potential of the solvent reservoir. The difference in solvent chemical potentials depends on the colloid-membrane interaction potential, which implies that the osmotic pressure is a state function of a state that itself is influenced by the membrane potential.

Statistical mechanics of transport processes in active fluids: Equations of hydrodynamics

The equations of hydrodynamics including mass, linear momentum, angular momentum, and energy are derived by coarse-graining the microscopic equations of motion for systems consisting of rotary dumbbells driven by internal torques. In deriving the balance of linear momentum, we find that the symmetry of the stress tensor is broken due to the presence of non-zero torques on individual particles. The broken symmetry of the stress tensor induces internal spin in the fluid and leads us to consider the balance of internal angular momentum in addition to the usual moment of momentum. In the absence of spin, the moment of momentum is the same as the total angular momentum. In deriving the form of the balance of total angular momentum, we find the microscopic expressions for the couple stress tensor that drives the spin field. We show that the couple stress contains contributions from both intermolecular interactions and the active forces. The presence of spin leads to the idea of balance of moment of inertia due to the constant exchange of particles in a small neighborhood around a macroscopic point. We derive the associated balance of moment of inertia at the macroscale and identify the moment of inertia flux that induces its transport. Finally, we obtain the balances of total and internal energy of the active fluid and identify the sources of heat and heat fluxes in the system.

Active Brownian equation of state: metastability and phase coexistence

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

Recovery of mechanical pressure in a gas of underdamped active dumbbells with Brownian noise

In contrast with a gas at thermodynamic equilibrium, the mean force exerted on a wall by a gas of active particles usually depends on the confining potential, thereby preventing a proper definition of mechanical pressure. In this paper, we investigate numerically the properties of a gas of underdamped self-propelled dumbbells subject to Brownian noise of increasing intensity, in order to understand how the notion of pressure is recovered as noise progressively masks the effects of self-propulsion and the system approaches thermodynamic equilibrium. The simulations performed for a mobile asymmetric wall separating two chambers containing an equal number of active dumbbells highlight some subtle and unexpected properties of the system. First, Brownian noise of moderate intensity is sufficient to let mean forces equilibrate for small values of the damping coefficient, while much stronger noise is required for larger values of the damping coefficient. Moreover, the displacement of the mean position of the wall upon increase of the intensity of the noise is not necessarily monotonous and may instead display changes of direction. Both facts actually reflect the existence of several mechanisms leading to the rupture of force balance, which tend to displace the mean position of the wall towards different directions and display different robustness against an increase of the intensity of Brownian noise. This work therefore provides a clear illustration of the fact that driving an autonomous system towards (or away from) thermodynamic equilibrium may not be a straightforward process, but may instead proceed through the variations of the relative weights of several conflicting mechanisms.

Microscopic derivation of the hydrodynamics of active-Brownian-particle suspensions

We derive the hydrodynamic equations of motion for a fluid of active particles described by underdamped Langevin equations that reduce to the active-Brownian-particle model, in the overdamped limit. The contraction into the hydrodynamic description is performed by locally averaging the particle dynamics with the nonequilibrium many-particle probability density, whose formal expression is found in the physically relevant limit of high friction through a multiple-time-scale analysis. This approach permits us to identify the conditions under which self-propulsion can be subsumed into the fluid stress tensor and thus to define systematically and unambiguously the local pressure of the active fluid.

Classical Nucleation Theory Description of Active Colloid Assembly

Nonaligning self-propelled particles with purely repulsive excluded volume interactions undergo athermal motility-induced phase separation into a dilute gas and a dense cluster phase. Here, we use enhanced sampling computational methods and analytic theory to examine the kinetics of formation of the dense phase. Despite the intrinsically nonequilibrium nature of the phase transition, we show that the kinetics can be described using an approach analogous to equilibrium classical nucleation theory, governed by an effective free energy of cluster formation with identifiable bulk and surface terms. The theory captures the location of the binodal, nucleation rates as a function of supersaturation, and the cluster size distributions below the binodal, while discrepancies in the metastable region reveal additional physics about the early stages of active crystal formation. The success of the theory shows that a framework similar to equilibrium thermodynamics can be obtained directly from the microdynamics of an active system, and can be used to describe the kinetics of evolution toward nonequilibrium steady states.

Active Particles with Soft and Curved Walls: Equation of State, Ratchets, and Instabilities

We study, from first principles, the pressure exerted by an active fluid of spherical particles on general boundaries in two dimensions. We show that, despite the nonuniform pressure along curved walls, an equation of state is recovered upon a proper spatial averaging. This holds even in the presence of pairwise interactions between particles or when asymmetric walls induce ratchet currents, which are accompanied by spontaneous shear stresses on the walls. For flexible obstacles, the pressure inhomogeneities lead to a modulational instability as well as to the spontaneous motion of short semiflexible filaments. Finally, we relate the force exerted on objects immersed in active baths to the particle flux they generate around them.