Arxiv:1908.03805V1 [Math.SP]

Home , Svetlana Jitomirskaya

arXiv:1908.03805v1 [math.SP] 10 Aug 2019 perturbative pc ieso.Tefitmlidmninllclzto a bandb for obtained was Chulaevsky-Dinaburg by localization multi-dimensional n methods fist a The technique becomes dimension. dynamical) bad, space therefore go (and sometimes one-dimensional do portantly, actually things o where c number more parameters (the significantly torus becomes 14 situation undelying [2, the the the dimension, references space of (e.g. and the universality dimension cially, 21] related the 19, and either ([4, transitions increases sides arithmetic coupling now sharp high are and and There small arguments. pertur both analyticity earlier on co the reducibilit by replaced theory powerful enhanced rich that where particularly remarkably 18] potentials, and non-perturbat 17, one-frequency well-developed first 8, dimensional most of 6, The development [4, the techniques. functions since Green’s years, 20 the last the in ntemlifeunycase, multi-frequency the in d a oapysm ftennprubtv ehd otetodime two the for to coupling methods high non-perturbative at the localization of some apply to way (1.2) where 5,weeh eeoe e oeflshm htalwdt circum to allowed that scheme powerful to [7] new restricted a that developed difficulties he where [5], (1.1) aeo general of case 1] oee,prubtv ehiushv o enmd owor to made been not have techniques perturbative However, [11]. h hoyo usproi ptrr ihaayi oetashss has potentials analytic with opetarors quasiperiodic of theory The Let e od n phrases. and words Key nti ae,w osdrtefloigoperators following the consider we paper, this In Let | v S n prtr on operators Abstract. | eara nltcfnto on function analytic real a be eaTelt oeao)mti on matrix (operator) Toeplitz a be =max := NESNLCLZTO O MULTI-FREQUENCY FOR LOCALIZATION ANDERSON 1 ≤ i k ≤ eut eoe eeal,fle[] otrwn wysalmeasure small away throwing so [4], false generally, become, results VTAAJTMRKY,WNA I,ADYNEGSHI YUNFENG AND LIU, WENCAI JITOMIRSKAYA, SVETLANA , d d | eetbihAdro oaiainfrgnrlanalytic general for localization Anderson establish We n Z i at nee infiatymr eea situation). general more significantly even an fact, (in i d | US-EIDCOEAOSON OPERATORS QUASI-PERIODIC for for nesnlclzto,ln-ag us-eidcoperat quasi-periodic long-range localization, Anderson arbitrary H n 1. ( ( = x > k nrdcinadmi results main and Introduction = ) n k k | 1 S k , = ,ee for even 1, S = , ( n d n 2 + . d , d , k n .Ti a xeddb orant arbitrary to Bourgain by extended was This 2. = · · · λv T .I hspprw xedBugi’ eutt the to result Bourgain’s extend we paper this In 2. = feunyoeaoson operators -frequency ′ ) b ≤ | ( where , , x n + d ℓ e ) 2 − d 1 nω ( ∈ ρ Z | .Bugi-odti-clgdvlpda developed Bourgain-Goldstein-Schlag 1. = n d Z ) b − δ satisfying, ) d nn = n . ′ | ′ i , P , =1 d > ρ n b i , 0, n ( b ′ i ∈ ∈ ℓ Z r o plcbei higher in applicable not are s 2 k N /yaia ehiusare techniques y/dynamical d ( feunyquasi-periodic -frequency Z , v ehd o oto of control for methods ive o 1 for d soa ae[] obtaining [7], case nsional Z aino eigenfunctions of bation for ) o-etraieresults non-perturbative csiy vnmr im- more Even ecessity. mlctd First, omplicated. en,goa hoy[1], theory global rein), crstecs fone- of case the ncerns rqece)o,espe- or, frequencies) f d 1].Hwvr fone if However, –16]). opoelocalization prove to k e rmtcadvances dramatic een r,sm-leri sets. semi-algebraic ors, etraie(KAM) perturbative y ≤ ettearithmetic the vent k i n arbitrary and 1 = ≤ d ). esof sets k non- = d 2 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI where x = (x , , x , , x , , x ) Tb, 11 ··· b11 ··· 1d ··· bdd ∈ nω = (n ω , , n ω , , n ω , , n ω ). 1 11 ··· 1 b11 ··· d 1d ··· d bdd Example 0. Taking bi = 1, i = 1, ..., d and the nearest neighbor Laplacian S we obtain operators considered in [5]. 3 Example 1. d = 2, b1 = 2, b2 = 1. v is a function on T . For x = (x1, x2, x3) and ω = (ω1, ω2, ω3), the operator becomes

(1.3) H = S + λv(x1 + n1ω1, x2 + n1ω2, x3 + n2ω3)δnn′ , where n = (n1, n2). k Example 2. b = kd, bi = k, i = 1, ..., d, f is a function on T , and v(x , , x , , x , , x )= f(x + + x , , x + + x ). 11 ··· k1 ··· 1d ··· kd 11 ··· 1d ··· k1 ··· kd Then the operator becomes

(1.4) H(x)= S + λf(x + nA)δnn′ , where x Tk, n Zd, and A is a d by k matrix of frequencies. This is the most general form of a d-dimensional∈ ∈ quasiperiodic operator with a k-dimensional phase space. The Aubry dual family has the form ˜ (1.5) H(x)= F + λs(x + Am)δnn′ , d k where x T , n Z , and F , S are Toeplitz operators with (n, n′) terms given by the n n′ Fourier coefficients∈ ∈ of, correspondingly f, s. The standard Laplacian is therefore dual to− the potential given by the sum of cosines, and the dual of a general analytic potential is a Toeplitz matrix as above. Remark 1. When considering families (1.4) with A restricted to a linear submanifold of d by k matrices of frequencies, one needs to take bi equal to the number of free variables in the ith row of A and adjust v accordingly. As such, the family considered in [5] can of course also be recast in this language: it corresponds to A restricted to diag(ω1, ..., ωd). We call x Tb the phase, ω Tb the frequency and λ 0 the coupling. Let ∈ ∈ ≥ xj := (x , , x ) Tbj (1 j d). 1j ··· bj j ∈ ≤ ≤ We assume v satisfies the following non-degeneracy condition: for any 1 j 1 j+1 d b b (x , , x − , x , x ) T − j , ··· ··· ∈ the function b 1 j 1 j+1 d T j θ v(x , , x − , θ, x , x ) ∋ 7→ ··· ··· is nonconstant. Denote by mes the Lebesgue measure. We say operator H satisfies Anderson localization if it has only pure point spectrum with exponentially decaying eigenfunctions. Theorem 1.1. Let H(x) be given by (1.2) with v satisfying the non-degeneracy condition. Then for any δ > 0, there is a λ0 = λ0(δ, v, ρ, b, d) > 0 such that the following statement holds: b b b for any λ λ0 and any x T , there exists Ω=Ω(x, λv, δ, ρ, b, d) T with mes(T Ω) δ such that for≥ ω Ω, H(x)∈satisfies Anderson localization. ⊂ \ ≤ ∈ 3

