Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2308.07431 · Physical Review Research

Similarity between a many-body quantum avalanche model and the ultrametric random matrix model

Jan Šuntajs, Miroslav Hopjan, Wojciech De Roeck, Lev Vidmar

Published Apr 5, 2024• 15 citations
arXivHTMLPDF

Abstract

In the field of ergodicity-breaking phases, it has been recognized that quantum avalanches can destabilize many-body localization at a wide range of disorder strengths. This has in particular been demonstrated by the numerical study of a toy model, sometimes simply called the “avalanche model” or the “quantum sun model” [], which consists of an ergodic seed coupled to a perfectly localized material. In this paper, we connect this toy model to a well-studied model in random matrix theory, the ultrametric ensemble. We conjecture that the models share the following features. (1) The location of the critical point may be predicted sharply by analytics. (2) On the localized site, both models exhibit Fock space localization. (3) There is a manifold of critical points. On the critical manifold, the eigenvectors exhibit nontrivial multifractal behavior that can be tuned by moving on the manifold. (4) The spectral statistics at criticality is intermediate between Poisson statistics and random matrix statistics, also tunable on the critical manifold. We confirm numerically these properties. Published by the American Physical Society 2024

DOI
10.1103/physrevresearch.6.023030

Related papers

How Random Are Ergodic Eigenstates of the Ultrametric Random Matrices and the Quantum Sun Model?

Tanay Pathak

We numerically study the extreme-value statistics of the Schmidt eigenvalues of reduced density matrices obtained from the ergodic eigenstates. We start by exploring the extreme value statistics of the ultrametric random matrices and then the related Quantum Sun Model, which is also a toy model of avalanche theory. It is expected that these ergodic eigenstates are purely random and thus possess random matrix theory-like features, and the corresponding eigenvalue density should follow the universal Marchenko-Pastur law. Nonetheless, we find deviations, specifically near the tail in both cases. Similarly, the distribution of maximum eigenvalue, after appropriate centering and scaling, should follow the Tracy-Widom distribution. However, our results show that, for both the ultrametric random matrix and the Quantum Sun model, it can be better described using the extreme value distribution. As the extreme value distribution is associated with uncorrelated or weakly correlated random variables, the results hence indicate that the Schmidt eigenvalues exhibit much weaker correlations compared to the strong correlations typically observed in Wishart matrices. Similar deviations are observed for the case of minimum Schmidt eigenvalues as well . Despite the spectral statistics, such as nearest neighbor spacing ratios, aligning with the random matrix theory predictions, our findings reveal that randomness is still not fully achieved. This suggests that deviations in extreme-value statistics offer a stringent test to probe the randomness of ergodic eigenstates and can provide deeper insights into the underlying structure and correlations in ergodic systems.

arXiv:2501.19244•Jan 31, 2025•0 citations

Many-body quantum chaos: analytic connection to random matrix theory

Pavel Kos, Marko Ljubotina, Tomaž Prosen

A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT). Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in the spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges. For single-particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc. R. Soc. A 400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of the spectral pair correlation function) from semiclassics has been completed by Müller et al. [Phys. Rev. Lett. 93, 014103 (2004)]. In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming to the forefront even for simple many-body quantum systems, such as locally interacting spin chains. Such systems display two universal types of behaviour which are termed the "many-body localized phase" and "ergodic phase." In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder. Here we provide a clear theoretical explanation for these observations. We compute K(t) in the leading two orders in t and show its agreement with RMT for nonintegrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin-1/2 models in a periodically kicking transverse field. In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t; hence, the leading-order RMT behavior K(t)≃2t is a consequence of translation and reflection symmetry of the Ising partition function.Received 5 February 2018Revised 12 April 2018DOI:https://doi.org/10.1103/PhysRevX.8.021062Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasQuantum chaosQuantum statistical mechanicsStatistical PhysicsNonlinear DynamicsCondensed Matter, Materials & Applied Physics

•Jun 13, 2018•144 citations

Ergodic inclusions in many-body localized systems

Luis Colmenárez, David J. Luitz, Wojciech De Roeck

We investigate the effect of ergodic inclusions in putative many-body localized systems. We consider the random field Heisenberg chain, which is many-body localized at strong disorder and we couple it to an ergodic bubble, modeled by a random matrix Hamiltonian. Recent theoretical work suggests that localized systems are unstable to ergodic bubbles, driving the delocalization transition. We tentatively confirm this by numerically analyzing the response of the on-site purities to the insertion of the bubble. For a range of intermediate disorder strengths, this response decays very slowly, or not at all, with increasing distance to the bubble. This suggests that at those disorder strengths, the system is actually delocalized in the thermodynamic limit. However, the signal is quite weak and artefacts in the numerics cannot be ruled out conclusively. Published by the American Physical Society 2024

Physical review. B./Physical review. B•Feb 28, 2024•11 citations

Fading ergodicity and quantum dynamics in random matrix ensembles

Rafał Świętek, Maksymilian Kliczkowski, Miroslav Hopjan +1 more

