advanced28 min studyLesson 343 in Physics

Diffusion and first-passage problems

Derive first-exit probabilities and mean times from backward equations, impose absorbing or reflecting boundaries, and separate hitting times from ordinary diffusion scales.

On this page

What you’ll understand

  • Use the one-dimensional diffusion generator for hitting questions.
  • Apply absorbing and reflecting boundary conditions to mean first-passage times.
  • Explain why a typical diffusion timescale need not be a finite mean hitting time.

A spreading density answers where particles are at a chosen time. A first-passage problem asks when a trajectory first reaches a target, which is a different question. Boundaries remove or reflect paths and can completely alter the answer. We will derive a compact backward equation, solve an interval problem by hand, and explain why estimating a hitting time by setting mean-square displacement equal to distance squared can give the right scale but the wrong statistical quantity.

Forward densities and backward questions

Consider unbiased one-dimensional Itô diffusion dx=2D dWtdx=\sqrt{2D}\,dW_t with constant positive D. Its density obeys ∂tp=D∂x2p\partial_tp=D\partial_x^2p. The same process has backward generator L=D∂x2\mathcal L=D\partial_x^2, acting on functions of the starting position. For a first exit from the interval between zero and L, the right-exit probability h obeys Dh′′=0Dh''=0, with h equal to zero at the left endpoint and one at the right. Integrating twice gives a linear function of the starting position.

This probability result is independent of D because multiplying the diffusion speed changes how fast a path explores the interval, but not the symmetry between left and right at a given starting point. It relies on zero drift and constant diffusion; a directional drift can change the hitting probabilities. Spatially varying diffusion additionally requires specifying its drift and stochastic interpretation before reusing the simple result. An absorbing boundary means the trajectory is stopped on its first encounter. Continuing to count paths after they cross a target would instead calculate ordinary transition probabilities, which include repeated returns and answer a different question.

Dh′′(x)=0,h(0)=0, h(L)=1,h(x)=xLDh''(x)=0,\qquad h(0)=0,\ h(L)=1,\qquad h(x)=\frac{x}{L}
h is the probability of reaching L before 0 for unbiased diffusion started inside the interval. It is dimensionless; D has length-squared per time units and both endpoints are absorbing.↔ Scroll wide equations horizontally.

The first-step argument produces the mean exit time

Let T(x) be the mean remaining time before absorption at either endpoint. After a short interval dt, the process has spent dt and has a new starting position. Averaging the Taylor expansion of the remaining time gives T(x)=dt+T(x)+DT′′(x)dtT(x)=dt+T(x)+DT''(x)dt to leading order. Cancelling the common term yields the backward Poisson equation DT′′=−1DT''=-1. Its negative curvature reflects a longer mean wait when starting farther from either absorbing boundary. Setting T to zero at both endpoints fixes the two integration constants.

A reflecting endpoint has a different condition. No probability is lost there, and for constant-D unbiased diffusion the corresponding backward condition is zero normal derivative of T. With reflection at zero and absorption at L, use T′(0)=0T'(0)=0 and T(L)=0T(L)=0. The resulting mean time is larger than when zero is absorbing, since trajectories visiting the left endpoint are retained instead of stopped. The boundary condition is part of the stochastic model, not a mathematical convenience that may be changed without changing the physical question.

DT′′(x)=−1,TAA(x)=x(L−x)2D,TRA(x)=L2−x22DDT''(x)=-1,\qquad T_{\mathrm{AA}}(x)=\frac{x(L-x)}{2D},\qquad T_{\mathrm{RA}}(x)=\frac{L^2-x^2}{2D}
AA denotes absorption at both 0 and L; RA denotes reflection at 0 and absorption at L. Both formulas are for constant unbiased diffusion and have units of time.↔ Scroll wide equations horizontally.

Survival tails determine whether a mean exists

The survival probability S(t) is the probability that absorption has not yet occurred. Its negative derivative gives the first-passage density, and the mean passage time equals the integral of S when that integral converges. In a finite interval with an absorbing endpoint, the survival tail decays sufficiently rapidly for the mean to be finite. On the infinite half-line, unbiased diffusion reaches a finite target with probability one, yet its mean hitting time diverges. Rare long excursions make probability-one arrival compatible with an infinite average wait.

