arxiv: v2 [math.sp] 1 Jun 2016

Transcription

1 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE OF HYPERBOLIC SURFACES EIENNE LE MASSON AN UOMAS SAHLSEN arxiv:65.57v [math.sp] Jun 6 Abstract. We present a quantum ergodicity theorem for fixed spectral window and sequences of compact hyperbolic surfaces converging to the hyperbolic plane in the sense of Benjamini and Schramm. his addresses a question posed by Colin de Verdière. Our theorem is inspired by results for eigenfunctions on large regular graphs by Anantharaman and the first-named author. It applies in particular to eigenfunctions on compact arithmetic surfaces in the level aspect, which connects it to a question of Nelson on Maass forms. he proof is based on a wave propagation approach recently considered by Brooks, Lindenstrauss and the first-named author on discrete graphs. It does not use any microlocal analysis, making it quite different from the usual proof of quantum ergodicity in the large eigenvalue limit. Moreover, we replace the wave propagator with renormalised averaging operators over discs, which simplifies the analysis and allows us to make use of a general ergodic theorem of Nevo. As a consequence of this approach, we only require bounded measurable observables.. Introduction.. Background and motivation. Let H be the hyperbolic plane, and X = Γ\H be a compact hyperbolic surface. Let be the Laplace-Beltrami operator on X. We consider an arbitrary orthonormal basis {ψ j } j N of L (X) consisting of eigenfunctions of, and denote by {λ j } j N the non-decreasing sequence of corresponding eigenvalues. he surface X is known to satisfy Quantum Ergodicity (QE), a spectral analog of the classical ergodic property, stating that N(λ) ψ j, a ψ j a d Vol λ j λ X, as λ, where a is any fixed continuous observable on X and the symbol means that we divide the integral by the total volume. Here N(λ) = {λ j λ} is the eigenvalue counting function (with multiplicity). More generally, one can replace multiplication by a in the inner product by a zero-order pseudodifferential operator A, and the a in the integral by the principal symbol σ A of A. he quantum ergodicity property implies that there exists a density subsequence (ψ jk ) such that the measures ψ jk d Vol converge weakly to the uniform measure. It was first shown to hold in this setting by Shnirelman and Zelditch [36, 39], and generalised to any Riemannian manifold with ergodic geodesic flow by Colin de Verdière [5]. Ergodicity of the geodesic flow is not enough to avoid having to extract a subsequence, as was shown by Hassell []. However, for negatively curved manifolds, Rudnick and Sarnak conjectured that the full sequence of eigenfunctions should equidistribute asymptotically [33]. his Quantum Unique Ergodicity (QUE) conjecture was proved in the setting of compact Mathematics Subject Classification. 8Q5, 374, F7. Key words and phrases. Quantum chaos, quantum ergodicity, Benjamini-Schramm convergence, short geodesics, hyperbolic dynamics, eigenfunctions of the Laplacian, Selberg transform, rate of mixing, mean ergodic theorem. E. Le Masson was partially supported by ERC advanced grants 6759 and 947 and the NSF grant for MSRI.. Sahlsten was partially supported by the ERC starting grant and Marie Sk lodowska-curie Individual Fellowship grant 6553.

2 EIENNE LE MASSON AN UOMAS SAHLSEN arithmetic hyperbolic surfaces of congruence type by Lindenstrauss, for joint eigenfunctions of the Laplacian and an algebra of Hecke operators, arising from the arithmetic structure of the surface ([8], see also [3] for a strengthening of this result using only one Hecke operator). By excluding the possibility of an escape of mass in the cusp, Soundararajan [37] completed Lindenstrauss s result for the non-compact case of the modular surface SL (Z)\H. An alternative approach was pursued by Anantharaman [4] (see also [6]), who proved entropy bounds for eigenfunctions on Riemannian manifolds with Anosov geodesic flow, excluding in particular the possibility of concentration of mass on closed geodesics. We will be interested here in a different point of view, where the eigenvalues stay bounded in a fixed interval, and we take instead a large scale geometric limit by considering sequences of hyperbolic surfaces converging to the plane. he convergence we consider is in the sense of Benjamini and Schramm, which we will define precisely. his kind of setting for the study of eigenfunctions has attracted attention separately in the field of modular forms and on discrete graphs. Let us review some of these results to motivate our theorem. In the field of modular forms, and in parallel to Lindenstrauss s results, a theory of QUE for holomorphic forms was developed (see [34] for a survey), culminating in the work of Holowinsky and Soundararajan [3]. Note that these results are concerned with holomorphic Hecke eigenforms and are proved for specific non-compact arithmetic surfaces Γ (q)\h, q N, where { ( ) } a b Γ (q) = SL c d (Z) : c mod q, that is, congruence coverings of the modular surface. In the setting of holomorphic forms, two parameters are available: the weight, which is roughly analogous to the Laplace eigenvalue, and the level q, which characterises the congruence covering. he results of Holowinsky and Soundararajan concern the limit of large weights, and it is natural to ask what happens for a fixed weight in the level aspect, that is when the level goes to infinity (which means in particular that Γ (q)\h converges to H as q in the sense of Benjamini and Schramm, see [, heorem 5.]). Building on the methods of Holowinsky and Soundararajan, QUE-type theorems in the level aspect have been proved recently by Nelson [3], and Nelson, Pitale and Saha [3]. However, some important elements are missing to adapt these methods to Maass forms, that is eigenfunctions of the Laplacian, as was noted by Nelson [3, Remark.7]. Independently of the theory of modular forms, there has been a growing interest in the problem of quantum chaos on discrete regular graphs since the pioneering work of Jakobson, Miller, Rivin, and Rudnick [5] (see also [35] for a survey on the subject). Regular graphs can be seen as discrete analogues of hyperbolic surfaces. In this setting the eigenvalues of the Laplacian are bounded, and the relevant asymptotics is that of large graphs, in the sense of the number of vertices going to infinity. A first result of delocalisation of eigenfunctions was shown by Brooks and Lindenstrauss []. Using a discrete version of microlocal calculus [7], Anantharaman and Le Masson then proved a general quantum ergodicity result for regular graphs [5]. Since then, alternative proofs of this result have been proposed by Brooks, Le Masson and Lindenstrauss [], and Ananthraman [3]. In the light of these results on graphs, one can ask if a quantum ergodicity theorem can be proved on manifolds for large spatial scale and fixed eigenvalues. his problem was suggested to the first-named author by Yves Colin de Verdière and we are addressing it in this paper. In the special case of compact arithmetic surfaces, our result gives in particular an equidistribution result for eigenfunctions in the level aspect, connecting it to the question of equidistribution of Maass forms. Perhaps surprisingly, the adaptation of the pseudo-differential methods of [5] in the continuous setting presents some important difficulties, and the proof of our theorem will instead be inspired by the wave propagation methods proposed in [] for the discrete Laplacian on graphs. By introducing a continuous analogue of these methods, we are able

