---
title: Run-and-Tumble Particles (RTPs) Overview
url: https://www.emergentmind.com/topics/run-and-tumble-particles-rtps
type: topic
---

# Run-and-Tumble Particles (RTPs) Overview

Searching arXiv for recent and foundational papers on run-and-tumble particles to ground the article in published work.
Run-and-tumble particles (RTPs) are active particles that self-propel at fixed speed along an internal orientation during persistent runs and then undergo stochastic reorientation events, or tumbles, that reset or randomize that orientation. In the minimal kinetic description, an RTP moves with propulsion velocity \(v{\bf u}\), with \({\bf u}\) reoriented by a Poisson process of rate \(\alpha\), or by analogous lattice updates in discrete-state models [1206.1805]. Across continuum, lattice, interacting, and confined settings, RTPs serve as a canonical nonequilibrium model for persistent transport, motility-induced phase separation (MIPS), anomalous collective diffusion, boundary accumulation, jamming, and rare-event kinetics. A recurrent theme in the literature is that RTPs are microscopically distinct from active Brownian particles (ABPs), yet often share the same leading coarse-grained hydrodynamics in isotropic torque-free settings, while departing sharply under strong confinement, anisotropy, or explicit internal-state structure [1504.07391].

## 1. Definition, microscopic dynamics, and canonical formulations

In the standard continuum formulation, the one-particle probability density \(\psi({\bf r},{\bf u},t)\) obeys a kinetic equation with self-propulsive advection, optional translational diffusion \(D_t\), optional rotational diffusion \(D_r\), and tumble loss/gain terms. For a pure RTP one sets \(D_r=0\), leaving tumbles at rate \(\alpha\) as the sole orientational decorrelation mechanism [1206.1805]. In one dimension, the basic microscopic law is \(dX/d\tau=v\sigma(\tau)\), where \(\sigma(\tau)\in\{-1,+1\}\) flips as a Poisson process, while in two dimensions the orientation is commonly parameterized by an angle \(\theta\) with \(\hat n=(\cos\theta,\sin\theta)\) and tumbles reset \(\theta\) uniformly on \([0,2\pi]\) [2012.08276; 2211.07337].

Lattice RTPs replace continuous orientations by discrete nearest-neighbor headings. On a two-dimensional square lattice, a representative model assigns each particle an internal orientation \({\bf \theta}_m\in\{0,\pi/2,\pi,3\pi/2\}\), enforces hard-core exclusion \(n_{\bf i}=0,1\), and combines oriented hopping, lateral hopping, and tumbling [2406.11726]. In that model, a particle hops along its current orientation with unit rate, to each of the other nearest neighbors with rate \(p\), and tumbles with rate \(\omega\). The parameter \(p\) acts as positional diffusion, and the paper emphasizes that this ingredient is essential for stable MIPS in that lattice family [2406.11726].

A useful generalization is the one-dimensional PDMP framework, where the state is \((x_t,\sigma_t)\in\mathbb R\times\Sigma\), deterministic motion is \(\partial_t x_t=v_\sigma(x_t)\), and the internal state jumps with rates \(\lambda_{\sigma\tilde\sigma}(x)\). This formulation covers RTPs in external potentials, finite-state internal dynamics, and resetting, and makes the invariant-measure problem amenable to local regularity analysis [2506.02274]. Another exact continuous-space PDMP treatment studies two interacting RTPs on a one-dimensional ring by reducing the problem to the interparticle distance \(r_t\in[0,L/2]\) plus an internal pair state, with jamming appearing as a boundary condition in configuration space [2306.00831].

These model classes share the same physical core: bounded-speed self-propulsion, orientational persistence, stochastic reorientation, and nonequilibrium steady states even in simple confining geometries. What varies across formulations is the representation of orientation, the role of exclusion and interaction, and whether one studies single-particle transport, few-body jamming, or collective phases.

## 2. Coarse-graining, hydrodynamics, and relation to active Brownian particles

