Convergence of spectral measures and eigenvalue rigidity

Similar documents
Concentration for Coulomb gases

Concentration for Coulomb gases and Coulomb transport inequalities

Eigenvalue variance bounds for Wigner and covariance random matrices

Concentration Inequalities for Random Matrices

Exponential tail inequalities for eigenvalues of random matrices

Lecture I: Asymptotics for large GUE random matrices

Lectures 6 7 : Marchenko-Pastur Law

Random Matrices: Beyond Wigner and Marchenko-Pastur Laws

Triangular matrices and biorthogonal ensembles

The norm of polynomials in large random matrices

Random regular digraphs: singularity and spectrum

Operator norm convergence for sequence of matrices and application to QIT

Second Order Freeness and Random Orthogonal Matrices

Markov operators, classical orthogonal polynomial ensembles, and random matrices

Determinantal point processes and random matrix theory in a nutshell

Bulk scaling limits, open questions

COMPLEX HERMITE POLYNOMIALS: FROM THE SEMI-CIRCULAR LAW TO THE CIRCULAR LAW

Distribution of Eigenvalues of Weighted, Structured Matrix Ensembles

The circular law. Lewis Memorial Lecture / DIMACS minicourse March 19, Terence Tao (UCLA)

Valerio Cappellini. References

Microscopic behavior for β-ensembles: an energy approach

1 Intro to RMT (Gene)

Random Matrix: From Wigner to Quantum Chaos

STAT C206A / MATH C223A : Stein s method and applications 1. Lecture 31

Random Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws. Symeon Chatzinotas February 11, 2013 Luxembourg

DISTRIBUTION OF EIGENVALUES OF REAL SYMMETRIC PALINDROMIC TOEPLITZ MATRICES AND CIRCULANT MATRICES

Quantum Chaos and Nonunitary Dynamics

Lectures 2 3 : Wigner s semicircle law

Concentration inequalities for non-lipschitz functions

Large Deviations for Random Matrices and a Conjecture of Lukic

Local law of addition of random matrices

From the mesoscopic to microscopic scale in random matrix theory

On asymptotics of large Haar distributed unitary matrices 1

Concentration inequalities: basics and some new challenges

Random matrices: Distribution of the least singular value (via Property Testing)

Universality of sine-kernel for Wigner matrices with a small Gaussian perturbation

Ginibre Random Field Hirofumi Osada (Kyushu University)

A Note on the Central Limit Theorem for the Eigenvalue Counting Function of Wigner and Covariance Matrices

Random Matrices: Invertibility, Structure, and Applications

Universality in Numerical Computations with Random Data. Case Studies.

DLR equations for the Sineβ process and applications

Invertibility of random matrices

RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS

Ferromagnets and the classical Heisenberg model. Kay Kirkpatrick, UIUC

From random matrices to free groups, through non-crossing partitions. Michael Anshelevich

Stochastic Differential Equations Related to Soft-Edge Scaling Limit

Products of Rectangular Gaussian Matrices

SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS

Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium

Fluctuations from the Semicircle Law Lecture 4

Random matrices: A Survey. Van H. Vu. Department of Mathematics Rutgers University

SPECTRAL MEASURE OF LARGE RANDOM HANKEL, MARKOV AND TOEPLITZ MATRICES

Wigner s semicircle law

Homogenization of the Dyson Brownian Motion

Comparison Method in Random Matrix Theory

Universality for random matrices and log-gases

Inhomogeneous circular laws for random matrices with non-identically distributed entries

Rigidity of the 3D hierarchical Coulomb gas. Sourav Chatterjee

Homework 1. Yuan Yao. September 18, 2011

ON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES

Local semicircle law, Wegner estimate and level repulsion for Wigner random matrices

ON THE CONVERGENCE OF THE NEAREST NEIGHBOUR EIGENVALUE SPACING DISTRIBUTION FOR ORTHOGONAL AND SYMPLECTIC ENSEMBLES

Random Fermionic Systems

Central Limit Theorems for linear statistics for Biorthogonal Ensembles

Eigenvalues and Singular Values of Random Matrices: A Tutorial Introduction

