Joint spectrum and large deviation principle for random
matrix products

arXiv:1702.06937v1 [math.PR] 22 Feb 2017

Cagri SERT
Abstract
The aim of this note is to announce some results about the probabilistic and deterministic asymptotic properties of linear groups. The first one is the analogue, for norms
of random matrix products, of the classical theorem of Cramér on large deviation principles (LDP) for sums of iid real random variables. In the second result, we introduce
a limit set describing the asymptotic shape of the powers S n = {g1 . . . . .gn | gi ∈ S}
of a subset S of a semisimple linear Lie group G (e.g. SL(d, R)). This limit set has
applications, among others, in the study of large deviations.

1
1.1

Large deviation principle for random matrix products
Introduction

Let G be a connected semisimple linear real algebraic group (e.g. SL(d, R)). A random
walk on G is a random process Yn = Xn . . . . .X1 where Xi ’s are independent and identically
distributed (iid) G-valued random variables. Starting from Bellman [8], Furstenberg-Kesten
[14] and Furstenberg [13], an important aim in the study of these non-commutative random walks was to establish the analogues of the classical limit theorems existing for the
iid real random variables. More precisely, one is interested in studying the probabilistic
limiting behavior of the (logarithms of the) norms log ||Yn || of random matrix products.
In fact, we will consider the slightly more general multi-norm given by the Cartan projection that we briefly describe before going on: let g be the Lie algebra of the group G,
a a Cartan subalgebra of g and a+ a chosen Weyl chamber. Let K be a maximal compact subgroup of G for which we have the Cartan decomposition G = K exp(a+ )K. Then
the map κ : G −→ a+ is well-defined by the following equality and is called the Cartan projection: for all g ∈ G, g = k exp(κ(g))u for some k, u ∈ K. As an example, in
the case of G = SL(d, R), the Cartan projection κ(g) of a matrix g ∈ SL(d, R) writes as
2
||∧d g||
k d
κ(g) = (log ||g||, log ||∧||g||g|| , . . . , log ||∧
d−1 g|| ), where ∧ R ’s are considered with their canonical Euclidean structures and ||.||’s denote the associated operator norms. The components
of κ(g) are the logarithms of the singular values of g.
The first limit theorem that was proven for random matrix products was the analogue
of the law of large numbers. Stating it in our setting, Furstenberg-Kesten’s
result [14] reads:
R
if µ is a probability measure on G with a finite first moment (i.e. ||κ(g)||µ(dg) < ∞ for
some norm ||.|| on a), then the µ-random walk Yn = Xn . . . . .X1 (i.e. Xi ’s are iid of law µ)
satisfies
1
a.s.
κ(Yn ) −→ ~λµ ∈ a
n→∞
n
where ~λµ can be defined by this and called the Lyapunov vector of µ. Nowadays, this
result is a rather straightforward corollary of Kingman’s subadditive ergodic theorem. A
1

1

LARGE DEVIATION PRINCIPLE FOR RANDOM MATRIX PRODUCTS

2

second important limit theorem that was established in increasing generality by Tutubalin
[26], Le Page [19], Goldsheid-Guivarc’h [16], and Benoist-Quint [6] [5] is the central limit
theorem (CLT). Benoist-Quint’s CLT reads: if µ is a probability measure on G of finite
second order moment and such that the support of µ generates a Zariski-dense semigroup
in G, then √1n (κ(Yn ) − n~λµ ) converges in distribution to a non- degenerate Gaussian law on
a. A feature of this result is the Zariski density assumption which also appears in our result
below. We note that the fact that the support S of the probability measure µ generates
a Zariski-dense semigroup can be read as: any polynomial that vanishes on ∪n≥1 S n also
vanishes on G (we recall that in R, a subset is Zariski dense if and only if it is infinite). Some
other limit theorems whose analogues have been obtained are the law of iterated logarithm
and local limit theorems for which we refer the reader to the nice books of Bougerol-Lacroix
[11] and more recently Benoist-Quint [5].
An essential and, up to our work, a rather incomplete aspect of these non-commutative
limit theorems is concerned with large deviations. The main result in this direction is that
of Le Page [19], (see also Bougerol [11]) and its extension by Benoist-Quint [5], stating the
exponential decay of large deviation probabilities off the Lyapunov vector. Before stating
this result, recall that a probability measure µ on
G is said to have a finite exponential
R ||κ(g)||
moment, if there exists α > 1 such that we have α
µ(dg) < ∞. We have
Theorem 1.1 (Le Page [19], Benoist-Quint [5]). Let G be as before, µ be a probability measure of finite exponential moment on G whose support generates a Zariski-dense semigroup
in G. Let Yn denote the nth step of the µ-random walk on G. Then, for all  > 0, we have
lim supn→∞ n1 log P(|| κ(Ynn ) − ~λµ || > ) < 0.