3 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 3 to dispense completely with any pseudo-differential calculus, which makes our proof quite different from the usual proof of the quantum ergodicity theorem. In contrast with the latter, where continuity of the observables is usually required, we only need L, or bounded measurable test functions depending on the case considered... Main results. o state our result, let us define the notion of Benjamini-Schramm convergence. his notion was introduced by Benjamini and Schramm as a model of distributional convergence of graphs [7], but it has a natural analogue in the continuous setting, where it was used in particular by Abert et al. [, ] and Bowen [9]. Let InjRad X (z) be the injectivity radius of X at z X and InjRad(X) = inf z X InjRad X (z) the minimal injectivity radius (see Section for a definition). We say that a sequence of compact hyperbolic surfaces X n = Γ n \H Benjamini-Schramm converges to the hyperbolic plane H if for any R >, Vol({z X n : InjRad lim Xn (z) < R}) =. (.) n Vol(X n ) he convergence means that the probability for a ball of radius R at a random point in X n to be isometric to a ball in H tends to as n. Note that this is automatically satisfied if InjRad(X n ). We will discuss some examples after the statement of the results. he theorem we prove is the following: heorem.. Fix a closed interval I (/4, + ). Suppose that (X n ) n N is a sequence of compact connected hyperbolic surfaces such that () (X n ) converges to H in the sense of Benjamini and Schramm; () here exists l min > such that for all n N, InjRad(X n ) l min. (3) here exists β > such that for all n N, the Laplacian on X n has a spectral gap β, i.e. the smallest eigenvalue of the Laplacian is greater than β; Let λ (n) λ (n)... be the sequence of eigenvalues of the Laplacian on X n, and {ψ (n) k } k N a corresponding orthonormal basis of eigenfunctions in L (X n ). For any uniformly bounded sequence of measurable functions (a n ) n N, we have N(X n, I) j : λ (n) j I ψ(n) j, a n ψ (n) a n d Vol Xn (.) when n +, where N(X n, I) is the number of eigenvalues in the interval I counted with multiplicity. he theorem means that on large hyperbolic surfaces, and for any fixed spectral interval, most of the eigenfunctions of the Laplacian approach the uniform measure in a weakconvergence sense. o be more precise, we need to discuss the type of sequences of test functions a n that allow such an interpretation. In order to get a true equidistribution, we want to ensure that cancellations occur between the terms ψ (n) j, a n ψ (n) j and an d Vol Xn. It means that the latter quantity should not tend to too fast compared to the rate of convergence of (.) when n +. From (.) we can deduce that cancellations will occur when the support of a n is of the order of Vol(X n ). For example, for any fixed constant < c <, the theorem says that we have equidistribution in a portion of the surface of area c Vol(X n ). By giving a quantitative estimate of the rate of convergence of (.), we will see that we actually have equidistribution on smaller scales. heorem. is a consequence of the following quantitative version, together with an asymptotic estimate of the number of eigenvalues in bounded intervals under Benjamini- Schramm convergence (see Section 9). heorem.. Let X = Γ\H be a compact connected hyperbolic surface. Let λ λ... be the sequence of eigenvalues of the Laplacian on X, and {ψ k } k N a corresponding orthonormal j

4 4 EIENNE LE MASSON AN UOMAS SAHLSEN basis of eigenfunctions in L (X). Fix a closed interval I (/4, + ), then there exists R I > such that for any a L (X) we have ψ a j, a ψ j a d Vol X I ϱ(β) R + e4r Vol({z X : InjRad l X (z) < R}) a min λ j I for all R > R I. Here ϱ(β) > is a constant depending only on the spectral gap of the Laplacian on X, and l min is the injectivity radius of X. In particular, if R < InjRad(X), we have ψ a j, a ψ j a d Vol X I ϱ(β) R, (.3) λ j I for any function a L (X). Here A I B means there exist a constant C > that only depends on I with A CB. As a consequence of this theorem, if we replace conditions () and () in heorem. with the assumption that the radius of injectivity go to infinity, then we get the same conclusion assuming only that the observables a n are in L and that a n = O(Vol(X n)). For simplicity, let us comment on the rate (.3) of the large injectivity radius case. Assuming a n = O(Vol(X n)), we obtain N(X n, I) j : λ (n) j I ψ(n) j, a n ψ (n) ( ) j a n d Vol Xn = O, R n where R n log(vol(x n )) is the injectivity radius. Note that this rate is analogous to the rate log(λ) obtained in the large eigenvalue case λ by Zelditch [4]. he rate in (.3) also gives equidistribution on smaller scales, for observables with support of area c Vol(X n )/R c Vol(X n )/ log(vol(x n )), where < c < is a fixed constant. In the large eigenvalue limit, the question of Quantum Ergodicity on small scales has been addressed in [, ]. In contrast with this case, because the quantitative estimate (.3) depends only on the L norm of the observable, we have no restriction on the size of the support of the observables and can choose observables with support on even smaller scales. In this case however, we can only deduce that there is no concentration of the mass of eigenfunctions at these very small scales, the rate is not strong enough to give equidistribution. he quantitative statements of heorem. show in particular that it would be possible to weaken the condition of uniformity of the spectral gap and shortest geodesic. However, the conditions we require are typical for hyperbolic surfaces, in the sense that they are satisfied almost surely in the large scale limit for a natural model of random surfaces []. he conditions are also satisfied for congruence coverings of compact arithmetic surfaces, where the level goes to infinity, making it possible to deduce immediately an equidistribution result for eigenfunctions in the level aspect in this setting (see for example [8] for an introduction to the arithmetic setting, and [, heorem 5.] for a proof of the Benjamini-Schramm convergence in this case). In general, any sequence of cocompact lattices (Γ n ) n N such that for all n N, Γ n+ Γ n, Γ n is normal in Γ, and n N Γ n = {id}, gives a sequence of surfaces X n = Γ n \H satisfying conditions () and () in heorem. (the radius of injectivity is tending to + in this case [6, heorem.])..3. Further discussions and generalisations. heorem. is a Quantum Ergodicity type theorem, in the sense that we average over a spectral interval. By analogy with the large eigenvalue case, it is natural to ask if a stronger property holds. Question.3. o sequences of surfaces converging in the sense of Benjamini-Schramm satisfy a QUE property for fixed eigenvalues intervals? hat is, if I (/4, + ) is a