The central coarse-grained result for isotropic RTPs is that orientational moments relax on a finite persistence time and can be adiabatically eliminated at long times and large scales. Decomposing \(\psi({\bf r},{\bf u},t)\) into scalar density \(\varphi\), polarization \({\bf p}\), and higher harmonics yields, to leading order in gradients,
\[
p_a=-\frac{1}{(d-1)D_r+\alpha}\nabla_a(v\varphi)+O(\nabla^2),
\]
and therefore a diffusion-drift density equation
\[
\partial_t\varphi=-\nabla\cdot{\bf J},\qquad {\bf J}={\bf V}\varphi-D\nabla\varphi,
\]
with
\[
D=\frac{v^2}{d[(d-1)D_r+\alpha]}+D_t,\qquad
{\bf V}=-\frac{v\nabla v}{d[(d-1)D_r+\alpha]}.
\]
For pure RTPs, \(D_r=0\), so the orientational relaxation rate is simply \(\alpha\) [1206.1805; 1504.07391].

This leads to the well-known ABP/RTP mapping
\[
\alpha \leftrightarrow (d-1)D_r,
\]
valid for isotropic motility parameters that depend on position or density but not on \({\bf u}\), and at leading order in the gradient expansion. Under that mapping, ABPs and RTPs share the same large-scale diffusion-drift equation, the same fluctuating hydrodynamics, and the same local MIPS criterion [1206.1805]. The equivalence is therefore hydrodynamic rather than microscopic: the full angular hierarchies differ term by term, but the leading conserved density dynamics depends only on the combined orientational relaxation rate \((d-1)D_r+\alpha\) [1504.07391].

The same coarse-graining yields a fluctuating density equation of the form
\[
\partial_t\rho=-\nabla\cdot\left({\bf V}[\rho]\rho-D[\rho]\nabla\rho+\sqrt{2D\rho}\,{\boldsymbol\Lambda}\right),
\]
with white noise \(\langle \Lambda_i({\bf r},t)\Lambda_j({\bf r}',t')\rangle=\delta_{ij}\delta({\bf r}-{\bf r}')\delta(t-t')\) [1206.1805]. When \(v\) is density dependent and local, \(v=v(\rho)\), this admits an equilibrium-like free-energy representation at leading order, with
\[
\mu=\ln(\rho v),
\]
and spinodal instability criterion
\[
\frac{dv}{d\rho}<-\frac{v}{\rho},
\]
which is the standard local MIPS condition in both RTP and ABP hydrodynamics [1504.07391].

The same papers are careful about where equivalence fails. Strong external potentials, confinement with size comparable to the run length, higher-order gradient terms, orientational torques, and explicitly \({\bf u}\)-dependent motilities all expose the discrete tumble dynamics and spoil the mapping [1206.1805; 1504.07391]. This distinction between microscopic nonequivalence and hydrodynamic equivalence organizes much of the RTP literature.

## 3. Single-particle transport, confinement, and stationary structure

In weak external potentials, RTPs admit an effective-equilibrium regime. For two-dimensional sedimentation under \(V_{\rm ext}(z)=\delta m g z\), the exact distal steady state is exponential,
\[
\rho(z)=\rho_0 e^{-\lambda z},\qquad
\lambda=\frac{2\alpha v_s}{v^2-v_s^2},
\]
with anisotropic angular distribution
\[
f(\theta)=\frac{1-v_s^2/v^2}{2\pi[1-2(v_s/v)\cos\theta+v_s^2/v^2]}.
\]
Expanding for \(v_s\ll v\) gives an effective temperature \(kT_{\rm eff}=v^2\zeta/(2\alpha)\), but the paper emphasizes that beyond weak forcing no single effective temperature survives and RTP sedimentation is genuinely nonequilibrium [1504.07391].

Strong confinement reveals more distinctively RTP physics. In a harmonic trap, ABPs and RTPs share the same weak-trap Gaussian profile, but for strong persistence \(\Phi=\kappa\tau\gg1\), RTPs travel across the trap along straight chords after tumbles at the boundary, leading to bulk density scaling \(\rho_{\rm bulk}\propto \alpha/\kappa\sim 1/\Phi\), much larger than for ABPs in the same regime [1504.07391]. In nested two-dimensional mazes, geometry alone separates RTPs from ABPs: ABPs escape outward from the center faster, whereas RTPs reach the center from the rim more easily because a tumble can instantly reverse direction and enter an inner opening. In the circular Matryoshka-like maze, ABPs accumulate in the outermost region at large persistence, while RTPs occupy all zones with nearly equal probability [1611.00191].

