Stochastic homogenization and ergodicity
Explain when a random stationary environment has a deterministic effective response and what random correctors actually require.
On this page
The question for today
What you’ll understand
- Distinguish stationarity from ergodicity and their different conclusions.
- State the stationary-gradient and sublinearity properties of a stochastic corrector.
- Calculate effective conductivities in nonergodic and ergodic one-dimensional examples.
Periodicity gives a cell that can be repeated exactly. Random media usually have no such finite cell, yet large regions can still have a reproducible effective response. The replacement is a probability law together with spatial translations. Stationarity makes the law consistent across locations; ergodicity makes translation-invariant macroscopic quantities deterministic. These assumptions must be separated. A material selected once from two homogeneous possibilities is stationary, but looking at a larger specimen never averages away that initial choice.
Translations organize the random environment
Let a probability space carry a measurable measure-preserving action τx of the translation group. A stationary coefficient field is represented as a(x,ω)=a(0,τxω), with the appropriate jointly measurable realization. Assume a is symmetric and satisfies fixed positive lower and upper ellipticity bounds almost surely, uniformly in location. Ergodicity means every event invariant under all translations has probability zero or one. It does not mean that coefficients at separated locations are independent, nor does it provide a numerical speed for spatial averaging.
For this stationary ergodic uniformly elliptic scalar divergence-form setting, a deterministic homogenized matrix governs almost every realization. On a fixed bounded Lipschitz physical domain, the zero-Dirichlet solutions for a(x/ε,ω) and fixed H⁻¹ load converge weakly in H¹ and strongly in L² to the homogenized solution. Without ergodicity, effective coefficients can remain random and depend on the invariant information in the environment. Qualitative convergence alone supplies neither a central limit theorem nor a universal rate in ε.
Corrector gradients are stationary; potentials need not be
A stochastic corrector solves −div a(p+∇φp)=0 in the whole space in the weak sense. Its gradient is stationary, has mean zero, and has finite second moment. The potential is defined only up to an additive constant. In dimensions where point evaluation is unavailable for H¹ functions, normalization should use a spatial average rather than a claimed value at the origin. The potential itself need not be a stationary square-integrable random field; treating that stronger property as automatic loses important low-dimensional behavior.
The useful large-scale property is sublinearity after subtracting a spatial mean: the mean-square oscillation of φp on a large ball, divided by the square of its radius, tends to zero almost surely. It ensures the rescaled potential correction becomes negligible at the macroscopic scale while its gradient still modifies the flux. Constructing this object involves the stationary potential-field space and ergodic arguments. The formula is not obtained by imposing periodic boundary conditions on one arbitrary random sample and declaring that finite sample exact.
One-dimensional resistance makes the role of ergodicity visible
In one dimension, a fixed flux produces local slope proportional to resistance r=1/a. For a stationary ergodic field, the spatial mean of resistance converges almost surely to its expectation, yielding the reciprocal of that expectation as the effective conductivity. If the environment has invariant random information, the spatial limit instead involves the corresponding conditional expectation. The arithmetic expectation of conductivity is generally the wrong series response, even when the random law has a very simple description.
For independent cell values on an integer lattice, stationarity initially refers to integer shifts. Adding an independent uniformly distributed origin phase in a unit cell produces a field stationary under all real translations; its translation action is ergodic when the underlying cell process is ergodic. This detail matters when applying a theorem formulated for continuous translations. Random finite-cell approximations can be useful computationally, but representativeness, boundary effects, and ensemble uncertainty must be analyzed separately from the existence of an infinite-volume effective coefficient.
Work it through
A stationary material whose effective coefficient stays random
Choose a random variable C equal to one or four with equal probability and set a(x,ω)=C(ω) for every x. Does stationarity imply a deterministic conductivity?
Use the identity translation action τxω=ω. The field is stationary because translating space changes no value or probability law.
The event C=1 is translation invariant and has probability one half. The action is therefore not ergodic; its invariant information includes the selected material.
In each realization the coefficient is already constant at every scale, so its effective conductivity is C itself. It is one in half the realizations and four in the others, rather than the ensemble arithmetic average 5/2 or the reciprocal mean resistance 8/5.
Stationarity permits a random effective response. Ergodicity removes invariant environmental choices and is the hypothesis behind the deterministic conclusion.
The living notebook
Explore the model
Interactive lab · markov
Compare mixing and invariant alternatives
Reset, set transition p=0.5 and q=0.5, then Step once; predict probability of state one 0.5. For a fresh run Reset, set p=0 and q=0, then Step once; predict it stays zero from the fixed state-zero initial law. The second chain admits every stationary mixture, while the first has a unique stationary law. The lab evolves a two-state distribution, not a spatial random coefficient or a sample trajectory.
- P(state 1) at step 0
- 0%
- Stationary P(state 1)
- 60%
- Chain behavior
- Monotone convergence
Exact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. 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 forA stationary law and the loss of invariant alternatives are distinct properties.
Pause & reason
Check your understanding
Your turn
Try it yourself
Let independent unit cells have conductivity one or four with equal probability, and add an independent uniform origin shift to obtain continuous stationarity. Compute the almost-sure one-dimensional effective conductivity and explain why the worked example does not share it.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Average resistance along the ergodic environment before taking its reciprocal; compare the invariant information in the two constructions.
Compare with the solution
The resistance is one or one quarter, so its expectation is (1+1/4)/2=5/8. The ergodic spatial average of resistance converges to 5/8, giving effective conductivity 8/5. Independent cell values supply genuine spatial averaging; the random origin makes the law invariant under real translations. In the worked example every cell equals the same random C, so observing more cells reveals no new independent variation and the invariant event C=1 persists. Its spatial mean resistance is 1/C, giving effective conductivity C rather than 8/5.
Keep the idea
Take it with you
- Stationarity describes translation-invariant laws; ergodicity makes invariant macroscopic responses deterministic.
- Stochastic corrector gradients can be stationary even when their potentials are not.
- Qualitative ergodic homogenization supplies no universal convergence or fluctuation rate.
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.