Edge spectra of Gaussian random symmetric matrices with correlated entries

Debapratim Banerjee,Soumendu Sundar Mukherjee,Dipranjan Pal
2024-09-18
Abstract:We study the largest eigenvalue of a Gaussian random symmetric matrix $X_n$, with zero-mean, unit variance entries satisfying the condition $\sup_{(i, j) \ne (i', j')}|\mathbb{E}[X_{ij} X_{i'j'}]| = O(n^{-(1 + \varepsilon)})$, where $\varepsilon > 0$. It follows from Catalano et al. (2024) that the empirical spectral distribution of $n^{-1/2} X_n$ converges weakly almost surely to the standard semi-circle law. Using a Füredi-Komlós-type high moment analysis, we show that the largest eigenvalue $\lambda_1(n^{-1/2} X_n)$ of $n^{-1/2} X_n$ converges almost surely to $2$. This result is essentially optimal in the sense that one cannot take $\varepsilon = 0$ and still obtain an almost sure limit of $2$. We also derive Gaussian fluctuation results for the largest eigenvalue in the case where the entries have a common non-zero mean. Let $Y_n = X_n + \frac{\lambda}{\sqrt{n}}\mathbf{1} \mathbf{1}^\top$. When $\varepsilon \ge 1$ and $\lambda \gg n^{1/4}$, we show that \[ n^{1/2}\bigg(\lambda_1(n^{-1/2} Y_n) - \lambda - \frac{1}{\lambda}\bigg) \xrightarrow{d} \sqrt{2} Z, \] where $Z$ is a standard Gaussian. On the other hand, when $0 < \varepsilon < 1$, we have $\mathrm{Var}(\frac{1}{n}\sum_{i, j}X_{ij}) = O(n^{1 - \varepsilon})$. Assuming that $\mathrm{Var}(\frac{1}{n}\sum_{i, j} X_{ij}) = \sigma^2 n^{1 - \varepsilon} (1 + o(1))$, if $\lambda \gg n^{\varepsilon/4}$, then we have \[ n^{\varepsilon/2}\bigg(\lambda_1(n^{-1/2} Y_n) - \lambda - \frac{1}{\lambda}\bigg) \xrightarrow{d} \sigma Z. \] While the ranges of $\lambda$ in these fluctuation results are certainly not optimal, a striking aspect is that different scalings are required in the two regimes $0 < \varepsilon < 1$ and $\varepsilon \ge 1$.
Probability,Mathematical Physics,Combinatorics,Statistics Theory
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the behavior of the largest eigenvalue and its fluctuation under different conditions in a Gaussian random symmetric matrix with correlated entries. Specifically, it studies the behavior of the largest eigenvalue \(\lambda_1(n^{-1/2}X_n)\) of the matrix \(n^{-1/2}X_n\) when the entries of the matrix \(X_n\) satisfy zero mean, unit variance, and the correlation between the entries satisfies \(\sup_{(i,j)\neq(i',j')}|E[X_{ij}X_{i'j'}]|=O(n^{-(1+\varepsilon)})\) (\(\varepsilon > 0\)). The main result of the paper shows that under the above conditions, the largest eigenvalue of \(n^{-1/2}X_n\) converges almost surely to 2. In addition, the paper also explores the fluctuation of the largest eigenvalue when the entries have non - zero means and obtains the fluctuation results in different parameter ranges. For example, when \(\varepsilon \geq 1\) and \(\lambda \gg n^{1/4}\), we have: \[ \sqrt{n} \left( \lambda_1(n^{-1/2}Y_n) - \left( \lambda + \frac{1}{\lambda} \right) \right) \overset{d}{\to} \sqrt{2}Z, \] where \(Z\) is a standard Gaussian variable. And when \(0 < \varepsilon < 1\) and \(\lambda \gg n^{\varepsilon/4}\), we have: \[ n^{\varepsilon/2} \left( \lambda_1(n^{-1/2}Y_n) - \left( \lambda + \frac{1}{\lambda} \right) \right) \overset{d}{\to} \sigma Z, \] where \(\sigma\) is a constant related to the variance of the matrix entries. These results reveal the fluctuation characteristics of the largest eigenvalue under different correlation decay rates.