Ultra high order cumulants and quantitative CLT for polynomials in Random Matrices
Zhigang Bao, Daniel Munoz George
TL;DR
The paper develops ultra high order cumulant bounds for traces of polynomials in self-adjoint random matrices from GUE/GOE and general Wigner ensembles, allowing deterministic matrices in the polynomials. By introducing quotient-graph combinatorics and genus-type expansions, the authors derive N-dependent bounds that remain meaningful even when r grows with N, revealing N^{2−r} scaling in Gaussian cases and N^{1−r/2} in general Wigner settings. These cumulant bounds fuel three quantitative CLTs for Tr P(X,D): a CLT with Cramér-type corrections, a Berry–Esseen bound, and a concentration inequality capturing Gaussian small-deviation tails and M-degree dependent large-deviation tails, where M is the polynomial degree. The results generalize second-order freeness by providing explicit rates and tail behavior for the full distribution of the trace, including polynomials in both random and deterministic matrices, with strong implications for fluctuations and concentration in non-invariant models. Overall, the work delivers a robust cumulant-based framework to obtain precise probabilistic laws for complex random-matrix polynomials in both Gaussian and non-Gaussian environments, with clear combinatorial and graphical underpinnings.
Abstract
From the study of the high order freeness of random matrices, it is known that the order $r$ cumulant of the trace of a polynomial of $N$-dimensional GUE/GOE is of order $N^{2-r}$ if $r$ is fixed. In this work, we extend the study along three directions. First, we also consider generally distributed Wigner matrices with subexponential entries. Second, we include the deterministic matrices into discussion and consider arbitrary polynomials in random matrices and deterministic matrices. Third, more importantly, we consider the ultra high order cumulants in the sense that $r$ is arbitrary, i.e., could be $N$ dependent. Our main results are the upper bounds of the ultra high order cumulants, for which not only the $N$-dependence but also the $r$-dependence become significant. These results are then used to derive three types of quantitative CLT for the trace of any given self-adjoint polynomial in these random matrix variables: a CLT with a Cramér type correction, a Berry-Esseen bound, and a concentration inequality which captures both the Gaussian tail in the small deviation regime and $M$-dependent tail in the large deviation regime, where $M$ is the degree of the polynomial. In contrast to the second order freeness which implies the CLT for linear eigenvalue statistics of polynomials in random matrices, our study on the ultra high order cumulants leads to the quantitative versions of the CLT.
