advanced29 min studyLesson 350 in Physics

Active matter and broken equilibrium relations

Derive enhanced diffusion from a persistent telegraph model, distinguish active transport from thermal mobility, and separate coarse-grained predictions from unresolved fuel dissipation.

On this page

What you’ll understand

  • Derive a persistent-velocity contribution to long-time diffusion.
  • Test the thermal Einstein relation without redefining temperature by assumption.
  • Distinguish exact model consequences, experimental observations, and unresolved thermodynamic information.

An active particle converts an energy supply into persistent motion. Its long-time trajectory can look diffusive, yet its spreading need not be tied to thermal drag by the equilibrium Einstein relation. This makes active matter a useful endpoint for the module: stationary distributions, noise models, response, and entropy all require their assumptions to be revisited. We will derive a simple persistent-motion prediction exactly within its model, then distinguish that derivation from claims about real swimmers and their fuel consumption.

A persistent direction is an extra state variable

Consider a one-dimensional telegraph model with direction σ=±1\sigma=\pm1 and speed v0. Direction flips from either sign to the other at rate alpha. This is explicitly a flip rate, not a tumble rate that may choose the same direction again. In the stationary equal-direction ensemble the correlation is ⟨σ(t)σ(0)⟩=e−2α∣t∣\langle\sigma(t)\sigma(0)\rangle=e^{-2\alpha|t|}. The factor two follows because either a forward or backward transition relaxes a direction imbalance. A different convention for random reorientation rates would produce a different parameter in the exponential.

Add passive mobility and thermal diffusion through dx=(v0σ+μF)dt+2Dth dWtdx=(v_0\sigma+\mu F)dt+\sqrt{2D_{\mathrm{th}}}\,dW_t, with independent thermal noise. For the passive bath, Dth=μkBTD_{\mathrm{th}}=\mu k_BT can still hold. The persistent propulsion is an additional nonequilibrium contribution rather than a change in that bath temperature. At zero force the ensemble mean displacement vanishes when the initial directions are equally probable, but each short trajectory can retain a preferred direction. Zero mean transport therefore does not make the propulsion equivalent to thermal white noise.

⟨va(t)va(0)⟩=v02e−2α∣t∣,τp=12α,Deff=Dth+∫0∞⟨va(t)va(0)⟩dt=Dth+v022α\langle v_{\mathrm a}(t)v_{\mathrm a}(0)\rangle=v_0^2e^{-2\alpha|t|},\qquad \tau_{\mathrm p}=\frac{1}{2\alpha},\qquad D_{\mathrm{eff}}=D_{\mathrm{th}}+\int_0^\infty\langle v_{\mathrm a}(t)v_{\mathrm a}(0)\rangle dt=D_{\mathrm{th}}+\frac{v_0^2}{2\alpha}
This is the long-time diffusion of a noninteracting one-dimensional particle with independent Poisson direction flips and independent thermal noise. Speed has units length/time, alpha has inverse-time units, and each D has length²/time units.↔ Scroll wide equations horizontally.

Short ballistic motion becomes enhanced diffusion

Integrating the velocity correlation twice gives the active contribution to displacement variance. At short times it grows as v02t2v_0^2t^2, expressing persistence, while the independent thermal part grows linearly. At times long compared with the persistence time, the active part becomes linear and adds v02/(2α)v_0^2/(2\alpha) to the effective diffusion coefficient. The noncommuting limits are instructive: setting alpha to zero first yields persistent ballistic motion, so the long-time diffusive formula cannot be used at zero flip rate.

Original experiments on chemically propelled colloids observed short-time directed motion and enhanced long-time spreading. That observation supports persistent active descriptions, but does not uniquely select this one-dimensional flip model or its numerical coefficient. Real swimmers rotate in two or three dimensions, interact with walls, and can change speed with fuel concentration. Their orientational correlation and the dimensional diffusion factors must be modeled or measured separately. The formula above is an exact deduction from stated ideal assumptions, not a parameter-free prediction for every experimental active particle.

⟨(Δx)2⟩=2Dtht+v02α[t−1−e−2αt2α](F=0),lim⁡t→∞⟨(Δx)2⟩2t=Deff\langle(\Delta x)^2\rangle=2D_{\mathrm{th}}t+\frac{v_0^2}{\alpha}\left[t-\frac{1-e^{-2\alpha t}}{2\alpha}\right]\quad(F=0),\qquad \lim_{t\to\infty}\frac{\langle(\Delta x)^2\rangle}{2t}=D_{\mathrm{eff}}
The initial direction is sampled from its stationary symmetric distribution and the initial position is fixed. The formula is one dimensional and assumes alpha > 0. Its small-time active term is ballistic rather than diffusive.↔ Scroll wide equations horizontally.

An effective diffusion does not define a full thermodynamics

In this simple model a weak uniform force changes the mean drift by mu times F without altering the propulsion statistics. The ratio Deff/μD_{\mathrm{eff}}/\mu therefore exceeds the bath thermal energy whenever propulsion is nonzero. Calling that ratio an effective temperature is a limited parametrization of one transport comparison. It does not establish Boltzmann statistics in arbitrary traps, a universal fluctuation–dissipation relation, or standard thermal work identities for all active protocols. Such claims require additional derivations and can fail across observables or timescales.

