Papers
Topics
Authors
Recent
Search
2000 character limit reached

Active Brownian Particle (ABP) Overview

Updated 24 November 2025
  • Active Brownian Particle (ABP) is a stochastic model for self-propelled particles exhibiting persistent ballistic motion disrupted by rotational diffusion.
  • ABP dynamics are modeled by overdamped Langevin equations that couple translation and orientation, facilitating studies of diffusion transitions and boundary effects.
  • The model reveals rich collective phenomena including motility-induced phase separation and glassy states, with implications for active matter research.

An active Brownian particle (ABP) is a paradigmatic stochastic model for self-propelled microswimmers, such as motile bacteria, synthetic Janus colloids, or cytoskeletal molecular motors. ABPs move via persistent ballistic motion disrupted by angular (rotational) diffusion, yielding non-equilibrium dynamics fundamentally distinct from classical Brownian particles. The ABP model provides a minimal yet versatile framework for the study of active matter, elucidating single-particle crosses from ballistic to effective diffusion, collective behaviors such as motility-induced phase separation, glassy states, and interactions with complex boundaries or environments.

1. Microscopic Dynamics and Generalized Langevin Description

A canonical ABP in d=2,3d=2,3 dimensions is governed by overdamped Langevin equations coupling translational and orientational degrees of freedom. In two dimensions, for particle ii with position ri\mathbf{r}_i and orientation θi\theta_i, the dynamics reads (Zhang et al., 2022, Basu et al., 2018):

r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}

where

  • v0v_0 is the self-propulsion speed,
  • ni=(cosθi,sinθi)\mathbf{n}_i = (\cos\theta_i, \sin\theta_i) defines the propulsion direction,
  • DtD_t (DrD_r): translational (rotational) diffusion coefficient,
  • ξ\xi, ii0 are independent Gaussian white noises,
  • ii1 encodes interparticle or external forces.

The persistence time is ii2; the persistence length ii3. For an isolated ABP (ii4), the mean-squared displacement (MSD) exhibits a crossover from ballistic to diffusive scaling. Explicit results for ii5 (Zhang et al., 2022, Basu et al., 2018):

ii6

At ii7, MSD ii8 (ballistic), while for ii9, rotational decorrelation yields ri\mathbf{r}_i0 with ri\mathbf{r}_i1. In three dimensions ri\mathbf{r}_i2 (Zhang et al., 2022, Wang et al., 2019).

The Peclet number, ri\mathbf{r}_i3, and ri\mathbf{r}_i4 control the balance of persistence and diffusion.

2. Short-Time and Large-Deviation Statistics

Short-time statistics of ABP positions, particularly for ri\mathbf{r}_i5, are governed by the colored, memory-carrying active noise (Majumdar et al., 2020, Basu et al., 2018). The resulting distributions of ri\mathbf{r}_i6, ri\mathbf{r}_i7 are strongly anisotropic and non-Gaussian:

  • ri\mathbf{r}_i8 is sharply peaked near ri\mathbf{r}_i9 with non-Gaussian tails,
  • θi\theta_i0 spreads as a random acceleration process, yielding broader (but still non-Brownian) marginals and anomalous first-passage exponents.

In the absence of translational diffusion, the ABP's position distribution has compact support, θi\theta_i1. Optimal fluctuation theory provides explicit rate functions for large deviations, revealing essential singularities as the boundary of this support is approached (Majumdar et al., 2020). These signatures of activity remain robust beyond the Markovian regime.

3. Confinement Effects and Boundary Accumulation

When subject to external potentials or geometric confinement, ABPs display steady-state distributions markedly distinct from passive particles (Zhang et al., 2022, Wang et al., 2019, Wagner et al., 2021). In harmonic traps,

  • For weak activity or low persistence (θi\theta_i2), the stationary distribution is approximately Boltzmann with θi\theta_i3.
  • As persistence grows, the ABP distribution develops pronounced non-Boltzmann features: spatial plateaus, boundary accumulation, and ultimately bimodal structure reflecting the run length θi\theta_i4 (Wang et al., 2019).

Interactions with hard or soft walls induce boundary accumulation and persistent currents. In the zero-noise, steady-state Smoluchowski equation limit, non-equilibrium fluxes are controlled by the persistence length θi\theta_i5 and the nature of the boundary condition. The normal flux at the boundary acts as a source for large-scale density and current patterns, inducing depletion or long-range interactions between obstacles ("active depletion") (Wagner et al., 2021). For boundaries with inhomogeneous or patterned activity profiles, the location-dependent propulsion speed induces steady-state density patterns that "mimic" the substrate shape; typically, the density is low in regions of high activity and peaks at sharp activity gradients (Mishra et al., 2022).

