Quantum chaotic systems exhibit certain universal statistical properties that closely resemble predictions from random matrix theory (RMT). With respect to observables, it has recently been conjectured that, when truncated to a sufficiently narrow energy window, their statistical properties can be described by an unitarily invariant ensemble, and testable criteria have been introduced, which are based on the scaling behavior of free cumulants. In this paper, we investigate the conjecture numerically in translationally invariant Heisenberg spin chains with spin quantum number $S =\frac{1}{2},1,\frac{3}{2}$. Combining a quantum-typicality-based numerical method with the exploitation of the system's symmetries, we study the spin current operator and find clear evidence of consistency with the proposed criteria in chaotic cases. Our findings further support the conjecture of the existence of RMT universality as manifest in the observable properties in quantum chaotic systems.
Random matrix theory universality of current operators in spin-S Heisenberg chains
Quantum chaotic systems exhibit certain universal statistical properties that closely resemble predictions from random matrix theory (RMT). With respect to observables, it has recently been conjectured that, when truncated to a sufficiently narrow energy window, their statistical properties can be described by an unitarily invariant ensemble, and testable criteria have been introduced, which are based on the scaling behavior of free cumulants. In this paper, we investigate the conjecture numerically in translationally invariant Heisenberg spin chains with spin quantum number S=21,1,23. Combining a quantum-typicality-based numerical method with the exploitation of the system's symmetries, we study the spin current operator and find clear evidence of consistency with the proposed criteria in chaotic cases. Our findings further support the conjecture of the existence of RMT universality as manifest in the observable properties in quantum chaotic systems.
A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT).Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges.For single particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc.R. Soc.London, A400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of spectral pair correlation function), from semiclassics has been completed by Müller et al. [Phys.Rev. Lett.93, 014103 (2004)].In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming at the forefront even for simple many-body quantum systems, such as locally interacting spin chains.Such systems display two universal types of behaviour which are termed as the 'many-body localized phase' and 'ergodic phase'.In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder.Here we provide a clear theoretical explanation for these observations.We compute K(t) in the leading two orders in t and show its agreement with RMT for non-integrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin 1/2 models in a periodically kicking transverse field.In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t, hence the leading order RMT behavior K(t)2t is a consequence of translation and reflection symmetry of the Ising partition function.
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.
We study the spectral properties of and spectral crossovers between different random matrix ensembles [Poissonian, Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE)] in correlated spin-chain systems, in the presence of random magnetic fields, and the scalar spin-chirality term, competing with the usual isotropic and time-reversal invariant Heisenberg term. We have investigated these crossovers in the context of the level-spacing distribution and the level-spacing ratio distribution. We use random matrix theory (RMT) analytical results to fit the observed Poissonian-to-GOE and GOE-to-GUE crossovers, and examine the relationship between the RMT crossover parameter $\ensuremath{\lambda}$ and scaled physical parameters of the spin-chain systems in terms of a scaling exponent. We find that the crossover behavior exhibits universality, in the sense that, it becomes independent of lattice size in the large Hamiltonian matrix dimension limit. Moreover, the scaling exponent obtained from such a finite size scaling analysis, seems to be quite robust and independent of the type of crossover considered or the specific spectral correlation measure used.
Random matrix theory (RMT) and semi-classical methods are both used to study the spectra of chaotic quantum systems. The former is usually applied to small distances within the spectrum residing on scales of the mean level spacing where the phenomenon of universality emerges. Through the quantum-classical correspondence provided by semi-classics this regime is, conversely, described by long time dynamical properties of the classical system. In the first half of this thesis we explore such a connection for the study of spectral properties of quantum graphs with dynamically broken time reversal invariance. In physical systems this invariance is often not fully broken and to emulate this we include a rank-1 perturbation on the quantum level resulting in a RMT model outside of the usual universality classes. As we further show, the outcome can depend on graph specific properties and the rank of the perturbation. For the second half of the thesis we consider an opposite limit: instead of long times we study short times, on which universality can not be expected, in a chain-like spin system. In these systems the number of spins takes on a similar role as time, therefore the short time behavior in long chains has remarkable similarities to long time dynamics in few-body system. For instance, the formulation of an evolution operator in spatial direction is possible. Exploring this spatial-time duality we address long range spectral statistic in many-body systems and, for the first time, resolve periodic orbits in a genuine many-body system from the traces of its quantum evolution.
A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT). Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in the spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges. For single-particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc. R. Soc. A 400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of the spectral pair correlation function) from semiclassics has been completed by Müller et al. [Phys. Rev. Lett. 93, 014103 (2004)]. In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming to the forefront even for simple many-body quantum systems, such as locally interacting spin chains. Such systems display two universal types of behaviour which are termed the "many-body localized phase" and "ergodic phase." In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder. Here we provide a clear theoretical explanation for these observations. We compute K(t) in the leading two orders in t and show its agreement with RMT for nonintegrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin-1/2 models in a periodically kicking transverse field. In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t; hence, the leading-order RMT behavior K(t)≃2t is a consequence of translation and reflection symmetry of the Ising partition function.Received 5 February 2018Revised 12 April 2018DOI:https://doi.org/10.1103/PhysRevX.8.021062Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasQuantum chaosQuantum statistical mechanicsStatistical PhysicsNonlinear DynamicsCondensed Matter, Materials & Applied Physics
The high complexity of many-body quantum dynamics means that essentially all analytical or numerical approaches either exploit special structure or are approximate in nature. One such approach—the memory function formalism—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, and the success of this approach depends on the accuracy of the fast space approximation. Using a formulation in operator Krylov space known as the recursion method, we prove the emergence of a universal random matrix description of the fast mode dynamics. This is captured by the “level- n Green’s function” G n ( z ) , which we show approaches universal scaling forms in the “fast limit” n → ∞ . Notably, this emergent universality can occur in both chaotic and nonchaotic systems, provided their spectral functions are sufficiently smooth. This universality of G n ( z ) turns out to be precisely analogous to the universality of eigenvalue correlations in random matrix theory (RMT), even though there is present in the Hamiltonian. Concretely, at finite z we show that G n ( z ) approaches the Wigner semicircle law, while if 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 of this universality, we give a new numerical method, the , for approximating spectral functions, including hydrodynamic transport data, from a finite number of Lanczos coefficients. Our proof involves a map to a Riemann-Hilbert problem which we solve using a steepest-descent-type method, rigorously controlled in the n → ∞ limit. Via the steepest-descent procedure, we are led to a related Coulomb gas optimization problem, and we discuss how a recent conjecture—the “Operator Growth Hypothesis”—implies that chaotic operator dynamics can generically be identified with the critical point of a confinement transition in this Coulomb gas. These results elevate the recursion method from a useful numerical technique to a theoretically principled framework with universal content.