We consider sample covariance matrices of the form $X^*X$, where $X$ is an $M times N$ matrix with independent random entries. We prove the isotropic local Marchenko-Pastur law, i.e. we prove that the resolvent $(X^* X - z)^{-1}$ converges to a multiple of the identity in the sense of quadratic forms. More precisely, we establish sharp high-probability bounds on the quantity $langle v, (X^* X - z)^{-1} w rangle - langle v,w rangle m(z)$, where $m$ is the Stieltjes transform of the Marchenko-Pastur law and $v, w in mathbb C^N$. We require the logarithms of the dimensions $M$ and $N$ to be comparable. Our result holds down to scales $Im z geq N^{-1+epsilon}$ and throughout the entire spectrum away from 0. We also prove analogous results for generalized Wigner matrices.

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.

We develop a new method for deriving local laws for a large class of random matrices. It is applicable to many matrix models built from sums and products of deterministic or independent random matrices. In particular, it may be used to obtain local laws for matrix ensembles that are \emph{anisotropic} in the sense that their resolvents are well approximated by deterministic matrices that are not multiples of the identity. For definiteness, we present the method for sample covariance matrices of the form Q:=TXX∗T∗, where T is deterministic and X is random with independent entries. We prove that with high probability the resolvent of Q is close to a deterministic matrix, with an optimal error bound and down to optimal spectral scales. As an application, we prove the edge universality of Q by establishing the Tracy-Widom-Airy statistics of the eigenvalues of Q near the soft edges. This result applies in the single-cut and multi-cut cases. Further applications include the distribution of the eigenvectors and an analysis of the outliers and BBP-type phase transitions in finite-rank deformations; they will appear elsewhere. We also apply our method to Wigner matrices whose entries have arbitrary expectation, i.e. we consider W+A where W is a Wigner matrix and a Hermitian deterministic matrix. We prove the anisotropic local law for W+A and use it to establish edge universality.