Tutorials › Stochastic Processes › Gaussian Increment Moments and Path Regularity

Path Properties of Brownian Motion · Tutorial 823 of 994

Gaussian Increment Moments and Path Regularity

Learn how Gaussian moments determine Brownian increment bounds and how their exponents guide moment-based regularity arguments.

Advanced 9 min read

What You'll Learn

  • Derive the even moments of a standard normal variable by integration by parts.
  • Compute even moments of Brownian increments, including the zero-length case.
  • Use moment scaling to bound the probability of a large increment.
  • Determine when a moment exponent is strong enough for regularity arguments.
  • Recognize why a single-increment estimate is not itself a path-continuity proof.

Why Increment Moments Matter

Brownian path regularity is studied by comparing values at nearby times. The distribution of a single increment gives more than its variance: it gives moments of every order. Higher moments penalize large increments more strongly, and their dependence on the time gap provides estimates that can be used in continuity and Hölder-regularity arguments. The crucial calculation is exact: an increment over a time interval of length \(h\) has its \(2m\)-th absolute moment proportional to \(h^m\).

Throughout, let \((\Omega,\mathcal{F},\mathbb{P})\) be a probability space and let \((B_t)_{t\ge0}\) be a real-valued standard Brownian motion with state space \(S=\mathbb{R}\) and continuous-time index \(t\ge0\). We use the defining increment law: for \(0\le s\le t\), \(B_t-B_s\) is Gaussian with mean \(0\) and variance \(t-s\); when \(t=s\), the increment is identically zero. The estimates below concern the distribution of increments and do not by themselves prove a path-regularity theorem.

Gaussian Moments and Brownian Scaling

Theorem (Even Gaussian moments): If \(Z\) has the standard normal distribution, then for every integer \(m\ge1\), \[ \mathbb{E}[Z^{2m}]=(2m-1)(2m-3)\cdots 3\cdot1=(2m-1)!!. \] By symmetry, every odd signed moment of \(Z\) is zero. For \(h\ge0\), a centered Gaussian variable \(Y\) with variance \(h\) therefore satisfies \(\mathbb{E}|Y|^{2m}=(2m-1)!!h^m\), including \(h=0\).

Proof. For \(Z\) standard normal, its density is \(\phi(x)=(2\pi)^{-1/2}e^{-x^2/2}\). Set \(I_m=\int_{-\infty}^{\infty}x^{2m}\phi(x)\,dx\), so \(I_m=\mathbb{E}[Z^{2m}]\). Since \(\phi'(x)=-x\phi(x)\), integration by parts gives, for \(m\ge1\), \[ I_m=-\int_{-\infty}^{\infty}x^{2m-1}\phi'(x)\,dx =-(x^{2m-1}\phi(x))\big|_{-\infty}^{\infty} +(2m-1)\int_{-\infty}^{\infty}x^{2m-2}\phi(x)\,dx. \] The boundary term is zero because a polynomial times \(e^{-x^2/2}\) tends to zero at both ends of the real line. Thus \(I_m=(2m-1)I_{m-1}\). Since \(I_0=1\), repeated application gives \(I_m=(2m-1)(2m-3)\cdots1\). The density is symmetric, so the integral of any odd power over the real line is zero. If \(h>0\), write \(Y=\sqrt{h}Z\) in distribution; then \(\mathbb{E}|Y|^{2m}=h^m\mathbb{E}|Z|^{2m}\). If \(h=0\), \(Y=0\) almost surely, and both sides are zero. \(\square\)

Applying this result to a Brownian increment gives the central formula. The absolute value can be omitted for an even power, but retaining it makes the nonnegativity of the moment explicit.

$$ \mathbb{E}|B_t-B_s|^{2m}=(2m-1)!!(t-s)^m,\qquad 0\le s\le t,\quad m=1,2,\ldots . $$

The moment order and time exponent are linked: a moment of order \(2m\) scales as the \(m\)-th power of the interval length. In particular, multiplying the time gap by \(c\ge0\) multiplies this moment by \(c^m\). This is the moment form of Brownian scaling; it follows directly from the Gaussian increment law.

Worked Example: The Second, Fourth, and Sixth Moments

