Spectra of Random Polynomial Matrices: the Petaloid Law

Abstract

We study the distribution of the zeros of $\det P_N(z)$ where $P_N(z)$ is a random monic polynomial matrix, i.e., $P_N(z)=z^{d}I-{\sum }_{j=0}^{d-1}{A}_{j,N}{z}^j$ for possibly coupled random matrices $A_{j,N}$, scaled to have entrywise variance $O(1/N)$. We provide general conditions under which this distribution almost-surely weakly converges as $N\to \infty$ to a deterministic measure, depending on just the variance and covariance of the entries of the $A_j$. This generalizes the circular, elliptic, and semicircle laws, which concern the special case of this question where $d=1$. Unlike those classical laws, these measures can combine nonuniform two-dimensional densities with singular components supported on curves, producing a variety of petal-shaped regions, inspiring our name ``the petaloid law''. We give explicit formulas for the limiting densities and supports. Under a Gaussianity assumption, we also show that there are almost surely no eigenvalues outside small neighborhoods of the limiting support.