4. Collective Behavior: MIPS, Glassy States, and Mode-Coupling Theory

At finite density, ABPs interact—commonly via short-range repulsion—generating a rich spectrum of collective phenomena:

  • Motility-Induced Phase Separation (MIPS): Above a threshold θi\theta_i6 and density, persistent self-propulsion leads to spontaneous liquid-gas coexistence, even in the absence of attractions. Phase boundaries can be derived from spinodal and binodal criteria based on coarse-grained field theory; e.g., effective diffusivity θi\theta_i7 vanishes at the spinodal θi\theta_i8 (Perez-Bastías et al., 17 Apr 2025).
  • Active Glass Transition: At high density and low persistence, ABPs can form nonergodic glassy states. Mode-Coupling Theory (MCT) captures the kinetic slowing-down and critical behavior. The θi\theta_i9-relaxation time diverges as

r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}0

and varies nontrivially with both r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}1 and r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}2 (Liluashvili et al., 2017). The active glass is softer than its passive analog and retains memory of initial orientational fluctuations within positional correlations, encoding the non-equilibrium character of the arrested state.

  • Effect of Confinement and Disorder: Quenched disorder, random obstacles, or porous environments localize ABPs, suppress large-scale phase separation, and enhance dynamical heterogeneity. MIPS may be suppressed, pinned, or rendered highly spatially heterogeneous under strong confinement (Moore et al., 2023).

5. ABPs in Complex and Viscoelastic Media: Anomalous Transport

ABPs embedded in complex environments such as polymer solutions or cross-linked gels experience anomalous subdiffusive dynamics due to viscoelastic feedback:

  • In star-polymer environments, an ABP exhibits subdiffusive mean-squared displacement scaling r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}3 with r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}4, decreasing as the Péclet number increases. For moderate r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}5, the motion matches fractional Brownian motion with Hurst exponent r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}6; for large r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}7, logarithmic subdiffusion and anti-correlated velocities emerge, exactly captured by a fractional Langevin equation with both thermal and athermal noise (Joo et al., 2020).
  • For ABPs in polymer solutions, the long-time diffusion coefficient r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}8 can be non-monotonic in particle size—larger particles may diffuse faster at high activity—a direct result of the competition between enhanced persistence length and the scale-dependent viscosity of the polymer mesh (Du et al., 2018). The active particle experiences an effective viscosity r˙i=v0ni+2Dtξi(t)+Fi/ζ θ˙i=2Drηi(t)\begin{aligned} & \dot{\mathbf{r}}_i = v_0 \, \mathbf{n}_i + \sqrt{2 D_t} \, \boldsymbol{\xi}_i(t) + \mathbf{F}_i/\zeta \ & \dot{\theta}_i = \sqrt{2 D_r} \, \eta_i(t) \end{aligned}9 larger than that of a passive probe under identical conditions.

6. Extensions: Anisotropic Particles, Noise-Induced Drifts, and Predictive Control

The minimal ABP model can be extended to include anisotropy and nonthermal noise:

  • Fluctuating point-applied propulsion leads to translation-rotation coupling, anisotropic diffusion tensors, and a noise-induced drift term that can persist even in the overdamped regime. These effects alter long-time diffusivity and steady-state properties, and can dominate over thermal contributions depending on particle geometry and noise statistics (Thiffeault et al., 2021).
  • In the context of engineering and optimal control, ABP systems can be steered using model-predictive control frameworks and deep-learning-based density forecasting, leveraging the multi-scale linkage between microscopic dynamics, continuum advection-diffusion PDEs, and control protocols acting through modulation of environmental or particle parameters. These approaches bridge single-particle physics to collective, pattern-forming behaviors in real time (Saremi et al., 7 Sep 2025).

7. First Passage and Survival Phenomena

ABPs exhibit fundamentally altered first-passage dynamics compared to passive particles, both in 1D and higher dimensions:

  • Asymptotic expansions for the mean first passage time (MFPT) in a 1D interval reveal leading v0v_00 corrections that can either decrease or increase MFPT depending on initial position and orientation bias, with symmetry-breaking appearing for non-uniform initial orientations. For unbiased orientations, the first non-trivial correction is v0v_01 and is always a universal reduction in MFPT (Iyaniwura et al., 2023).
  • In higher dimensions, partially absorbing or sticky walls induce non-trivial angular first-passage processes for the orientation, and the MFPT is linked to splitting probabilities and angular sojourn statistics. Non-Markovian (encounter-based) killing can be incorporated through higher-order Fokker-Planck equations for joint position–orientation–occupation-time statistics (Bressloff, 2023).

References

Definition Search Book Streamline Icon: https://streamlinehq.com
References (17)

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Active Brownian Particle (ABP).