5 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 5 compact interval, then under the assumptions of heorem., do we have max λ j I ψ(n) j, a n ψ (n) j a n d Vol Xn as (X n ) Benjamini-Schramm converges to H? As far as we know this problem is open also in a weaker form where instead of taking the maximum over an interval we choose any sequence (ψ (n) j n ) of eigenfunctions and ask if, a n ψ (n) a n d Vol Xn. ψ(n) j n j n Let us just mention again that an analog of this question for holomorphic forms has been addressed in [3] and [3]. he methods of the proof of heorem. and heorem. are relatively elementary. hey require only some kind of wave propagation and mixing properties of the dynamics. In this sense they have a potential to be generalised to a large range of different settings, as is already apparent from their use on graphs. We adapted the techniques to the specific case of hyperbolic surfaces. In particular we replace the wave propagator with renormalised averaging operators on discs of growing radius. his allows us to use an ergodic theorem of Nevo that becomes particularly natural in this context (see Section 6). wo extensions of the result seem attainable with similar methods: the case of Maass forms on congruence coverings of the modular surface, and general observables in phase space. Indeed, the compactness of the surfaces is essentially used to ensure that the spectrum is discrete, and Nevo s ergodic theorem is valid in phase space. A more challenging problem would be to adapt the methods to the case of variable curvature or more general homogeneous spaces, in particular higher dimension. In this case it is unclear how the choice of propagator and the estimate on its spectral action would generalise. We believe our result shows how the study of the discrete graph model can give new insights on the continuous case. here are also open problems on graphs that might benefit from the understanding of the manifold case, in particular the connection between non-regularity of the graphs and variable curvature. If the history of QUE is any indication, this result is just the beginning of a new approach to the understanding of eigenfunctions of the Laplacian. he differences with the usual setting are yet to be explored, and we hope our methods can provide in particular new insights to the large eigenvalue theory. One of the striking differences is the minimal regularity required for the observables: is it a characteristic of the large scale setting or is it a strength of our method that could be transported to the large eigenvalue case? In the other direction, is it possible to find a proof of our theorem using the usual microlocal analysis approach? he large eigenvalue version of the Quantum Ergodicity theorem only requires ergodicity of the geodesic flow, but we use in a crucial way some mixing properties of the dynamics and a sufficient decay of correlations. Our result is just a first step. It is difficult to predict the difficulty of the QUE problem in the large scale setting compared to the large eigenvalue one. We hope that this question will bring new perspectives and ideas to the field. Organisation of the paper. he paper is organised as follows. Section contains some preliminaries on harmonic analysis on hyperbolic surfaces. In Section 3 we give the general idea of the proof of the main theorems. Section 4 introduces a wave propagation operator and gives the proof of heorems. and. assuming two fundamental estimates: a Hilbert- Schmidt norm estimate and an estimate on the spectral action of the propagator. o prove the first estimate we introduce a general Hilbert-Schmidt norm estimate adapted to the case of Benjamini-Schramm convergence in Section 5. hen we invoke in Section 6 the main tool

6 6 EIENNE LE MASSON AN UOMAS SAHLSEN from ergodic theory we will be using, that is a mean ergodic theorem in L for averaging operators. he proof of the Hilbert-Schmidt norm estimate is then carried out in 7. he spectral action estimate is done in Section 8. Finally Section 9 introduces an asymptotic estimate of the number of eigenvalues in a fixed bounded interval which can be seen as an analog of the Weyl law for the Benjamini-Schramm convergence setting.. Elements of harmonic analysis and hyperbolic geometry In this section, we give some definitions and introduce elements of harmonic analysis on hyperbolic surfaces that we will use in the proof. For more background on the geometry and spectral theory of hyperbolic surfaces we refer to the books [4, 4, 8]... Hyperbolic surfaces. he hyperbolic plane is identified with the upper-half plane H = {z = x + iy C y > }, equipped with the hyperbolic Riemannian metric ds = dx + dy y. We will denote by d(z, z ) the distance between two points z, z H. he hyperbolic volume is given by dx dy dµ(z) = y. he group of isometries of H is identified with PSL(, R), the group of real matrices of determinant modulo ± id, acting by Möbius transformations (( ) ) ( ) a b a b PSL(, R), z H z = az + b c d c d cz + d. A hyperbolic surface can be seen as a quotient X = Γ\H of H by a discrete subgroup Γ PSL (R). We denote by a fundamental domain associated to Γ. If we fix z H, an example of a fundamental domain is given by the set = {z H d(z, z) < d(z, γz) for any γ Γ {± id}}. he injectivity radius on the surface X = Γ\H at a point z is given by InjRad X (z) = min{d(z, γz) γ Γ {id}}. hus InjRad X (z) gives the largest R > such that B X (z, R) is isometric to a ball of radius R in the hyperbolic plane. Let g PSL(, R), we define the translation operator g, such that for any function f on H g f(z) = f(g z). We will generally see a function f on a hyperbolic surface X = Γ\H as a Γ-invariant function f : H C, γ f(z) = f(γ z) = f(z) for all γ Γ. he integral of the function on the surface is then equal to the integral of the invariant function over any fundamental domain f(z) dµ(z).

7 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 7.. Geodesic flow. he tangent bundle of H can be identified with H C. he hyperbolic metric gives the following inner product for two tangent vectors (z, re iθ ) and (z, r e iθ ) on the tangent plane z H re iθ, r e iθ z = r r Im (z) cos(θ θ). As a consequence, the map (z, θ) H S (z, Im (z) e iθ ) H C, where S = R/πZ, identifies H S with the unit tangent bundle. he group PSL(, R) acts on the tangent bundle via the differential of its action on H. It is well known (see for example [6]) that this action induces a homeomorphism between PSL(, R) and the unit tangent bundle of H, such that the action of PSL(, R) on itself by left multiplication corresponds to the action of PSL(, R) on the unit tangent bundle. We denote by ϕ t : H S H S the geodesic flow associated to H. It is invariant with respect to the Liouville measure dµ dθ, where dθ is the Lebesgue measure on S. Via the identification H S PSL(, R), the geodesic flow is equal to the multiplication on the right by the diagonal subgroup ( ) e t/ ϕ t (g) = g e t/, g G, t R. For a hyperbolic surface Γ\H, the unit tangent bundle is identified with Γ\ PSL(, R), and via this identification the geodesic flow will be given simply by ( ) e t/ ϕ t (Γg) = Γg e t/..3. Polar coordinates. We will make use of the polar coordinates on H. Let z H be an arbitrary point. For any point z H different from z, there is a unique geodesic of length r going from z to z. Using the geodesic flow, it means that there is a unique θ S and r (, ) such that z is the projection of ϕ r (z, θ) on the first coordinate. he change of variable z (r, θ) is called polar coordinates. he induced metric is ds = dr + sinh r dθ, and the hyperbolic volume in these coordinates is given by dµ(r, θ) = sinh r dr dθ..4. Invariant integral operators and Selberg transform. In the coordinates z = x + iy, the Laplacian on H is the differential operator ( ) = y x + y. A fundamental property of the Laplacian is that it commutes with isometries. We have for any g PSL(, R), g = g. he Laplacian can therefore be seen as a differential operator on any hyperbolic surface Γ\H. We say that a kernel K : H H C is invariant under the diagonal action of Γ if for any γ Γ we have K(γ z, γ w) = K(z, w), (z, w) H H. For any Γ-invariant function f, such a kernel defines an integral operator A on the surface X by the formula Af(z) = K(z, w)f(w) dµ(w) = K(z, γw)f(w) dµ(w), z. H γ Γ

