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

TL;DR

This paper derives and verifies local stress expressions in active Brownian particle systems, showing that active stress is due to momentum transport and establishing an equation of state for such inhomogeneous active fluids.

Contribution

The paper introduces a novel local stress expression for active Brownian particles and demonstrates its validity and physical origin through simulations.

Findings

Local bulk stress equals wall stress in confined systems.

Active stress results from momentum transport, not boundary effects.

Equation of state exists for spherical ABPs.

Abstract

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.

Equations

$$m\ddot{{\mathbf{r}}_i}+\gamma \dot{{\mathbf{r}}_i}(t)=\gamma v_0{\mathbf{e}}_i(t)+{\mathbf{F}}_i(t)+{\Gamma }_i(t),$$

$$⟨{\Gamma }_i(t)⟩=0,⟨{\Gamma }_{i\alpha }(t){\Gamma }_{j\beta }(t^{\prime })⟩=2\gamma k_{B}T{\delta }_{\alpha \beta }{\delta }_{ij}\delta (t-t^{\prime }).$$

$$\dot{{\mathbf{e}}_i}(t)={\mathbf{e}}_i(t)\times {\mathbf{\eta }}_i(t)-2D_R{\mathbf{e}}_i$$

$$⟨{\eta }_{i\alpha }(t)⟩=0,⟨{\eta }_{i\alpha }(t){\eta }_{j\beta }(t^{\prime })⟩={\gamma }_R{\delta }_{\alpha \beta }{\delta }_{ij}\delta (t-t^{\prime }).$$

$$⟨{\mathbf{v}}_i^a(t)\cdot {\mathbf{v}}_i^a(0)⟩=v_0^2⟨{\mathbf{e}}_i(t)\cdot {\mathbf{e}}_i(0)⟩=v_0^2{e}^{-{\gamma }_{R}t}.$$

$${\mathbf{F}}_i={\mathbf{F}}_i^w+\text{sideset}{}^{\prime }{\sum }_{j=1}^N{\mathbf{F}}_{ij}.$$

$$\frac{d}{dt}\left[m⟨{\dot{r}}_{i\alpha }{r}_{i\alpha }⟩+\frac{\gamma }{2}⟨r_{i\alpha }^2⟩\right]=m⟨{\dot{r}}_{i\alpha }^2⟩+\gamma ⟨v_{i\alpha }^a{r}_{i\alpha }⟩+⟨F_{i\alpha }^w{r}_{i\alpha }⟩+\text{sideset}{}^{\prime }{\sum }_{j=1}^N⟨F_{ij\alpha }{r}_{i\alpha }⟩,$$

$$V{\sigma }_{\alpha \alpha }^e=\sum _{i=1}^N⟨F_{i\alpha }^w{S}_{i\alpha }⟩,$$

$$V{\sigma }_{\alpha \alpha }^i=-\sum _{i=1}^{N}m⟨{\dot{r}}_{i\alpha }^2⟩-\sum _{i=1}^N\gamma ⟨v_{i\alpha }^a{r}_{i\alpha }⟩-\sum _{i=1}^N⟨F_{i\alpha }^w(r_{i\alpha }-S_{i\alpha })⟩-\frac{1}{2}\sum _{i=1}^N\text{sideset}{}^{\prime }{\sum }_{j=1}^N⟨F_{ij\alpha }{r}_{ij\alpha }⟩,$$

$$V{\sigma }_{\alpha \alpha }^s=-\sum _{i=1}^N\gamma ⟨v_{i\alpha }^a{r}_{i\alpha }⟩.$$

$$\sum _{i=1}^{N}m⟨{\dot{r}}_i^2{\Lambda }_i⟩+\sum _{i=1}^{N}m⟨{\dot{r}}_i{r}_i{\dot{\Lambda }}_i⟩+\sum _{i=1}^N⟨F_i{r}_i{\Lambda }_i⟩+\sum _{i=1}^N⟨{\Gamma }_i{r}_i{\Lambda }_i⟩+\gamma \sum _{i=1}^N⟨v_i^a{r}_i{\Lambda }_i⟩=0,$$

$$\sum _{i=1}^N⟨e_i{\dot{r}}_i{\Lambda }_i⟩+\sum _{i=1}^N⟨e_i{r}_i{\dot{\Lambda }}_i⟩-{\gamma }_R\sum _{i=1}^N⟨e_i{r}_i{\Lambda }_i⟩=0,$$

$$\frac{d}{dt}⟨{\dot{r}}_i{r}_i{\Lambda }_i⟩$$

$$\sum _{i=1}^{N}m⟨{\dot{r}}_i^2{\Lambda }_i⟩+\sum _{i=1}^{N}m⟨{\dot{r}}_i{r}_i{\dot{\Lambda }}_i⟩+\sum _{i=1}^N⟨F_i{r}_i{\Lambda }_i⟩+\frac{\gamma }{{\gamma }_R}\left[\sum _{i=1}^N⟨v_i^a{\dot{r}}_i{\Lambda }_i⟩+\sum _{i=1}^N⟨v_i^a{r}_i{\dot{\Lambda }}_i⟩\right]=0.$$

