Skip to main content
Have a personal or library account? Click to login
Global Central Limit Theorems for Stationary Markov Chains Cover

Global Central Limit Theorems for Stationary Markov Chains

By:   
Open Access
|May 2025

Full Article

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 = Sx 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 L0pm=fLpm :Sfdm=0  .

We assume m ergodic for P, which means (by one of the equivalent definitions) that if fL2(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 fL02m satisfies the annealed CLT if in (Ω, ℙm) we have

1nk=1nfXk𝒟𝒩0, σ2,where𝒩0,0:=δ0.

We say that a real centered fL02m satisfies the L2-normalized CLT if

1σnfk=1nfXk𝒟𝒩0,1,
provided σnf:=k=1nfXkL2m>0 for sufficiently large n ∈ ℕ.

We denote by P also the Markov operator defined as

Pfx:=SfyPx,dy
for every bounded measurable f and every xS. 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 = fLp 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 fLp(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 fL2(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

supx,ySPkx,APky,A<δ,AΣ.
This condition implies uniform geometric ergodicity: supxPn(x, ·) − mTVn for some M > 0 and 0 < ρ < 1. But the latter condition implies ‖PnE → 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 fL2(m).

Theorem 1 (M. Rosenblatt, [24]).

IfPnE2 → 0, then every centered fL2(m) satisfies the annealed CLT.

Rosenblatt proved that his condition is equivalent to ρ-mixing of the chain, and gave examples that it yields neither ‖PnE → 0 nor ‖PnE1 → 0, although each of these conditions implies it; but ‖PnE2 → 0 if and only ifPnEp → 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, fL1(𝕋, m), m the normalized Haar (Lebesgue) measure. It is shown in [11] that if limkμ^k=0 , that is, the Fourier transform of μ vanishes at infinity (i.e. μ is Rajchman), then ‖PnE2 → 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 1nk=1nTk converges in the operator norm. The limit is a projection onto Fix(T) := {f𝒳 : Tf = f} corresponding to the decomposition X=FixTITX¯ . A contraction T is uniformly ergodic if and only if (IT)𝒳 is closed in 𝒳 ([20]).

When P is uniformly ergodic in L2(m), we have L02m=IPL2m=I PL02m . (Recall that L02m:=fL2:Ef=0 ). If ‖PnE2 → 0, then P is uniformly ergodic on L2(m); moreover, the spectral radius rP|L02m<1 , meaning P has a spectral gap in the complex L02m .

Theorem 2 (Gordin-Lifshits, [15]).

Let P be a Markov operator with invariant probability measure m, and assume that P is ergodic.

If f ∈ (IP)L2(m), then f satisfies the annealed CLT, with

σ2=σf2:=limn1nk=1nfXk22=g2Pg2,

where f = (IP)g with gL02m .

When σf2>0 (which is the case when PP is ergodic), f satisfies also the L2-normalized CLT, which follows from a theorem of Slutsky ([25]) (see [8, p. 254]).

By [7], f ∈ (IP)L2(m) if and only if supn k=1nPkf2< .

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 fL02m satisfies the annealed CLT.

Note that uniform ergodicity does not necessarily imply Harris recurrence.

Problem 1.

Let P be a Markov operator with invariant probability measure m, and assume that P is ergodic. If every fL02m satisfies the annealed CLT, does it follow that P is uniformly ergodic in L2(m)?

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 fL02m we have supn1 1nk=1nfXkL2m2< .

  • (iii) For every fL02m we have supn1 1nk=1nPkf2< .

  • (iv) For every fL02m we have supn1 k=1nPkf,f< .

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), thenPnEp → 0.

The proof primarily relies on positivity and ergodicity.

Lemma 6.

If PP is ergodic, then for every fL02m 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, ‖(PP)n f2 → 0 for every fL02m if and only if PP is ergodic.

Proof

We assume that PP is ergodic. Let 𝒦 be the unitary space of P:

K:=gL2m:Png2=P*ng2=g2foreveryn1.
Clearly Pg22=g22 if and only if P*Pg,g=g22 . Hence, by the Cauchy-Schwarz inequality, g𝒦 implies PPg = g, and the ergodicity of PP implies that 𝒦 contains only the constant functions. Any f centered is therefore orthogonal to 𝒦, and by [13] both Pn f → 0 and Pn 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 PP is symmetric positive semi-definite in the complex L2(m), so its spectrum is a subset of [0, 1]. If PP is ergodic, then for centered fL2(m) we have ‖(PP)n f2 → 0 by the spectral theorem.

Conversely, if ‖(PP)n f2 → 0 for every centered fL2(m), then obviously PP is ergodic.

Lemma 7.

Let the shift θ be totally ergodic on (Ω, 𝒜,m), which is the case when PP is ergodic. If f ≠ 0 belongs to L02m , then σn(f) > 0 for every n ≥ 1.

Proof

By stationarity of the chain (Xn), σn(f) = 0 implies

k=0n1fXkL2m=0,
so
fX0θnfX0=fXn+k=0n1fXkfX0=k=1nfXk=k=0n1fXkθ=0.
By ergodicity of θn, f(X0) is a constant, which is zero since f is centered.

4. Global central limit theorems

Theorem 8.

Let P be a Markov operator with invariant probability measure m. If PP is ergodic and P is uniformly ergodic, thenPnE2 → 0, and every centered 0 ≠ fL2(m) satisfies a non-degenerate annealed CLT and the L2-normalized CLT.

Moreover, if 0 ≠ fL3(m) is centered, then

(1)
suptmk=1nfXkσfnt12πtex2/2dx=O1n.

Proof

Ergodicity of PP implies ergodicity of P, by Lemma 6. The assumption of uniform ergodicity implies that every fL02m is of the form f = (IP)g with gL2(m) centered.

Fix 0 ≠ f = (IP)g with gL2(m) centered. By the Gordin-Lifshits CLT, the annealed CLT holds for f, with variance of the limit expressed as

σ2=σf2=limn1nk=1nfXkL2m2=g22Pg22,gL02m.
Hence σf = 0 if and only if PPg = g. If σf = 0, then g is constant by the ergodicity of PP. 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 PP implies weak mixing of P (Lemma 6), so uniform ergodicity yields ‖PnE2 → 0 (by Theorem 5). For 0 ≠ fL3(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 ≠ fL2(m) satisfies a non-degenerate annealed CLT if and only if PP is ergodic.

Proof

When PP is ergodic Theorem 8 applies. For the converse, if PPg = g for non-constant gL2(m), then PP(gEg) = gEg, and f = (IP)(gEg) ≠ 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. PP = PP. IfPnE2 → 0, then PP is ergodic, and Theorem 8 applies.

Proof

Let PPg = gL2(m). Since PE = E, normality yields ‖gEg2 = ‖(PP)n gEg2 = ‖PnPngPnEg2 ≤ ‖PngEg2 → 0.

Example

In general, ‖PnE2 → 0 does not imply that PP is ergodic.

Let us define P on S := {1, 2, 3} by the matrix 1212000112120 . The invariant probability vector is 13,13,13 , and P is given by the adjoint matrix. P has no non-trivial invariant sets, its only unimodular eigenvalue is 1, but PP is not ergodic.

Problem 2.

If a Markov operator P is ergodic, and every centered nonzero fL2(m) satisfies a non-degenerate annealed CLT, doesPnE2 → 0?

Note that PP 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. 1nk=1nPkEp0 ); hence (by Corollary 3) every centered fL2(m) satisfies the annealed CLT.

Example (A hyperbounded Markov operator).

Let (S, m) be the unit circle with normalized Lebesgue measure. Let 0 ≤ gL2(m) with g dm = 1, and define P by P f = gf. Then m is invariant, P is ergodic and normal in L2(m). Since ‖P f2 = ‖gf2 ≤ ‖g2f1 for fL1(m), P maps L1(m) into L2(m).

Proposition 12 (Becker, [2])).

A power-bounded operator T (i.e. supn0Tn‖ < ∞) on a Banach space 𝒳 is uniformly ergodic if and only if for every fITX¯ the series ∑n1 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)PnE2 → 0.

  • (iii) For every fL02m the series k=1Pkf,f converges.

  • (iv) For every fL02m we have n=1Pnf22< .

  • (v) There exists 1 ≤ p < ∞ such that for every fL0pm there exists r > 1 with n=1Pnfpr< .