Take \(s=\tfrac12\) and \(t=2\). The interval length is \(h=t-s=2-\tfrac12=\tfrac32\). For \(m=1\), the second absolute moment is \[ \mathbb{E}|B_2-B_{1/2}|^2=1\cdot\tfrac32=\tfrac32. \] For \(m=2\), the fourth moment is \[ \mathbb{E}|B_2-B_{1/2}|^4=3\left(\tfrac32\right)^2 =3\cdot\tfrac94=\tfrac{27}4. \] For \(m=3\), the sixth moment is \[ \mathbb{E}|B_2-B_{1/2}|^6=15\left(\tfrac32\right)^3 =15\cdot\tfrac{27}8=\tfrac{405}8. \]

Check: Dividing each moment by its corresponding power of \(h=\tfrac32\) gives \(1\), \(3\), and \(15\), the standard-normal second, fourth, and sixth moments. Thus the scale factors are consistent with the Gaussian law.

Worked Example: Standardizing an Increment

Let \(s=1\) and \(t=\tfrac{10}9\), so \(h=t-s=\tfrac19>0\). The standardized increment \(Z=(B_{10/9}-B_1)/\sqrt{1/9}\) has the standard normal distribution. Its eighth moment is \(7!!=7\cdot5\cdot3\cdot1=105\). Equivalently, the eighth moment of the unstandardized increment is \[ \mathbb{E}|B_{10/9}-B_1|^8 =105\left(\tfrac19\right)^4 =\frac{105}{6561} =\frac{35}{2187}. \]

Check: Standardizing multiplies the increment by \(3\), so its eighth absolute moment is multiplied by \(3^8=6561\). Indeed, \(6561\cdot\tfrac{35}{2187}=3\cdot35=105\), as required for a standard normal variable.

From Moments to Increment-Probability Bounds

Proposition (Moment bound for a large increment): For \(0\le s<t\), an integer \(m\ge1\), and \(a>0\), Brownian increments satisfy \[ \mathbb{P}(|B_t-B_s|>a)\le \frac{(2m-1)!!(t-s)^m}{a^{2m}}. \] The inequality is also valid with the event \(|B_t-B_s|\ge a\) on the left.

Proof. Set \(Y=|B_t-B_s|^{2m}\), which is nonnegative. If \(|B_t-B_s|>a\), then \(Y>a^{2m}\). Therefore \[ a^{2m}\mathbf{1}_{\{|B_t-B_s|>a\}}\le |B_t-B_s|^{2m}. \] Taking expectations and using the even-moment formula gives \[ a^{2m}\mathbb{P}(|B_t-B_s|>a) \le \mathbb{E}|B_t-B_s|^{2m} =(2m-1)!!(t-s)^m. \] Division by \(a^{2m}>0\) proves the strict-event bound. Since the same pointwise inequality holds with the event \(|B_t-B_s|\ge a\), it proves the non-strict version as well. \(\square\)

For regularity arguments, the threshold is often chosen in relation to the interval length. Fix \(\alpha\ge0\) and \(c>0\), and take \(a=c(t-s)^\alpha\). Writing \(h=t-s>0\), the proposition yields \[ \mathbb{P}(|B_t-B_s|>c h^\alpha) \le \frac{(2m-1)!!}{c^{2m}}h^{m-2m\alpha}. \] Thus the exponent of \(h\) is \(m(1-2\alpha)\). If one compares increments on a time grid with spacing \(h\), there are on the order of \(h^{-1}\) neighboring intervals in a fixed time range. The exponent exceeds \(1\) precisely when the single-interval bound decays faster than this count grows. This exponent comparison motivates moment-based regularity methods, but it is not by itself a proof of a uniform pathwise bound.

Worked Example: A Threshold Bound at a Chosen Scale

Use the sixth moment, so \(m=3\) and \((2m-1)!!=5\cdot3\cdot1=15\). Choose \(\alpha=\tfrac14\), \(c=1\), and an interval of length \(h=\tfrac1{64}\). Then \[ \mathbb{P}(|B_t-B_s|>h^{1/4}) \le 15h^{3-6(1/4)} =15h^{3/2} =15\left(\tfrac1{64}\right)^{3/2} =\frac{15}{512}. \]