$${\mathbf{p}}_i^a=\frac{\gamma }{{\gamma }_R}{\mathbf{v}}_i^a$$

$$\Delta V{\sigma }_{\alpha \alpha }^e=$$

$$\Delta V{\sigma }_{\alpha \alpha }^i=$$

$$\Delta V{\sigma }_{\alpha \alpha }^{if}=-\frac{1}{2}\sum _{i=1}^N\text{sideset}{}^{\prime }{\sum }_{j=1}^N⟨{\lambda }_{ij}{r}_{ij}{F}_{ij}⟩$$

$$\Delta V{\sigma }_{\alpha \alpha }^s=-\frac{\gamma }{{\gamma }_R}\sum _{i=1}^N⟨v_{i\alpha }^a{\dot{r}}_{i\alpha }{\Lambda }_{i\alpha }⟩,$$

$$\sum _{i=1}^N⟨p_{i\alpha }^a{r}_{i\alpha }{\dot{\Lambda }}_{i\alpha }⟩=-\frac{\gamma }{{\gamma }_R}\sum _{i=1}^N⟨v_{i\alpha }^a{\dot{r}}_{i\alpha }{\Lambda }_{i\alpha }⟩.$$

$$V{\sigma }_{\alpha \alpha }^i=-N{k}_{B}T-\sum _{i=1}^N\gamma ⟨v_{i\alpha }^a{r}_{i\alpha }⟩-\sum _{i=1}^N⟨F_{i\alpha }^w(r_{i\alpha }-S_{i\alpha })⟩,$$

$$\gamma ⟨v_{i\alpha }^a{r}_{i\alpha }⟩=\frac{\gamma }{d{\gamma }_R}{v}_0^2+\frac{1}{{\gamma }_R}⟨F_{i\alpha }^w{v}_{i\alpha }^a⟩,$$

$$V{\sigma }_{\alpha \alpha }^i=-N{k}_{B}T-\frac{\gamma }{3{\gamma }_R}N{v}_0^2-\frac{1}{{\gamma }_R}\sum _{i=1}^N⟨F_{i\alpha }^w{v}_{i\alpha }^a⟩-\sum _{i=1}^N⟨F_{i\alpha }^w(r_{i\alpha }-L_{\alpha })⟩$$

$${\sigma }_{\alpha \alpha }^{id}=-\frac{\gamma N}{3{\gamma }_{R}V}{v}_0^2.$$

$$\Delta V{\sigma }_{\alpha \alpha }^i=-\sum _{i=1}^{N}m⟨{\dot{r}}_{i\alpha }^2{\Lambda }_{i\alpha }⟩-\frac{\gamma }{{\gamma }_R}\sum _{i=1}^N⟨v_{i\alpha }^a{\dot{r}}_{i\alpha }{\Lambda }_{i\alpha }⟩.$$

$${\dot{r}}_{i\alpha }=\frac{1}{m}{\int }_{-\infty }^t{e}^{-\gamma (t-t^{\prime })/m}\left[\gamma v_{i\alpha }^a(t^{\prime })+F_{i\alpha }(t^{\prime })+{\Gamma }_{i\alpha }(t^{\prime })\right]d{t}^{\prime },$$

$$⟨{\dot{r}}_{i\alpha }^2⟩=\frac{k_{B}T}{m}+\frac{v_0^2}{d(1+m{\gamma }_R/\gamma )},⟨v_{i\alpha }^a{\dot{r}}_{i\alpha }⟩=\frac{v_0^2}{d(1+m{\gamma }_R/\gamma )}$$

$$\Delta V{\sigma }_{\alpha \alpha }^i=-N_{\Delta V}\left(k_{B}T+\frac{\gamma }{d{\gamma }_R}{v}_0^2\right),$$

$$\Delta V{\sigma }_{\alpha \alpha }^i=-\frac{\gamma }{{\gamma }_R}\sum _{i=1}^{N_{\Delta V}}⟨(v_{i\alpha }^a{)}^2{\Lambda }_{i\alpha }⟩=-N_{\Delta V}\frac{\gamma }{3{\gamma }_R}{v}_0^2$$

$$\frac{{\sigma }_{zz}^i{\sigma }^2}{\gamma v_0}=-\frac{{\rho }_{\Delta V}{\sigma }^3}{6}Pe,$$

$$\frac{\Delta \rho }{\rho }=\frac{1}{1+0.7L_z/l_p},$$

Full text

Introduction

