Disordered quantum many-body systems pose one of the central challenges in condensed matter physics and quantum information science, as their dynamics are generally intractable for classical computation. Many-body localization (MBL), hypothesized to evade thermalization indefinitely under strong disorder, exemplifies this difficulty. Here, we study the stability of MBL in two dimensions using ultracold atoms in optical lattices with variable system sizes up to 24 24 sites, well beyond the classically simulable regime. Using the imbalance as a probe, we trace the long-time dynamics under two distinctive disorder potentials: quasiperiodic and random disorder. For random disorder, the MBL crossover point shifts to higher disorder strength with increasing system size, consistent with the avalanche scenario. In contrast, with quasiperiodic disorder, we observe no clear system size dependence, suggesting possible stability of MBL in two dimensions.
Hur et al. (Thu,) studied this question.