Recent work has proposed fading ergodicity as a mechanism for many-body ergodicity breaking. Here, we show that two paradigmatic random matrix ensembles -- the Rosenzweig-Porter model and the ultrametric model -- fall within the same universality class of ergodicity breaking when embedded in a many-body Hilbert space of spins-1/2. By calibrating the parameters of both models via their Thouless times, we demonstrate that the matrix elements of local observables display similar statistical properties, allowing us to identify the fractal phase of the Rosenzweig-Porter model with the fading-ergodicity regime. This correspondence is further supported through the analyses of quantum-quench dynamics of local observables, their temporal fluctuations and power spectra, and survival probabilities. Our findings reveal that local observables thermalize within the fading-ergodicity regime on timescales shorter than the Heisenberg time, thus providing a unified framework for understanding ergodicity breaking across these distinct models.

cond-mat.stat-mechCross-listedarXiv:2603.23616•Physical review. B./Physical review. B•Mar 26, 2026•1 citation

Many-Body-Localization Transition : strong multifractality spectrum at the critical point towards a non-ergodic delocalized phase

Cécile Monthus

For short-ranged disordered quantum models in one dimension, the Many-Body-Localization is analyzed via the adaptation to the Many-Body context [M. Serbyn, Z. Papic and D.A. Abanin, PRX 5, 041047 (2015)] of the Thouless point of view on the Anderson transition : the question is whether a local interaction between two long chains is able to reshuffle completely the eigenstates (Delocalized phase with a volume-law entanglement) or whether the hybridization between tensor states remains limited (Many-Body-Localized Phase with an area-law entanglement). The central object is thus the level of Hybridization induced by the matrix elements of local operators, as compared with the difference of diagonal energies. The multifractal analysis of these matrix elements of local operators is used to analyze the corresponding statistics of resonances. Our main conclusion is that the critical point is characterized by the Strong-Multifractality Spectrum $f(0 \leq \alpha \leq 2)=\frac{\alpha}{2}$, well known in the context of Anderson Localization in spaces of effective infinite dimensionality, where the size of the Hilbert space grows exponentially with the volume. As a consequence, the level statistics at criticality is the closest possible to the Poisson statistics of the Localized phase, and the Delocalized phase near the critical point is not ergodic.

•Mar 15, 2016•2 citations

Unconventional thermalization of a localized chain interacting with an ergodic bath

Konrad Pawlik, Nicolas Laflorencie, Jakub Zakrzewski

The study of many-body localized (MBL) phases intrinsically links spectral properties with eigenstate characteristics: localized systems exhibit Poisson level statistics and area-law entanglement entropy, while ergodic systems display volume-law entanglement and follow random matrix theory predictions, including level repulsion. Here, we introduce the interacting Anderson Quantum Sun model, which significantly deviates from these conventional expectations. In addition to standard localized and ergodic phases, we identify a regime that exhibits volume-law entanglement coexisting with intermediate spectral statistics. We also identify another nonstandard regime marked by Poisson level statistics, sub-volume entanglement growth, and rare-event-dominated correlations, indicative of emerging ergodic instabilities. These results highlight unconventional routes of ergodicity breaking and offer fresh perspectives on how Anderson localization may be destabilized.

