Stability of mobility edges in disordered interacting systems
Many-body localization provides a mechanism to avoid thermalization in isolated interacting quantum systems. The breakdown of thermalization may be complete, when all eigenstates in the many-body spectrum become localized, or partial, when the so-called many-body mobility edge separates localized and delocalized parts of the spectrum. Previously, De Roeck \textit{et al.}[arXiv:1506.01505] suggested a possible instability of the many-body mobility edge in energy density. The local ergodic regions -- so called "bubbles" -- resonantly spread throughout the system, leading to delocalization. In order to study such instability mechanism, in this work we design a model featuring many-body mobility edge in \emph{particle density}: the states at small particle density are localized, while increasing the density of particles leads to delocalization. Using numerical simulations with matrix product states we demonstrate the stability of MBL with respect to small bubbles in large dilute systems for experimentally relevant timescales. In addition, we demonstrate that processes where the bubble spreads are favored over processes that lead to resonant tunneling, suggesting a possible mechanism behind the observed stability of many-body mobility edge. We conclude by proposing experiments to probe particle density mobility edge in Bose-Hubbard model.