1. Introduction
Let P = P(x, A) be a Markov transition probability function on a general state space (S, Σ), with invariant probability measure m (i.e. m(·) = ∫S P(x, ·)dm(x)). Let Ω := Sℕ be the space of trajectories with σ-algebra 𝒜 := Σ⊗ℕ, and let ℙx be the probability measure on 𝒜 governing the chain with transition probability function P and initial distribution δx. The probability of the chain with initial distribution m is then ℙm = ∫S ℙx dm(x). By invariance of m, ℙm is shift invariant on (Ω, 𝒜). Let Xn be the projection of Ω on the nth coordinate. Then (Xn) on (Ω, 𝒜, ℙm) is a stationary Markov chain with state space S.
For 1 ≤ p < ∞ we denote by Lp(m) the Banach space {f : S → ℝ : ∫S | f |p dm < ∞}, and put .
We assume m ergodic for P, which means (by one of the equivalent definitions) that if f ∈ L2(m) satisfies f(x) = ∫S f(y)P(x, dy) m-a.e., then f is constant a.e. Then the chain is ergodic too, i.e. the shift θ on (Ω, 𝒜, ℙm), defined by θ(Xn)n∈ℕ = (Xn+1)n∈ℕ, is ergodic.
We say that a real centered satisfies the annealed CLT if in (Ω, ℙm) we have
We say that a real centered satisfies the L2-normalized CLT if
provided for sufficiently large n ∈ ℕ.We denote by P also the Markov operator defined as
for every bounded measurable f and every x ∈ S. By invariance of m, P extends to all L1(m) functions, and is a contraction of all Lp(m) spaces, 1 ≤ p ≤ ∞, meaning that it does not increase the norm of functions in these spaces. As previously mentioned, ergodicity implies that Pf = f ∈ Lp holds only for f constant. We denote by Pn the n-fold composition of the operator P, and by Ef := ∫S f dm the expectation (with respect to the probability measure m) of f ∈ Lp(m), p ≥ 1.Following the early work of Doeblin, many efforts were made to identify conditions on an ergodic Markov operator P with invariant measure m which would ensure that every centered f ∈ L2(m) satisfies the annealed CLT – an L2-global annealed CLT for the chain.
2. History
Nagaev ([21]) used the following condition of Dobrushin: there exist k ∈ ℕ and δ < 1
This condition implies uniform geometric ergodicity: supx ‖Pn(x, ·) − m‖TV ≤ Mρn for some M > 0 and 0 < ρ < 1. But the latter condition implies ‖Pn − E‖∞ → 0, which turns out to be equivalent to Doeblin’s condition; see [24, p. 213]. Ibragimov ([18]) used a strong mixing condition (φ-mixing), which also turns out to imply Doeblin’s condition. Davydov ([9], [10]) constructed a positive recurrent aperiodic chain with countable state space such that the CLT fails for some centered f ∈ L2(m).Theorem 1 (M. Rosenblatt, [24]).
If ‖Pn − E‖2 → 0, then every centered f ∈ L2(m) satisfies the annealed CLT.
Rosenblatt proved that his condition is equivalent to ρ-mixing of the chain, and gave examples that it yields neither ‖Pn − E‖∞ → 0 nor ‖Pn − E‖1 → 0, although each of these conditions implies it; but ‖Pn − E‖2 → 0 if and only if ‖Pn − E‖p → 0 for some (every) 1 < p < ∞. Importantly, Rosenblatt’s condition does not necessarily imply Harris recurrence, see an example below.
Example (Random walks on the unit circle 𝕋).
Let μ be a probability measure on 𝕋, and define the convolution operator Pf = μ ∗ f, f ∈ L1(𝕋, m), m the normalized Haar (Lebesgue) measure. It is shown in [11] that if , that is, the Fourier transform of μ vanishes at infinity (i.e. μ is Rajchman), then ‖Pn − E‖2 → 0. When μ is Rajchman with all its powers singular with respect to Lebesgue measure, P is not Harris recurrent.
A contraction T on a Banach space 𝒳 is called uniformly ergodic if converges in the operator norm. The limit is a projection onto Fix(T) := {f ∈ 𝒳 : Tf = f} corresponding to the decomposition . A contraction T is uniformly ergodic if and only if (I−T)𝒳 is closed in 𝒳 ([20]).
When P is uniformly ergodic in L2(m), we have . (Recall that ). If ‖Pn − E‖2 → 0, then P is uniformly ergodic on L2(m); moreover, the spectral radius , meaning P has a spectral gap in the complex .
Theorem 2 (Gordin-Lifshits, [15]).
Let P be a Markov operator with invariant probability measure m, and assume that P is ergodic.
If f ∈ (I − P)L2(m), then f satisfies the annealed CLT, with
where f = (I − P)g with .
When (which is the case when P∗P is ergodic), f satisfies also the L2-normalized CLT, which follows from a theorem of Slutsky ([25]) (see [8, p. 254]).
By [7], f ∈ (I − P)L2(m) if and only if .
Theorem 1 now follows from Corollary 3 below.
Corollary 3.
Let P be a Markov operator with invariant probability measure m, and assume that P is uniformly ergodic in L2(m) with limit equal to E. Then every satisfies the annealed CLT.
Note that uniform ergodicity does not necessarily imply Harris recurrence.
3. Some ergodic properties
Theorem 4 (Derriennic-Lin, [11]).
Let P be a Markov operator with invariant probability measure m, and assume P is ergodic. Then the following conditions are equivalent:
(i) P is uniformly ergodic in L2(m).
(ii) For every we have .
(iii) For every we have .
(iv) For every we have .
Note that P is a contraction also of each complex Lp(m) space, 1 ≤ p ≤ ∞, and it is uniformly ergodic in the complex Lp(m) iff it is uniformly ergodic in the real Lp(m). A similar statement holds also for norm convergence of Pn.
Theorem 5.
Let P be a Markov operator with invariant probability measure m. If P is uniformly ergodic on Lp(m), 1 ≤ p < ∞, and is weakly mixing on the complex Lp(m) (the only unimodular eigenvalue of P is 1), then ‖Pn − E‖p → 0.
The proof primarily relies on positivity and ergodicity.
Lemma 6.
If P∗P is ergodic, then for every we have Pn f → 0 weakly in L2(m); thus the shift θ on (Ω, 𝒜, ℙm) is weakly mixing, hence totally ergodic (all powers θk are ergodic). Moreover, ‖(P∗P)n f ‖2 → 0 for every if and only if P∗P is ergodic.
Proof
We assume that P∗P is ergodic. Let 𝒦 be the unitary space of P:
Clearly if and only if . Hence, by the Cauchy-Schwarz inequality, g ∈ 𝒦 implies P∗Pg = g, and the ergodicity of P∗P implies that 𝒦 contains only the constant functions. Any f centered is therefore orthogonal to 𝒦, and by [13] both Pn f → 0 and P∗n f → 0 weakly in L2(m). Thus P is weakly mixing.The weak mixing of P implies that the shift θ is weakly mixing; see [1, Section 2].
The operator P∗P is symmetric positive semi-definite in the complex L2(m), so its spectrum is a subset of [0, 1]. If P∗P is ergodic, then for centered f ∈ L2(m) we have ‖(P∗P)n f ‖2 → 0 by the spectral theorem.
Conversely, if ‖(P∗P)n f ‖2 → 0 for every centered f ∈ L2(m), then obviously P∗P is ergodic.
4. Global central limit theorems
Theorem 8.
Let P be a Markov operator with invariant probability measure m. If P∗P is ergodic and P is uniformly ergodic, then ‖Pn − E‖2 → 0, and every centered 0 ≠ f ∈ L2(m) satisfies a non-degenerate annealed CLT and the L2-normalized CLT.
Moreover, if 0 ≠ f ∈ L3(m) is centered, then
Proof
Ergodicity of P∗P implies ergodicity of P, by Lemma 6. The assumption of uniform ergodicity implies that every is of the form f = (I − P)g with g ∈ L2(m) centered.
Fix 0 ≠ f = (I − P)g with g ∈ L2(m) centered. By the Gordin-Lifshits CLT, the annealed CLT holds for f, with variance of the limit expressed as
Hence σf = 0 if and only if P∗Pg = g. If σf = 0, then g is constant by the ergodicity of P∗P. Since g is centered, σf = 0 implies g = 0, so f = 0.By Lemma 6 the shift is totally ergodic, so Lemma 7 yields σn(f) > 0 for n ≥ 1. Thus, for centered f ≠ 0 we have n−1/2σn(f) → σf > 0, so the annealed CLT implies the L2-normalized CLT, by Slutsky’s theorem [25].
Ergodicity of P∗P implies weak mixing of P (Lemma 6), so uniform ergodicity yields ‖Pn − E‖2 → 0 (by Theorem 5). For 0 ≠ f ∈ L3(m) centered σf > 0 as shown above, and (1) holds by [17].
Corollary 9.
Let P be a Markov operator with invariant probability measure m, and assume that P is ergodic and uniformly ergodic. Every centered 0 ≠ f ∈ L2(m) satisfies a non-degenerate annealed CLT if and only if P∗P is ergodic.
Proof
When P∗P is ergodic Theorem 8 applies. For the converse, if P∗Pg = g for non-constant g ∈ L2(m), then P∗P(g − Eg) = g − Eg, and f = (I − P)(g − Eg) ≠ 0 satisfies the CLT with σf = 0.
Proposition 10.
Let P be a Markov operator with invariant probability measure m, and assume that P is normal in L2(m), i.e. P∗P = PP∗. If ‖Pn − E‖2 → 0, then P∗P is ergodic, and Theorem 8 applies.
Proof
Let P∗Pg = g ∈ L2(m). Since P∗E = E, normality yields ‖g − Eg‖2 = ‖(P∗P)n g − Eg‖2 = ‖P∗nPng − P∗nEg‖2 ≤ ‖Png − Eg‖2 → 0.
Example
In general, ‖Pn − E‖2 → 0 does not imply that P∗P is ergodic.
Let us define P on S := {1, 2, 3} by the matrix . The invariant probability vector is , and P∗ is given by the adjoint matrix. P has no non-trivial invariant sets, its only unimodular eigenvalue is 1, but P∗P is not ergodic.
Problem 2.
If a Markov operator P is ergodic, and every centered nonzero f ∈ L2(m) satisfies a non-degenerate annealed CLT, does ‖Pn − E‖2 → 0?
Note that P∗P is ergodic (proof of Corollary 9), so P is weakly mixing.
Below we present a sufficient “moment improving” condition for uniform ergodicity (called hyperboundedness); this condition is sometimes easy to check.
Theorem 11 (Glück, [14]).
Let P be a Markov operator with invariant probability measure m, assumed to be ergodic. Assume that for some 1 ≤ s < r < ∞ we have P Ls(m) ⊂ Lr(m). Then P is uniformly ergodic in all Lp(m) spaces, 1 < p < ∞ (i.e. ); hence (by Corollary 3) every centered f ∈ L2(m) satisfies the annealed CLT.
Example (A hyperbounded Markov operator).
Let (S, m) be the unit circle with normalized Lebesgue measure. Let 0 ≤ g ∈ L2(m) with ∫ g dm = 1, and define P by P f = g ∗ f. Then m is invariant, P is ergodic and normal in L2(m). Since ‖P f ‖2 = ‖g ∗ f ‖2 ≤ ‖g‖2‖f ‖1 for f ∈ L1(m), P maps L1(m) into L2(m).
Proposition 12 (Becker, [2])).
A power-bounded operator T (i.e. supn≥0 ‖Tn‖ < ∞) on a Banach space 𝒳 is uniformly ergodic if and only if for every the series ∑n≥1 n−1Tnf converges in 𝒳.
Proposition 13.
Let P be a Markov operator with invariant probability measure m, assumed to be ergodic. Then the following conditions are equivalent:
(i) The Markov chain is ρ-mixing1.
(ii) ‖Pn − E‖2 → 0.
(iii) For every the series converges.
(iv) For every we have .
(v) There exists 1 ≤ p < ∞ such that for every there exists r > 1 with .
If either of the above conditions holds, then the annealed CLT holds for every . The variance of the limiting normal distribution is
Proof
The equivalence of (i) and (ii) is by [24, p. 207].
By [11, Proposition 3.1], condition (ii) is equivalent to the existence of ρ < 1 and M > 0 such that ‖Pn − E‖2 ≤ M ρn for n ≥ 1. This yields (iii) and (iv).
(iii) implies uniform ergodicity, by Theorem 4. By [13, Lemma 2.1], (iii) implies Pn f → 0 weakly in L2(m) for every ; hence P is weakly mixing. Now (ii) holds by Theorem 5.
Obviously (iv) implies (v) with p = 2.
If (v) holds, then for every centered f ∈ Lp(m), Hölder’s inequality, applied with s = r/(r − 1), yields
Hence the series is convergent in Lp-norm when f ∈ Lp(m) is centered. By Becker’s Proposition 12, P is then uniformly ergodic in Lp(m). Since condition (v) implies that P has no unimodular eigenvalues, we have ‖Pn − E‖p → 0 (by Theorem 5), and by [24, Theorem VII.4.1] (ii) holds.Finally, (ii) implies the CLT statement by Theorem 1. By Theorem 2 the variance of the limit is limn→∞ σn(f)2/n.
5. α-mixing
Rosenblatt in [24] introduced a certain “strong mixing” condition, now called α-mixing, and proved that for the stationary chain generated by P with invariant probability measure m, α-mixing is equivalent to
The above supremum is bounded by ‖Pn − E‖2, so ρ-mixing implies α-mixing. Clearly α-mixing implies ‖Pn g − Eg‖2 → 0 for every g ∈ L2(m), hence total ergodicity of the shift θ.A stationary Markov chain which is Harris recurrent and aperiodic is α-mixing; see [4, Section 3.2].
Theorem 15.
Let P be a Markov operator with invariant probability measure m, and assume that the chain is α-mixing. If every satisfies the L2-normalized CLT, then P∗P is ergodic, every satisfies a non-degenerate annealed CLT, and ‖Pn − E‖2 → 0.
Proof
By Proposition 14 P∗P is ergodic, so the shift is totally ergodic. Hence for , σn(f) > 0 for every n ≥ 1, by Lemma 7.
Let γ ∈ (0, 1) be fixed. Fix , and put σn = σn(f). Since the chain is α-mixing, the stationary sequence {f(Xj)} is also α-mixing. By a result in [19], the L2-normalized CLT implies that there exists a function L(t), t > 0, slowly varying at ∞, such that . By a property of slowly varying functions, we obtain . Then
The above convergence holds for every . Denoting ϵ = (1 − γ)/2, we apply it to f = g − Eg, g ∈ L2(m), to obtain By the Banach-Steinhaus theorem, the norms are bounded, so . Thus P is uniformly ergodic. Theorem 8 yields ‖Pn − E‖2 → 0 and the non-degenerate annealed CLT for every .Theorem 16.
Let P be an ergodic Markov operator with invariant probability measure m. Then the following conditions are equivalent:
(i) ‖Pn − E‖2 → 0 and P∗P is ergodic.
(ii) The chain is α-mixing and every satisfies the L2-normalized CLT.
(iii) Every satisfies a non-degenerate annealed CLT and the L2-normalized CLT.
Proof
(i) implies (ii) follows from Theorem 8 and the fact that ρ-mixing implies α-mixing (combined with Proposition 13).
(ii) implies (i): Indeed, P∗P is ergodic by Proposition 14, and ‖Pn − E‖2 → 0 by Theorem 15.
(i) implies (iii) by Theorem 8.
(iii) implies (i): First of all, P∗P is ergodic by Proposition 14. Further, fix . We shall prove that is bounded. For the sake of contradiction, suppose it is not bounded. Then there is an increasing sequence {nk}k∈ℕ such that converges to zero, whence
The left-hand side of (2) converges in distribution to 𝒩(0, 1) by the assumption of the L2-normalized CLT for f; the right-hand side converges to 𝒩(0, 0), by the assumed annealed CLT for f and Slutsky’s theorem, leading to a contradiction. Hence is bounded for every . By Theorem 4, P is uniformly ergodic. By Proposition 14, P∗P is ergodic, so P is weakly mixing by Lemma 6, and then ‖Pn − E‖2 → 0 by Theorem 5.Problem 3.
Assume that P is a Makov operator with invariant probability measure m such that
and assume that every non-zero satisfies the L2-normalized CLT. Does it follow that P is uniformly ergodic in L2(m)?If yes, then ‖Pn − E‖2 → 0 by Theorem 5, since P is weakly mixing by the strong convergence of Pn. Note that the assumption implies that P∗P is ergodic, by Proposition 14.
By Theorem 15, the answer is yes for P which is Harris recurrent and aperiodic.
Example (P not uniformly ergodic with (P∗P) ergodic).
Let Q be ergodic with invariant probability measure m which is not uniformly ergodic. For ε ∈ (0, 1) define P = Pε := εI +(1 − ε)Q. We shall prove that P∗P is ergodic. Clearly m is invariant also for P and for P∗P. For A ∈ Σ we have
If P∗P1A = 1A a.e., then for almost every x ∉ A the above summands are zero, so in particular Q1A ≤ 1A a.e. Since m is invariant, Q1A = 1A, and A is trivial by the ergodicity of Q. By definition (I − P)L2(m) = (I − Q)L2(m), so when Q is not uniformly ergodic (I − P)L2(m) is not closed; hence P is not uniformly ergodic.6. Geometric ergodicity
Definition.
A Markov operator P with invariant probability measure m is called geometrically ergodic if, for some ρ < 1,
Geometric ergodicity implies aperiodic Harris recurrence and α-mixing, with the α-mixing coefficients α(n) converging to 0 exponentially fast; see [4, Section 3.2].
Theorem 17 (Doukhan-Massart-Rio, [12]).
Let Σ be countably generated and let P be a geometrically ergodic Markov operator. Then any centered f with ∫ | f |2 log+ | f | dm < ∞ satisfies the annealed CLT.
Theorem 18 (Roberts-Tweedie, [23]).
Let Σ be countably generated, and let P be a Harris positive recurrent Markov chain. If ‖Pn − E‖2 → 0, then P is geometrically ergodic.
Note that ‖Pn − E‖2 → 0 does not necessarily imply Harris recurrence; therefore Harris recurrrence must be assumed.
Note.
The converse may fail – in [3] and [16] are examples of P geometrically ergodic with some centered f ∈ L2(m) which does not satisfy the annealed CLT, so limn→∞ ‖Pn − E‖2 > 0.
Theorem 19.
Let P be a Markov operator with invariant probability measure m, and assume that P is normal in L2(m). Then ‖Pn − E‖2 → 0 if (and only if) the α-mixing coefficients converge to zero (at least) exponentially fast.
Bradley ([5]) proved the theorem when P is symmetric.
In general, if P is geometrically ergodic, then P is Harris aperiodic and the α-mixing coefficients converge to zero exponentially fast. We do not know if a Harris aperiodic P whose α–mixing coefficients converge to zero exponentially fast is geometrically ergodic.
Corollary 20.
Let Σ be countably generated. If a Markov operator P is geometrically ergodic, and is additionally normal in L2(m), then
The symmetric case is in [22]. For S countable Corollary 20 is established in [26].
Remarks.
P in Theorem 19 need not be Harris recurrent.
When Σ is countably generated and P is Harris recurrent and normal in L2(m), Theorems 18 and 19 yield that exponential decay to 0 of α(n), geometric ergodicity and ρ-mixing are equivalent.
In Bradley’s and Häggström’s examples P is geometrically ergodic, and every centered f ∈ Lp(m), p > 2, satisfies the CLT, by Theorem 17; however, P does not have a spectral gap in Lp(m), i.e. limn→∞ ‖Pn − E‖p > 0, since otherwise it would imply limn→∞ ‖Pn − E‖2 = 0 ([24]), and so the CLT for every centered f ∈ L2(m). By Corollary 20, P in such examples cannot be normal in L2(m).
The examples of Bradley and Häggström show that without normality Theorem 19 fails, although we have geometric ergodicity.
Problem 4.
Let P be a Harris aperiodic Markov chain, and suppose that every centered f such that ∫ | f |2 log+ | f | dm < ∞ satisfies the annealed CLT. Does this imply that P is geometrically ergodic? (Is a converse of Theorem 17 true?).
Dedecker informed the author that an example of Bradley ([6]) exhibits P Harris recurrent which is not geometrically ergodic, such that every , p > 2, satisfies the annealed CLT. In Problem 4 we (necessarily) assume more, i.e. that the annealed CLT is satisfied by a strictly larger subset of .
Acknowledgements
The author is grateful to the referee for his careful reading of the manuscript and for his detailed comments which improved the presentation.
The author thanks the Department of Mathematics of the University of Silesia in Katowice for its warm hospitality.