Uniform mixing makes the transition operator an contraction on mean-zero functions
ProvedMarkovChainCLT.integral_sq_iterKernel_le_of_tvDistA uniformly ergodic chain contracts . If the -step kernel satisfies the uniform total-variation bound , then for every square-integrable with ,
Why this is the bridge from bounded to square-integrable observables. For a bounded , uniform ergodicity gives directly, and everything — the Poisson equation, the martingale approximation, the central limit theorem — follows easily. For a merely square-integrable that estimate is unavailable, and the sup-norm must be replaced by the norm. Interpolating naively fails; what works is the Cauchy–Schwarz bound against total variation, which costs a square root of the mixing rate but keeps the norm on the right-hand side.
Since can be made as small as one likes by taking large (uniform ergodicity gives ), choosing makes a strict contraction of the mean-zero subspace , so that decays geometrically. This is what makes the Neumann series converge in — solving the Poisson equation for square-integrable data — and what yields the summable covariances behind the variance bound for partial sums.
Proof. For -almost every the measure integrates (because by invariance), so the square-root total-variation inequality applies with , . Using ,
Squaring and using ,
and integrating in against — where invariance turns back into — gives the factor .
import Definitions.Def_MarkovErgodicity import Definitions.Def_MarkovIterKernel import Definitions.Def_TotalVariationDist import Mathlib.Probability.Kernel.Invariance import Mathlib.MeasureTheory.Integral.Bochner.Set open Filter Function MeasurableSpace MeasureTheory ProbabilityTheory open MarkovChainCLT open scoped ENNReal NNReal Topology
theorem MarkovChainCLT.integral_sq_iterKernel_le_of_tvDist {X : Type*} [MeasurableSpace X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hinv : Kernel.Invariant P π) (N : ℕ) (ρ : ℝ) (hρ0 : 0 ≤ ρ)
(hρ : ∀ x, tvDist (iterKernel P N x) π ≤ ρ)
(h : X → ℝ) (hh : Measurable h) (hL2 : Integrable (fun x => (h x) ^ 2) π)
(hmean : ∫ x, h x ∂π = 0) :
∫ x, (∫ y, h y ∂(iterKernel P N x)) ^ 2 ∂π ≤ 4 * ρ * ∫ x, (h x) ^ 2 ∂π := by sorry