Check: Since \((\tfrac1{64})^{1/2}=\tfrac18\), we have \((\tfrac1{64})^{3/2}=\tfrac1{64\cdot8}=\tfrac1{512}\). The bound is therefore \(15/512\), which lies between \(0\) and \(1\) and is a valid probability bound. The power of \(h\) is \(3/2>1\), the feature of the estimate that can be useful when controlling many short increments.

Choosing the Moment Order

The Gaussian law supplies moments of every even order, so the moment order can be chosen to match a desired threshold exponent. If the target threshold is proportional to \(h^\alpha\), the exponent in the probability bound is \(m(1-2\alpha)\). For a fixed \(\alpha<\tfrac12\), the factor \(1-2\alpha\) is positive. One can therefore choose an integer \(m\) large enough that \(m(1-2\alpha)>1\). At \(\alpha=\tfrac12\), this exponent is \(0\) for every \(m\); for \(\alpha>\tfrac12\), it is negative. This describes the strength of this particular moment-bound method, not a proof of the sharp path regularity of Brownian motion.

Worked Example: Selecting an Order for a Given Exponent

Suppose the threshold is \(h^{1/3}\), so \(\alpha=\tfrac13\). With \(m=2\), the fourth-moment bound has exponent \[ m(1-2\alpha)=2\left(1-\tfrac23\right)=\tfrac23, \] which is not greater than \(1\). Choose instead \(m=4\), corresponding to the eighth moment. Its coefficient is \(7!!=7\cdot5\cdot3\cdot1=105\), and the exponent is \[ m(1-2\alpha)=4\left(1-\tfrac23\right)=\tfrac43. \] Consequently, for any \(0\le s<t\) and \(h=t-s>0\), \[ \mathbb{P}(|B_t-B_s|>h^{1/3}) \le 105h^{4/3}. \]

Check: Substituting \(m=4\) and \(a=h^{1/3}\) directly into the proposition gives \(105h^4/(h^{1/3})^8=105h^{4-8/3}=105h^{4/3}\), the same bound. The exponent exceeds \(1\), whereas the fourth-moment choice gave only \(2/3\).

Common Pitfall: Scaling the Variance Incorrectly

A frequent error is to use the time gap itself as the scale of a fourth moment. The variance of \(B_t-B_s\) is \(h=t-s\), but a fourth moment scales as the square of that variance, not as the variance alone. For example, if \(h=2\), the incorrect calculation \(3h=3\cdot2=6\) gives the wrong fourth moment. The correct calculation is \[ \mathbb{E}|B_t-B_s|^4=3h^2=3\cdot2^2=12. \] The error is visible by standardizing: dividing the correct fourth moment by \(h^2=4\) gives \(3\), the fourth moment of a standard normal variable; dividing the incorrect value \(6\) by \(4\) does not.

There is a second distinction to keep in view: an estimate for one fixed pair \(s,t\) is not automatically a statement holding simultaneously for every pair of times. Path regularity requires control across many intervals and, ultimately, a uniform argument. The moment formula provides the quantitative input; additional reasoning is needed to convert that input into a pathwise conclusion. It is also important to use \(h=t-s\ge0\), not \(t+s\): the increment variance depends on the length of the interval between the times.

Takeaway: For every integer \(m\ge1\), Brownian increments satisfy \(\mathbb{E}|B_t-B_s|^{2m}=(2m-1)!!(t-s)^m\). Markov's inequality turns this into a large-increment bound, and choosing \(m\) determines the decay exponent available for continuity and Hölder arguments.

Check Your Understanding

In each question, use the interval length \(h=t-s\) and verify any numerical moment by substituting it into the Gaussian moment formula.

  1. For \(t-s=\tfrac14\), compute \(\mathbb{E}|B_t-B_s|^4\).
  2. For \(t-s=3\), compute the sixth moment of the increment.
  3. Using the fourth moment, what power of \(h\) appears in the bound for \(\mathbb{P}(|B_t-B_s|>h^{1/4})\)?
  4. For threshold \(h^{2/5}\), find an integer \(m\) for which the exponent \(m(1-2\alpha)\) is greater than \(1\), where \(\alpha=\tfrac25\).
  5. Why does an estimate for each fixed pair of times not, by itself, prove a pathwise bound that holds for all pairs?