From a Symmetric Walk to a Brownian Path
A simple symmetric random walk moves one unit to the right or left at each step, with equal probabilities. At a fixed number of steps, its typical displacement is of the order of the square root of that number. Donsker’s theorem turns this scaling observation into a statement about entire paths: after rescaling time and space together, the walk approaches Brownian motion on any fixed time interval.
Work on a probability space \((\Omega,\mathcal{F},\mathbb{P})\). Let \((\xi_k)_{k\ge1}\) be independent random variables with \(\mathbb{P}(\xi_k=1)=\mathbb{P}(\xi_k=-1)=\tfrac12\), and define \(S_0=0\) and \(S_m=\sum_{k=1}^m \xi_k\) for \(m\ge1\). The walk \((S_m)_{m\ge0}\) has state space \(\mathbb{Z}\) and discrete time \(m\). For each positive integer \(n\), define a continuous-time, continuous-path process on a fixed interval \([0,T]\) by linearly interpolating the walk between its step times and scaling space by \(1/\sqrt n\):
Here \(\lfloor nt\rfloor\) is the integer part of \(nt\). At a grid time \(t=m/n\), this gives \(W_n(t)=S_m/\sqrt n\). Between grid times, it joins adjacent scaled walk values by a straight line. The time scaling means that \(n\) walk steps correspond to one unit of rescaled time.
The Donsker Result and Its Hypotheses
To apply the theorem, check its hypotheses for the particular increments, not just the shape of the plotted paths. The simple symmetric increments are identically distributed and independent by construction. Their mean and second moment are
Thus the increments are centered and have variance \(1\). They also have finite second moment, as required. Consequently, for every fixed finite \(T\), the scaled, linearly interpolated simple symmetric walk converges in distribution in \(C[0,T]\) to standard Brownian motion. The conclusion is about the random path, not merely the value at one time.
The unit variance is what fixes the Brownian scale. If increments instead had variance \(\sigma^2>0\), the unadjusted space scaling \(1/\sqrt n\) would give limiting covariance \(\sigma^2\min(s,t)\), rather than the standard Brownian covariance \(\min(s,t)\). That broader setting belongs to the later treatment of finite-variance increments; here the increments have variance exactly \(1\).
Checks on the Scaling
Two elementary calculations help interpret the theorem. First, interpolation changes the stepwise walk by at most one step’s scaled size. Second, the covariance of the interpolated process approaches the Brownian covariance. These facts do not replace Donsker’s theorem: convergence of covariances alone would not establish convergence of path distributions. They do provide direct checks that the interpolation and scaling have the expected size and dependence.
Proof. Write \(r=nt-\lfloor nt\rfloor\), so \(0\le r<1\). The difference between the interpolated and stepwise values is \(r\xi_{\lfloor nt\rfloor+1}/\sqrt n\). Since \(|\xi_k|=1\), its absolute value is \(r/\sqrt n\le1/\sqrt n\).
For the covariance, it is enough to take \(s\le t\), since covariance is symmetric. Write \(ns=a+\alpha\) and \(nt=b+\beta\), where \(a,b\) are nonnegative integers and \(0\le\alpha,\beta<1\). The process value at time \(s\) has coefficient \(1/\sqrt n\) on each of \(\xi_1,\ldots,\xi_a\) and coefficient \(\alpha/\sqrt n\) on \(\xi_{a+1}\). The corresponding coefficients at \(t\) are \(1/\sqrt n\) through \(\xi_b\) and \(\beta/\sqrt n\) on \(\xi_{b+1}\). Distinct increments have covariance zero, and each increment has variance \(1\). If \(a<b\), the shared coefficients give covariance \((a+\alpha)/n=s\). If \(a=b\), the covariance is \((a+\alpha\beta)/n\), so its difference from \(s=(a+\alpha)/n\) is \(\alpha(1-\beta)/n\). In this case \(\alpha\le\beta\), and \(0\le\alpha(1-\beta)\le1\), so the difference is at most \(1/n\). In both cases the covariance differs from \(\min(s,t)=s\) by at most \(1/n\), which tends to zero. \(\square\)
The bound on interpolation is uniform over the whole time interval: \(\sup_{0\le t\le T}|W_n(t)-U_n(t)|\le1/\sqrt n\). It explains why joining the grid values does not create a large extra fluctuation. The covariance calculation checks the variance scale and the dependence between times; the theorem supplies the stronger path-distribution conclusion.
Worked Example: Interpolating a Four-Step Path
Take \(n=4\) and consider the particular increment sequence \((\xi_1,\xi_2,\xi_3,\xi_4)=(1,1,-1,1)\). Compute the interpolated scaled path at its grid times and at \(t=\tfrac38\), with \(T=1\).
The partial sums are \(S_0=0\), \(S_1=1\), \(S_2=2\), \(S_3=1\), and \(S_4=2\). Since \(\sqrt n=2\), the values of \(W_4\) at \(t=0,\tfrac14,\tfrac12,\tfrac34,1\) are, respectively, \(0,\tfrac12,1,\tfrac12,1\). At \(t=\tfrac38\), \(nt=\tfrac32\), so \(\lfloor nt\rfloor=1\) and \(nt-\lfloor nt\rfloor=\tfrac12\). Therefore
Check: \(t=\tfrac38\) is halfway between grid times \(\tfrac14\) and \(\tfrac12\), where the scaled path has values \(\tfrac12\) and \(1\). Their midpoint is \((\tfrac12+1)/2=\tfrac34\), agreeing with the interpolation formula.
Worked Example: The Endpoint at Time One
For \(n=4\), calculate the exact distribution, mean, and variance of \(W_4(1)\). Then check how the endpoint calculation extends to general \(n\).
At time \(1\), \(W_4(1)=S_4/2\). There are \(2^4=16\) equally likely sign sequences. The possible sums \(-4,-2,0,2,4\) have counts \(1,4,6,4,1\), respectively, so the scaled endpoint takes values \(-2,-1,0,1,2\) with probabilities \(\tfrac1{16},\tfrac4{16},\tfrac6{16},\tfrac4{16},\tfrac1{16}\). Symmetry gives mean zero. Its variance is
For general \(n\), \(W_n(1)=S_n/\sqrt n\). Independence and centering give \(\mathbb{E}S_n=0\) and \(\operatorname{Var}(S_n)=\sum_{k=1}^n\operatorname{Var}(\xi_k)=n\). Hence \(\mathbb{E}W_n(1)=0\) and \(\operatorname{Var}(W_n(1))=n/n=1\) for every \(n\). The endpoint’s characteristic function also shows the normal scaling: independence gives \(\mathbb{E}e^{iuW_n(1)}=(\cos(u/\sqrt n))^n\), and for fixed real \(u\), \(n\log\cos(u/\sqrt n)\to-u^2/2\). Thus this characteristic function tends to \(e^{-u^2/2}\), the characteristic function of a standard normal random variable.
Check: At \(n=4\), the displayed probabilities sum to \((1+4+6+4+1)/16=1\), and the computed variance is \(1\), matching the general formula \(n/n=1\).
Worked Example: Comparing Two Times
For \(n=4\), calculate the covariance of \(W_4(\tfrac12)\) and \(W_4(1)\), and compare it with the covariance of Brownian motion at those times.
Both times are grid times, so \(W_4(\tfrac12)=S_2/2\) and \(W_4(1)=S_4/2\). Since \(S_4=S_2+\xi_3+\xi_4\), and \(\xi_3,\xi_4\) are independent of \(S_2\) with mean zero, \(\operatorname{Cov}(S_2,S_4)=\operatorname{Var}(S_2)\). The first two increments are independent, centered, and each has variance \(1\), so \(\operatorname{Var}(S_2)=1+1=2\). Therefore
Standard Brownian motion has covariance \(\operatorname{Cov}(B_s,B_t)=\min(s,t)\), so \(\operatorname{Cov}(B_{\tfrac12},B_1)=\tfrac12\). Check: The covariance calculation for the walk is exactly \(\tfrac12\), the same value as the Brownian covariance at these two times.
Why Path Convergence Is Stronger Than an Endpoint Limit
The endpoint calculation concerns only \(W_n(1)\). Even convergence of the joint distributions at every finite collection of times would not, by itself, rule out irregular behavior between those times. Donsker’s theorem establishes convergence as random elements of \(C[0,T]\); its proof requires both control of finite-dimensional distributions and tightness of the path laws. The tightness argument is a substantive part of the earlier proof structure, not a consequence of the variance calculation alone.
Path convergence is useful when the quantity of interest depends on a whole trajectory, such as the maximum over an interval or the time spent above a level. For example, the map \(f\mapsto\sup_{0\le t\le T}f(t)\) is continuous in the uniform norm: for any \(f,g\in C[0,T]\),
Thus, by the continuous mapping theorem, Donsker convergence also gives convergence in distribution of the path maxima to \(\sup_{0\le t\le T}B_t\). This conclusion concerns the maximum of the interpolated path; interpolation makes that maximum a well-defined continuous-path functional.
Common scaling error. Dividing the walk by \(n\), rather than \(\sqrt n\), makes the fluctuations vanish. For instance, at time \(1\) with \(n=16\), the wrongly scaled endpoint \(S_{16}/16\) has variance \(16/16^2=1/16\), whereas the correctly scaled endpoint \(S_{16}/\sqrt{16}=S_{16}/4\) has variance \(16/16=1\). In general, the incorrect variance is \(n/n^2=1/n\to0\); it cannot yield a nondegenerate standard Brownian endpoint. The square-root scaling is dictated by the accumulated variance of \(n\) independent unit-variance increments.
Check Your Understanding
Use the definitions and calculations above to answer these questions.
- For the simple symmetric increment, calculate its mean, variance, and fourth moment. Which of these values is needed for the stated form of Donsker’s theorem?
- For \(n=9\), what is the space-scaling factor, and at which rescaled time does the walk’s ninth step occur?
- Explain why finite-dimensional convergence alone does not establish convergence of the whole interpolated path.
- For \(n=4\), calculate \(\operatorname{Var}(W_4(\tfrac12))\) directly from \(S_2\).
- What is the maximum possible difference between the interpolated and stepwise scaled paths, uniformly over \(0\le t\le T\), for a fixed \(n\)?