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:2002.00682

Probing the randomness of ergodic states: extreme-value statistics in the ergodic and many-body-localized phases

Rajarshi Pal, Arul Lakshminarayan

Published Feb 3, 2020• 0 citations
arXivHTMLPDF

Abstract

The extreme-value statistics of the entanglement spectrum in disordered spin chains possessing a many-body localization transition is examined. It is expected that eigenstates in the metallic or ergodic phase, behave as random states and hence the eigenvalues of the reduced density matrix, commonly referred to as the entanglement spectrum, are expected to follow the eigenvalue statistics of a trace normalized Wishart ensemble. In particular, the density of eigenvalues is supposed to follow the universal Marchenko-Pastur distribution. We find deviations in the tails both for the disordered XXZ with total $S_z$ conserved in the half-filled sector as well as in a model that breaks this conservation. A sensitive measure of deviations is provided by the largest eigenvalue, which in the case of the Wishart ensemble after appropriate shift and scaling follows the universal Tracy-Widom distribution. We show that for the models considered, in the metallic phase, the largest eigenvalue of the reduced density matrix of eigenvector, instead follows the generalized extreme-value statistics bordering on the Fisher-Tipett-Gumbel distribution indicating that the correlations between eigenvalues are much weaker compared to the Wishart ensemble. We show by means of distributions conditional on the high entropy and normalized participation ratio of eigenstates that the conditional entanglement spectrum still follows generalized extreme value distribution. In the deeply localized phase we find heavy tailed distributions and Lévy stable laws in an appropriately scaled function of the largest and second largest eigenvalues. The scaling is motivated by a recently developed perturbation theory of weakly coupled chaotic systems.

DOI
10.48550/arxiv.2002.00682

Related papers

Tracking the many-body localized to ergodic transition via extremal statistics of entanglement eigenvalues

Abhisek Samanta, Kedar Damle, Rajdeep Sensarma

Some interacting disordered many-body systems are unable to thermalize when the quenched disorder becomes larger than a threshold value. Although several properties of nonzero energy density eigenstates (in the middle of the many-body spectrum) exhibit a qualitative change across this many-body localization (MBL) transition, many of the commonly used diagnostics only do so over a broad transition regime. Here we provide evidence that the transition can be located precisely even at modest system sizes by sharply defined changes in the distribution of extremal eigenvalues of the reduced density matrix of subsystems. In particular, our results suggest that $p*={lim}_{{\ensuremath{\lambda}}_{2}\ensuremath{\rightarrow}ln{(2)}^{+}}{P}_{2}({\ensuremath{\lambda}}_{2})$, where ${P}_{2}({\ensuremath{\lambda}}_{2})$ is the probability distribution of the second lowest entanglement eigenvalue ${\ensuremath{\lambda}}_{2}$, behaves as an ``order parameter'' for the MBL phase: $p*>0$ in the MBL phase, while $p*=0$ in the ergodic phase with thermalization. Thus, in the MBL phase, there is a nonzero probability that a subsystem is entangled with the rest of the system only via the entanglement of one subsystem qubit with degrees of freedom outside the region. In contrast, this probability vanishes in the thermal phase.

arXiv:2001.10198•Physical review. B./Physical review. B•Sep 2, 2020•11 citations

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

Statistics of correlation functions in the random Heisenberg chain

Luis A. Colmenarez, Paul McClarty, Masud Haque +1 more

Ergodic quantum many-body systems satisfy the eigenstate thermalization hypothesis (ETH). However, strong disorder can destroy ergodicity through many-body localization (MBL) – at least in one dimensional systems – leading to a clear signal of the MBL transition in the probability distributions of energy eigenstate expectation values of local operators. For a paradigmatic model of MBL, namely the random-field Heisenberg spin chain, we consider the full probability distribution of eigenstate correlation functions across the entire phase diagram. We find gaussian distributions at weak disorder, as predicted by pure ETH. At intermediate disorder – in the thermal phase – we find further evidence for anomalous thermalization in the form of heavy tails of the distributions. In the MBL phase, we observe peculiar features of the correlator distributions: a strong asymmetry in S_i^z S_{i+r}^z S i z S i + r z correlators skewed towards negative values; and a multimodal distribution for spin-flip correlators. A quantitative quasi-degenerate perturbation theory calculation of these correlators yields a surprising agreement of the full distribution with the exact results, revealing, in particular, the origin of the multiple peaks in the spin-flip correlator distribution as arising from the resonant and off-resonant admixture of spin configurations. The distribution of the S_i^zS_{i+r}^z S i z S i + r z correlator exhibits striking differences between the MBL and Anderson insulator cases.

arXiv:1906.10701•SciPost Physics•Nov 21, 2019•34 citations

Extreme value statistics and eigenstate thermalization in kicked quantum chaotic spin-$1/2$ chains

Tanay Pathak, Masaki Tezuka