The agents of active matter perpetually convert internal energy or energy from the environment into systematic translational motion and are, thus, usually far from thermal equilibrium [1, 2, 3, 4, 5]. Their motility gives rise to remarkable phenomena such as collective motion, activity-induced phase transitions [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 5, 17], or wall accumulation [18, 19, 20, 21]. Activity is associated with the generation of active stresses, which are responsible for the persistent nonequilibrium state of a system [22, 23, 24]. Such stresses are particularly interesting, since there is no equilibrium counterpart and thus, they usually cannot be described by equilibrium thermodynamics. An expression for the active pressure in systems of active Brownian particles (ABPs), based on a generalization of the virial including the active force, has been introduced recently [22, 23]. Subsequently, expressions for the active pressure have been derived in various ways for both, the overdamped dynamics of ABPs [25, 26, 27, 28] as well as in presence of inertia [29, 30].

Mechanical stress/pressure as force per area can be defined even far from equilibrium. The fundamental question is then, how this mechanical pressure is related to nonequilibrium thermodynamics. As a step toward such a relation, the existence of an equation of state for systems of ABPs has been addressed [24, 28, 31]. Various simulation and theoretical studies confirm that indeed such a relation, including the active pressure, exists for spherical ABPs [23, 24, 22, 31, 32, 25, 33]. However, for (anisotropic) active particles which experience a torque, pressure is in general not a state function anymore [28, 34].

A pressure equation of state implies equivalence of wall and local pressure at any point inside the homogeneous relevant volume. Exactly this identity has been questioned in Ref.  27. Typically, confined or periodic systems are considered and relations between wall forces, hence the wall pressure, and virial, including the active-force virial, are determined [23, 28, 22, 25], but no local (bulk) pressure is calculated. Aside from Ref.  27, local pressure and the existence of an equation of state has been addressed in Ref.  35.

In this article, we derive an instantaneous expression for the local stress (pressure) in a system of spherical active Brownian particles using the virial approach. Thereby, we follow the strategy of Ref.  36 for passive systems. A volume $\Delta V$ within a larger volume $V$ is considered and an expression is derived for the pressure inside $\Delta V$. Considering the equations of motions of the ABPs inside the volume $\Delta V$, including inertia, we obtain the well-known virial expressions for the pressure of a passive system. In addition, however, we find activity contributions. On the one hand, we obtain a virial expression for the active stress (swim stress) of particles inside $\Delta V$, which differs from the swim stress discussed so far, and, on the other hand, we find an original momentum-flux contribution of APBs to stress, which we denote as active momentum. The notation rests upon the similarity to the classical linear momentum of a passive particle. It is precisely this term, which distinguishes our studies from those of Ref.  27. In general, we agree with the arguments presented in Ref.  27, however, the momentum flux across a plane in the system, similar to the momentum flux by inertia, yields an extra contribution to the local stress which not considered in Ref.  27. In fact, in Ref.  35 an “active impulse, defined as the mean momentum a particle will receive on average from the substrate in the future”, is introduce, resembling our instantaneous active momentum, but the meaning is different. By computer simulations of non-interacting ABPs and ABPs interacting by a Lennard-Jones potential, we demonstrate that the derived expressions for the wall and bulk stress are identical. This highlights the existence of a pressure equation state for spherical ABP systems [22, 24, 35]. Moreover, it opens up the possibility to calculate stresses locally, even in inhomogeneous systems, which we confirm by calculating the local stress in a phase separated system.

Stress in active systems

Model

We consider a system of $N$ ABPs of diameter $\sigma$ in a three-dimensional volume $V$. Their equations of motion, including the inertia contribution, are given by

[TABLE]

where ${\bm{r}}_{i}$, $\dot{\bm{r}_{i}}$, and ${\bm{F}}_{i}$ ($i=1,\ldots,N$) denote the particle positions, velocities, and forces, respectively, $m$ is the mass, $\gamma$ the translational friction coefficient, and ${\bm{\varGamma}}_{i}$ are Gaussian white-noise random forces, with the moments

[TABLE]

Here, $T$ denotes the temperature, $k_{B}$ the Boltzmann constant, and $\alpha,\beta\in\{x,y,z\}$ the Cartesian coordinate directions. A particle $i$ is propelled with constant velocity $v_{0}$ along its orientation vector ${\bm{e}}_{i}$. The latter changes in a diffusive manner according to

[TABLE]

within the Ito interpretation of the multiplicative noise process [37]; $\bm{\eta}_{i}$ is a Gaussian and Markovian stochastic process, with the moments

[TABLE]

The damping factor $\gamma_{R}$ is related to the rotational diffusion coefficient, $D_{R}$, via $\gamma_{R}=2D_{R}$. Alternatively, the velocity equation of motion of an active Ornstein-Uhlenbeck particle can be considered [38, 39], where aside from the orientation also the magnitude of the propulsion velocity is changing. In any case, the correlation function for the active velocity $\bm{v}_{i}^{a}=v_{0}\bm{e}_{i}$ is [38, 25]

[TABLE]

