Is the localization transition observed for interacting, strongly disordered bosons in two-dimensional optical lattices a genuine dynamical phase transition between a thermalizing and a many-body localized (MBL) phase, or a sharp finite-time and finite-size crossover? This question, posed by the first site-resolved two-dimensional MBL experiment (Choi et al., Science 2016), remains unresolved. Its full resolution requires establishing whether a non-ergodic localized phase survives the thermodynamic and infinite-time limit in two dimensions despite avalanche instabilities, and, if it does, determining the transition's critical properties — dynamical critical exponent, scaling laws, and Griffiths subdiffusive regime — together with a direct eigenstate-level characterization and a quantification of how weak environmental coupling limits the observable lifetime of localization. A decade of theory, numerics, and experiments has produced conflicting evidence in both directions, and no component of this characterization is settled.
Background
Statistical mechanics rests on the expectation that an isolated, generic quantum many-body system serves as its own heat bath: after long evolution, local observables reach thermal values, as formalized by the eigenstate thermalization hypothesis (ETH), which states that individual high-energy eigenstates already encode thermal equilibrium. Strong quenched disorder can break this paradigm. Building on Anderson localization of single particles, Basko, Aleiner and Altshuler (2006) argued that interactions alone do not restore transport at strong disorder: interacting disordered systems can enter a many-body localized (MBL) phase, in which transport is blocked, local memory of an initial state persists forever, and eigenstates violate ETH. In one dimension this picture is supported by extensive numerics, by a phenomenology of local integrals of motion (LIOMs) — quasi-local conserved operators that label the localized phase, with dephasing between them producing a characteristic logarithmic-in-time growth of entanglement entropy — and by ultracold-atom experiments that measure the 'imbalance', i.e. the persistence of an initially prepared density pattern, and even entanglement-related correlations (Lukin et al., 2019). In one dimension, the MBL transition is described by scaling theories in which a diverging length scale is accompanied by a diverging timescale; the associated dynamical critical exponent, governing the power-law divergence of relaxation times as the transition is approached, and subdiffusive Griffiths regimes — slow transport caused by rare, locally weakly disordered regions — are central elements of the one-dimensional theory (Abanin, Altman, Bloch, Serbyn, 2019). Dimensionality is decisive: for noninteracting particles, all states are localized in two dimensions for arbitrarily weak disorder, but interactions dramatically change the problem, and analytical and exact numerical tools are essentially confined to one dimension because the Hilbert space grows exponentially with system size. The experiment of Choi et al. (2016) realized a disordered two-dimensional Bose–Hubbard model (interacting bosons hopping in an optical lattice with a two-dimensional random onsite disorder potential, generated by projecting computer-generated random patterns with a digital mirror device, with a short correlation length set by the imaging numerical aperture) at scales far beyond classical simulation, and observed, through site-resolved imbalance measurements at evolution times of roughly (with the tunneling time), a transition from a thermalizing to a localized regime, with evidence for a diverging length scale and a critical disorder strength that shifts when the interaction energy is reduced. However, the authors emphasized the open issues: an experiment at finite times cannot distinguish a true asymptotic phase transition from a sharp crossover; longer-time and larger-scale data would be needed to extract dynamical critical behavior; eigenstate properties (entanglement growth, LIOMs, site-resolved density correlations, Griffiths effects) were not probed; noise-induced heating, coupling to higher bands, and atom loss bound the accessible times and could eventually relax the quasi-steady state; and higher-dimensional MBL theory was (and remains) incomplete. The problem acquired a sharper theoretical edge with the avalanche mechanism of De Roeck and Huveneers (2017): a rare thermal inclusion in a localized background can, in dimensions , resonantly absorb neighboring localized regions, thermalize them, and thereby grow without bound — an instability that would preclude any true MBL phase in two dimensions, leaving only long-lived transient localization. Subsequent numerical work (Potirniche, Banerjee, Altman, 2019) supported this scenario, while stabilization mechanisms have also been proposed (Foo et al., 2023, showed that Gaussian — super-exponential — localization of LIOMs, natural in the confining potentials of experiments, pushes the avalanche critical dimension from to ), and independent two-dimensional numerics found MBL signatures in strongly disordered two-dimensional models (Bose–Hubbard, Heisenberg, and quantum-Monte-Carlo-amenable constructions) (Wahl, Pal, Simon, 2019; Decker, Kennes, Karrasch, 2022; Swain et al., 2025). Two-dimensional experiments have since measured slow relaxation in quasiperiodic disorder (Bordia et al., 2017), bath-induced delocalization (Rubio-Abadal et al., 2019), and finite-size scaling of the imbalance on optical-lattice and superconducting-simulator platforms (Hur et al., 2025; Li et al., 2025), with random and quasiperiodic disorder showing opposite system-size trends. The question of whether the observed two-dimensional 'MBL' is a phase or a crossover, and what its critical properties are, is therefore a precise and actively contested open problem.
Problem Statement
For interacting quantum many-body systems with strong quenched disorder in two spatial dimensions — as realized by bosons in a disordered two-dimensional optical lattice (the two-dimensional disordered Bose–Hubbard model) — determine whether the observed transition between a thermalizing and a localized regime is a genuine dynamical phase transition, i.e. whether a many-body localized phase with persistent non-ergodic dynamics survives the thermodynamic and infinite-time limit, or whether it is a finite-time, finite-size crossover that eventually thermalizes, as predicted by avalanche-instability arguments based on rare thermal inclusions. This overarching question comprises the following jointly linked components, all originating in the source experiment's outlook: (1) establish whether the quasi-steady localized regime (observed at evolution times of order , with the lattice tunneling time, and at fixed system sizes) remains stable as time and system size grow, including how the answer depends on the disorder statistics (random versus quasiperiodic) and on the interaction strength, and whether the apparent critical disorder drifts without bound; (2) if the transition exists, determine each of its critical properties — the diverging length and timescales, the dynamical critical exponent governing relaxation times near criticality, and the presence and extent of Griffiths-type rare-region effects and subdiffusive transport; (3) if the transition exists, characterize the localized phase at the eigenstate level in two dimensions (logarithmic entanglement-entropy growth, local integrals of motion, and density–density correlations beyond averaged profiles); and (4) quantify how weak coupling to the environment and experimental imperfections (such as heating, low-frequency noise, atom loss, coupling to higher bands) bound the localization lifetimes that any real two-dimensional system can exhibit, separating genuine asymptotic physics from observationally limited transients.
References
- Jae-yoon Choi, Sebastian Hild, Johannes Zeiher, Peter Schauß, Antonio Rubio-Abadal, Tarik Yefsah, Vedika Khemani, David A. Huse, Immanuel Bloch, Christian Gross (2016). Exploring the many-body localization transition in two dimensions. Science 352, 1547–1552. DOI: 10.1126/science.aaf8834, https://doi.org/10.1126/science.aaf8834 (preprint arXiv:1604.04178, https://arxiv.org/abs/1604.04178)
- D. M. Basko, I. L. Aleiner, B. L. Altshuler (2006). Metal–insulator transition in a weakly interacting many-electron system with localized single-particle states. Annals of Physics 321, 1126–1205. DOI: 10.1016/j.aop.2005.11.014, https://doi.org/10.1016/j.aop.2005.11.014
- Dmitry A. Abanin, Ehud Altman, Immanuel Bloch, Maksym Serbyn (2019). Colloquium: Many-body localization, thermalization, and entanglement. Reviews of Modern Physics 91, 021001. DOI: 10.1103/RevModPhys.91.021001, https://doi.org/10.1103/RevModPhys.91.021001 (preprint arXiv:1804.11065, https://arxiv.org/abs/1804.11065)
- Wojciech De Roeck, François Huveneers (2017). Stability and instability towards delocalization in many-body localization systems. Physical Review B 95, 155129. DOI: 10.1103/PhysRevB.95.155129, https://doi.org/10.1103/PhysRevB.95.155129 (preprint arXiv:1608.01815, https://arxiv.org/abs/1608.01815)
- Ionut-Dragos Potirniche, Sumilan Banerjee, Ehud Altman (2019). Exploration of the stability of many-body localization in . Physical Review B 99, 205149. DOI: 10.1103/PhysRevB.99.205149, https://doi.org/10.1103/PhysRevB.99.205149 (preprint arXiv:1805.01475, https://arxiv.org/abs/1805.01475)
- Thorsten B. Wahl, Arijeet Pal, Steven H. Simon (2019). Signatures of the many-body localized regime in two dimensions. Nature Physics 15, 164–169. DOI: 10.1038/s41567-018-0339-x, https://doi.org/10.1038/s41567-018-0339-x (preprint arXiv:1711.02678, https://arxiv.org/abs/1711.02678)
- Kevin S. C. Decker, Dante M. Kennes, Christoph Karrasch (2022). Many-body localization and the area law in two dimensions. Physical Review B 106, L180201. DOI: 10.1103/PhysRevB.106.L180201, https://doi.org/10.1103/PhysRevB.106.L180201 (preprint arXiv:2106.12861, https://arxiv.org/abs/2106.12861)
- Nyayabanta Swain, Ho-Kin Tang, Darryl Chuan Wei Foo, Brian J. J. Khor, Gabriel Lemarié, Fakher F. Assaad, Pinaki Sengupta, Shaffique Adam (2025). Engineering many-body quantum Hamiltonians with nonergodic properties using quantum Monte Carlo. Physical Review B 111, 224201. DOI: 10.1103/PhysRevB.111.224201, https://doi.org/10.1103/PhysRevB.111.224201 (preprint arXiv:2106.08587, https://arxiv.org/abs/2106.08587)
- D. C. W. Foo, N. Swain, P. Sengupta, G. Lemarié, S. Adam (2023). Stabilization mechanism for many-body localization in two dimensions. Physical Review Research 5, L032011. DOI: 10.1103/PhysRevResearch.5.L032011, https://doi.org/10.1103/PhysRevResearch.5.L032011 (preprint arXiv:2202.09072, https://arxiv.org/abs/2202.09072)
- Hyunsoo Ha, Alan Morningstar, David A. Huse (2023). Many-Body Resonances in the Avalanche Instability of Many-Body Localization. Physical Review Letters 130, 250405. DOI: 10.1103/PhysRevLett.130.250405, https://doi.org/10.1103/PhysRevLett.130.250405
- Pranjal Bordia, Henrik Lüschen, Sebastian Scherg, Sarang Gopalakrishnan, Michael Knap, Ulrich Schneider, Immanuel Bloch (2017). Probing slow relaxation and many-body localization in two-dimensional quasiperiodic systems. Physical Review X 7, 041047. DOI: 10.1103/PhysRevX.7.041047, https://doi.org/10.1103/PhysRevX.7.041047 (preprint arXiv:1704.03063, https://arxiv.org/abs/1704.03063)
- Antonio Rubio-Abadal, Jae-yoon Choi, Johannes Zeiher, Simon Hollerith, Jun Rui, Immanuel Bloch, Christian Gross (2019). Many-body delocalization in the presence of a quantum bath. Physical Review X 9, 041014. DOI: 10.1103/PhysRevX.9.041014, https://doi.org/10.1103/PhysRevX.9.041014 (preprint arXiv:1805.00056, https://arxiv.org/abs/1805.00056)
- Alexander Lukin, Matthew Rispoli, Robert Schittko, M. Eric Tai, Adam M. Kaufman, Soonwon Choi, Vedika Khemani, Julian Léonard, Markus Greiner (2019). Probing entanglement in a many-body–localized system. Science 364, 256–260. DOI: 10.1126/science.aau0818, https://doi.org/10.1126/science.aau0818 (preprint arXiv:1805.09819, https://arxiv.org/abs/1805.09819)
- Julian Léonard, Sooshin Kim, Matthew Rispoli, Alexander Lukin, Robert Schittko, Joyce Kwan, Eugene Demler, Dries Sels, Markus Greiner (2023). Probing the onset of quantum avalanches in a many-body localized system. Nature Physics 19, 481–485. DOI: 10.1038/s41567-022-01887-3, https://doi.org/10.1038/s41567-022-01887-3 (preprint arXiv:2012.15270, https://arxiv.org/abs/2012.15270)
- Yunyan Yao, Liang Xiang, Zexian Guo, Zehang Bao, Yong-Feng Yang, Zixuan Song, Haohai Shi, Xuhao Zhu, Feitong Jin, Jiachen Chen, Shibo Xu, Zitian Zhu, Fanhao Shen, Ning Wang, Chuanyu Zhang, Yaozu Wu, Yiren Zou, Pengfei Zhang, Hekang Li, Zhen Wang, Chao Song, Chen Cheng, Rubem Mondaini, H. Wang, J. Q. You, Shi-Yao Zhu, Lei Ying, Qiujiang Guo (2023). Observation of many-body Fock space dynamics in two dimensions. Nature Physics 19, 1459–1465. DOI: 10.1038/s41567-023-02133-0, https://doi.org/10.1038/s41567-023-02133-0 (preprint arXiv:2211.05803, https://arxiv.org/abs/2211.05803)
- J. Šuntajs, J. Bonča, T. Prosen, L. Vidmar (2020). Quantum chaos challenges many-body localization. Physical Review E 102, 062144. DOI: 10.1103/PhysRevE.102.062144, https://doi.org/10.1103/PhysRevE.102.062144 (preprint arXiv:1905.06345, https://arxiv.org/abs/1905.06345)
- Jan Šuntajs, Tomaž Prosen, Lev Vidmar (2023). Localization challenges quantum chaos in the finite two-dimensional Anderson model. Physical Review B 107, 064205. DOI: 10.1103/PhysRevB.107.064205, https://doi.org/10.1103/PhysRevB.107.064205 (preprint arXiv:2212.10625, https://arxiv.org/abs/2212.10625)
- Junhyeok Hur, Joey Li, Byungjin Lee, Kiryang Kwon, Minseok Kim, Samgyu Hwang, Sumin Kim, Yong Soo Yu, Amos Chan, Thorsten Wahl, Jae-yoon Choi (2025). Stability of many-body localization in two dimensions. arXiv:2508.20699, https://arxiv.org/abs/2508.20699
- Tian-Ming Li, Zheng-Hang Sun, Yun-Hao Shi, Zhen-Ting Bao, Yong-Yi Wang, Jia-Chi Zhang, Yu Liu, Cheng-Lin Deng, Yi-Han Yu, Zheng-He Liu, Chi-Tong Chen, Li Li, Hao Li, Hao-Tian Liu, Si-Yun Zhou, Zhen-Yu Peng, Yan-Jun Liu, Ziting Wang, Yue-Shan Xu, Kui Zhao, Yang He, Da'er Feng, Jia-Cheng Song, Cai-Ping Fang, Junrui Deng, Mingyu Xu, Yu-Tao Chen, Bozhen Zhou, Gui-Han Liang, Zhong-Cheng Xiang, Guangming Xue, Dongning Zheng, Kaixuan Huang, Zheng-An Wang, Haifeng Yu, Piotr Sierant, Kai Xu, Heng Fan (2025). Many-body delocalization with a two-dimensional 70-qubit superconducting quantum simulator. arXiv:2507.16882, https://arxiv.org/abs/2507.16882
- C. A. Trugenberger (2025). Infinite symmetry prevents disorder-induced localization in two dimensions. Physical Review B 111, 245101. DOI: 10.1103/PhysRevB.111.245101, https://doi.org/10.1103/PhysRevB.111.245101