We present a generalization of the method of the local relaxation flow to establish the universality of local spectral statistics of a broad class of large random matrices. We show that the local distribution of the eigenvalues coincides with the local statistics of the corresponding Gaussian ensemble provided the distribution of the individual matrix element is smooth and the eigenvalues ${x_j}_{j=1}^N$ are close to their classical location ${\gamma_j}_{j=1}^N$ determined by the limiting density of eigenvalues. Under the scaling where the typical distance between neighboring eigenvalues is of order 1/N, the necessary apriori estimate on the location of eigenvalues requires only to know that $\E |x_j - \gamma_j |^2 \le N^{-1-\e}$ on average. This information can be obtained by well established methods for various matrix ensembles. We demonstrate the method by proving local spectral universality for Wishart matrices.

Consider $N\times N$ hermitian or symmetric random matrices $H$ with independent entries, where the distribution of the $(i,j)$ matrix element is given by the probability measure $\nu_{ij}$ with zero expectation and with variance $\sigma_{ij}^2$. We assume that the variances satisfy the normalization condition $\sum_{i} \sigma^2_{ij} = 1$ for all $j$ and that there is a positive constant $c$ such that $c\le N \sigma_{ij}^2 \le c^{-1}$. We further assume that the probability distributions $\nu_{ij}$ have a uniform subexponential decay. We prove that the Stieltjes transform of the empirical eigenvalue distribution of $H$ is given by the Wigner semicircle law uniformly up to the edges of the spectrum with an error of order $ (N \eta)^{-1}$ where $\eta$ is the imaginary part of the spectral parameter in the Stieltjes transform. There are three corollaries to this strong local semicircle law: (1) Rigidity of eigenvalues: If $\gamma_j =\gamma_{j,N}$ denotes the {\it classical location} of the $j$-th eigenvalue under the semicircle law ordered in increasing order, then the $j$-th eigenvalue $\lambda_j$ is close to $\gamma_j$ in the sense that for any $\xi>1$ there is a constant $L$ such that \[\mathbb P \Big (\exists \, j : \; |\lambda_j-\gamma_j| \ge (\log N)^L \Big [ \min \big (\, j, N-j+1 \, \big) \Big ]^{-1/3} N^{-2/3} \Big) \le C\exp{\big[-c(\log N)^{\xi} \big]} \] for $N$ large enough. (2) The proof of the {\it Dyson's conjecture} \cite{Dy} which states that the time scale of the Dyson Brownian motion to reach local equilibrium is of order $N^{-1}$. (3) The edge universality holds in the sense that the probability distributions of the largest (and the smallest) eigenvalues of two generalized Wigner ensembles are the same in the large $N$ limit provided that the second moments of the two ensembles are identical.

Let $\widetilde{X}_{M\times N}$ be a rectangular data matrix with independent real-valued entries $[\widetilde{x}_{ij}]$ satisfying $\mathbb {E}\widetilde{x}_{ij}=0$ and $\mathbb {E}\widetilde{x}^2_{ij}=\frac{1}{M}$, $N,M\to\infty$. These entries have a subexponential decay at the tails. We will be working in the regime $N/M=d_N,\lim_{N\to\infty}d_N\neq0,1,\infty$. In this paper we prove the edge universality of correlation matrices ${X}^{\dagger}X$, where the rectangular matrix $X$ (called the standardized matrix) is obtained by normalizing each column of the data matrix $\widetilde{X}$ by its Euclidean norm. Our main result states that asymptotically the $k$-point ($k\geq1$) correlation functions of the extreme eigenvalues (at both edges of the spectrum) of the correlation matrix ${X}^{\dagger}X$ converge to those of the Gaussian correlation matrix, that is, Tracy-Widom law, and, thus, in particular, the largest and the smallest eigenvalues of ${X}^{\dagger}X$ after appropriate centering and rescaling converge to the Tracy-Widom distribution. The asymptotic distribution of extreme eigenvalues of the Gaussian correlation matrix has been worked out only recently. As a corollary of the main result in this paper, we also obtain that the extreme eigenvalues of Gaussian correlation matrices are asymptotically distributed according to the Tracy-Widom law. The proof is based on the comparison of Green functions, but the key obstacle to be surmounted is the strong dependence of the entries of the correlation matrix. We achieve this via a novel argument which involves comparing the moments of product of the entries of the standardized data matrix to those of the raw data matrix. Our proof strategy may be extended for proving the edge universality of other random matrix ensembles with dependent entries and hence is of independent interest.

We consider the ensemble of adjacency matrices of Erd{\H o}s-R\'enyi random graphs, i.e.\ graphs on $N$ vertices where every edge is chosen independently and with probability $p \equiv p(N)$. We rescale the matrix so that its bulk eigenvalues are of order one. Under the assumption $p N \gg N^{2/3}$, we prove the universality of eigenvalue distributions both in the bulk and at the edge of the spectrum. More precisely, we prove (1) that the eigenvalue spacing of the Erd{\H o}s-R\'enyi graph in the bulk of the spectrum has the same distribution as that of the Gaussian orthogonal ensemble; and (2) that the second largest eigenvalue of the Erd{\H o}s-R\'enyi graph has the same distribution as the largest eigenvalue of the Gaussian orthogonal ensemble. As an application of our method, we prove the bulk universality of generalized Wigner matrices under the assumption that the matrix entries have at least $4 + \epsilon$ moments.

Consider $N\times N$ Hermitian or symmetric random matrices $H$ where the distribution of the $(i,j)$ matrix element is given by a probability measure $\nu_{ij}$ with a subexponential decay. Let $\sigma_{ij}^2$ be the variance for the probability measure $\nu_{ij}$ with the normalization property that $\sum_{i} \sigma^2_{ij} = 1$ for all $j$. Under essentially the only condition that $c\le N \sigma_{ij}^2 \le c^{-1}$ for some constant $c>0$, we prove that, in the limit $N \to \infty$, the eigenvalue spacing statistics of $H$ in the bulk of the spectrum coincide with those of the Gaussian unitary or orthogonal ensemble (GUE or GOE). We also show that for band matrices with bandwidth $M$ the local semicircle law holds to the energy scale $M^{-1}$.