The forces $\bm{F}_{i}$ include contributions from (short-range) interactions with confining walls, $\bm{F}_{i}^{w}$, as well as possible pair-wise ABP interactions, $\bm{F}_{ij}$, such that (the prime at the sum indicates that the index $j=i$ is excluded)

[TABLE]

Global stress in confined system

We examine now ABPs confined in a rectangle (2D) or cuboid (3D) with walls located at $\pm L_{\alpha}/2$. The virial expression for the stress follows by multiplying the respective individual Cartesian equations of Eqs. (1) by $r_{i\alpha}$, which yields [40, 25]

[TABLE]

where $\left\langle\ldots\right\rangle$ denotes either a time or an ensemble average. Due to confinement, the left-hand side of Eq. (7) vanishes. With the definition of an external stress $\sigma_{\alpha\alpha}^{e}$ as average total force exerted on a wall,

[TABLE]

where $S_{i\alpha}=\pm L_{\alpha}/2$ denotes the location of the wall with which particle $i$ interacts, Eq. (7) yields the equivalent internal stress

[TABLE]

with $\bm{r}_{ij}=\bm{r}_{i}-\bm{r}_{j}$, as sum over all particles of the confined systems. Hence, we denote this expression as global (wall) stress. Note that by definition $\sigma_{\alpha\alpha}^{e}=\sigma_{\alpha\alpha}^{i}$. The contribution with the active force is denoted as swim stress [22, 41, 25]

[TABLE]

The pressure itself follows as $p^{i/e}=-(\sum_{\alpha}\sigma_{\alpha\alpha}^{i/e})/3$. Aside from the wall term with the forces $\bm{F}_{i}^{w}$, the stress tensor (9) is identical with that derived in Ref.  25. The wall term captures the finite range of the wall force and vanishes in the limit of vanishing wall interaction range.

Local virial and local stress

In order to determine the local stress at the position $\bm{r}$ within the volume $V$, we consider a cuboidal subvolume $\Delta V=\Delta L_{x}\times\Delta L_{y}\times\Delta L_{z}$ centered around $\bm{r}$, with the side lengths $\Delta L_{\alpha}$ (cf. Fig. 1). ABPs inside of $\Delta V$ interact with each other as well as ABPs outside. Moreover, ABPs enter and leave $\Delta V$ in the course of time. The internal stress follows from the suitable virial expression of ABPs inside of $\Delta V$. To derive the virial expression, as before, we multiply the equations (1) and (3) by $r_{i\alpha}$, an additional factor $\Lambda_{i\alpha}(\bm{r})$, and sum over all particles. The factor $\Lambda_{i\alpha}$ determines the volume $\Delta V$—$\Lambda_{i\alpha}(\bm{r})$ is unity when particle $i$ is within $\Delta V$ in the coordinate direction $\alpha$ and zero otherwise [36]. As a result, averaging yields

[TABLE]

with the Ito interpretation of the colored noise process for $\bm{e}$ [42, 25], and the relations

[TABLE]

Here, and in the following, we often suppress the index $\alpha$ for compactness of the expressions. Due to confinement, the left-hand sides of Eqs. (13) vanish. Similarly, the averages $\left\langle\dot{r}_{i}r_{i}\Lambda_{i}\right\rangle=\langle r_{i}^{2}\dot{\Lambda}_{i}\rangle=0$. Inserting the third term of Eq. (12) into Eq. (11), we obtain the virial expression

[TABLE]

Here, we neglect terms with the averages $\langle\Gamma_{i}r_{i}\Lambda_{i}\rangle$; these averages vanish as long as inertia (cf. Eq. (1)) is taken into account. For an overdamped dynamics, $m=0$ in Eq. (1), $\langle\Gamma_{i}r_{i}\Lambda_{i}\rangle$ yields the thermal contribution to the stress [25].

For passive systems with $v_{0}=0$, Eq. (14) reduces to the local virial expression in thermal equilibrium presented in Ref.  36. Compared to the global virial [25], additional terms with the derivative of $\Lambda_{i}$, $\dot{\Lambda}_{i}$, appear due to particles entering and leaving the local volume $\Delta V$. Most remarkably, activity yields a similar contribution as inertia, namely $\gamma\sum_{i}\langle v_{i\alpha}^{a}r_{i\alpha}\dot{\Lambda}_{i}\rangle/\gamma_{R}$, compared to $\sum_{i}\langle m\dot{r}_{i\alpha}r_{i\alpha}\dot{\Lambda}_{i}\rangle$. This demonstrates that transport of ABPs across a plane yields a contribution to stress similarly to the momentum $\bm{p}_{i}^{m}=m\dot{\bm{r}}_{i}$ of a (passive) particle, in contrast to previous statements [27]. Hence, we can consider

[TABLE]

as an instantaneous active momentum. Formally, an analogous momentum (impulse) has been introduced in Ref.  35 “as momentum the active particle will receive on average from its active force in the future”, thus, the meaning is rather different.

