advanced27 min studyLesson 341 in Physics

Langevin forces and white-noise models

Connect coarse-grained random forces to Wiener increments, track noise units, and distinguish the Itô and Stratonovich interpretations of multiplicative noise.

On this page

What you’ll understand

  • Normalize white noise through the variance of its integrated increments.
  • Identify when overdamped motion and an equilibrium Einstein relation apply.
  • Convert scalar multiplicative-noise drift between Itô and Stratonovich conventions.

A Langevin equation replaces unresolved microscopic collisions by drag and fluctuating forcing. It is a model of an observed slow variable, not a claim that molecular forces literally have zero correlation time. Its usefulness depends on a separation of timescales and on correct normalization. We will state the noise convention explicitly, connect it to measurable displacement statistics, and identify the extra decision required when noise strength depends on the particle's position.

White noise is defined by integrated kicks

For a one-dimensional inertial particle, write m dv=(−ζv+F)dt+2ζkBT dWtm\,dv=(-\zeta v+F)dt+\sqrt{2\zeta k_BT}\,dW_t. Here zeta is a friction coefficient in kilograms per second, F is a resolved force, and WtW_t is a standard Wiener process with independent increments of zero mean and variance dt. The random force itself is a generalized process, formally proportional to the time derivative of W. Its pointwise value and variance are not ordinary finite numbers. Measurements always average its effect over a nonzero interval.

The ideal delta correlation is appropriate when the bath correlation time is short compared with resolved times. At very short times molecular memory and particle inertia can matter; a colored-noise model or generalized friction kernel may be needed. The thermal amplitude above is tied to equilibrium bath temperature and friction. Choosing it independently of friction describes a different model and generally loses the thermal equilibrium distribution. Zero mean random force does not imply zero accumulated variance: independent kicks cancel in average displacement while their squared effects add.

⟨η(t)η(s)⟩=δ(t−s),ΔW=∫tt+Δtη(s) ds,⟨ΔW⟩=0,⟨(ΔW)2⟩=Δt\langle\eta(t)\eta(s)\rangle=\delta(t-s),\qquad \Delta W=\int_t^{t+\Delta t}\eta(s)\,ds,\qquad \langle\Delta W\rangle=0,\quad \langle(\Delta W)^2\rangle=\Delta t
Eta is unit white noise with units s⁻¹ᐟ²; a Wiener increment has units s¹ᐟ². The delta correlation is understood under integration, not as a finite force variance at one instant.↔ Scroll wide equations horizontally.

The overdamped limit retains fluctuations

At times long compared with m/ζm/\zeta, a suitable overdamped description is dx=μFdt+2D dWtdx=\mu Fdt+\sqrt{2D}\,dW_t, with mobility μ=1/ζ\mu=1/\zeta. Mobility has units of velocity divided by force and D has units of length squared per time. For an isothermal equilibrium thermal bath and constant mobility, D=μkBTD=\mu k_BT. This Einstein relation is an assumption justified by thermal consistency, not a definition valid for every fluctuating environment. An actively driven bath can violate it even when its position trajectories appear diffusive.

For constant force and diffusion, the mean displacement grows linearly with time, while the centered variance is twice D times time in one dimension. The uncentered mean-square displacement also includes the square of the mean drift. Confusing these quantities leads to an apparent diffusion coefficient that increases with observation time. In several independent Cartesian directions each contributes its own variance; isotropic three-dimensional diffusion therefore has mean-square displacement six D times time after drift is subtracted. The factor two is already present in each coordinate's noise normalization.

dx=μF dt+2D dWt,⟨Δx⟩=μFt,Var⁡(Δx)=2Dt,D=μkBT for a thermal equilibrium bathdx=\mu F\,dt+\sqrt{2D}\,dW_t,\qquad \langle\Delta x\rangle=\mu Ft,\qquad \operatorname{Var}(\Delta x)=2Dt,\qquad D=\mu k_BT\ \text{for a thermal equilibrium bath}
The displacement formulas assume constant F, mu and D, a fixed initial position, and resolved times beyond inertial relaxation. The noise amplitude sqrt(2D) has units length divided by sqrt(time).↔ Scroll wide equations horizontally.

Multiplicative noise requires a calculus convention

Unless stated otherwise this module writes stochastic differential equations in the Itô convention. For dx=aI(x)dt+b(x)dWtdx=a_I(x)dt+b(x)dW_t, the coefficient of each increment is evaluated at its left endpoint. A Stratonovich equation uses a midpoint prescription and is marked with a circle on its stochastic differential. For constant b both conventions give the same paths and drift. For varying b they describe the same process only after the drift is changed by aI=aS+bb′/2a_I=a_S+bb'/2, with the derivative taken with respect to x.

