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
The question for today
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 with constant positive D. Its density obeys . The same process has backward generator , acting on functions of the starting position. For a first exit from the interval between zero and L, the right-exit probability h obeys , 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.
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 to leading order. Cancelling the common term yields the backward Poisson equation . 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 and . 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.
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 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.
Work it through
First exit from an absorbing channel
A particle diffuses on a channel from 0 to , with both ends absorbing and . Starting at , find its right-exit probability and mean time to either end.
Both boundaries stop the process, so use the two-absorbing formulas rather than a reflecting-wall result. The right-exit probability is .
The mean exit time is . The numerator is in square micrometres and the denominator in square micrometres per second.
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 particle exits at the right end with probability one-third, and its mean time to reach either absorbing end is 2 seconds.
The living notebook
Explore the model
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.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /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.
Pause & reason
Check your understanding
Your turn
Try it yourself
Now let zero be reflecting and absorbing, with . Starting at , 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 . The reflecting condition sets , and absorption at L gives . Thus . Its derivative is , 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.
Keep the idea
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.