In this paper, we consider the universality of the local eigenvalue statistics of random matrices. Our main result shows that these statistics are determined by the first four moments of the distribution of the entries. As a consequence, we derive the universality of eigenvalue gap distribution and$k$-point correlation, and many other statistics (under some mild assumptions) for both Wigner Hermitian matrices and Wigner real symmetric matrices.

The universality for the eigenvalue spacing statistics of generalized Wigner matrices was established in our previous work \cite{EYY} under certain conditions on the probability distributions of the matrix elements. A major class of probability measures excluded in \cite{EYY} are the Bernoulli measures. In this paper, we extend the universality result of \cite{EYY} to include the Bernoulli measures so that the only restrictions on the probability distributions of the matrix elements are the subexponential decay and the normalization condition that the variances in each row sum up to one. The new ingredient is a strong local semicircle law which improves the error estimate on the Stieltjes transform of the empirical measure of the eigenvalues from the order $(N \eta)^{-1/2}$ to $(N \eta)^{-1}$. Here $\eta$ is the imaginary part of the spectral parameter in the definition of the Stieltjes transform and $N$ is the size of the matrix.

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.