Emergent random matrix universality in quantum operator dynamics

The memory function description of quantum operator dynamics involves a carefully chosen split into 'fast' and 'slow' modes. An approximate model for the fast modes can then be used to solve for Green's functions $G(z)$ of the slow modes. One common approach is to form a Krylov operator basis $\{O_m \}_{m\ge 0}$, and then approximate the dynamics restricted to the 'fast space' $\{O_m \}_{m\ge n}$ for $n\gg 1$. Denoting the fast space Green's function by $G_n (z)$, in this work we prove that $G_n (z)$ exhibits universality in the $n \to \infty$ limit, approaching different universal scaling forms in different regions of the complex $z$-plane. This universality of $G_n (z)$ is precisely analogous to the universality of eigenvalue correlations in random matrix theory (RMT). At finite $z$, we show that $G_n (z)$ approaches the Wigner semicircle law, the same as the average bulk resolvent for the Gaussian Unitary Ensemble. When $G(z)$ is the Green's function of certain hydrodynamical variables, we show that at low frequencies $G_n (z)$ is instead governed by the Bessel universality class from RMT. As an application we give a new numerical method, the spectral bootstrap, for approximating spectral functions from Lanczos coefficients. Our proof is complex analytic in nature, involving a Riemann-Hilbert problem which we solve using a steepest-descent-type method, controlled in the $n\to\infty$ limit. This proof technique assumes some regularity conditions on the spectral function, and we comment on their possible connections to the eigenstate thermalization hypothesis. We are led via steepest-descent to a related Coulomb gas ensemble, and we discuss how a recent conjecture--the 'Operator Growth Hypothesis'--is linked to a confinement transition in this Coulomb gas. We then explain the implications of this criticality for the computational resources required to estimate transport coefficients.