We study models of interacting fermions in one dimension to investigate the crossover from integrability to nonintegrability, i.e., quantum chaos, as a function of system size. Using exact diagonalization of finite-sized systems, we study this crossover by obtaining the energy level statistics and Drude weight associated with transport. Our results reinforce the idea that for system size $L\ensuremath{\rightarrow}\ensuremath{\infty}$ nonintegrability sets in for an arbitrarily small integrability-breaking perturbation. The crossover value of the perturbation scales as a power law $\ensuremath{\sim}{L}^{\ensuremath{-}\ensuremath{\eta}}$ when the integrable system is gapless. The exponent $\ensuremath{\eta}\ensuremath{\approx}3$ appears to be robust to microscopic details and the precise form of the perturbation. We conjecture that the exponent in the power law is characteristic of the random matrix ensemble describing the nonintegrable system. For systems with a gap, the crossover scaling appears to be faster than a power law.
We consider the geometry of quantum states associated with random matrix Hamiltonians belonging to ensembles that exhibit an integrable-to-chaotic transition in terms of the nearest-neighbour energy level spacing distribution, focusing on the
β-Gaussian ensembles with generic Dyson index. For the tridiagonal Gaussian
β-ensemble, which shows chaos-to-integrability transition with varying Dyson index, we first calculate an analytical expression for the ensemble-averaged fidelity susceptibility for a two-by-two matrix representation of the Hamiltonian for generic values of the Dyson index, and show that it diverges as the total Hamiltonian, which is the sum of a diagonal matrix with independent elements and a tridiagonal
β-matrix, goes over to the integrable phase. Next, for large-dimensional matrices, we numerically compute the fidelity susceptibility and the quantum metric tensor, respectively, for a one-parameter and a two-parameter class of Hamiltonians that we construct from the
β-ensemble, and find scaling relations of these quantities for chaotic and integrable phases. We also consider variations of the
β-ensemble, such as the one that preserves the rotational invariance of the ensemble, and compute the relevant ensemble-averaged fidelity susceptibility to show that it has similar features as the tridiagonal ensemble in the integrable as well as the chaotic phase, thereby establishing the universality of these properties. Finally, the presence of non-vanishing non-diagonal elements, which arises due to the rotational non-invariance of the
β-ensembles, is a specific feature of the corresponding metric tensor, and we use this to illuminate the difference between the quantum state space geometry of these ensembles and the rotationally invariant ones, such as the classical Gaussian ensembles.
Quantum chaotic systems are conjectured to display a spectrum whose fine-grained features (gaps and correlations) are well described by random matrix theory (RMT). We propose and develop a complementary version of this conjecture: quantum chaotic systems display a Lanczos spectrum whose local means and covariances are well described by RMT. To support this proposal, we first demonstrate its validity in examples of chaotic and integrable systems. We then show that for Haar-random initial states in RMTs the mean and covariance of the Lanczos spectrum suffice to produce the full long-time behavior of general survival probabilities including the spectral form factor, as well as the spread complexity. In addition, for initial states with continuous overlap with energy eigenstates, we analytically find the long-time averages of the probabilities of Krylov basis elements in terms of the mean Lanczos spectrum. This analysis suggests a notion of eigenstate complexity, the statistics of which differentiate integrable systems and classes of quantum chaos. Finally, we clarify the relation between spread complexity and the universality classes of RMT by exploring various values of the Dyson index and Poisson distributed spectra.
We study the transition between integrable and chaotic behavior in dissipative open quantum systems, exemplified by a boundary driven quantum spin chain. The repulsion between the complex eigenvalues of the corresponding Liouville operator in radial distance s is used as a universal measure. The corresponding level spacing distribution is well fitted by that of a static two-dimensional Coulomb gas with harmonic potential at inverse temperature β∈[0,2]. Here, β=0 yields the two-dimensional Poisson distribution, matching the integrable limit of the system, and β=2 equals the distribution obtained from the complex Ginibre ensemble, describing the fully chaotic limit. Our findings generalize the results of Grobe, Haake, and Sommers, who derived a universal cubic level repulsion for small spacings s. We collect mathematical evidence for the universality of the full level spacing distribution in the fully chaotic limit at β=2. It holds for all three Ginibre ensembles of random matrices with independent real, complex, or quaternion matrix elements.
We demonstrate the connection between an operator's matrix element distribution and entangling power via numerical simulations of random, pseudo-random, and quantum chaotic operators. Creating operators with a random distribution of matrix elements is more difficult than creating operators that reproduce other statistical properties of random matrices. Thus, operators that fulfill many random matrix statistical properties may not generate states of high multi-partite entanglement. To quantify the randomness of various statistical distributions and, by extension, entangling power, we use properties of interpolating ensembles that transition between integrable and random matrix ensembles.
Our current understanding of quantum chaos in many-body quantum systems hinges on the random matrix theory (RMT) behavior of eigenstates and their energy level statistics. Although RMT has been remarkably successful in describing "coarse" features of many-body quantum Hamiltonians in chaotic regimes, such as the Wigner-Dyson level spacing statistics or the volume-law behavior of eigenstate entanglement entropy, it remains a challenge to describe their "finer" features, particularly those arising from spatial locality. Here, we show that we can accurately describe the statistical behavior of eigenstate ensembles in many-body Hamiltonians by using pure random states with physical constraints that capture the essential features of the Hamiltonian, specifically spatial locality and symmetries. We demonstrate our approach on local spin Hamiltonians with a scalar U(1) charge. By constructing ensembles of constrained random states that account for two commuting scalar charges playing the role of energy and magnetization, we describe the patterns of entanglement of midspectrum eigenstates beyond their average volume-law behavior, including O(1) corrections and fluctuations, analytically and numerically. When defining the correspondence between quantum chaotic eigenstates in many-body Hamiltonians and RMT ensembles, our work highlights the important role played by spatial locality in describing universal features beyond the volume-law behavior.