Random commuting matrices
John E. McCarthy
TL;DR
This work develops a systematic framework to study random $d$-tuples of $n\times n$ matrices constrained by algebraic relations, with a focus on commuting tuples and their eigenvalue distributions. It extends classical random matrix theory to multi-matrix settings by introducing unitarily invariant ensembles on algebraic sets such as ${\frak C}^d_n$ and its Hermitian and non-Hermitian variants, deriving a multidimensional Ginibre-type formula for eigenvalue densities and establishing large-$n$ limits via equilibrium measures for logarithmic potentials. In the Hermitian case, the eigenvalues of suitably scaled random tuples converge to explicit equilibrium measures supported on balls or spheres, yielding projections that follow semicircular laws for low $d$ and concentration near the origin for higher $d$. In the non-Hermitian regime, the irreducible case admits a structured decomposition that leads to explicit densities in low dimensions (notably $n=2$) and reveals attraction phenomena for $d\ge3$, contrasting with repulsion in the Hermitian setting. The results connect random matrix theory with potential theory, free probability, and multivariable operator theory, and open the door to further explorations of spectral laws under algebraic constraints and their implications for related von Neumann algebras and group-theoretic constructions.
Abstract
We define a random commuting $d$-tuple of $n$-by-$n$ matrices to be a random variable that takes values in the set of commuting $d$-tuples and has a distribution that is a rapidly decaying continuous weight on this algebraic set. In the Hermitian case, we characterize the eigenvalue distribution as $n$ tends to infinity. In the non-Hermitian case, we get a formula that holds if the set is irreducible. We show that there are qualitative differences between the single matrix case and the several commuting matrices case.