1.2

Statement of main result

In our first main result, under the usual Zariski density assumption, we prove the matrix analogue of a classical theorem (see below) about large deviations for iid real random
variables. Let X be a topological space and F be a σ- algebra on X.
Definition 1.2. A sequence Zn of X-valued random variables is said to satisfy a large
deviation principle (LDP) with rate function I : X −→ [0, ∞], if for every measurable
subset R of X, we have
− inf I(x) ≤ lim inf
◦

x∈R

n→∞

1
1
log P(Zn ∈ R) ≤ lim sup log P(Zn ∈ R) ≤ − inf I(x)
n
n→∞ n
x∈R

◦

where, R denotes the interior and R the closure of R.
With this definition,
the classical Cramér-Chernoff theorem says that the sequence of
P
averages Yn = n1 ni=1 Xi of real iid random variables of finite exponential moment satisfies
an LDP with a proper convex rate function I, given by the convex conjugate (FenchelLegendre transform) of the logarithmic moment generating function of Xi ’s. Our first main
result reads
Theorem 1.3. Let G be a connected semisimple linear real algebraic group and µ be a
probability measure of finite exponential moment on G, whose support generates a Zariski
dense sub- semigroup of G. Then, the sequence random variables n1 κ(Yn ) satisfies an LDP
with a proper convex rate function I : a −→ [0, ∞] assuming a unique zero on the Lyapunov
vector ~λµ of µ.

1

LARGE DEVIATION PRINCIPLE FOR RANDOM MATRIX PRODUCTS

3

Remark 1.4. 1. Without any moment assumptions on µ, we also obtain a weaker result
which is an analogue of a result of Bahadur [7] for iid real random variables.
2. Under a stronger exponential moment condition (sometimes called finite super-exponential
moment), by exploiting convexity of I, we are able to identify the rate function I with the
convex conjugate of a limiting logarithmic moment generating function of the random variables n1 κ(Yn ).
3. It follows by convexity of I that the effective support DI := {x ∈ a | I(x) < ∞} of the
rate function I is a convex subset of a+ . We also show that this set DI depends only on
the support S of µ (i.e. not on the particular mass distribution µ on S) and we identify
DI with a Hausdorff limiting set of a deterministic construction, namely the joint spectrum
J(S) of S (see Theorem 2.3 below).
4. We conjecture that a similar LDP holds for the Jordan projection λ : G → a+ in place
of κ (see Def. in paragraph 2.1).

1.3

The idea of the proof

The fact that I(~λµ ) = 0 follows from the definition of ~λµ and I. That this zero is unique
is merely the translation of the deep result of Le Page (Theorem 1.1). That I is proper
is a rather straightforward consequence of Chernoff’s estimates and that one can identify
I as a Fenchel-Legendre transform (2. of Remark 1.4) follows from its convexity using
standard techniques, namely the Varadhan’s integral lemma (see Theorem 4.3.1. in [12])
and elementary properties of convex conjugation of functions. Moreover, the convexity of
the rate function is proven using ideas similar to those used to prove its existence. Therefore,
in this note, we will focus on the proof of the existence of an LDP. For simplicity, we shall
assume that the measure µ is compactly supported. We make use of the following general
fact:
Theorem 1.5 (see Theorem 4.1.11 in [12]). Let X be a topological space endowed with its
Borel σ-algebra βX , and Zn be a sequence of X- valued random variables taking values in a
compact subset of X. Let A be a base of open sets for the topology of X. For each x ∈ X,
define:
Ili (x) := sup − lim inf
A∈A
x∈A

