LAUSR

We have recently shown that the logarithmic growth of the entanglement entropy following a quantum quench in a many-body localized phase is accompanied by a slow growth of the number entropy ${S}_{N}\ensuremath{\sim}ln\phantom{\rule{0.16em}{0ex}}ln\phantom{\rule{0.16em}{0ex}}t$. Here we provide an in-depth numerical study of ${S}_{N}(t)$ for the disordered Heisenberg chain and show that this behavior is not transient and persists even for very strong disorder. Calculating the truncated R\'enyi number entropy ${S}_{N}^{(\ensuremath{\alpha})}(t)={(1\ensuremath{-}\ensuremath{\alpha})}^{\ensuremath{-}1}ln{\ensuremath{\sum}}_{n}{p}^{\ensuremath{\alpha}}(n)$ for $\ensuremath{\alpha}\ensuremath{\ll}1$ and $p(n)g{p}_{c}$---which is sensitive to large number fluctuations occurring with low probability---we demonstrate that the particle number distribution $p(n)$ in one half of the system has a continuously growing tail. This indicates a slow but steady increase in the number of particles crossing between the partitions in the interacting case and is in sharp contrast to Anderson localization for which we show that ${S}_{N}^{(\ensuremath{\alpha}\ensuremath{\rightarrow}0)}(t)$ saturates for any cutoff ${p}_{c}g0$. We show, furthermore, that the growth of ${S}_{N}$ is not the consequence of rare states or rare regions but rather represents typical behavior. These findings indicate that the interacting system is never fully localized even for very strong but finite disorder.