Local limits in p-adic random matrix theory

We study the distribution of singular numbers of products of certain classes of p-adic random matrices, as both the matrix size and number of products go to infinity simultaneously. In this limit, we prove convergence of the local statistics to a new random point configuration on Z, defined explicitly in terms of certain intricate mixed q-series/exponential sums. This object may be viewed as a nontrivial p-adic analogue of the interpolating distributions of Akemann-Burda-Kieburg, which generalize the sine and Airy kernels and govern limits of complex matrix products. Our proof uses new Macdonald process computations and holds for matrices with iid additive Haar entries, corners of Haar matrices from GL(N)(Z(p)), and the p-adic analogue of Dyson Brownian motion studied by the author (https://arxiv.org/pdf/2309.02865).