Equation (14) includes volume and boundary contributions. Thereby, terms with $\dot{\Lambda}_{i}$ are pure boundary terms. Following the procedure of Ref.  36, the following external, $\sigma_{\alpha\alpha}^{e}$, and internal, $\sigma_{\alpha\alpha}^{i}$, stress tensor components are identified,

[TABLE]

Here, we assume that the volume $\Delta V$ is sufficiently far from the surface, and, hence, short-range wall forces do not contribute to the local stress in $\Delta V$. The external stress tensor $\sigma_{\alpha\alpha}^{e}$ includes contributions by momentum flux across the boundary, with the active flux $\gamma\sum_{i=1}^{N}\langle v^{a}_{i}r_{i}\dot{\Lambda}_{i}\rangle/\gamma_{R}$, as well as boundary force terms; the third term on the right-hand side of Eq. (16) accounts for stress contributions by forces between particles $i$ inside and $j$ outside of $\Delta V$, and the fourth term for force contributions, when all particles are outside of $\Delta V$, but the connecting line insects the volume $\Delta V$; $\lambda_{ij}$ denotes the fraction of the line connecting particle $i$ and $j$ inside of the volume $\Delta V$ (cf. Fig. 1) [36]. Note, by suppression of the averages, instantaneous stresses are obtained.

The internal stress tensor $\sigma_{\alpha\alpha}^{i}$ (17) includes kinetic contributions from inertia and activity of particles, which are inside the volume of interest, as well as interparticle forces. The interparticle-force term

[TABLE]

includes contributions from particles both, inside of $\Delta V$, here $\lambda_{ij}=1$—this is the standard volume term—, as well as combinations, where a particle is inside and all others, within the finite interaction range, are outside, or both are outside, but their connecting line intersects the volume $\Delta V$, here $\lambda_{ij}<1$. The latter are boundary terms in analogy to the wall contributions in Eq. (9) of the global system with finite-range wall forces [40].

Remarkably, a new expression for the swim stress is obtained

[TABLE]

denote as local swim stress in the following, different from the “traditional” swim stress of Eq. (10). Most importantly, the swim stress (10) vanishes in $\Delta V$, and, hence, it is not the adequate expression for the calculation of the local active stress. The vanishing expression (10) in a subvolume lead to the conclusion in Ref.  27 that the local pressure in an ideal active gas is identical with the passive-ideal-gas pressure when accounting for the respective local density. As a consequence, it was concluded that the active pressure is a boundary effect, because active forces were not contributing to momentum flux. Our calculations contradict this conclusions. First of all, we find an active flux, Eq. (15), and, secondly, another expression for the local active stress, Eq. (19). The vanishing average of Eq. (10) in $\Delta V$ together with Eq. (12) shows that the contribution by the active momentum is the negative of the active virial, i.e.,

[TABLE]

However, if we extend the volume $\Delta V$ to the whole volume $V$, i.e., $\Lambda_{i}\equiv 1$, Eq. (3) in the Ito formulation yields $\left\langle v^{a}_{i\alpha}\dot{r}_{i\alpha}\right\rangle=\gamma_{R}\left\langle v^{a}_{i\alpha}r_{i\alpha}\right\rangle$, and Eqs. (9) and (17) become identical. Note that the wall force has then to be taken into account in Eq. (17). Thus, our Eq.  (19) is a more general expression for the swim stress, since it applies locally as well as globally.

Global stress of ideal ABP gas

For an ideal gas of ABPs, the interparticle forces are zero, i.e., $\bm{F}_{ij}=0$, and Eq. (9) becomes

[TABLE]

in the overdamped limit, $m=0$. According to Refs.  25, 24, the swim stress (10) can be written as

[TABLE]

where $d$ is the dimension, which leads to the expression

[TABLE]

in 3D. The first term on the right-hand side is the thermal contribution to stress, and the second term the ideal gas swim stress [25, 24, 22, 43, 44]

[TABLE]

Noteworthy, the wall interaction appears twice in Eq. (23). The term proportional to the active velocity reflects and accounts for wall accumulation, whereas the other term captures finite-range effects of the wall force.

Local tress of ideal ABP gas

For zero interparticle forces, Eq. (17) becomes

[TABLE]

With the formal solution of Eq. (1),

[TABLE]

we find

[TABLE]

in the stationary state. In the evaluation of these expressions, all wall contributions vanish as long as we assume that the volume $\Delta V$ is located far enough from a wall such that $e^{-\gamma(t-t^{w}_{i})}\ll 1$, where $t_{i}^{w}$ is the time of the last encounter of an ABP with a wall before entering $\Delta V$. Thus, we obtain

[TABLE]

