Active Brownian particle

From HandWiki
Revision as of 17:28, 6 February 2024 by Dennis Ross (talk | contribs) (update)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

An active Brownian particle (ABP) is a model of self-propelled motion in a dissipative environment.[1][2] It is a nonequilibrium generalization of a Brownian particle. The self-propulsion results from a force that acts on the particle's center of mass and points in the direction of an intrinsic body axis (the particle orientation).[2] It is common to treat particles as spheres, though other shapes (such as rods) have also been studied.[3][4] Both the center of mass and the direction of the propulsive force are subjected to white noise, which contributes a diffusive component to the overall dynamics. In its simplest version, the dynamics is overdamped and the propulsive force has constant magnitude, so that the magnitude of the velocity is likewise constant (speed-up to terminal velocity is instantaneous).

The term active Brownian particle usually refers to this simple model and its straightforward extensions, though some authors have used it for more general self-propelled particle models.[4][5]

Equations of motion

Mathematically, an active Brownian particle is described by its center of mass coordinates [math]\displaystyle{ \mathbf{r} }[/math] and a unit vector [math]\displaystyle{ \hat{\mathbf{n}} }[/math] giving the orientation. In two dimensions, the orientation vector can be parameterized by the 2D polar angle [math]\displaystyle{ \theta }[/math], so that [math]\displaystyle{ \hat{\mathbf{n}} = (\cos \theta, \sin \theta) }[/math]. The equations of motion in this case are the following stochastic differential equations:

[math]\displaystyle{ \begin{align} \dot{\mathbf{r}} &= v_0 \hat{\mathbf{n}} - (m \xi)^{-1} \nabla V(\mathbf{r}) + \sqrt{2 D_t} \, \boldsymbol{\eta}_{\text{trans}}(t) \\ \dot{\theta} &= \sqrt{2 D_r} \, \eta_{\text{rot}}(t). \end{align} }[/math]

where

[math]\displaystyle{ \begin{align} \langle \eta_{\text{rot}}(t)\rangle &= 0; \qquad \langle \eta_{\text{rot}}(t) \eta_{\text{rot}}(t')\rangle = \delta(t - t') \\ \langle \boldsymbol{\eta}_{\text{trans}}(t)\rangle &= \boldsymbol{0}; \qquad \langle \boldsymbol{\eta}_{\text{trans}}(t) \boldsymbol{\eta}^{\intercal}_{\text{trans}}(t') \rangle = \mathbf{I} \delta(t-t') \end{align} }[/math]

with [math]\displaystyle{ \mathbf{I} }[/math] the 2×2 identity matrix. The terms [math]\displaystyle{ \boldsymbol{\eta}_{\text{trans}}(t) }[/math] and [math]\displaystyle{ \eta_{\text{rot}}(t) }[/math] are translational and rotational white noise, which is understood as a heuristic representation of the Wiener process. Finally, [math]\displaystyle{ V(\mathbf{r}) }[/math] is an external potential, [math]\displaystyle{ m }[/math] is the mass, [math]\displaystyle{ \xi }[/math] is the friction, [math]\displaystyle{ v_0 }[/math] is the magnitude of the self-propulsion velocity, and [math]\displaystyle{ D_t }[/math] and [math]\displaystyle{ D_r }[/math] are the translational and rotational diffusion coefficients.[6]

The dynamics can also be described in terms of a probability density function [math]\displaystyle{ f(\mathbf{r},\theta, t) }[/math], which gives the probability, at time [math]\displaystyle{ t }[/math], of finding a particle at position [math]\displaystyle{ \mathbf{r} }[/math] and with orientation [math]\displaystyle{ \theta }[/math]. By averaging over the stochastic trajectories from the equations of motion, [math]\displaystyle{ f(\mathbf{r},\theta, t) }[/math] can be shown to obey the following partial differential equation:

[math]\displaystyle{ \frac{\partial f}{\partial t} + v_0 \hat{n} \cdot \nabla f = (m \xi)^{-1} \nabla \cdot (\nabla V(\mathbf{r}) \, f) + D_r \frac{\partial^2 f}{\partial \theta^2} + D_t \nabla^2 f }[/math]

Behavior

For an isolated particle far from boundaries, the combination of diffusion and self-propulsion produces a stochastic (fluctuating) trajectory that appears ballistic over short length scales and diffusive over large length scales. The transition from ballistic to diffusive motion is defined by a characteristic length [math]\displaystyle{ \ell = v_0/D_r }[/math], called the persistence length.[1]

In the presence of boundaries or other particles, more complex behavior is possible. Even in the absence of attractive forces, particles tend to accumulate at boundaries. Obstacles placed within a bath of active Brownian particles can induce long-range density variations and nonzero currents in steady state.[7][8]

Sufficiently concentrated suspensions of active Brownian particles phase separate into a dense and dilute regions.[9][10] The particles' motility drives a positive feedback loop, in which particles collide and hinder each other's motion, leading to further collisions and particle accumulation.[1] At a coarse-grained level, a particle's effective self-propulsion velocity decreases with increased density, which promotes clustering. In the more general context of self-propelled particle models, this behavior is known as motility-induced phase separation.[9] It is a type of athermal phase separation because it occurs even if the particles are spheres with hard-core (purely repulsive) interactions.

Variations

A variant of active Brownian motion involves complete directional reversals in addition to rotational diffusion. This movement pattern is seen in bacteria like Myxococcus xanthus, Pseudomonas putida, Pseudoalteromonas haloplanktis, Shewanella putrefaciens, and Pseudomonas citronellolis.[11]

See also

Notes

Sources