Communication: Green-Kubo approach to the average swim speed in active Brownian systems

Microscopic theory for nonequilibrium correlation functions in dense active fluids

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

Oscillatory Force Autocorrelations in Equilibrium Odd-Diffusive Systems

The force autocorrelation function (FACF), a concept of fundamental interest in statistical mechanics, encodes the effect of interactions on the dynamics of a tagged particle. In equilibrium, the FACF is believed to decay monotonically in time, which is a signature of slowing down of the dynamics of the tagged particle due to interactions. Here, we analytically show that in odd-diffusive systems, which are characterized by a diffusion tensor with antisymmetric elements, the FACF can become negative and even exhibit temporal oscillations. We also demonstrate that, despite the isotropy, the knowledge of FACF alone is not sufficient to describe the dynamics: the full autocorrelation tensor is required and contains an antisymmetric part. These unusual properties translate into enhanced dynamics of the tagged particle quantified via the self-diffusion coefficient that, remarkably, increases due to particle interactions.

Tracer dynamics in crowded active-particle suspensions

We discuss the dynamics of active Brownian particles (ABPs) in crowded environments through the mean-squared displacement (MSD) of active and passive tracer particles in both active and passive host systems. Exact equations for the MSD are derived using a projection operator technique, extending to dense systems the known solution for a single ABP. The interaction of the tracer particle with the host particles gives rise to strong memory effects. Evaluating these approximately in the framework of a recently developed mode-coupling theory for active Brownian particles (ABP-MCT), we discuss the various dynamical regimes that emerge: While self-propelled motion gives rise to super-diffusive MSD, at high densities, this competes with an interaction-induced sub-diffusive regime. The predictions of the theory are shown to be in good agreement with results obtained from an event-driven Brownian dynamics (ED-BD) simulation scheme for the dynamics of two-dimensional active Brownian hard disks.

How Activity Landscapes Polarize Microswimmers without Alignment Forces

Active-particle suspensions exhibit distinct polarization-density patterns in activity landscapes, even without anisotropic particle interactions. Such polarization without alignment forces is at work in motility-induced phase separation and betrays intrinsic microscopic activity to mesoscale observers. Using stable long-term confinement of a single thermophoretic microswimmer in a dedicated force-free particle trap, we examine the polarized interfacial layer at a motility step and confirm that it does not exert pressure onto the bulk. Our observations are quantitatively explained by an analytical theory that can also guide the analysis of more complex geometries and many-body effects.

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.

Inertial effects on the Brownian gyrator

The recent interest into the Brownian gyrator has been confined chiefly to the analysis of Brownian dynamics both in theory and experiment despite the applicability of general cases with definite mass. Considering mass explicitly in the solution of the Fokker-Planck equation and Langevin dynamics simulations, we investigate how inertia can change the dynamics and energetics of the Brownian gyrator. In the Langevin model, the inertia reduces the nonequilibrium effects by diminishing the declination of the probability density function and the mean of a specific angular momentum, j_{θ}, as a measure of rotation. Another unique feature of the Langevin description is that rotation is maximized at a particular anisotropy while the stability of the rotation is minimized at a particular anisotropy or mass. Our results suggest that the Langevin dynamics description of the Brownian gyrator is intrinsically different from that with Brownian dynamics. In addition, j_{θ} is proven to be essential and convenient for estimating stochastic energetics such as heat currents and entropy production even in the underdamped regime.

Transport coefficients in dense active Brownian particle systems: mode-coupling theory and simulation results

We discuss recent advances in developing a mode-coupling theory of the glass transition (MCT) of two-dimensional systems of active Brownian particles (ABPs). The theory describes the structural relaxation close to the active glass in terms of transient dynamical density correlation functions. We summarize the equations of motion that have been derived for the collective density-fluctuation dynamics and those for the tagged-particle motion. The latter allow to study the dynamics of both passive and active tracers in both passive and active host systems. In the limit of small wave numbers, they give rise to equations of motion describing the mean-squared displacements (MSDs) of these tracers and hence the long-time diffusion coefficients as a transport coefficient quantifying long-range tracer motion. We specifically discuss the case of a single ABP tracer in a glass-forming passive host suspension, a case that has recently been studied in experiments on colloidal Janus particles. We employ event-driven Brownian dynamics (ED-BD) computer simulations to test the ABP-MCT and find good agreement between the two for the MSD, provided that known errors in MCT already for the passive system (i.e., an overestimation of the glassiness of the system) are accounted for by an empirical mapping of packing fractions and host-system self-propulsion forces. The ED-BD simulation results also compare well to experimental data, although a peculiar non-monotonic mapping of self-propulsion velocities is required. The ABP-MCT predicts a specific self-propulsion dependence of the Stokes-Einstein relation between the long-time diffusion coefficient and the host-system viscosity that matches well the results from simulation. An application of ABP-MCT within the integration-through transients framework to calculate the density-renormalized effective swim velocity of the interacting ABP agrees qualitatively with the ED-BD simulation data at densities close to the glass transition and quantitatively for the full density range only after the mapping of packing fractions employed for the passive system.

Stochastic resetting of active Brownian particles with Lorentz force