8 8 EIENNE LE MASSON AN UOMAS SAHLSEN he function k : H H C given by k(z, w) = γ Γ K(z, γw) is such that k(γz, γ w) = k(z, w) for any γ, γ Γ, which defines the kernel on the surface. A special case of invariant kernels is given by radial kernels. Let k : [, + ) C, then K(z, w) = k(d(z, w)), (z, w) H H is an invariant kernel. For k : [, + ) C, the Selberg transform S(k) of k is obtained as the Fourier transform S(k)(s) = + e isu g(u) du of the Abel transform g(u) = + k(ϱ) sinh ϱ dϱ. u cosh ϱ cosh u For a function h : R C, the Selberg transform is inverted using the inverse Fourier transform and the formula g(u) = π + e isu h(s) ds k(ϱ) = + g(u) du. π cosh u cosh ϱ ϱ Eigenfunctions of the Laplacian are eigenfunctions of all operators of convolution by a radial kernel and the eigenvalues are given precisely by the Selberg transform. Proposition.. Let X = Γ\H be a hyperbolic surface. Let k : [, + ) C be a bounded measurable function with compact support. If ψ λ is an eigenfunction of the Laplacian on X of eigenvalue λ, then it is an eigenfunction of the radial integral operator A associated to k. hat is, Aψ λ (z) = k(d(z, w))ψ λ (w)dµ(w) = h(s λ )ψ λ (z), where the eigenvalue h(s λ ) is given by the Selberg transform of the kernel k: h(s λ ) = S(k)(s λ ), and s λ C is defined by the equation λ = 4 + s λ. 3. Strategy of the proof Proving our theorem is equivalent to estimating ψ j, a ψ j λ j I for an observable a L, such that a dµ =. () For this purpose we introduce in Section 4 a wave propagation operator P t that behaves to some extent like a quantum evolution operator (preserving the quantum average ψ j, a ψ j of an observable a). Although not unitary, we choose it such that there exists a constant C I > depending only on the choice of the spectral interval I, so that ψ j, a ψ j C I ψ j, P t a P t dt ψ j. λ j I λ j I Studying the action of P t on eigenfunctions and proving this inequality is the object of Section 8. Note that we also choose P t to have characteristics of a wave propagator,

9 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 9 namely finite speed of propagation. his is because we will work locally, at a scale where the surfaces look like the hyperbolic plane around most of the points, and we need to prevent the propagation from going beyond this scale. () We then use the fact that the diagonal elements of a sufficiently regular operator in the orthonormal basis given by the eigenfunctions is bounded by the Hilbert-Schmidt norm of this operator: ψ j, λ j I P t a P t dt ψ j P t a P t dt We give a way to estimate the Hilbert-Schmidt norm on a surface using the kernel of the operator on the plane in Section 5. For this we divide between the points with a small radius of injectivity, that are controlled by the Benjamini-Schramm convergence, and points with a radius of injectivity greater than, the maximal distance of propagation. (3) We then study the kernel of the operator P t a P t. he kernel of this operator consists of averages of the observable a over intersection of balls of radius t. o study these averages we introduce in Section 6 a mean ergodic theorem of Nevo for averages over sets of increasing volume. Using the decay estimate given by this theorem and depending on the spectral gap, we would find that the norm of the operator P t a P t can be bounded by a quantity independent of t: P t a P t HS a. he idea is then to show that we have an orthogonality property of the type P t a P t dt = P t a P t HS dt such that P t a P t dt HS HS HS = P t a P t HS dt a. (3.) whenever is less than the radius of injectivity. Assuming for simplicity that the radius of injectivity goes to infinity, we obtain the type of estimate that we want. In reality, the orthogonality property is not obtained exactly in terms of the variable t, and appears in the form of a change of variable lemma. he computation is carried out in Section 7. (4) In (3.), the L norm of the observable a can be of the order of the area of the surface. We thus need to normalise our expression by dividing by this area. he link between the area and the number of eigenvalues in the interval I is provided by an analog of Weyl s law, that we prove in Section Wave propagation and reduction of the proof of heorem. Let X = Γ\H be a hyperbolic surface and H a fundamental domain. We define the operator P t u(z) := u(w) dµ(w), z H, cosh t B(z,t) for any function u : H R and t. he operator P t can be seen as a regularised version of the wave propagator, or a renormalised averaging operator over balls of radius t. It is self-adjoint in L (H) and it is an integral operator associated to the radial Γ-invariant kernel K t (z, z ) := k t (d(z, z )), z, z H H,

10 EIENNE LE MASSON AN UOMAS SAHLSEN where k t (ϱ) := cosh t {ϱ t}, ϱ. Fix now a test function a L (H). We still denote by a the operator of multiplication by the function a. he kernel of the conjugation P t a P t is given by [P t a P t ](z, z ) = a(w) dµ(w), z, z H. (4.) cosh t B(z,t) B(z,t) Note in particular that [P t a P t ](z, z ) = whenever d(z, z ) > t. We want to estimate the Hilbert-Schmidt norm of the time averages, P t a P t dt. HS Reducing the Quantum Ergodicity theorem to an estimate on the Hilbert-Schmidt norm of time averages appears in the proof of the large eigenvalue version of the theorem [4]. he way we estimate this average differs however from the large eigenvalue case. Another difference is that we use the L -normalised averaging operator on discs P t instead of the standard propagator of the wave equation. We prove the following bound: Proposition 4.. For any bounded a L (H), and for any >, we have P t ap t dt a HS ϱ(β) + e4 Vol({z X : InjRad l X (z) < }) a, min where ϱ(β) = + 4β with < β /4 the spectral gap of the Laplacian on X. We also need to know how the propagator P t acts on eigenfunctions of the Laplacian. his is given in the the following proposition: Proposition 4.. he Selberg transform of the kernel k t is given by the integral t h t (s) = cos(su) cosh t cosh u du. cosh t Moreover, for any fixed compact interval I (/4, ), there exist constants C I, I > such that for all s R with λ = /4 + s I, and for all I we have h t (s) dt C I. From these propositions we can prove heorem.: Proof of heorem.. Assume that a d Vol X =. Since P t is an operator of convolution with a radial kernel k t, we have by Proposition. that the action of the operator P t on an eigenfunction ψ j with eigenvalue λ j is given by P t ψ j = h t (s j )ψ j, where h t is the Selberg transform of the kernel k t and λ j = /4 + s j. We deduce from this and the symmetry of P t that ( ) ψ j, P t a P t dt ψ j = h t (s j ) dt ψ j, aψ j. j:λ j I j:λ j I We know by Proposition 4. that there is a constant C I depending only on the interval I such that inf j:λ j I h t (s j ) dt C I.