arXiv:2507.18286•Physical review. B./Physical review. B•Apr 13, 2026•2 citations
In the field of ergodicity-breaking phases, it has been recognized that quantum avalanches can destabilize many-body localization at a wide range of disorder strengths. This has in particular been demonstrated by the numerical study of a toy model, sometimes simply called the “avalanche model” or the “quantum sun model” [], which consists of an ergodic seed coupled to a perfectly localized material. In this paper, we connect this toy model to a well-studied model in random matrix theory, the ultrametric ensemble. We conjecture that the models share the following features. (1) The location of the critical point may be predicted sharply by analytics. (2) On the localized site, both models exhibit Fock space localization. (3) There is a manifold of critical points. On the critical manifold, the eigenvectors exhibit nontrivial multifractal behavior that can be tuned by moving on the manifold. (4) The spectral statistics at criticality is intermediate between Poisson statistics and random matrix statistics, also tunable on the critical manifold. We confirm numerically these properties. Published by the American Physical Society 2024
We numerically study the extreme-value statistics of the Schmidt eigenvalues of reduced density matrices obtained from the ergodic eigenstates. We start by exploring the extreme value statistics of the ultrametric random matrices and then the related Quantum Sun Model, which is also a toy model of avalanche theory. It is expected that these ergodic eigenstates are purely random and thus possess random matrix theory-like features, and the corresponding eigenvalue density should follow the universal Marchenko-Pastur law. Nonetheless, we find deviations, specifically near the tail in both cases. Similarly, the distribution of maximum eigenvalue, after appropriate centering and scaling, should follow the Tracy-Widom distribution. However, our results show that, for both the ultrametric random matrix and the Quantum Sun model, it can be better described using the extreme value distribution. As the extreme value distribution is associated with uncorrelated or weakly correlated random variables, the results hence indicate that the Schmidt eigenvalues exhibit much weaker correlations compared to the strong correlations typically observed in Wishart matrices. Similar deviations are observed for the case of minimum Schmidt eigenvalues as well . Despite the spectral statistics, such as nearest neighbor spacing ratios, aligning with the random matrix theory predictions, our findings reveal that randomness is still not fully achieved. This suggests that deviations in extreme-value statistics offer a stringent test to probe the randomness of ergodic eigenstates and can provide deeper insights into the underlying structure and correlations in ergodic systems.
A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT). Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in the spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges. For single-particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc. R. Soc. A 400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of the spectral pair correlation function) from semiclassics has been completed by Müller et al. [Phys. Rev. Lett. 93, 014103 (2004)]. In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming to the forefront even for simple many-body quantum systems, such as locally interacting spin chains. Such systems display two universal types of behaviour which are termed the "many-body localized phase" and "ergodic phase." In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder. Here we provide a clear theoretical explanation for these observations. We compute K(t) in the leading two orders in t and show its agreement with RMT for nonintegrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin-1/2 models in a periodically kicking transverse field. In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t; hence, the leading-order RMT behavior K(t)≃2t is a consequence of translation and reflection symmetry of the Ising partition function.Received 5 February 2018Revised 12 April 2018DOI:https://doi.org/10.1103/PhysRevX.8.021062Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasQuantum chaosQuantum statistical mechanicsStatistical PhysicsNonlinear DynamicsCondensed Matter, Materials & Applied Physics
We investigate the effect of ergodic inclusions in putative many-body localized systems. We consider the random field Heisenberg chain, which is many-body localized at strong disorder and we couple it to an ergodic bubble, modeled by a random matrix Hamiltonian. Recent theoretical work suggests that localized systems are unstable to ergodic bubbles, driving the delocalization transition. We tentatively confirm this by numerically analyzing the response of the on-site purities to the insertion of the bubble. For a range of intermediate disorder strengths, this response decays very slowly, or not at all, with increasing distance to the bubble. This suggests that at those disorder strengths, the system is actually delocalized in the thermodynamic limit. However, the signal is quite weak and artefacts in the numerics cannot be ruled out conclusively. Published by the American Physical Society 2024
Recent work has proposed fading ergodicity as a mechanism for many-body ergodicity breaking. Here, we show that two paradigmatic random matrix ensembles -- the Rosenzweig-Porter model and the ultrametric model -- fall within the same universality class of ergodicity breaking when embedded in a many-body Hilbert space of spins-1/2. By calibrating the parameters of both models via their Thouless times, we demonstrate that the matrix elements of local observables display similar statistical properties, allowing us to identify the fractal phase of the Rosenzweig-Porter model with the fading-ergodicity regime. This correspondence is further supported through the analyses of quantum-quench dynamics of local observables, their temporal fluctuations and power spectra, and survival probabilities. Our findings reveal that local observables thermalize within the fading-ergodicity regime on timescales shorter than the Heisenberg time, thus providing a unified framework for understanding ergodicity breaking across these distinct models.
For short-ranged disordered quantum models in one dimension, the Many-Body-Localization is analyzed via the adaptation to the Many-Body context [M. Serbyn, Z. Papic and D.A. Abanin, PRX 5, 041047 (2015)] of the Thouless point of view on the Anderson transition : the question is whether a local interaction between two long chains is able to reshuffle completely the eigenstates (Delocalized phase with a volume-law entanglement) or whether the hybridization between tensor states remains limited (Many-Body-Localized Phase with an area-law entanglement). The central object is thus the level of Hybridization induced by the matrix elements of local operators, as compared with the difference of diagonal energies. The multifractal analysis of these matrix elements of local operators is used to analyze the corresponding statistics of resonances. Our main conclusion is that the critical point is characterized by the Strong-Multifractality Spectrum f(0≤α≤2)=α2f(0 \leq \alpha \leq 2)=\frac{\alpha}{2}f(0≤α≤2)=2α​, well known in the context of Anderson Localization in spaces of effective infinite dimensionality, where the size of the Hilbert space grows exponentially with the volume. As a consequence, the level statistics at criticality is the closest possible to the Poisson statistics of the Localized phase, and the Delocalized phase near the critical point is not ergodic.
The study of many-body localized (MBL) phases intrinsically links spectral properties with eigenstate characteristics: localized systems exhibit Poisson level statistics and area-law entanglement entropy, while ergodic systems display volume-law entanglement and follow random matrix theory predictions, including level repulsion. Here, we introduce the interacting Anderson Quantum Sun model, which significantly deviates from these conventional expectations. In addition to standard localized and ergodic phases, we identify a regime that exhibits volume-law entanglement coexisting with intermediate spectral statistics. We also identify another nonstandard regime marked by Poisson level statistics, sub-volume entanglement growth, and rare-event-dominated correlations, indicative of emerging ergodic instabilities. These results highlight unconventional routes of ergodicity breaking and offer fresh perspectives on how Anderson localization may be destabilized.