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:2110.05120 · Physical review. B./Physical review. B

Approaching the Thouless energy and Griffiths regime in random spin systems by singular value decomposition

Wen-Jia Rao

Published Feb 22, 2022• 9 citations
arXivHTMLPDF

Abstract

We employ singular value decomposition (SVD) to study the eigenvalue spectra of random spin systems. By SVD, eigenvalue spectrum is decomposed into orthonormal modes ${W}_{k}$ with weight ${\ensuremath{\lambda}}_{k}$. We show that the scree plot (${\ensuremath{\lambda}}_{k}$ with respect to $k$) in the ergodic phase contains two branches that both follow power law ${\ensuremath{\lambda}}_{k}\ensuremath{\sim}{k}^{\ensuremath{-}\ensuremath{\alpha}}$ but with different exponents $\ensuremath{\alpha}$. By evaluating ${W}_{k}$, we verified the part of ${\ensuremath{\lambda}}_{k}$ with $k>{k}_{\text{Th}}$ is universal that follows random matrix theory, where ${k}_{\text{Th}}$ is related to the Thouless energy. We further demonstrate that $\ensuremath{\alpha}$ corresponds only to the exponential part of the level spacing distribution while being insensitive to the level repulsion, or equivalently the system's symmetry. Consequently, $\ensuremath{\alpha}$ gives an underestimation for the many-body localization transition point, which suggests a nonergodic behavior that may be attributed to the Griffiths regime.

DOI
10.1103/physrevb.105.054207

Related papers

Probing the metallic energy spectrum beyond the Thouless energy scale using singular value decomposition

Richard Berkovits

Disordered quantum systems feature an energy scale known as the Thouless energy. For energy ranges below this scale, the properties of the energy spectrum can be described by random matrix theory. Above this scale a different behavior sets in. For a metallic system it was shown long ago by Altshuler and Shklovskii [Sov. Phys. JETP 64, 127 (1986)] that the number variance should increase as a power law with power dependent on only the dimensionality of the system. Although tantalizing hints at this behavior were seen in previous numerical studies, it is quite difficult to verify this prediction using the standard local unfolding methods. Here we use a different unfolding method, i.e., singular value decomposition, and establish a connection between the power law behavior of the scree plot (the singular values ranked by their amplitude) and the power law behavior of the number variance. Thus, we are able to numerically verify Altshuler and Shklovskii's prediction for disordered three-, four-, and five-dimensional single-electron Anderson models on square lattices in the metallic regime. The same method could be applied to systems such as the Sachdev-Ye-Kitaev model and various interacting many-body models for which the many-body localization occurs. It was recently reported that such systems exhibit a Thouless energy, and analyzing the spectrum's behavior on larger scales is of much current interest.

arXiv:2107.04353•Physical review. B./Physical review. B•Aug 23, 2021•17 citations

Large-scale behavior of energy spectra of the quantum random antiferromagnetic Ising chain with mixed transverse and longitudinal fields

Richard Berkovits

In recent years, it became clear that the metallic regime of systems that exhibit a many-body localization (MBL) behavior shows properties that are quite different than the vanilla metallic region of the single-particle Anderson regime. Here we show that the large-scale energy spectrum of a canonical microscopical model featuring MBL displays a nonuniversal behavior at intermediate scales, which is distinct from the deviation from universality seen in the single-particle Anderson regime. The crucial step in revealing this behavior is a global unfolding of the spectrum performed using the singular value decomposition (SVD) which takes into account the sample to sample fluctuations of the spectra. The spectrum properties may be observed directly in the singular value amplitudes via the scree plot, or by using the SVD to unfold the spectra and then perform a number of states variance calculation. Both methods reveal an intermediate scale of energies that follow super-Poissonian statistics.

arXiv:2201.12749•Physical review. B./Physical review. B•Mar 15, 2022•10 citations

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

Rajarshi Pal, Arul Lakshminarayan

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.

arXiv:2002.00682•Feb 3, 2020•0 citations

