Understanding the timescales associated with relaxation to equilibrium in closed quantum many-body systems is one of the central focuses in the study of their nonequilibrium dynamics. At late times, these relaxation processes exhibit universal behavior, emerging from the inherent randomness of chaotic Hamiltonians. In this work, we investigate a disordered spin- 1 / 2 XX ladder—an experimentally realizable model known for its diffusive dynamics—to explore the connection between transport properties and spectral measures derived solely from the energy levels of the system via these relaxation timescales. We begin by analyzing the spectral form factor, which yields the time when the system begins to follow the random matrix theory (RMT) statistics, known as the RMT time. We then determine the Thouless times—the average times for a local excitation to diffuse across the entire finite system—through the linear-response theory for both spin and energy transport. Our numerical results confirm that the RMT time scales quadratically with system size and upper bounds the Thouless times. Interestingly, we also find that, unlike other nonintegrable models, spin diffusion proceeds faster than energy diffusion.
Understanding the timescales associated with relaxation to equilibrium in closed quantum many-body systems is one of the central focuses in the study of their nonequilibrium dynamics. At late times, these relaxation processes exhibit universal behavior, emerging from the inherent randomness of chaotic Hamiltonians. In this work, we investigate a disordered spin- 1 / 2 XX ladder—an experimentally realizable model known for its diffusive dynamics—to explore the connection between transport properties and spectral measures derived solely from the energy levels of the system via these relaxation timescales. We begin by analyzing the spectral form factor, which yields the time when the system begins to follow the random matrix theory (RMT) statistics, known as the RMT time. We then determine the Thouless times—the average times for a local excitation to diffuse across the entire finite system—through the linear-response theory for both spin and energy transport. Our numerical results confirm that the RMT time scales quadratically with system size and upper bounds the Thouless times. Interestingly, we also find that, unlike other nonintegrable models, spin diffusion proceeds faster than energy diffusion.
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.
While the high-temperature spin diffusion in spin chains with random local fields has been the subject of numerous studies concerning the phenomenon of many-body localization (MBL), the energy diffusion in the same models has been much less explored. We demonstrate that energy diffusion is faster at weak random fields but becomes essentially equal at strong fields; hence, both diffusions determine the slowest relaxation time scale (Thouless time) in the system. Numerically reachable finite-size systems reveal the anomalously large distribution of diffusion constants with respect to actual field configurations. Despite the exponential-like dependence of diffusion on field strength, the results for sensitivity to twisted boundary conditions are incompatible with the Thouless criterion for localization and the presumed transition to MBL, at least for numerically reachable sizes. In contrast, we find indications of the scenario of subdiffusive transport, particularly in the dynamical diffusivity response.
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.
Phenomena in quantum chaotic systems such as spectral fluctuations are known to be described remarkably well by the random matrix theory (RMT) in small energy intervals and over large time scales. However, in exchange for analytical tractability, RMT has been constructed to leave out an important element of realistic physical systems { the notion of locality. Since quantum information propagates at a finite speed in a spatially-extended quantum chaotic system, there exists time scales before which RMT is effective. How do we describe the non-equilibrium quantum dynamics at intermediate timescales which are much longer than the local relaxation time but too short for RMT to apply? What are the signatures of many-body quantum chaotic systems beyond the description of RMT? In this thesis, we address these questions by analysing minimal models that incorporate a notion of locality and at the same time retain some analytical tractability of RMT. More specifically, we study Floquet quantum circuits that are random in space but periodic in time. These circuits consist of local and random quantum gates that act on a q-state `spin' at each site. After a brief review of classical chaos, quantum chaos, and quantum circuits in Chapter 1, we present three main results. In Chapter 2, we formulate a diagrammatical approach for computing ensemble averages of observables over random unitary circuits, and derive exact results on spectral statistics and dynamical observables for a Floquet circuit in the large-q limit. In Chapter 3, we further investigate the spectral statistics in spatiallyextended many-body quantum chaotic systems by computing the spectral form factor K(t) analytically and numerically. We show that K(t) follows the RMT description at times longer than the Thouless time tTh and identify the scaling of tTh at the large-q limit of a Floquet circuit. In Chapter 4, we discuss eigenstate correlations for spatially-extended many-body quantum chaotic systems, in terms of the statistical properties of matrix elements of local observables. While the eigenstate thermalization hypothesis (ETH) is known to give an excellent description of these quantities, the phenomenon of scrambling and the butterfly effect imply structure beyond ETH. We determine analytically and verify numerically the universal form of this structure at long distances and small eigenvalue separations for Floquet systems.
A major open question in studies of nonequilibrium quantum dynamics is the identification of the time scales involved in the relaxation process of isolated quantum systems that have many interacting particles. We demonstrate that long time scales can be analytically found by analyzing dynamical manifestations of spectral correlations. Using this approach, we show that the Thouless time tTh and the relaxation time tR increase exponentially with system size. We define tTh as the time at which the spread of the initial state in the many-body Hilbert space is complete and verify that it agrees with the inverse of the Thouless energy. tTh marks the point beyond which the dynamics acquire universal features, while relaxation happens later when the evolution reaches a stationary state. In chaotic systems, ${t}_{\text{Th}}\ensuremath{\ll}{t}_{\text{R}}$, while for systems approaching a many-body localized phase, ${t}_{\text{Th}}\ensuremath{\rightarrow}{t}_{\text{R}}$. Our analytical results for tTh and tR are obtained for the survival probability, which is a global quantity. We show numerically that the same time scales appear also in the evolution of the spin autocorrelation function, which is an experimental local observable. Our studies are carried out for realistic many-body quantum models. The results are compared with those for random matrices.
Quantum transport in disordered systems poses intriguing fundamental questions about the interplay of disorder, interactions, and decoherence, with important implications for nanoscale energy transfer and quantum information transfer. Here, we investigate the emergence of multiple transport timescales in the dissipative dynamics of a spin impurity coupled to a small, spatially disordered network of spins. Using a two-dimensional tight-binding model with dipolar interactions and local dephasing, we demonstrate that geometric heterogeneity leads to hierarchical coupling strengths and pronounced separation of dynamical timescales. By analyzing different metrics for dynamics, we identify distinct relaxation timescales associated with cluster-level equilibration and global equilibration. A minimal three-site model reveals the physical origin of the longest timescale: strong internal hybridization generates an effective detuning that suppresses transfer to other weakly coupled sites, yielding a parametrically enhanced relaxation time in the weak-dephasing regime. We corroborate this picture with nonequilibrium steady-state transport calculations and simulations of disordered spin configurations, demonstrating orders-of-magnitude slowing of relaxation when hierarchical couplings are present. Our results highlight the central role of geometry and connectivity in spin networks and open quantum systems in general, and provide experimentally relevant predictions for relaxation times in small spin baths.