The coarse-grained propulsion variable also omits chemical fuel and the mechanism that turns fuel into force. Position and direction trajectories alone may not determine the total physical entropy production. A full energetic model must identify reservoirs, fuel reactions, time-reversal rules, and hidden dissipation. Interacting active particles can display collective organization and phase separation in suitable models, but those phenomena do not follow from the single-particle formula alone. Separate exact consequences of a specified stochastic model from empirical observations and broader open questions about real active materials.

Enhanced spreading with unchanged passive mobility

A telegraph swimmer has v0=1.2 μm s−1v_0=1.2\,\mathrm{\mu m\,s^{-1}}, α=0.6 s−1\alpha=0.6\,\mathrm{s^{-1}}, Dth=0.20 μm2 s−1D_{\mathrm{th}}=0.20\,\mathrm{\mu m^2\,s^{-1}}, and μ=50 μm (pN s)−1\mu=50\,\mathrm{\mu m\,(pN\,s)^{-1}}. Its bath has kBT=0.004 pN μmk_BT=0.004\,\mathrm{pN\,\mu m}. Find persistence time, effective diffusion, and the ratio Deff/mu.

  1. The direction correlation time is τp=1/(2α)=1/1.2=5/6 s\tau_{\mathrm p}=1/(2\alpha)=1/1.2=5/6\,\mathrm s. The passive Einstein relation is consistent because μkBT=50(0.004)=0.20 μm2 s−1\mu k_BT=50(0.004)=0.20\,\mathrm{\mu m^2\,s^{-1}}.

  2. The active diffusion contribution is v02/(2α)=1.44/1.2=1.20 μm2 s−1v_0^2/(2\alpha)=1.44/1.2=1.20\,\mathrm{\mu m^2\,s^{-1}}. Adding the passive value gives Deff=1.40 μm2 s−1D_{\mathrm{eff}}=1.40\,\mathrm{\mu m^2\,s^{-1}}.

  3. The transport ratio is Deff/μ=1.40/50=0.028 pN μmD_{\mathrm{eff}}/\mu=1.40/50=0.028\,\mathrm{pN\,\mu m}, seven times the actual bath thermal energy. This comparison does not measure a sevenfold increase of the bath temperature or determine chemical heat dissipation.

The result

The persistence time is about 0.833 s and long-time diffusion is 1.40 square micrometres per second. The equilibrium Einstein relation fails for total active spreading despite remaining valid for the passive bath contribution.

Explore the model

Interactive lab

Interactive lab · markov

Faster direction switching shortens memory

Reset and set BOTH Transition p = 0.2 and Transition q = 0.2, then Step once; predict state-1 probability 0.2. For a fresh comparison, Reset and set p = 0.4 and q = 0.4 and Step once; predict 0.4. Both approach equal stationary weights, but the second loses its initial state-0 bias faster.

Faster direction switching shortens memoryExact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. Current readings: P(state 1) at step 0: 0%; Stationary P(state 1): 60%; Chain behavior: Monotone convergence. All parameters are adjustable with the labeled controls below.Two-state Markov distribution after 0 steps: state zero probability 1, state one probability 0.Initial distribution (1, 0) · step 0p = 0.3q = 0.201100%0%1−p1−qP(0)P(1)Monotone convergence
P(state 0)P(state 1)Stationary distribution
P(state 1) at step 0
0%
Stationary P(state 1)
60%
Chain behavior
Monotone convergence
t = 0 s

Exact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThis is a discrete two-state probability model, not a continuous telegraph trajectory or active swimmer. Its transition probabilities are not physical flip rates, and it has no position, fuel, diffusion, or entropy-production readout.

Check your understanding

What does a measured Deff/mu larger than kBT imply by itself?

Try it yourself

A one-dimensional telegraph particle has v0=3 μm s−1v_0=3\,\mathrm{\mu m\,s^{-1}}, Dth=0.1 μm2 s−1D_{\mathrm{th}}=0.1\,\mathrm{\mu m^2\,s^{-1}}, and initial flip rate α=1.5 s−1\alpha=1.5\,\mathrm{s^{-1}}. Find Deff and persistence time before and after doubling alpha, with speed and thermal diffusion unchanged.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Only the active term and persistence time scale inversely with flip rate. The passive diffusion contribution remains fixed.

Compare with the solution

Initially the active contribution is 9/(2×1.5)=3 μm2 s−19/(2\times1.5)=3\,\mathrm{\mu m^2\,s^{-1}}, so Deff=3.1 μm2 s−1D_{\mathrm{eff}}=3.1\,\mathrm{\mu m^2\,s^{-1}} and τp=1/3 s\tau_{\mathrm p}=1/3\,\mathrm s. At α=3 s−1\alpha=3\,\mathrm{s^{-1}}, the active contribution is 1.5, total diffusion is 1.6 μm2 s−11.6\,\mathrm{\mu m^2\,s^{-1}}, and persistence time is 1/6 s1/6\,\mathrm s. Total diffusion does not exactly halve because its passive 0.1 contribution is unchanged. These comparisons concern separately specified steady models, not a continuously switched protocol.

Take it with you

  • Persistent active velocity adds a correlation-integral contribution to long-time diffusion.
  • A diffusion–mobility ratio alone does not establish a universal thermodynamic temperature.
  • Hidden propulsion mechanisms require extra information before total entropy production can be assigned.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.