Several exact one-dimensional and two-dimensional stationary-state problems further illustrate how persistence reshapes invariant measures. For a 1D RTP in an inhomogeneous medium with telegraphic noise \(\sigma(t)\) and position-dependent flip rates
\[
R_1(x)=\left(\frac{|x|}{l}\right)^\alpha[\gamma_1\theta(x)+\gamma_2\theta(-x)],\qquad
R_2(x)=\left(\frac{|x|}{l}\right)^\alpha[\gamma_2\theta(x)+\gamma_1\theta(-x)],
\]
the particle has a normalizable steady state on the infinite line when \(\gamma_1>\gamma_2\), with exact density
\[
P_\alpha^{\rm st}(x)\propto
\exp\!\left[-\frac{\gamma_1-\gamma_2}{v(\alpha+1)l^\alpha}|x|^{\alpha+1}\right].
\]
This gives exponential confinement for \(\alpha=0\) and Gaussian confinement for \(\alpha=1\), generated purely by inhomogeneous tumbling rather than an external potential [2004.11041].

Under stochastic resetting in two dimensions, a Poisson reset at rate \(\alpha\) to the origin with orientation randomization creates a nonequilibrium stationary state with qualitatively different radial and Cartesian marginals. The radial stationary density tends to a constant, \(\mathcal P_\alpha^s(r)\to \alpha/v_0\) as \(r\to0\), whereas the stationary \(x\)-marginal diverges logarithmically,
\[
P_\alpha^s(x)=-\frac{\alpha}{\pi}\ln|x|+O(1),
\]
and both decay at large distance with the same exponent \(\sqrt{\alpha^2+2\alpha\gamma}/v_0\) [2009.09891]. In one-dimensional PDMP language, the regularity of stationary RTP measures under confinement can be characterized without an explicit closed-form density: for a wide class of 1D RTPs in a potential, the invariant density is continuous at high tumble rates and develops divergences at force-balance points at low tumble rates, a phenomenon termed shape transition [2506.02274].

These results collectively show that RTP stationary structure is governed not only by external potentials, but by a broader competition among bounded propulsion, tumbling, inhomogeneous switching, and boundary conditions. A plausible implication is that singular accumulation at force-balance points is a generic nonequilibrium signature of persistent bounded-speed motion, rather than a peculiarity of any single model.

## 4. Interactions, MIPS, clustering, and criticality

For interacting RTPs, the central coarse-grained mechanism behind MIPS is density-suppressed motility. In the hydrodynamic theory with local \(v(\rho)\), particles accumulate where they move more slowly; if \(dv/d\rho<-v/\rho\), this positive feedback destabilizes the homogeneous state and drives phase separation even without attractive forces [1206.1805; 1504.07391]. Microscopic simulations with programmed \(v(\rho)\) and with explicit WCA repulsions show that RTPs and ABPs have closely matching phase diagrams and coexistence densities under the mapping \(\alpha=(d-1)D_r\), supporting the practical usefulness of the hydrodynamic equivalence for MIPS [1504.07391].

A more microscopic collective-diffusion picture emerges in interacting lattice RTPs with hard-core exclusion. On a \(d\)-dimensional periodic lattice, the bulk diffusion coefficient \(D(\rho,\gamma)\) obeys a scaling form
\[
D(\rho,\gamma)=D^{(0)}\mathcal F(\psi),\qquad \psi=\frac{\rho a v}{\gamma},
\]
controlled by the competition between persistence length \(l_p=v/\gamma\) and mean gap \(\langle g\rangle\sim1/(\rho a)\) [2209.11995]. In the strong-persistence regime, the diffusivity varies as \(D\propto \rho^{-\alpha}\) over a broad density range, with \(\alpha\) crossing over from \(0\) at low density to \(2\) at high density, leading to nonlinear diffusion and anomalous density spreading [2209.11995].

On a two-dimensional square lattice with oriented hopping and positional diffusion \(p\), the percolation transition of occupied clusters provides a geometric route to characterize MIPS. At small tumble rate \(\omega\), increasing \(p\) produces a re-entrant sequence
\[
\text{mixed}\to \text{phase-separated/percolating}\to \text{mixed},
\]
because moderate positional diffusion stabilizes macroscopic clustering while too much diffusion destroys it [2406.11726]. At \(\omega=0.020\), Binder-cumulant crossings locate two transitions at \(p_c=0.029\) and \(p_c=0.235\); at \(\omega=0.022\), the corresponding values are \(p_c=0.0414\) and \(p_c=0.191\) [2406.11726].