11 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE We thus have by Proposition 4. that ψ j, aψ j ( ψ j, C I λ j I I λ j I When we divide by the volume we obtain P t a P t dt HS ) P t a P t dt ψ j I a β + e4r l min Vol({z X : InjRad X (z) < R}) a ψ j, aψ j a I Vol(X) β + e4 Vol({z X : InjRad X (z) < }) a l min Vol(X). s j I hen Lemma 9. on the relationship between the number N(X, I) of eigenvalues λ j I and the volume Vol(X) yields the claim of heorem. o complete the proof of the main theorems, we need to prove Proposition 4. and Proposition 4.. In order to prove the first proposition, we start with two preliminary sections. We first give a way to estimate the Hilbert-Schmidt norm of the operators we are considering that depends on the radius of injectivity. hen we introduce a mean ergodic theorem in L for averaging operators over sets of increasing volume. 5. Hilbert-Schmidt norm and injectivity radius he Hilbert-Schmidt norm of a bounded operator A on X is given by the trace of A A A HS = r(a A), where A is the adjoint ot A. If K is the kernel of A, then A HS = K(z, z ) dµ(z)dµ(z ). he following lemma tells us how to compute the Hilbert-Schmidt norm of integral operators with Γ-invariant kernels on the surface X = Γ\H in terms of the injectivity radius, connecting it to the Benjamini-Schramm convergence. Lemma 5.. Let K : H H R be invariant under the diagonal action of Γ. he integral operator A on X defined by the kernel has the following Hilbert-Schmidt norm: A HS = K(z, γ w) dµ(z)dµ(w). γ Γ If K(z, w) = whenever d(z, w) R, we have A HS K(z, w) dµ(z)dµ(w) H + er l min Vol{z X : InjRad X (z) < R} sup (z,w) H where l min denotes the length of the shortest closed geodesic on Γ\H. K(z, w), (5.) Proof. We denote by (R) the point in the fundamental domain with radius of injectivity greater than R: (R) = {z : InjRad X (z) R},

12 EIENNE LE MASSON AN UOMAS SAHLSEN and by (R) c the complement of this set in. We then split the integral into two parts, and use that in the first part, the sum over Γ is reduced to one term. A HS = K(z, γ w) dµ(z)dµ(w) + K(z, γ w) dµ(z)dµ(w) γ Γ (R) γ Γ (R) c For the second sum, we use that it contains at most e R /l min terms. Indeed, we know that the volume of the ball of radius R is equal to cosh R and the minimal distance between two lattice points is l min, the number of lattice points is thus bounded by (cosh(r + l min ) )/(cosh l min ), which is bounded by e R /l min. We get using Cauchy-Schwarz inequality A HS H K(z, w) dµ(z)dµ(w) + er l min (R) c he second term on the right-hand side is bounded by e R l min K(z, γ w) dµ(z)dµ(w). γ Γ sup K(z, w) Vol(B H (z, R)) dµ(z) (z,w) H (R) c er l min Vol{z X : InjRad X (z) < R} sup (z,w) H K(z, w) 6. Mean ergodic theorem for averaging operators Our main tool from ergodic theory is given by an equidistribution result for averaging operators over general sets. Let us first introduce the problem in its general setting. Let (Y, ν) be a probability space, and G a group equipped with left-invariant Haar measure dg, and a measure-preserving action on Y. For a collection of measurable sets F t G we are interested in the behaviour of the averaging operators defined by π Y (F t )f(x) = F t F t f(g x) dg, f L (Y ), x Y, as t +. An extensive study of this question has been published by Gorodnik and Nevo [8] (see also the survey [9]). We will use in particular a theorem of Nevo, written in the following form in [9, heorem 4.] heorem 6. (Nevo [3, heorem ]). If G is a connected simple Lie group and the measure-preserving action on the probability space (Y, ν) has a spectral gap, then there exist C, θ > such that for any family F t G, t, of measurable sets of positive measure, we have π Y (F t )f f dµ C F t θ f L (Y,ν) Y L (Y,ν) for any f L (Y, ν), where we denote by F t the measure of the set F t. he constant C depends only on G and θ depends only on the spectral gap. For a family F t with volume increasing in t, heorem 6. gives a rate at which the averages π Y (F t )f converge to Y f dµ in L (Y ), and this rate only depends on the spectral gap and volume of F t. Let us insist on the fact that no further assumption on F t is needed. We will apply this theorem in the next section with G = PSL(, R), Y = Γ\ PSL(, R), and for suitable sets F t arising from the wave propagators P t, to bound the Hilbert-Schmidt norms of the time averages P tap t dt.

13 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 3 7. Bounding the Hilbert-Schmidt norm o prove Proposition 4. using Lemma 5. applied to the kernel (4.) of P t a P t, it is sufficient to bound the quantity [P t ap t ](z, z ) dt dµ(z ) dµ(z) H = cosh(t) a(w) dµ(w) dt dµ(z) dµ(z ). H B(z,t) B(z,t) In order to use the L estimate for ergodic averages given in heorem 6., we will rewrite the previous expression as the L norm of an averaging operator acting on the function a. For this purpose we introduce a change of variable lemma that allows us to go from the double integral over the surface to an integral over the unit tangent bundle. For R > write B R = {(z, z ) H d(z, z ) < R}. We define the mapping Φ : B R S (, R) by Φ(z, z ) = (m(z, z ), θ(z, z ), d(z, z )), where m(z, z ) is the middle of the geodesic segment [z, z ], θ(z, z ) is the direction of a unit vector in m(z, z ) tangent to this segment, and d(z, z ) is the geodesic distance between z and z. his gives us the following change of variable. Lemma 7.. For any f : H S [, ) C invariant under the action of Γ, i.e. such that γ Γ, f(γ (z, θ), r) = f(z, θ, r), we have B R f(m(z, z ), θ(z, z ), d(z, z )) dµ(z) dµ(z ) = R sinh(r) f(z, θ, r) dθ dµ(z) dr. S Figure. Lemma 7. allows us to go from the double integration over middle points to an integral on the unitary tangent bundle. Proof. We fix the variable z, and we use for z H the polar coordinates centred in z, z (θ, r), where r = r(z, z ) = d(z, z ) and θ = θ(z, z ) is the direction of the unit vector at m(z, z ) tangent to the geodesic (z, z ), see Figure. As we have dµ(z ) = sinh r drdθ, this change of variable gives us R f(m(z, z ), θ(z, z ), d(z, z )) dµ(z)dµ(z ) = f(ϕ r/ (z, θ), r) sinh(r) dµ(z) dθ dr, B R S with ϕ r/ the geodesic flow at time r/.