Freeness and the Transpose

Dependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.

NUMERICAL SOLUTION OF DYSON BROWNIAN MOTION AND A SAMPLING SCHEME FOR INVARIANT MATRIX ENSEMBLES

The Matrix Dyson Equation in random matrix theory

Random matrix theory and log-correlated Gaussian fields

Differential Equations for Dyson Processes

Susceptible-Infective-Removed Epidemics and Erdős-Rényi random

EECS 598: Statistical Learning Theory, Winter 2014 Topic 11. Kernels

Approximation of BSDEs using least-squares regression and Malliavin weights

Matrix estimation by Universal Singular Value Thresholding

Large deviations of the top eigenvalue of random matrices and applications in statistical physics

Lecture 13 October 6, Covering Numbers and Maurey s Empirical Method

SEA s workshop- MIT - July 10-14

Free Probability and Random Matrices: from isomorphisms to universality

NONCOLLIDING BROWNIAN MOTIONS AND HARISH-CHANDRA FORMULA

Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½

BULK SPECTRUM: WIGNER AND MARCHENKO-PASTUR THEOREMS

The following definition is fundamental.

Large deviations and stochastic calculus for large random matrices

Consistency of Modularity Clustering on Random Geometric Graphs

Statistical mechanics of random billiard systems

Preface to the Second Edition...vii Preface to the First Edition... ix

Random Toeplitz Matrices

Lecture III: Applications of Voiculescu s random matrix model to operator algebras

MATH 505b Project Random Matrices

ORTHOGONAL POLYNOMIALS

ELEMENTS OF PROBABILITY THEORY

The largest eigenvalues of the sample covariance matrix. in the heavy-tail case

Math Solutions to homework 5

Concentration inequalities and the entropy method

Stein s Method: Distributional Approximation and Concentration of Measure

Research Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance

Determinantal Processes And The IID Gaussian Power Series

arxiv: v5 [math.na] 16 Nov 2017

arxiv: v3 [math-ph] 29 Aug 2017

Transcription:

Convergence of spectral measures and eigenvalue rigidity Elizabeth Meckes Case Western Reserve University ICERM, March 1, 2018

Macroscopic scale: the empirical spectral measure

Macroscopic scale: the empirical spectral measure Suppose that M is an n n random matrix with eigenvalues λ 1,..., λ n. The empirical spectral measure µ of M is the (random) measure µ := 1 n n δ λk. k=1

Macroscopic scale: the empirical spectral measure Suppose that M is an n n random matrix with eigenvalues λ 1,..., λ n. The empirical spectral measure µ of M is the (random) measure µ := 1 n n δ λk. k=1

Wigner s Theorem

Wigner s Theorem For each n N, let {Y i } 1 i, {Z ij } 1 i<j be independent collections of i.i.d. random variables, with EY 1 = EZ 12 = 0 EZ 2 12 = 1 EY 2 1 <. Let M n be the symmetric random matrix with diagonal entries Y i and off-diagonal entries Z ij or Z ji. The empirical spectral measure µ n of 1 n M n is close, for large n, to the semi-circular law: 1 4 x 2π 2 1 x 2 dx.

Other examples

Other examples The circular law (Ginibre): The empirical spectral measure of a large random matrix with i.i.d. Gaussian entries is approximately uniform on a disc.

Other examples The circular law (Ginibre): The empirical spectral measure of a large random matrix with i.i.d. Gaussian entries is approximately uniform on a disc. The classical compact groups (Diaconis Shahshahani): The empirical spectral measure of a uniform random matrix in O (n), U (n), Sp (2n) is approximately uniform on the unit circle when n is large.

Other examples Truncations of random unitary matrices (Petz Reffy): Let U m be the upper-left m m block of a uniform random matrix in U (n), and let α = m n.

