Active Brownian Particles

Updated 15 February 2026

Active Brownian Particles are self-propelled agents converting energy into directed motion with persistent propulsion and rotational diffusion.
They display unique non-equilibrium phenomena such as motility-induced phase separation, enhancing diffusivity and leading to dense-dilute phase coexistence.
Rigorous models employing Langevin dynamics and continuum theories enable exploration of structural ordering, glassy dynamics, and boundary effects in both synthetic and biological systems.

Active Brownian Particles (ABPs) are a fundamental class of self-propelled, non-equilibrium agents that convert energy into directed motion, with their orientation subject to Brownian rotational diffusion. Distinguished from classical Brownian particles by a persistent propulsion driven along an internal axis, ABPs serve as the canonical model for a broad range of synthetic and biological active matter systems, including colloidal Janus particles, motile bacteria, and chemically-powered microswimmers. The field has seen the development of a robust mathematical framework encompassing stochastic particle models, kinetic and hydrodynamic continuum theories, and systematic thermodynamic approaches.

1. Mathematical Formulation and Stochastic Dynamics

The microscopic dynamics of an ABP are governed by overdamped Langevin equations coupling self-propulsion and rotational diffusion. In $d$ dimensions, the position $\mathbf{r}_i$ and orientation unit vector $\mathbf{p}_i$ of particle $i$ evolve as

$\dot{\mathbf{r}}_i = v_0 \mathbf{p}_i + \sqrt{2D_t}\, \boldsymbol{\eta}_i (t)$

$\dot{\mathbf{p}}_i = \mathbf{p}_i \times \boldsymbol{\xi}_i(t)$

where $v_0$ is the propulsion velocity, $D_t$ is the translational diffusion coefficient, and $\boldsymbol{\eta}_i(t)$ , $\boldsymbol{\xi}_i(t)$ are independent, zero-mean Gaussian white noises with $\langle\eta_{i\alpha}(t)\eta_{j\beta}(t')\rangle = \delta_{ij}\delta_{\alpha\beta}\delta(t-t')$ and similarly for $\xi$ (Zhang et al., 2022). The orientation dynamics describe rotational diffusion on the unit sphere, with persistence time $\tau_r=1/(2D_r)$ or $1/D_r$ in two dimensions.

This minimal model extends to interacting systems by the inclusion of pairwise potentials (e.g., Weeks–Chandler–Andersen) or hard-sphere exclusions, yielding

$\dot{\mathbf{r}}_i = v_0 \mathbf{p}_i + \frac{1}{\gamma}\sum_{j\neq i} \mathbf{F}_{ij} + \sqrt{2D_t}\, \boldsymbol{\eta}_i(t)$

and corresponding orientation equations (Stenhammar et al., 2013, Reichert et al., 2020).

ABPs differ fundamentally from passive Brownian particles: their self-propulsion violates detailed balance at the single-particle level, with the distribution of displacement non-Gaussian at short times, and persistent motion yielding enhanced long-time diffusivity. For a single ABP in 2D, the mean-squared displacement is

$\langle |\Delta \mathbf{r}(t)|^2 \rangle = 4D_t t + \frac{2v_0^2}{D_r^2}\left[e^{-D_r t} -1 + D_r t\right]$

yielding an effective diffusion constant $D_\text{eff} = D_t + v_0^2/2D_r$ (Zhang et al., 2022, Wang et al., 2019).

2. Collective Phenomena: Motility-Induced Phase Separation and Nonequilibrium Thermodynamics

ABPs exhibit emergent phenomena absent in equilibrium systems. Purely repulsive ABPs can undergo motility-induced phase separation (MIPS), where propulsion creates effective attractions, leading to coexistence of dense and dilute phases reminiscent of a gas–liquid transition (Stenhammar et al., 2013, Yang et al., 2022). The mechanism stems from crowding-induced slowdowns: as local density increases, persistent self-propelled particles become caged, favoring further accumulation.

Continuum theories derived via coarse-graining lead to non-equilibrium Cahn–Hilliard-type equations for the density $\rho(\mathbf{r},t)$ : $\partial_t \rho = -\nabla \cdot \left[ -D(\rho)\rho \nabla \mu(\rho) + \sqrt{2D(\rho)\rho}\, \boldsymbol{\Lambda} \right]$ with an effective chemical potential $\mu(\rho) = \ln \rho + \ln v(\rho)$ and mobility $D(\rho)$ (Stenhammar et al., 2013). The inhomogeneous stationary state violates detailed balance at gradient level, with coarse-grained free energy functionals including non-equilibrium contributions (Chakraborti et al., 2016).

Key features:

Spinodal and binodal densities: Phase coexistence controlled by the sign of $f''(\rho)$ of the effective free-energy density, with critical activity $P_c$ and density $\rho_c$ closely matching simulations.
Equation of state and large deviations: An additivity principle enables construction of non-equilibrium analogs of chemical potential, compressibility, and free energy, yielding accurate predictions for subsystem fluctuations (Chakraborti et al., 2016).
Coarsening kinetics: Domain sizes follow $L(t)\sim t^\alpha$ with $\alpha \approx 0.27$ , slightly subdiffusive compared to passive coarsening (Stenhammar et al., 2013).

Recent work generalizes this framework to coupled density–polarization “two-field” models, enabling explicit description of wetting, interfacial phenomena, and polarization heterogeneity at boundaries (Perez-Bastías et al., 17 Apr 2025).

3. Structure, Dynamics, and Velocity Correlations

Even in the absence of aligning or attractive interactions, ABPs display nontrivial order and spatiotemporal coherence:

Structural ordering: In the dense phase, local hexagonal (crystalline-like) order is quantified using bond-orientational parameters such as $\psi_i = (1/N_i) \sum_{j \in \mathrm{nn}(i)} \exp[i 6 \theta_{ij}]$ (Yang et al., 2022).
Velocity correlations and coherence: Although thermal noise dominates instantaneously, averaging particle velocities over an intermediate lag time ( $\Delta t^*\sim 0.075\tau$ ) reveals regions of significant velocity alignment, spatially scaling with structural clusters and persisting over the cage retention timescale. This coherence is accessible only via appropriate temporal averaging, which cancels the uncorrelated thermal contributions and reveals the persistent, correlated component of motion (Yang et al., 2022).
Glassy dynamics: At high densities, ABPs display dynamical arrest analogous to glass transitions. Mode-coupling theory for ABPs gives a unified approach for passive/active tracer relaxation, correctly predicting fluidization by activity and distinctive non-equilibrium signatures such as time-reversal asymmetry in angular correlations (Reichert et al., 2020).

4. Effects of Confinement, Boundaries, and Substrate Inhomogeneity

Boundary conditions and substrate patterning have pronounced effects on ABP behavior:

Wall attraction, absorption, and sticky boundaries: ABPs accumulate at confining walls due to persistent propulsion. Explicit treatment of partially absorbing or sticky boundaries yields analytical first-passage and mean trapping-time statistics as combined angular diffusion and bulk return processes, with the mean absorption time controlled by motility, wall geometry, and absorption rate (Bressloff, 2023).
Morphological transitions at porous walls: ABPs confined by porous boundaries undergo wetting–dewetting transitions from connected dense layers to droplet states as wall porosity or activity increases. This active wetting can be interpreted via an effective wall–liquid interfacial tension, which changes sign with control parameters (Das et al., 2020).
Mimicking substrate patterns: ABPs on substrates with spatially heterogeneous activity fields develop steady-state density and polarization profiles that directly mirror the underlying substrate inhomogeneity. Notably, regions of high activity are depleted in density, with accumulation at interfaces, confirming the capacity of ABPs to “read out” patterns embedded in motility landscapes purely via their stochastic dynamics (Mishra et al., 2022).
Chemotactic/anti-chemotactic behavior: In time-stationary confined gradients, ABPs exhibit “anti-chemotaxis,” accumulating away from fuel sources. However, in non-stationary, bursty fuel environments, their distribution transiently shifts toward sources, enabling functional chemotactic behavior even without explicit gradient sensing (Merlitz et al., 2019).

5. Advanced Theoretical Developments: Field Theory, Anisotropy, and Control

Contemporary directions extend ABP theory in multiple dimensions:

Field-theoretic approaches: Janssen–De Dominicis path integral methods cast the ABP process into a field theory over density and polarization variables, allowing systematic perturbation theory and RG analysis. This yields exact results for mean-square displacement and density in potentials, with explicit nonequilibrium corrections (Zhang et al., 2022).
Anisotropic and fluctuating-propulsion ABPs: Modeling of ABPs subject to nonthermal, off-center, or anisotropic propulsion extends conventional models. These generalizations generate novel noise-induced drifts and translational–rotational coupling in the diffusion tensor, which persist even in the overdamped limit and can dominate swimmer motion in practical settings (Thiffeault et al., 2021).
Control and machine learning of ABP collectives: Multiscale approaches integrate particle-based simulations, continuum spectral PDE solvers, PDE-constrained model predictive control, and neural surrogates (Conv1D-LSTM-attention hybrids) to forecast and steer ABP density fields in real time. This allows high-fidelity tracking of prescribed patterns and rapid evaluation of control policies (Saremi et al., 7 Sep 2025).
Kinetic theory at moderate densities: Boltzmann–Lorentz kinetic equations quantify the reduction of effective swim speed and diffusion due to collisions, predicting a density-dependent effective diffusion coefficient $D(\rho) = V_{\rm eff}^2(\rho)/(2 D_r)$ up to area fractions $\eta \lesssim 0.42$ (Soto, 2 Jan 2025).

6. Applications, Benchmarks, and Experimental Connections

ABP models have led to several demonstrable technological and biological applications:

Size-selective demixing and cargo delivery: Field-driven methods employing time-dependent activity waves can achieve separation of ABP mixtures by size, as the resonance conditions (optimal wave velocity and spatial frequency) vary strongly with particle diameter. This enables reversible, fully controlled demixing and targeted delivery in active colloidal suspensions (Merlitz et al., 2018).
Contagion modeling: By embedding infection and recovery dynamics in ABP simulations, one derives first-principles expressions for contagion rates, critical densities, and spatiotemporal epidemic patterns without free parameters. This approach elucidates the microscopic origins of SI/SIR kinetics and reveals features (e.g., symmetric epidemic peaks) absent from classic mean-field models (Norambuena et al., 2020).
Validation via simulation and experiment: Extensive MD, Brownian, and event-driven simulations confirm continuum and thermodynamic theory predictions, with measured binodals, coarsening exponents, and subsystem distributions in close agreement. Quantitative tests of ABP diffusion, phase coexistence, and wetting transitions are feasible with colloidal Janus particles and other microswimmer platforms (Stenhammar et al., 2013, Reichert et al., 2020, Das et al., 2020, Wang et al., 2019).

Active Brownian Particles form a versatile and mathematically tractable paradigm for non-equilibrium statistical mechanics, combining rich phenomenology (phase separation, glassiness, pattern formation, chemotaxis-like competition) with rigorous kinetic, hydrodynamic, and thermodynamic frameworks. Their theoretical and computational tractability, combined with experimental realizability, continues to drive progress in the broader field of active matter.