The critical behavior is unusual. Along the critical line in the \((\omega,p)\) plane, the exponent ratios vary continuously, for example from \((1/\nu,\beta/\nu,\gamma/\nu)=(1.43,0.14,1.72)\) at \((p_c,\omega_c)=(0.235,0.020)\) toward the two-dimensional \(Z_2\)-percolation values near small \((p,\omega)\), while the scaling function \(U_4=F(\xi_2/L)\) remains invariant and matches the Ising-cluster percolation function [2406.11726]. The same group extended this picture to interacting RTPs with explicit nearest-neighbor attraction \(J\): at fixed \(p=0.05\) and \(\rho=1/2\), moderate attraction can suppress a motility-induced phase-separated state, while stronger attraction restores it, producing a reentrant transition in the \((J,\omega)\) plane [2506.05264]. The critical line again shows continuously varying exponents, but Binder-cumulant curves plotted against \(\xi_2/L\) collapse onto the same function as equilibrium interacting percolation, motivating the language of Ising-like super universality [2506.05264].

The geometric viewpoint is important here. These papers distinguish the percolation transition, defined by the emergence of a macroscopic connected dense cluster, from the MIPS transition, defined by density segregation into high- and low-density regions. Their claim is that the two occur at the same critical line, but with exponents related analogously to Ising versus Ising-cluster percolation [2406.11726; 2506.05264]. This suggests that RTP phase separation can be simultaneously understood as a nonequilibrium thermodynamic instability and as a geometric connectivity transition.

## 5. Crowding, jamming, collective transport, and current fluctuations

Crowding changes the meaning of persistence. For a single RTP moving in a lattice of hard-core obstacles, the long-time diffusivity can become nonmonotonic in the tumbling probability \(\alpha\). In the discrete-time model of “Optimized Diffusion of Run-and-Tumble Particles in Crowded Environments,” the mean free run time between obstacle encounters is
\[
\langle \tau_r\rangle=\frac{1}{\rho}
\]
for fixed obstacles, independent of \(\alpha\), by a generalized Kac theorem [1711.05209]. The mean trapping time, however, diverges as \(\alpha\to0\), so maximal persistence is no longer optimal. For static obstacles, the diffusivity has an optimum at \(\alpha_m\propto \rho\), with \(\mathcal D_m\propto1/\rho\), whereas obstacle mobility softens or suppresses this optimum [1711.05209].

At the pair level, jamming can be solved exactly in continuous space. For two RTPs on a one-dimensional ring, the stationary state falls into two universality classes determined by whether probability flows entering and leaving jamming configurations respect a detailed internal-state symmetry [2306.00831]. In the detailed-jamming class, the invariant measure is a uniform bulk density plus a Dirac mass at the jammed boundary \(r=0\); in the global-jamming class, the jammed boundary injects an internal-state skewness that relaxes through the bulk via a matrix ODE, producing exponential or catenary-like stationary profiles [2306.00831]. The paper identifies active global balance as the nonequilibrium analog of ordinary global balance, and interprets jamming as a boundary source in configuration space rather than merely a contact event [2306.00831].

Collective current fluctuations in interacting one-dimensional RTP systems exhibit another layer of hydrodynamic universality. For hardcore RTPs on a ring and for a long-ranged lattice-gas surrogate, the bond-current variance \(\langle Q_i^2(T,L)\rangle\) shows three regimes: model-dependent short-time behavior, a universal intermediate-time subdiffusive growth \(\sim\sqrt{T}\), and a long-time diffusive growth \(\sim T\) [2309.02896]. After scaling by the bulk diffusion coefficient \(D(\rho,\gamma)\), mobility \(\chi(\rho,\gamma)\), and system size \(L\),
\[
\frac{D\,\langle Q_i^2(T,L)\rangle}{2\chi L}=\mathcal W\!\left(\frac{DT}{L^2}\right),
\]
both models collapse onto the same scaling function \(\mathcal W(y)\), with asymptotics \(\mathcal W(y)\sim\sqrt{y/\pi}\) for \(y\ll1\) and \(\mathcal W(y)\sim y\) for \(y\gg1\) [2309.02896]. The same work finds exponentially decaying spatial current correlations with correlation length \(\xi\sim\tau_p^{1/2}\) at large persistence, and a negative temporal current tail \(\mathcal C_0^{JJ}(t)\sim -t^{-3/2}\) underlying the \(\sqrt{T}\) law [2309.02896].

