In this paper, we prove a necessary and sufficient condition for Tracy-Widom law of Wigner matrices. Consider $N \times N$ symmetric Wigner matrices $H$ with $H_{ij} = N^{-1/2} x_{ij}$, whose upper right entries $x_{ij}$ $(1\le i< j\le N)$ are $i.i.d.$ random variables with distribution $\mu$ and diagonal entries $x_{ii}$ $(1\le i\le N)$ are $i.i.d.$ random variables with distribution $\wt \mu$. The means of $\mu$ and $\wt \mu$ are zero, the variance of $\mu$ is 1, and the variance of $\wt \mu $ is finite. We prove that Tracy-Widom law holds if and only if $\lim_{s\to \infty}s^4\p(|x_{12}| \ge s)=0$. The same criterion holds for Hermitian Wigner matrices.

Large dimensional Gram type matrices are common objects in high-dimensional statistics and machine learning. In this paper, we study the limiting distribution of the edge eigenvalues for a general class of high-dimensional Gram type random matrices, including separable sample covariance matrices, sparse sample covariance matrices, bipartite stochastic block model and random Gram matrices with general variance profiles. Specifically, we prove that under (almost) sharp moment conditions and certain tractable regularity assumptions, the edge eigenvalues, i.e., the largest few eigenvalues of non-spiked Gram type random matrices or the extremal bulk eigenvalues of spiked Gram type random matrices, satisfy the Tracy-Widom distribution asymptotically. Our results can be used to construct adaptive, accurate and powerful statistics for high-dimensional statistical inference. In particular, we propose data-dependent statistics to infer the number of signals under general noise structure, test the one-sided sphericity of separable matrix, and test the structure of bipartite stochastic block model. Numerical simulations show strong support of our proposed statistics. The core of our proof is to establish the edge universality and Tracy-Widom distribution for a rectangular Dyson Brownian motion with regular initial data. This is a general strategy to study the edge statistics for high-dimensional Gram type random matrices without exploring the specific independence structure of the target matrices. It has potential to be applied to more general random matrices that are beyond the ones considered in this paper.

In this paper, we prove a necessary and sufficient condition for Tracy-Widom law of Wigner matrices. Consider N×N symmetric Wigner matrices H with Hij=N−1/2xij, whose upper right entries xij (1≤i<j≤N) are i.i.d. random variables with distribution μ and diagonal entries xii (1≤i≤N) are i.i.d. random variables with distribution $\wt \mu$. The means of μ and $\wt \mu$ are zero, the variance of μ is 1, and the variance of $\wt \mu$ is finite. We prove that Tracy-Widom law holds if and only if $\lim_{s\to \infty}s^4\p(|x_{12}| \ge s)=0$. The same criterion holds for Hermitian Wigner matrices.

We analyse the spectrum of additive finite-rank deformations of $N \times N$ Wigner matrices $H$. The spectrum of the deformed matrix undergoes a transition, associated with the creation or annihilation of an outlier, when an eigenvalue $d_i$ of the deformation crosses a critical value $\pm 1$. This transition happens on the scale $|d_i| - 1 \sim N^{-1/3}$. We allow the eigenvalues $d_i$ of the deformation to depend on $N$ under the condition $|\abs{d_i} - 1| \geq (\log N)^{C \log \log N} N^{-1/3}$. We make no assumptions on the eigenvectors of the deformation. In the limit $N \to \infty$, we identify the law of the outliers and prove that the non-outliers close to the spectral edge have a universal distribution coinciding with that of the extremal eigenvalues of a Gaussian matrix ensemble. A key ingredient in our proof is the \emph{isotropic local semicircle law}, which establishes optimal high-probability bounds on the quantity $< v,[(H - z)^{-1} - m(z) 1] w >$, where $m(z)$ is the Stieltjes transform of Wigner's semicircle law and $v, w$ are arbitrary deterministic vectors.

We consider ASEP on a bounded interval and on a half‐line with sources and sinks. On the full line, Bertini and Giacomin in 1997 proved convergence under weakly asymmetric scaling of the height function to the solution of the KPZ equation. We prove here that under similar weakly asymmetric scaling of the sources and sinks as well, the bounded interval ASEP height function converges to the KPZ equation on the unit interval with Neumann boundary conditions on both sides (different parameter for each side), and likewise for the half‐line ASEP to KPZ on a half‐line. This result can be interpreted as showing that the KPZ equation arises at the triple critical point (maximal current / high density / low density) of the open ASEP.

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.

