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

Transition-path sampling for run-and-tumble particles

We elaborate and validate a generalization of the renowned transition-path-sampling algorithm for a paradigmatic model of active particles, namely, the run-and-tumble particles. Notwithstanding the nonequilibrium character of these particles, we show how the consequent lack of the microscopical reversibility property, which is usually required by transition-path sampling, can be circumvented by identifying reasonable backward dynamics with a well-defined path-probability density. Our method is then applied to characterize the structure and kinetics of rare transition pathways undergone by run-and-tumble particles having to cross a potential barrier in order to find a target.

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.

From a microscopic solution to a continuum description of active particles with a recoil interaction in one dimension

We consider a model system of persistent random walkers that can jam, pass through each other, or jump apart (recoil) on contact. In a continuum limit, where particle motion between stochastic changes in direction becomes deterministic, we find that the stationary interparticle distribution functions are governed by an inhomogeneous fourth-order differential equation. Our main focus is on determining the boundary conditions that these distribution functions should satisfy. We find that these do not arise naturally from physical considerations, but they need to be carefully matched to functional forms that arise from the analysis of an underlying discrete process. The interparticle distribution functions, or their first derivatives, are generically found to be discontinuous at the boundaries.

Critical behavior of active Brownian particles: Connection to field theories

We explore the relation between active Brownian particles, a minimal particle-based model for active matter, and scalar field theories. Both show a liquid-gas-like phase transition toward stable coexistence of a dense liquid with a dilute active gas that terminates in a critical point. However, a comprehensive mapping between the particle-based model parameters and the effective coefficients governing the field theories has not been established yet. We discuss conflicting recent numerical results for the critical exponents of active Brownian particles in two dimensions. Starting from the intermediate effective hydrodynamic equations, we then present a construction for a scalar order parameter for active Brownian particles that yields the active model B+. We argue that a crucial ingredient is the coupling between density and polarization in the particle current. The renormalization flow close to two dimensions exhibits a pair of perturbative fixed points that limit the attractive basin of the Wilson-Fisher fixed point, with the perspective that the critical behavior of active Brownian particles in two dimensions is governed by a strong-coupling fixed point different from Wilson-Fisher and not necessarily corresponding to Ising universality.

Pair-distribution function of active Brownian spheres in two spatial dimensions: Simulation results and analytic representation

We investigate the full pair-distribution function of a homogeneous suspension of spherical active Brownian particles interacting by a Weeks-Chandler-Andersen potential in two spatial dimensions. The full pair-distribution function depends on three coordinates describing the relative positions and orientations of two particles, the Péclet number specifying the activity of the particles, and their mean packing density. This five-dimensional function is obtained from Brownian dynamics simulations. We discuss its structure taking into account all of its degrees of freedom. In addition, we present an approximate analytic expression for the product of the full pair-distribution function and the interparticle force. We find that the analytic expression, which is typically needed when deriving analytic models for the collective dynamics of active Brownian particles, is in good agreement with the simulation results. The results of this work can thus be expected to be helpful for the further theoretical investigation of active Brownian particles as well as nonequilibrium statistical physics in general.

Predictive local field theory for interacting active Brownian spheres in two spatial dimensions

We present a predictive local field theory for the nonequilibrium dynamics of interacting active Brownian particles with a spherical shape in two spatial dimensions. The theory is derived by a rigorous coarse-graining starting from the Langevin equations that describe the trajectories of the individual particles. For high accuracy and generality of the theory, it includes configurational order parameters and derivatives up to infinite order. In addition, we discuss possible approximations of the theory and present reduced models that are easier to apply. We show that our theory contains popular models such as Active Model B+  as special cases and that it provides explicit expressions for the coefficients occurring in these and other, often phenomenological, models. As a further outcome, the theory yields an analytical expression for the density-dependent mean swimming speed of the particles. To demonstrate an application of the new theory, we analyze a simple reduced model of the lowest nontrivial order in derivatives, which is able to predict the onset of motility-induced phase separation of the particles. By a linear stability analysis, an analytical expression for the spinodal corresponding to motility-induced phase separation is obtained. This expression is evaluated for the case of particles interacting repulsively by a Weeks-Chandler-Andersen potential. The analytical predictions for the spinodal associated with these particles are found to be in very good agreement with the results of Brownian dynamics simulations that are based on the same Langevin equations as our theory. Furthermore, the critical point predicted by our analytical results agrees excellently with recent computational results from the literature.

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.

Response of active Brownian particles to boundary driving

We computationally study the behavior of underdamped active Brownian particles in a sheared channel geometry. Due to their underdamped dynamics, the particles carry momentum a characteristic distance away from the boundary before it is dissipated into the substrate. We correlate this distance with the persistence of particle trajectories, determined jointly by their friction and self-propulsion. Within this characteristic length, we observe counterintuitive phenomena stemming from the interplay of activity, interparticle interactions, and the boundary driving. Depending on the values of friction and self-propulsion, interparticle interactions can either aid or hinder momentum transport. More dramatically, in certain cases we observe a flow reversal near the wall, which we correlate with an induced polarization of the particle self-propulsion directions. We rationalize these results in terms of a simple kinetic picture of particle trajectories.

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.

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.

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.