The isolated one-dimensional Heisenberg model with static random magnetic fields has become paradigmatic for the analysis of many-body localization. Here, we study the dynamics of this system initially prepared in a highly-excited nonstationary state. Our focus is on the probability for finding the initial state later in time, the so-called survival probability. Two distinct behaviors are identified before equilibration. At short times, the decay is very fast and equivalent to that of clean systems. It subsequently slows down and develops a power-law behavior with an exponent that coincides with the multifractal dimension of the eigenstates.

We study numerically and analytically the quench dynamics of isolated many-body quantum systems. Using full random matrices from the Gaussian orthogonal ensemble, we obtain analytical expressions for the evolution of the survival probability, density imbalance, and out-of-time-ordered correlator. They are compared with numerical results for a one-dimensional-disordered model with two-body interactions and shown to bound the decay rate of this realistic system. Power-law decays are seen at intermediate times, and dips below the infinite time averages (correlation holes) occur at long times for all three quantities when the system exhibits level repulsion. The fact that these features are shared by both the random matrix and the realistic disordered model indicates that they are generic to nonintegrable interacting quantum systems out of equilibrium. Assisted by the random matrix analytical results, we propose expressions that describe extremely well the dynamics of the realistic chaotic system at different time scales.

We study one-dimensional lattices of interacting spins-$\frac{1}{2}$ and show that the effects of quenching the amplitude of a local magnetic field applied to a single site of the lattice can be comparable to the effects of a global perturbation applied instantaneously to the entire system. Both quenches take the system to the chaotic domain, the energy distribution of the initial states approaches a Breit-Wigner shape, the fidelity (Loschmidt echo) decays exponentially, and thermalization becomes viable.

A main feature of a chaotic quantum system is a rigid spectrum where the levels do not cross. We discuss how the presence of level repulsion in lattice many-body quantum systems can be detected from the analysis of their time evolution instead of their energy spectra. This approach is advantageous to experiments that deal with dynamics, but have limited or no direct access to spectroscopy. Dynamical manifestations of avoided crossings occur at long times. They correspond to a drop, referred to as correlation hole, below the asymptotic value of the survival probability and to a bulge above the saturation point of the von Neumann entanglement entropy and the Shannon information entropy. By contrast, the evolution of these quantities at shorter times reflects the level of delocalization of the initial state, but not necessarily a rigid spectrum. The correlation hole is a general indicator of the integrable–chaos transition in disordered and clean models and as such can be used to detect the transition to the many-body localized phase in disordered interacting systems. This article is part of the themed issue ‘Breakdown of ergodicity in quantum systems: from solids to synthetic matter’.