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.
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 bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization at least up to time that scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a nonmonotonic dependence of the dynamical exponent on the electric field. As the field increases toward a critical value proportional to the Hubbard interaction strength, transport slows down, becoming subdiffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting superdiffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.
We investigate the ergodicity-to-localization transition in interacting fermion systems subjected to a spatially uniform electric field. For that, we employ the recently proposed equations (TFE), a type of continuous unitary flow equations. This enables us to iteratively determine an approximate diagonal basis of the quantum many-body system. We present improvements to the method, which achieves good accuracy at small to intermediate interaction strengths, even in the absence of an electric field or disorder. Then, we examine two quantities that reveal the fate of Stark MBL in 1D and 2D. First, we investigate the structure of the resulting basis to determine the crossover between ergodic and localized regimes with respect to electric field strength. Second, we simulate long-time dynamics at infinite temperature. Our results in 1D show a localization transition at nonzero field for finite interaction that vanishes with increasing system size leading to localization at infinitesimally small field even in the presence of interactions. In 2D, we find less clear signatures of localization and strong finite-size effects. We establish that the TFE work accurately up to intermediate times but cannot capture higher order effects in interaction strength that lead to delocalization at longer times in finite-size Stark MBL systems.
We consider spinless fermions on a finite one-dimensional lattice, interacting via nearest-neighbor repulsion and subject to a strong electric field. In the noninteracting case, due to Wannier-Stark localization, the single-particle wave functions are exponentially localized even though the model has no quenched disorder. We show that this system remains localized in the presence of interactions and exhibits physics analogous to models of conventional many-body localization (MBL). In particular, the entanglement entropy grows logarithmically with time after a quench, albeit with a slightly different functional form from the MBL case, and the level statistics of the many-body energy spectrum are Poissonian. We moreover predict that a quench experiment starting from a charge-density wave state would show results similar to those of Schreiber et al. [Science 349, 842 (2015)SCIEAS0036-807510.1126/science.aaa7432].
In one-dimensional (1D) disorder-free interacting systems, a sufficiently strong linear potential can induce localization of the many-body eigenstates, a phenomenon dubbed as Stark many-body localization (MBL). In this paper, we investigate the fate of Stark MBL in 1D spinless fermions systems with long-range interactions, specifically focusing on the role of interaction strength. We obtain the Stark MBL phase diagrams by computing the mean gap ratio and many-body inverse participation ratio at half-filling. We show that, for short-range interactions, there is a qualitative symmetry between the limits of weak and strong interactions. However, this symmetry is absent in the case of long-range interactions, where the system is always Stark many-body localized at strong interactions, regardless of the linear potential strength. Furthermore, we study the dynamics of imbalance and entanglement with various initial states using time-dependent variational principle (TDVP) numerical methods. We reveal that the dynamical quantities display a strong dependence on the initial conditions, which suggests that the Hilbert-space fragmentation precludes thermalization. Our results demonstrate the robustness of Stark MBL even in the presence of long-range interactions and offer an avenue to explore MBL in disorder-free systems with long-range interactions.
Novel paradigms of strong ergodicity breaking have recently attracted significant attention in condensed matter physics. Understanding the exact conditions required for their emergence or breakdown not only sheds more light on thermalization and its absence in closed quantum many-body systems, but it also has potential benefits for applications in quantum information technology. A case of particular interest is many-body localization whose conditions are not yet fully settled. Here, we prove that spin chains symmetric under a combination of mirror and spin-flip symmetries and with a non-degenerate spectrum show finite spin transport at zero total magnetization and infinite temperature. We demonstrate this numerically using two prominent examples: the Stark many-body localization system (Stark-MBL) and the symmetrized many-body localization system (symmetrized-MBL). We provide evidence of delocalization at all energy densities and show that delocalization persists when the symmetry is broken. We use our results to construct two localized systems which, when coupled, delocalize each other. Our work demonstrates the dramatic effect symmetries can have on disordered systems, proves that the existence of exact resonances is not a sufficient condition for delocalization, and opens the door to generalization to higher spatial dimensions and different conservation laws.
The principles of ergodicity and thermalization constitute the foundation of statistical mechanics, positing that a many-body system progressively loses its local information as it evolves. Nevertheless, these principles can be disrupted when thermalization dynamics lead to the conservation of local information, as observed in the phenomenon known as many-body localization. Quantum spin chains provide a fundamental platform for exploring the dynamics of closed interacting quantum many-body systems. This study explores the dynamics of a spin chain with S≥1/2 within the Majumdar-Ghosh model, incorporating a non-uniform magnetic field and single-ion anisotropy. Through the use of exact numerical diagonalization, we unveil that a nearly constant-gradient magnetic field suppress thermalization, a phenomenon termed Stark many-body localization (SMBL), previously observed in S=1/2 chains. Furthermore, our findings reveal that the sole presence of single-ion anisotropy is sufficient to prevent thermalization in the system. Interestingly, when the magnitudes of the magnetic field and anisotropy are comparable, they compete, favoring delocalization. Despite the potential hindrance of SMBL by single-ion anisotropy in this scenario, it introduces an alternative mechanism for localization. Our interpretation, considering local energetic constraints and resonances between degenerate eigenstates, not only provides insights into SMBL but also opens avenues for future experimental investigations into the enriched phenomenology of disordered free localized S≥1/2 systems.