Quantum chaos in supersymmetric quantum mechanics: an exact diagonalization study
We use exact diagonalization to study energy level statistics and out-of-time-order correlators (OTOCs) for the simplest supersymmetric extension $\hat{H}_S = \hat{H}_B \otimes I + \hat{x}_1 \otimes σ_1 + \hat{x}_2 \otimes σ_3$ of the bosonic Hamiltonian $\hat{H}_B = \hat{p}_1^2 + \hat{p}_2^2 + \hat{x}_1^2 \, \hat{x}_2^2$. For a long time, this bosonic Hamiltonian was considered one of the simplest systems which exhibit dynamical chaos both classically and quantum-mechanically. Its structure closely resembles that of spatially compactified pure Yang-Mills theory. Correspondingly, the structure of our supersymmetric Hamiltonian is similar to that of spatially compactified supersymmetric Yang-Mills theory, also known as the BFSS model. We present numerical evidence that a continuous energy spectrum of the supersymmetric model leads to monotonous growth of OTOCs down to the lowest temperatures, which is also expected for the BFSS model from holographic duality. This growth is saturated by low-energy eigenstates with effectively one-dimensional wave functions and a completely non-chaotic energy level distribution. We observe a sharp boundary separating these low-energy states from the bulk of chaotic high-energy states. Our data suggests, although with a limited confidence, that at low temperatures the OTOC growth might be exponential over a finite range of time, with the corresponding Lyapunov exponent scaling linearly with temperature. In contrast, the gapped low-energy spectrum of the bosonic Hamiltonian leads to oscillating OTOCs at low temperatures without any signatures of exponential growth. We also find that the OTOCs for the bosonic Hamiltonian are never sufficiently close to the classical Lyapunov distance. On the other hand, the OTOCs for the supersymmetric system agree with the classical limit reasonably well over a finite range of temperatures and evolution times.