with $N_{\Delta V}$ the number of ABPs in the volume $\Delta V$. The contributions with the Stokes number $m\gamma_{R}/\gamma$ cancel, and the result for the overdamped dynamics is obtained, as already discussed in Ref.  30. Evidently, the local stress is independent of the wall force. However, due to wall accumulation, the ABP density is not necessarily homogeneous over the whole volume, in particular, it is reduced in the bulk. Hence, walls can affect the bulk stress and pressure.

For an overdamped and athermal dynamics, Eq. (25) becomes

[TABLE]

in 3D by inserting $\dot{r}_{i\alpha}=v^{a}_{i\alpha}$ from Eq. (1), Hence, stress is given by the square of the (active) velocity, similar to the case of a passive system.

Evidently, Eqs. (23) and (28) become identical in the limit $L_{\alpha}\to\infty$, i.e., an infinite volume corresponding to the thermodynamic limit in case of a passive system. Then, the local particle density $N_{\Delta V}/\Delta V$ is equal to the density $N/V$, since finite-size effects and the inhomogeneities by surface accumulation become negligible compared to the overall particle number. However, the equality of the two expressions (23) and (25) applies for all surface separations as long as $(L_{\alpha}-\Delta L_{\alpha})/l_{p}\gg 1$, where $l_{p}=v_{0}/\gamma_{R}$ is the persistence length of the active motion. The surface forces in Eq. (23) yield additional contributions to the ideal active gas stress. Associated bulk-density variations change the respective values of the local stress (28).

Figure 2 shows simulation results for the local swim stress of non-interacting ABPs in the subvolume $\Delta V$ according to Eq. (19) as function of the center position $r_{z}$ of $\Delta V$, where $\Delta L_{x}=\Delta L_{y}=100\sigma$ and $\Delta L_{z}=0.05L_{z}$. Activity is characterized by the Péclet number $Pe=v_{0}/\sigma D_{R}$, where $Pe=20$ in Fig. 2. In addition, the global wall stress according to Eq. (21) ($k_{B}T=0$) is presented. Evidently, the local stress is independent of the location of the subvolume $\Delta V$ in the bulk of the system and it agrees with the wall stress within $1\%$. Our results indicate that stress in an active system is not a wall effect. According to Eq. (29), the local stress is proportional to the square of the active velocity of the ABPs within $\Delta V$, similar to the dependence on the kinetic energy in passive systems. Moreover, the simulations confirm the presence of an active momentum flux with the momentum of Eq. (15).

The dependence of the stress on the activity is presented in Fig. 3(a). The local swim stress (19) coincides perfectly with the wall stress (23) ($k_{B}T=0$) of the system. The scaled stress

[TABLE]

where $\rho_{\Delta V}=N_{\Delta V}/\Delta V$ is the bulk density, decreases monotonically with increasing Péclet number, but deviates from the linear dependence of an unconfined (periodic boundary conditions) ideal ABP gas for large $Pe$.[45, 46, 27] The stress reduction is a consequence of an enhanced ABP surface accumulation with increasing $Pe$ and, correspondingly, a reduction of the bulk density. Here, we like to emphasize that the stresses calculated via Eqs. (8) and (9) perfectly agree with each other.

Figure 3(b) shows the relative density difference $\Delta\rho/\rho=(\rho-\rho_{\Delta V})/\rho$ with respect to the mean density $\rho=N/V$ as function of the Péclet number (squares). The simulation data, with the bulk density $\rho_{\Delta V}$, can well be described by the relation

[TABLE]

which approaches zero in the limits $L_{z}\to\infty$ and $Pe\to 0$, and the bulk density becomes equal to the average density $\rho$ as for an passive ideal gas. Such a dependence on $L_{z}/l_{p}$ has been obtained in Ref.  27 for the fraction of adsorbed ABPs at a wall. However, the numerical factor is somewhat different, we find $0.7$ compared to $0.6$ in Ref.  27. The difference could be related to the different adopted dimensionalities, $d=3$ here versus $d=2$ in Ref.  27. Moreover, Fig. 3(b) shows that the wall virial $\sum_{i}\left\langle F_{i\alpha}^{w}v^{a}_{i\alpha}\right\rangle/\gamma_{R}$ (23) divided by the ideal ABP swim stress (24) is identical with the relative density difference $\Delta\rho/\rho$ and, thus, exhibits the same dependence on the Péclet number (the finite-range wall-force contribution is negligible). Of course, this is expected by the equivalence of the stress expressions (21) and (25).

Figure 4 illustrates the dependence of the stress on the wall separation $L_{z}$ for the Péclet number $Pe=40$. Again, simulation results for wall and local stress agree very well. With increasing system size, the wall effect decreases and the stress approaches the value of the ideal bulk stress. The relation between the local density and the system size is displayed in the inset of Fig. 4, which is well described by Eq. (31). The dependence of the ideal active gas stress on the Péclet number and system size is solely described by the respective dependence of the density $\rho_{\Delta V}$ on these quantities, since $\sigma_{zz}^{i}=\sigma_{zz}^{id}\rho_{\Delta V}/\rho$. This is confirmed in Figs. 3(a) and 4.