n→∞

1
log P(Zn ∈ A)
n

and

1
Ils (x) := sup − lim sup log P(Zn ∈ A)
n→∞ n
A∈A
x∈A

Suppose that for all x ∈ X, we have Ili (x) = Ils (x). Then, the sequence Zn satisfies an
LDP with rate function I given by I(x) := Ili (x) = Ils (x).
In view of this theorem, to prove the existence of the LDP in Theorem 1.3, we have
to show that the equality Ili = Ils is satisfied. A first remark is if κ was an additive
mapping (i.e. κ(gh) = κ(g) + κ(h)), this would follow rather easily from the independence
of random walk increments and uniform continuity of κ. A further important remark is that
in fact a weaker form of additivity (i.e. ||κ(gh) − κ(g) − κ(h)|| is uniformly bounded for all
g, h ∈ supp(µ)) is sufficient to insure this equality. A key result of Benoist [2], [3] shows that
this weak form of additivity is satisfied in a given (r, )-Schottky semigroup. This already
finishes the proof in the case when µ is supported on such a semigroup. For the general case,
we need an argument showing that we can restrict the random walk on Schottky semigroups
with no loss in the exponential rate of decay probabilities (as in Lemma 1.6). This is done

2

JOINT SPECTRUM

4

by using a result of Abels-Margulis-Soifer [1] about proximal elements in Zariski dense
semigroups (which in turn uses a result of Benoist-Labourie [9] and Prasad [20]) together
with the uniform continuity of the Cartan projection and then a simple partitioning and
pigeonhole argument. A key step in this proof is the following Lemma 1.6.
Recall that an element g in PGL(d, R) is called -proximal if it has a unique eigenvalue of
maximal modulus and the ratio of its first two singular values is at least 1 . It is called (r, )proximal if additionally the top eigenvector is at least r away from projective hyperplane
spanned by the other generalized eigenspaces. An element g in G is called (r, )-loxodromic
if it is (r, )-proximal in each of the rk(G) fundamental proximal representations of G.
Abels-Margulis-Soifer show that for a Zariski dense semigroup Γ in G, there exists r > 0
such that for every  > 0, one can find a finite subset F ⊂ Γ with the property that for all
γ ∈ Γ, there exists f ∈ F such that γ.f is (r, )-loxodromic. We denote by Γ the semigroup
generated by the support of µ.
Lemma 1.6. Let 0 <  < r = r(Γ). There exist a compact set C = C(Γ, ) ⊂ a, a natural
number i0 = i0 (, Γ, µ), and a constant d1 = d1 (, Γ, µ) > 0 such that for all n0 ∈ N and
subset R ⊂ a+ , there exists a natural number n1 ≥ n0 with n1 − n0 ≤ i0 such that we have
P(κ(Yn1 ) ∈ R + C and Yn1 is (r, )-loxodromic) ≥ d1 .P(κ(Yn0 ) ∈ R)
The next step in the proof consists in observing that one can further restrict the random
walk to an (r, )-Schottky semigroup, again losing only a uniformly bounded proportion of
the probability (as in Lemma 1.6). By doing so, we reduce the situation to a random walk
on a semigroup on which the Cartan projection κ(.) is almost additive and hence we can
conclude as we sketched in the beginning of the argument.
Remark 1.7. In a further work in preparation [25], using essentially the same method as
in the proof of Theorem 1.3, we show that the LDP holds for the average word length of
random walks on Gromov hyperbolic groups.

2
2.1

Joint spectrum
Introduction

