A random matrix model towards the quantum chaos transition conjecture

Consider $D$ independent random systems that are modeled by independent $N\times N$ complex Hermitian Wigner matrices. Suppose they are lying on a circle and the neighboring systems interact with each other through a deterministic matrix $A$. We prove that in the asymptotic limit $N\to \infty$, the whole system exhibits a quantum chaos transition when the interaction strength $\|A\|_{HS}$ changes. More precisely, when $\|A\|_{HS}\ge N^{\varepsilon}$, we prove that the bulk eigenvalue statistics match those of a $DN\times DN$ GUE asymptotically and each bulk eigenvector is equally distributed among all $D$ subsystems with probability $1-o(1)$. These phenomena indicate quantum chaos of the whole system. In contrast, when $\|A\|_{HS}\le N^{-\varepsilon}$, we show that the system is integrable: the bulk eigenvalue statistics behave like $D$ independent copies of GUE statistics asymptotically and each bulk eigenvector is localized on only one subsystem. In particular, if we take $D\to \infty$ after the $N\to \infty$ limit, the bulk statistics of the whole system converge to a Poisson point process under the $DN$ scaling.