Some of our findings seem to be rather evident as soon as the validity of the Eq. (29) is taken for granted. However, this expression follows from our local stress with the local swim stress (19), and not the stress expression for the system with walls. Specifically, Ref.  27 discusses wall-induced effects only. Hence, by our studies we emphasize the existence of a bulk stress and pressure and its identity with the respective wall values.

Interacting ABPs

Simulation results for the local interparticle force (18), local active (19), and global stress of ABP systems with excluded-volume interactions are depicted in Fig 5(a) as function of the packing fractions $\phi$. Again, we obtain excellent agreement between the various stress-tensor contributions calculated both, globally and locally. Consistent with previous studies, the swim stress contribution decreases with increasing density, whereas the interparticle force contribution monotonically increases over the considered concentration range [25, 16, 47]. Above $\phi\approx 0.4$, the contribution due to interparticle forces dominates the overall stress. Figure 5(b) compares the local swim and interparticle stresses with those in the overall system for various wall distances and concentrations. Remarkably, the local swim stress (bullets) is independent of the wall separation, except for $\phi=0.25$. In this case, the volume $\Delta V$ is not sufficiently far away from the walls, and due to the comparably large persistence length $l_{p}=5\sigma$, the ABP velocity is affected by wall interactions (cf. Eq. (26)). The effect vanishes with increasing concentration, since the persistence length decreases due to enhanced particle collisions at higher concentrations. The local swim stress in $\Delta V$ (19) can be separated into two contributions by inserting the athermal and overdamped equation of motion (Eq. (1)), namely

[TABLE]

Hence, the reduction of the relative local swim stress in Fig. 5(a) is a consequence of the correlations between the interparticle forces and the propulsion direction. The effect increases with increasing concentration. We like to emphasize that density variations in the bulk are negligible small for the considered wall separations in Fig. 5(b). In fact, surface accumulation is lower in a system of self-avoiding ABPs compared to an ideal ABP gas. The global swim stress (10) includes, in addition to the terms in Eq. (32), a wall contribution as shown in Ref.  25. It is the latter part, which implies a system-size dependence of the global swim stress (open circles in Fig. 5(b)). The effect gradually vanishes with increasing system size. In contrast, the global interparticle-force stress is independent of system size, which can be expected for a system at constant average density $\rho$. However, the local interparticle-force stress depends on $L_{z}$. Since the total local and global stresses are equal, the system-size dependence of the local interparticle-force stress and that of the global swim stress are equal.

Figure 5(b) clearly reveals differences between the global swim (Eq. (9)) and local swim stress (Eq. (17)), specifically for finite system sizes. However, in the thermodynamic limit of infinitely large systems, the two definitions yield identical results. The total stress in a volume $V$ (global) or $\Delta V$ (local) is not affected by the difference, since other contributions by interparticle forces yield identical differences and, hence, compensate for the disparity in the swim stresses.

Local stress in inhomogeneous systems of ABPs

The stress in phase-separated, inhomogeneous ABP systems has been calculated before by various methods [48, 49]. The calculation of the global stress in Ref.  48 is performed with a particle-based model (ABPs) in a system with periodic boundary conditions. However, no stress normal to the interface has been obtained. In Ref.  49, a continuum description based on a generalization of the Cahn-Hilliard equation is applied. Here, we illustrate the suitability of our particle-based expression for the calculation of local stresses for a phase-separated system of ABPs confined between two walls. The same geometry is considered as before, but with the dimensions $L_{x}=L_{y}=25\sigma$. The Péclet number $Pe=80$ is sufficiently large to yield phase separation into crystalline-like layers adjacent to walls and a fluid phase in the center, as is displayed in Fig. 6(a) and (b). In order to calculate the local stress, we consider local volumes of width $\Delta L_{z}=4\sigma$, i.e., we average over four crystalline layers. As shown in Fig. 6(c), the stress normal to the walls and the interface ($z$ direction) is constant over the whole system. Moreover, the local stress agrees with the stress exerted on the wall. In the various regimes, stress is dominated either by the swim stress (Eq. (19)) (central fluid part) or the force contribution to stress (Eq. (18)) (crystalline wall part). Specifically in the crossover regime, where both stresses contribute, the local stress is equal to the global stress (wall stress). This result confirms that the derived expression is suitable for calculating local stress even in inhomogeneous systems.

Summary and conclusions

