Exact cospectrality probabilities for uniform random matrices
Alexander Van Werde
TL;DR
The paper derives an exact, non-asymptotic formula for the probability that conjugating a random symmetric matrix by a fixed orthogonal matrix preserves integrality, expressed via Smith ideals of the scaled orthogonal matrix. It then applies the formula to rational cospectrality, obtaining precise low-dimensional results: for $n=2$, $\mathbb{E}[N_2(\ell)]=2^{r+3}/\ell^2$ under a specific prime-factor condition, and for $n=3$, $\mathbb{E}[N_3(\ell)]=48/\ell^2 \prod_{p\mid\ell}(1+1/p)$, with non-monotone dependence on the denominator. The method hinges on Smith normal form invariants over Dedekind domains and a translation-invariant averaging argument modulo ideals, offering a unified approach that extends beyond the specific ensemble. The results yield nontrivial bounds on the probability of rational cospectrality at fixed denominators and illuminate how arithmetic fluctuations influence spectral-conjugacy phenomena in low dimensions, with potential extensions to Hermitian and non-symmetric matrices.
Abstract
We study the conjugation action of orthogonal matrices on symmetric random matrices. Given a fixed orthogonal matrix over an algebraic number field and a random matrix with entries sufficiently uniform in the ring of integers, we wonder what the probability is that the conjugate is again integral. Our main result establishes an exact formula for this probability in terms of the Smith ideals associated to the orthogonal matrix. As an illustrative application, we establish exact formulas for the expected number of rational orthogonal matrices that preserve the integrality of a random matrix for every fixed denominator in dimensions two and three. Notably, the dependence on the denominator turns out to be non-monotone due to number-theoretic fluctuations. We also prove bounds on the probability of rational cospectrality with bounded but arbitrarily large denominator.