These transport results correct a common misconception that persistence always improves exploration. In free space, decreasing the tumble rate typically increases the diffusivity; in crowded media or under exclusion, persistence also enhances trapping, jamming, or boundary sticking, so optimal transport often occurs at intermediate reorientation rates [1711.05209; 2309.02896].

## 6. Exact asymptotics, branching variants, and rare trajectories

A large body of RTP work concerns exact or asymptotic results for propagators and stationary measures. For a free two-dimensional RTP with tumble rate \(\gamma\), the late-time propagator can be expanded as a series in \((\gamma t)^{-1}\); the leading term is Gaussian with diffusivity
\[
D_{RT}=\frac{v_0^2}{2\gamma},
\]
while each subleading correction satisfies an inhomogeneous diffusion equation whose source depends on lower-order terms [2211.07337]. The position density in the diffusive sector takes the form
\[
P_r(r,t)=\frac{1}{\sqrt{\pi D_{RT}t}}\exp\!\left(-\frac{r^2}{4D_{RT}t}\right)\sum_{k=0}^\infty (\gamma t)^{-k} q_{2k}\!\left(\frac{r}{\sqrt{4D_{RT}t}}\right),
\]
with explicit polynomials \(q_2,q_4,q_6\) computed recursively [2211.07337]. This makes precise how persistence leaves non-Gaussian fingerprints on top of the asymptotic diffusion law.

Branching or source-driven variants produce different long-time structures. In “Run-and-tumble particles on a line with a fertile site,” crossing the origin activates offspring production during the current run, with creation rate modulated by a fertility kernel \(\Theta\) [2012.08276]. The density at the fertile site grows exponentially,
\[
R(t)\sim \rho e^{\chi t},
\]
with \(\chi\) determined solely by the fertility kernel and fertility rate through the renewal equation \(\tilde\Xi(\chi)=1\), while the normalized density converges to a stationary profile with a local minimum at the source,
\[
\mathcal N''(0)=2\chi(\chi+2)>0
\]
[2012.08276]. This is a solvable branching extension of 1D RTP dynamics rather than a conservative transport model.

Rare-event sampling for RTPs has also become technically accessible. “Transition-path sampling for Run-and-Tumble particles” constructs a trajectory-space Monte Carlo scheme for an irreversible RTP model by defining a backward dynamics with a well-defined path probability and deriving the appropriate Metropolis acceptance rule [2411.12368]. Applied to barrier crossing in a two-dimensional double-well potential, the method shows that intermittent run-and-tumble dynamics can yield shorter transition-path times than both passive diffusion and purely active motion, because tumbles help the particle escape long wall-surfing trajectories and favor more direct reactive channels [2411.12368].

Finally, mean-field interacting RTPs with nonconvex attraction-repulsion can have a far richer stationary structure than equilibrium Brownian systems. For one-dimensional RTPs interacting via
\[
W(r)=-\frac{k_0}{2}r^2+\frac{g}{4}r^4,
\]
the large-\(N\) stationary density has compact support and undergoes a transition from connected to disconnected support at a renormalized \(k_c=3/2^{2/3}\), with multivalued self-consistency below a critical tumbling rate \(\gamma_c=0.1787369\ldots\) and even a continuous family of asymmetric stationary states in the disconnected-support phase [2510.07212]. This is a particularly sharp example of non-uniqueness generated by bounded active noise.

Run-and-tumble particles thus occupy a distinctive place in nonequilibrium statistical mechanics. They are simple enough to admit exact propagators, solvable PDMP reductions, and controlled hydrodynamics, yet rich enough to display confinement without potentials, jamming-induced singularities, MIPS, re-entrant criticality, anomalous collective diffusion, and trajectory-level optimization. A plausible summary of the literature is that RTPs are not merely a minimal active-particle toy model: they are a unifying framework in which persistence, bounded propulsion, and stochastic reorientation can be studied across microscopic, hydrodynamic, geometric, and thermodynamic levels [1206.1805; 1504.07391].

Source: https://www.emergentmind.com/topics/run-and-tumble-particles-rtps