Remarks (1) In particular, this holds for all operators (1.4) with arbitrary k, d and any non-constant analytic function f,1 which is an important building block in the proof of absolutely continuous spectrum for operators (1.5) [9] and was the key initial motivation for our work. (2) We note that our phase space dimension b satisﬁes b d since b = d b , b 1. ≥ i=1 i i ≥ This is essential for our arguments. As shown in Example 2, generalP quasiperiodic operators always have b d. However, operators (1.2) with b < d, for example ≥ Vn1,n2 = v(x + n1ω, x + n2ω), also appear naturally, e.g. in the study of interacting particles, and our proof does not apply in this setting. A localization result for a model with b = 1, d = 2 was recently obtained by Bourgain-Kachkovskiy [10]. (3) Previous multidimensional/multifrequency localization results [5, 7] were not only restricted to k = d , but also done only for the nearest neighbor Laplacian, i.e.

S(n, n′)= δ n n′ ,1. The extension to general S as in (1.1) is motivated by the Aubry duality purposes| − | in [9]. Localization for long-range operators (general S) was previ- ously obtained for k = 1 in [11] and, nonperturbatively, for k = d = 1 in [4, 8]. The main scheme of our proof is definitely adapted from Bourgain [5]. However, while our result is significantly more general and more technically complex, our argument can also be viewed as both a clarification and at the same time streamlining of [5]. Indeed, our proof, while including more detail and hopefully increasing the readability, is only shorter than the corresponding part of Bourgain’s. This is due to several important technical improvements that we add to Bourgain’s scheme. One important highlight is that, in the process of deterministic multi-scale analysis proceeding from scale N1 to N2, a chain of scales between N1 and N2 has always been used in the past work, [5, 7]. Here, instead of gluing “good” Green’s function at multiple scales between N1 and N2 to establish the “goodness” of Green’s function at scale N2, we find a way to directly use the “good” Green’s function at scale N1 + subexponential bound of the norm to prove the “goodness” of Green’s function at scale N2. Another issue we want to highlight is that the k = d = 2 analysis of [7] required dealing with many different types of elementary regions, something that would be prohibitively difficult to carry out in higher dimensions. In dealing with higher dimensions in [5] Bourgain significantly reduces the allowed elementary regions. This comes at the price of some complications in dealing with the lattice points at the boundary of the elementary regions, which Bourgain claims can be carried out, but provides no detail. We use the same (slightly corrected) type of restriction on the elementary regions but believe this issue is not entirely trivial and tackling it requires a certain modification of the procedure, which we provide in full detail. Non-perturbative proofs of localization for d > 1 are in a sense a version of deterministic multi-scale analysis. The latter is a powerful method originally developed for random operators by Fröhlich and Spencer [12], that crucially relied on independence and Wegner’s Lemma that is effectively dependent on rank-one perturbations. For the deterministic version, difficulties with lack of independence/rank one perturbations are circumvented by the semi-algebraic sets considerations and subharmonicity arguments [4]. The non-perturbative proofs consist of two parts. First, one needs to obtain measure and complexity estimates for phases/frequencies with exponential off-diagonal decay and subexponential upper bounds for the matrix elements of the Green’s function for box-restricted operators for a given energy. From this, localization follows through elimination of energy via an argument involving complexity bounds on semi-algebraic

1As in Remark 1, the non-degenracy condition on v leads to additional non-degeneracy conditions on f if the number of free variables in a certain row of the submanifold is bounded by 1. In particular, for A restricted to diag(ω1,..., ωd), as in [5], the required non-degeneracy condition is exactly as in [5]. 4 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI sets. The second part is by now rather standard and follows the reasonably short argument in [5] essentially verbatim. In fact, it is the first part that presents the main difficulty associated with higher dimensions. Thus we focus only on the first, single energy, part here. This is also where the key difficulty in extending [7] and the key difference between [5] and [7] lies. One needs to guarantee a sublinear upper bound on the number of times the ergodic trajectory hits certain forbidden regions of given measure/algebraic complexity, without further detail on the structure of those fordbidden regions. A key argument in [7] is a Lemma that does guarantee it for k = b = 2 under an explicit arithmetic condition on the frequencies. Roughly, it means that too many points on the trajectory of rotation close to an algebraic curve of a bounded degree would force it to oscillate more than the degree allows. However, this statement is not extendable to d 3. In [5] Bourgain instead developed a way to restrict to suitable frequencies already for the first≥ step, which turned out to be a very robust approach that we also develop here. Besides the elimination of energy argument, we do not include detailed proofs of two further statements very similar to those in [5], and with proofs presented there in a very clear way. The proofs that are similar to Bourgain’s that we do present either have certain novelty or contain important technical clarifications. In Section 2 we introduce the main concepts and also list the above mentioned results for which we do not present detailed proofs. One such concept is “property P at scale N”- essentially, the single energy statement one wants to establish for all large scales, that allows to streamline certain formulations. Section 3 is devoted to the main multi-scale argument: property P at scales N, N c implies property P at an interval of subexponentially large scales, Theorem 3.6. In section 4 we take care of the initial scale and give a very short argument to obtain the final single energy estimate, Theorem 4.1, from Theorem 3.6. In the appendix we prove a several variables matrix-valued Cartan estimate (Lemma 3.4 used in the proof of Theorem 3.6), that follows Bourgain’s one-variable argument in [4] but also uses high- dimensional Cartan sets estimates of [13].

2. Preparations

2.1. Notation. For any x Rd1 and X Rd1+d2 , denote the x-section of X by ∈ ⊂ X(x) := y Rd2 : (x, y) X . { ∈ ∈ } ˜ Tb 1 j 1 j+1 d Tb bj Let b = maxi bi. For any x and 1 j d, let xj¬ = (x , , x − , x , x ) − . ∈ ≤ ≤ l ··· ··· ∈ For x = (x1, x2, , xl), y = (y1, y2, , yl) R , let x y = maxi xi yi . For Λ ,Λ Z···d, we introduce ··· ∈ | − | | − | 1 ⊂ d diam(Λ) = sup n n′ , dist(m, Λ) = inf m n (m R ), n,n Λ | − | n Λ | − | ∈ ′∈ ∈ and dist(Λ1, Λ) = inf dist(n, Λ). n Λ ∈ 1 We also use as ℓ2 norm of the matrix. For convenience, in the following, we study 1 k·k 1 operator λ− H(x). We always assume λ > 1. Since the spectra of λ− H(x) are bounded by C(S, v), we can further assume E is bounded.

d 2.2. Green’s functions and elementary regions. For Λ Z , let RΛ be the restriction operator, i.e., (R ξ)(n) = ξ(n) for n Λ, and (R ξ)(n) =⊂ 0 for n / Λ. Denote by H = Λ ∈ Λ ∈ Λ RΛHRΛ and the Green’s functions 1 1 G (E; x) = (R (λ− H E + i0)R )− . Λ Λ − Λ We will also write GΛ when there is no ambiguity. Clearly,

