Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Eigenvectors distribution and quantum unique ergodicity for deformed Wigner matrices

Eigenvectors distribution and quantum unique ergodicity for deformed Wigner matrices We analyze the distribution of eigenvectors for mesoscopic, mean-field 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 specific variance. For a well spread initial spectrum, this variance profile 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 flow 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 profil 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 fixé de coordonnées. La preuve se base sur des lois locales a priori données dans [15, 35, 37] et le flot 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-field 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 first proved for Wigner matrices in [33,50] under a matching condition on the first 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-field 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-field type, which also interpolates between delocalized and localized (or partially delocalized) phases, but al- ways with GOE/GUE statistics. It is defined as a perturbation of a potential, consisting of a deterministic diagonal matrix, by a mean field 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 specific, 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 first 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 flow 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 first appeared in [15] for generalized Wigner matrices. It is stated as a high- probability bound showing that eigenvectors are asymptotically flat 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 different 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-field models. In [16], a stronger form of probabilistic quantum unique ergodicity has been actually proved, showing that the eigenvectors mass is asymptotically flat 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 flat but concentrates onto an explicit and deterministic profile with a quantitative error. The key ingredient for this analysis is the Bourgade-Yau eigenvector moment flow [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 difficulties 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 first defined in [35]. Definition 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 fixed 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 first 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, define then for any small positive ω, ω −ω T = N η , N r . ω ? For the spectral domain, take first 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. Definition 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 define 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 diffusion 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 first define 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 define 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). Define the quantiles (γ ) of this measure by i,t 06i6N i,t ρ (x)dx = . −∞ 4 We can also now define 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 define 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 Definition 1.1 and m ∈ N. Let t ∈ T and I ⊂ A be a deterministic (N-dependent) set of m elements. Let W as in Definition 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 define 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 defined up to a phase, we first need to define 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 Definition 1.2 and W as in Definition 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 Definition (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 sufficiently 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 definition or from a Gaussian divisible ensemble. Definition 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 flow 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 define 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 Definition 1.6 and W as ω t in (1.1). Define 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 fluctuations around this deterministic profile. 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 first 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 definition. Definition 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 defined 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. Definition 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 flow was first computed in different 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 fixed matrix H , then 0 0 0 the solution to the dynamics from Definition 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 off-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 flow. Take q ∈ R such that kqk = 1 a fixed direction onto which we will project our eigen- H H vectors. For u , . . . , u the eigenvectors of the matrix (1.15), define 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 configuration 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 configuration ξ. i,j Define ξ to be the configuration 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 configuration ξ with ξ = j and η = 0 if k k i k l l ∈/ {i , . . . , i }. 1 m Define 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 flow that f λ,s undergoes. Theorem 1.12 ([15, Theorem 3.1]). Suppose that u is the solution of the symmetric Dyson vector flow (1.17) and f (ξ) is given by (1.20) with the polynomials Q . Then it satisfies λ,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 flow with a configu- ration of 7 particles. Now that we have the expression of the eigenvector moment flow, we can give an heuristic for the apparition of a Cauchy profile 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, define 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 differential 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. Define 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 configuration of n particles on N sites, define 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)}, define p(e) = p , P (G) = p(e) and finally 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 configuration ξ 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 flow (1.17) and F (ξ) is given by (1.27). Then it satisfies 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 flow 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 different 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 finish 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 flow 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 different 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 ). Define 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 first state a weak local law on the resolvent entries and the Stieltjes transform. The result then follows from a fluctuation −1/2 −1 averaging lemma in order to go from the scale (Nη) to (Nη) . We first give the definition 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 defined 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 Definition 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 Definition 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 off-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 fixed 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 first 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 final 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 . Defining, for a fixed 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 first 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 definition of ≺ and that ϑ is as small as we want and the smoothness of m (z). We finish the proof by definition 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 define 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 Definition 2.8. Let ε, ϑ > 0 and κ ∈ (0, 1). An eigenvalue configuration λ 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 Definition 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 first 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 flow 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 defined 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 first introduced in [25]. Letting 1  `  N be a parameter that we will choose later, we then define 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 finite speed of propagation in the following sense: define the following distance on the set of configurations 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 configurations 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 configurations of n particles and a configuration ξ outside of the support of h in the sense that d(ξ,η) > N ` for any configuration η 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 Definition 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 defined 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 configurations of n particles and a configuration ξ supported outside of the support of h in the sense that d(ξ,η) 6 N ` for any configuration η 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 difference 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 configuration 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 definition 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 configuration 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 configuration ξ 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 finite 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 configurations close to each other, ε 2 namely such that d(ξ,η) 6 N `. Since we have `  Nt, by (3.14), σ (q,ξ, η) varies slowly and we can finally 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 finite 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 first that, by definition 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 final result with the bound |∂ m (z)| 6 C (log N )/t. t t Let ξ ⊂ A be a fixed configuration, 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 define the following flattening 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 first prove in the case of one particle with the dynamics (3.2). We first 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 configuration ξ in this case. However, we need to know that the maximum stays in that window, that is why we first flatten 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 finite speed of propagation (3.1) and we should choose Nt  ` > Nτ for the range cut-off. We will take an explicit value at the end of the proof. The different 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 profile plotted in the spectral dimension: the projection q is fixed and we plot the profile 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 “flat” 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 first 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 different 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 finite 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 final 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 finite speed of propagation, k is in the interval [[k − u, k + u]]. Indeed, if it is not then the difference 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 first 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 differentiable 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 infimum. 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 first define 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 specific window. More precisely define, with ξ being the t 0 22 configuration 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 , define z = γ + iη and let ξ a r k r r configuration 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 configuration 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 first 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 configuration 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 first sum by definition 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 flattened and averaged moment. First define ξ = max g (ξ), (3.54) m τ N (ξ)=n and let k , . . . , k be the positions of the particles of the configuration ξ 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 infimum 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 finite speed of propagation holds. Consider now a deterministic set of indices I ⊂ [[1, N ]]. Note that in the definition 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 fix 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 define, 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 define F (ξ) = U (0, s)Φ (ξ) s B t where B is defined 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 configuration 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 Definition 1.1 and W a Wigner matrix given by 0 −1+a −a Definition 1.2. Let τ ∈ T := [N , N t], then if u , . . . , u are the eigenvectors of the a 1 N matrix H = D + tW + τ GOE, define 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 specific τ and optimize all our different parameters when using the reverse heat flow technique to remove this small Gaussian component. One of these parameters will be a cut-off 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 defined respectively in (3.3) and (3.4). Lemma 3.1 will help us localize the dynamics onto a small set of configurations. Now the following lemma says that most of 26 the information of the dynamics is given by the short-range, bounding the difference between B and S . It is analogous to Lemma 3.5 in [16]. First define (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 configuration ξ 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 defined 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]. Define the following flattening operator. For f a function on configurations of n particles and ξ such a configuration, 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 definition of the flattening 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 finite 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 configuration of n > m s 1 1 m m particles supported in J . out We have, by definition 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 first sum in (3.63), we will first restrict it to the sites k such that there are no particles i ,k in the configuration ξ, so that we will have η ˜ = 1. Note that, if we have η ˜ 6= 0 then by 27 i ,k ˜ ˜ p definition 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 definition 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 different types of perfect matchings. Define 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 first 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, first define 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 first 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 definition of p and of σ and defined 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), first 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 first 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 finite 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. Define the following averaging operator. For f a function on configurations and ξ a configuration, ε 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 definition of the flattening 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. Define now ξ ξ ξ the following dynamics ∂ Γ = S (s)Γ , 0 6 s 6 τ s s s (3.73) Γ (ξ) = (Av Φ )(ξ). 0 t 30 Now, if one takes a configuration ξ supported on J , it suffices 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 first term by using Lemma 3.9 and the second term using the finite speed of propagation. Indeed, since ξ is supported on J , (Id − Av)Φ vanishes for any in t configuration 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 definition of F . However, if −C this supremum is greater than N , then by the finite 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 first 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 first term, we will use the finite 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 definition 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 definition 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 }, first note that there can only be In 2n`N such a, then one can see that we can use the finite 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 defined in (2.13). For (3.78), we will use Lemma 3.9. Indeed, first 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 first term, see first 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 definition 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 definition 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 defined 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 first 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 define the following sequence of sets of indices, defined 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 definition 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). Define, 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 first 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 confirmed 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 first 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 definition 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 off-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 first 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 definition 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 first inequality, the fact that |g | 6 Ct in the second and (4.9) again in the final 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 first 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 finishing the bound by Lemma 4.1. Again, now that we have any polynomial of fixed 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 defined 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 first 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 first 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 first need to use −1+ξ local laws which hold down to mesoscopic scales η = N . Now for the third derivative of F , first 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 first 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 definition 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 off-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 final 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 first 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 flow 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 definition 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 flow 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 Definition 1.6. We first 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 flow 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 Definition 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 fixed 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), first 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 finish 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 definition 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 defined 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 final 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 diffusion 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 fluctuation 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 flow 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, Diffusions 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 flow, Invent. Math. 185 (2011), no. 1, 75–119. [24] L. Erdős and K. Schnelli, Universality for random matrix flows 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. Ioffe, 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. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Mathematics arXiv (Cornell University)

Eigenvectors distribution and quantum unique ergodicity for deformed Wigner matrices

Benigni, Lucas
Mathematics , Volume 2020 (1711) – Nov 19, 2017
Download PDF
46 pages

Loading next page...
 
/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

Abstract

We analyze the distribution of eigenvectors for mesoscopic, mean-field 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 specific variance. For a well spread initial spectrum, this variance profile 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 flow 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 profil 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 fixé de coordonnées. La preuve se base sur des lois locales a priori données dans [15, 35, 37] et le flot 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-field 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 first proved for Wigner matrices in [33,50] under a matching condition on the first 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-field 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-field type, which also interpolates between delocalized and localized (or partially delocalized) phases, but al- ways with GOE/GUE statistics. It is defined as a perturbation of a potential, consisting of a deterministic diagonal matrix, by a mean field 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 specific, 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 first 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 flow 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 first appeared in [15] for generalized Wigner matrices. It is stated as a high- probability bound showing that eigenvectors are asymptotically flat 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 different 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-field models. In [16], a stronger form of probabilistic quantum unique ergodicity has been actually proved, showing that the eigenvectors mass is asymptotically flat 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 flat but concentrates onto an explicit and deterministic profile with a quantitative error. The key ingredient for this analysis is the Bourgade-Yau eigenvector moment flow [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 difficulties 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 first defined in [35]. Definition 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 fixed 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 first 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, define then for any small positive ω, ω −ω T = N η , N r . ω ? For the spectral domain, take first 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. Definition 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 define 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 diffusion 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 first define 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 define 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). Define the quantiles (γ ) of this measure by i,t 06i6N i,t ρ (x)dx = . −∞ 4 We can also now define 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 define 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 Definition 1.1 and m ∈ N. Let t ∈ T and I ⊂ A be a deterministic (N-dependent) set of m elements. Let W as in Definition 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 define 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 defined up to a phase, we first need to define 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 Definition 1.2 and W as in Definition 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 Definition (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 sufficiently 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 definition or from a Gaussian divisible ensemble. Definition 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 flow 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 define 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 Definition 1.6 and W as ω t in (1.1). Define 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 fluctuations around this deterministic profile. 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 first 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 definition. Definition 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 defined 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. Definition 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 flow was first computed in different 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 fixed matrix H , then 0 0 0 the solution to the dynamics from Definition 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 off-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 flow. Take q ∈ R such that kqk = 1 a fixed direction onto which we will project our eigen- H H vectors. For u , . . . , u the eigenvectors of the matrix (1.15), define 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 configuration 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 configuration ξ. i,j Define ξ to be the configuration 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 configuration ξ with ξ = j and η = 0 if k k i k l l ∈/ {i , . . . , i }. 1 m Define 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 flow that f λ,s undergoes. Theorem 1.12 ([15, Theorem 3.1]). Suppose that u is the solution of the symmetric Dyson vector flow (1.17) and f (ξ) is given by (1.20) with the polynomials Q . Then it satisfies λ,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 flow with a configu- ration of 7 particles. Now that we have the expression of the eigenvector moment flow, we can give an heuristic for the apparition of a Cauchy profile 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, define 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 differential 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. Define 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 configuration of n particles on N sites, define 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)}, define p(e) = p , P (G) = p(e) and finally 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 configuration ξ 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 flow (1.17) and F (ξ) is given by (1.27). Then it satisfies 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 flow 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 different 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 finish 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 flow 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 different 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 ). Define 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 first state a weak local law on the resolvent entries and the Stieltjes transform. The result then follows from a fluctuation −1/2 −1 averaging lemma in order to go from the scale (Nη) to (Nη) . We first give the definition 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 defined 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 Definition 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 Definition 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 off-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 fixed 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 first 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 final 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 . Defining, for a fixed 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 first 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 definition of ≺ and that ϑ is as small as we want and the smoothness of m (z). We finish the proof by definition 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 define 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 Definition 2.8. Let ε, ϑ > 0 and κ ∈ (0, 1). An eigenvalue configuration λ 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 Definition 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 first 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 flow 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 defined 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 first introduced in [25]. Letting 1  `  N be a parameter that we will choose later, we then define 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 finite speed of propagation in the following sense: define the following distance on the set of configurations 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 configurations 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 configurations of n particles and a configuration ξ outside of the support of h in the sense that d(ξ,η) > N ` for any configuration η 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 Definition 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 defined 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 configurations of n particles and a configuration ξ supported outside of the support of h in the sense that d(ξ,η) 6 N ` for any configuration η 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 difference 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 configuration 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 definition 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 configuration 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 configuration ξ 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 finite 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 configurations close to each other, ε 2 namely such that d(ξ,η) 6 N `. Since we have `  Nt, by (3.14), σ (q,ξ, η) varies slowly and we can finally 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 finite 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 first that, by definition 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 final result with the bound |∂ m (z)| 6 C (log N )/t. t t Let ξ ⊂ A be a fixed configuration, 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 define the following flattening 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 first prove in the case of one particle with the dynamics (3.2). We first 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 configuration ξ in this case. However, we need to know that the maximum stays in that window, that is why we first flatten 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 finite speed of propagation (3.1) and we should choose Nt  ` > Nτ for the range cut-off. We will take an explicit value at the end of the proof. The different 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 profile plotted in the spectral dimension: the projection q is fixed and we plot the profile 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 “flat” 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 first 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 different 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 finite 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 final 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 finite speed of propagation, k is in the interval [[k − u, k + u]]. Indeed, if it is not then the difference 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 first 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 differentiable 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 infimum. 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 first define 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 specific window. More precisely define, with ξ being the t 0 22 configuration 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 , define z = γ + iη and let ξ a r k r r configuration 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 configuration 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 first 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 configuration 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 first sum by definition 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 flattened and averaged moment. First define ξ = max g (ξ), (3.54) m τ N (ξ)=n and let k , . . . , k be the positions of the particles of the configuration ξ 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 infimum 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 finite speed of propagation holds. Consider now a deterministic set of indices I ⊂ [[1, N ]]. Note that in the definition 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 fix 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 define, 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 define F (ξ) = U (0, s)Φ (ξ) s B t where B is defined 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 configuration 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 Definition 1.1 and W a Wigner matrix given by 0 −1+a −a Definition 1.2. Let τ ∈ T := [N , N t], then if u , . . . , u are the eigenvectors of the a 1 N matrix H = D + tW + τ GOE, define 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 specific τ and optimize all our different parameters when using the reverse heat flow technique to remove this small Gaussian component. One of these parameters will be a cut-off 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 defined respectively in (3.3) and (3.4). Lemma 3.1 will help us localize the dynamics onto a small set of configurations. Now the following lemma says that most of 26 the information of the dynamics is given by the short-range, bounding the difference between B and S . It is analogous to Lemma 3.5 in [16]. First define (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 configuration ξ 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 defined 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]. Define the following flattening operator. For f a function on configurations of n particles and ξ such a configuration, 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 definition of the flattening 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 finite 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 configuration of n > m s 1 1 m m particles supported in J . out We have, by definition 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 first sum in (3.63), we will first restrict it to the sites k such that there are no particles i ,k in the configuration ξ, so that we will have η ˜ = 1. Note that, if we have η ˜ 6= 0 then by 27 i ,k ˜ ˜ p definition 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 definition 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 different types of perfect matchings. Define 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 first 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, first define 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 first 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 definition of p and of σ and defined 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), first 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 first 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 finite 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. Define the following averaging operator. For f a function on configurations and ξ a configuration, ε 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 definition of the flattening 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. Define now ξ ξ ξ the following dynamics ∂ Γ = S (s)Γ , 0 6 s 6 τ s s s (3.73) Γ (ξ) = (Av Φ )(ξ). 0 t 30 Now, if one takes a configuration ξ supported on J , it suffices 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 first term by using Lemma 3.9 and the second term using the finite speed of propagation. Indeed, since ξ is supported on J , (Id − Av)Φ vanishes for any in t configuration 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 definition of F . However, if −C this supremum is greater than N , then by the finite 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 first 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 first term, we will use the finite 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 definition 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 definition 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 }, first note that there can only be In 2n`N such a, then one can see that we can use the finite 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 defined in (2.13). For (3.78), we will use Lemma 3.9. Indeed, first 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 first term, see first 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 definition 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 definition 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 defined 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 first 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 define the following sequence of sets of indices, defined 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 definition 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). Define, 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 first 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 confirmed 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 first 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 definition 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 off-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 first 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 definition 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 first inequality, the fact that |g | 6 Ct in the second and (4.9) again in the final 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 first 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 finishing the bound by Lemma 4.1. Again, now that we have any polynomial of fixed 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 defined 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 first 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 first 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 first need to use −1+ξ local laws which hold down to mesoscopic scales η = N . Now for the third derivative of F , first 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 first 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 definition 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 off-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 final 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 first 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 flow 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 definition 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 flow 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 Definition 1.6. We first 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 flow 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 Definition 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 fixed 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), first 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 finish 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 definition 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 defined 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 final 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 diffusion 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 fluctuation 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 flow 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, Diffusions 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 flow, Invent. Math. 185 (2011), no. 1, 75–119. [24] L. Erdős and K. Schnelli, Universality for random matrix flows 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. Ioffe, 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.

Journal

Mathematics – arXiv (Cornell University)

Published: Nov 19, 2017

There are no references for this article.