14 4 EIENNE LE MASSON AN UOMAS SAHLSEN Now we note that dµ(z) dθ is the Liouville measure, which is invariant under the action of the geodesic flow (z, θ) ϕ r/ (z, θ). We thus have R f(m(z, z ), θ(z, z ), d(z, z )) dµ(z)dµ(z ) = f(z, θ, r) sinh(r) dµ(z) dθ dr. B R We are now ready to prove the proposition. Proof of Proposition 4.. We have [P t ap t ](z, z ) dt dµ(z ) dµ(z) H = cosh(t) a(w) dµ(w) dt dµ(z) dµ(z ) H B(z,t) B(z,t) For every r (, ), we define the family of sets F t (r) PSL(, R) such that for any point (z, θ) H S, the set F t (r) (z, θ) is the lift to the unit tangent bundle of the intersection of two balls of radius t, where the two centres are given by the projection on the surface of ϕ r/ (z, θ) and ϕ r/ (z, θ). S Figure. he volume of the sets F t (r) used in the proof of Proposition 4. can be controlled by the volume of the balls B(z, t r/) and B(z, ϱ), where cosh ϱ = by the hyperbolic version of Pythagoras theorem. cosh t cosh(r/) he volume of both of these balls is O(e t r/ ). We write for any function f L (Γ\H) π(f t (r))f(z, θ) = F t (r) F t(r) f(g (z, θ)) dg,

15 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 5 and we have, using the change of variable Lemma 7., [P t ap t ](z, z ) dt dµ(z ) dµ(z) H = sinh r cosh(t) F t (r) π(f t (r))a(z, θ) dt dµ(z) dθ dr, S r/ ( ) sinh r cosh(t) F t (r) π(f t (r))a L ( S ) dt dr. r/ Where we used Minkowski s integral inequality to obtain the last line. We then apply heorem 6. with G = PSL(, R) and Y = Γ\ PSL(, R), identifying the latter to S. We obtain that there is a constant ϱ(β) > depending only on the spectral gap such that the previous quantity is bounded up to a uniform constant by ( sinh r cosh(t) F t (r) ϱ(β) a dt) dr. r/ We then remark that the volume of the intersection of balls of radius t with centres at a distance r from each other is O(e t r/ ), which gives us an estimate of F t (r). Hence we have [P t ap t ](z, z ) dt dµ(z ) dµ(z) H ( sinh r e r/ e ϱ(β)(t r/) a dt) dr r/ a ϱ(β) dr a ϱ(β). 8. Spectral action of the propagator In this section, we study the spectral action of the propagator P t and prove Proposition 4.. We first give the expression for the Selberg transform we will need: Lemma 8.. he Selberg transform of the kernel k t (ϱ) = {ϱ t} cosh t is given by the integral h t (s) = t cos(su) Proof. Applying the Selberg transform to the kernel k t, we have with g(u) = h t (s) = + cosh u du. (8.) cosh t e isu g(u) du t sinh ϱ dϱ = cosh t cosh u ( t,t) (u). cosh t u cosh ϱ cosh u cosh t Using the fact that cosh( u) = cosh(u) we obtain (8.).

16 6 EIENNE LE MASSON AN UOMAS SAHLSEN Now to obtain Proposition 4., we will need that for a fixed compact interval I (/4, ), we want to find constants C I, I > such that for all λ I with λ = /4 + s and I we have h t (s) dt C I. (8.) o simplify the notation we will denote also by I the interval of the parameter s and write s I. Figure 3. Illustration of the Selberg transform h t (s) of the kernel k t (ϱ) with k =. he constant c(i) is constructed from the integral up to time t (in Lemma 8.3) using the compactness of I. he remaining part from t to t may subtract mass from this, but we can control this uniformly over k due to the uniform Lipschitz continuity of t h t (s) in Lemma 8.. he idea of the proof of (8.) is to show first that the function t h t (s) satisfies a Lipschitz regularity condition (see Lemma 8. below), and then use the fact that uniformly over some increasing sequence (t k ) k N such that t k, there is a given positive lower bound for h t (s), depending only on s (see Lemma 8.3 below). See Figure 3 for an illustration of the proof. Let us first prove the Lipschitz regularity of t h t (s): Lemma 8.. he function t h t (s) is uniformly Lipschitz on (, ) with constant independent of s. Proof. Let us assume that t >. We are going to prove that the derivative of t h t (s) is bounded uniformly in s. he derivative is given by t t h t (s) = sinh t cosh u cosh(t) 3/ cos(su) du cosh t cosh u t cosh u du. cosh t cosh t cosh u Fix < δ < t. We divide the integral into two parts: on the interval (δ, t) we have the inequality cosh u coth δ sinh u = cosh δ sinh u, sinh δ

17 and QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 7 coth δ t sinh u = coth δ cosh t cosh δ cosh t cosh t cosh u cosh t δ coth δ. On the interval (, δ) we have the bound, using that t > δ cosh u du cosh t cosh t cosh u δ cosh δ cosh t cosh t cosh δ ( ) / cosh δ cosh δ. We thus obtain a bound for t h t (s) that is uniform in s and t for any choice of δ, and we deduce that the function h t (s) is Lipschitz in t (, + ) uniformly in s by the mean value theorem. Now let us construct the desired sequence (t k ) k N, for which we have uniform bounds for h t (s) in t and over s I. For s > and k N write t k := πk s. Lemma 8.3. Given a bounded interval I, there exists a constant c(i) > and k (I) N such that for any k k (I) and any s I we have h tk (s) < c(i). Proof. First of all, for any j =,..., k we have tj cos(su) cosh u tj du = sin(su) t j cosh t k s cosh t k t j since the map u (sinh u)/ cosh u cosh t k [, t k ). hus it is enough to find c(i) > such that tk t k cos(su) sinh u du < cosh u cosh t k is continuous and monotonically increasing on cosh u cosh t k du < c(i). (8.3) he idea is to obtain a bound that is independent of k and s I. Recall that t k = πk s. Let us do a change of variable v := u t k which maps [t k, t k ] onto [, π/s]. hen as cos(s(v + t k )) = cos(sv + π(k )) = cos(sv) for all v [, π/s], we have tk t k cos(su) For v [, π/s] denote hen as t k = t k π s, we have f k (v) = cosh u cosh t k du = f k (v) := π s ev+t k + e (v+t k ) e t k + e t k = cos(sv) cosh(v + t k ) cosh t k. cosh(v + t k ) cosh t k dv. e v π s + e (v+t k 4π + e t k Hence we have that the sequence of functions (f k ) k N converges (as k ) pointwise to the limit function f(v) := e v π s. s ).