In this second part of this note, we define the notion of joint spectrum of a bounded subset
of a semisimple Lie group. We then relate this notion to the effective support DI = {x ∈
a | I(x) < ∞} of the rate function from Theorem 1.3. Recall the definition of the Jordan
projection λ : G −→ a+ : if g = ge gh gu is the Jordan decomposition of g with ge elliptic, gh
hyperbolic and gu is unipotent, then λ(g) is defined as κ(gh ). Let S be a bounded subset
of G and S n = {g1 . . . . .gn | gi ∈ S} denote its nth power. We are interested in the following
questions: do the sequences of bounded subsets of a+ , n1 κ(S n ) and n1 λ(S n ) have a limit in
the Hausdorff topology ? If yes, are the limits the same and can one describe these limit
sets ?

2.2

Statement of the main results

Regarding the above questions, we show the following:
Theorem 2.1. Let G be a connected semisimple linear real algebraic group and S a bounded
subset of G generating a Zariski dense sub-semigroup. Then,

2

JOINT SPECTRUM

5

1. The following Hausdorff limits exist, and we have the equality:
1
1
κ(S n ) = lim λ(S n )
n→∞ n
n→∞ n
This common limit will be denoted as J(S) and called the joint spectrum of S.
lim

2. J(S) is a compact convex subset of a+ with non-empty interior.
Remark 2.2. 1. It is not hard to see that if two subsets S and S 0 of G generate the same
Zariski dense semigroup Γ in G, then the projective images of J(S) and J(S 0 ) are the same.
Therefore the corresponding cone in a+ only depends on Γ. It turns out that this is precisely
the Benoist limit cone of Γ [3].
2. For a Banach algebra (B, ||.||) and a bounded subset S of B, denote the joint spectral
radius of S by r(S) := limn→∞ n1 log(supx∈S n ||x||). When B is finite dimensional, this does
not depend on the particular norm on B. Let now ρ : G → SL(V ) be an irreducible rational
representation of G and let χρ ∈ a∗ denote the highest weight of ρ. Then we have the equality
supz∈J(S) χρ (z) = r(ρ(S)). This allows us to derive a multi-dimensional generalization of
Berger-Wang identity [10].
We now relate the joint spectrum of the support of a probability measure µ with the
effective support DI of the rate function I for the LDP of n1 κ(Yn ), where Yn denotes as
usual the µ-random walk on G.
Theorem 2.3. Let G and µ be as in Theorem 1.3, and suppose moreover that the support S
of µ is bounded. Let I be the rate function given by Theorem 1.3 and DI = {x ∈ a | I(x) <
∞} be the effective support of I. We then have
◦

◦

1. DI = J(S) and DI = J(S). If S is moreover finite, then DI = J(S).
2. The Lyapunov vector ~λµ ∈ a+ of µ belongs to the interior of J(S).
For the point 2. of this theorem, we note that for a probability measure as in Theorem
1.3, the fact that ~λµ belongs to the interior of a+ was obtained by Guivarc’h-Raugi [18]
and Goldsheid-Margulis [15]. Our result gives a more precise location for ~λµ in case µ is,
moreover, boundedly supported. Our method extends also to the case where µ only has a
finite exponential moment to show that ~λµ belongs to the Benoist cone of the semigroup
generated by the support of µ.
The tools that go into the proof of Theorem 2.1 and Theorem 2.3 are mostly similar to
those used in the proof of Theorem 1.3. We use an additional tool to prove the fact that the
joint spectrum is of non-empty interior: while this result can be directly deduced from the
properties of Benoist cone proved in [3], we adopt an indirect approach and use the central
limit theorem of Goldsheid-Guivarc’h [16] and Guivarc’h [17] that we combine with AbelsMargulis-Soifer result [1] and Benoist’s estimates [2]. This in turn allows us to derive point
2. of the previous theorem and hence also part of the point 2. of Theorem 2.1.
Remark 2.4. The notion of joint spectrum plays also an important role in another work
in preparation [23], where we study a new exponential counting function for a finite subset
S in a group G as before.

Acknowledgements
These results are part of author’s doctoral thesis realized under the supervision of Emmanuel
Breuillard in Université Paris-Sud during 2013-2016. The author would like to take the

opportunity to thank him for numerous discussions and for suggesting the notion of joint
spectrum. The author also thanks to WWU Münster where part of this work was conducted.

