advanced34 min studyLesson 804 in Mathematics

Mean ergodic theorem and unitary operators

Prove norm convergence through invariant vectors and telescoping coboundaries, while keeping unitary and merely isometric dynamics distinct.

On this page

What you’ll understand

  • State the unitary mean ergodic theorem as strong Hilbert-space convergence.
  • Derive the invariant-vector and coboundary decomposition.
  • Compute a finite averaging operator and an explicit oscillatory average.

The pointwise ergodic theorem follows individual trajectories almost everywhere. The mean ergodic theorem instead studies averages as vectors in a Hilbert space, proving convergence in norm. Its unitary version has a short geometric proof: invariant vectors remain unchanged, while differences of a vector and its shift telescope. The remaining vectors are limits of such differences. In a probability-preserving invertible system the limiting projection is again conditional expectation on invariant information, so the two theorems identify the same object through different convergence statements.

Averages of a unitary operator

Let U be a unitary operator on a real or complex Hilbert space. Average the first N powers of U, beginning with identity. The von Neumann mean ergodic theorem says these averages applied to any fixed vector converge in Hilbert norm to its orthogonal projection onto the fixed-vector subspace. This is strong operator convergence, meaning convergence for each vector, rather than convergence in operator norm. The theorem does not supply one uniform rate simultaneously for every unit vector.

For an invertible bimeasurable probability-preserving transformation, composing a square-integrable observable with T gives a unitary Koopman operator. Its fixed vectors are the square-integrable invariant observables. The orthogonal projection onto them is conditional expectation on the invariant sigma-algebra. A noninvertible preserving map still gives an isometry, but surjectivity and hence unitarity may fail. There is a broader mean ergodic theorem for contractions, but the unitary proof below is stated under the actual unitary hypothesis.

AN=1N∑k=0N−1Uk,ANf⟶PFix⁡(U)f in Hilbert norm,U unitaryA_N=\frac1N\sum_{k=0}^{N-1}U^k,\qquad A_Nf\longrightarrow P_{\operatorname{Fix}(U)}f\text{ in Hilbert norm},\qquad U\text{ unitary}
For Koopman operators in the invertible probability setting, this limit is the L2 conditional expectation on invariant information. Norm convergence and pointwise convergence are distinct claims.↔ Scroll wide equations horizontally.

Coboundaries telescope

Call a vector of the form g minus Ug a coboundary. Its average equals g minus the Nth shift of g, divided by N, because all intermediate terms cancel. Unitarity bounds its norm by twice the norm of g divided by N. Every fixed vector is orthogonal to every coboundary. Conversely, a vector orthogonal to the range of identity minus U is fixed by the adjoint; for a unitary operator this is equivalent to being fixed by U.

Thus the Hilbert space splits orthogonally into fixed vectors and the closure of the coboundary range. The averaging operators have norm at most one, so convergence on actual coboundaries extends by approximation to that closure. This proves the theorem. Closure cannot generally be omitted: a vector can be approximated by coboundaries without being one with a controlled transfer vector g. Consequently the explicit inverse-N bound for one coboundary does not become a universal inverse-N rate for arbitrary observables.

H=Fix⁡(U)⊕ran⁡(I−U)‾,AN(g−Ug)=g−UNgN,∥AN(g−Ug)∥≤2∥g∥N\mathcal H=\operatorname{Fix}(U)\oplus\overline{\operatorname{ran}(I-U)},\qquad A_N(g-Ug)=\frac{g-U^Ng}{N},\qquad \|A_N(g-Ug)\|\le\frac{2\|g\|}{N}
The orthogonal complement calculation uses unitarity. The bound applies to a specified coboundary with a known transfer vector, not uniformly to every vector in the closure.↔ Scroll wide equations horizontally.

Oscillatory modes can average out without mixing

If a nonzero vector is an eigenvector with a unit-modulus eigenvalue different from one, its average is a scalar geometric sum times that vector. The sum divided by N tends to zero, even though the individual shifted vectors keep their original norm forever. This is a clean example of norm convergence of Cesàro averages without convergence of individual iterates. A circle rotation supplies exactly such eigenmodes, and its preserved oscillations will obstruct weak mixing in a later lesson.

On a finite cycle, averaging over a full period projects exactly onto the constant functions on that cycle. With several cycles, the projection is constant separately on each one rather than globally. These finite examples also clarify why a spatial norm includes the probability weights: the squared L2 distance is an average over starting states, not an unweighted Euclidean sum unless the measure was chosen that way. A convergence statement should specify both its norm and its underlying invariant law.

A two-step average on a four-cycle

On a uniform four-cycle, let Uf(j)=f(j+1(mod4))Uf(j)=f(j+1\pmod4) and f=(1,3,−1,1)f=(1,3,-1,1). Compute the two-step average, its L2 distance to the invariant projection, and the full-period average.

  1. The invariant projection is the constant mean one. The two-step average (f+Uf)/2(f+Uf)/2 is (2,1,0,1)(2,1,0,1), whose difference from the projection is (1,0,−1,0)(1,0,-1,0).

  2. Under uniform probability the squared L2 distance is (1+0+1+0)/4=1/2(1+0+1+0)/4=1/2, so the distance is 1/21/\sqrt2. Four successive shifts include every coordinate, giving A4f=(1,1,1,1)A_4f=(1,1,1,1) exactly.

  3. The centered vector f−1=(0,2,−2,0)f-1=(0,2,-2,0) equals g−Ugg-Ug for g=(0,0,−2,0)g=(0,0,-2,0), whose L2 norm is one. The coboundary bound at N equal to two is one, consistent with the actual distance 1/21/\sqrt2.

The result

Finite-period averaging is an explicit orthogonal projection. The norm uses the normalized invariant probability and the coboundary estimate has a checkable, nonsharp constant.

Explore the model

Interactive lab

Interactive lab · waves

Two phase contributions can cancel partially

Pause and Reset; set both phase difference 120120 degrees and amplitude ratio 11. Predict resultant amplitude ∣1+e2πi/3∣=1|1+e^{2\pi i/3}|=1. Reset and set both phase 00 degrees and ratio 11, predicting amplitude two. The wave model adds two equal-frequency contributions; compute longer geometric sums on paper to model an eigenmode's Cesàro average. It does not implement a general unitary operator or certify norm convergence. Reset clears both parameters and time; edits alone retain the clock.

Two phase contributions can cancel partiallyTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). 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 forAveraging can suppress an oscillatory component even though each shifted component retains its original norm.

Check your understanding

What does strong convergence in the mean ergodic theorem assert?

Try it yourself

For rotation by one third on the circle, take f(x)=e2πixf(x)=e^{2\pi ix}. Compute the L2 norm of its five-term average and compare it with the geometric-sum upper bound.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The eigenvalue is a primitive cube root omega. Reduce five consecutive powers using the sum of the first three roots.

Compare with the solution

The average is f(x)(1+ω+ω2+1+ω)/5=f(x)(1+ω)/5f(x)(1+\omega+\omega^2+1+\omega)/5=f(x)(1+\omega)/5. Since ∥f∥2=1\|f\|_2=1 and ∣1+ω∣=1|1+\omega|=1, its norm is one fifth. The general geometric bound is 2/(5∣1−ω∣)=2/(53)2/(5|1-\omega|)=2/(5\sqrt3), larger than one fifth. This observable averages to zero although the rational rotation is not ergodic on Lebesgue probability.

Take it with you

  • The unitary mean theorem gives strong Hilbert-norm convergence of averages.
  • Fixed vectors survive and coboundaries telescope.
  • Averaging out an eigenmode does not imply mixing or ergodicity of the whole system.

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.