Signatures of a critical point in the many-body localization transition

Ángel L. Corps, Rafael A. Molina, Armando Relaño

Disordered interacting spin chains that undergo a many-body localization transition are characterized by two limiting behaviors where the dynamics are chaotic and integrable. However, the transition region between them is not fully understood yet. We propose here a possible finite-size precursor of a critical point that shows a typical finite-size scaling and distinguishes between two different dynamical phases. The kurtosis excess of the diagonal fluctuations of the full one-dimensional momentum distribution from its microcanonical average is maximum at this singular point in the paradigmatic disordered J_1 J 1 - J_2 J 2 model. For system sizes accessible to exact diagonalization, both the position and the size of this maximum scale linearly with the system size. Furthermore, we show that this singular point is found at the same disorder strength at which the Thouless and the Heisenberg energies coincide. Below this point, the spectral statistics follow the universal random matrix behavior up to the Thouless energy. Above it, no traces of chaotic behavior remain, and the spectral statistics are well described by a generalized semi-Poissonian model, eventually leading to the integrable Poissonian behavior. We provide, thus, an integrated scenario for the many-body localization transition, conjecturing that the critical point in the thermodynamic limit, if it exists, should be given by this value of disorder strength.

arXiv:2010.08759•SciPost Physics•May 17, 2021•20 citations

Thermalization and Thouless Energy in the Ergodic Side of the Many-Body Localization Transition

Ángel L. Corps, Rafael A. Molina, Armando Relaño

Within the Thouless energy scale, many-body Hamiltonians behave like random matrices. We study the Thouless energy of the disordered Heisenberg quantum spin chain. This model has become the standard model for the study of many-body localization and the transition to an ergodic phase for weak disorder. We focus on the ergodic side of the Many-Body Localization transition. The Thouless energy extracted from long range spectral statistics and from the power spectrum of the fluctuations of the full momentum distribution coincide. This scale behaves non-monotonically and is peaked at an intermediate value of the disorder. An analysis of the scaling behavior of the Thouless energy and its comparison with the number of populated eigenstates after a quench suggests that this scale is not large enough to guarantee thermalization in the thermodynamic limit.

•Jan 12, 2020•1 citation

Statistical and Dynamical Properties in the Energy Spectra of XXZ Spin Chains

Kazue Kudo