We consider a general class of symmetric or Hermitian random band matrices $H=(h_{xy})_{x,y \in \llbracket 1,N\rrbracket^d}$ in any dimension $d\ge 1$, where the entries are independent, centered random variables with variances $s_{xy}=\mathbb E|h_{xy}|^2$. We assume that $s_{xy}$ vanishes if $|x-y|$ exceeds the band width $W$, and we are interested in the mesoscopic scale with $1\ll W\ll N$. Define the {\it{generalized resolvent}} of $H$ as $G(H,Z):=(H - Z)^{-1}$, where $Z$ is a deterministic diagonal matrix with entries $Z_{xx}\in \mathbb C_+$ for all $x$. Then we establish a precise high-probability bound on certain averages of polynomials of the resolvent entries. As an application of this fluctuation averaging result, we give a self-contained proof for the delocalization of random band matrices in dimensions $d\ge 2$. More precisely, for any fixed $d\ge 2$, we prove that the bulk eigenvectors of $H$ are delocalized in certain averaged sense if $N\le W^{1+\frac{d}{2}}$. This improves the corresponding results in \cite{HeMa2018} under the assumption $N\ll W^{1+\frac{d}{d+1}}$, and in \cite{ErdKno2013,ErdKno2011} under the assumption $N\ll W^{1+\frac{d}{6}}$. For 1D random band matrices, our fluctuation averaging result was used in \cite{PartII,PartI} to prove the delocalization conjecture and bulk universality for random band matrices with $N\ll W^{4/3}$.

We consider the ensemble of adjacency matrices of Erdős–Rényi random graphs, that is, 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. We prove that, as long as $pN\to\infty$ (with a speed at least logarithmic in $N$), the density of eigenvalues of the Erdős–Rényi ensemble is given by the Wigner semicircle law for spectral windows of length larger than $N^{-1}$ (up to logarithmic corrections). As a consequence, all eigenvectors are proved to be completely delocalized in the sense that the $\ell^{\infty}$-norms of the $\ell^{2}$-normalized eigenvectors are at most of order $N^{-1/2}$ with a very high probability. The estimates in this paper will be used in the companion paper [Spectral statistics of Erdős–Rényi graphs II: Eigenvalue spacing and the extreme eigenvalues (2011) Preprint] to prove the universality of eigenvalue distributions both in the bulk and at the spectral edges under the further restriction that $pN\gg N^{2/3}$.

We prove an estimate on the smallest singular value of a multiplicatively and additively deformed random rectangular matrix. Suppose $n\le N \le M \le dN$ for some constant $d\ge 1$. Let $X$ be an $M\times n$ random matrix with independent and identically distributed entries, which have zero mean, unit variance and arbitrarily high moments. Let $T$ be an $N\times M$ deterministic matrix with comparable singular values $c\le s_{N}(T) \le s_{1}(T) \le c^{-1}$ for some constant $c>0$. Let $A$ be an $N\times n$ deterministic matrix with $\|A\|=O(\sqrt{N})$. Then we prove that for any $\epsilon>0$, the smallest singular value of $TX-A$ is larger than $N^{-\epsilon}(\sqrt{N}-\sqrt{n-1})$ with high probability. If we assume further the entries of $X$ have subgaussian decay, then the smallest singular value of $TX-A$ is at least of the order $\sqrt{N}-\sqrt{n-1}$ with high probability, which is an essentially optimal estimate.

We introduce a class of M×M sample covariance matrices Q which subsumes and generalizes several previous models. The associated population covariance matrix Σ=EQ is assumed to differ from the identity by a matrix of bounded rank. All quantities except the rank of Σ−IM may depend on M in an arbitrary fashion. We investigate the principal components, i.e.\ the top eigenvalues and eigenvectors, of Q. We derive precise large deviation estimates on the generalized components ⟨w,ξi⟩ of the outlier and non-outlier eigenvectors ξi. Our results also hold near the so-called BBP transition, where outliers are created or annihilated, and for degenerate or near-degenerate outliers. We believe the obtained rates of convergence to be optimal. In addition, we derive the asymptotic distribution of the generalized components of the non-outlier eigenvectors. A novel observation arising from our results is that, unlike the eigenvalues, the eigenvectors of the principal components contain information about the \emph{subcritical} spikes of Σ. The proofs use several results on the eigenvalues and eigenvectors of the uncorrelated matrix Q, satisfying EQ=IM, as input: the isotropic local Marchenko-Pastur law established in [9], level repulsion, and quantum unique ergodicity of the eigenvectors. The latter is a special case of a new universality result for the joint eigenvalue-eigenvector distribution.