Tutorials › Stochastic Processes › Donsker Convergence for the Simple Symmetric Walk

Brownian Motion: Definition and Construction · Tutorial 816 of 994

Donsker Convergence for the Simple Symmetric Walk

See why the linearly interpolated simple symmetric walk, scaled diffusively, converges to standard Brownian motion on every fixed time interval.

Advanced 10 min read

What You'll Learn

  • Define the polygonally interpolated, diffusively scaled simple symmetric walk.
  • Verify the increment assumptions required by Donsker’s theorem.
  • Interpret convergence in distribution in the space of continuous paths.
  • Check endpoint and covariance calculations for the scaled walk.
  • Identify why linear interpolation and square-root scaling matter.

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\):

$$ W_n(t)=\frac{S_{\lfloor nt\rfloor}+(nt-\lfloor nt\rfloor)\xi_{\lfloor nt\rfloor+1}}{\sqrt n},\qquad 0\le t\le T. $$

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.

Definition: Convergence in distribution of \(W_n\) to a continuous-path process \(B\) on \([0,T]\) means weak convergence of their probability laws on \(C[0,T]\), the space of continuous real-valued functions with the uniform norm. Equivalently, \(\mathbb{E}[F(W_n)]\to\mathbb{E}[F(B)]\) for every bounded continuous function \(F:C[0,T]\to\mathbb{R}\).

The Donsker Result and Its Hypotheses

Theorem (Donsker’s invariance principle, unit-variance form): Let \((\xi_k)_{k\ge1}\) be independent and identically distributed real random variables with \(\mathbb{E}\xi_1=0\) and \(\mathbb{E}[\xi_1^2]=1\). Form the linearly interpolated processes \(W_n\) using their partial sums and the scaling above. For every fixed \(T<\infty\), \(W_n\) converges in distribution in \(C[0,T]\) to standard Brownian motion \((B_t)_{0\le t\le T}\). This is a standard functional limit theorem; its proof structure and tightness argument are treated in the earlier material on Donsker’s invariance principle.

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

$$ \mathbb{E}\xi_k=\tfrac12(1)+\tfrac12(-1)=0,\qquad \mathbb{E}[\xi_k^2]=\tfrac12(1^2)+\tfrac12((-1)^2)=1. $$

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.

Proposition: Let \(U_n(t)=S_{\lfloor nt\rfloor}/\sqrt n\) be the stepwise version of the scaled walk, considered on \([0,T]\). At every \(t\in[0,T]\), \(|W_n(t)-U_n(t)|\le 1/\sqrt n\). Also, for \(s,t\in[0,T]\), \(\operatorname{Cov}(W_n(s),W_n(t))\to\min(s,t)\) as \(n\to\infty\).

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

$$ W_4(\tfrac38)=\frac{S_1+\tfrac12\xi_2}{2} =\frac{1+\tfrac12(1)}{2} =\tfrac34. $$

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

$$ \operatorname{Var}(W_4(1)) =\tfrac1{16}(4)+\tfrac4{16}(1)+\tfrac6{16}(0) +\tfrac4{16}(1)+\tfrac1{16}(4) =1. $$

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

$$ \operatorname{Cov}(W_4(\tfrac12),W_4(1)) =\frac{\operatorname{Cov}(S_2,S_4)}{4} =\frac{2}{4} =\tfrac12. $$

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]\),

$$ \left|\sup_{0\le t\le T}f(t)-\sup_{0\le t\le T}g(t)\right| \le \sup_{0\le t\le T}|f(t)-g(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.

  1. 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?
  2. For \(n=9\), what is the space-scaling factor, and at which rescaled time does the walk’s ninth step occur?
  3. Explain why finite-dimensional convergence alone does not establish convergence of the whole interpolated path.
  4. For \(n=4\), calculate \(\operatorname{Var}(W_4(\tfrac12))\) directly from \(S_2\).
  5. What is the maximum possible difference between the interpolated and stepwise scaled paths, uniformly over \(0\le t\le T\), for a fixed \(n\)?