We analyze the dynamics of an initially trapped cloud of interacting quantum particles on a lattice under a linear (Stark) potential. We reveal a dichotomy: initially trapped interacting systems possess features typical of both many-body-localized and thermalizing systems. We consider both fermions (t-V model) and bosons (Bose-Hubbard model). For the zero and infinite interaction limits, both systems are integrable: we provide analytic solutions in terms of the moments of the initial cloud shape and clarify how the recurrent dynamics (many-body Bloch oscillations) depends on the initial state. Away from the integrable points, we identify and explain the timescale at which Bloch oscillations decohere.
We analyze the dynamics of an initially trapped cloud of interacting quantum particles on a lattice under a linear (Stark) potential. We reveal a dichotomy: initially trapped interacting systems possess features typical of both many-body-localized and thermalizing systems. We consider both fermions (t-V model) and bosons (Bose-Hubbard model). For the zero and infinite interaction limits, both systems are integrable: we provide analytic solutions in terms of the moments of the initial cloud shape and clarify how the recurrent dynamics (many-body Bloch oscillations) depends on the initial state. Away from the integrable points, we identify and explain the timescale at which Bloch oscillations decohere.
Variational solutions of the Boltzmann equation usually rely on the concept of linear response. We extend the variational approach for tight-binding models at high entropies to a regime far beyond linear response. We analyze both weakly interacting fermions and incoherent bosons on a lattice. We consider a case where the particles are driven by a constant force, leading to the well-known Bloch oscillations, and we consider interactions that are weak enough not to overdamp these oscillations. This regime is computationally demanding and relevant for ultracold atoms in optical lattices. We derive a simple theory in terms of coupled dynamic equations for the particle density, energy density, current, and heat current, allowing for analytic solutions. As an application, we identify damping coefficients for Bloch oscillations in the Hubbard model at weak interactions and compute them for a one-dimensional toy model. We also approximately solve the long-time dynamics of a weakly interacting, strongly Bloch-oscillating cloud of fermionic particles in a tilted lattice, leading to a subdiffusive scaling exponent.
In this work we demonstrate that nonrandom mechanisms that lead to single-particle localization may also lead to many-body localization, even in the absence of disorder. In particular, we consider interacting spins and fermions in the presence of a linear potential. In the noninteracting limit, these models show the well-known Wannier-Stark localization. We analyze the fate of this localization in the presence of interactions. Remarkably, we find that beyond a critical value of the potential gradient these models exhibit nonergodic behavior as indicated by their spectral and dynamical properties. These models, therefore, constitute a class of generic nonrandom models that fail to thermalize. As such, they suggest new directions for experimentally exploring and understanding the phenomena of many-body localization. We supplement our work by showing that by using machine-learning techniques the level statistics of a system may be calculated without generating and diagonalizing the Hamiltonian, which allows a generation of large statistics.
Many-body localization, an emergent quantum phenomenon, challenges the fundamental principles of thermodynamics and has attracted significant attention. In this study, we demonstrate Stark many-body localization in a one-dimensional Bose gas under the continuous limit, a scenario that extends well beyond the conventional tight-binding model commonly employed in most discussions. This localization occurs within a background harmonic trap without any random disorder, where the potential gradient effectively halts thermalization. Contrary to prior studies on Stark localization, our findings reveal that an increasing gradient does not necessarily lead to greater localization. Instead, we observe a complex transition within this continuous interacting system. Another surprising result is that a higher initial density, which typically correlates with a higher collisional rate and thermalization rate, actually induces stronger localization in this one-dimensional model. Additionally, we observe a density plateau with zero compressibility and low-density fluctuations, which may signal an insulating behavior associated with prethermalization during nonequilibrium dynamics.
The phenomenon of many-body Stark localization of bosons in tilted optical lattice is studied. Despite the fact that no disorder is necessary for Stark localization to occur, it is very similar to well-known many-body localization (MBL) in sufficiently strong disorder. Not only the mean gap ratio reaches Poissonian value as characteristic for localized situations, but also the eigenstates reveal multifractal character as in standard MBL. Stark localization enables a coexistence of spacially separated thermal and localized phases in the harmonic trap similarly to fermions. Stark localization may also lead to spectacular trapping of particles in a reversed harmonic field which naively might be considered as an unstable configuration.
Many-body localization for a system of bosons trapped in a one-dimensional lattice is discussed. Two models that may be realized for cold atoms in optical lattices are considered. The model with a random on-site potential is compared with previously introduced random interactions model. While the origin and character of the disorder in both systems is different they show interesting similar properties. In particular, many-body localization appears for a sufficiently large disorder as verified by a time evolution of initial density wave states as well as using statistical properties of energy levels for small system sizes. Starting with different initial states, we observe that the localization properties are energy-dependent which reveals an inverted many-body localization edge in both systems (that finding is also verified by statistical analysis of energy spectrum). Moreover, we consider computationally challenging regime of transition between many body localized and extended phases where we observe a characteristic algebraic decay of density correlations which may be attributed to subdiffusion (and Griffiths-like regions) in the studied systems. Ergodicity breaking in the disordered Bose–Hubbard models is compared with the slowing-down of the time evolution of the clean system at large interactions.
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.