(2.1) Gn+Λ(x)= GΛ(x + nω). 5

We denote by QN an elementary region of size N centered at 0, which is one of the following regions, Q = [ N, N]d N − or Q = [ N, N]d n Zd : n ς 0,1 i d , N − \{ ∈ i i ≤ ≤ } d where for i = 1, 2, , d, ςi <, >, and at least two ςi are not . 0 ··· ∈{ ∅} ∅ Denote by N the set of all elementary regions of size N centered at 0. Let N be the set of all translatesE of elementary regions, namely, E

N := n + QN n Zd,Q 0 . E { } ∈ N ∈EN 2.3. Semi-algebraic sets. Definition 2.1 (Chapter 9, [4]). A set Rn is called a semi-algebraic set if it is a finite union of sets defined by a finite numberS ⊂ of polynomial equalities and inequalities. More precisely, let P1, , Ps R[x1, , xn] be a family of real polynomials whose degrees are bounded by d.{ A (closed)··· } semi-algebraic ⊂ ··· set is given by an expression S (2.2) = x Rn : P (x)ς 0 , S { ∈ ℓ jℓ } [j ℓ\ ∈Lj where j 1, , s and ςjℓ , ,= . Then we say that has degree at most sd. In fact, the degreeL ⊂{ of ···which} is denoted∈ {≥ by deg(≤ }), means the smallestS sd over all representations as in (2.2). S S In [5], Bourgain proved a result for eliminating several variables. Lemma 2.2 (Lemma 1.18, [5]). Let [0, 1]d+r be a semi-algebraic set of degree B and such that S ⊂ mes( (y)) < η for y [0, 1]r. S ∀ ∈ Then the set

d2r (x , , x r ) [0, 1] : (x ) =  1 ··· 2 ∈ S i 6 ∅  1 \i 2r  ≤ ≤ is semi-algebraic of degree at most BC and measure at most 

r r(r 1)/2 BC ηd− 2− − , where C = C(d, r) > 0. Another important fact is the following decomposition Lemma for semi-algebraic sets in the product spaces.

Lemma 2.3 ([4, 5]). Let [0, 1]d=d1+d2 be a semi-algebraic set of degree deg( )= B and mes ( ) η, where S ⊂ S d S ≤ 1 (2.3) log B log , ≪ η with 1 η d ǫ. ≤ Then there is a decomposition of as S = S S1 ∪ S2 6 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI such that the projection of on [0, 1]d1 has small measure S1 C(d) mes (Proj d ) B ǫ, d1 [0,1] 1 S1 ≤ and has the transversality property S2 C(d) 1 1 mes ( ) B ǫ− η d , d2 L ∩ S2 ≤ where is any d -dimensional hyperplane in [0, 1]d s.t., L 2 max Proj (ej ) <ǫ, 1 j d1 | L | ≤ ≤ where we denote by e , , e the coordinate vectors in Rd1 . 1 ··· d1 We then have r Lemma 2.4. Suppose that ωi Rli (i = 1, 2, 3, , r) and l = l . Let [0, 1]lJ be a ∈ ··· i=1 i S ⊂ semi-algebraic set of degree B and such that P mes( ) < η. S For ω = (ω1, , ωr) [0, 1]l and n = (n , n , , n ) Zr, deﬁne ··· ∈ 1 2 ··· r ∈ nω = (n ω1, n ω2, , n ωr). 1 2 ··· r 1 J 1 r Let , , − Z be ﬁnite sets with the following property N ··· N ⊂ C min ns > (B max ms ) , 1 s r 1 s r ≤ ≤ | | ≤ ≤ | | i i 1 where n , m − (2 i J 1), where C = C(J, l). Assume also ∈ N ∈ N ≤ ≤ − 1 (2.4) max n C < . n J 1 | | η ∈N − Then l (i) i (1) (J 1) lJ C mes( ω [0, 1] : n s.t., (ω, n ω, , n − ω) mod Z ) B δ, { ∈ ∃ ∈ N ··· ∈ S} ≤ where 1 δ− = min min ns . n 1 1 s r | | ∈N ≤ ≤ Proof. The proof follows from Lemmas 2.2 and 2.3 just as the proof of Lemma 1.20 in [5] follows from the corresponding Lemma 1.18 and property (1.5) of semi-algebraic sets in [5].

Definition 2.5. We say (E, x) is (¯ρ, N) good, if for any Q 0 , N ∈EN (2.5) G (E; x) e√N , k QN k≤ ρ¯ n n′ N (2.6) G (E; x)(n, n′) e− | − | for n n′ . | QN |≤ | − |≥ 10 Definition 2.6. We say Green’s function satisfies property P with parameters (γ,¯ρ) at size N if there is a semi-algebraic set Ω = Ω (λv, ρ, b, d) Tb with deg(Ω ) N 4d such N N ⊂ N ≤ that the following statement is true: for any ω ΩN and E R, there exists a set XN = X (λv, ρ, b, d, ω, E) Tb such that ∈ ∈ N ⊂ N γ (2.7) sup mes(XN (xj¬)) e− , b b 1 j d,x T − j ≤ ≤ ≤ j¬∈ and for any x not in XN , (E, x) is (¯ρ, N) good. 7

Theorem 2.7. There exist small positive constants c3 < c4 < 1, where c3 and c4 depend on b, d such that the following statements are true. Let c = c3 and c = c2/2. Fix a large 1 4˜b 2 1 c 2/c1 2 N2 number N1. Let N2 = N1 and N3 = e . Suppose the Green’s functions satisfy property

P at size N1 with parameters (c1,¯ρ), and corresponding semi-algebraic sets ΩN1 . Then there 4d c3 exists a semi-algebraic set Ω3 ΩN1 with deg(Ω3) N3 and mes((ΩN1 Ω3) N3− such ⊂ ≤ c \ c ≤ that, if ω Ω , then for any E R and x Tb, there exists N 3

3. Resolvent identities and Cartan’s Lemma d 1 Let Λ1,Λ2 Z and Λ1 Λ2 = . LetΛ=Λ1 Λ2. Suppose that RΛ(λ− H(x) E)RΛ 1 ⊂ ∩ ∅ ∪ − and R (λ− H(x) E)R , i = 1, 2 are invertible. Then Λi − Λi 1 G = G + G λ− (G + G )(H H H )G . Λ Λ1 Λ2 − Λ1 Λ2 Λ − Λ1 − Λ2 Λ If m Λ and n Λ, we have ∈ 1 ∈ 1 ρ n′ n′′ (3.1) G (m, n) G (m, n) χ (n)+ λ− e− | − | G (m, n′) G (n′′, n) . | Λ | ≤ | Λ1 | Λ1 | Λ1 || Λ | n ΛX,n Λ ′∈ 1 ′′∈ 2 We remind

Lemma 3.1 (Schur test). Suppose A = Aij is a symmetric matrix. Then

A sup Aij . k k≤ i | | Xj We now prove 2 ρ Lemma 3.2. Let M0 (log N) ,¯ρ [ 2 , ρ] and M1 N. Let diam(Λ) 2N +1. Suppose that for any n Λ, there exists≥ some W ∈= W (n) with≤ M M M ≤such that n W Λ, ∈ ∈EM 0 ≤ ≤ 1 ∈ ⊂ dist(n,Λ W ) M and \ ≥ 2 (3.2) G (E; x) 2e√M , k W (n) k≤ ρ¯ n n′ M (3.3) G (E; x)(n, n′) 2e− | − | for n n′ . | W (n) |≤ | − |≥ 10 We assume further that N is large enough so that

3ρ ∞ ρ 1 √M d M d j 1 (3.4) sup 2λ− e (2M + 1) e− 20 (2j + 1) e− 2 . M M M ≤ 2 0≤ ≤ 1 Xj=0 8 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI

Then G (E; x) 4(2M + 1)de√M1 . k Λ k≤ 1 Proof. For simplicity, we drop the dependence on E and x. Under the assumption of (3.4), it is easy to check that for all M M M , 0 ≤ ≤ 1 1 d √M+ ρ M ρ n n 1 (3.5) 2λ− (2M + 1) e 10 e− 2 | − 2| . n Λ ≤ 2 X2∈ n n M | 2− |≥ 2 By (3.2) and (3.3), one has

ρ¯ √M+ M ρ¯ n n′ (3.6) G (n, n′) 2e 10 e− | − |. | W (n) |≤ For each n Λ, applying (3.1) with Λ1 = W (n), one has (3.7) ∈ 1 ρ n n G (n, n′) G (n, n′) χ (n′)+ λ− e− | 1− 2| G (n, n ) G (n , n′) . | Λ | ≤ | W (n) | W (n) | W (n) 1 || Λ 2 | n W (n) 1X∈ n Λ W (n) 2∈ \ By (3.6) and the fact that W (n) (2M + 1)d, one has | |≤ 1 √M+ ρ¯ M ρ¯ n n ρ n n G (n, n′) G (n, n′) χ (n′)+2λ− e 10 e− | − 1|e− | 1− 2| G (n , n′) | Λ | ≤ | W (n) | W (n) | Λ 2 | n W (n) 1X∈ n Λ W (n) 2∈ \ 1 d √M+ ρ M ρ n n G (n, n′) χ (n′)+2λ− (2M + 1) e 10 e− 2 | − 2| G (n , n′) ≤ | W (n) | W (n) | Λ 2 | n ΛXW (n) 2∈ \ ρ ρ 1 d √M+ M n n2 (3.8) GW (n)(n, n′) χW (n)(n′)+2λ− (2M + 1) e 10 e− 2 | − | GΛ(n2, n′) . ≤ | | n Λ | | X2∈ n n M | 2− |≥ 2 where the last inequality holds by the assumption dist(n,Λ W (n)) M . \ ≥ 2 Summing over n′ Λ in (3.8) and noticing (3.5) yields ∈ d √M 1 (3.9) sup G (n, n′) 2(2M + 1) e 1 + sup G (n , n′) . | Λ | ≤ 1 2 | Λ 2 | n Λ nXΛ n2 Λ nXΛ ∈ ′∈ ∈ ′∈ Now the lemma follows from Lemma 3.1.

1 ρ 4ρ 2d 2 Theorem 3.3. Let diam(Λ) 2N +1 and diam(Λ1) N . Let M0 (log N) ,¯ρ [ 2 , 5 ]. Suppose that for any n Λ Λ≤, there exists some W =≤ W (n) with≥ M M such∈ that ∈ \ 1 ∈ EM ≥ 0 n W , dist(n,Λ Λ W ) M , W Λ Λ and ∈ \ 1\ ≥ 2 ⊂ \ 1 G (E; x) e√M , k W k≤ ρ¯ n n′ M G (E; x)(n, n′) e− | − | for n n′ . | W |≤ | − |≥ 10 Suppose that G (E; x) e√N . k Λ k≤ Then

O(1) (¯ρ 1/2 ) n n′ − − M | − | N G (E; x)(n, n′) e 0 for n n′ . | Λ |≤ | − |≥ 10 9

1 Proof. As usual, we drop the dependence on E and x. Suppose n n′ N d + 1. Obviously, | − |≥ one of n and n′ is not in Λ1. By the self-adjointness of Green’s functions, we can assume n / Λ . ∈ 1 Applying (3.1) with Λ1 = W = W (n), one has

1 ρ n1 n2 (3.10) GΛ(n, n′) λ− e− | − | GW (n, n1) GΛ(n2, n′) . | |≤ n W | || | X1∈ n Λ W 2∈ \ It implies (since λ> 1)

ρ n n G (n, n′) e− | 1− 2| G (n, n ) G (n , n′) | Λ | ≤ | W 1 || Λ 2 | n W , nXn M 1 1∈ | 1− |≤ 10 − n Λ W 2∈ \ ρ n n + e− | 1− 2| G (n, n ) G (n , n′) | W 1 || Λ 2 | n W ,Xn n M 1∈ | 1− |≥ 10 n Λ W 2∈ \ √M ρ n n e e− | 1− 2| G (n , n′) ≤ | Λ 2 | n W , nXn M 1 1∈ | 1− |≤ 10 − n Λ W 2∈ \ ρ n n ρ¯ n n + e− | 1− 2|e− | − 1| G (n , n′) | Λ 2 | n W ,Xn n M 1∈ | 1− |≥ 10 n Λ W 2∈ \ √M ρ¯ n n e e− | − 2| G (n , n′) ≤ | Λ 2 | n W , nXn M 1 1∈ | 1− |≤ 10 − n Λ W 2∈ \ ρ¯ n n + e− | − 2| G (n , n′) | Λ 2 | n W ,Xn n M 1∈ | 1− |≥ 10 n Λ W 2∈ \ O(1) (¯ρ ) n n2 2d − − √M | − | (3.11) (2N + 1) sup e 0 GΛ(n2, n′) , ≤ n Λ W | | 2∈ \ where the third inequality holds because ofρ ¯ 4 ρ and n n M . ≤ 5 | − 2|≥ 2 1 2 n n′ N Iterating (3.11) until n2 n′ N 2 (but at most | − | times), we have n n′ , | − |≤ M0 | − |≥ 10 O(1) 1 O( n n ) 2 | − ′| (¯ρ )( n n′ N ) √ M − − √M0 | − |− N GΛ(n, n′) (2N + 1) 0 e e | | ≤ O(1) (¯ρ ) n n′ − − 1/2 | − | e M0 . ≤

Lemma 3.4 (Several variables matrix-valued Cartan estimate). Let T (x) be a self-adjoint N N matrix function of a parameter x [ δ, δ]J (J Z+) satisfying the following conditions: × ∈ − ∈ (i) T (x) is real analytic in x [ δ, δ]J and has a holomorphic extension to ∈ − CJ δ,δ1 = x = (xi)1 i J : sup xi δ, sup xi δ D ≤ ≤ ∈ 1 i J |ℜ |≤ 1 i J |ℑ |≤ ≤ ≤ ≤ ≤ satisfying

(3.12) sup T (x) B1, B1 1. x k k≤ ≥ ∈Dδ,δ 10 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI

(ii) For all x [ δ, δ]J , there is subset V [1, N] with ∈ − ⊂ V M, | |≤ and 1 (3.13) (R[1,N] V T (x)R[1,N] V )− B2, B2 1. k \ \ k≤ ≥ (iii) J 1 3J J J J J (3.14) mes x [ δ, δ] : T − (x) B 10− J − δ (1 + B )− (1 + B )− . { ∈ − k k≥ 3}≤ 1 2 Let 10M (3.15) 0 <ǫ (1 + B + B )− . ≤ 1 2 Then 1 1/J log ǫ− J 1 1 J c M log(B +B ) (3.16) mes x [ δ/2, δ/2] : T − (x) ǫ− Cδ e− 2 3 , n ∈ − k k≥ o ≤ where C = C(J, B1), c = c(J, B1) > 0. Proof. The proof is similar to that of the case J = 1 in Chapter 14 of [4] (see also Remark 3 there). We use the higher dimensional Cartan sets techniques of [13]. For convenience, we give the details in the Appendix.

Theorem 3.5. Under the assumptions of Theorem 2.7, let ω ΩN2 ΩN1 . We assume for j Tb 1 c3 c4 ¯∈ ∩ ¯ some x = (x , xj¬) , there exist N [ 4 N3 , N3 ] and Λ Λ N with diam(Λ) 1 ∈ ∈ ⊂ ∈ E ≤ 10N 10d such that, for any k Λ Λ¯, there exists some k W ,W Λ Λ¯ such that ∈ \ ∈ ∈ EN1 ⊂ \ dist(k,Λ Λ¯ W ) N1 , and x + kω mod Zb / X . Let \ \ ≥ 2 ∈ N1 b j ρN √N Y = y R j : y x e− 1 , G (E; (y, x¬)) e . { ∈ | − |≤ k Λ j k≥ } Then N 1/3bj (3.17) mes(Y ) e− . ≤ 1 Tb1 Tb b1 Proof. Without loss of generality, we assume j = 1. Fix x and x1¬ − . ρN 1 ∈ ∈ Let be the e− 1 neighbourhood of x in the complex plane, i,e., D b ρN 1 ρN = z C 1 : z e− 1 , z x e− 1 . D { ∈ |ℑ |≤ |ℜ − |≤ } 0 By the assumption of Theorem 3.5, one has for all k Λ Λ¯ and QN , ∈ \ 1 ∈EN1 √N1 (3.18) GQ (E; x + kω) e , k N1 k≤ ρ¯ n n′ N1 (3.19) GQ (E; x + kω)(n, n′) e− | − | for n n′ . | N1 |≤ | − |≥ 10 2 0 By standard perturbation arguments , (3.18) and (3.19), we have for any y , QN , ∈ D 1 ∈EN1 and k Λ Λ,¯ ∈ \ 1 √N1 (3.20) GQ (E; (x + y, x¬)+ kω) 2e , k N1 1 k≤ 1 ρ¯ n n′ N1 (3.21) GQ (E; (x + y, x¬)+ kω)(n, n′) 2e− | − | for n n′ . | N1 1 |≤ | − |≥ 10 Substituting Λ with Λ Λ¯ in Lemma 3.2, one has for any y , \ ∈ D 1 2√N1 (3.22) GΛ Λ¯ (E; (x + y, x1¬)) e . k \ k≤

2See e.g. the proof of Theorem 4.3. 11

We want to use Lemma 3.4. For this purpose, let 1 1 ρN T (y)= λ− H ((x + y, x¬)) E, J = b , δ = e− 1 . Λ 1 − 1 Now we are in the position to check the assumptions of Lemma 3.4. Obviously, B1 = O(1) since λ> 1 and E is bounded. Let V = Λ.¯ By (3.22), one has (3.23) M = Λ¯ 30dN 1/10, B = e2√N1 . | |≤ 2 By the fact that the Green’s functions satisfy property P and (2.7), one has that both (2.5) b N c1 and (2.6) hold at scale N for all y except a set of y T 1 with measure less than e− 2 . It 2 ∈ implies both (2.5) and (2.6) holds at scale N2 for all x + kω with k N3 except a set of d N c1 | | ≤ measure less than (2N3 + 1) e− 2 . Applying Lemma 3.2 with M0 = M1 = N2 and (2.1), one has 1 d √N 2√N T − (y) 4(2N + 1) e 2 4e 2 =: B , k k≤ 2 ≤ 3 c1 b1 d N except on a set of y T with measure less than (2N3 + 1) e− 2 . 2 ∈ c1 Since N2 = N1 , direct computation shows that c1 3b b1 b1 b b N /2 10− 1 b− δ (1 + B )− 1 (1 + B )− 1 e− 2 . 1 1 1 2 ≥ This veriﬁes (iii) in Lemma 3.4. √N For ǫ = e− , by (3.23), one has 10M ǫ< (1 + B1 + B2)− . By (3.16) of Lemma 3.4,

1/b1 √N c 1/3b − N N1/10 N 1 (3.24) mes(Y ) Ce 2 e− . ≤ ≤

Theorem 3.6. Let c1, c2, c3, c4, N1, N2, N3,Ω3 be given by Theorem 2.7, so in particular, Green’s functions satisfy property P at N1, N2 with parameters (c1,¯ρ). Then for all N3 2 O(1) ≤ N N3 , Green’s functions satisfy property P at size N with parameters (c1,¯ρ 1/2 ) and ≤ − N1 Ω =Ω Ω , where O(1) only depends on d. N 3 ∩ N2 2 0 Proof. We ﬁx N [N3, N3 ] and QN N . Let ω ΩN3 . ∈ ∈E ∈ c c For any n Q , replacing x with x + nω in Theorem 2.7, there exits N 3 < N¯

N˜ (3.28) n Λ Q , dist(n, Q Λ ) ∈ new ⊂ N N \ new ≥ 2 12 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI and

1 (3.29) Diam(Λ¯ ) 4N˜ 10d . new ≤

Also for any k Λ Λ¯ , there exists some W Λ Λ¯ such that ∈ new\ new EN1 ∋ ⊂ new\ new

N (3.30) dist(k,Λ Λ¯ W ) 1 . new\ new\ ≥ 2

We split the proof into three cases. Case 1: n +Λ QN . In this case, let Λnew = n + Λ and Λ¯ new = n + Λ.¯ See Case 1 of Fig.1. Case 2: (n + Λ)⊂ (Zd Q ) is non-empty and dist(n + Λ,¯ ∂Q ) 2N . See Case 2 of Fig.1. ∩ \ N N ≥ 1 In this case, let Λ¯ new = n + Λ¯ (the red square). By shrinking n + Λ a little bit, we can obtain proper Λnew (n + Λ) QN satisfying (3.28). Since dist(n + Λ,¯ ∂QN ) 2N1, we can also guarantee (3.30)⊂ holds. ∩ ≥

Λnew

Case 2 Case 1 Fig.1

Case 3: (n + Λ) (Zd Q ) is non-empty and dist(n + Λ,¯ ∂Q ) 2N . In this case, making ∩ \ N N ≤ 1 (n + Λ)¯ QN possibly larger, we obtain Λ¯ new QN . We can also make sure for any k Q Λ¯ ∩, there exists some W Q Λ¯ ⊂ ∈ N \ new ∈EN1 ⊂ N \ new

N (3.31) dist(k, Q Λ¯ W ) 1 . N \ new\ ≥ 2

See Fig.2. For B, Λ¯ = (n + Λ)¯ Q (the red part). For A and C, (n + Λ)¯ Q is the red new ∩ N ∩ N part, and Λ¯ new is union of the red part and the blue part. Shrinking n + Λ, we can obtain proper Λnew satisfying (3.28). This implies (3.30) by (3.31). 13

Λnew Λnew

A B

Λnew

Fig.2: Case 3

Tbj 2bj ρN1 ρN1 Fix xj¬. Divide into e cubes of size e− . Applying Theorem 3.5 in each cube, ((3.27), (3.29) and (3.30) ensure we can use Theorem

3.5), there exists a set YN˜ (xj¬) such that

1 3bj 2bj ρN1 N˜ (3.32) mes(Y ˜ (x¬)) e e− , N j ≤ j j and for x = (x , x¬) with x / Y ˜ (x¬), j ∈ N j j √N˜ (3.33) G (E; (x , x¬)) e . k Λnew j k≤ Setting M0 = N1,Λ=Λnew and Λ1 = Λ¯ new in Theorem 3.3 ((2.5), (2.6), (3.27), (3.29), (3.30) and (3.33) ensure we can use Theorem 3.3), we have for such x,

O(1) (¯ρ 1/2 ) n n′ ˜ − − N | − | N (3.34) G (E; x)(n, n′) e 1 for n n′ . | Λnew |≤ | − |≥ 10 Let

(3.35) BN (xj¬)= YN˜ (xj¬). 1 N c3 [N˜ N c4 4 3 ≤ ≤ 3 b b By (3.32), (3.35) and since c = c /4˜b, one has for any j and x¬ T − j , 1 3 j ∈ N c1 (3.36) mes(B (x¬)) e− . N j ≤ j 1 c3 c4 Suppose x / B (x¬). Applying Λ = Q , M = N and M = N in Lemma 3.2 since ∈ N j N 0 4 3 1 3 N [N , N 2] ((3.26), (3.28), (3.33) and (3.34) ensure the assumption of Lemma 3.2), one has ∈ 3 3 c d N c4 √N (3.37) G (E; x) 4(2N 4 + 1) e√ 3 e . k QN k≤ 3 ≤ 0 1 c3 Applying Λ = QN N , M0 = 4 N3 and Λ1 = in Theorem 3.3, by (3.33), (3.34) and (3.37), we have ∈E ∅

O(1) (¯ρ 1/2 ) n n′ − − N | − | N (3.38) G (E; x)(n, n′) e 1 for n n′ . | QN |≤ | − |≥ 10 14 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI

Let O(1) X = x Tb : (E, x) is not (¯ρ , N) good ,Ω =Ω Ω . N { ∈ − 1/2 } N 3 ∩ N2 N1 The theorem follows from (3.38), (3.37) and (3.36).

4. Large deviation theorem for Green’s functions and proof of Theorem 1.1 The main result of this section is the following large deviation theorem (LDT) for Green’s functions. Theorem 4.1 (LDT). There exist constants γ = γ(b, d) (0, 1), N = N (v, ρ, b, d) and ∈ 0 0 λ0 = λ0(v, ρ, b, d), such that for all N N0 and λ λ0, the Green’s functions satisfy property ρ ≥ ≥ P with parameters (γ, 2 ) at size N, and the corresponding semi-algebraic set ΩN satisfying b mes(T N N ΩN ) 0, \ ∩ ≥ 0 → as λ . → ∞ Compactness arguments and Theorem 8 in [20] immediately imply Lemma 4.2 (Lojasiewicz type Lemma). For E R, δ > 0, deﬁne ∈ X := x Tb : v(x) E <δ . { ∈ | − | } Then there are constants C(v), a(v) > 0 such that a(v) (4.1) sup mes(X(xj¬)) C(v)δ . b b 1 j d,x T − j ≤ ≤ ≤ j¬∈ Theorem 4.3. Let X be as in Lemma 4.2 and (4.2) X := x : x + nω mod Zb X . N ∈ n[N | |≤ Then we have d a(v) (4.3) sup mes(XN (xj¬)) C(v)(2N + 1) δ . b b 1 j d,x T − j ≤ ≤ ≤ j¬∈ Moreover, if 1 d (4.4) λ 2δ− (2N + 1) , ≥ then for any x / X , ω Tb, we have for Q 0 , ∈ N ∈ N ∈EN 1 (4.5) G (E; x) 2δ− , k QN k≤ 1 ρ n n′ (4.6) G (E; x)(n, n′) 2δ− e− | − |. | QN |≤ Proof. The bound (4.3) follows from Lemma 4.2 immediately. 0 1 Let x / XN and ﬁx QN N . Let A = λ− RQN SRQN , with the kinetic term S being given ∈ ∈E 1 by (1.1). Let B be diagonal part of the restriction of λ− H E on Q , namely, − N B = R (v(x + nω)δ Eδ )R . QN nn′ − nn′ QN By (4.2), one has min v(x + nω) E δ. n QN | − |≥ ∈ It leads to 1 1 (4.7) B− δ− . k k≤ 15

ρ n n′ Since S(n, n′) e− | − | for all n, n′, by Lemma 3.1 again, one has for N N(ρ, d), | |≤ ≥

1 ρ n n′ 1 d (4.8) A λ− sup e− | − | λ− (2N + 1) . k k≤ n QN ≤ ∈ n′XQN ∈ By (4.4),

1 1 AB− . k k≤ 2 Combining with (4.7) and (4.8), we have the following Neumann series expansion

1 1 s (4.9) G = B− ( AB− ) . QN − Xs 0 ≥ Thus one has

1 1 1 (4.10) G B− 2δ− . k QN k ≤ || k1 AB 1 ≤ − || − k It implies (4.5). In particular, (4.6) is also true for n = n′. For n = n′, by (1.1), (4.9) and the fact that B is diagonal, we have 6

1 s s ρ n k1 ρ k1 k2 ks 1 n′ G (n, n′) B− λ− δ− e− | − |− | − |−···| − − | | QN | ≤ k k Xs 1 k ≥ N | i|≤ 1 ρ n n′ sd s s δ− e− | − | (2N + 1) λ− δ− ≤ Xs 1 ≥ 1 ρ n n′ 2δ− e− | − |, ≤ where the last inequality holds by (4.4).

Proof of Theorem 4.1. Denote the relation between N1 and N3 in Theorem 2.7 by f, i.e., xc1 (n) f(x)= e . Let N0 = N0(v, ρ, b, d) be suﬃciently large. Denote by f (x) the nth iteration of f, namely, f (n)(x)= f(f(f( x ))). Let g(x)= f 2(x). Clearly, g(x) f(x+1) for large x. ··· ··· ≥ 1 N¯ 1/2 By letting δ = 2 e− and Theorem 4.3, since c1 < 1/2, the Green’s functions satisfy 4ρ ¯ Tb N¯ 1/2 ¯ d property P with parameters (c1, 5 ) for N0 N N and ΩN = if λ 4e (2N + 1) . ≤ ≤ 2 ≥ Theorem 3.6 allows us to proceed from scales N, N c1 to scales [f(N), g(N)]. Since we want to cover all scales, our initial step will consist of property P at the interval of scales [N1, f(N1)]. For this reason, we need to take N1 = log log λ. Initial step: For large λ, the Green’s functions satisfy property P with parameters (c1, ρ0) Tb 4ρ for all N0 N g(log log λ) and ΩN = , where ρ0 = 5 . Let ≤ ≤

i 4ρ O(1) (4.11) ρ = . i 5 − f (j)(log log λ)1/2 Xj=1

Applying Theorem 3.6 to N1 = log log λ, log log λ + 1, log log λ + 2, , f(log log λ), the Green’s functions satisfy property P with parameters (c , ρ ) for all g···(log log λ) N 1 1 ≤ ≤ 16 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI g(f(log log λ)) since g(x) f(x + 1). Moreover, ≥ f(log log λ) f(log log λ) 1 mes  ΩN  1 ≥ − f(N)c3 N=log\ log λ N=logX log λ   f(log log λ) 1 (4.12) 1 . ≥ − N 5 N=logX log λ

(2) Applying Theorem 3.6 to N1 = f(log log λ), f(log log λ)+1, f(log log λ)+2, , f (log log λ), the Green’s functions satisfy property P with parameters (c , ρ ) for all g(f(log··· log λ)) N 1 2 ≤ ≤ g(f (2)(log log λ)). Moreover,

f (2)(log log λ) f (2)(log log λ) 1 mes  ΩN  1 . ≥ − N 5 N=f(log\ log λ)+1 N=f(logX log λ)+1   By induction, we have the Green’s functions satisfy property P with parameters (c1, ρi) for (i 1) (i) g(f − (log log λ)) N g(f (log log λ)), i = 1, 2, . Moreover, ≤ ≤ ··· f (i) (log log λ) f (i) (log log λ) 1 (4.13) mes  ΩN  1 . ≥ − N 5 N=f (i 1)\(log log λ)+1 N=f (i 1)X(log log λ)+1  −  − Now Theorem 4.1 follows from (4.11) and (4.13).

Proof of Theorem 1.1. With Theorem 4.1 at hand, the proof Theorem 1.1 is rather standard. We refer the readers to [4, Section 3] or [7, Section 6] for details.

Appendix A. In the following, we will prove the several variables matrix-valued Cartan estimate, i.e., Lemma 3.4. The proof is similar to that in [3, 4]. Before going to the details, we recall some useful lemmas. The ﬁrst result is the standard Schur’s complement theorem. For convenience, we include a proof here. Lemma A.1. Let T be the matrix

T1 T2 T = t , T2 T3 where T is an invertible n n matrix , T is an n k matrix and T is a k k matrix. Let 1 × 2 × 3 × t 1 S = T T T − T . 3 − 2 1 2 Then T is invertible if and only if S is invertible, and

1 1 1 2 1 (A.1) S− T − C(1 + T − ) (1 + S− ), k k≤k k≤ k 1 k k k where C depends only on T . k 2k Proof. It is easy to check that

T1 T2 I 0 IT2 T1 0 (A.2) T = t = t 1 . T2 T3 T2 T1− I 0 S 0 I 17

It implies T is invertible if and only if S is invertible and also the second inequality of (A.1). By (A.2), one has 1 1 1 1 T1 0 − IT2 − I 0 − T − = 1 1 0 I 0 S T2− T1− I 1 1 T1− 0 I T2S− I 0 = − 1 1 1 0 I 0 S− T − T − I − 2 1 ⋆ ⋆ = 1 . ⋆ S− implying the first inequality of (A.1). We then introduce the higher dimensional Cartan sets Lemma of Goldstein-Schlag [13]. We denote by (z, r) the standard disk on C of center z and radius r> 0. D Lemma A.2. [13, Lemma 2.15] Let f(z , , z ) be an analytic function defined in a ploydisk 1 ··· J = (zi,0,1/2) and φ = log f . Let sup φ(z) M, m φ(z0), z0 = (z1,0, , zJ,0). P 1 i J D | | z ≤ ≤ ··· ≤Q≤ ∈P Given F 1, there exists a set such that ≫ B⊂P (A.3) φ(z) >M C(J)F (M m), for z (z ,1/4) , − − ∀ ∈ D i,0 \B 1 Yi J ≤ ≤ and J F 1/J (A.4) mes( R ) C(J)e− . B ∩ ≤ Proof of Lemma 3.4. The proof is similar to that of Proposition 14.1 in [4] in case J = 1 and Lemma 1.43 in [3] without explicit bounds. In the following proof, C = C(B1, J) and c = c(B1, J). Let 2 1 1 1 µ = 10− J − δ(1 + B1)− (1 + B2)− . Fix x [ δ/2, δ/2]J 0 ∈ − and consider T (z) with z x0 = sup zi x0,i < µ. Thanks to Cauchy’s estimate and | − | 1 i J | − | ≤ ≤ (3.12), one obtains for z x <µ, | − 0| 4B ∂ T (z) 1 , i = 1, 2, , J, k zi k≤ δ ··· which implies 4JB1µ 1 1 T (z) T (x ) 25− (1 + B )− . k − 0 k≤ δ ≤ 2 From the assumption (ii) of Lemma 3.4, we can find V = V (x0) so that V M and (3.13) is satisfied. Denote by V c = [1, N] V . Thus using the standard Neumann| series|≤ argument and (3.13), one has \ 1 (A.5) (R c T (z)R c )− 2B for z x <µ. k V V k≤ 2 | − 0| We define for z x <µ the analytic self-adjoint function | − 0| 1 (A.6) S(z)= R T (z)R R T (z)R c (R c T (z)R c )− R c T (z)R . V V − V V V V V V Then by (A.5) and (A.6), we have (A.7) S(z) 3B2B . k k≤ 1 2 18 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI

Recalling Lemma A.1, if S(z) is invertible, so is T (z) and by (A.1), 1 1 2 1 (A.8) S− (z) C T − (z) CB (1 + S− (z) ). k k≤ k k≤ 2 k k For x RJ , one has ∈ M 1 M (A.9) S(x) det S(x) = λ S− (x) − . || || ≥ | | | |≥k k λ σY(S(x)) ∈ By (A.7), one has M 1 2 M 1 S(x) − (3B1 B2) (A.10) S− (x) k k . k k≤ det S(x) ≤ det S(x) | | | | Let φ(z) = log det S(x + µz) , z < 1. | 0 | | | Then by (A.9) and (A.7),

(A.11) sup φ(z) CM log B2. z <1 ≤ | | By (3.14) and the deﬁnition of µ, there is some x with x x < µ/10 such that 1 | 0 − 1| 1 (A.12) T − (x ) B . k 1 k≤ 3 1 Hence by (A.8), S− (x ) CB , and from (A.9), k 1 k≤ 3 (A.13) φ(a) CM log B , ≥− 3 x x where a = 1− 0 , so a < 1/10. Let µ | | = (a ,1/2). P D i 1 Yi J ≤ ≤ Then one has sup φ(z) CM log B2, φ(a) CM log B3. z ≤ ≥− ∈P Applying Lemma A.2 and recalling (A.3), (A.4), for any F 1, there is some set ≫ B ⊂ (ai,1/4) with 1 i J D ≤Q≤ (A.14) φ(z) CFM log(B + B ) for z (a ,1/4) , ≥− 2 3 ∈ D i \B 1 Yi J ≤ ≤ and J F 1/J (A.15) mes( R ) Ce− . B ∩ ≤ For 0 <ǫ< 1, let c log ǫ F = − . M log(B2 + B3) Then by (A.14) and (A.15), mes x RJ : x x < µ/4 and det S(x) ǫ ∈ | − 1| | |≤ = µJ mes x RJ : x a < 1/4 and φ(x) log ǫ ∈ | − | ≤ J F 1/J Cµ e− . ≤ Since x x < µ/10, we have | 0 − 1| 1 1/J log ǫ− J J c M log(B +B ) (A.16) mes x R : x x < µ/8 and det(S(x)) ǫ Cµ e− 2 3 . ∈ | − 0| | |≤ ≤ 19

Recalling (A.8), (A.10) and (3.15), one has for x x < µ/8 and det S(x) ǫ, | − 0| | |≥ 1 2 1 2 M 2 (A.17) T − (x) C(1 + B )(1 + ǫ− (3B B ) ) Cǫ− . k k≤ 2 1 2 ≤ Covering [ δ , δ ]J by cubes of side µ/4, and combining (A.16) and (A.17), one has − 2 2 1 1/J log ǫ− J 1 2 J c M log(B +B ) mes x [ δ/2, δ/2] : T − (x) ǫ− Cδ e− 2 3 . n ∈ − k k≥ o ≤

Acknowledgments We are grateful to Jean Bourgain for his encouragement. This research was supported by NSF DMS-1401204, DMS-1901462, and DMS-1700314.

References [1] A. Avila. Global theory of one-frequency Schrödinger operators. Acta Math., 215(1):1–54, 2015. [2] A. Avila, J. You, and Q. Zhou. Sharp phase transitions for the almost Mathieu operator. Duke Math. J., 166(14):2697–2718, 2017. [3] J. Bourgain. Estimates on Green’s functions, localization and the quantum kicked rotor model. Ann. of Math. (2), 156(1):249–294, 2002. [4] J. Bourgain. Green’s function estimates for lattice Schrödinger operators and applications, volume 158 of Annals of Mathematics Studies. Princeton University Press, Princeton, NJ, 2005. [5] J. Bourgain. Anderson localization for quasi-periodic lattice Schrödinger operators on Zd, d arbitrary. Geom. Funct. Anal., 17(3):682–706, 2007. [6] J. Bourgain and M. Goldstein. On nonperturbative localization with quasi-periodic potential. Ann. of Math. (2), 152(3):835–879, 2000. [7] J. Bourgain, M. Goldstein, and W. Schlag. Anderson localization for Schrödinger operators on Z2 with quasi-periodic potential. Acta Math., 188(1):41–86, 2002. [8] J. Bourgain and S. Jitomirskaya. Absolutely continuous spectrum for 1D quasiperiodic operators. Invent. Math., 148(3):453–463, 2002. [9] J. Bourgain, S. Jitomirskaya, and L. Parnovski. Absolutely continuous spectrum for multidimensional quasiperiodic operators. In preparation. [10] J. Bourgain and I. Kachkovskiy. Anderson localization for two interacting quasiperiodic particles. Geom. Funct. Anal., 29(1):3–43, 2019. [11] V. A. Chulaevsky and E. I. Dinaburg. Methods of KAM-theory for long-range quasiperiodic operators on Zν . Pure point spectrum. Comm. Math. Phys., 153(3):559–577, 1993. [12] J. Fröhlich and T. Spencer. Absence of diffusion in the Anderson tight binding model for large disorder or low energy. Comm. Math. Phys., 88(2):151–184, 1983. [13] M. Goldstein and W. Schlag. Fine properties of the integrated density of states and a quantitative separation property of the Dirichlet eigenvalues. Geom. Funct. Anal., 18(3):755–869, 2008. [14] S. Jitomirskaya and W. Liu. Universal hierarchical structure of quasiperiodic eigenfunctions. Ann. of Math. (2), 187(3):721–776, 2018. [15] S. Jitomirskaya and W. Liu. Universal reflective-hierarchical structure of quasiperiodic eigenfunctions and sharp spectral transition in phase. arXiv preprint arXiv:1802.00781, 2018. 20 SVETLANAJITOMIRSKAYA,WENCAILIU,ANDYUNFENGSHI

[16] S. Jitomirskaya and S. Zhang. Quantitative continuity of singular continuous spectral measures and arithmetic criteria for quasiperiodicSchr¨odinger operators. arXiv preprint arXiv:1510.07086, 2015. [17] S. Y. Jitomirskaya. Anderson localization for the almost Mathieu equation: a nonperturbative proof. Comm. Math. Phys., 165(1):49–57, 1994. [18] S. Y. Jitomirskaya. Metal-insulator transition for the almost Mathieu operator. Ann. of Math. (2), 150(3):1159–1175, 1999. [19] C. A. Marx and S. Jitomirskaya. Dynamics and spectral theory of quasi-periodic Schr¨odinger-type operators. Ergodic Theory Dynam. Systems, 37(8):2353–2393, 2017. [20] D. H. Phong, E. M. Stein, and J. A. Sturm. On the growth and stability of real-analytic functions. Amer. J. Math., 121(3):519–554, 1999. [21] J. You. Quantitative almost reducibility and its applications. In International Congress of Mathematicians, Rio de Janeiro, page 1916, 2018.

(S. Jitomirskaya) Department of Mathematics, University of California, Irvine, CA, USA E-mail address: [email protected]

(W. Liu) Department of Mathematics, University of California, Irvine, CA, USA E-mail address: [email protected]

Current address: Department of Mathematics, Texas A&M University, College Station, TX, USA

(Y. Shi) School of Mathematical Sciences, Peking University, Beijing 100871, China E-mail address: [email protected]