No Arabic abstract
We study the frequency dependence of the optical conductivity $text{Re} , sigma(omega)$ of the Heisenberg spin-$1/2$ chain in the thermal and near the transition to the many-body localized phase induced by the strength of a random $z$-directed magnetic field. Using the method of dynamical quantum typicality, we calculate the real-time dynamics of the spin-current autocorrelation function and obtain the Fourier transform $text{Re} , sigma(omega)$ for system sizes much larger than accessible to standard exact-diagonalization approaches. We find that the low-frequency behavior of $text{Re} , sigma(omega)$ is well described by $text{Re} , sigma(omega) approx sigma_text{dc} + a , |omega|^alpha$, with $alpha approx 1$ in a wide range within the thermal phase and close to the transition. We particularly detail the decrease of $sigma_text{dc}$ in the thermal phase as a function of increasing disorder for strong exchange anisotropies. We further find that the temperature dependence of $sigma_text{dc}$ is consistent with the existence of a mobility edge.
Many-body localized systems in which interactions and disorder come together defy the expectations of quantum statistical mechanics: In contrast to ergodic systems, they do not thermalize when undergoing nonequilibrium dynamics. What is less clear, however, is how topological features interplay with many-body localized phases as well as the nature of the transition between a topological and a trivial state within the latter. In this work, we numerically address these questions, using a combination of extensive tensor network calculations, specifically DMRG-X, as well as exact diagonalization, leading to a comprehensive characterization of Hamiltonian spectra and eigenstate entanglement properties.
The entanglement spectrum of the reduced density matrix contains information beyond the von Neumann entropy and provides unique insights into exotic orders or critical behavior of quantum systems. Here, we show that strongly disordered systems in the many-body localized phase have power-law entanglement spectra, arising from the presence of extensively many local integrals of motion. The power-law entanglement spectrum distinguishes many-body localized systems from ergodic systems, as well as from ground states of gapped integrable models or free systems in the vicinity of scale-invariant critical points. We confirm our results using large-scale exact diagonalization. In addition, we develop a matrix-product state algorithm which allows us to access the eigenstates of large systems close to the localization transition, and discuss general implications of our results for variational studies of highly excited eigenstates in many-body localized systems.
We numerically study both the avalanche instability and many-body resonances in strongly-disordered spin chains exhibiting many-body localization (MBL). We distinguish between a finite-size/time MBL regime, and the asymptotic MBL phase, and identify some landmarks within the MBL regime. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition, in infinite-length systems, occurs much deeper in the MBL regime than has been suggested by most previous studies. Our other landmarks involve system-wide resonances. We find that the effective matrix elements producing eigenstates with system-wide resonances are enormously broadly distributed. This means that the onset of such resonances in typical samples occurs quite deep in the MBL regime, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regime in to three subregimes: (i) at strongest disorder, typical samples do not have any eigenstates that are involved in system-wide many-body resonances; (ii) there is a substantial intermediate regime where typical samples do have such resonances, but the pair of eigenstates with the minimum spectral gap does not; and (iii) in the weaker randomness regime, the minimum gap is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain non-resonant and the system thus still appears MBL in many respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i).
The many-body localised (MBL) to thermal crossover observed in exact diagonalisation studies remains poorly understood as the accessible system sizes are too small to be in an asymptotic scaling regime. We develop a model of the crossover in short 1D chains in which the MBL phase is destabilised by the formation of many-body resonances. The model reproduces several properties of the numerically observed crossover, including an apparent correlation length exponent $ u=1$, exponential growth of the Thouless time with disorder strength, linear drift of the critical disorder strength with system size, scale-free resonances, apparent $1/omega$ dependence of disorder-averaged spectral functions, and sub-thermal entanglement entropy of small subsystems. In the crossover, resonances induced by a local perturbation are rare at numerically accessible system sizes $L$ which are smaller than a emph{resonance length} $lambda$. For $L gg sqrt{lambda}$, resonances typically overlap, and this model does not describe the asymptotic transition. The model further reproduces controversial numerical observations which Refs. [v{S}untajs et al, 2019] and [Sels & Polkovnikov, 2020] claimed to be inconsistent with MBL. We thus argue that the numerics to date is consistent with a MBL phase in the thermodynamic limit.
The collective and quantum behavior of many-body systems may be harnessed to achieve fast charging of energy storage devices, which have been recently dubbed quantum batteries. In this paper, we present an extensive numerical analysis of energy flow in a quantum battery described by a disordered quantum Ising chain Hamiltonian, whose equilibrium phase diagram presents many-body localized (MBL), Anderson localized (AL), and ergodic phases. We demonstrate that i) the low amount of entanglement of the MBL phase guarantees much better work extraction capabilities than the ergodic phase and ii) interactions suppress temporal energy fluctuations in comparison with those of the non-interacting AL phase. Finally, we show that the statistical distribution of values of the optimal charging time is a clear-cut diagnostic tool of the MBL phase.