18 8 EIENNE LE MASSON AN UOMAS SAHLSEN hanks to the compactness of I, this convergence is uniform over s I as we will see now. Write I = [a, b], fix s I and v [, π/s]. As s a, t k = πk/s and v we obtain hus f k (v) f(v) and + e (v+t k 4π s ) + e t k f k (v) + e 4( k)π a. e v π s ( + e 4( k)π a ). By the mean value theorem, we can find ξ [ e v π s f(v) e v π s ( + e 4( k)π a ) = Since e v π s hus we have proved we have by the choice of ξ that e v π s ξ e 4( k)π a for all s I and v. herefore π s cos(sv) sup s I sup s I π a π a e ( k)π a. Now for any s > write e 4( k)π a e 4( k)π a χ [,π/s] (v) f k (v) f(v) e ( k)π a cosh(v + t k ) cosh t k χ [,π/s] (v) f k (v) f(v) dv c(s) := π s ( + e 4( k)π a ), e v π s ] satisfying π ev s ξ e 4( k)π a. = e ( k)π a. π s dv cos(sv) e v π s dv cos(sv) e v π s dv. hen c(s) >. Indeed, a change of variable x := sv and integration by parts gives π s cos(sv) e v π s dv = π cos(x) e x π s dx = π s s sin(x) e x π s e x π s so the positivity follows as x e x π s / e x π s is strictly increasing on [, π). Hence, as I is compact and s c(s) is continuous, we have that the infimum c(i) := inf c(s) >. s I Since π a e ( k)π a as k and this sequence only depends on I, the above computation yields k (I) N such that for all k k (I) and all s I we have that π s cos(sv) cosh(v + t k ) cosh t k dv + c(s) c(i). Since c(i) c(s) for all s I, we have for any s I that π s cos(sv) cosh(v + t k ) dv < c(s) + c(i) c(i). cosh t k dx

19 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE 9 Multiplying both sides by yields the claim by the definition of h tk (s). Now we are ready to prove Proposition 4.. Proof of Proposition 4.. Let > and s I. By Lemma 8.3, we know that for any k k(i), h tk (s) c(i). We thus have a sequence of times (t k ) k k (I) on which h tk (s) is bounded uniformly from below. Using the uniform Lipschitz continuity of t h t (s) on (, + ) we can find intervals around each of these points on which the function is also uniformly bounded from below, so that we have a lower bound on a set of positive measure. More precisely, we can find an interval J around and k k (I) such that for any k k, h t (s) c(i) on the translated interval J + t k. enote N A(s, ) = J + t k, k=k where N is chosen such that Note that and we have in particular and using that N s 4π, we obtain πn/s < π(n + )/s. A(s, ) (N k ) J >, ( ) N h t (s) k dt J c(i) h t (s) dt ( s 4π k ) J c(i). hus for I large enough there exists a constant C I > such that for any I which completes the proof. h t (s) dt C I, 9. Asymptotic number of eigenvalues In this section, we prove an analog of Weyl s law for eigenvalues in a fixed bounded interval and sequences of surfaces converging in the sense of Benjamini-Schramm. his allows us to normalise by the number of eigenvalues in heorem.. his estimate of the number of eigenvalues is a consequence of the results of [6, 7] in the case of sequences of coverings, and more generally [] for Benjamini-Schramm convergence. We provide a proof for completeness in the special case we are considering. Lemma 9.. For any sequence of compact hyperbolic surfaces X n = Γ n \H converging to the plane in the sense of Benjamini and Schramm, and any interval I (, + ), we have the following asymptotics: N(X n, I) I Vol(X n ), where N(X n, I) is the number of eigenvalues in the interval I.

20 EIENNE LE MASSON AN UOMAS SAHLSEN o establish this result, we use Selberg s pre-trace formula. For a proof of this theorem and a more general treatment of Selberg s trace formula, see for example [9, 4, 8]. Given h : C C, the invariant kernel k(z, z ) = k h (z, z ) associated to h is defined by k(z, z ) = k h (z, z ) := π d(z,z ) (ĥ) (ϱ) cosh ϱ cosh d(z, z ) dϱ. his kernel is simply obtained from the inverse Selberg transform of h. heorem 9. (Selberg s pre-trace formula). Let X be a compact hyperbolic surface. enote by (ψ j ) an orthonormal basis of eigenfunctions of the Laplacian in L (X) and by (λ j ) the corresponding non decreasing sequence of eigenvalues. Suppose h : C C satisfies: () h is analytic on the strip Im z σ for some σ > /; () h is even, that is, h( z) = h(z); (3) h(z) ( + Re z ) δ for some fixed δ >, uniformly for all z in the strip Im z σ. hen the following formula converges absolutely: j= h(s j ) ψ j (z) = 4π R h(ϱ) tanh(πϱ)ϱ dϱ + γ Γ(X) {id} uniformly in z H, where k(z, z ) is the invariant kernel associated to h. k(z, γ z) We want to replace h in the previous formula with the indicator function of the bounded interval we are considering. However, as we need h to be analytic, we will use an approximation of the indicator by analytic functions. For this purpose we can use Beurling-Selberg extremal functions. We have the following result (see for example [38] for explicit formulas) heorem 9.3 (Beurling-Selberg extremal functions). Given an interval I R and m N there exists entire functions h m, H m : C C such that () the functions h m, H m minorise and majorise the indicator function of I: () the functions h m, H m satisfy h m (x) I (x) H m (x), x R; h m χ I L (R) = H m χ I L (R) = m ; (3) the Fourier transforms ĥm, Ĥm are compactly supported in [ m, m]. Now that we have the main tools, we can prove Lemma 9.. Proof of Lemma 9.. Let m, n N and write X n (m) := {z X n : InjRad Xn (z) m}. We apply Selberg s pre-trace formula (heorem 9.) for the surface X n and the functions h m and H m from heorem 9.3 and obtain h m (s j ) ψ j (z) = h m (ϱ) tanh(πϱ)ϱ dϱ + k m (z, γ z) 4π and j= j= H m (s j ) ψ j (z) = 4π R R H m (ϱ) tanh(πϱ)ϱ dϱ + Where the associated invariant kernels k m and K m are given by k m (z, z ) = π d(z,z ) γ Γ n {id} γ Γ n {id} (ĥm) (ϱ) cosh ϱ cosh d(z, z ) dϱ K m (z, γ z).