Spectral properties of XXZ spin chains are studied intensively. The integrable XXZ spin chain is one of the most important quantum integrable systems, and its Hamiltonian is given by H = J L ∑ j=1 (S j S x j+1 + S y j S y j+1 +∆S z jS z j+1). (1) Here, L is the number of sites; S j = 1 2 σ j and σ α j (α = x, y, z) is the Pauli matrix on the jth site. Equation (1) corresponds to the Heisenberg spin chain at ∆ = 1, to the XY model at ∆ = 0, and to the Ising spin chain when ∆ → ∞. In this thesis, we study on XXZ spin chains with next-nearest-neighbor (NNN) coupling and/or with a magnetic field as well as Eq. (1). First, we discuss level crossings in the energy spectra of the integrable XXZ spin chain. Some of the level crossings are counterexamples of a theorem by von Neumann and Wigner: two energy levels with the same symmetry rarely cross each other in the flow of the spectrum. Multiple degeneracies appear at the particular crossing points where the symmetry of the system is enhanced in the spectral flow. The eigenvalues and the eigenvectors are obtained in an analytical form for 3-site and 4-site chains. The spectral flow of a 6-site chain is shown by numerical calculation. Not only level crossings but also “level bifurcations” are found in the flow of eigenvalues for the inhomogeneous transfer matrix of the XXZ spin chain. The “level bifurcations” are branches such as bifurcation or multiple branches in the flow of eigenvalues and not observed for the homogeneous transfer matrix. The observations are obtained by numerical calculation. We also confirm the appearance of some of the branches through the Bethe ansatz method. The behavior of some branching points are analyzed in the flow of eigenvalues for a particular inhomogeneous case. Level statistics is discussed for the XXZ spin chains with NNN coupling or with a random magnetic field. Evaluating the level-spacing distribution, the spectral rigidity and the number variance, we confirm the following conjecture: the level-spacing distribution should be described by the Poisson distribution for an integrable model and by the Wigner distribution for a nonintegrable model. The XXZ spin chains with the NNN coupling (J2) or with a random magnetic field (h) are nonintegrable, so that mostly the level statistics shows Wigner behavior. However, the level statistics even for a nonintegrable model sometimes shows Poisson-like behavior because of the finite-size effects. For the NNN coupled XXZ spin chain, the behavior of level statistics comes close to the Poissonian around J2 = 0 and/or ∆ = 1. As h/J becomes large, Poisson-like behavior appears for the level statistics of the XXZ spin chain under random magnetic field. When h/J 0, however, Poisson-like behavior appears for small L. In addition to the Poisson-like behavior, unexpected non-Wigner behavior of level statistics can appear due to mixed symmetry. When we discuss level statistics, we have to consider the whole symmetry of a system and perform desymmetrization. The NNN coupled XXZ chains are invariant under spin rotation around the z axis, translation, reflection. However, we have to consider spin reversal symmetry only when S z tot = 0. If spin reversal symmetry is not considered for S z tot = 0, the behavior of level statistics is described by almost the numerical average of Poisson and Wigner behavior. To observe the dynamical properties of the energy spectra, we investigate energy diffusion for NNN coupled XXZ spin chains in the presence of a periodically oscillating magnetic field. Diffusion coefficients (D) are found to obey the power law with respect to both the field strength (B0) and driving frequency (ω): D ∝ (B0ω) . We observe that β = 2 for “linear response regime” and β = 1 for “non-perturbative regime”. In the linear response regime, the diffusion coefficient can be calculated using the Kubo formula. On the other hand, B0ω is so large in the non-perturbative regime that the perturbation theory fails. The diffusion coefficient is smaller than that predicted by the Kubo formula and the perturbation theory fails in the non-perturbative regime. The ranges of the linear response and the non-perturbative regimes depend on J2. When J2 is small, we observe large-amplitude oscillations, which is attributed to a non-diffusive and ballistic nature of the underlying energy diffusion. In this way, we investigate the energy spectra of finite quantum XXZ spin chains in several different aspects. We hope that the results of the thesis will contribute to the development of the physics of low-dimensional systems and mesoscopic systems.

