No Arabic abstract
Using numerically exact methods we study transport in an interacting spin chain which for sufficiently strong spatially constant electric field is expected to experience Stark many-body localization. We show that starting from a generic initial state, a spin-excitation remains localized only up to a finite delocalization time, which depends exponentially on the size of the system and the strength of the electric field. This suggests that bona fide Stark many-body localization occurs only in the thermodynamic limit. We also demonstrate that the transient localization in a finite system and for electric fields stronger than the interaction strength can be well approximated by a Magnus expansion up-to times which grow with the electric field strength.
We study the dynamics of an interacting quantum spin chain under the application of a linearly increasing field. This model exhibits a type of localization known as Stark many-body localization. The dynamics shows a strong dependence on the initial conditions, indicating that the system violates the conventional (strong) eigenstate thermalization hypothesis at any finite gradient of the field. This is contrary to reports of a numerically observed ergodic phase. Therefore, the localization is crucially distinct from disorder-driven many-body localization, in agreement with recent predictions on the basis of localization via Hilbert-space shattering.
We theoretically study the response of a many-body localized system to a local quench from a quantum information perspective. We find that the local quench triggers entanglement growth throughout the whole system, giving rise to a logarithmic lightcone. This saturates the modified Lieb-Robinson bound for quantum information propagation in many-body localized systems previously conjectured based on the existence of local integrals of motion. In addition, near the localization-delocalization transition, we find that the final states after the local quench exhibit volume-law entanglement. We also show that the local quench induces a deterministic orthogonality catastrophe for highly excited eigenstates, where the typical wave-function overlap between the pre- and post-quench eigenstates decays {it exponentially} with the system size.
We introduce techniques for analysing the structure of quantum states of many-body localized (MBL) spin chains by identifying correlation clusters from pairwise correlations. These techniques proceed by interpreting pairwise correlations in the state as a weighted graph, which we analyse using an established graph theoretic clustering algorithm. We validate our approach by studying the eigenstates of a disordered XXZ spin chain across the MBL to ergodic transition, as well as the non-equilibrium dyanmics in the MBL phase following a global quantum quench. We successfully reproduce theoretical predictions about the MBL transition obtained from renormalization group schemes. Furthermore, we identify a clear signature of many-body dynamics analogous to the logarithmic growth of entanglement. The techniques that we introduce are computationally inexpensive and in combination with matrix product state methods allow for the study of large scale localized systems. Moreover, the correlation functions we use are directly accessible in a range of experimental settings including cold atoms.
Recent developments in matrix-product-state (MPS) investigations of many-body localization (MBL) are reviewed, with a discussion of benefits and limitations of the method. This approach allows one to explore the physics around the MBL transition in systems much larger than those accessible to exact diagonalization. System sizes and length scales that can be controllably accessed by the MPS approach are comparable to those studied in state-of-the-art experiments. Results for 1D, quasi-1D, and 2D random systems, as well as 1D quasi-periodic systems are presented. On time scales explored (up to $t approx 300$ in units set by the hopping amplitude), a slow, subdiffusive transport in a rather broad disorder range on the ergodic side of the MBL transition is found. For 1D random spin chains, which serve as a standard model of the MBL transition, the MPS study demonstrates a substantial drift of the critical point $W_c(L)$ with the system size $L$: while for $L approx 20$ we find $W_c approx 4$, as also given by exact diagonalization, the MPS results for $L = 50$--100 provide evidence that the critical disorder saturates, in the large-$L$ limit, at $W_c approx 5.5$. For quasi-periodic systems, these finite-size effects are much weaker, which suggests that they can be largely attributed to rare events. For quasi-1D ($dtimes L$, with $d ll L$) and 2D ($Ltimes L$) random systems, the MPS data demonstrate an unbounded growth of $W_c$ in the limit of large $d$ and $L$, in agreement with analytical predictions based on the rare-event avalanche theory.
The Loschmidt echo, defined as the overlap between quantum wave function evolved with different Hamiltonians, quantifies the sensitivity of quantum dynamics to perturbations and is often used as a probe of quantum chaos. In this work we consider the behavior of the Loschmidt echo in the many body localized phase, which is characterized by emergent local integrals of motion, and provides a generic example of non-ergodic dynamics. We demonstrate that the fluctuations of the Loschmidt echo decay as a power law in time in the many-body localized phase, in contrast to the exponential decay in few-body ergodic systems. We consider the spin-echo generalization of the Loschmidt echo, and argue that the corresponding correlation function saturates to a finite value in localized systems. Slow, power-law decay of fluctuations of such spin-echo-type overlap is related to the operator spreading and is present only in the many-body localized phase, but not in a non-interacting Anderson insulator. While most of the previously considered probes of dephasing dynamics could be understood by approximating physical spin operators with local integrals of motion, the Loschmidt echo and its generalizations crucially depend on the full expansion of the physical operators via local integrals of motion operators, as well as operators which flip local integrals of motion. Hence, these probes allow to get insights into the relation between physical operators and local integrals of motion, and access the operator spreading in the many-body localized phase.