Many-body localization and thermalization in disordered Hubbard chains

We study the many-body localization transition in one-dimensional Hubbard chains using exact diagonalization and quantum chaos indicators. We also study dynamics in the delocalized (ergodic) and localized phases and discuss thermalization and eigenstate thermalization, or the lack thereof, in such systems. Consistently within the indicators and observables studied, we find that ergodicity is very robust against disorder, namely, even in the presence of weak Hubbard interactions the disorder strength needed for the system to localize is large. We show that this robustness might be hidden by finite size effects in experiments with ultracold fermions.