This is why the scale ℓ2/D\ell^2/D should not automatically be called a mean first-passage time. It estimates when typical diffusive displacement becomes comparable to a distance ell. Actual means depend on domain size, drift, boundary type, and the tail of the stopping-time distribution. For narrow targets in higher dimensions, geometry matters as well. In experiments finite observation windows censor long trajectories, which can make a divergent or very large theoretical mean appear finite unless the unobserved survival tail is considered.

First exit from an absorbing channel

A particle diffuses on a channel from 0 to L=6 μmL=6\,\mathrm{\mu m}, with both ends absorbing and D=2 μm2 s−1D=2\,\mathrm{\mu m^2\,s^{-1}}. Starting at x=2 μmx=2\,\mathrm{\mu m}, find its right-exit probability and mean time to either end.

  1. Both boundaries stop the process, so use the two-absorbing formulas rather than a reflecting-wall result. The right-exit probability is h=x/L=2/6=1/3h=x/L=2/6=1/3.

  2. The mean exit time is x(L−x)/(2D)=2(6−2)/(2×2)=2 sx(L-x)/(2D)=2(6-2)/(2\times2)=2\,\mathrm s. The numerator is in square micrometres and the denominator in square micrometres per second.

  3. Check the endpoints: the time formula vanishes at either boundary and the right-exit probability runs from zero to one. Reflection symmetry gives the same mean time at x = 4 micrometres, but a right-exit probability of two-thirds there.

The result

The particle exits at the right end with probability one-third, and its mean time to reach either absorbing end is 2 seconds.

Explore the model

Interactive lab

Interactive lab · heat

Diffusion damps shorter spatial scales faster

Reset and set BOTH Spatial mode n = 1 and Diffusivity κ = 0.4; Step twice to t = 1 and predict amplitude exp(−0.4). For a fresh comparison, Reset and set n = 2 and κ = 0.4; Step twice and predict exp(−1.6). The second mode decays four times faster in its exponent.

Diffusion damps shorter spatial scales fasterThe exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Current readings: Current amplitude: 1; Spatial L² norm: 1.772; Decay rate · κn²: 3.6 /s. All parameters are adjustable with the labeled controls below.Periodic heat equation mode 3 at time 0: the amplitude is 1.0π/2π3π/22π−1−0.500.51x · periodic intervalTemperature perturbation uExact mode · n = 3
Current temperature modeInitial sineCurrent amplitude
Current amplitude
1
Spatial L² norm
1.772
Decay rate · κn²
3.6 /s
t = 0 s

The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. 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 forThe heat instrument is one analytic periodic sine mode, not a probability density. It has no absorbing target, reflecting boundary, sampled trajectory, or first-passage readout.

Check your understanding

What follows when unbiased diffusion on a half-line reaches a finite target with probability one?

Try it yourself

Now let zero be reflecting and L=4 μmL=4\,\mathrm{\mu m} absorbing, with D=1.5 μm2 s−1D=1.5\,\mathrm{\mu m^2\,s^{-1}}. Starting at x=1 μmx=1\,\mathrm{\mu m}, derive the mean passage time and check both boundary conditions.

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

Give me a hint

Integrate DT double-prime = −1 twice. Reflection fixes the first derivative at zero rather than the value of T there.

Compare with the solution

Integration gives T(x)=−x2/(2D)+C1x+C2T(x)=-x^2/(2D)+C_1x+C_2. The reflecting condition sets C1=0C_1=0, and absorption at L gives C2=L2/(2D)C_2=L^2/(2D). Thus T(1)=(16−1)/(2×1.5)=5 sT(1)=(16-1)/(2\times1.5)=5\,\mathrm s. Its derivative is −x/D-x/D, which vanishes at zero, while T vanishes at x = L. Setting T(0) to zero instead would incorrectly replace reflection by absorption and shorten the predicted stopping time.

Take it with you

  • First-passage questions stop trajectories at their first target encounter.
  • Backward equations require boundary conditions matching the actual stopping rule.
  • A typical diffusion scale and an average hitting time are not interchangeable.

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.