We consider a realistic nonequilibrium protocol, where a quantum system in thermal equilibrium is suddenly subjected to an external force. Due to this force, the system is driven out of equilibrium and the expectation values of certain observables acquire a dependence on time. Eventually, upon switching off the external force, the system unitarily evolves under its own Hamiltonian and, as a consequence, the expectation values of observables equilibrate towards specific constant long-time values. Summarizing our main results, we show that, in systems which violate the eigenstate thermalization hypothesis (ETH), this long-time value exhibits an intriguing dependence on the strength of the external force. Specifically, for weak external forces, i.e., within the linear response regime, we show that expectation values thermalize to their original equilibrium values, despite the ETH being violated. In contrast, for stronger perturbations beyond linear response, the quantum system relaxes to some nonthermal value which depends on the previous nonequilibrium protocol. While we present theoretical arguments which underpin these results, we also numerically demonstrate our findings by studying the real-time dynamics of two low-dimensional quantum spin models.

The eigenstate thermalization hypothesis explains the emergence of the thermodynamic equilibrium in isolated quantum many-body systems by assuming a particular structure of the observable's matrix elements in the energy eigenbasis. Schematically, it postulates that off-diagonal matrix elements are random numbers and the observables can be described by random matrix theory (RMT). To what extent a RMT description applies, more precisely at which energy scale matrix elements of physical operators become truly uncorrelated, is, however, not fully understood. We study this issue by introducing a novel numerical approach to probe correlations between matrix elements for Hilbert-space dimensions beyond those accessible by exact diagonalization. Our analysis is based on the evaluation of higher moments of operator submatrices, defined within energy windows of varying width. Considering nonintegrable quantum spin chains, we observe that matrix elements remain correlated even for narrow energy windows corresponding to timescales of the order of thermalization time of the respective observables. We also demonstrate that such residual correlations between matrix elements are reflected in the dynamics of outof-time-ordered correlation functions.

No abstract