We have derived expressions for the stress (pressure) of a system of confined active Brownian particles, specifically, generally valid expressions for the local stress in a subvolume $\Delta V(\bm{r})$ centered around a point $\bm{r}$. In addition to previously derived contributions to the stress of active systems confined between walls [25], the contribution of the finite-range surface interaction is taken into account in the calculation of the global stress. For the stress in a local volume ($\Delta V$), we find two expressions by following the strategy of Ref.  36. On the one hand, stress is given by momentum transfer across a hypothetical plane, including force as well as kinetic (momentum) contributions. Here, we find a motility contribution by introducing an active momentum in analogy to the kinetic momentum of a particle of mass $m$. On the other hand, we identify a virial expression for the active stress (local swim stress) of particles inside of $\Delta V$, different from the known swim stress [22, 41, 25]. In fact, in Ref.  35 a formally similar expression has been defined, however, with an active impulse of very different meaning. Computer simulations show that our obtained local stress expression agrees quantitatively with the corresponding global stress expression of an ABP system confined between walls for both, non-interacting ABPs as well as ABPs interacting via a Lennard-Jones potential. This underlines the existence of an equation of state for a system of spherical ABPs.

Our results reach beyond previous considerations and discussions on the stress in active systems. Initially, swim pressure has been introduce via Clausius’ virial [50], taking active forces into account [22, 23]. Based on the equations of motion of ABPs for an infinite system, the swim pressure of an ideal active gas has then been proposed [22]. Subsequently, derivations of active pressure expressions have been presented for ABPs confined between walls [25, 26, 27] or in periodic systems [25, 27]. Alternatively, pressure expressions have been determined via the Fokker-Planck equation of an ABP system [28], or the equations of motion for the density field [35]. A derivation of an expression for a local pressure has be attempted in Ref.  27, with the conclusion that a local pressures exists, but the bulk value differs from the mechanical pressure exerted on confining walls. In case of an ideal gas of ABPs, the bulk pressure would be given by $\rho(\bm{r})k_{B}T$, with the position dependent density $\rho(\bm{r})$ and, hence, would vanish in an athermal system. Hence, active pressure is considered as a boundary effect [27]. In contrast, the present calculations provide a general expression for the local stress and pressure and thus, mark a substantial advancement. The derived local stress/pressure expressions are more general and reduce to the known expressions in the limit of the local volume being equal to the overall volume. This is emphasized by the fact that the local swim stress (19) is finite, while the “traditional” swim stress (10) vanishes in $\Delta V$. Most importantly, our calculations demonstrated that active stress/pressure is not a wall (boundary) effect, as stated in Ref.  27, but is caused by active momentum transport in analogy to momentum transfer in a passive system. Hence, our results allow us to resolve apparent contradictions in the interpretation of local stress and consequently in the existence of equations of state in active systems.

As we demonstrated, local stress or pressure can be calculated in inhomogeneous ABP systems. Extensions to local shear stresses are straightforward and, thus, surface tensions of active systems can be determined.

Methods

We perform simulations in three dimensions considering the overdamped Eq. (1), $m=0$, without thermal fluctuations. The translational equations of motion (1) are solved via the Ermak-McCammon algorithm [51]. The equations of motion (3) for the orientation vectors are solved using the scheme described in Ref.  25.

Pair-wise interparticle interactions are taken into account by the repulsive Lennard-Jones potential

[TABLE]

where $r_{ij}=|\bm{r}_{ij}|$ and the interaction strength is set to $\epsilon/k_{B}T=Pe$; $Pe$ is Péclet number $Pe$ defined as

[TABLE]

This yields the same particle overlap independent of the Péclet number [52]. The ABPs are confined between two walls parallel to the $xy$-plane of the Cartesian reference frame located at $S_{iz}=\pm L_{z}/2$ (cf. Fig. 1). Parallel to the walls, periodic boundary conditions are applied. The purely repulsive ABP-wall interaction is described by the Lennard-Jones potential of Eq. (35) with $r_{ij}$ replaced by $z_{i}-S_{iz}$.

We determine the stress $\sigma_{zz}$ along the $z$-axis in two ways. Firstly, we calculate the global stress of the whole system according to Eq. (9) ($m=0$). Secondly, the local stress in a volume $\Delta V$ is calculated via Eq. (17). Thereby, the box dimension $L_{z}$ normal to the walls is chosen such that it is significantly larger than the persistence length $l_{p}=v_{0}/2D_{R}=Pe\sigma/2$ of an ABP, in particular, we consider wall separations $L_{z}\geqslant 5l_{p}$. The volume $\Delta V$ of width $\Delta L_{z}$ is located in the bulk of the system along the $z$-axis (cf. Fig. 1). If not indicated otherwise, the dimensions of the volume $V$ are $L_{x}=L_{y}=40\sigma$, $L_{z}=100\sigma$, and $\Delta L_{z}=0.2L_{z}$. In Fig. 1, $r_{z}$ denotes the position of the center of the slit within the interval $L_{z}/2<r_{z}<L_{z}/2$. The number density of ABPs in $V$ is $\rho\sigma^{3}=0.3$.

Acknowledgements

Financial support by the Deutsche Forschungsgemeinschaft (DFG) within the priority program SPP 1726 “Microswimmers—from Single Particle Motion to Collective Behaviour” is gratefully acknowledged.