Other examples Truncations of random unitary matrices (Petz Reffy): Let U m be the upper-left m m block of a uniform random matrix in U (n), and let α = m n. For large n, the empirical spectral measure of U m is close to the measure with density f α (z) = { 2(1 α), α(1 z 2 ) 2 0 < z < α; 0, otherwise.

Other examples Truncations of random unitary matrices (Petz Reffy): Let U m be the upper-left m m block of a uniform random matrix in U (n), and let α = m n. For large n, the empirical spectral measure of U m is close to the measure with density f α (z) = { 2(1 α), α(1 z 2 ) 2 0 < z < α; 0, otherwise. α = 4 5 α = 2 5 Figures from Truncations of random unitary matrices, Życzkowski Sommers, J. Phys. A, 2000

Other examples Brownian motion on U (n) (Biane): Let {U t } t 0 be a Brownian motion on U (n); i.e., a solution to du t = U t dw t 1 2 U tdt, with U 0 = I and W t a standard B.M. on u(n). There is a deterministic family of measures {ν t } t 0 on the unit circle such that the spectral measure of U t converges weakly almost surely to ν t.

Other examples Brownian motion on U (n):

Levels of randomness

Levels of randomness Let µ n be the (random) spectral measure of an n n random matrix, and let ν be some deterministic measure which supposedly approximates µ n.

Levels of randomness Let µ n be the (random) spectral measure of an n n random matrix, and let ν be some deterministic measure which supposedly approximates µ n. The annealed case: The ensemble-averaged spectral measure is Eµ n : fd(eµ n ) := E fdµ n.

Levels of randomness Let µ n be the (random) spectral measure of an n n random matrix, and let ν be some deterministic measure which supposedly approximates µ n. The annealed case: The ensemble-averaged spectral measure is Eµ n : fd(eµ n ) := E fdµ n. One may prove that Eµ n ν, possibly via explicit bounds on d(eµ n, ν) in some metric d(, ).

Levels of randomness The quenched case:

Levels of randomness The quenched case: Convergence weakly in probability or weakly almost surely: for any bounded continuous test function f, P a.s. f dµ n f dν or f dµ n f dν.

Levels of randomness The quenched case: Convergence weakly in probability or weakly almost surely: for any bounded continuous test function f, P a.s. f dµ n f dν or f dµ n f dν. The random variable d(µ n, ν): Look for ɛ n such that with high probability (or even probability 1), d(µ n, ν) < ɛ n.

Microscopic scale: eigenvalue rigidity

Microscopic scale: eigenvalue rigidity In many settings, eigenvalues concentrate strongly about predicted locations.

Microscopic scale: eigenvalue rigidity In many settings, eigenvalues concentrate strongly about predicted locations. 1.0 0.5-1.0-0.5 0.5 1.0-0.5-1.0

The eigenvalues of U m for m = 1, 5, 20, 45, 80, for U a realization of a random 80 80 unitary matrix.

Theorem (E. M. M. Meckes) Let 0 θ 1 < θ 2 < < θ n < 2π be the eigenvalue angles of U p, where U is a Haar random matrix in U (n). For each j and t > 0, P [ θ j 2πj N > 4π ] N t t 4 exp min ( 2 p log N p ) + 1, t.

Concentration of empirical spectral measures 2-D Coulomb gases

Concentration of empirical spectral measures 2-D Coulomb gases Coulomb transport inequality (Chafaï Hardy Maïda): Consider the 2-D Coulomb gas model with Hamiltonian H n (z 1,..., z n ) = j k log z j z k + n n V (z j ); let µ V denote the equilibrium measure. There is a constant C V such that j=1 d BL (µ, µ V ) 2 W 1 (µ, µ V ) 2 C V [E V (µ) E V (µ V )], where E V is the modified energy functional E V (µ) = E(µ) + Vdµ.

Truncations of random unitary matrices Let U be distributed according to Haar measure in U (n) and let n 1 m n. Let U m denote the top-left m m block of m U. The eigenvalue density of U m is given by 1 c n,m 1 j<k m z j z k 2 m j=1 ( 1 m n z j 2) n m 1 dλ(z1 ) dλ(z n ), which corresponds to a two-dimensional Coulomb gas with external potential n m 1 m log ( 1 m n Ṽ n,m (z) = z 2) n. z < m ; n, z m.

Truncations of random unitary matrices Theorem (M. Lockwood) Let µ m,n be the spectral measure of the top-left m m block of n m U, where U is a random n n unitary matrix and 1 m n 2 log(n). Let α = m n, and let ν α have density then g α (z) = { 2(1 α), (1 α z 2 ) 2 0 < z < 1; 0, otherwise. P[d BL (µ m,n, ν α ) > r] e Cαm2 r 2 +2m[log(m)+C α] + e cn, { } 1 where C α = min log(α 1 ), 1 and C α { log( 1 α ), α 0; log(1 α), α 1.

Concentration of empirical spectral measures Ensembles with concentration properties

Concentration of empirical spectral measures Ensembles with concentration properties If M is an n n normal matrix with spectral measure µ M and f : C R is 1-Lipschitz, it follows from the Hoffman-Wielandt inequality that M fdµ M is a 1 n -Lipschitz function of M.

Concentration of empirical spectral measures Ensembles with concentration properties If M is an n n normal matrix with spectral measure µ M and f : C R is 1-Lipschitz, it follows from the Hoffman-Wielandt inequality that M fdµ M is a 1 n -Lipschitz function of M. = For any reference measure ν, is 1 n -Lipschitz M W 1 (µ M, ν)

Concentration of empirical spectral measures Many random matrix ensembles satisfy the following concentration property: Let F : S M N R be 1-Lipschitz with respect to H.S.. Then [ F(M) ] P EF(M) > t Ce cnt2.

Concentration of empirical spectral measures Many random matrix ensembles satisfy the following concentration property: Let F : S M N R be 1-Lipschitz with respect to H.S.. Then [ F(M) ] P EF(M) > t Ce cnt2. Some Examples:

Concentration of empirical spectral measures Many random matrix ensembles satisfy the following concentration property: Let F : S M N R be 1-Lipschitz with respect to H.S.. Then [ F(M) ] P EF(M) > t Ce cnt2. Some Examples: GUE; Wigner matrices in which the entries satisfy a quadratic transportation cost inequality with constant c N.

Concentration of empirical spectral measures Many random matrix ensembles satisfy the following concentration property: Let F : S M N R be 1-Lipschitz with respect to H.S.. Then [ F(M) ] P EF(M) > t Ce cnt2. Some Examples: GUE; Wigner matrices in which the entries satisfy a quadratic transportation cost inequality with constant Wishart (sort of) c N.

Concentration of empirical spectral measures Many random matrix ensembles satisfy the following concentration property: Let F : S M N R be 1-Lipschitz with respect to H.S.. Then [ F(M) ] P EF(M) > t Ce cnt2. Some Examples: GUE; Wigner matrices in which the entries satisfy a quadratic transportation cost inequality with constant Wishart (sort of) c N. Haar measure and heat kernel measure on the compact classical groups: SO (N), U (N), SU (N), Sp (2N)

Concentration of empirical spectral measures Many random matrix ensembles satisfy the following concentration property: Let F : S M N R be 1-Lipschitz with respect to H.S.. Then [ F(M) ] P EF(M) > t Ce cnt2. Some Examples: GUE; Wigner matrices in which the entries satisfy a quadratic transportation cost inequality with constant Wishart (sort of) c N. Haar measure and heat kernel measure on the compact classical groups: SO (N), U (N), SU (N), Sp (2N) Ensembles with matrix density e N Tr(u(M)), with u (x) c > 0.

Typical vs. average

Typical vs. average In ensembles with the concentration property, W 1 (µ n, ν), this means P[W 1 (µ n, ν) > EW 1 (µ n, ν) + t] Ce cn2 t 2.

Typical vs. average In ensembles with the concentration property, W 1 (µ n, ν), this means P[W 1 (µ n, ν) > EW 1 (µ n, ν) + t] Ce cn2 t 2. To show W 1 (µ n, ν) is typically small, it s enough to show that EW 1 (µ n, ν) is small.

Average distance to average One approach: consider the stochastic process X f := fdµ n E fdµ n.

Average distance to average One approach: consider the stochastic process X f := fdµ n E fdµ n. Under the concentration hypothesis, {X f } f satisfies a sub-gaussian increment condition: P [ X f X g > t ] 2e cn 2 t 2 f g 2 L. Dudley s entropy bound together with approximation theory, truncation arguments, etc., can lead to a bound on ( ) EW 1 (µ n, Eµ n ) = E sup X f. f L 1

Theorem (M. Melcher) Let µ N t be the spectral measure of U t, where {U t } t 0 is a Brownian motion on U (n) with U 0 = I.

Theorem (M. Melcher) Let µ N t be the spectral measure of U t, where {U t } t 0 is a Brownian motion on U (n) with U 0 = I. 1. For any t, x > 0, ( P W 1 (µ N t, µ N t ) > c ( ) t 1/3 N 2 + x) 2e N2 x 2 t.

Theorem (M. Melcher) Let µ N t be the spectral measure of U t, where {U t } t 0 is a Brownian motion on U (n) with U 0 = I. 1. For any t, x > 0, ( P W 1 (µ N t, µ N t ) > c ( ) t 1/3 N 2 + x) 2e N2 x 2 t. 2. There are constants c, C such that for T 0 and x c T 2/5 log(n) P ( N 2/5, sup W 1 (µ N t, ν t ) > x 0 t T ) ( ) T C x 2 + 1 e N2 x 2 T. In particular, with probability one for N sufficiently large sup W 1 (µ N t, ν t ) c T 2/5 log(n) 0 t T N 2/5.

Second approach: Eigenvalue rigidity!

Second approach: Eigenvalue rigidity! The set of eigenvalues of many types of random matrices are determinantal point processes with symmetric kernels:

Second approach: Eigenvalue rigidity! The set of eigenvalues of many types of random matrices are determinantal point processes with symmetric kernels: GUE U (N) Complex Ginibre j=0 K N (x, y) n 1 h j (x)h j (y)e (x2 +y 2 ) 2 R 1 π N 1 j=0 N 1 j=0 (zw) j j! Λ e ij(x y) [0, 2π) e ( z 2 + w 2 ) 2 { z = 1}

The gift of determinantal point processes

The gift of determinantal point processes Theorem (Hough/Krishnapur/Peres/Virág) Let K : Λ Λ C be the kernel of a determinantal point process, and suppose the corresponding integral operator is self-adjoint, nonnegative, and locally trace-class. For D Λ, let N D denote the number of particles of the point process in D. Then d N D = ξ k, where {ξ k } is a collection of independent Bernoulli random variables. k

Concentration of the counting function

Concentration of the counting function Since N D is a sum of i.i.d. Bernoullis, Bernstein s inequality applies: ( { }) P [ N D EN D > t] 2 exp min where σd 2 = Var N D. t 2 4σ 2 D, t, 2

Rigidity of individual eigenvalues

Rigidity of individual eigenvalues If U is a random unitary matrix, then U has eigenvalues {e iθ k } N k=1, for 0 θ 1 < θ 2 < < θ N < 2π.

Rigidity of individual eigenvalues If U is a random unitary matrix, then U has eigenvalues {e iθ k } N k=1, for 0 θ 1 < θ 2 < < θ N < 2π. We define the predicted locations to be {e 2πik N } N k=1.

Rigidity of individual eigenvalues If U is a random unitary matrix, then U has eigenvalues {e iθ k } N k=1, for 0 θ 1 < θ 2 < < θ N < 2π. We define the predicted locations to be {e 2πik N } N k=1. concentration of N [0,θ] concentration of e iθ k about e ik N

Rigidity of individual eigenvalues If U is a random unitary matrix, then U has eigenvalues {e iθ k } N k=1, for 0 θ 1 < θ 2 < < θ N < 2π. We define the predicted locations to be {e 2πik N } N k=1. concentration of N [0,θ] concentration of e iθ k about e ik N EW 1 (µ N, ν) C log(n) + 1, N where ν is the uniform distribution on S 1 C.

Co-authors Kathryn Lockwood (Ph.D. student, CWRU): truncations of random unitary matrices Tai Melcher (UVA): Brownian motion on U (n) Mark Meckes (CWRU): most of the rest