Stability of many-body localization in two dimensions

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\times 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.