It is often expected (and assumed) for a quantum chaotic system that the presence of correlated eigenvalues implies that all the other properties as dictated by random matrix theory are satisfied. We demonstrate using the spin-$1/2$ kicked field Ising model that this is not necessarily true. We study the properties of eigenvalues of the reduced density matrix for this model, which constitutes the entanglement spectrum. It is shown that the largest eigenvalue does not follow the expected Tracy--Widom distribution even for the large system sizes considered. The distribution instead follows the extreme value distribution of Weibull type. Furthermore, we also show that such deviations do not lead to drastic change in the thermalization property of this system by showing that the models satisfy the diagonal and off-diagonal eigenstate thermalization hypothesis. Finally, we study the spin-spin autocorrelation function and numerically show that it has the characteristic behavior for chaotic systems: it decreases exponentially and saturates to a value at late time that decreases with system size.

arXiv:2505.16525•May 22, 2025•0 citations

Universal Power-Law Entanglement Spectrum in Many-Body Localized Phases

Maksym Serbyn, Alexios A. Michailidis, Dmitry A. Abanin +1 more

Entanglement spectrum of the reduced density matrix contains information beyond the von Neumann entropy and provides unique insights into exotic orders or critical behavior of electronic and quantum spin systems. Here we show that strongly-disordered systems in the many-body localized phase have universal power-law entanglement spectra, arising from the presence of extensively many local integrals of motion. The power-law entanglement spectrum distinguishes many-body localized systems from ergodic systems, as well as from ground states of gapped integrable models or systems in the vicinity of scale-invariant critical points. We confirm our results using large-scale exact diagonalization. In addition, we develop a matrix-product state algorithm which allows us to access the eigenstates of large systems close to the localization transition, and discuss general implications of our results for variational studies of highly excited eigenstates in many-body localized systems.

•May 18, 2016•0 citations

Statistics of correlation functions in the random Heisenberg chain

Colmenarez, Luis, Paul McClarty, Masudul Haque +1 more

Ergodic quantum many-body systems satisfy the eigenstate thermalization hypothesis (ETH). However, strong disorder can destroy ergodicity through many-body localization (MBL) - at least in one dimensional systems - leading to a clear signal of the MBL transition in the probability distributions of energy eigenstate expectation values of local operators. For a paradigmatic model of MBL, namely the random-field Heisenberg spin chain, we consider the full probability distribution of eigenstate correlation functions across the entire phase diagram. We find gaussian distributions at weak disorder, as predicted by pure ETH. At intermediate disorder - in the thermal phase - we find further evidence for anomalous thermalization in the form of heavy tails of the distributions. In the MBL phase, we observe peculiar features of the correlator distributions: a strong asymmetry in (SiSi+rz)-S-z correlators skewed towards negative values; and a multimodal distribution for spin-flip correlators. A quantitative quasi-degenerate perturbation theory calculation of these correlators yields a surprising agreement of the full distribution with the exact results, revealing, in particular, the origin of the multiple peaks in the spin-flip correlator distribution as arising from the resonant and off-resonant admixture of spin configurations. The distribution of the (SiSi+rz)-S-z correlator exhibits striking differences between the MBL and Anderson insulator cases.