If either of the above conditions holds, then the annealed CLT holds for every fL02m . The variance of the limiting normal distribution is

σf2=f22+2k=1Pkf,f.

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 ‖PnE2M ρ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 fL02m ; 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 fLp(m), Hölder’s inequality, applied with s = r/(r − 1), yields

n=1Pnfpnn=11ns1sn=1Pnfpr1r<.
Hence the series n=1Pnfn is convergent in Lp-norm when fLp(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 ‖PnEp → 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.

Proposition 14.

Let P be a Markov operator with invariant probability measure m. If every 0fL02m satisfies the L2-normalized CLT, then PP is ergodic. Consequently (Lemma 6 and Theorem 5), if P is uniformly ergodic, thenPnE2 → 0.

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

4αn:=supfdm=0Pnf1f0asn.
The above supremum is bounded by ‖PnE2, so ρ-mixing implies α-mixing. Clearly α-mixing implies ‖Pn gEg2 → 0 for every gL2(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 0fL02m satisfies the L2-normalized CLT, then PP is ergodic, every 0fL02m satisfies a non-degenerate annealed CLT, andPnE2 → 0.

Proof

By Proposition 14 PP is ergodic, so the shift is totally ergodic. Hence for 0fL02m , σn(f) > 0 for every n ≥ 1, by Lemma 7.

Let γ ∈ (0, 1) be fixed. Fix 0fL02m , 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 σn2=nLn . By a property of slowly varying functions, we obtain nγ+1σn2=nγLn0 . Then

1nγ+1/2k=1nPkf21nγ+1/2k=1nfXkL2m=nγ+1/2σn0.
The above convergence holds for every fL02m . Denoting ϵ = (1 − γ)/2, we apply it to f = gEg, gL2(m), to obtain
nε1nk=1nPkgEg2=1nγ+1/2k=1nPkgEg2CgngL2m.
By the Banach-Steinhaus theorem, the norms nε1nk=1nPkE2 are bounded, so 1nk=1nPkE2Knε0 . Thus P is uniformly ergodic. Theorem 8 yields ‖PnE2 → 0 and the non-degenerate annealed CLT for every 0fL02m .

Theorem 16.

Let P be an ergodic Markov operator with invariant probability measure m. Then the following conditions are equivalent:

  • (i)PnE2 → 0 and PP is ergodic.

  • (ii) The chain is α-mixing and every 0fL02m satisfies the L2-normalized CLT.

  • (iii) Every 0fL02m 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, PP is ergodic by Proposition 14, and ‖PnE2 → 0 by Theorem 15.

(i) implies (iii) by Theorem 8.

(iii) implies (i): First of all, PP is ergodic by Proposition 14. Further, fix 0fL02m . We shall prove that σnf/n is bounded. For the sake of contradiction, suppose it is not bounded. Then there is an increasing sequence {nk}k such that nk/σnkf converges to zero, whence

(2)
1σnkfj=1nkfXj=nkσnkf1nkj=1nkfXj.
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 σnf/n is bounded for every fL02m . By Theorem 4, P is uniformly ergodic. By Proposition 14, PP is ergodic, so P is weakly mixing by Lemma 6, and then ‖PnE2 → 0 by Theorem 5.

Problem 3.

Assume that P is a Makov operator with invariant probability measure m such that

limnPngEg2=limnP*ngEg2=0foreverygL2m,
and assume that every non-zero fL02m satisfies the L2-normalized CLT. Does it follow that P is uniformly ergodic in L2(m)?

If yes, then ‖PnE2 → 0 by Theorem 5, since P is weakly mixing by the strong convergence of Pn. Note that the assumption implies that PP 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 (PP) 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 PP is ergodic. Clearly m is invariant also for P and for PP. For A ∈ Σ we have

P*P1A=ε21A+ε1εQ*1A+Q1A+(1ε)2Q*Q1A.
If PP1A = 1A a.e., then for almost every xA 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 (IP)L2(m) = (IQ)L2(m), so when Q is not uniformly ergodic (IP)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,

Mx:=supnρnPnx,mTV<a.e.

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. IfPnE2 → 0, then P is geometrically ergodic.

Note that ‖PnE2 → 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 fL2(m) which does not satisfy the annealed CLT, so limn→∞PnE2 > 0.

Theorem 19.

Let P be a Markov operator with invariant probability measure m, and assume that P is normal in L2(m). ThenPnE2 → 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

PnE20.

The symmetric case is in [22]. For S countable Corollary 20 is established in [26].

Remarks.

  1. P in Theorem 19 need not be Harris recurrent.

  2. 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 fLp(m), p > 2, satisfies the CLT, by Theorem 17; however, P does not have a spectral gap in Lp(m), i.e. limn→∞PnEp > 0, since otherwise it would imply limn→∞PnE2 = 0 ([24]), and so the CLT for every centered fL2(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 fL0pm , 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 L02m .

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.

DOI: https://doi.org/10.2478/amsil-2025-0011 | Journal eISSN: 2391-4238 (formerly 0860-2107) | Journal ISSN: 0860-2107
Language: English
Page range: 177 - 189
Submitted on: Feb 6, 2025
Accepted on: Apr 23, 2025
Published on: May 20, 2025
Published by: University of Silesia in Katowice, Institute of Mathematics
In partnership with: Paradigm Publishing Services

© 2025 Michael Lin, published by University of Silesia in Katowice, Institute of Mathematics
This work is licensed under the Creative Commons Attribution 4.0 License.