This drift conversion is not an optional numerical correction. Keeping the same displayed drift while switching interpretations generally changes the stochastic process and its stationary density. The ordinary chain rule holds in Stratonovich form; Itô's formula instead includes the quadratic-variation term. The correct interpretation of a reduced physical model follows from how fast degrees of freedom were eliminated, including the noise and friction timescales. Neither notation is universally more physical. A complete model declares the interpretation and any required noise-induced drift before computing predictions.

Separate drift from spreading

An overdamped particle starts at x = 0 with constant drift μF=2 μm s−1\mu F=2\,\mathrm{\mu m\,s^{-1}} and D=0.5 μm2 s−1D=0.5\,\mathrm{\mu m^2\,s^{-1}}. At t = 4 s, find its mean displacement, centered variance, and uncentered mean-square displacement.

  1. The deterministic drift contributes ⟨x⟩=(2)(4)=8 μm\langle x\rangle=(2)(4)=8\,\mathrm{\mu m}. This is the ensemble mean, not the position of every trajectory.

  2. The Wiener increments contribute Var⁡(x)=2Dt=2(0.5)(4)=4 μm2\operatorname{Var}(x)=2Dt=2(0.5)(4)=4\,\mathrm{\mu m^2}, corresponding to a standard deviation of 2 μm2\,\mathrm{\mu m} around the drifting mean.

  3. Use ⟨x2⟩=Var⁡(x)+⟨x⟩2\langle x^2\rangle=\operatorname{Var}(x)+\langle x\rangle^2. The result is 4+64=68 μm24+64=68\,\mathrm{\mu m^2}. Dividing this uncentered quantity by twice the elapsed time would incorrectly mix transport with diffusion.

The result

The mean is 8 micrometres, the centered variance is 4 square micrometres, and the mean-square displacement from the origin is 68 square micrometres.

Explore the model

Interactive lab

Interactive lab · sampling

Independent averages narrow without changing their mean

For a fresh run, Reset and set BOTH Success probability p = 0.5 and Independent samples n = 20. Predict variance 0.0125 for the sample mean. Then Reset and set p = 0.5 and n = 80; predict variance 0.003125 and half the standard deviation. Inspect exact readouts or Step.

Independent averages narrow without changing their meanExact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Current readings: Mean sample proportion: 0.5; Sample-mean variance: 0.0125; Uniform CDF error bound: ≤ 0.1062. All parameters are adjustable with the labeled controls below.Exact distribution of the mean of 20 independent Bernoulli variables with success probability 0.5; playback highlights outcome 0.00.250.50.75100.0580.1150.1730.23Sample mean · k / nProbability massNormal bin masses are an approximationk = 0 · exact probability = 9.54e-7
Exact binomial massNormal interval massHighlighted outcome / mean
Mean sample proportion
0.5
Sample-mean variance
0.0125
Uniform CDF error bound
≤ 0.1062
t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). 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 an exact binomial sampling distribution, not simulated Wiener noise or resampled trajectories. Playback highlights possible outcomes; it does not generate a bath realization.

Check your understanding

What happens when multiplicative noise is changed from Stratonovich to Itô notation for the same scalar process?

Try it yourself

An isothermal particle has ζ=2 pN s μm−1\zeta=2\,\mathrm{pN\,s\,\mu m^{-1}}, kBT=0.004 pN μmk_BT=0.004\,\mathrm{pN\,\mu m}, and constant force F = 4 pN. Find its drift speed, D, and the centered displacement standard deviation over 0.1 s.

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

Give me a hint

Invert the friction coefficient for mobility, apply the thermal Einstein relation, and use the square root of twice D times time.

Compare with the solution

The mobility is μ=0.5 μm (pN s)−1\mu=0.5\,\mathrm{\mu m\,(pN\,s)^{-1}}, giving drift μF=2 μm s−1\mu F=2\,\mathrm{\mu m\,s^{-1}}. Thermal diffusion is D=μkBT=0.002 μm2 s−1D=\mu k_BT=0.002\,\mathrm{\mu m^2\,s^{-1}}. The centered variance over 0.1 s is 2(0.002)(0.1)=0.0004 μm22(0.002)(0.1)=0.0004\,\mathrm{\mu m^2} and its standard deviation is 0.02 μm0.02\,\mathrm{\mu m}. These values assume the interval exceeds inertial and bath-memory times; no such overdamped prediction is implied at arbitrarily short times.

Take it with you

  • White noise is normalized by integrated increments with variance proportional to time.
  • Centered spreading and drift contribute different terms to mean-square displacement.
  • Multiplicative noise requires a declared interpretation and a consistent drift conversion.

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.