•Jan 1, 2019•33 citations
The extreme-value statistics of the entanglement spectrum in disordered spin chains possessing a many-body localization transition is examined. It is expected that eigenstates in the metallic or ergodic phase, behave as random states and hence the eigenvalues of the reduced density matrix, commonly referred to as the entanglement spectrum, are expected to follow the eigenvalue statistics of a trace normalized Wishart ensemble. In particular, the density of eigenvalues is supposed to follow the universal Marchenko-Pastur distribution. We find deviations in the tails both for the disordered XXZ with total SzS_zSz​ conserved in the half-filled sector as well as in a model that breaks this conservation. A sensitive measure of deviations is provided by the largest eigenvalue, which in the case of the Wishart ensemble after appropriate shift and scaling follows the universal Tracy-Widom distribution. We show that for the models considered, in the metallic phase, the largest eigenvalue of the reduced density matrix of eigenvector, instead follows the generalized extreme-value statistics bordering on the Fisher-Tipett-Gumbel distribution indicating that the correlations between eigenvalues are much weaker compared to the Wishart ensemble. We show by means of distributions conditional on the high entropy and normalized participation ratio of eigenstates that the conditional entanglement spectrum still follows generalized extreme value distribution. In the deeply localized phase we find heavy tailed distributions and Lévy stable laws in an appropriately scaled function of the largest and second largest eigenvalues. The scaling is motivated by a recently developed perturbation theory of weakly coupled chaotic systems.
Some interacting disordered many-body systems are unable to thermalize when the quenched disorder becomes larger than a threshold value. Although several properties of nonzero energy density eigenstates (in the middle of the many-body spectrum) exhibit a qualitative change across this many-body localization (MBL) transition, many of the commonly used diagnostics only do so over a broad transition regime. Here we provide evidence that the transition can be located precisely even at modest system sizes by sharply defined changes in the distribution of extremal eigenvalues of the reduced density matrix of subsystems. In particular, our results suggest that $p*={lim}_{{\ensuremath{\lambda}}_{2}\ensuremath{\rightarrow}ln{(2)}^{+}}{P}_{2}({\ensuremath{\lambda}}_{2})$, where ${P}_{2}({\ensuremath{\lambda}}_{2})$ is the probability distribution of the second lowest entanglement eigenvalue ${\ensuremath{\lambda}}_{2}$, behaves as an ``order parameter'' for the MBL phase: p∗>0p*>0p∗>0 in the MBL phase, while p∗=0p*=0p∗=0 in the ergodic phase with thermalization. Thus, in the MBL phase, there is a nonzero probability that a subsystem is entangled with the rest of the system only via the entanglement of one subsystem qubit with degrees of freedom outside the region. In contrast, this probability vanishes in the thermal phase.
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.
Ergodic quantum many-body systems satisfy the eigenstate thermalization hypothesis (ETH). However, strong disorder can destroy ergodicity through many-body localization (MBL) – at least in one dimensional systems – leading to a clear signal of the MBL transition in the probability distributions of energy eigenstate expectation values of local operators. For a paradigmatic model of MBL, namely the random-field Heisenberg spin chain, we consider the full probability distribution of eigenstate correlation functions across the entire phase diagram. We find gaussian distributions at weak disorder, as predicted by pure ETH. At intermediate disorder – in the thermal phase – we find further evidence for anomalous thermalization in the form of heavy tails of the distributions. In the MBL phase, we observe peculiar features of the correlator distributions: a strong asymmetry in S_i^z S_{i+r}^z S i z S i + r z correlators skewed towards negative values; and a multimodal distribution for spin-flip correlators. A quantitative quasi-degenerate perturbation theory calculation of these correlators yields a surprising agreement of the full distribution with the exact results, revealing, in particular, the origin of the multiple peaks in the spin-flip correlator distribution as arising from the resonant and off-resonant admixture of spin configurations. The distribution of the S_i^zS_{i+r}^z S i z S i + r z correlator exhibits striking differences between the MBL and Anderson insulator cases.
Extreme value statistics and eigenstate thermalization in kicked quantum chaotic spin-1/21/21/2 chains
It is often expected (and assumed) for a quantum chaotic system that the presence of correlated eigenvalues implies that all the other properties as dictated by random matrix theory are satisfied. We demonstrate using the spin-1/21/21/2 kicked field Ising model that this is not necessarily true. We study the properties of eigenvalues of the reduced density matrix for this model, which constitutes the entanglement spectrum. It is shown that the largest eigenvalue does not follow the expected Tracy--Widom distribution even for the large system sizes considered. The distribution instead follows the extreme value distribution of Weibull type. Furthermore, we also show that such deviations do not lead to drastic change in the thermalization property of this system by showing that the models satisfy the diagonal and off-diagonal eigenstate thermalization hypothesis. Finally, we study the spin-spin autocorrelation function and numerically show that it has the characteristic behavior for chaotic systems: it decreases exponentially and saturates to a value at late time that decreases with system size.
Entanglement spectrum of the reduced density matrix contains information beyond the von Neumann entropy and provides unique insights into exotic orders or critical behavior of electronic and quantum spin systems. Here we show that strongly-disordered systems in the many-body localized phase have universal power-law entanglement spectra, arising from the presence of extensively many local integrals of motion. The power-law entanglement spectrum distinguishes many-body localized systems from ergodic systems, as well as from ground states of gapped integrable models or systems in the vicinity of scale-invariant critical points. We confirm our results using large-scale exact diagonalization. In addition, we develop a matrix-product state algorithm which allows us to access the eigenstates of large systems close to the localization transition, and discuss general implications of our results for variational studies of highly excited eigenstates in many-body localized systems.
Ergodic quantum many-body systems satisfy the eigenstate thermalization hypothesis (ETH). However, strong disorder can destroy ergodicity through many-body localization (MBL) - at least in one dimensional systems - leading to a clear signal of the MBL transition in the probability distributions of energy eigenstate expectation values of local operators. For a paradigmatic model of MBL, namely the random-field Heisenberg spin chain, we consider the full probability distribution of eigenstate correlation functions across the entire phase diagram. We find gaussian distributions at weak disorder, as predicted by pure ETH. At intermediate disorder - in the thermal phase - we find further evidence for anomalous thermalization in the form of heavy tails of the distributions. In the MBL phase, we observe peculiar features of the correlator distributions: a strong asymmetry in (SiSi+rz)-S-z correlators skewed towards negative values; and a multimodal distribution for spin-flip correlators. A quantitative quasi-degenerate perturbation theory calculation of these correlators yields a surprising agreement of the full distribution with the exact results, revealing, in particular, the origin of the multiple peaks in the spin-flip correlator distribution as arising from the resonant and off-resonant admixture of spin configurations. The distribution of the (SiSi+rz)-S-z correlator exhibits striking differences between the MBL and Anderson insulator cases.