Access the full text.

Sign up today, get DeepDyve free for 14 days.

Mathematics
, Volume 2020 (1711) – Nov 19, 2017

/lp/arxiv-cornell-university/eigenvectors-distribution-and-quantum-unique-ergodicity-for-deformed-5xsee0239I

- ISSN
- 0246-0203
- eISSN
- ARCH-3343
- DOI
- 10.1214/20-AIHP1060
- Publisher site
- See Article on Publisher Site

We analyze the distribution of eigenvectors for mesoscopic, mean-ﬁeld perturbations of diagonal matrices, in the bulk of the spectrum. Our results apply to a generalized N × N Rosenzweig– Porter model. We prove that the eigenvector entries are asymptotically Gaussian with a speciﬁc variance. For a well spread initial spectrum, this variance proﬁle universally follows a heavy-tailed Cauchy distribution. In the case of smooth entries, we also obtain a strong form of quantum unique ergodicity in the form of a strong concentration inequality for the mass of eigenvectors on a given set of coordinates. The proof relies on a priori local laws for this model as given in [15, 35, 37], and the eigenvector moment ﬂow from [15, 16]. Résumé Nous analysons la distribution des vecteurs propres de perturbations mésoscopiques de matrices diagonales à l’intérieur du spectre. Nos résultats s’appliquent a un modèle généralisé de Rosenzweig– Porter. Nous prouvons que les entrées des vecteurs propres sont asymptotiquement gaussiennes avec une variance explicite. Lorsque le spectre initial est bien étalé, ce proﬁl de variance suit de manière universelle une distribution de Cauchy à queue lourde. Lorsque les entrées sont lisses, nous obtenons aussi une forme forte d’unique ergodicité quantique sous la forme d’une inégalité de concentration sur la masse des vecteurs propres sur un domaine ﬁxé de coordonnées. La preuve se base sur des lois locales a priori données dans [15, 35, 37] et le ﬂot des moments des vecteurs propres de [15, 16] Keywords: Deformed random matrices; Eigenvector distribution; Quantum unique ergodicity; Dyson Brownian motion MSC: Primary 60B20; secondary 58J51 Contents 1 Introduction 1 2 Local laws 10 3 Short time relaxation 14 4 Approximation by a Gaussian divisible ensemble 35 5 Proofs of main results 43 1. Introduction In the study of large interacting quantum systems, Wigner conjectured that empirical results on energy levels are well approximated by statistics of eigenvalues of large random matrices. We refer to [39] for an overview and the mathematical formalization of the conjecture. This vision has not been shown for realistic correlated quantum systems but is regarded to hold for numerous models. For instance, the Bohigas−Giannoni−Schmit conjecture in quantum chaos [10] connects eigenvalues distributions in the semiclassical limit to the Gaudin distribution for GOE statistics. These statistics also conjecturally appear for random Schrödinger operators [6] in the delocalized phase. Most of these hypotheses are unfortunately far from being proved arXiv:1711.07103v3 [math.PR] 3 Nov 2020 with mathematical rigor. It is, however, possible to study systems given by large random matrices. One of the most important models of this type is the Wigner ensemble, random Hermitian or symmetric matrices whose elements are, up to the symmetry, independent and identically distributed zero-mean unit variance random variables. For this ensemble, local statistics of the spectrum only depend on the symmetry class and not on the laws of the elements (see [12, 23, 25, 27, 49]). The Wigner−Dyson−Mehta conjecture was solved for numerous, more general mean-ﬁeld models such as the generalized Wigner matrices, random matrices for which the laws of the matrix elements can have distinct variances (see [26] and references therein). The statistics of eigenvectors were not present in Wigner’s original study but localization, or delocalization, has been broadly studied in random matrix theory. For Wigner matrices, it has been shown in [22] that eigenvectors are completely delocalized in the following sense: denoting u , . . . , u the L -normalized eigenvectors of an N × N Wigner matrix, we have 1 N with very high probability, C (log N ) sup|u (α)| 6 √ . α N −C Thus, eigenvectors cannot concentrate onto a set of size smaller than N (log N ) . See also [53] for optimal bounds in some cases of Wigner matrices or [28] for similar estimates for gen- eralized Wigner matrices and [45] for an improved bound which also holds for non-Hermitian matrices. In the GOE and GUE cases, the distribution of the matrix is orthogonally invariant and eigenvectors are distributed according to the Haar measure on the orthogonal group. In particular, the entries of bulk eigenvectors are asymptotically normal: Nu (α) −−−−→ N, N→∞ where N is a standard Gaussian random variable. Asymptotic normality was ﬁrst proved for Wigner matrices in [33,50] under a matching condition on the ﬁrst four moments of the entries using comparison theorems introduced in [27, 49]. These conditions were later removed in [15] where asymptotic normality holds for generalized Wigner matrices. Beyond mean-ﬁeld models, conjectures of interest, for example for band matrices, are still yet to be proved. A sharp transition is conjectured to occur when the band width W cross the critical value √ √ N. For W N, eigenvectors are expected to be exponentially localized on O(W ) sites and eigenvalue statistics are Poisson, while for W N eigenvectors would be completely delocalized and one would get Wigner-Dyson-Mehta statistics for the eigenvalues. For the most recent works on this subject see [42, 46] for localization results, [16] for delocalization results, [48] for another transition occurring at the edge of the spectrum and [11] for a recent review on the subject. In this paper, we consider a generalized Rosensweig-Porter model, of mean-ﬁeld type, which also interpolates between delocalized and localized (or partially delocalized) phases, but al- ways with GOE/GUE statistics. It is deﬁned as a perturbation of a potential, consisting of a deterministic diagonal matrix, by a mean ﬁeld noise, given by a Wigner random matrix, scaled by a parameter t. This model follows two distinct phase transitions. When t 1/N, eigenvalue statistics coincide with t = 0 and eigenvectors are localized on O(1) sites [51], while when t 1, local statistics fall in the Wigner-Dyson-Mehta universality class [37] with fully delocalized eigenvectors [36]. For 1/N t 1, it has been shown in [35] that eigen- value statistics are in the Wigner-Dyson-Mehta universality class and in [52] that eigenvectors are not completely delocalized when the noise is Gaussian. In this intermediate phase, also called the bad metal regime (see [29] or [41] for instance), eigenstates are partially delocalized over Nt sites, a diverging number as N grows but a vanishing fraction of the eigenvector coordinates. The existence of this regime for more intricate models is only conjectured or even debated in the physics literature though progress has been made recently, for instance for the Anderson model on the Bethe lattice and regular graph in [34]. Our results give the asymptotic distribution of the eigenvectors for the Rosensweig-Porter model with a complete understanding of the bad metal regime. We show that bulk eigenvec- tors are asymptotically Gaussian with a speciﬁc, explicit variance depending on the initial potential, the parameter t and the position in the spectrum. For a well-spread initial con- dition, this variance is heavy-tailed and follows a Cauchy distribution. This shape was ﬁrst 2 unearthed in a non-rigorous way in [2,3] for W a matrix from the Gaussian ensembles, where the Gaussian distribution of eigenvectors (Corollary 1.4) was conjectured. Note that eigen- vector dynamics was also considered in [1] and used for denoising matrices in [18, 19]. In the case of Gaussian entries, the eigenvector distribution has been exhibited in the physics literature in [29] using the resolvent ﬂow and in [41] using supersymmetry techniques. Another strong form of delocalization of eigenfunctions is quantum ergodicity. It has been proved for the Laplace-Beltrami operator for a wide class of manifolds by Shnirel’man [47], Colin de Verdière [20] and Zelditch [54] but also for regular graphs by Ananthamaran- Le Masson [4]. In [44], Rudnick-Sarnack conjectured a stronger form of delocalization for eigenfunctions of the Laplacian called the quantum unique ergodicity. More precisely, de- note (φ ) the eigenfunctions of the Laplace operator on any negatively curved compact k k>1 Riemannian manifold M, they then supposedly become equidistributed with respect to the volume measure μ in the following sense: for any open set A ⊂ M Z Z |φ | dμ −−−−→ dμ. k→∞ A A This convergence has been established for arithmetic surfaces (see [30, 31, 38]). A probabilistic form of quantum unique ergodicity exists for eigenvectors of large random matrices. It ﬁrst appeared in [15] for generalized Wigner matrices. It is stated as a high- probability bound showing that eigenvectors are asymptotically ﬂat in the following way: let (u ) be the eigenvectors of a N×N generalized Wigner matrix, then for any k ∈ [[1, N ]], k 16k6N for any deterministic N-dependent set I ∈ [[1, N ]] such that |I| → +∞ and any δ > 0, −ε N 1 N P u (α) − > δ 6 |I| N δ α∈I for some ε > 0 using the Bienaymé–Chebyshev inequality. Similar high-probability bounds were proved for diﬀerent models of random matrices such as d-regular random graphs in [7], or band matrices in [13,16]. In these last papers on band matrices, it was seen that quantum unique ergodicity is a useful property to study non mean-ﬁeld models. In [16], a stronger form of probabilistic quantum unique ergodicity has been actually proved, showing that the eigenvectors mass is asymptotically ﬂat with overwhelming probability (the probability decreases faster than any polynomial). Our result adapts the method introduced in [16] to show a strong deformed quantum unique ergodicity for eigenvectors of a class of deformed Wigner matrices. Indeed, the probability mass is not ﬂat but concentrates onto an explicit and deterministic proﬁle with a quantitative error. The key ingredient for this analysis is the Bourgade-Yau eigenvector moment ﬂow [15], a multi-particle random walk in a random environment given by the trajectories of the eigen- values. This method was used for generalized Wigner matrices [15] and sparse random graphs [14], and both settings correspond to equilibrium or close to equilibrium situations. Our main contribution consists in treating the non-equilibrium case, which implies additional diﬃculties made explicit in the next section. 1.1. Main Results Consider a deterministic diagonal matrix D = diag(D , . . . , D ). The eigenvalues (or diago- 1 N nal entries) need to be regular enough on a window of size r in the following way ﬁrst deﬁned in [35]. Deﬁnition 1.1. Let η and r be two N-dependent parameters satisfying 0 0 0 −1 −ε ε −ε N 6 η 6 N , N η 6 r 6 N ? ? for some ε > 0. A deterministic diagonal matrix D is said to be (η , r)-regular at E if there ? 0 exists c > 0 and C > 0 independent of N such that for any E ∈ [E − r, E + r] and D D 0 0 η 6 η 6 10, we have c 6 Im m (E + iη) 6 C , D D D 3 where m is the Stieltjes transform of D: 1 1 m (z) = . N D − z k=1 We want to study the perturbation of such a diagonal matrix, notably the eigenvectors, by a mesoscopic Wigner random matrix. We will now suppose that D is (η , r)-regular at E , a ﬁxed energy point. Letting 0 < κ < 1, we will denote in the rest of the paper the spectral window as I = [E − (1− κ)r, E + (1− κ)r]. 0 0 We remove a certain window of energy to avoid any possible complications at the edge. We will also use the following domains: the ﬁrst domain is used for the size of our deformation while the other is the spectral domain in which we perform our analysis. We need the perturbation to be mesoscopic but smaller than the energy window size r, deﬁne then for any small positive ω, ω −ω T = N η , N r . ω ? For the spectral domain, take ﬁrst t ∈ T and note that we will consider only Im(z) := η smaller than t but most results such as local laws holds up to macroscopic η. Let ϑ > 0 be an arbitrarily small constant and ϑ,κ κ ϑ −ϑ D = z = E + iη : E ∈ I , N /N 6 η 6 N t . r r Hereafter is our assumptions on our Wigner matrix. Deﬁnition 1.2. A Wigner matrix W is a N×N Hermitian/symmetric matrix satisfying the following conditions (i) The entries (W ) are independent. i,j 16i6j6N 2 −1 (ii) For all i, j, E[W ] = 0 and E[|W | ] = N . i,j ij (iii) For every p ∈ N, there exists a constant C such that NW 6 C . p ij p Let W be a Wigner matrix and deﬁne the following t-dependent matrix for t ∈ T W = D + tW. (1.1) The eigenvectors of D are exactly the vectors of the canonical basis since the matrix is diagonal. However, if t were of order one instead of being in T , the local statistics of W ω t would become universal and would be given by local statistics from the Gaussian ensemble. In particular, the eigenvectors would be completely delocalized [37]. Our model consists in looking at the diﬀusion of the eigenvectors on the canonical basis after a mesoscopic perturbation. Our main result is that the coordinates of bulk eigenvectors are time and position dependent Gaussian random variables. Before stating our result, we ﬁrst deﬁne the asymptotic distribution of the eigenvalues of the matrix W which is the free convolution of the semicircle law (coming from W ) and the empirical distribution of D. We will deﬁne this distribution through its Stieltjes transform m (z) as the solution to the following self- consistent equation 1 1 m (z) = . (1.2) N D − z − tm (z) α t α=1 It is known that this equation has a unique solution with positive imaginary part and is the Stieltjes transform of a measure with density denoted by ρ (see [9] for more details). Deﬁne the quantiles (γ ) of this measure by i,t 06i6N i,t ρ (x)dx = . −∞ 4 We can also now deﬁne the set of indices in the spectral window corresponding to the indices such that the corresponding classical locations lies in the energy window I , κ κ A = {i ∈ [[1, N ]], γ ∈ I } . r i,t r We can now state our main results, denoting u (t), . . . , u (t) the L -normalized eigenvectors 1 N of W (we will often omit the t-dependence for u). We will deﬁne the following quantity, for −1 N η t and η t r: q t 2 α σ (q, k, η) = . (1.3) (D − γ − tRe m (γ + iη)) + (tIm m (γ + iη)) α k ,t t k ,t t k ,t i i i α=1 It is the asymptotic deterministic variance of our eigenvector projections. By taking q to be a canonical basis vector e , we see that u (α) has a variance of the form α k 1 t 2 2 N (D − γ − tRe m (γ )) + (tIm m (γ )) α k ,t t k ,t t k ,t i i i For regularly spaced D ’s, this is heavy-tailed with Cauchy shape . It localizes the α 2 2 N (x +t ) entries onto a subset of indices of size Nt N: a fraction of the eigenvector coordinates vanishing as N grows. Such a partial localization appears in [52] for W GOE-distributed. 0 0 Theorem 1.3. (Gaussianity of bulk eigenvectors) Fix κ ∈ (0, 1), ω ∈ (0, ε /10) where ε is as in Deﬁnition 1.1 and m ∈ N. Let t ∈ T and I ⊂ A be a deterministic (N-dependent) set of m elements. Let W as in Deﬁnition 1.2 and W as in (1.1). Write I = {k , . . . , k }, t 1 m take a deterministic q ∈ R such that kqk = 1, and deﬁne for i ∈ [[1, m]], 2 2 σ (q, k ) := lim σ (q, k , η) (1.4) i i t t η↓0 Then we have |hq, u i| −−−−→ (|N |) in the symmetric case, (1.5) k i i i=1 N→∞ σ (q, k ) i=1 2N (1) (2) |hq, u i| −−−−→ (|N + iN |) in the Hermitian case (1.6) i i i i=1 σ (q, k ) N→∞ i=1 (1) (2) in the sense of convergence of moments, where all N , N and N are independent Gaus- i i sian random variables with variance 1. The convergence is uniform in over the choice of sets I ⊂ [[1, N ]] of size m. One can deduce joint weak convergence of eigenvector entries from the previous conver- N−1 gence of moments because q is arbitrary in S (see [15, Section 5.3]). However, since the eigenvectors are deﬁned up to a phase, we ﬁrst need to deﬁne the following equivalence iω relation: u ∼ v if and only if u = ±v in the symmetric case and u = e v for some ω ∈ R in the Hermitian case. Corollary 1.4. Let κ ∈ (0, 1) and m ∈ N, let W as in Deﬁnition 1.2 and W as in Deﬁnition 1.1. Then for any deterministic k ∈ A and J ⊂ [[1, N ]] such that |J| = m we have s ! u (α) −−−−→ (N ) in the symmetric case, (1.7) k i 2 i=1 σ (e , k) N→∞ α∈J s ! 2N (1) (2) u (α) −−−−→ N + iN in the Hermitian case (1.8) i i 2 2 σ (e , k) N→∞ i=1 α∈J (1) (2) in the sense of convergence of moments modulo ∼, where all N , N and N are indepen- i i dent Gaussian random variables with variance 1. In more precise terms, for any polynomial P in m variables there exists δ depending on P such that, for N large enough, " ! !# h i m −δ sup E P u (α) − E P (N ) 6 N k j 2 j=1 σ (e , k) J⊂[[1,N ]],|J|=m α κ α∈J k∈A 5 where is taken uniformly at random in the set {−1, 1}. The convergence in the Hermitian case is similar by taking uniform on the circle. This result states that the entries of bulk eigenvectors are asymptotically independent Gaussian random variables with variance σ which answers a conjecture from [3, Section 3.2], stated in the more restrictive case where W is GOE. The asymptotic normality of the eigenvectors gives the following weak form of quantum unique ergodicity. Corollary 1.5. (Weak Quantum Unique Ergodicity) Let W as in 1.2 and W as in Deﬁnition (1.1). There exists ϑ > 0 such that for any c > 0, there exists C > 0 such that the following holds : for any I ⊂ [[1, N ]] and k ∈ A , we have X X Nt 1 2 2 −ϑ −1 P |u (α)| − σ (e , k) > c 6 C (N +|I| ). (1.9) k α |I| N α∈I α∈I This high probability bound is not the strongest form of quantum unique ergodicity one can obtain for random matrices. Indeed, if we consider the Gaussian ensembles for which the eigenbasis is Haar-distributed on the orthogonal group and each eigenvectors is uniformly distributed on the sphere, one can get that for any ε and D positive constants, for any 1 6 k 6 N, |I| I 2 ε −D P |u (α)| − > N 6 N for any N suﬃciently large N N α∈I In this paper, we will obtain a similar overwhelming probability bound on the probability mass of a single eigenvector with an explicit error for a more restrictive model of matrices: deformed random matrix with smooth entries given by the following deﬁnition or from a Gaussian divisible ensemble. Deﬁnition 1.6. A smooth Wigner matrix W is a N × N Hermitian/symmetric matrix with the following conditions (i) The matrix entries (W ) are independent and identically distributed random ij 16i6j6N −1/2 variables following the distribution N ν where ν has mean zero and variance 1. −Θ(x) (ii) The distribution ν has a positive density ν(x) = e such that for any j, there are constants C and C such that 0 1 (j) 2 C |Θ (x)| 6 C (1 + x ) (1.10) (iii) The tail of the distribution ν has a subexponential decay. In other words, there exists C and q two positive constants such that 1 dν(x) 6 C exp(−y ) (1.11) |x|>y We need the smoothness assumptions on W in order to use the reverse heat ﬂow techniques from [21,23]. Indeed, our result is an overwhelming probability bound on the eigenvectors of W . We think, however, that this property holds for a larger matrix ensembles and that the smoothness property is simply technical. In the following, since the eigenvectors are concentrated on Nt sites, it is relevant to deﬁne the following notation for any set (which can be N-dependent) A, denote |A| A = ∧ 1. Nt Indeed, having errors involving A allows us to get bounds improving for |A| 6 Nt but still holding for |A| Nt. Theorem 1.7. Let κ ∈ (0, 1), ω a small positive constant, for instance ω < ε /10. Let t ∈ T , I ⊂ [[1, N ]] be a deterministic (N-dependent) set, W as in Deﬁnition 1.6 and W as ω t in (1.1). Deﬁne now b b I I 2 2 Ξ = and σ (α, k) := σ (α, k, η ) with Nη = . (1.12) 0 0 t t 1/3 2 (Nt) Ξ 6 κ Then we have, for any ε > 0 (small) and D > 0 (large) and for k, ` ∈ A with k 6= `, in the symmetric case X X 2 2 ε −D P u (α) − σ (α, k) + u (α)u (α) > N Ξ 6 N (1.13) k k ` α∈I α∈I and in the Hermitian case, X X 2 2 ε −D P |u (α)| − σ (α, k) + u (α)u ¯ (α) > N Ξ 6 N . (1.14) k k ` 2N α∈I α∈I Remark 1.8. The choice of η depends on our proof and is the one we should take to optimize our error Ξ. However, this error and the choice of η do not seem optimal since we actually expect to have some form of Gausian ﬂuctuations around this deterministic proﬁle. 1.2. Method of Proof Our proof is based on the three-step strategy from [21,23] (see [26] for a recent book presenting this method). The ﬁrst step is to have an optimal, local control of the spectral elements of the matrix ensemble given by a local law on the resolvent. The second step is to obtain the wanted result for a relaxation of the model by a small Gaussian perturbation. Finally, the third and last step consists of removing this Gaussian part. We will give the proof of Theorem 1.3, Corollary 1.5, Theorem 1.7 only in the symmetric case and refer the reader to [15, 16] for the tools needed in the Hermitian case. First step: local laws for our model. In [35], Landon-Yau showed a local law for the Dyson Brownian motion with a diagonal initial condition at all times. This result gives us an averaged local law on the Stieltjes transform but also an entrywise anisotropic local law for the resolvent. Since we want to consider any projection of the eigenvectors, we will also need a local law on the quadratic form hq, G(z)qi. This control of the resolvent for mesoscopic perturbation has been showed in [14]. Note that these results were done in the Gaussian case but can easily be generalized to the Wigner case with the right assumptions on moments. Second step: short time relaxation. The second step consists of perturbing W by a small Gaussian component. We will obtain this perturbed model by making W undergo the Dyson Brownian motion given by the following deﬁnition. Deﬁnition 1.9. Here is our choice of Dyson Brownian motion. Let B be a N × N symmetric matrix such that B for i < j and B / 2 are independent ij ii standard brownian motions. The N × N symmetric Dyson Brownian motion with initial condition H is deﬁned as H = H + √ B . (1.15) s 0 s We also give the dynamics followed by the eigenvalues and the eigenvectors of such ma- trices. Deﬁnition 1.10. Let λ be in the simplex Σ = {λ < ··· < λ }, u be an orthogonal 0 N 1 N 0 N × N matrix, and B as in (1.15). Consider the dynamics dB 1 ds kk dλ = √ + , (1.16) N λ − λ N k ` `6=k X X 1 dB 1 ds k` du = √ u − u (1.17) k ` k λ − λ 2N (λ − λ ) N k ` k ` `6=k `6=k with initial condition (λ ,u ), 0 0 This eigenvector ﬂow was ﬁrst computed in diﬀerent contexts such as [5] for GOE or GUE matrices, [17] for real Wishart processes and [40] for Brownian motion on ellipsoids. 7 Remark 1.11. If λ and u are the eigenvalues and eigenvectors of a ﬁxed matrix H , then 0 0 0 the solution to the dynamics from Deﬁnition 1.10 have the same distribution for any time s as the eigenvalues and eigenvectors of H = H + sGOE s 0 with GOE being a matrix from the normalized Gaussian Orthogonal Ensemble (in the sense that its oﬀ-diagonal entries have variance 1/N). In this paper, taking W to be such a matrix in (1.1), we study the eigenvectors of the Dyson Brownian motion with a diagonal initial condition after a mesoscopic time. −1 We will then need to study the eigenvectors of H for a small N τ t. The convergence of joint moments of eigenvectors projections will be obtained by the maximum principle technique introduced in [15]. It is based on analyzing the dynamics followed by these moments. We will now recall notations and results on this eigenvector moment ﬂow. Take q ∈ R such that kqk = 1 a ﬁxed direction onto which we will project our eigen- H H vectors. For u , . . . , u the eigenvectors of the matrix (1.15), deﬁne z (s) = Nhq, u (s)i. (1.18) Now for m ∈ [[1, N ]], denote by j , . . . , j positive integers and let i , . . . , i in [[1, N ]] be 1 m 1 m distinct indices. We will consider the following normalized polynomials Y Y j ,...,j 2j −1 1 m l Q = z a(2j ) where a(n) = k (1.19) i ,...,i i 1 m l l=1 k≤n,k odd 2n Note that a(2n) = E[N ] with N a standard Gaussian random variables. Consider a conﬁguration of particles ξ : [[1, N ]] → N where ξ := ξ(j) is seen as the number of particles at the site j. We denote N (ξ) = ξ the total number of particles in the conﬁguration ξ. i,j Deﬁne ξ to be the conﬁguration obtained by moving one particle from i to j. If there is i,j no particle in i then ξ = ξ. It is clear that we can map {(i , j ), . . . , (i , j )} with distinct 1 1 m m i ’s and positive j ’s summing to an n > 0 to a conﬁguration ξ with ξ = j and η = 0 if k k i k l l ∈/ {i , . . . , i }. 1 m Deﬁne now, given this map, h i j ,...,j 1 m f (ξ) := E Q |λ , (1.20) λ,s i ,...,i 1 m a n-th joint moment of the coordinates of u . The conditioning here is on the full path of eigenvalues from 0 to ∞. The next theorem gives the eigenvector moment ﬂow that f λ,s undergoes. Theorem 1.12 ([15, Theorem 3.1]). Suppose that u is the solution of the symmetric Dyson vector ﬂow (1.17) and f (ξ) is given by (1.20) with the polynomials Q . Then it satisﬁes λ,s s the equation i,j 2ξ (1 + 2ξ ) f (ξ )− f (ξ) i j λ,s λ,s ∂ f (ξ) = . (1.21) s λ,s N (λ − λ ) i j i6=j N( ) i 2 2 N( ) i N1 N(ii1) 1 2 N1 N Figure 1: Example of the symmetric eigenvector moment ﬂow with a conﬁgu- ration of 7 particles. Now that we have the expression of the eigenvector moment ﬂow, we can give an heuristic for the apparition of a Cauchy proﬁle in the variance (1.3). Indeed, the single particle case m = 1 gives us the variance of an entry of an eigenvector. To understand this result, consider 8 the diagonal entries of the matrix D to be the quantiles of the semicircle law for instance, it is then interesting to consider the following continuous dynamics, deﬁne the operator K acting on smooth functions on [−2, 2] as f (x)− f (y) (Kf )(x) = dρ(y). (1.22) (x− y) −2 The diﬀerential equation ∂ f = Kf can be seen as a deterministic and continuous equivalent of (1.21) because of the rigidity property of the Dyson Browian motion eigenvalues. We then get the following lemma from [12] Lemma 1.13 ([12]). Let f be smooth with all derivatives uniformly bounded. For any x, y ∈ (−2, 2), denote x = 2 cos θ, y = 2 cos φ with θ, φ ∈ (0, π). Then −tK e f (x) = p (x, y)f (y)dρ(y) (1.23) where the kernel is given by −t 1− e p (x, y) := . (1.24) 2 2 i(θ+φ) −t/2 i(θ−φ) −t/2 e − e e − e Now at our small time-scale, we have p (x, y) . 2 2 (x− y) + t Hence (1.3) where m = 1 can be considered as a result of stochastic homogenization in a non-equilibrium setting when we consider the dynamics (1.21) in the bulk. For Theorem 1.7, we will study another observable which follows the same dynamics as in Theorem 1.12. This new observable has been analyzed in [16] to obtain universality for a class of band matrices. Deﬁne now the centered eigenvectors overlaps for symmetric matrices, p = u (α)u (α), i 6= j ∈ I, (1.25) ij i j α∈I p = u (α) − C , i ∈ I (1.26) ii i 0 α∈I where u are the eigenvectors of H and C is any constant in the sense that it does not s 0 depend on i but can depend on N. Now for ξ a conﬁguration of n particles on N sites, deﬁne the following set V = {(i, a), 1 6 i 6 N, 1 6 a 6 2η }. ξ i The set V will be a set of vertices. Consider now G the set of perfect matchings on V . For ξ ξ any edge on G, e = {(i, a), (j, b)}, deﬁne p(e) = p , P (G) = p(e) and ﬁnally ij e∈E (G) F (ξ) = E P (G) λ (1.27) λ,s M(ξ) G∈G where M(ξ) = (2ξ )!!, with (2m)!! being the number of perfect matchings of the com- i=1 plete graph on 2m vertices. Note that this quantity depend on the eigenvalues trajectories λ. i i 1 i2 3 i1 i i3 (a) A conﬁguration ξ with N (ξ) = 6 particles (b) An example of a perfect matching G ∈ G with 2 2 P (G) = p p p p i i i i i i i i 1 2 2 2 2 3 3 3 The previous quantity follows the same dynamics (1.21) as f (ξ). λ,s 9 Theorem 1.14 ([16, Theorem 2.6]). Suppose that u is the solution of the symmetric Dyson vector ﬂow (1.17) and F (ξ) is given by (1.27). Then it satisﬁes the equation λ,s i,j 1 2ξ (1 + 2ξ ) F (ξ )− F (ξ) i j λ,s λ,s ∂ F (ξ) = . (1.28) s λ,s N (λ − λ ) i j i6=j Third step: invariance of local statistics. The third and last step will be to obtain the result for the matrix W without any Gaussian part. We will do so using a variant of the dynamical method introduced in [15, Appendix A] which will show the continuity of resolvent statistics along the trajectory. This method was also used in [32] to study sparse matrices. We will see −1 that we need to take τ of order larger than N but smaller than t/N to use the continuity argument. In order to remove the Gaussian component for Theorem 1.7, we will use the reverse heat ﬂow and will need smoothness of the entries of our original matrix. Note that a moment matching scheme between the two matrix ensembles, which holds for any time −1 N τ 1, could be used to obtain the invariance of local statistics. However it can be of later interest to obtain the continuity estimate up to time t/N. The next section will state the local laws proved in diﬀerent papers ([35] and [14]). The third section is dedicated to prove Theorem 1.3, Corollary 1.5 and Theorem 1.7 for a short time relaxation of the matrix W . We use a maximum principle on f or F , a basic tool for t s s the analysis of parabolic equations. However, since we want a local result, remember that the variance depends on the position of the spectrum, we need to localize the maximum principle. We ﬁnish with an induction on the number of particles in the multi-particle random walk or a mutli-scale argument. For the third step, we will need a continuity result for the Dyson Brownian motion which will be shown in Subsection 4.1 and we will give the reverse heat ﬂow technique in Subsection 4.2. We will then conclude by combining the three steps in Section Acknowledgments. The author would like to kindly thank his advisor Paul Bourgade for many insightful and helpful discussions about this work and anonymous referees for many helpful comments on this article. 2. Local laws In this section, we focus on the diﬀerent local laws result for W . These local laws are high probability bounds, for simplicity we will now introduce the following notation for stochastic domination. For X = (X (u), N ∈ N, u ∈ U ) , Y = (Y (u), N ∈ N, u ∈ U ) . N N N N two families of nonnegative random variables depending on N (note that U can also depend on N), we will say that X is stochastically dominated by Y uniformly in ω, and write X ≺ Y , if for all τ > 0 and D > 0 we have τ −D sup P (X (u) > N Y (u)) 6 N N N u∈U for N large enough. If we have |X| ≺ Y for some family X, we will write X = O (Y ). Deﬁne now the resolvent of W and its normalized trace, N N X X |u ihu | 1 1 1 k k −1 G(z) = (W − z) = , ς (z) = Tr G(z) = (2.1) t t λ − z N N λ − z k k k=1 k=1 and denote G (z) the (i, j) entry of the resolvent matrix. In the rest of the section we will ij omit the dependence in t of the resolvent since we are not looking at its dynamics. 10 2.1. Anisotropic local law for deformed Wigner matrices An averaged local law was proved in [35]. The proof relies on Schur’s complement formula, large deviations bounds and interlacing formula in order to ﬁrst state a weak local law on the resolvent entries and the Stieltjes transform. The result then follows from a ﬂuctuation −1/2 −1 averaging lemma in order to go from the scale (Nη) to (Nη) . We ﬁrst give the deﬁnition of the limiting Stieltjes transform as the solution m (z) such that Im m (z) > 0 on the upper t t half plane of the following equation N N X X 1 1 1 m (z) = = g (t, z) (2.2) t α N D − z − tm (z) N α t α=1 α=1 where we deﬁned g (t, z) := D − t− m (z) α t We will also need the following lemma on the Stieltjes transform claiming that its imagi- nary part is of order one. ϑ,κ Lemma 2.1 ([35, Lemma 7.2]). Let ϑ, ω > 0 small constants and κ ∈ (0, 1). Take z ∈ D , for N large enough, the following bounds holds for t ∈ T , c 6 Im m (z) 6 C. Moreover ct 6 |D − z − tm (z)| 6 C. α t Note that the constants above do not depend on any parameter. Here is the averaged local law taken from [35]. Theorem 2.2 ( [35, Theorem 3.3]). Let W be as in Deﬁnition 1.1, ϑ > 0 and κ ∈ (0, 1), |ς (z)− m (z)| ≺ (2.3) t t Nη ϑ,κ uniformly in z = E + iη ∈ D The proof of Theorem 2.2 also gives the following entrywise local law from [35] also properly stated in [14]. Theorem 2.3 ([14, Theorem 2.4]). Let W be as in Deﬁnition 1.1 and ϑ > 0, κ ∈ (0, 1). ϑ,κ Uniformly in z = E + iη ∈ D , we have for the diagonal entries |G (z)− g (t, z)| ≺ |g (t, z)| , (2.4) αα α α Nη and for the oﬀ-diagonal entries |G (z)| ≺ min{|g (t, z)|,|g (t, z)|}. (2.5) αβ α β Nη In order to study hq, u i, we will need the following local law for hq, G(z)qi proved in [14], Theorem 2.4 ([14, Theorem 2.1]). Let ϑ > 0, κ ∈ (0, 1) and q a L -normalized vector of R , we have N N X X 2 2 hq, G(z)qi− q g (t, z) ≺ √ Im q g (t, z) (2.6) α α α α Nη α=1 α=1 ϑ,κ uniformly in z = E + iη ∈ D This theorem also gives us control of the resolvent as a bilinear form by polarization. We will give the proof of this corollary for completeness. 11 2 N Corollary 2.5. Let ϑ > 0, κ ∈ (0, 1), let v and w two L -normalized vectors of R , we have ! ! N N N X X X 2 2 hv, G(z)wi− v w g (t, z) ≺ √ Im v g (t, z) Im w g (t, z) (2.7) α α α α α α α Nη α=1 α=1 α=1 ϑ,κ uniformly in z = E + iη ∈ D Proof. Let μ ∈ R, a parameter ﬁxed later. Consider h(v + μw) , G (v + μw)i = hv, Gvi + μ hw, Gwi + 2μhv, Gwi, (2.8) by linearity and symmetry of the resolvent G. On one hand, using Theorem 2.4 on the ﬁrst two terms of the right hand side of (2.8), we get the equation N N X X 2 2 2 h(v + μw) , G (v + μw)i = 2μhv, Gwi + v g (t, z) + μ w g (t, z) α α α α α=1 α=1 !!! N N X X 2 2 2 +O √ Im v g (t, z) + μ w g (t, z) . (2.9) ≺ α α α α Nη α=1 α=1 On the other hand, using Theorem 2.4 on the left hand side of (2.8), we obtain h(v + μw) , G (v + μw)i = (v + μw ) g (t, z) α α α α=1 !! +O √ Im (v + μw ) g (t, z) . (2.10) ≺ α α α Nη α=1 Finally, combining (2.9) and (2.10) and choosing Im w v g (t, z) α α α α=1 μ = , Im w g (t, z) α=1 α we get the ﬁnal result. We will also need the following rigidity result from [35]. Theorem 2.6 ([35, Theorem 3.5]). Let ω > 0 be a small constant and κ ∈ (0, 1). For any t ∈ T , we have |λ − γ | ≺ k k,t uniformly in k ∈ A This control of the resolvent allows us to give an upper bound for the moments of the N−1 eigenvectors of W . Deﬁning, for a ﬁxed q ∈ S , " # Y 2ξ ( Nhq, u i) ϕ (ξ) = E λ (2.11) a(2η ) k=1 with u , . . . , u the eigenvectors of W and λ its eigenvalues, we have the following corollary. 1 N t Corollary 2.7. Let κ ∈ (0, 1), ϑ > 0 and ξ : A → N 2ξ 2 ϕ (ξ) ≺ σ (q, k, η) =: σ (q,ξ, η). (2.12) t t k=1 −1+ϑ −ϑ uniformly in N 6 η 6 N t. 12 −1+ϑ κ Proof. Fix any ϑ ∈ (0, 1) such that η = N t and k ∈ A , we have the following ﬁrst high probability bound with z = λ + iη k k √ 2 Nhq, u i η √ 2 k Nhq, u i = 6 N Imhq, G(z )qi k k 2 2 η (λ − Re(z )) + η k k ! !! N N 2 2 X X q 1 q α α = N Im +O Im . D − z − tm (z ) Nη D − z − tm (z ) α k t k α k t k α=1 α=1 We can then write, √ 2 X α 2 Nhq, u i ≺ NηIm ≺ σ (q, k, η). D − z − tm (z ) α k t k α=1 where we used the deﬁnition of ≺ and that ϑ is as small as we want and the smoothness of m (z). We ﬁnish the proof by deﬁnition of ϕ . t λ,t The local law gives us a strong control on both the eigenvalues and eigenvectors of our matrix ensemble. As W will undergo the Dyson Brownian motion, these quantities will still be controlled through a local law up to a small error coming from the time of the relaxation. We will now deﬁne the event of good eigenvalue paths (λ(t + s)) ) in the sense that all s∈(0,τ ) the estimates and bound from the previous section such as local laws from Theorems 2.2, 2.4 and 2.6 hold. First denote the resolvent of H and its normalized trace by 1 1 1 −1 G(s, z) = (H − z) and ς (s, z) = Tr H = . (2.13) s s N N λ (t + s)− z k=1 Deﬁnition 2.8. Let ε, ϑ > 0 and κ ∈ (0, 1). An eigenvalue conﬁguration λ is good is the following holds with overwhelming probability conditioning on λ(τ ) = λ for N large enough, ε −1 ϑ,κ 1. sup |ς (s, z)− m (z)| 6 N (Nη) uniformly in z ∈ D . t+s 06s6τ r P P N N 2 2ε −1/2 2 2. sup |hq, G(s, z)qi − q g (t + s, z)| 6 N (Nη) Im q g (t + s, z) α α 06s6τ α=1 α α=1 α ϑ,κ uniformly in z ∈ D . ε −1 κ 3. sup |λ (t + s)− γ | 6 N N uniformly in i ∈ A i i,t+s 06s6τ r Note that by the considerations in this section and a continuity argument, we see that good eigenvalue paths occur with overwhelming probability so that we can condition on having such a path in the following section. Note that we will make W undergo the Dyson Brownian motion for a small time τ so that the classical location γ and the deterministic counterpart i,t to resolvent entries g have small variations. This is the statement of the next lemma. Lemma 2.9. Let ε, ϑ > 0 and κ ∈ (0, 1). Conditionally on a good eigenvalue path as in Deﬁnition 2.8, we have the following control of the resolvent, N τ sup |ς (s, z)− m (z)| 6 + , Nη t 06s6τ 2ε 1 N τ sup G(s, z) − 6 √ + |g (t, z)| αα α D − z − tm (z) Nη t α t 06s6τ ϑ,κ uniformly in z ∈ D . For the rigidity results for eigenvalues we have sup |λ (t + s)− γ | 6 N + τ k k,t 06s6τ uniformly in k ∈ A . 13 ε Proof. The ﬁrst result comes from the fact that |∂ m (z)| 6 N /t which can be deduced from t t the time evolution of m from [35, Lemma 7.6] ∂ m (z) = ∂ (m (z)(m (z) + z)) t t z t t combined with the estimates |∂ m (z)| 6 C/t and |m (z)| 6 log N. The other error term z t t comes from the local law holding for a good eigenvalue path. The proof of the second bound comes from the following simple identity m (z) + t∂ m (z) N t t t ∂ g (t, z) = which gives |∂ g (t, z)| 6 |g (t, z)|. t α t α α (D − z − tm (z)) t i t using the fact that |m (z)| 6 C log N and |∂ m (z)| 6 N /t. For the rigidity estimate, we t t t combine the estimate |∂ γ | 6 C log N which can be found in [35, Lemma 7.6] as well with t i,t the rigidity coming from the good eigenvalue path. 3. Short time relaxation In this section, we are going to prove Theorems 1.3 and 1.7 for the Dyson Brownian motion starting from W using maximum principle. Note that in this section we will omit the subscript λ for simplicity. Recall the dynamics of the eigenvector moment ﬂow with n particles for H with H = W . τ 0 t i,j 1 2η (1 + 2η ) f (ξ )− f (ξ) i j τ τ ∂ f (ξ) = =: (B f )(ξ), τ τ τ τ N (λ (t + τ )− λ (t + τ )) i j (3.1) i6=j f (ξ) = ϕ (ξ) 0 t with ϕ (ξ) is deﬁned in (2.11). where we noted λ (t + τ ) the eigenvalues of H . Note that in t i τ the case of a single particle in k, we can write the dynamics 2 f (j)− f (k) τ τ ∂ f (k) = , τ τ N (λ (t + τ )− λ (t + τ )) (3.2) j k j=1 f (k) = ϕ (k). 0 t We cut the dynamics into two parts : the short range where most of the information will be and the long range. This decomposition in this context was ﬁrst introduced in [25]. Letting 1 ` N be a parameter that we will choose later, we then deﬁne i,j 2η (1 + 2η ) f (ξ )− f (ξ) i j τ (S (τ )f )(ξ) = , (3.3) N (λ − λ ) i j |j−k|6` i,j 1 2η (1 + 2η ) f (ξ )− f (ξ) i j τ τ (L (τ )f )(k) = . (3.4) N (λ − λ ) i j |j−k|>` Denote by U (s, τ ) the semigroup associated with S from time s to τ : ∂ U (s, τ ) = S (τ )U (s, τ ) (3.5) τ S S for any s 6 τ. We will denote in the same way U . It has been proved in [14, 15] that the parabolic short range dynamics has a ﬁnite speed of propagation in the following sense: deﬁne the following distance on the set of conﬁgurations with n particles d(η,ξ) = |x − y | (3.6) k k k=1 14 where (x , . . . , x ) are the positions of the particles in nondecreasing order of η and y of 1 n α ξ. The following lemma then states that if two conﬁgurations are far from each other, the short-range dynamics started at one and evaluated at the other is exponentially small with high probability. Lemma 3.1 ([14, Corollary 3.3]). Choose ` > Nτ, let ε > 0 be a small constant and κ ∈ (0, 1). Conditioning on a good eigenvalue path (λ(s)) , uniformly, for any function h on conﬁgurations of n particles and a conﬁguration ξ outside of the support of h in the sense that d(ξ,η) > N ` for any conﬁguration η inside the support of h, we have 0 n −cN sup U (s, s )h(ξ) 6 N khk e (3.7) S ∞ 06s6s 6t for any D > 0. In this section, we condition on an event occurring with overwhelming probability so that we can state the results deterministically. We state it here as the following lemma. Lemma 3.2. Let ω, a, b, ϑ and ε be small positive constants and D as in Deﬁnition 1.1. −1+a −a 1+b −b Let t ∈ T , and W as in (1.1). Let κ ∈ (0, 1), τ ∈ [N , N t] and ` ∈ [τN , N t]. ω t The dynamics (H ) induces a measure on the space of eigenvalues and eigenvectors s 06s6τ (λ(t + s), u(t + s)) . The event A of trajectories deﬁned by the following holds with 06s6τ overwhelming probability: ε −1 −1 ϑ,κ 1. sup |ς (s, z)− m (z)| 6 N (Nη) + τt uniformly in z ∈ D . 06s6τ r P P N N 2 2ε −1/2 −1 2 2. sup |hq, G(s, z)qi− q g (t, z)| 6 (N (Nη) + τt )Im q g (t, z) α α 06s6τ α α α=1 α=1 ϑ,κ uniformly in z ∈ D . ε −1 κ 3. sup |λ (t + s)− γ | 6 N (N + τ ) uniformly in i ∈ A . i i,t 06s6τ r 4. For any function h on conﬁgurations of n particles and a conﬁguration ξ supported outside of the support of h in the sense that d(ξ,η) 6 N ` for any conﬁguration η inside the support of h we have 0 n −cN sup U (s, s )h(ξ) 6 N khk e S ∞ 06s6s 6t The following lemma gives us a bound on the diﬀerence between the short-range and long-range dynamics basically stating that most of the information lies in the short-range dynamics. Lemma 3.3. Fix ` Nt and consider ξ to be a conﬁguration of n particles supported on A then for any eigenvalue paths (λ(t + s),u(t + s)) in A, we have for any ϑ > 0 and r 06s6τ −1+ϑ −ϑ N 6 η 6 N `/N, Nτ (n+4)ε 2 |(U (0, τ )− U (0, τ ))ϕ (ξ)| 6 N σ (q,ξ, η). (3.8) B S t −1+ϑ −ϑ Proof. Let η ∈ [N , N t]. Using Duhamel’s formula we can write |(U (0, τ )− U (0, τ ))ϕ (ξ)| = U (s, τ )L (s)f (ξ)ds B S t S s Now, by deﬁnition of the operator L (s) we have that j,k f (ξ )− f (ξ) s s L (s)f (ξ) = 2η (1 + 2η ) s j k N (λ − λ ) j k j,k,|j−k|>` so that we can bound, using Corollary 2.7 since ξ is supported on A , nε 2 k,j N σ (q,ξ, η) +|f (ξ )| |L (s)f (ξ)| 6 2n(1 + 2n) . N (λ − λ ) j k j,k:|j−k|>` 15 Now, say the conﬁguration is supported on p sites denoted (k , . . . , k ), then one can write 1 p X X nε 2 k ,j N σ (q,ξ, η) +|f (ξ )| |L (s)f (ξ)| 6 C . (3.9) s n N (λ − λ ) j k i=1 j:|j−k |>` Now, since ξ is supported on A , we have that by Corollary 2.7 and denoting ξ \ k the conﬁguration ξ where we removed a particule from the site k , k ,j (n−1)ε 2 f (ξ ) 6 N σ (q,ξ \ k , η)f (j). (3.10) s i s Consider now η = 2 `/N for q = [0,dlog (N/`)e], then we can bound dlog (N/`)e X X X 2 f (j) 1 hq, u iη s j q 2 + 2 N (λ − λ ) η (λ − λ ) η j i q k j p i q j:|j−k |>` j:|j−k |>` q=0 i i dlog (N/`)e 6 Imhq, G(s, λ + iη )qi. (3.11) k q 2 ` q=0 We can now use the anisotropic local law since λ lies in the spectral window and since we are on the event A, Imhq, G(s, λ + iη )qi 6 Imhq, G(s, λ + iη)qi 6 N Im m (λ + iη)σ (q, k , η) k q k t k t i i i i 2ε 2 6 N σ (q, k , η). (3.12) where we used the fact that η `/N 6 η and that m (z) is bounded in the spectral window. q t Combining the estimates (3.11) and (3.12) we obtain that f (j) N 2ε 2 6 N σ (q, k , η). N (λ − λ ) ` j k j:|j−k |>` Injecting this bound in (3.9) with (3.10), we obtain that (n+2)ε 2 |L (s)ϕ (ξ)| 6 N σ (q,ξ, η) where we used the fact that, by the same argument as in (3.11), 1 N . . N (λ − λ ) ` j k |j−k |>` Now, we can use that U is a contraction combined with the ﬁnite speed of propagation from Lemma 3.1 so that we can write that ε (n+3)ε 2 |U (s, τ )L (s)f (ξ)| 6 N sup |L (s)f (η)| 6 N sup σ (q,η, η). S s s ε ` ε η:d(η,ξ)6N ` η:d(η,ξ)6N ` (3.13) However, since we have that η is close to ξ we can use regularity of σ (q,η, η). Indeed, if one looks at the function q t ψ(x) = 2 2 (D − x− tRe m (x)) + (tIm m (x)) α t t α=1 Then we have that 2 2 q t(2(1 + t∂ Re m (x))(D − x− tRe m (x))− 2t ∂ Im m (x)Im m (x) x t α t x t t ∂ ψ(x) = 2 2 [(D − x− tRe m (x)) + (tIm m (x)) ] α t t α=1 So that we can obtain the bound using the fact that |∂ m (x)| 6 N /t, x t |∂ ψ(x)| 6 ψ(x). 16 2 We can use this bound in order to obtain the following variation formula for σ (q, k, η), d(ξ,η) 2 2 ε 2 σ (q,ξ, η)− σ (q,η, η) 6 N σ (q,ξ, η). (3.14) t t t Nt In (3.13), one can see that the supremum is only taken over conﬁgurations close to each other, ε 2 namely such that d(ξ,η) 6 N `. Since we have ` Nt, by (3.14), σ (q,ξ, η) varies slowly and we can ﬁnally bound (n+4)ε 2 |U (s, τ )L (s)f (ξ)| 6 N σ (q,ξ, η) S s The two previous lemmas will be very useful tools to prove Theorem 1.4. Indeed, the ﬁnite speed of propagation in Lemma 3.1 allows us to localize our problem but is a property of the short range dynamics, Lemma 3.3 then tells us that most of the information of the global dynamics is in this short range part. 3.1. Analysis of the moment observable To prove Theorem 1.3 we will prove the following intermediary proposition, Proposition 3.4. Conditionally on (λ,u) ∈ A, let κ ∈ (0, 1), ε > 0 and n be an integer. If N−1 κ q ∈ S , for any ξ : A → N such that N (ξ) = n, there exists a p depending on n such that we have 1/3 1 τ 2 −ε pε 2 −ε f (ξ) = σ (q,ξ, N τ ) +O N √ + σ (q,ξ, N τ ) . (3.15) t t Nτ where σ (q, k, τ ) is given by (1.3). The 1/3 exponent that we give here in the error is not optimal. We are not able to reach an optimal error because of the strong dichotomy we do between the short range and the long range dynamics and also the localization technique where more parameters need to be tuned. Using a multi-scale partition of the dynamics could improve the error term. Note also −ε the choice of the parameter η corresponds to N τ which optimize our error term. Let ε > 0 be a small constant. Recall that t ∈ T and τ t. First, the following lemma gives us a local law for f , in the case of a single particle, deduced from the isotropic local law for W in Theorem 2.4. ϑ,κ Lemma 3.5. For z ∈ D , we have ! ! N N X X 1 f (k) q τ α Im = Im N λ (t + τ )− z D − z − tm (z) k α t α=1 k=1 !! X 2 1 τ q ε α +O N √ + Im . Nη t D − z − tm (z) α t α=1 Proof. See ﬁrst that, by deﬁnition of f , we have 1 f (k) = E hq, G (z)qi λ N λ (t + τ )− z k=1 √ (d) H −1 where G := (H −z) is the resolvent of H . Now, the law of H is D + tW + τ GOE = τ τ τ 0 0 D + t + τW for some W a Wigner matrix. We can use Theorem 2.4 for this matrix and write H α Im hq, G qi = Im D − z − (t + τ )m (z) α t+τ α=1 !! X 2 1 q ε α +O N √ Im Nη D − z − (t + τ )m (z) α t+τ α=1 17 with m (z) the solution with positive imaginary part of the following self-consistent equa- t+τ tion, 1 1 m (z) = . t+τ N D − z − (t + τ )m (z) α t+τ α=1 ϑ,κ Note that we have, for z ∈ D , 0 < Im(m (z)) 6 C for some constant C so that, by a t+τ ϑ,κ Taylor expansion in τ t (remember that η t for z ∈ D ), we obtain ! ! N N 2 2 X X q q α α Im = Im D − z − (t + τ )m (z) D − z − tm (z) α t+τ α t+τ α=1 α=1 !! τ q ε α +O N Im . t D − z − tm (z) α t α=1 We get the ﬁnal result with the bound |∂ m (z)| 6 C (log N )/t. t t Let ξ ⊂ A be a ﬁxed conﬁguration, we want to use a maximum principle on a window centered around ξ of size w according to the distance (3.6). Since we make a small pertur- bation τ t, in order to notice the dynamics in this window, we need to have w Nτ. Furthermore, we want to look in the part of the spectrum where the eigenvector will be of typical size 1/ Nt, to localize the dynamics in this small part of the spectrum. We then need to take w Nt. We deﬁne the following ﬂattening and averaging operator, for a > 0 f (ξ) if d(ξ,ξ ) 6 a, Flat f (ξ) = (3.16) σ (q,ξ ) if d(ξ,ξ ) > a 0 w and (Av f )(ξ) = (Flat f )(ξ)da. (3.17) 0 ξ0 w/2 Notice that for every ξ, there exists a ∈ [0, 1] such that Av f (ξ) = a f (ξ) + (1− a )σ (q,ξ ). (3.18) ξ ξ ξ 0 w t 3.1.1 Proof of Proposition 3.4 : Case of a single particle To show the result (3.15) by induction, we will ﬁrst prove in the case of one particle with the dynamics (3.2). We ﬁrst prove under the hypotheses as in Proposition 3.4 in the case where n = 1, 1/3 1 τ 2 ε 2 f (k) = σ (q, k, η) +O N + σ (q, k, η) . (3.19) t t Nτ To do so, we want to use a localized maximum principle centered around k which is the position of the particle for the conﬁguration ξ in this case. However, we need to know that the maximum stays in that window, that is why we ﬁrst ﬂatten and average f (k) and use it as an initial condition for the dynamics (1.21). We will then make the short range dynamics work on f (k) during a time τ t. Since we use the short range dynamics for a time τ we will be able to use the ﬁnite speed of propagation (3.1) and we should choose Nt ` > Nτ for the range cut-oﬀ. We will take an explicit value at the end of the proof. The diﬀerent parameters and scaling is illustrated in Figure 3. Consider g , 1 g (j)− g (k) τ τ ∂ g (k) = with g (k) = (Av ϕ )(k). τ τ 0 k t N (λ (t + τ )− λ (t + τ )) j k |j−k|6` First note that, in order to prove (3.19), it is enough to show that ` Nτ 1 τ w Nη 2 2 g (k) = σ (q, k, η) +O + + √ + + + σ (q, k, η) (3.20) τ t t w ` Nη η Nt ` 18 Nt hq; u i N ` w Nt Figure 3: We see here a sketch of the variance proﬁle plotted in the spectral dimension: the projection q is ﬁxed and we plot the proﬁle as a function of the eigenvector u . We then localize the dynamics onto the small window plot- ted here: the window is small enough so that the eigenvector can be seen as “ﬂat” but large enough so that the short-range dynamics will not involve indices outside of this window. −1 where the parameter η, the spectral resolution, will be chosen so that N η t. Indeed we have, |f (k)− g (k)| = |(U (0, τ )ϕ )(k)− (U (0, τ )Av ϕ )(k)| (3.21) τ τ B t S k t = [(U (0, τ )− U (0, τ ))ϕ ] (k) + [U (0, τ )(Id− Av )ϕ ](k). (3.22) B S t S k t By Lemma 3.3, we get Nτ (n+4)ε 2 [(U (0, τ )− U (0, τ ))ϕ ](k) 6 N σ (q, k, η). (3.23) B S t Now, for the second term in (3.22). Since (ϕ − Av ϕ )(k) = 0 for k ∈ [[k −w/2, k +w/2]], t k t w w and taking w `N , looking at k ∈ [[k − w/3, k + w/3]] for instance, Lemma 3.1 tells us w w that the term is exponentially small. Thus, we obtain Nτ (n+4)ε 2 f (k) = g (k) +O N σ (q, k, η) . τ τ We ﬁrst prove the following equation which can be seen as an averaged version of the result. We now show the following lemma which is analogous to [15, Lemma 7.3] (k ) ϑ,κ Lemma 3.6. For k ∈ [[k − u, k + u]], set z = γ + iη ∈ D 0 w w k ,t 0 r 1 g (k) (k ) 2 Im − Im(m (z ))σ (q, k , η) t w (k ) N λ (t)− z k=1 ` Nτ 1 τ w √ (n+6)ε 2 6 N + + + + + η σ (q, k , η). (3.24) w ` Nη η Nt Proof. First, decompose the left hand side term into three diﬀerent terms : 1 (U (0, τ )Av ϕ )(k)− (Av U (0, τ )ϕ )(k) S k t k S t w w Im (3.25) (k ) N λ − z k=1 1 (Av U (0, τ )ϕ )(k)− Av f (k) k S t k τ w w + Im (3.26) (k ) N λ − z k=1 1 Av f (k) k τ w (k ) 2 + Im − Im(m (z ))σ (q, k , η) . (3.27) t w (k ) N λ − z k=1 To bound (3.25), we write (U (0, τ )Av ϕ )(k)− (Av U (0, τ )ϕ )(k) = U (ϕ )(k)da (3.28) S k t k S t t w w k w/2 19 with a a a U (ϕ )(k) = (U (0, τ )Flat ϕ )(k)− (Flat U (0, τ )ϕ )(k). t S t S t k k k w w w Look now at what happens around k − a, the other boundary of the window k + a can w w be bounded exactly the same way. By ﬁnite speed of propagation from Lemma 3.1, for k < k − a− `N , we easily get ε/2 a a n −N U (0, τ )Flat ϕ (k) = Flat (U (0, τ )ϕ ) (k) +O N e . S t S t k k w w The same equality is true for k > k − a + `N using the same argument. ε ε ∞ For k − a− `N 6 k 6 k − a + `N , since the operator U is bounded in ` , we have w w S a ε 2 U (ϕ )(k) 6 2 sup ϕ (j) 6 N σ (q, k , η) (3.29) t t w k t j:|j−k |6N ` where we used Corollary 2.7 combined with (3.14). Now, in the integrand in (3.28) there is a set of measure at most N ` which is not exponentially small which gives, combined with Theorem 2.2, that we can bound (3.25) by ` ` ε (k ) 2 ε 2 (3.25) 6 N Im m (z ))σ (q, k ) 6 N σ (q, k , η), (3.30) t w w t t w w (k ) where we used that Im(m (z )) is bounded in the spectral window from Lemma 2.1. To bound (3.26), noting that f = U (0, τ )ϕ , τ B t |(Av U (0, τ )ϕ )(k)− (Av U (0, τ )ϕ )(k)| 6 |[(U (0, τ )− U (0, τ )) ϕ ] (k)| k S t k B t S B t w w Nτ 5ε 2 6 N σ (q, k, η) where we applied Lemma 3.3. Thus we have CNτ 1 1 Nτ 5ε 2 5ε 2 (3.26) 6 N σ (q, k, η)Im 6 N σ (q, k, η) (3.31) t t (k ) ` N λ − z ` k=1 where we used the averaged local law from Theorem 2.2 and the fact that in the spectral (k ) window we have Im m (z ) 6 C. (k ) To bound (3.27), we want to use the local laws. Recalling that z = γ + iη, then X X 1 (Av U (0, τ )ϕ )(k) 1 (Av U (0, τ )ϕ )(k) k B t k B t w w Im = Im (k ) (k ) 0 0 N λ − z N λ − z k k k=1 |k−k |6N η +O ησ (q, k , η) . where we used the fact that similarly to the proof of Lemma 3.3, we have that, for any threshold ` , 1 (Av f )(k) Nη k τ ε 2 Im = O N σ (q, k , η) (3.32) (k ) N ` λ − z k 0 |k−k |>` 0 0 20 √ Taking ` = N η gives here the bound. If we use the notation (3.18), we obtain 1 (Av f )(k) k τ Im = (3.33) (k ) N λ − z k=1 1 a f (k) + (1− a )σ (q, k , η) √ k τ k w t 2 = Im +O ησ (q, k , η) (3.34) (k ) √ λ − z |k−k |6N η ! ! N N X X 1 f (k) 1 1 = a Im + (1− a )σ (q, k , η)Im (3.35) k k w 0 0 t (k ) (k ) 0 0 N λ − z N λ − z k k k=1 k=1 1 (a − a )f (k) + (a − a )σ (q, k , η) √ k k τ k k w 0 0 t + Im +O ησ (q, k , η) . (k ) N λ − z √ k |k−k |6N η (3.36) Note now that q η α (k ) 2 Im = Im(m (z )) + σ (q, k , η), t 0 (k ) (k ) 0 0 D − z − tm (z ) t α t α=1 so that, by Lemma 3.5 and Theorem 2.2, we obtain (k ) 2 2 (k ) 0 0 (3.35) = a Im(m (z ))σ (q, k , η) + (1− a )σ (q, k , η)Im(m (z )) k t 0 k w t 0 t 0 t 1 1 τ √ ε 2 +O N + + + η σ (q, k , η) Nη t Nη Note that we do not keep the error η/t as we take η τ so that τ/t is of larger order. Now using the deterministic bound (3.14), we obtain that since k ∈ [[k − w, k + w]], w w 1 τ w √ (k ) 2 ε 2 (3.35) = Im(m (z ))σ (q, k , η) +O N √ + + + η σ (q, k , η) t w w t t t Nt Nη C|i−k| Finally, with the elementary property |a −a | 6 , we get that (3.36) 6C ησ (q, k , η). i k w Putting these estimates together, we obtain 1 τ w √ ε 2 (3.27) = O N √ + + + η σ (q, k , η) (3.37) Nη t Nt Combining (3.30),(3.31) and (3.37), we get the ﬁnal result 1 g (k) (k ) 2 Im − Im(m (z ))σ (q, k , η) t w (k ) N λ − z k=1 ` Nτ 1 τ w √ ε 2 = O N + + √ + + + η σ (q, k , η) . w ` Nη t Nt Now, we just need to prove that (3.20). Let k be the index such that g (k ) = m τ m max g (k) and z = λ + iη. If we have that k τ k 2 −10 g (k )− σ (q, k , η) 6 N , (3.38) τ m m −10 there is nothing to prove. Now if the left hand side is greater than N , by ﬁnite speed of propagation, k is in the interval [[k − u, k + u]]. Indeed, if it is not then the diﬀerence in m w w 21 (3.38) would be exponentially small. We then have 1 g (j)− g (k ) τ τ m ∂ g (k ) = τ τ m N (λ − λ ) j k |j−k |6` j6=k X X 1 ηg (j) g (k ) η τ τ m 6 − 2 2 2 2 Nη (λ − λ ) + η Nη (λ − λ ) + η j k j k m m |j−k |6` |j−k |6` m m j6=k j6=k m m 1 1 g (k) g (k ) N Nη 1 τ τ m 6 Im − Im(m (z)) +O + σ (q, k , η) t m η N λ − z η η ` Nη k=1 6 σ (q, k , η)− g (k ) (3.39) w τ m N Nη 1 ` Nτ 1 τ w √ 2 2 +O + σ (q, k , η) + + + √ + + + η σ (q, k , η) m w t t η ` Nη w ` Nη t Nt (3.40) where we used in the ﬁrst inequality that g (k ) is the maximum, in the second inequality τ m 1+ε 2 that extending the sum to all j adds an error (N η/`σ (q, k , η)) using (3.32) and Theorem 2.2 combined with the estimate from Lemma 2.1. Finally in the last inequality we used (3.24), c 6 Im(m (z)) 6 C in the spectral window and that the rigidity errors that appears from changing λ into γ from Theorem 2.6 are smaller than the other terms. Injecting our k k ,t m m variance σ (q, k , η) in (3.39) which does not depend of τ, 2 2 ∂ g (k )− σ (q, k , η) 6 − g (k )− σ (q, k , η) τ τ m w τ m w t t Nη ` Nτ 1 τ w √ N +O + + + + + + η σ (q, k , η) . ` w ` t Nt η Nη Thus, writing S = g (k )− σ (q, k ) we get τ τ m t w C Nη ` Nτ 1 τ w √ N ∂ S 6 − S +O + + + √ + + + η σ (q, k , η) . τ τ τ w η ` w ` Nη t Nt η Note that S is not necessary diﬀerentiable as the maximum is not necesarily unique for instance, but one can get the result by instead considering S − S s τ ∂ S = lim sup . τ τ s− τ s→τ Using Gronwall’s lemma, we have if η τ, Nη ` Nτ 1 τ w √ ε 2 −C S = O N + + + √ + + + η σ (q, k , η) + N τ w ` w ` Nη t Nt for any C. We can do the same reasoning with the inﬁmum. Finally taking the following parameters 1/3 −ε 2 η = N τ, w = Nτ (Nt) , ` = Nτu (3.41) we get the result for a single particle. Note that with these choices of parameters we have −1 the correct relations: N η τ ` w Nt. 3.1.2 Proof of Proposition 3.4: Case of n particles In the previous part of the proof, we looked only at the second moment E Nhq, u (t)i |λ which corresponds to a single particle in the site k. Now, we will do the proof of (3.15) by induction on the number of particles. We can ﬁrst deﬁne the same objects as the single particle case: we will consider the short range dynamics for a small time τ t with initial condition an average of the eigenvectors moment of W localized onto a speciﬁc window. More precisely deﬁne, with ξ being the t 0 22 conﬁguration with n particles that lies at the center of our window of size w in the sense of the distance (3.6), ij 1 g (ξ )− g (ξ) τ τ ∂ g (ξ) = , (3.42) τ τ N (λ − λ ) i j |i−j|6` i6=j g (ξ) = (Av f ) (ξ). (3.43) 0 ξ t By the same reasoning as for the one particle case, using Lemmas 3.3 and 3.1 with n particles, we get CNτ (n+4)ε 2 |f (ξ)− g (ξ)| 6 N σ (q,ξ, η). (3.44) τ τ To reason by induction on the number of particles, we need to show the following equation, κ (k ) similar to (3.24) in the case of one particle. For k ∈ A , deﬁne z = γ + iη and let ξ a r k r r conﬁguration of n particles with at least one particle in k , we need to show k ,k X r 1 g ξ 2 2 Im − a f (ξ \ k )σ (q, k , η) + (1− a )σ (q,ξ , η) ξ τ r t r ξ 0 (k ) N λ (t + τ )− z k=1 ` Nτ 1 τ (n+4)ε 2 =O N + + √ + σ (q,ξ , η) . w ` Nη t (3.45) where ξ \ k denote the conﬁguration where we removed one particle in k from ξ. r r We apply the same decomposition in three terms as in the single particle case, the ﬁrst two terms can be bounded the same way and we can bound the left hand side of (3.45) by k ,k (Av f ) ξ ξ τ Im (k ) N λ − z k=1 (k ) 2 (k ) 2 r r − a f (ξ \ k )Im(m (z ))σ (q, k , η) + (1− a )Im(m (z ))σ (q,ξ , η) ξ t r t r ξ t 0 t t ` Nτ (n+4)ε 2 +O N + σ (q,ξ , η) (3.46) u ` Now we need to see that, k ,k k ,k 2 r r Av f (ξ ) = a k ,kf (ξ ) + (1− a k ,k )σ (q,ξ , η) ξ τ r τ r 0 0 ξ ξ t k ,k 2 = a f (ξ ) + (1− a )σ (q,ξ , η) ξ τ ξ 0 k ,k 2 + (a k ,k − a )f (ξ ) + (a − a k ,k )σ (q,ξ , η) r ξ τ ξ r 0 ξ ξ t k ,k 2 d(ξ ,ξ)σ (q,ξ , η) k ,k 2 t = a f (ξ ) + (1− a )σ (q,ξ , η) +O . ξ τ ξ 0 We can use the same decomposition into |k − k | 6 N η and the averaged local law from Theorem 2.2 to get ! ! N N k ,k r k ,k X X r 1 Av f ξ 1 f (ξ ) ξ τ w τ Im = a Im (k ) (k ) r r N λ − z N λ − z k k k=1 k=1 (k ) 2 2 + (1− a )Im(m (z ))σ (q,ξ , η) +O σ (q,ξ , η) . (3.47) ξ t 0 0 t t Nη Consider the sum in the right hand side, recall that there is at least one particle in k and denote k , . . . , k with m 6 n, the sites where there is at least one particle in the conﬁguration 1 m (k ) ξ. Recall that z = γ + iη, k ,t N m k ,k k ,k nε X r X r X 1 ηf (ξ ) 1 ηf (ξ ) N τ τ 2 k ,k r i = +O σ (q,ξ , η) 2 2 2 2 N (γ − λ ) + η N (γ − λ ) + η Nη k ,t k k ,t k r r k=1 i=1 k∈{ / k ,...k } 1 m (3.48) 23 where we used Corollary 2.7 on the indices we removed from the sum. Now we have the following equality for the ﬁrst sum by deﬁnition of f , k ,k X r 1 ηf (ξ ) 2 2 N (γ − λ ) + η k ,t k k∈{ / k ,...k } 1 m 2j 0 2(j −1) r 2 Y X ηz z z r j = E λ (3.49) 2 2 a(2j ) a(2(j − 1)) N ((γ − λ ) + η ) r r k ,t j j∈{ / k ,...,k } 16r 6m 1 m r 6=r By Lemma 3.5, we have, N m 2 ε X X X ηz 1 1 ηz N k 2 = +O σ (q, k , η) , 2 2 2 2 N (γ − λ ) + η N (γ − λ ) + η Nη k ,t j k ,t k r r j∈{ / k ,...,k } k=1 j=1 1 m (3.50) N m 2 ε X X 1 z N k 2 = Im +O σ (q, k , η) , (3.51) (k ) N λ − z Nη j=1 k=1 1 τ N (k ) 2 ε 2 2 = Im(m (z ))σ (q, k ) +O N + σ (q, k ) + σ (q, k , η) . t r r j t t t t Nη Nη j=1 (3.52) Combining (3.52) and (3.49) and using the bounds on the variations of σ (3.14), we get k ,k X r 1 ηf (ξ ) (k ) 2 = Im(m (z ))σ (q, k , η) f (ξ \ k ) t t r τ r 2 2 N (γ − λ ) + η k ,t k k=1 1 τ ε 2 +O N √ + σ (q,ξ, η) (3.53) Nη Finally, combining (3.53) and (3.47), we have k ,k 1 Av f ξ ξ τ (k ) 2 Im = a f (ξ \ k )Im(m (z ))σ (q, k , η) ξ t r t t r (k ) N λ − z k=1 1 τ (k ) 2 ε 2 + (1− a )Im(m (z ))σ (q,ξ , η) +O N √ + σ (q,ξ, η) ξ t 0 t t Nη which, combined with (3.46), gives us (3.45). We now follow the same proof as in the case of one particle : we state a maximum principle on the ﬂattened and averaged moment. First deﬁne ξ = max g (ξ), (3.54) m τ N (ξ)=n and let k , . . . , k be the positions of the particles of the conﬁguration ξ with m 6 n. We 1 m m 24 are going to use our induction hypothesis in the maximum principle inequalities by (3.45). i,j C g (ξ )− g (ξ ) τ τ m ∂ g (ξ ) 6 τ τ m N (λ − λ ) i j |i−j|6l i6=j k ,j X X r X C 1 ηg (ξ ) g (ξ ) η τ τ m 6 − 2 2 2 2 N η (λ − λ ) + η η (λ − λ ) + η j k j k r r r=1 |j−k |6l |j−k |6l r r j6=k j6=k r r m N k ,j 2 X X r C 1 g (ξ ) Nη σ (q,ξ , η) τ m m (k ) t 6 Im − g (ξ )Im s (z ) +O τ m τ (k ) η N ` η λ − z r=1 j=1 (k ) 2 (k ) 2 r r 6 a f (ξ \ k )Im(m (z ))σ (q, k , η) + (1− a )Im(m (z ))σ (q,ξ , η) ξ τ m r t r ξ t 0 m t m t r=1 (3.55) (n+4)ε C N 1 τ ` Nτ Nη (k ) 2 − g (ξ )Im(m (z )) +O √ + + + + σ (q,ξ , η) . τ m t m η η t w ` ` Nη r=1 Now, we use the induction assumption on f (ξ \ k ) which is a (n− 1)th moment and τ m r obtain 1/3 1 τ 2 2 2 f (ξ \ k )σ (q, k , η) = σ (q,ξ , η) +O √ + σ (q,ξ , η) . (3.56) τ m r t r m m t t Nτ Besides, we can easily see that, since d(ξ ,ξ ) 6 2w, from (3.14), w m N w 2 2 2 2 2 σ (q,ξ , η) − σ (q,ξ , η) 6 σ (q,ξ , η). (3.57) 0 m m t t t Nt Now, injecting (3.56) and (3.57) in (3.55), we get 2 2 ∂ g (ξ )− σ (q,ξ , η) 6 − g (ξ )− σ (q,ξ , η) τ τ m 0 τ m 0 t t r=1 (n+6)ε 1/3 Nτ ` 1 τ Nη w 1 τ N √ √ +O + + + + + + + σ (q,ξ , η) . ` w t ` Nt t η Nη Nτ Doing the same reasoning as in the proof for one particle, we get, by applying Gronwall’s lemma, g (ξ ) = σ (q,ξ , η) τ m m 1/3 Nτ ` 1 τ Nη w 1 τ +O + + √ + + + + √ + σ (q,ξ , η) ` w Nη t ` Nt t Nτ We can again do the same reasoning with the inﬁmum and choosing the parameters as in (3.41) the claim from Proposition 3.4 follows. 3.2. Analysis of the perfect matching observable In this subsection we will again condition on the event A of good eigenvalue paths where the local laws and the ﬁnite speed of propagation holds. Consider now a deterministic set of indices I ⊂ [[1, N ]]. Note that in the deﬁnition of the centered overlaps p , we can only center ii by a constant not depending on i. However in Theorem 1.7, one can see that the expectation of the probability mass of the ith-eigenvector on I clearly depends on i. Thus, we will need to localize our perfect matching observables onto a window of size w chosen later and show that these p are, up to an error depending on w, centered around the same constant. The ii size of the window w will be taken so that Nτ w Nt similarly to the previous section. More precisely, we will ﬁx an integer i ∈ A and consider the set of indices A (i ) = {i ∈ [[1, N ]], γ ∈ [γ − (1− κ)w, γ + (1− κ)w]} 0 t,i t,i t,i w 0 0 25 so that we will take for our centered diagonal overlaps X X 2 2 p = u (α) − C with C = σ (e , i ), (3.58) ii i 0 0 α 0 α∈I α∈I the overlaps for i 6= j will not change. First consider these overlaps for the matrix W and deﬁne, similarly to the previous subsection Φ (ξ) = E P (G) λ . M(ξ) G∈G This quantity corresponds the perfect matching observable for our initial matrix W and we make it undergo the dynamics (1.28) so that we deﬁne F (ξ) = U (0, s)Φ (ξ) s B t where B is deﬁned in (3.1). We now prove the result from Theorem 1.7 for a Gaussian divisible ensemble for p . We will need the following technical lemma allowing us to bound i i 0 0 the p by the perfect matching observables. ij Lemma 3.7 ([16, Lemma 3.6]). Take an even integer n, there exists C > 0 depending on n such that for any i < j and any time s we have n (1) (2) (3) E [p (s) |λ] 6 C F (ξ ) + F (ξ ) + F (ξ ) (3.59) ij s s s (1) (2) where ξ is the conﬁguration of n particles in the site i, ξ n particles in the site j, and (3) ξ an equal number of particles between the site i and j. We will also use repeatedly the following bound on the eigenvectors which comes from Corollary 2.7 2 ε u (α) 6 N I. k∈I The purpose of this section is to prove Theorem 1.7 for another matrix ensemble: a deformed Wigner matrix perturbed by a small Gaussian component. More precisely, we state it as the following theorem. Theorem 3.8. Consider a and ω two small positive constants and κ ∈ (0, 1), take t ∈ T , D a deterministic diagonal matrix given by Deﬁnition 1.1 and W a Wigner matrix given by 0 −1+a −a Deﬁnition 1.2. Let τ ∈ T := [N , N t], then if u , . . . , u are the eigenvectors of the a 1 N matrix H = D + tW + τ GOE, deﬁne the error I τ Ξ(τ ) = + I , Nτ we have, for any k, ` ∈ A with k 6= ` and any ε > 0 and D > 0, X X 2 2 ε −D P u (α) − σ (α, k, τ ) + u (α)u (α) > N Ξ(τ ) 6 N . k k ` α∈I α∈I As one can see in the statement of Theorem 3.8, we will need the small time τ corre- sponding to the size of the Gaussian perturbation to be of order smaller than t so that the eigenvalues and eigenvectors barely changed during that time. We will later choose a speciﬁc τ and optimize all our diﬀerent parameters when using the reverse heat ﬂow technique to remove this small Gaussian component. One of these parameters will be a cut-oﬀ for the dy- namics as in Subsection 3.1. Indeed, in order to use a maximum principle on the dynamics, we will split it in the same way: a short-range dynamics with generator S that will contain most of the information and a long-range part with generator L we need to control as in [25] where S and L are deﬁned respectively in (3.3) and (3.4). Lemma 3.1 will help us localize the dynamics onto a small set of conﬁgurations. Now the following lemma says that most of 26 the information of the dynamics is given by the short-range, bounding the diﬀerence between B and S . It is analogous to Lemma 3.5 in [16]. First deﬁne (u,v) S = sup F (ξ). ξ⊂I,u6s6v ε κ Lemma 3.9. For any intervals J ⊂ A (i ) and J = {i, d(i, J ) 6 N `} ⊂ A (i ) since in w 0 out in 0 we will take ` Nw, any conﬁguration ξ such that N (ξ) = n supported on J and any in −1 N τ w we have n−1 n−2 Nτ w n I n (0,τ ) (0,τ ) (0,τ ) |((U (0, τ )− U (0, τ ))Φ (ξ)| 6 N S + I S + S B S t J J J out out out ` t ` (3.60) where F is the perfect matching observable deﬁned in 1.27. This bound is used in this form so we can obtain information on a box in space by extracting information from a larger box. Iterating this bound will give us Theorem 3.8. Proof. We will follow the proof from [16]. Deﬁne the following ﬂattening operator. For f a function on conﬁgurations of n particles and ξ such a conﬁguration, f (ξ) if ξ ⊂ {i, d(i, J ) 6 a} , in (Flat f )(ξ) = (3.61) 0 otherwise. We make the functions vanish outside of a certain interval. We use now Duhamel’s formula and write ((U (0, τ )− U (0, τ )Φ ) (ξ) = U (s, τ )L (s)F (ξ)ds. S B t S s Now, see that, by deﬁnition of the ﬂattening operator and the fact that ξ is supported on J , in d(Supp (L (s)F − Flat (L (s)F )) ,ξ) > N `. s N ` s With this bound, we can use the ﬁnite speed of propagation from Lemma 3.1 and obtain, using that U is a contraction in ` , −cN /2 |U (s, τ )L (s)F (ξ)| 6 max (Flat (L (s)F )) (ξ) +O e . S s N ` s ξ⊂J out ˜ ˜ Thus we need to control |(L (s)F )(ξ)| for ξ = {(i , j ), . . . (i , j )} a conﬁguration of n > m s 1 1 m m particles supported in J . out We have, by deﬁnition of L , i ,k X X ˜ p X X F (ξ ) 1 ˜ ˜ L (s)F (ξ) 6 C + C F (ξ) . s n n s 2 2 N (λ − λ ) N (λ − λ ) i k i k p p 16p6m 16p6m |i −k|>` |i −k|>` p p For the second term in the previous inequality, we can use rigidity estimates from Theorem 2.6 and a dyadic decomposition and see that 1 N 6 N , (3.62) N (λ − λ ) ` i k k,|i −k|>` so that we have the bound i ,k X X ˜ p F (ξ ) N s (0,τ ) ε ε L (s)F (ξ) 6 N + N S . (3.63) 2 out N (λ − λ ) ` i k 16p6m k,|i −k|>` For the ﬁrst sum in (3.63), we will ﬁrst restrict it to the sites k such that there are no particles i ,k in the conﬁguration ξ, so that we will have η ˜ = 1. Note that, if we have η ˜ 6= 0 then by 27 i ,k ˜ ˜ p deﬁnition of ξ supported on J , ξ is also supported on J . This gives us the bound out out i ,k i ,k X ˜ p X ˜ p X F (ξ ) F (ξ ) 1 s s (0,τ ) 6 + S . 2 out 2 N (λ − λ ) N (λ 2 N (λ − λ ) i k i −λ ) i k p p k p k,|i −k|>` k,|i −k|>` k,|i −k|>` p p p η ˜ =0 η ˜ 6=0 k k X i ,k ˜ p F (ξ ) N s (0,τ ) 6 + CN S 2 2 out N (λ − λ ) ` i k k,|i −k|>` η ˜ =0 where in the second inequality we used the rigidity of the eigenvalues, which gives us (λ − 2 2 λ ) > C (`/N ) for |i −k| > `, and the fact that there is at most m sites k such that η ˜ 6= 0. k p k By deﬁnition of the perfect matching observables, we can write i ,k X ˜ p X X p(e) F (ξ ) s e∈E (G) = C (n)E λ . 2 2 N (λ − λ ) N (λ − λ ) i k i k p p G∈G k,|i −k|>` k,|i −k|>` i ,k p p p η ˜ =0 η ˜ =0 k k In order to control this term, we will consider two diﬀerent types of perfect matchings. Deﬁne the following partition of G into two subsets (1) G = {G ∈ G , {(k, 1), (k, 2)} ∈ E (G)} , (3.64) (2) G = {G ∈ G , {(k, 1), (k, 2)} ∈/ E (G)} . (3.65) i k i i i k i i 1 2 3 1 2 3 (1) (2) (a) A perfect matching from G (b) A perfect matching from G i ,k i ,k 2 2 ξ ξ (1) We will begin by bounding the contribution from (3.64). Note ﬁrst that, for G ∈ G , we have p(e) = p × Q ((p(e)) ) (3.66) kk 1 e∈E (G) e∈E (G) with Q ((p(e)) ) = p(e). 1 e∈E (G) e∈E (G)\{(k,1),(k,2)} See also that Q is a monic monomial of degree n − 1 so that we can use Lemma 3.7 and obtain n−1 (0,τ ) E Q ((p(e)) ) λ 6 C sup |F (ξ)| 6 C S . (3.67) 1 s e∈E (G) out 06s6τ G∈G ξ\i ξ⊂J ,N (ξ)=n−1 out Combining (3.66) and (3.67), we now only need to bound X X p p N kk kk = +O . 2 2 2 N (λ − λ ) N (λ − λ ) ` i k i k p p k,|k−i |>l, k,|k−i |>` p p η ˜ =0 In order to bound the sum from the right hand side of the previous equation, ﬁrst deﬁne the −ε following functions, for |z − λ | 6 N `/N, kk f (z) = , N (z − λ ) − + k, γ ∈/[E ,E ] t,k 1 2 kk g(z) = N (z − λ ) − + − + k, γ ∈/[E ,E ]∪[E ,E ] t,k 1 1 2 2 28 − + − ε ε ε where E = γ , E = γ , E = γ , E = γ , E = γ , 1 t,i −` t,i −`−N t,i −`+N 2 t,i +` t,i +`−N p 1 p 1 p p 2 p E = γ . Let also Γ be the rectangle with vertices E ± i`/N and E ± i`/N. We t,i +`+N 1 2 2 p therefore want to bound, up to a N term, kk = ∂ f (z) N (λ − λ ) z=λ i k i p p k,|k−i |>` −ε ` Now consider C the circle centered in λ with radius N , then by Cauchy’s formula, we i i p p can write, 1 f (z) ∂ f (λ ) = dξ. z i 2iπ (z − λ ) C p By using another Cauchy integral formula on the contour Γ for f and seeing that for λ , z int inside the contour and λ outside of the contour we have, by a residue calculus, ext Z Z dξ dξ = = 0 (ξ − λ )(ξ − λ ) int i (ξ − λ )(ξ − λ ) Γ p Γ ext i . Note that we also have Z Z 1 p 1 p kk kk = O 2π (ξ − λ ) (ξ − λ ) 2π (ξ − λ ) (ξ − λ ) i k γ p γ i k so that we can write Z Z 1 g(ξ) Img(ξ) ∂ f (λ ) = dξ = O dξ . (3.68) z i 2 2 2π (ξ − λ ) (ξ − λ ) i i Γ p Γ p We will ﬁrst control the part of the contour closest to the real axis. Consider Γ = {z = E + iη ∈ Γ, |η| < N /N}, as in [16] and bounding p by I, we obtain kk Img(ξ) IN dξ = O N . 2 2 (ξ − λ ) ` Γ p Now for the rest of the contour, note that we can add the missing eigenvalues to the total sum in g up to adding an error of order N I/`. Finally we just have to bound N p kk Im |dξ| ` N (ξ − λ ) Γ\Γ k=1 X X N 1 Imm (ξ) 1 = Im G (ξ)− Img (t, z ) |dξ| (3.69) αα α 0 ` N Imm (z ) N Γ\Γ t 0 α∈I α∈I where we used the deﬁnition of p and of σ and deﬁned z = γ + iη , with η t is the kk t 0 t,i 0 0 center of our window of size w t with positive imaginary part. Now, using the entrywise local law from Theorem 2.3 and expanding between z and ξ since |ξ − z | 6 w, we have 0 0 X X 1 Imm (ξ) 1 Im G (ξ)− Img (t, z ) αα α 0 N Imm (z ) N t 0 α∈I α∈I X X 1 1 Imm (ξ) 6 Im (G (ξ)− g (t, ξ)) + Img (t, ξ)− Img (t, z ) αα α α α 0 N N Imm (z ) t 0 α∈I α∈I I w ε ε 6 N p + N I . N|Imξ| Injecting this bound in the contour integral, we get the following bound. w I ε ε (3.69) 6 N I + N √ . 29 Finally, putting all the contributions together and coming back to (3.68), we obtain N I w |∂ f (λ )| 6 N + I . (3.70) z i ` ` t (2) Consider now the contribution from (3.65), ﬁrst see that for G ∈ G , there exists q and q in {1, . . . m}, such that p(e) = p p × Q ((p(e)) ) ki ki 2 e∈E (G) q q 1 2 e∈E (G) with Q a monic monomials of degree n− 2. Then using Lemma 3.7, we can bound " # ! (n−2) (0,τ ) E Q ((p(e)) ) λ = O S . 2 e∈E (G) out G∈G Besides, we can bound the term with the cross-edges in the following way X X X p p N N N i k i q q 1 2 2 2 2 2 6 (p + p ) = u (α) + u (α) 6 I . q q i k i k 1 2 q q 2 2 2 2 1 2 N (λ − λ ) ` ` ` i k k,|k−i |>` k=1 α∈I η ˜ =0 Putting everything together, we get n−1 n−2 Nτ w n I n (0,τ ) (0,τ ) (0,τ ) |((U (0, τ )− U (0, τ ))Φ (ξ)| 6 N S + I S + S B S t J J J out out out ` t ` which is exactly the result wanted. We will ﬁrst prove the following proposition in order to deduce Theorem 1.7. Proposition 3.10. For any ε and N large enough the following holds. For any intervals κ −ε J ⊂ A and J ⊂ {i, d(i, J ) 6 N Nw} we have in out in n−1 ` Nτ 1 I w (0,τ ) (0,τ ) (0,τ ) ε ε S 6 N + + S + N √ + I (S ) J J J in out out Nw ` Nτ t Nτ n−2 I Nτ (0,τ + N + I (S ) . (3.71) 2 out Nτ ` Proof. We will use the short-range dynamics and its ﬁnite speed of propagation property in order to localize the maximum principle in J . We will then use the local laws and Lemma in 3.9 in order to get a Gronwall type bound. Deﬁne the following averaging operator. For f a function on conﬁgurations and ξ a conﬁguration, ε X 3N Av(f ) = Flat (f ). (3.72) Nw 1 Nw 2 Nw <a< ε ε 3 3 N N Note that, if ξ is not included in J , by deﬁnition of the ﬂattening operator, (Av(f ))(ξ) = 0. out The purpose of these operators is to change the initial condition in order to remove the particles far from the initial interval J . Note also that we can write, for any ξ, in Av(f )(ξ) = a f (ξ), with a ∈ [0, 1]. Note that we have the elementary bound |a − a | 6 CN /Nw. Deﬁne now ξ ξ ξ the following dynamics ∂ Γ = S (s)Γ , 0 6 s 6 τ s s s (3.73) Γ (ξ) = (Av Φ )(ξ). 0 t 30 Now, if one takes a conﬁguration ξ supported on J , it suﬃces to show the bound in in Proposition 3.10 for Γ. Indeed |F (ξ)− Γ (ξ)| 6 |(U (0, τ )Φ − U (0, τ )Φ ) (ξ)| +|U (0, τ )(Φ − AvΦ )(ξ)| τ τ B t S t S t t n−1 n−2 CNτ w n I n (0,τ ) (0,τ ) (0,τ ) ε ε 6 N S + I S + S + exp(−cN ) J J J out out out ` t ` where we bounded the ﬁrst term by using Lemma 3.9 and the second term using the ﬁnite speed of propagation. Indeed, since ξ is supported on J , (Id − Av)Φ vanishes for any in t conﬁguration supported on J . Note that we can use Lemma 3.9 since we will take `/N w. out In the rest of the proof, we will prove the bound from Proposition 3.10 for Γ . If we already have for some C > 0, −C Γ (η ) := sup Γ (ξ) 6 N τ m τ ξ, N (ξ)=n then we have nothing to prove by the argument above and the deﬁnition of F . However, if −C this supremum is greater than N , then by the ﬁnite speed of propagation of S , we know 3Nw that η will be supported in, for instance, i, d(i, J ) 6 . m in ε 4N Consider now, a parameter η that we will choose later and denote also m 6 n the number of sites with at least a particle and j , . . . , j those sites. Then, we can write 1 m k,j 2η (1 + 2η ) Γ η − Γ (η ) X m,k m,j τ m τ m ∂ Γ (η ) = (3.74) τ τ m N (λ (t + τ )− λ (t + τ )) k j 0<|j−k|6` k,j η(Γ η − Γ (η )) X τ m τ m 6 (3.75) 2 2 Nη (λ − λ ) + η k j 16p6m k, 0<|j −k|6` k,j Γ η X m X C 1 1 6 − Im − Γ (η ) Im (3.76) τ m Nη λ − z Nη z − λ k j j k p p 16p6m 16p6m k, 0<|j −k|6` k, 0<|j −k|<` p p with z = λ +iη. For the second term, see that for p ∈ [[1, m]], if we choose η to be smaller j j p p than `/N, #{k, 0 < |j − k| 6 `} > C` > CNη > C # k, |λ − λ | 6 η > C Nη p j k where we used, in the two last inequalities, the rigidity of the eigenvalues. Now we can write X X 1 η Im > > CN. 2 2 z − λ (λ − λ ) + η j k j k p p 16p6m 16p6m k, 0<|j −k|6` k, 0<|λ −λ |6η p j k We now need to control the ﬁrst term in (3.76). To do so, we will split it in the three following terms: j ,k j ,k p p (U (0, τ )AvΦ ) (η )− (AvU (0, τ )Φ ) (η ) S t m S t m Im (3.77) N (z − λ ) j k k, 0<|j −k|6` j ,k j ,k p p (AvU (0, τ )Φ ) (η )− (AvU (0, τ )Φ ) (η ) m m S t B t +Im (3.78) N (z − λ ) j k k, 0<|j −k|6` j ,k (AvU (0, τ )Φ ) (η ) B t m +Im (3.79) N (z − λ ) j k k, 0<|j −k|6` To bound the ﬁrst term, we will use the ﬁnite speed of propagation property of S from Lemma 3.1. Indeed, we can write 3N j ,k j ,k j ,k p p p (U (0, τ )AvΦ ) (η )− (AvU (0, τ )Φ ) (η ) = U η S t m S t m a m Nw 1 Nw 2 Nw <a< ε ε 3 N 3 N 31 with U = U (0, τ )Flat Φ − Flat U (0, τ )Φ . a S a t a S t j ,k Fix a and consider three cases, if η is supported on {i, d(i, J ) > a + N `} then by In deﬁnition of Flat we have Flat U (0, τ )Φ = 0, a S t and by Lemma 3.1 we have cN |U (0, τ )Flat Φ | 6 exp − . S a t j ,k Now if η is supported on {i, d(i, J ) 6 a− N `} then again by deﬁnition of Flat , m In a j ,k j ,k p p Flat (U (0, τ )Φ ) η = (U (0, τ )Φ ) η . a S t m S t m Thus cN j ,k j ,k p p U η 6 U (0, τ ) (Φ − Flat Φ ) η 6 exp − . a m S t a t m j ,k Finally, if η is supported on {i, d(i, J ) 6 a + `N }, ﬁrst note that there can only be In 2n`N such a, then one can see that we can use the ﬁnite speed of propagation if we remove particle away from a at distance 2`N for instance, then j ,k U η 6 |Flat U (0, τ )Φ | +|Flat U (0, τ )Flat ε Φ | a a S t a S a+2`N t +|Flat U (0, τ ) (Φ − Flat )| a S t a+2`N cN 6 kFlat Φ k +kFlat k + exp − a t a+2`N ∞ ∞ cN (0,τ ) 6 2S + exp − , out where we used that U is a contraction in kk . Finally we can bound (3.77), S ∞ ` cN (0,τ ) (3.77) 6 N S + exp − (3.80) out Nw 2 where we used the fact that 1 1 Im 6 ς τ, z 6 N (3.81) N z − λ j k k,0<|j −k|<` where the Stieltjes transform ς is deﬁned in (2.13). For (3.78), we will use Lemma 3.9. Indeed, ﬁrst note that in the short-range regime, the set of k such that |j − k| 6 ` is included in 2κ A . Then we can bound n−1 n−2 Nτ w n I n (0,τ ) (0,τ ) (0,τ ) (3.78) 6 N S + I S + S (3.82) J J J out out out ` t ` where we used the fact that Av is a contraction and (3.81). Finally, in order to bound the third term (3.79), we will use the local law for Φ . First write j ,k j ,k p p (AvU (0, τ )Φ ) η a j ,kF η X m X p m B t τ Im = Im N (z − λ ) N (z − λ ) j k j k p p k, 0<|j −k|6` k, 0<|j −k|6` p p j ,k j ,k p p a F η a j ,k − a F η X X p η τ m η τ m m m = Im + Im N (z − λ ) N (z − λ ) j k j k p p k, 0<|j −k|6` k, 0<|j −k|6` p p j ,k a F η X η τ m (0,τ ) = Im +O N 1 a − a j ,k S . (3.83) |j −k|6` η p m J η out N (z − λ ) j k k, 0<|j −k|6` 32 But we have the bound, for |j − k| 6 `, j ,k N d η ,η ε m m N ` a − a j ,k 6 6 . η p Nw Nw In order to bound the ﬁrst term, see ﬁrst that we can remove the contributions of k ∈ {j , . . . , j } writing 1 p j ,k j ,k p p a F η a F η ε X η τ m X η τ m m m (0,τ ) Im = Im +O S . out N (z − λ ) N (z − λ ) Nη j k j k p p k, 0<|j −k|6` k, 0<|j −k|6` p p k∈{ / j ,...,j } 1 p Recall now the deﬁnition of Φ from (1.27) and write, j ,k F η X τ m X X p(e) e∈E (G) = E λ . (3.84) j ,k N (z − λ ) N (z − λ ) j k j k p p M η m G∈G k, 0<|j −k|6` k, 0<|j −k|6` j ,k p p p k∈{ / j ,...,j } k∈{ / j ,...,j } 1 p 1 p First consider the contribution of (3.64) in the sum in (3.84), denote e = {(k, 1), (k, 2)} and write X X p(e) X X Y e∈E (G) kk = p(e) N (z − λ ) N (z − λ ) j k j k p p (1) (1) k, 0<|j −k|6` k, 0<|j −k|6` e∈E (G)\{e } p p k G∈G G∈G j ,k j ,k p p k∈{ / j ,...,j } η k∈{ / j ,...,j } η 1 p 1 p m m X Y X kk = p(e) . N (z − λ ) j k G∈G e∈E (G) k, 0<|j −k|6` η \j p m p k∈{ / j ,...,j } 1 p (3.85) To control the last term in (3.85), we can use the local law. First write, X X X b p p p I kk kk kk = + +O N , (3.86) N (z − λ ) N (z − λ ) N (z − λ ) Nη j k j k j k p p p k=1 k, 0<|j −k|6` k,|j −k|>` p p k∈{ / j ,...,j } 1 p where we used the bound |p | 6 N I. Now, recall the deﬁnition of p from (3.58) so that kk kk we have X X X Imm (z ) p 1 t j 1 kk p Im = Im G (z )− Img (t, z ) αα j α 0 N (z − λ ) N Imm (z ) N j k t 0 k=1 α∈I α∈I X X 1 1 Imm (z ) t j 6 Im G (z )− g (t, z ) + Img (t, z )− Img (t, z ) αα j α j α j α 0 p p p N N Imm (z ) t 0 α∈I α∈I I w 6 N √ + I Nη t where z := γ + iη. Note that we used the fact that η is supported in J , so that 0 t,i m out |z − z | = |γ − λ | 6 N w. We can then use Lemma 3.7 and bound 0 j t,i j p 0 p X Y n−1 (0,τ ) p(e) = O sup |Φ (ξ)| = O (S ) out ξ⊂J ,N (ξ)=n−1 out G∈G e∈E (G η \j m p Now, consider the contribution of (3.65) in the sum from (3.84). Note that for any graph (2) G in G , there exists q and q in {1, . . . , m}, a ∈ {1 . . . , η } and b ∈ {1, . . . , η }, such that i i 0 ξ q 33 0 e := {(k, 1), (q, a)} and e := {(k, 2), (q , b)} are edges in G. We can then write q q X X p(e) X X Y p p i ,k i ,k e∈E (G) q 0 = p(e) , N (z − λ ) N (z − λ ) j k j k p p (2) (2) k, 0<|j −k|6` k, 0<|j −k|6` e∈E (G)\{e ,e } p p q 0 G∈G G∈G j ,k j ,k p p η η m m X X Y X p p i ,k i 0,k q = p(e) , N (z − λ ) j k 0 0 p q,q =1 q,q e∈E (G) k, 0<|j −k|6` G∈H η \j m p (3.87) q,q where we deﬁned the set of graphs H to be the set of perfect matching of the complete graph on the set of vertices V where we removed a single particle at the site i and i . Note ξ q q q,q that for any graph G ∈ H , p(e) is a monomial of degree n− 2. e∈E (G) η \j m p Now, we can bound the imaginary part of the second sum in (3.87), X X b p p C N I i ,k i 0,k q q 2 2 Im 6 p + p = O . (3.88) i ,k i 0,k q N (z − λ ) Nη Nη j k k=1 k, 0<|j −k|6` k∈{ / j ,...,j } 1 p For the last inequality, we used the following identity on eigenvectors X X 2 2 p = u (α) i,k k=1 α∈I and that for any ε > 0,using the entrywise local law from Theorem 2.3 on a diagonal entry of the resolvent, 2 −1+ε −1+ε u (α) 6 N Im G(τ, λ + iN ) 6 k k α,α Nt Again, we can bound the other term from (3.87) using Lemma 3.7, n−2 X Y (0,τ ) p(e) = O N sup |Φ (ξ)| = O S out ξ⊂J ,N (ξ)=n−2 0 out q,q e∈E (G) G∈H η \j m p Finally, putting all these estimates together, we get the Gronwall-type inequality, 1 N ` Nτ 1 (0,τ ) ∂ Γ (η ) 6 − Γ (η ) +O + + S τ τ m τ m out η η Nw ` Nη ! ! !! b b n−1 n−2 I w I Nτ (0,τ ) (0,τ n n + + I (S ) + + (S ) (3.89) J J out 2 out t Nη ` Nη In order to get a proper bound using Gronwall’s lemma, we need to take η τ but to get the best estimates possible, we also have to take η as large as possible. Hence, considering −ε η = N τ we have the bound ` Nτ 1 I w n−1 (0,τ ) (0,τ ) (0,τ ) ε 2ε S 6 N + + S + N + I (S ) J J J in out out Nw ` Nτ t Nτ n−2 I Nτ (0,τ + N + I (S ) (3.90) 2 out Nτ ` which gives the Proposition 3.10. Now that we have the bound from Proposition 3.10, we are able to get a bound on the p using Lemma 3.7. To do so, we can use a sequence of set of indices with decreasing size ij and apply recursively Proposition 3.10. We will also need to choose the right parameters `, w and τ. 34 −1+δ Proof of Theorem 3.8. Consider ﬁrst any ε small enough, such that if we write t = N −1 (recall that t ∈ T so that t N ) we have ε < 4δ/3. and a large D > 0. Then we can take the following parameters : w = N τ and ` = N τw, (3.91) −1 note that we have the right bounds between these parameters: N τ `/N w t, and deﬁne the following sequence of sets of indices, deﬁned implicitely, J = A (i ), 0 0 −ε J = {i : d(i, J ) 6 N Nw}. i i+1 From Proposition 3.10 we have the following bound holding with overwhelming probability, I N τ n−1 1 n−2 (0,τ ) (0,τ ) (0,τ ) (0,τ ) −ε/2 b b n n S 6 N S + √ + I (S ) + I (S ) . J J J J i+1 i i i t Nτ Nτ Now see that as long as we have (0,τ ) 1/n 3ε/2 (S ) > CN Ξ(τ ) with Ξ given in (1.12), we obtain the recursive bound (0,τ ) (0,τ ) −ε/2 S 6 N S . J J i+1 i −1 But if we take a very large i so that the previous bound cannot hold, for instance i = d3ε e, then it means that for such a i we have the bound (0,τ ) 1/n 3ε/2 (S ) 6 CN Ξ. Now using the deﬁnition of p , we have for i ∈ A (i ), ii 0 1 w 2 2 u (α) − σ (α, i) 6 |p | + I 6 |p | + Ξ(τ ). (3.92) i ii ii N t α∈I Finally, using Markov’s inequality, taking for instance n = b3D/εc, and using Lemma 3.7 to (0,τ ) bound the p by S we have ij ε −D P (|p | +|p | > N Ξ(τ )) 6 N . (3.93) ii ij The result then follows from combining (3.92) and (3.93). 4. Approximation by a Gaussian divisible ensemble 4.1. Continuity of the Dyson Brownian motion In Subsection 3.1 we showed that the moments of the eigenvectors of the matrix H are asymptotically those of a Gaussian random variable with variance σ . Now, since τ is a small −1 time, recall that N τ t, we can use the continuity of the Dyson Brownian motion to show that H and H = W have the same local statistics. In order to state a proper τ 0 t continuity lemma we need to have a dynamics with constant second moments and vanishing expectation. First see that the variance of the centered model is h i E (W − D ) = . t,ij ij Consider, for 0 6 s 6 τ, the following variance-preserving dynamics on symmetric matrices. dB 1 ˜ ˜ d H (s)− D = − H (s)− D ds, (4.1) 2t H (0) = W = D + tW. (4.2) The following lemma gives us a continuity argument between H (τ ) and W . It is similar to Lemma A.1 in [15] or Lemma 4.3 in [32]. We will later use this lemma on the resolvent entries. 35 Lemma 4.1. Denote ∂ = ∂ . Take F a smooth function of the matrix entries satisfying αβ ˜ αβ X N| H (s)− D | αβ ˜ ˜ E sup +| H (s)− D | ∂ F (θH ) 6 M (4.3) αβ αβ N t θ, 06s6τ α6β where (θH ) = θ H with θ = 1 for {k, l} =6 {α, β} and θ ∈ [0, 1]. Then αβ αβ αβ kl αβ ˜ ˜ E[F (H (τ ))]− E[F (H (0))] = O (τ ) M. Proof. By Itô’s formula we have 1 N ˜ ˜ ˜ ˜ ∂ E[F (H (s))] = − E[ H (s)− D ∂ F (H )]− E[∂ F (H )]. s αβ s s αβ αβ 2N t α6β Using Taylor expansions, we can write, forgetting the dependence in time for clarity, ˜ ˜ ˜ ˜ E[ H − D ∂ F (H )] = E[ H − D ∂ F ] + E[ H − D ∂ F ] αβ αβ ˜ ˜ αβ H =D H =D αβ αβ αβ αβ αβ αβ αβ ˜ ˜ +O E sup H − D ∂ F (θH ) . αβ αβ t 3 2 3 ˜ ˜ = ∂ F +O E sup H − D ∂ F (θH ) . αβ H =D αβ αβ αβ αβ and 2 2 3 ˜ ˜ ˜ E[∂ F (H )] = E[∂ F ] +O E sup H − D ∂ F (θH ) . s ˜ s αβ αβ αβ H =D αβ αβ αβ Putting everything together the claim follows. This continuity property of the Dyson Brownian motion gives us a control over the eigen- ˜ ˜ values and eigenvectors of H (0) and H . 2m Corollary 4.2. Let κ ∈ (0, 1) and m ∈ N. Let Θ : R → R be a smooth function satisfying (k) −C sup |Θ (x)(1 +|x|) < ∞, 2m k∈[[0,5]],x∈R for some C > 0. Denote u ˜ (s), . . . , u ˜ (s) the eigenvectors of H (s) associated with the eigen- 1 N ˜ ˜ values λ (s), . . . , λ (s). Deﬁne, for a small a the time domain 1 N " # N t 0 −a T = , N , N N then for any τ ∈ T , there exists p > 0 depending on Θ, a, κ and r such that ˜ ˜ H H 2 −p s 0 ˜ sup E − E Θ N (λ − γ ), hq, u ˜ i 6 N . k k,t k κ σ (q, k, η) I⊂I ,|I|=m,kqk=1 t k∈I (4.4) Proof. We prove this corollary using the continuity estimate from Lemma 4.1. To do so, we use the techniques introduced in [33] in order to change estimates of the resolvent below −1 microscopic scales, in other words control G(E + iη) for η N , into estimates on eigen- vectors. Indeed, it has been shown that such a control of the resolvent combined with an estimate on the number of eigenvalues in a very small interval allows us via integrating the resolvent over such an interval to gain estimates on eigenvectors. We can then split the proof into two results we need to show: (i) A level repulsion estimate on the eigenvalues for both matrix ensembles of the following form: for E ∈ I and a small ξ > 0 there exists d > 0 such that −1−ξ −1−ξ −ξ−d P λ ∈ [E − N , E + N ] > 2 6 N . 36 (ii) Comparison of the resolvent below microscopic scales: for any smooth function F of −1−ξ polynomial growth, there exists a c > 0 and a ξ > 0 such that for all N < η < t, H W −c τ t sup E − E F hq, G(z )qi 6 N kqk =1, Im q g (t, z ) α k κ α=1 E ,...,E ∈I 1 m k=1 (4.5) where z = E + iη. k k We ﬁrst prove (i) for the eigenvalues of W . This property can be deduced from gap universality, note that gap universality for Gaussian perturbation of size t ∈ T has been shown in [35] (a stronger level repulsion estimate can also be found in [35, Section 5]). In [35, Subsection 2.4], Landon−Yau explains that they can deduce universality for deformed Wigner ensembles. However, they state the result for an initial condition such that r t. It has been conﬁrmed by the authors that it is a simple typographical error and should be read as r t. We will nonetheless give an idea of the proof of (i) for the sake of completeness. As said earlier, we will ﬁrst apply Lemma 4.1 to −1 κ −1−ξ ˜ ˜ F (H ) = Tr(H − z) for z in z = E + iη, E ∈ I , N < η < t s s for ξ > 0 arbitrarily small. See that in Lemma 4.1, we need to bound a functional of the form F (θH ) for θ a perturbation of two entries of the matrix. Since we will only need bounds such as Theorem 2.4 or Corollary 2.5 with a possible N for any small positive ε window, such bounds still hold for the perturbated matrix (see [26, Section 15]). So we will explain ˜ ˜ the bound for the third derivative of F applied directly to H . Note that by deﬁnition of H, we have |(H − D) | ≺ so that we can bound the left hand side of 4.3 by αβ 1 t |∂ F (H )|. (4.6) αβ N N α6β −1 Taking the third derivative of F with respect to an entry, we obtain, writing G = (H − z) for simplicity XX ∂ F (H ) = − G G G G s γα β α β α β γ αβ 1 1 2 2 3 3 γ=1 α,β where {α , β } = {α, β} for ` = 1, 2, 3. To bound the sum in the previous equation, we will ` ` need the following high probability bound for the oﬀ-diagonal entries of the resolvent which follows directly from Theorem 2.3 |G (z)| ≺ √ |g (t, z)g (t, z)|. (4.7) αβ α β Nη −1 Note that these bounds holds for η N , we will ﬁrst consider such η. Finally we can bound (4.6), r r 2ξ X X XX 1 t N t 1/2 ∂ F (H ) ≺ |g | (|g g g g g g |) . (4.8) s γ α β α β α β αβ 1 1 2 2 3 3 N N N N α6β α6β γ=1 α,β Now, from Lemma 7.5 of [35], we have |g (t, z)| 6 C log N (4.9) γ=1 where the constant C only depend on D our diagonal matrix. Besides, in the last product of (4.8), there are, by deﬁnition of α and β, three occurrences of g and three occurrences of 37 g . Thus, r r 2 2ξ 2ξ X X CN log N t CN log N t 3/2 (4.8) 6 |g (t, z)g (t, z)| 6 |g (t, z)| α β α N N Nt N α6β α=1 (4.10) 2ξ 3 6 N log N (4.11) −1 where we used (4.9) in the ﬁrst inequality, the fact that |g | 6 Ct in the second and (4.9) again in the ﬁnal inequality. In order to go below microscopic scales, recall that we used local laws that holds down to mesoscopic scales, we can use the following identity, for y 6 η, 1 η 1 Im TrG(E + iy) 6 Im TrG(E + iη) . N y N Finally, using Lemma 4.1, we get, 1 N H W 5ξ −e τ t sup E − E TrG(z) 6 N τ 6 N (4.12) N t E∈I for some e > 0 by taking τ ∈ T for a > 5ξ. We can easily generalize this result to a product of trace, indeed taking F (H ) = F with F = TrG(z ), s k k k k=1 we can take the third derivative and write m m X Y X X Y 3 3 2 ∂ F = ∂ F F + 3 ∂ F ∂ F F k k αβ k k αβ αβ αβ 1 2 k =1 k6=k k =1 k 6=k k6=k ,k 1 1 1 2 1 1 2 X X X Y + ∂ F ∂ F ∂ F F . (4.13) αβ k αβ k αβ k k 1 2 3 k =1 k 6=k k 6=k ,k k6=k ,k ,k 1 2 1 3 1 2 1 2 3 Then, using the ﬁrst and second derivative of F , X X ∂ F = G G , αβ k kα βk k=1 α, β {α,β}={i,j} X X ∂ F = G G G , k kα β α β k αβ 1 1 2 2 k=1 α,β {α ,β }={i,j} k k we can bound (4.13) in a similar way and ﬁnishing the bound by Lemma 4.1. Again, now that we have any polynomial of ﬁxed degree, we can also extend to any smooth function F with polynomial growth. ˜ ˜ Now, a consequence of these uniform bounds in Re(z) between H = W and H for 0 t τ τ ∈ T for some small a gives us a comparison of the gap distribution between these two matrix ensembles (see [27] for instance). Namely, there exists c > 0 such that for any O a smooth test function of n variables and any index i such that γ ∈ I , we have for N large i,t enough and i , . . . , i indices such that i 6 N , 1 n k h i (N ) (N ) W H −c t τ 1 E − E O Nρ (γ )(λ − λ ), . . . , Nρ (γ )(λ − λ ) 6 N . i,t i i,i+i i,t i i+i t 1 t n (4.14) But H is a matrix with a small Gaussian component following the conditions of [35], so that we have, for this matrix ensemble, gap universality. Hence, combining this gap universality with the continuity of the Green’s function, we obtain gap universality of the 38 √ matrix ensemble D + tW . In other words, there exists c > 0 such that, taking the same assumptions as for (4.14), we can write, h i (N ) (N ) E O Nρ (γ )(λ − λ ), . . . , Nρ (γ )(λ − λ ) i,t i i,i+i i,t i i+i t 1 t n h i GOE (N ) (N ) −c −E O Nρ (μ )(λ − λ ), . . . , Nρ (μ )(λ − λ ) 6 N (4.15) i i i,i+i i i i+i sc 1 sc n where ρ is the density of Wigner’s semicircular law and μ its quantiles deﬁned by sc i p i ρ (x) = 4− x 1 , dρ (E) = . (4.16) sc [−2,2] sc −∞ Finally, (4.15) combined with Theorem 2.6 gives us the level repulsion estimate (i) for the matrix W , indeed consider E ∈ I , and ` the index such that |γ − E| 6 min |γ − E| , `,t k,t k∈I then, for any ε ˜ > 0, we have −1−ξ −1−ξ W −1−ξ P i, λ ∈ [E − N , E + N ] > 2 6 P |λ − λ | < N i k k+1 ε˜ |k−`|6N GOE −1−ξ −c +ε ˜ 6 P (|λ − λ | < N + N k k+1 ε˜ |k−`|6N −2ξ+ε ˜ −c +ε ˜ 6 N + N −ξ−δ 6 N for some δ > 0 by taking ε ˜ and ξ > 0 small enough. Note that we used rigidity in the ﬁrst inequality, gap universality in the second and a level repulsion estimate for GOE matrix which for instance can be found in [25] . In order to get the resolvent estimate (ii), we will use Lemma 4.1. To do so, we will ﬁrst explain how to get the bound M for −1 ˜ ˜ F (H ) = hq, (H − z) qi Im q g (t, z) γ=1 for z ∈ C down to below microscopic scales. To get the right bound, we will ﬁrst need to use −1+ξ local laws which hold down to mesoscopic scales η = N . Now for the third derivative of F , ﬁrst write X X |∂ F (H )| = q G G G G q (4.17) s a aα β α β α β b b αβ P 1 1 2 2 3 3 Im q g (t, z) γ 16a,b6N α,β γ=1 γ where {α , β } = {i, j} for k = 1, 2, 3. In order to bound the four terms coming up in the k k previous equation we will need Corollary 2.5. Writing (2.6) for v = q and w = e , we obtain v u X t 2 hq, Ge i = q g (t, z) +O √ Im q g (t, z) Im (g (t, z)) . α α α ≺ γ α Nη γ=1 Note that since we want a bound holding down to microscopic scales, the error terms has −1−ξ to be taken into account. In particular, we will consider η & N so that we can bound −1/2 ξ/2 every (Nη) by N . In the following computations, we will not bound the errors coming cross terms for simplicity, they can be bounded in a similar way. 39 We can divide the sum in (4.17) in three parts. The ﬁrst case consists in {β , α } = 1 2 {β , α } = {α, β}. In this case, note that, necessarily, {α , β } = {α, β} and write 2 3 1 3 q G G G G q = hq, Ge iG G he , Gqi a aα β α β α β b b α β α β α β 1 1 2 2 3 3 1 1 2 2 3 3 16a,b6N ≺ min (|g (t, z)|, |g (t, z)|) hq, Ge ihe , Gqi α β α β 1 3 Nη 2ξ ≺ N min (|g (t, z)|, |g (t, z)|) |q g (t, z)q g (t, z)| (4.18) α β α α β β ! ! + min (|g (t, z)|, |g (t, z)|) Im q g (t, z) |g (t, z)g (t, z)| . α β γ α β γ=1 (4.19) Putting the leading order (4.18) in the sum of (4.6), we have the bound t 1 (4.18) (4.20) N 2 N Im q g (t, z) γ 16α6β6N γ=1 γ 2ξ t N 6 |q q ||g (t, z)g (t, z)| α β α β N Im q g (t, z) γ 16i6j6N γ=1 γ (4.21) 2ξ t CN 2 2 2 6 q + q |g (t, z)g (t, z)| α β P α β N Im q g (t, z) γ 16α6β6N γ=1 γ (4.22) N N 2ξ X X t CN 2 2 2 6 q |g (t, z)| |g (t, z)| . α β P α N Im q g (t, z) γ α=1 β=1 γ=1 γ (4.23) −1 Note that by deﬁnition of g (t, z) = (D − z − tm (z)) , the fact that Imm (z) 1 and α α t t η 6 t, we can write |g (t, z)| Im (g (t, z)) . (4.24) α α Besides we also have from (4.9), N N X X C CN |g (t, z)| 6 |g (t, z)| 6 log N. (4.25) α α t t α=1 α=1 Injecting now (4.24) and (4.25) in (4.20), we get the bound r r 2ξ t 1 t N N (4.18) 6 C log N (4.26) N 2 N N t N Im q g (t, z) γ 16α<β6N γ=1 γ 3ξ N N 6 . (4.27) Nt t Looking now at the error term (4.19) and injecting it in the sum (4.6), we obtain r r 2ξ X X t 1 t N 3/2 (4.19) 6 |g (t, z)g (t, z)| α β N N N N Im q g (t, z) γ 16α<β6N 16α<β6N γ=1 r 2 2ξ t N 6 |g (t, z)| N Nt α=1 3ξ 6 N . (4.28) 40 The second case are the terms where one term is diagonal and the other is an oﬀ-diagonal term. More precisely the set α and β such that β = α and β 6= α or β 6= α and β = α . 1 2 2 3 1 2 2 3 Note that necessarily, in that case, α = β . For instance consider the term 1 3 hq, Ge iG G he , Gqi = hq, Ge iG G he , Gqi. (4.29) α β α β α β i jj αβ i 1 1 2 2 3 3 Putting all the leading terms from Theorem 2.3, (4.7) and (2.7), we obtain the bound 2 2 hq, Ge iG G he , Gqi ≺ |q | |g (t, z)| |g (t, z)|√ min(|g (t, z)|,|g (t, z)|) (4.30) α ββ αβ α α α β α β Nη + Im q g (t, z) |g (t, z)g (t, z)||g (t, z)| min(|g (t, z)|,|g (t, z)|) γ α β β α β 3/2 (Nη) γ=1 (4.31) ! ! 2ξ 2 2 2 3/2 6 N q |g (t, z)g (t, z)| + Im q g (t, z) |g (t, z)g (t, z)| . α β γ α β α γ γ=1 (4.32) Then injecting the bounds (4.24) and (4.25) in the sum of (4.6), one gets 2ξ 3ξ t N CN N 2 2 q |g g | 6 (4.33) α β P α N Nt t N Im q g (t, z) γ 16α<β6N γ=1 and for the second term, r r 2ξ t N N 3/2 3ξ |g (t, z)g (t, z)| 6 N . (4.34) α β N N t 16α<β6N The ﬁnal case consists of α and β such that {{β , α },{β , α }} = {{α, α},{β, β}}. Note 1 2 2 3 that, in this case, we necessarily have α 6= β . For instance, consider the term 1 3 hq, Ge iG G he , Gqi = hq, Ge iG G he , Gqi. (4.35) α β α β α β α ββ αα β 1 1 2 2 3 3 Again, taking the leading terms from the local laws from Threom 2.3 and Corollary 2.5, 2 ξ 3/2 2 hq, Ge iG G he , Gqi ≺ |q q ||g (t, z)g (t, z)| + N |g (t, z)g (t, z)| Im q g (t, z) . α ββ αα β α β α β α β γ γ=1 (4.36) Then using similar bounds as the ﬁrst case one gets r r r t 1 1 N N 2ξ (4.36) ≺ N + . (4.37) N 2 Nt t t N Im q g (t, z) γ 16α<β6N γ=1 γ Finally, putting together (4.27), (4.28), (4.33), (4.34) and (4.37), we get the bound, for −1+ξ η = N . r r t 1 N 3 3ξ ∂ F (H ) ≺ N . (4.38) αβ N N t 16α<β6N In order to get a bound for η below microscopic scales, we can use the following inequality, for any y 6 η, which can be found in [27, Section 8], |hv, G(E + iy)wi| 6 C log N Imhv, G(E + iη)wi. This bounds allows us to get below microscopic scales for F and its derivatives since they only κ −1−ξ involves such quantity as hv, G(E + iy)wi. Thus, uniformly in E ∈ I and N 6 η 6 t, we have 5ξ M = O N . (4.39) −5ξ t Using now Lemma 4.1, we can make W undergo the dynamics H up to a time τ N t s with ξ arbitrarily small in order to get the right bound. 41 For a product of resolvent entries, one can do similar computations and bounds. Indeed consider m > 0, and ˜ ˜ ˜ F (H ) = F (H ) with F (H ) = hq, G(z )qi, s k s k s k k=1 then one can write the third derivative of F as (4.13) and using the fact that ∂ F = − hq, Ge ihe , Gqi, (4.40) αβ k α β {α,β}={i,j} ∂ F = hq, Ge iG he , Gqi (4.41) k α β ,α β αβ 1 1 2 2 α,β where {α , β } = {i, j} and using the same type of bounds as for (4.38), we obtain the result i i (4.5) since the extension to any smooth function with polynomial growth is also clear. 4.2. Reverse heat ﬂow In Subsection 3.2, we showed Theorem 3.8, which corresponds to our main result for the −1 matrix H = W + τ GOE for N τ t with a general Wigner matrix W in the τ t deﬁnition of W . Thus, the overwhelming probability bound holds for the eigenvectors of this matrix H giving us a strong form of quantum unique ergodicity for the deformed Gaussian divisible ensemble. In order to remove the small Gaussian component in the matrix, we will use the reverse heat ﬂow technique from [21,23] which allows us to obtain an error as small as we want in total variation between two matrix ensembles. In order to use this technique, we need the smoothness assumption on the matrix W given by Deﬁnition 1.6. We ﬁrst introduce some notation for this section. As before, we let ν denote the distribution of the entries of W , and ϕ denote the density of ν with respect to ρ, the Gaussian distribution with mean zero and variance one, that is, dν = ϕd%. The reverse heat ﬂow technique gives the existence of a probability distribution ν ˜ for any s small enough such that making ν ˜ undergo the Ornstein-Uhlenbeck process of s s generator 1 ∂ x ∂ A := − 2 ∂x 2 ∂x approaches the distribution ν in total variation. This process on all the matrix entries induces the Dyson Brownian motion process on the eigenvalues. Thus the following proposition tells us that there exists a distribution of a matrix from the Gaussian divisible process of the form √ √ f f W = 1− s W + sGOE that approximates as close as polynomially possible a smooth Wigner matrix W . The precise statement is written in the following proposition. Proposition 4.3 ([23]). Let K be a positive integer and ν = ϕ% a distribution smooth in the sense that it follows the conditions (ii) and (iii) of Deﬁnition 1.6. Then there exists s a small positive constant depending on K such that for any 0 < s 6 s , there exists a K K probability density ψ with mean zero and variance one such that we have the inequality sA K e ψ − ϕ d% 6 Cs (4.42) for some positive constant C depending only on K. Besides we also have the inequality for the joint probability of all matrix entries in the following sense, ⊗N 2 2 sA ⊗N ⊗N 2 K e ψ − ϕ d% 6 CN s (4.43) −ε Now, see that this proposition holds for any ﬁxed K so that, taking s = N for some small ε we can choose a large K only depending on ε (and not on N) so that we can obtain any polynomial bound between the two matrix ensembles. This property allows us to get overwhelming probability bounds on the eigenvectors since the total variation distance of the distribution of the eigenvector entries is smaller than the total variation distance between the joint probability of the matrix entries. 42 5. Proofs of main results −1 Now that we have the result for the Gaussian divisible ensemble H with N τ t by Section 3, combining it with the continuity argument from the last subsection, we are able to prove Theorem 1.3 and Corollary 1.5. These two results are a consequence of the following proposition showing the convergence of moments for the eigenvectors of W . Proposition 5.1. Let κ ∈ (0, 1) and m an integer, for a set of indices I ⊂ A , such that |I| = m, we have for any deterministic unit vector q = q , h i N m 2 2 E P |hq, u i| −−−−−→ E P N (5.1) 2 k=1 σ (q, k) N−→∞ k∈I with (N ) a family of independent normal random variables. k k See now the proof of Theorem 1.3 and Corollary 1.5 given by Proposition 3.4. Proof of Theorem 1.3. Proposition 5.1 exactly gives us that the joint moments of the renor- malized eigenvectors converge to those of independent normal random variables which is the result of Theorem 1.3. Proof of Corollary 1.5. By Proposition 3.4 and Corollary 4.2, we have the following inequal- ity, for some ε > 0, −ε 1 N 2 2 E |u (α) | = σ (α, k) +O . (5.2) N Nt By Markov’s inequality, we can write 2 2 X X X X Nt 1 N t 1 2 2 2 2 P |u (α)| − σ (α, k) > c 6 E |u (α)| − σ (α, k) , k k t t 2 2 |A| N c |A| N α∈A α∈A α∈A α∈A 2 2 N t 6 (A− 2B + C) (5.3) 2 2 c |A| We now need to evaluate the three terms in the last inequality using (5.2), ﬁrst we have ! ! 2 2 −ε 2 X X X 1 2 N |A| 2 2 4 A := E |u (α)| = σ (α, k) + σ (α, k) +O . t t 2 2 2 2 N N N t α∈A α∈A α∈A Likewise, " # X X X −ε 2 1 1 N |A| 2 2 2 2 B := σ (β, N )E |u (α)| = σ (α, k)σ (β, k) +O . t t t 2 2 2 N N N t β∈A α∈A α,β∈A Finally, C is just a deterministic term, C := σ (α, k) . α∈A Putting all three terms together, we get −ε 2 −ε 2 2 N |A| |A| N |A| A− 2B + C = σ (α, k) +O 6 C + (5.4) 2 2 2 2 2 2 2 N N t N t N t α∈A The claim then follows from injecting the last inequality in (5.3). We ﬁnish now with the proof of Proposition 5.1. 43 0 Proof of Proposition 5.1. By Corollary 4.2, we know that for some τ ∈ T there exists ε > 0 such that N N 2 H 2 −ε E P |hq, u i| − E P |hq, u i| 6 N (5.5) 2 2 σ (q, k) σ (q, k) t t k∈I k∈I and by Proposition 3.4, we know that, recalling the deﬁnition of H , for some τ t there exists a ε > 0 such that h i N m 0 τ 2 2 −ε E P |hq, u i| − E P N 6 N . (5.6) 2 k=1 σ (q, k) ) k∈I Now, we need to see that H deﬁned in (4.1) has the same law as H for some τ t. τ τ Note that we can write the law of the entries of H as d τ 1 − (ij) 2t t H (τ ) = D + e tW + t 1− e √ N , (5.7) ij ij ij (ij) where N is a family of independent standard Gaussian random variables. Doing the i6j scaling 0 −τ/t t = te = O(t), (5.8) τ = t 1− e = O(τ ) one can write √ √ ˜ 0 0 H = D + t W + τ GOE. Finally, we can apply Proposition 3.4 to H so that (5.6) applies and combining it with (5.5) we get the convergence of moments for the eigenvectors of W . Combining Theorem 3.8 and Proposition 3.4, we are now able to prove Theorem 1.7. Proof of Theorem 1.7. Let ε and D two positive constants and consider s = τ/t,. There exists then a large K, which does not depend on N, such that by Proposition 4.3 there exists a matrix W such that the total variation distance between the distribution of W and −D 1− sW + NsGOE is smaller than N . Denote u , . . . , u the L −normalized eigenvectors of W = D + tW and u ˜ , . . . , u ˜ the 1 N t 1 N p √ f f normalized eigenvectors of W (s) = D + t(1− s) W + tsGOE. Now, since we have in the overwhelming probability bound (1.13) the N degree of liberty, we can do the scaling t = t(1− s) as s 1 and still get (1.13) for the deformed Gaussian divisible ensemble W (s), thus one can write 2 2 ε P u (α) − σ (α, k) > N Ξ(τ ) α∈I ! ! X X 2 2 ε/2 2 ε/2 6 P u ˜ (α) − σ (α, k) > N Ξ(τ ) + P u (α) − u ˜ (α) > N Ξ(τ ) k k k α∈I α∈I −D 6 N where for the last inequality we used the quantum unique ergodicity proved in Theorem 3.8 for the deformed Gaussian divisible ensemble of which u ˜ are the eigenvectors and Proposition 4.3 in order. Now, in order to get the error Ξ we now need to optimize the error 1/3 b b I τ I t Ξ(τ ) = + I = = Ξ with τ = . 0 0 1/3 t N Nτ (Nt) We can do the same thing for the quantity u (α)u (α) and get the ﬁnal result. k l α∈I 44 References [1] R. Allez and J.-P. Bouchaud, Eigenvector dynamics: general theory and some applications, Physical Review E 86 (2012), no. 4, 046202. [2] , Eigenvector dynamics under free addition, Random Matrices: Theory and Applications 3 (2014), no. 03. [3] R. Allez, J. Bun, and J.-P. Bouchaud, The eigenvectors of Gaussian matrices with an external source, arXiv preprint (2014). [4] N. Anantharaman and E. Le Masson, Quantum ergodicity on large regular graphs, Duke Math. J. 164 (2015), no. 4, 723–765. [5] G. W. Anderson, A. Guionnet, and O. Zeitouni, An introduction to random matrices, Cambridge Studies in Advanced Mathematics, vol. 118, Cambridge University Press, Cambridge, 2010. [6] P. Anderson, Absences of diﬀusion in certain random lattices, Phys. Rev. (1958), 1492-1505. [7] R. Bauerschmidt, J. Huang, and H.-T. Yau, Local Kesten–Mckay law for random regular graphs, Comm. Math. Phys. (2019), to appear. [8] F. Benaych-Georges and A. Knowles, Lectures on the local semicircle law for Wigner matrices, arXiv preprint (2016). [9] P. Biane, On the free convolution with a semi-circular distribution, Indiana Univ. Math. J. 46 (1997), no. 3, 705–718. [10] O. Bohigas, M.-J. Giannoni, and C. Schmit, Characterization of chaotic quantum spectra and universality of level ﬂuctuation laws, Phys. Rev. Lett. 52 (1984), no. 1, 1–4. [11] P. Bourgade, Random band matrices, Proceedings ICM-2018 (2018). [12] P. Bourgade, L. Erdős, H.-T. Yau, and J. Yin, Fixed energy universality for generalized Wigner matrices, Comm. Pure Appl. Math. 69 (2016), no. 10, 1815–1881. [13] , Universality for a class of random band matrices, Adv. Theo. Math. Phys. 21 (2017), 739-800. [14] P. Bourgade, J. Huang, and H.-T. Yau, Eigenvector statistics of sparse random matrices, Electron. J. Probab. 22 (2017), Paper No. 64, 38. [15] P. Bourgade and H.-T. Yau, The eigenvector moment ﬂow and local quantum unique ergodicity, Comm. Math. Phys. 350 (2017), no. 1, 231–278. [16] P. Bourgade, H.-T. Yau, and J. Yin, Random band matrices in the delocalized phase, I: Quantum unique ergodicity and universality, arXiv preprint (2018). [17] M.-F. Bru, Diﬀusions of perturbed principal component analysis, J. Multivariate Anal. 29 (1989), no. 1, 127–136. [18] J. Bun, R. Allez, J.-P. Bouchaud, and M. Potters, Rotational Invariant Estimator for General Noisy Matrices, IEEE Trans. Information Theory 62 (2016), no. 12, 7475–7490. [19] J. Bun, J.-P. Bouchaud, and M. Potters, Overlaps between eigenvectors of correlated random matrices, Physical Review E 98 (2018), no. 5. [20] Y. Colin de Verdière, Ergodicité et fonctions propres du laplacien, Comm. Math. Phys. 102 (1985), no. 3, 497–502. [21] L. Erdős, S. Péché, J. A. Ramírez, B. Schlein, and H.-T. Yau, Bulk universality for Wigner matrices, Comm. Pure Appl. Math. 63 (2010), no. 7, 895–925. [22] L. Erdős, B. Schlein, and H.-T. Yau, Local semicircle law and complete delocalization for Wigner random matrices, Comm. Math. Phys. 287 (2009), no. 2, 641–655. [23] , Universality of random matrices and local relaxation ﬂow, Invent. Math. 185 (2011), no. 1, 75–119. [24] L. Erdős and K. Schnelli, Universality for random matrix ﬂows with time-dependent density, Ann. Inst. Henri Poincaré Probab. Stat. 53 (2017), no. 4, 1606–1656. [25] L. Erdős and H.-T. Yau, Gap universality of generalized Wigner and β-ensembles, J. Eur. Math. Soc. (JEMS) 17 (2015), no. 8, 1927–2036. [26] , A dynamical approach to random matrix theory, Courant Lecture Notes in Mathematics, vol. 28, Courant Institute of Mathematical Sciences, New York; American Mathematical Society, Providence, RI, [27] L. Erdős, H.-T. Yau, and J. Yin, Bulk universality for generalized Wigner matrices, Probab. Theory Related Fields 154 (2012), no. 1-2, 341–407. [28] , Rigidity of eigenvalues of generalized Wigner matrices, Adv. Math. 229 (2012), no. 3, 1435–1515. [29] D. Facoetti, P. Vivo, and G. Biroli, From non-ergodic eigenvectors to local resolvent statistics and back: A random matrix perspective, EPL (Europhysics Letters) 115 (2016), no. 4, 47003. [30] R. Holowinsky, Sieving for mass equidistribution, Ann. of Math. (2) 172 (2010), no. 2, 1499–1516. [31] R. Holowinsky and K. Soundararajan, Mass equidistribution for Hecke eigenforms, Ann. of Math. (2) 172 (2010), no. 2, 1517–1528. 45 [32] J. Huang, B. Landon, and H.-T. Yau, Bulk universality of sparse random matrices, J. Math. Phys. 56 (2015), no. 12, 123301, 19. [33] A. Knowles and J. Yin, Eigenvector distribution of Wigner matrices, Probab. Theory Related Fields 155 (2013), no. 3-4, 543–582. [34] V. E. Kravtsov, B. L. Altshuler, and L. B. Ioﬀe, Non-ergodic delocalized phase in Anderson model on Bethe lattice and regular graph, Ann. Phys. 389 (2017). [35] B. Landon and H.-T. Yau, Convergence of local statistics of Dyson Brownian motion, Comm. Math. Phys. 355 (2017), no. 3, 949–1000. [36] J. O. Lee and K. Schnelli, Local deformed semicircle law and complete delocalization for Wigner matrices with random potential, J. Math. Phys. 54 (2013), no. 10, 103504, 62. [37] J. O. Lee, K. Schnelli, B. Stetler, and H.-T. Yau, Bulk universality for deformed Wigner matrices, Ann. Probab. 44 (2016), no. 3, 2349–2425. [38] E. Lindenstrauss, Invariant measures and arithmetic quantum unique ergodicity, Ann. of Math. (2) 163 (2006), no. 1, 165–219. [39] M. L. Mehta, Random matrices, 3rd ed., Pure and Applied Mathematics (Amsterdam), vol. 142, Else- vier/Academic Press, Amsterdam, 2004. [40] J. R. Norris, L. C. G. Rogers, and D. Williams, Brownian motions of ellipsoids, Trans. Amer. Math. Soc. 294 (1986), no. 2, 757–765. [41] K. Truong and A. Ossipov, Eigenvectors under a generic perturbation: Non-perturbative results from the random matrix approach, EPL (Europhysics Letters) 116 (2016), no. 3, 37002. [42] R. Peled, J. Schenker, M. Shamis, and S. Sodin, On the Wegner Orbital Model, Int. Math. Res. Not. IMRN 4 (2019), 1030–1058. [43] N. Rosenzweig and C. E. Porter, "Repulsion of Energy Levels" in Complex Atomic Spectra, Phys. Rev. 120 (1960), no. 5, 1698–1714. [44] Z. Rudnick and P. Sarnak, The behaviour of eigenstates of arithmetic hyperbolic manifolds, Comm. Math. Phys. 161 (1994), no. 1, 195–213. [45] M. Rudelson and R. Vershynin, Delocalization of eigenvectors of random matrices with independent entries, Duke Math. J. 164 (2015), no. 13, 2507–2538. [46] J. Schenker, Eigenvector localization for random band matrices with power law band width, Comm. Math. Phys. 290 (2009), no. 3, 1065–1097. [47] A. I. Šnirel’man, Ergodic properties of eigenfunctions, Uspehi Mat. Nauk 29 (1974), no. 6(180), 181–182. [48] S. Sodin, The spectral edge of some random band matrices, Ann. of Math. (2) 172 (2010), no. 3, 2223– [49] T. Tao and V. Vu, Random matrices: universality of local eigenvalue statistics, Acta Math. 206 (2011), no. 1, 127–204. [50] , Random matrices: universal properties of eigenvectors, Random Matrices Theory Appl. 1 (2012), no. 1, 1150001, 27. [51] P. von Soosten and S. Warzel, The phase transition in the ultrametric ensemble and local stability of Dyson Brownian motion, Electron. J. Probab. 23 (2018), 24 pp. [52] , Non-ergodic delocalization in the Rosenzweig–Porter model, Letters in Mathematical Physics (2018). [53] V. Vu and K. Wang, Random weighted projections, random quadratic forms and random eigenvectors, Random Structures Algorithms 47 (2015), no. 4, 792–821. [54] S. Zelditch, Uniform distribution of eigenfunctions on compact hyperbolic surfaces, Duke Math. J. 55 (1987), no. 4, 919–941.

Mathematics – arXiv (Cornell University)

**Published: ** Nov 19, 2017

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.