21 QUANUM ERGOICIY AN BENJAMINI-SCHRAMM CONVERGENCE and K m (z, z ) = π d(z,z ) (Ĥm) (ϱ) cosh ϱ cosh d(z, z ) dϱ. Note that, as the Fourier transforms ĥm and Ĥm are compactly supported in [ m, m], so are their derivatives, and as a consequence both k m (z, z ) and K m (z, z ) vanish whenever d(z, z ) > m. When we integrate over z we thus obtain j= h m (s j ) = Vol(X n) 4π R h m (ϱ) tanh(πϱ)ϱ dϱ + X n(m) γ Γ n {id} k m (z, γ z) dµ(z). We have X n(m) γ Γ n {id} which gives ( h m (s j ) Vol(X n) 4π j= k m (z, γ z) dµ(z) em Vol(X n (m)) sup k m (z, z ), l min z,z X n R h m (ϱ) tanh(πϱ)ϱ dϱ em Vol(X n (m)) ) sup k m (z, z ) l min Vol(X n ) z,z and for m large enough we can remove the absolute value from the left-hand side of the inequality, as the series is non-negative. his is because of the property () for Beurling-Selberg extremal functions: it yields that h m χ I in L as m so there is a subsequence h mk χ I pointwise, so dominated convergence theorem and Benjamini- Schramm convergence yields the claim. Let us define c I = 4π I tanh(πϱ)ϱ dϱ. For any ε >, there exist m such that h m (ϱ) tanh(πϱ)ϱ dϱ c I < ε 4π. R Using the Benjamini-Schramm convergence of X n to H, we can then find n large enough such that e m Vol(X n (m)) l min Vol(X n ) sup k m (z, z ) < ε z,z X n. Note that sup z,z k m (z, z ) is bounded from above by a number independent of n: sup k m (z, z ) (ĥm) (ϱ) dϱ < z,z X n 4π cosh ϱ so this choice of n is possible. herefore, for n large enough we have N(I, X n ) Vol(X n )(c I + ε). Similarly, using the function H m we obtain that for any ε >, we have N(I, X n ) Vol(X n )(c I + ε) for n large enough. We deduce from the two previous inequalities that N(I, X n ) c I Vol(X n ) when n +.

22 EIENNE LE MASSON AN UOMAS SAHLSEN Acknowledgements We thank Yves Colin de Verdière for suggesting the problem. We are grateful to Nalini Anantharaman, Elon Lindenstrauss and Nicolas de Saxcé for useful discussions and comments. Moreover, we thank the enthusiastic community in the School of Mathematics at the University of Bristol, where the majority of the research was carried out, in particular, Alex Gorodnik, homas Jordan, Jens Marklof, Abhishek Saha, Roman Schubert and Corinna Ulcigrai for several helpful conversations about this work. Finally we thank Zeev Rudnick and Steve Zelditch for useful comments on the first version of the manuscript. References [] Miklos Abert, Nicolas Bergeron, Ian Biringer, sachik Gelander, Nikolay Nikolov, Jean Raimbault, and Iddo Samet. On the growth of L -invariants for sequences of lattices in lie groups. Preprint, 5 (latest), arxiv: [] Miklos Abert, Nicolas Bergeron, Ian Biringer, sachik Gelander, Nikolay Nikolov, Jean Raimbault, and Iddo Samet. On the growth of Betti numbers of locally symmetric spaces. C. R. Math. Acad. Sci. Paris, 349(5-6):83 835,. [3] Nalini Anantharaman. Quantum ergodicity on large graphs. Preprint, 5, arxiv: [4] Nalini Anantharaman. Entropy and the localization of eigenfunctions. Ann. of Math. (), 68(): , 8. [5] Nalini Anantharaman and Etienne Le Masson. Quantum ergodicity on large regular graphs. uke Math. J., 64(4):73 765, 5. [6] Nalini Anantharaman and Stéphane Nonnenmacher. Half-delocalization of eigenfunctions for the Laplacian on an Anosov manifold. Ann. Inst. Fourier (Grenoble), 57(7):465 53, 7. Festival Yves Colin de Verdière. [7] Itai Benjamini and Oded Schramm. Recurrence of distributional limits of finite planar graphs. Electron. J. Probab., 6:no. 3, 3 pp. (electronic),. [8] Nicolas Bergeron. he spectrum of hyperbolic surfaces. ranslated from the French by Farrell Brumley. Les Ulis: EP Sciences; Cham: Springer, 6. [9] Lewis Bowen. Cheeger constants and L -Betti numbers. uke Math. J., 64(3):569 65, 5. [] Robert Brooks and Eran Makover. Random construction of Riemann surfaces. J. ifferential Geom., 68(): 57, 4. [] Shimon Brooks, Etienne Le Masson, and Elon Lindenstrauss. Quantum ergodicity and averaging operators on the sphere. IMRN, to appear, 5, arxiv: [] Shimon Brooks and Elon Lindenstrauss. Non-localization of eigenfunctions on large regular graphs. Israel J. Math., 93(): 4, 3. [3] Shimon Brooks and Elon Lindenstrauss. Joint quasimodes, positive entropy, and quantum unique ergodicity. Invent. Math., 98():9 59, 4. [4] Peter Buser. Geometry and spectra of compact Riemann surfaces. Modern Birkhäuser Classics. Birkhäuser Boston, Inc., Boston, MA,. Reprint of the 99 edition. [5] Yves Colin de Verdière. Ergodicité et fonctions propres du laplacien. Comm. Math. Phys., (3):497 5, 985. [6] avid L. de George and Nolan R. Wallach. Limit formulas for multiplicities in L (Γ\G). Ann. of Math. (), 7():33 5, 978. [7] avid L. egeorge and Nolan R. Wallach. Limit formulas for multiplicities in L (Γ\G). II. he tempered spectrum. Ann. of Math. (), 9(3): , 979. [8] Alexander Gorodnik and Amos Nevo. he ergodic theory of lattice subgroups, volume 7 of Annals of Mathematics Studies. Princeton University Press, Princeton, NJ,. [9] Alexander Gorodnik and Amos Nevo. Quantitative ergodic theorems and their number-theoretic applications. Bull. Amer. Math. Soc. (N.S.), 5():65 3, 5. [] Xiaolong Han. Small scale quantum ergodicity in negatively curved manifolds. Nonlinearity, 8(9): , 5. [] Andrew Hassell. Ergodic billiards that are not quantum unique ergodic. Ann. of Math. (), 7():65 69,. With an appendix by the author and Luc Hillairet. [] Hamid Hezari and Gabriel Rivière. L p norms, nodal sets, and quantum ergodicity. Adv. Math., 9: , 6. [3] Roman Holowinsky and Kannan Soundararajan. Mass equidistribution for Hecke eigenforms. Ann. of Math. (), 7():57 58,.