Equilibrium properties of a system of passive diffusing particles in an external magnetic field are unaffected by Lorentz force. In contrast, active Brownian particles exhibit steady-state phenomena that depend on both the strength and the polarity of the applied magnetic field. The intriguing effects of the Lorentz force, however, can only be observed when out-of-equilibrium density gradients are maintained in the system. To this end, we use the method of stochastic resetting on active Brownian particles in two dimensions by resetting them to the line x = 0 at a constant rate and periodicity in the y direction. Under stochastic resetting, an active system settles into a nontrivial stationary state which is characterized by an inhomogeneous density distribution, polarization and bulk fluxes perpendicular to the density gradients. We show that whereas for a uniform magnetic field the properties of the stationary state of the active system can be obtained from its passive counterpart, novel features emerge in the case of an inhomogeneous magnetic field which have no counterpart in passive systems. In particular, there exists an activity-dependent threshold rate such that for smaller resetting rates, the density distribution of active particles becomes non-monotonic. We also study the mean first-passage time to the x axis and find a surprising result: it takes an active particle more time to reach the target from any given point for the case when the magnetic field increases away from the axis. The theoretical predictions are validated using Brownian dynamics simulations.

Photo-bioconvection: towards light control of flows in active suspensions

The persistent motility of individual constituents in microbial suspensions represents a prime example of the so-called active matter systems. Cells consume energy, exert forces and move, overall releasing the constraints of equilibrium statistical mechanics of passive elements and allowing for complex spatio-temporal patterns to emerge. Moreover, when subject to physico-chemical stimuli their collective behaviour often drives large-scale instabilities of a hydrodynamic nature, with implications for biomixing in natural environments and incipient industrial applications. In turn, our ability to exert external control of these driving stimuli could be used to govern the emerging patterns. Light, being easily manipulable and, at the same time, an important stimulus for a wide variety of microorganisms, is particularly well suited to this end. In this paper, we will discuss the current state, developments and some of the emerging advances in the fundamentals and applications of light-induced bioconvection with a focus on recent experimental realizations and modelling efforts. This article is part of the theme issue 'Stokes at 200 (part 2)'.

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.

Nondiffusive fluxes in a Brownian system with Lorentz force

The Fokker-Planck equation provides a complete statistical description of a particle undergoing random motion in a solvent. In the presence of Lorentz force due to an external magnetic field, the Fokker-Planck equation picks up a tensorial coefficient, which reflects the anisotropy of the particle's motion. This tensor, however, cannot be interpreted as a diffusion tensor; there are antisymmetric terms which give rise to fluxes perpendicular to the density gradients. Here, we show that for an inhomogeneous magnetic field these nondiffusive fluxes have finite divergence and therefore affect the density evolution of the system. Only in the special cases of a uniform magnetic field or carefully chosen initial condition with the same full rotational symmetry as the magnetic field can these fluxes be ignored in the density evolution.

Linear Response Theory and Green-Kubo Relations for Active Matter

We address the question of how interacting active systems in a nonequilibrium steady state respond to an external perturbation. We establish an extended fluctuation-dissipation theorem for active Brownian particles (ABP), which highlights the role played by the local violation of detailed balance due to activity. By making use of a Markovian approximation we derive closed Green-Kubo expressions for the diffusivity and mobility of ABP and quantify the deviations from the Stokes-Einstein relation. We compute the linear response function to an external force using unperturbed simulations of ABP and compare the results with the analytical predictions of the transport coefficients. Our results show the importance of the interplay between activity and interactions in the departure from equilibrium linear response.

Isotropic-nematic transition of self-propelled rods in three dimensions

Using overdamped Brownian dynamics simulations we investigate the isotropic-nematic (IN) transition of self-propelled rods in three spatial dimensions. For two well-known model systems (Gay-Berne potential and hard spherocylinders) we find that turning on activity moves to higher densities the phase boundary separating an isotropic phase from a (nonpolar) nematic phase. This active IN phase boundary is distinct from the boundary between isotropic and polar-cluster states previously reported in two-dimensional simulation studies and, unlike the latter, is not sensitive to the system size. We thus identify a generic feature of anisotropic active particles in three dimensions.

Mode-coupling theory for active Brownian particles

We present a mode-coupling theory (MCT) for the slow dynamics of two-dimensional spherical active Brownian particles (ABPs). The ABPs are characterized by a self-propulsion velocity v_{0} and by their translational and rotational diffusion coefficients D_{t} and D_{r}, respectively. Based on the integration-through-transients formalism, the theory requires as input only the equilibrium static structure factors of the passive system (where v_{0}=0). It predicts a nontrivial idealized-glass-transition diagram in the three-dimensional parameter space of density, self-propulsion velocity, and rotational diffusivity that arise because at high densities, the persistence length of active swimming ℓ_{p}=v_{0}/D_{r} interferes with the interaction length ℓ_{c} set by the caging of particles. While the low-density dynamics of ABPs is characterized by a single Péclet number Pe=v_{0}^{2}/D_{r}D_{t}, close to the glass transition the dynamics is found to depend on Pe and ℓ_{p} separately. At fixed density, increasing the self-propulsion velocity causes structural relaxation to speed up, while decreasing the persistence length slows down the relaxation. The active-MCT glass is a nonergodic state that is qualitatively different from the passive glass. In it, correlations of initial density fluctuations never fully decay, but also an infinite memory of initial orientational fluctuations is retained in the positions.

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.

Brownian systems with spatially inhomogeneous activity

We generalize the Green-Kubo approach, previously applied to bulk systems of spherically symmetric active particles [J. Chem. Phys. 145, 161101 (2016)JCPSA60021-960610.1063/1.4966153], to include spatially inhomogeneous activity. The method is applied to predict the spatial dependence of the average orientation per particle and the density. The average orientation is given by an integral over the self part of the Van Hove function and a simple Gaussian approximation to this quantity yields an accurate analytical expression. Taking this analytical result as input to a dynamic density functional theory approximates the spatial dependence of the density in good agreement with simulation data. All theoretical predictions are validated using Brownian dynamics simulations.