Why Sum Squared Increments?
A Brownian path is continuous, but its small increments do not behave like those of a smooth curve. To quantify this difference, divide a time interval into short pieces, square each Brownian increment, and add the squares. The resulting sum is called a quadratic variation sum. As the pieces become shorter, the sum converges in probability to the length of the time interval.
This result is a basic reason Brownian motion is a useful noise model in stochastic differential equations: accumulated squared increments remain substantial even as the time steps shrink. The quantity that survives is elapsed time, not zero. The focus here is the precise partition result and its calculations, rather than later rules for manipulating stochastic integrals.
The result concerns any sequence of deterministic partitions \(\pi_k\) for which \(|\pi_k|\to 0\). The partitions need not be equally spaced or nested. Standard Brownian motion has independent increments, and each increment over an interval of length \(h\) has distribution \(N(0,h)\). Those are the facts needed for the derivation.
Deriving the Mean and Variance
Write \(\Delta_i=t_i-t_{i-1}\) and \(X_i=W(t_i)-W(t_{i-1})\). The increments \(X_i\) are independent, with \(X_i\sim N(0,\Delta_i)\). For a centered normal random variable with variance \(v\), its second moment is \(v\) and its fourth moment is \(3v^2\). Thus \(E[X_i^2]=\Delta_i\), and $$ \operatorname{Var}(X_i^2)=E[X_i^4]-(E[X_i^2])^2=3\Delta_i^2-\Delta_i^2=2\Delta_i^2. $$ The independence of the increments also gives independence of their squares. Therefore the sum has expectation equal to the sum of the interval lengths, while its variance is the sum of the individual variances.
The expectation is exactly the elapsed time for every partition, not merely in the limit. The variability of the sum is what shrinks with the mesh. Since each \(\Delta_i\leq|\pi|\), $$ \sum_{i=1}^{m}\Delta_i^2\leq |\pi|\sum_{i=1}^{m}\Delta_i=|\pi|(b-a). $$ This bound gives a direct proof of the theorem.
Because \(E[Q_\pi]=b-a\), the mean-square error is \(E[(Q_\pi-(b-a))^2]=\operatorname{Var}(Q_\pi)\).
The variance formula and the mesh bound imply \(\operatorname{Var}(Q_\pi)=2\sum_i\Delta_i^2\leq2|\pi|(b-a)\).
Along \(\pi_k\), the bound tends to zero because \(|\pi_k|\to0\). Thus \(E[(Q_{\pi_k}-(b-a))^2]\to0\), which proves convergence in mean square.
For completeness, mean-square convergence implies convergence in probability here by Markov's inequality applied to the nonnegative squared error. For any \(\varepsilon>0\), $$ P\bigl(|Q_{\pi_k}-(b-a)|>\varepsilon\bigr) \leq \frac{E[(Q_{\pi_k}-(b-a))^2]}{\varepsilon^2} \leq \frac{2|\pi_k|(b-a)}{\varepsilon^2}\longrightarrow 0. $$ This is precisely convergence in probability. The variance bound also provides a usable finite-partition error bound. \(\square\)
Equal-Spaced Partitions
For an equal-spaced partition of \([0,T]\) into \(n\) intervals, each interval has length \(T/n\). The increment on interval \(i\) can be written \(X_i=\sqrt{T/n}\,Z_i\), where the \(Z_i\) are independent \(N(0,1)\) random variables. Consequently, $$ Q_n=\sum_{i=1}^{n}X_i^2=\frac{T}{n}\sum_{i=1}^{n}Z_i^2. $$ The sum of squares of \(n\) independent standard normal variables has a chi-square distribution with \(n\) degrees of freedom. Hence \(Q_n\) has the distribution \((T/n)\chi_n^2\), with mean \(T\) and variance \(2T^2/n\).
Worked Example: Four Equal Intervals on a Two-Unit Horizon
Take \(T=2\) and \(n=4\). Each interval has length \(2/4=0.5\), so $$ Q_4=0.5\sum_{i=1}^{4}Z_i^2\sim 0.5\chi_4^2. $$
The mean is \(0.5\times4=2\), matching the elapsed time. A chi-square variable with four degrees of freedom has variance \(2(4)=8\), so \(\operatorname{Var}(Q_4)=0.5^2\times8=2\). The general formula gives the same answer: \(2T^2/n=2(2^2)/4=2\).
For a tolerance of 1 time unit, Chebyshev's inequality gives $$ P(|Q_4-2|\geq1)\leq\frac{\operatorname{Var}(Q_4)}{1^2}=2. $$ Since a probability cannot exceed 1, this bound provides no useful assurance at this coarse resolution. It does not mean the probability equals 2; it shows that this particular inequality is too weak to be informative here.
Worked Example: A Finer Equal-Spaced Approximation
Now take \(T=4\), \(n=100\), and ask how likely the sum is to differ from 4 by at least 1. The interval length is \(4/100=0.04\), and $$ E[Q_{100}]=4,\qquad \operatorname{Var}(Q_{100})=\frac{2(4^2)}{100}=0.32. $$
Chebyshev's inequality yields $$ P(|Q_{100}-4|\geq1)\leq\frac{0.32}{1^2}=0.32. $$ The exact distribution is \(Q_{100}=0.04\chi_{100}^2\). The bound is valid but conservative; it guarantees only that the probability is no larger than 0.32. The point of refining the partition is visible in the variance: compared with \(n=25\), for which the variance would be \(2(16)/25=1.28\), using 100 intervals reduces it by a factor of four.
Worked Example: An Unequal Partition and a General Error Bound
Consider a partition of \([0,1]\) whose successive interval lengths are \(0.1,0.1,0.2,0.1,0.2,0.1,0.2\). They sum to \(1\), and the mesh is \(0.2\). The expected squared-increment sum is therefore \(1\). Its variance is $$ 2\sum_i\Delta_i^2 =2(0.01+0.01+0.04+0.01+0.04+0.01+0.04) =2(0.16)=0.32. $$
For a deviation of at least \(1\), Chebyshev's inequality gives \(P(|Q_\pi-1|\geq1)\leq0.32\). The mesh-based bound from the proof gives the same upper bound here: \(2|\pi|(b-a)/1^2=2(0.2)(1)=0.4\), which is valid but looser. In general the mesh bound need not equal the exact variance; it is a simple bound that works for every deterministic partition.
What the Result Does—and Does Not—Say
The convergence statement is probabilistic: for every fixed positive tolerance, the probability of an error larger than that tolerance tends to zero as the mesh shrinks. It does not say that every finite partition gives exactly the elapsed time. Even though \(E[Q_\pi]=b-a\) at every resolution, individual sums vary; their variance quantifies that variation.
It is also important to distinguish the squared-increment sum from the square of the total increment. With two increments \(x_1\) and \(x_2\), the endpoint square is $$ (x_1+x_2)^2=x_1^2+x_2^2+2x_1x_2. $$ It includes both squared contributions as well as a cross-product, but combines them rather than retaining the separate squared increments. For example, if \(x_1=1\) and \(x_2=2\), the endpoint square is \(9\), while the squared-increment sum is \(1^2+2^2=5\). Thus the endpoint square is not a substitute for the partition sum. For Brownian increments, the endpoint square and the squared-increment sum can even have the same expectation, without being the same quantity.
The result assumes deterministic partitions. The proof uses the fact that the interval lengths are fixed and the corresponding Brownian increments are independent normal variables with specified variances. A partition chosen using observed path values is not covered by this argument. The result also does not, by itself, make a claim about convergence along every path. Those are different questions from the convergence-in-probability statement proved here.
Diagnostics and Practical Interpretation
For an equal-spaced calculation, the variance \(2T^2/n\) is a quick diagnostic: if a reported uncertainty does not decrease as \(n\) grows, check the interval length, the squared-increment calculation, and whether the increments have actually been treated as Brownian increments. For an irregular partition, calculate \(2\sum_i\Delta_i^2\); using \(2|\pi|(b-a)\) is a safe upper bound, not generally the exact variance.
In a pharmacokinetic model, this result gives mathematical context for Brownian-driven process variation between observation times. It does not imply that observed concentration differences themselves are Brownian increments: the state dynamics, sampling schedule, and measurement error all matter. Nor is quadratic variation an ordinary slope or endpoint displacement. For a smooth deterministic trajectory with bounded derivative, increments are typically proportional to interval length, so the sum of their squares tends to zero; Brownian increments have variance proportional to interval length, producing a nonzero limiting sum instead.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- For an equal-spaced partition of \([0,3]\) into 30 intervals, calculate the mean and variance of the squared-increment sum.
- A deterministic partition of \([0,2]\) has mesh \(0.1\). Give the mesh-based upper bound on the variance.
- Why is \((W(2)-W(0))^2\) not the squared-increment sum for a partition that includes an intermediate time?
- State the condition on a sequence of deterministic partitions that is sufficient for the theorem.
Solutions. (1) The mean is \(T=3\). The variance is \(2T^2/n=2(9)/30=0.6\). (2) The bound is \(2|\pi|(b-a)=2(0.1)(2)=0.4\). (3) The endpoint square combines squared increments with cross-products; for two increments it is \(x_1^2+x_2^2+2x_1x_2\), whereas the partition sum retains \(x_1^2\) and \(x_2^2\) separately and omits that cross-product. (4) The partitions must be deterministic and their meshes must tend to zero; they need not be equally spaced or nested.