Degenerate case :
ProvedMarkovChainCLT.stationary_clt_of_asymptoticVariance_nonposThe degenerate branch of the Markov-chain CLT: zero asymptotic variance.
Let be a Harris ergodic Markov kernel with invariant probability measure , let , write and for the centered functional along the stationary chain, and suppose the covariance series converges absolutely, so that the asymptotic variance
is well defined. If , then
Jones (2004) states Theorems 5–8 under the standing assumption , which is the only case of practical interest: it is what makes the limit a genuine Gaussian and the Monte Carlo error bars meaningful. This statement is the complementary branch, and it is exactly what a formalization needs in order to conclude the mission's unconditional SatisfiesCLT predicate, which permits a degenerate limit (convergence to the point mass at ).
Mathematically the content is that under stationarity and absolute convergence of the covariance series one has ; when this forces in , hence in probability, hence in distribution to . Note always holds under these hypotheses, so "" is really the single case ; it is phrased as an inequality only so that it is the exact logical complement of the hypothesis carried by the mission's other CLT statements, letting the two branches be combined without a side condition.
Formalization Note seqAsymptoticVariance is the mission's ; .toNNReal sends the degenerate value to , so the limit law gaussianReal 0 0 is the Dirac mass at . Indexing starts at , matching the mission's convention.
import Definitions.Def_MarkovErgodicity import Definitions.Def_MarkovChainPathMeasure import Definitions.Def_MixingCoefficients open MeasureTheory ProbabilityTheory Filter open scoped ENNReal NNReal Topology ProbabilityTheory
theorem MarkovChainCLT.stationary_clt_of_asymptoticVariance_nonpos {X : Type*} [MeasurableSpace X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hP : HarrisErgodic P π) (f : X → ℝ) (hf : Measurable f) (hL2 : MemLp f 2 π)
(hsum : Summable (fun k : ℕ => ∫ ω, (f (ω 0) - ∫ x, f x ∂π) *
(f (ω (k + 1)) - ∫ x, f x ∂π) ∂(chainMeasure P π)))
(hzero : seqAsymptoticVariance (chainMeasure P π)
(fun i ω => f (ω i) - ∫ x, f x ∂π) ≤ 0) :
TendstoInDistribution
(fun (n : ℕ) (ω : ℕ → X) => (Real.sqrt n)⁻¹ *
∑ i ∈ Finset.range n, (f (ω i) - ∫ x, f x ∂π))
atTop (id : ℝ → ℝ) (fun _ => chainMeasure P π)
(gaussianReal 0 (seqAsymptoticVariance (chainMeasure P π)
(fun i ω => f (ω i) - ∫ x, f x ∂π)).toNNReal) := by sorry