We study the behaviour of two-time correlation functions at late times for finite system sizes considering observables whose (one-point) average value does not depend on energy. In the long time limit, we show that such correlation functions display a ramp and a plateau determined by the correlations of energy levels, similar to what is already known for the spectral form factor. The plateau value is determined, in absence of degenerate energy levels, by the fluctuations of diagonal matrix elements, which highlights differences between different symmetry classes. We show this behaviour analytically by employing results from Random Matrix Theory and the Eigenstate Thermalisation Hypothesis, and numerically by exact diagonalization in the toy example of a Hamiltonian drawn from a Random Matrix ensemble and in a more realistic example of disordered spin glasses at high temperature. Importantly, correlation functions in the ramp regime do not show self-averaging behaviour, and, at difference with the spectral form factor the time average does not coincide with the ensemble average.
We study the behaviour of two-time correlation functions at late times for finite system sizes considering observables whose (one-point) average value does not depend on energy. In the long time limit, we show that such correlation functions display a ramp and a plateau determined by the correlations of energy levels, similar to what is already known for the spectral form factor. The plateau value is determined, in absence of degenerate energy levels, by the fluctuations of diagonal matrix elements, which highlights differences between different symmetry classes. We show this behaviour analytically by employing results from Random Matrix Theory and the Eigenstate Thermalisation Hypothesis, and numerically by exact diagonalization in the toy example of a Hamiltonian drawn from a Random Matrix ensemble and in a more realistic example of disordered spin glasses at high temperature. Importantly, correlation functions in the ramp regime do not show self-averaging behaviour, and, at difference with the spectral form factor the time average does not coincide with the ensemble average.
We numerically study the spectral statistics of open quantum many-body systems (OQMBS) as signatures of quantum chaos (or the lack thereof), using the dissipative spectral form factor (DSFF), a generalization of the spectral form factor to complex spectra. We show that the DSFF of chaotic OQMBS generically displays the quadratic ramp-plateau behaviour of the Ginibre ensemble from random matrix theory, in contrast to the linear ramp-plateau behaviour of the Gaussian ensemble in closed quantum systems. Furthermore, in the presence of many-body interactions, such RMT behaviour emerges only after a time scale τdev, which generally increases with system size for sufficiently large system size, and can be identified as the non-Hermitian analogue of the many-body Thouless time. The universality of the random matrix theory behavior is demonstrated by surveying twelve models of OQMBS, including random Kraus circuits (quantum channels) and random Lindbladians (Liouvillians) in several symmetry classes, as well as Lindbladians of paradigmatic models such as the Sachdev-Ye-Kitaev (SYK), XXZ, and the transverse field Ising models. We devise an unfolding and filtering procedure to remove variations of the averaged density of states which would otherwise hide the universal RMT-like signatures in the DSFF for chaotic OQMBS. Beyond chaotic OQMBS, we study the spectral statistics of non-chaotic OQMBS, specifically the integrable XX model and a system in the many-body localized (MBL) regime in the presence of dissipation, which exhibit DSFF behaviours distinct from the ramp-plateau behaviour of random matrix theory. Lastly, we study the DSFF of Lindbladians with the Hamiltonian term set to zero, i.e. only the jump operators are present, and demonstrate that the results of RMT universality and scaling of many-body Thouless time survive even without coherent evolution.
The fact that generic isolated many-body quantum systems thermalize is understood using the eigenstate thermalization hypothesis (ETH). In recent years, there has been much interest in the behavior of the ETH spectral functions, which characterize the smooth dependence of the variance of the off-diagonal matrix elements of observables on the associated energy and frequency, and whose low-frequency part contains information about the long-time dynamics. In finite systems described by the ETH, the spectral functions are expected to exhibit plateaus below a characteristic frequency ωETH. In this regime, the statistics of the matrix elements of observables are expected to be described by random matrix theory. Related frequencies that have been studied in the literature are ωSFF, which controls the onset of random-matrix behavior in the spectral form factor, and the transport frequencies ωtr, which are derived from transport coefficients. However, a direct quantitative comparison of these frequencies is lacking. Using exact diagonalization, we conduct such a comparison for the spectral functions of current operators in clean and disordered quantum spin ladders with diffusive energy and spin transport. We find clear evidence for the expected low-frequency plateaus in the ETH spectral functions. For the accessible system sizes, ωSFF is consistent with the extent of the plateaus, while the transport frequencies are systematically larger and lie in the nonuniversal regime of the spectral functions. Our findings highlight the need to better understand the origin of these quantitative differences.
A bstract We present a comprehensive analytical study of a variation of the eigenvector ensemble initially proposed by Deutsch [1, 2] for the foundations of the Eigenstate Thermalization Hypothesis (ETH). This ensemble, called the C -ensemble, incorporates additional system-dependent information, enabling the study of complex quantum systems beyond the universal predictions of Random Matrix Theory (RMT). Specifically, we focus on how system-specific details influence late-time behavior in correlation functions, such as the spectral form factor, and how explicit Hamiltonian corrections not captured by RMT can be included. We demonstrate the consistency of this ensemble with respect to the universal (Haar) results by showing that it defines a unitary 1-design for arbitrary systems and for strongly chaotic systems it becomes an approximated 2-design. Universal expressions for two- and four-point ensemble-averaged correlation functions are derived, revealing how system-dependent information is spectrally decoupled. Furthermore, we show that for small energy windows, the correlation functions defined by this ensemble reduce to the predictions made by the ETH.
We investigate the appearance of quantum chaos in a single many-body wave function by analyzing the statistical properties of the eigenvalues of its reduced density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$ of a spatial subsystem A. We find that (i): the spectrum of the density matrix is described by so-called Wishart random matrix theory, which (ii): exhibits besides level repulsion, spectral rigidity, and universal spectral correlations between eigenvalues separated by distances ranging from one up to many mean level spacings, which we investigate. We use these universal spectral characteristics of the reduced density matrix as a definition of chaos in the wave function. A simple and precise characterization of such universal correlations in a spectrum is a segment of strictly linear growth at sufficiently long times, recently called the ``ramp,'' of the spectral form factor which is the Fourier transform of the correlation function between a pair of eigenvalues. It turns out that Wishart and standard random matrix theory have the same universal ``ramp.'' Specifically, here numerical results for the spectral form factor of the density matrix of generic nonintegrable many-body systems, such as one-dimensional quantum Ising and Floquet spin models, are found to exhibit a universal ``ramp'' identical to that appearing for a ``random pure state'' (``Page state,'' or ``Haar state''). The density matrix of the latter is precisely the Wishart random matrix, the reduced density matrix of a completely random wave function. In addition, we study the development of chaos in the wave function by letting an initial direct product state evolve under the unitary time evolution. We find that the universal spectral correlations as manifested by the ``ramp'' set in as soon as the entanglement entropy begins to grow, and first develop for the eigenvalues at the top of the spectrum of the density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$, subsequently spreading over the entire spectrum at later times. Finally, we study a prethermalized regime described by a generalized Gibbs ensemble, which develops in a rapidly driven Floquet model at intermediate times. We find that the prethermalized regime exhibits no chaos, as evidenced by the absence of a ``ramp'' in the spectral form factor of the density matrix, while the universal spectral correlations start to develop when the prethermalized regime finally relaxes at late times to the fully thermalized (infinite temperature) chaotic regime.
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.
Abstract We argue that in a large class of disordered quantum many-body systems, the late time dynamics of time-dependent correlation functions is captured by random matrix theory, specifically the energy eigenvalue statistics of the corresponding ensemble of disordered Hamiltonians. We find that late time correlation functions approximately factorize into a time-dependent piece, which only depends on spectral statistics of the Hamiltonian ensemble, and a time-independent piece, which only depends on the data of the constituent operators of the correlation function. We call this phenomenon “spectral decoupling”, which signifies a dynamical onset of random matrix theory in correlation functions. A key diagnostic of spectral decoupling isk-invariance, which we refine and study in detail. Particular emphasis is placed on the role of symmetries, and connections betweenk-invariance, scrambling, and OTOCs. Disordered Pauli spin systems, as well as the SYK model and its variants, provide a rich source of disordered quantum many-body systems with varied symmetries, and we studyk-invariance in these models with a combination of analytics and numerics.