Geometric drift gives a total-variation rate proportional to (Meyn-Tweedie Thm 15.0.1)
ProvedMarkovChainCLT.ergodicWithRate_of_geoDriftConditionLet be a Harris ergodic Markov chain with transition kernel and invariant probability distribution . Suppose a measurable function satisfies the geometric drift condition towards a measurable small set : is integrable under every and, for some and some ,
Then the chain converges geometrically in total variation with a rate constant proportional to : there are and with
This is the direction "drift (5) geometrically ergodic" of Theorem 15.0.1 of Meyn and Tweedie (1993), restricted from the -norm to the total-variation norm, and it is the content of the source's Remark 1 that under the drift condition "we can take " in the rate bound (3).
Formalization Note Harris ergodicity is the mission's total-variation encoding (HarrisErgodic), which supplies the -irreducibility and aperiodicity needed for convergence to ; the conclusion is the platform's ErgodicWithRate with constant and rate , exactly the shape consumed by the mission's geometric ergodicity predicate. The -algebra of the state space is assumed countably generated (MeasurableSpace.CountablyGenerated X), the standing assumption of Meyn and Tweedie (1993, Section 3.1) on which the existence of small sets (their Theorem 5.2.2) rests; the mission's Theorem 1(i) milestone carries the same hypothesis.
import Definitions.Def_MarkovErgodicity import Definitions.Def_MarkovDriftMinorization open MeasureTheory ProbabilityTheory Filter open scoped ENNReal NNReal Topology ProbabilityTheory
theorem MarkovChainCLT.ergodicWithRate_of_geoDriftCondition {X : Type*} [MeasurableSpace X]
[MeasurableSpace.CountablyGenerated X]
(P : Kernel X X) [IsMarkovKernel P] (π : Measure X) [IsProbabilityMeasure π]
(hP : HarrisErgodic P π)
(V : X → ℝ) (hV : Measurable V) (hV1 : ∀ x, 1 ≤ V x)
(C : Set X) (hC : MeasurableSet C) (hsmall : IsSmallSet P C)
(d b : ℝ) (hd : 0 < d) (hdrift : GeoDriftCondition P V d b C) :
∃ R ρ : ℝ, 0 ≤ R ∧ 0 ≤ ρ ∧ ρ < 1 ∧
ErgodicWithRate P π (fun x => R * V x) (fun n => ρ ^ n) := by sorry