References
[1] Abels, H., Margulis, G.A. and Soifer, G.A., 1995. Semigroups containing proximal
linear maps. Israel journal of mathematics, 91(1-3), pp.1-30.
[2] Benoist, Y., 1996. Actions propres sur les espaces homogènes réductifs. Annals of
mathematics, pp.315-347.
[3] Benoist, Y., 1997. Propriétés asymptotiques des groupes linéaires. Geometric &
Functional Analysis GAFA, 7(1), pp.1-47.
[4] Benoist, Y., 2000, Propriétés asymptotiques des groupes linéaires II. Advanced Studies
Pure Math. 26 (2000) pp.33-48.
[5] Benoist, Y. and Quint, J.F., 2013. Random walks on reductive groups. to appear.
[6] Benoist, Y. and Quint, J.F., 2016. Central limit theorem for linear groups. Ann.
Probab, 44(2), pp.1308-1340.
[7] Bahadur, R.R., 1971. Some limit theorems in statistics (No. 04; QA276. A1, B3 1971.).
Philadelphia: Society for industrial and applied mathematics.
[8] Bellman, R., 1954. Limit theorems for non-commutative operations. I. Duke Mathematical Journal, 21(3), pp.491-500.
[9] Benoist, Y., and Labourie, F., 1993. Sur les difféomorphismes d’Anosov affines à feuilletages stable et instable différentiables. Inventiones mathematicae, 111(1), pp.285-308.
[10] Berger, M.A. and Wang, Y., 1992. Bounded semigroups of matrices. Linear Algebra
and its Applications, 166, pp.21-27.
[11] Bougerol, P. and Lacroix J., 1985. Products of random matrices with applications to
Schrödinger operators. Progress in Probability and Statistics, 8. Birkhäuser Boston,
Inc., Boston, MA, xii+283 pp.
[12] Dembo, A. and Zeitouni, O. Large deviations techniques and applications (Vol. 38).
Springer Science & Business Media. 2009
[13] Furstenberg, H. Noncommuting random products. Transactions of the American
Mathematical Society, 1963, 108(3), pp.377-428.
[14] Furstenberg, H. and Kesten, H. Products of random matrices. The Annals of Mathematical Statistics, 1960, 31(2), pp.457-469.
[15] Goldsheid, I. and Margulis, G. Lyapunov Indices of a Product of Random Matrices,
Russian Math. Surveys 44, 1989, pp.11-81.
[16] Goldsheid, I.Y. and Guivarc’h, Y. Zariski closure and the dimension of the Gaussian
low of the product of random matrices. I. Probability theory and related fields, 1996,
105(1), pp.109-142.
[17] Guivarc’h, Y. On the spectrum of a large subgroup of a semisimple group. Journal of
modern dynamics, 2008, 2(1), p.15.
[18] Guivarc’h, Y. and Raugi, A. Frontière de Furstenberg, propriétés de contraction
et théorèmes de convergence. Probability Theory and Related Fields, 1985, 69(2),
pp.187-242.
[19] Le Page, E. Théorèmes limites pour les produits de matrices aléatoires. In Probability
measures on groups. Springer Berlin Heidelberg, 1982, pp. 258-303.

REFERENCES

7

[20] Prasad, G. R-regular elements in Zariski dense subgroups. The Quarterly Journal of
Mathematics, 1994, 45(4), pp.542-545.
[21] Sert C. Joint spectrum and large deviation principle for random matrix products,
Phd. thesis, Université Paris-Sud, 2016.
[22] Sert C. Large deviation principle for random matrix products (in preparation).
[23] Sert C. Growth indicator on semisimple linear groups (in preparation).
[24] Sert C. Joint spectrum on semisimple linear groups (in preparation).
[25] Sert C. Large deviation principle and growth indicator for Gromov hyperbolic groups
(in preparation).
[26] Tutubalin, V.N. A central limit theorem for products of random matrices and some of
its applications. In Symposia Mathematica Vol. 21, 1977, pp.101-116.
E-mail address: cagri.sert@math.ethz.ch
Telephone number: +41 44 632 3737