•Jan 1, 2005•2 citations
We employ singular value decomposition (SVD) to study the eigenvalue spectra of random spin systems. By SVD, eigenvalue spectrum is decomposed into orthonormal modes Wk{W}_{k}Wk​ with weight \ensuremathλk{\ensuremath{\lambda}}_{k}\ensuremathλk​. We show that the scree plot (\ensuremathλk{\ensuremath{\lambda}}_{k}\ensuremathλk​ with respect to kkk) in the ergodic phase contains two branches that both follow power law \ensuremathλk\ensuremath∼k\ensuremath−\ensuremathα{\ensuremath{\lambda}}_{k}\ensuremath{\sim}{k}^{\ensuremath{-}\ensuremath{\alpha}}\ensuremathλk​\ensuremath∼k\ensuremath−\ensuremathα but with different exponents \ensuremathα\ensuremath{\alpha}\ensuremathα. By evaluating Wk{W}_{k}Wk​, we verified the part of \ensuremathλk{\ensuremath{\lambda}}_{k}\ensuremathλk​ with k>kThk>{k}_{\text{Th}}k>kTh​ is universal that follows random matrix theory, where kTh{k}_{\text{Th}}kTh​ is related to the Thouless energy. We further demonstrate that \ensuremathα\ensuremath{\alpha}\ensuremathα corresponds only to the exponential part of the level spacing distribution while being insensitive to the level repulsion, or equivalently the system's symmetry. Consequently, \ensuremathα\ensuremath{\alpha}\ensuremathα gives an underestimation for the many-body localization transition point, which suggests a nonergodic behavior that may be attributed to the Griffiths regime.
Disordered quantum systems feature an energy scale known as the Thouless energy. For energy ranges below this scale, the properties of the energy spectrum can be described by random matrix theory. Above this scale a different behavior sets in. For a metallic system it was shown long ago by Altshuler and Shklovskii [Sov. Phys. JETP 64, 127 (1986)] that the number variance should increase as a power law with power dependent on only the dimensionality of the system. Although tantalizing hints at this behavior were seen in previous numerical studies, it is quite difficult to verify this prediction using the standard local unfolding methods. Here we use a different unfolding method, i.e., singular value decomposition, and establish a connection between the power law behavior of the scree plot (the singular values ranked by their amplitude) and the power law behavior of the number variance. Thus, we are able to numerically verify Altshuler and Shklovskii's prediction for disordered three-, four-, and five-dimensional single-electron Anderson models on square lattices in the metallic regime. The same method could be applied to systems such as the Sachdev-Ye-Kitaev model and various interacting many-body models for which the many-body localization occurs. It was recently reported that such systems exhibit a Thouless energy, and analyzing the spectrum's behavior on larger scales is of much current interest.
In recent years, it became clear that the metallic regime of systems that exhibit a many-body localization (MBL) behavior shows properties that are quite different than the vanilla metallic region of the single-particle Anderson regime. Here we show that the large-scale energy spectrum of a canonical microscopical model featuring MBL displays a nonuniversal behavior at intermediate scales, which is distinct from the deviation from universality seen in the single-particle Anderson regime. The crucial step in revealing this behavior is a global unfolding of the spectrum performed using the singular value decomposition (SVD) which takes into account the sample to sample fluctuations of the spectra. The spectrum properties may be observed directly in the singular value amplitudes via the scree plot, or by using the SVD to unfold the spectra and then perform a number of states variance calculation. Both methods reveal an intermediate scale of energies that follow super-Poissonian statistics.
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.
Disordered interacting spin chains that undergo a many-body localization transition are characterized by two limiting behaviors where the dynamics are chaotic and integrable. However, the transition region between them is not fully understood yet. We propose here a possible finite-size precursor of a critical point that shows a typical finite-size scaling and distinguishes between two different dynamical phases. The kurtosis excess of the diagonal fluctuations of the full one-dimensional momentum distribution from its microcanonical average is maximum at this singular point in the paradigmatic disordered J_1 J 1 - J_2 J 2 model. For system sizes accessible to exact diagonalization, both the position and the size of this maximum scale linearly with the system size. Furthermore, we show that this singular point is found at the same disorder strength at which the Thouless and the Heisenberg energies coincide. Below this point, the spectral statistics follow the universal random matrix behavior up to the Thouless energy. Above it, no traces of chaotic behavior remain, and the spectral statistics are well described by a generalized semi-Poissonian model, eventually leading to the integrable Poissonian behavior. We provide, thus, an integrated scenario for the many-body localization transition, conjecturing that the critical point in the thermodynamic limit, if it exists, should be given by this value of disorder strength.
Within the Thouless energy scale, many-body Hamiltonians behave like random matrices. We study the Thouless energy of the disordered Heisenberg quantum spin chain. This model has become the standard model for the study of many-body localization and the transition to an ergodic phase for weak disorder. We focus on the ergodic side of the Many-Body Localization transition. The Thouless energy extracted from long range spectral statistics and from the power spectrum of the fluctuations of the full momentum distribution coincide. This scale behaves non-monotonically and is peaked at an intermediate value of the disorder. An analysis of the scaling behavior of the Thouless energy and its comparison with the number of populated eigenstates after a quench suggests that this scale is not large enough to guarantee thermalization in the thermodynamic limit.
Spectral properties of XXZ spin chains are studied intensively. The integrable XXZ spin chain is one of the most important quantum integrable systems, and its Hamiltonian is given by H = J L ∑ j=1 (S j S x j+1 + S y j S y j+1 +∆S z jS z j+1). (1) Here, L is the number of sites; S j = 1 2 σ j and σ α j (α = x, y, z) is the Pauli matrix on the jth site. Equation (1) corresponds to the Heisenberg spin chain at ∆ = 1, to the XY model at ∆ = 0, and to the Ising spin chain when ∆ → ∞. In this thesis, we study on XXZ spin chains with next-nearest-neighbor (NNN) coupling and/or with a magnetic field as well as Eq. (1). First, we discuss level crossings in the energy spectra of the integrable XXZ spin chain. Some of the level crossings are counterexamples of a theorem by von Neumann and Wigner: two energy levels with the same symmetry rarely cross each other in the flow of the spectrum. Multiple degeneracies appear at the particular crossing points where the symmetry of the system is enhanced in the spectral flow. The eigenvalues and the eigenvectors are obtained in an analytical form for 3-site and 4-site chains. The spectral flow of a 6-site chain is shown by numerical calculation. Not only level crossings but also “level bifurcations” are found in the flow of eigenvalues for the inhomogeneous transfer matrix of the XXZ spin chain. The “level bifurcations” are branches such as bifurcation or multiple branches in the flow of eigenvalues and not observed for the homogeneous transfer matrix. The observations are obtained by numerical calculation. We also confirm the appearance of some of the branches through the Bethe ansatz method. The behavior of some branching points are analyzed in the flow of eigenvalues for a particular inhomogeneous case. Level statistics is discussed for the XXZ spin chains with NNN coupling or with a random magnetic field. Evaluating the level-spacing distribution, the spectral rigidity and the number variance, we confirm the following conjecture: the level-spacing distribution should be described by the Poisson distribution for an integrable model and by the Wigner distribution for a nonintegrable model. The XXZ spin chains with the NNN coupling (J2) or with a random magnetic field (h) are nonintegrable, so that mostly the level statistics shows Wigner behavior. However, the level statistics even for a nonintegrable model sometimes shows Poisson-like behavior because of the finite-size effects. For the NNN coupled XXZ spin chain, the behavior of level statistics comes close to the Poissonian around J2 = 0 and/or ∆ = 1. As h/J becomes large, Poisson-like behavior appears for the level statistics of the XXZ spin chain under random magnetic field. When h/J 0, however, Poisson-like behavior appears for small L. In addition to the Poisson-like behavior, unexpected non-Wigner behavior of level statistics can appear due to mixed symmetry. When we discuss level statistics, we have to consider the whole symmetry of a system and perform desymmetrization. The NNN coupled XXZ chains are invariant under spin rotation around the z axis, translation, reflection. However, we have to consider spin reversal symmetry only when S z tot = 0. If spin reversal symmetry is not considered for S z tot = 0, the behavior of level statistics is described by almost the numerical average of Poisson and Wigner behavior. To observe the dynamical properties of the energy spectra, we investigate energy diffusion for NNN coupled XXZ spin chains in the presence of a periodically oscillating magnetic field. Diffusion coefficients (D) are found to obey the power law with respect to both the field strength (B0) and driving frequency (ω): D ∝ (B0ω) . We observe that β = 2 for “linear response regime” and β = 1 for “non-perturbative regime”. In the linear response regime, the diffusion coefficient can be calculated using the Kubo formula. On the other hand, B0ω is so large in the non-perturbative regime that the perturbation theory fails. The diffusion coefficient is smaller than that predicted by the Kubo formula and the perturbation theory fails in the non-perturbative regime. The ranges of the linear response and the non-perturbative regimes depend on J2. When J2 is small, we observe large-amplitude oscillations, which is attributed to a non-diffusive and ballistic nature of the underlying energy diffusion. In this way, we investigate the energy spectra of finite quantum XXZ spin chains in several different aspects. We hope that the results of the thesis will contribute to the development of the physics of low-dimensional systems and mesoscopic systems.