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:2401.08031 · Physical Review Research

Operator growth in disordered spin chains: Indications for the absence of many-body localization

A. Weisse, Robert Gerstner, Jesko Sirker

Published Jun 11, 2025• 5 citations
arXivHTMLPDF

Abstract

We consider the spreading of a local operator A in one-dimensional many-body systems with Hamiltonian H by calculating the k -fold commutator [ H , [ H , [ . . . , [ H , A ] ] ] ] . We derive bounds for the operator norm of this commutator in free and interacting systems with and without disorder thus directly connecting the with questions of localization. We analytically show, in particular, that an almost-factorial growth of the operator norm—as recently proven for the random Ising model and strongly suggested here to also hold for the Heisenberg model with random fields—is inconsistent with an exponential localization of A . Assuming there exists a quasilocal unitary U which maps H onto an effective Hamiltonian H ̃ = U H U † = ∑ n E n τ n z + ∑ i , j J i j τ i z τ j z + ⋯ , we show that A ̃ = U A U † is a quasilocal operator which, in contrast to the Anderson case, indeed in the general many-body case, leading to an almost-factorial growth of the commutator norm. Therefore, either the unitary U in many-body systems with maximal norm growth does not exist, and such systems are always ergodic, or unusual nonergodic phases described by H ̃ exist which violate the operator growth hypothesis and in which local operators spread over the entire lattice, implying that transport will eventually set in. To investigate this issue further, we concentrate on the chain with random magnetic fields. We analytically and symbolically verify our general results for the noninteracting Anderson and Aubry-André models. For the case, the symbolic calculations are consistent with a maximal norm growth. Furthermore, we find no indication of a weakened exponential localization of A , expected for strong disorder and low commutator orders if the unitary U does exist. Finally, we study the differences between the interacting and noninteracting cases when trying to perturbatively construct U by consecutive Schrieffer-Wolff transformations. While it is straightforward to show that this construction converges in the Anderson case, we find no indications for a convergence in the interacting case, suggesting that U does not exist and that many-body localization is absent.

DOI
10.1103/wgss-nt8t

Related papers

Length-resolved Operator Growth and Path-Entropy Obstructions to Many-Body Localization

Jesko Sirker

For the disordered Ising chain with transverse and longitudinal fields, where couplings and fields are drawn from strictly positive distributions, Cao~\cite{Cao} has shown that the moments $μ_{2k} = \|[H,σ^z_0]^{(k)}\|_2^2$ grow almost factorially, $μ_{2k}^{1/(2k)}\sim k/\ln k$, and thus asymptotically at the maximal allowed rate. We generalize this result by resolving the operator norm in support length and show that the weight at length $\ell_k \sim k/\ln k$ already exhibits almost factorial growth, $\|[H,σ^z_0]^{(k)}_{\ell_k}\|_2 \gtrsim (k/\ln k)^k$. This implies maximal spatial delocalization of local operators and, in particular, rules out dynamical locality---the strongest form of many-body localization---at any disorder strength. We further establish rigorously a finite-size crossover scale $L\sim (W/J)^2$, where $W$ is the disorder and $J$ the coupling strength. For $L\lesssim (W/J)^2$ numerical studies only access a pre-asymptotic regime. Finally, we identify a structural path-entropy obstruction to perturbative LIOM constructions, based on the almost factorial branching of operator content and independent of resonance effects; the same mechanism strongly suggests ballistic real-time operator spreading, so sub-ballistic or localized dynamics would require a presently unidentified cancellation principle acting on almost factorially many disorder-dependent paths with random amplitudes.

arXiv:2606.07774•Jun 5, 2026•0 citations

Explicit construction of local conserved operators in disordered many-body systems

Thomas E. O’Brien, Dmitry A. Abanin, Guifré Vidal +1 more

The presence and character of local integrals of motion---quasilocal operators that commute with the Hamiltonian---encode valuable information about the dynamics of a quantum system. In particular, strongly disordered many-body systems can generically avoid thermalization when there are extensively many such operators. In this work, we explicitly construct local conserved operators in one-dimensional spin chains by directly minimizing their commutator with the Hamiltonian. We demonstrate the existence of an extensively large set of local integrals of motion in the many-body localized phase of the disordered XXZ spin chain. These operators are shown to have exponentially decaying tails, in contrast to the ergodic phase where the decay is (at best) polynomial in the size of the subsystem. We study the algebraic properties of localized operators and confirm that in the many-body localized phase, they are well described by ``dressed'' spin operators.

arXiv:1608.03296•Physical review. B./Physical review. B•Oct 24, 2016•58 citations

Typical growth behavior of the out-of-time-ordered commutator in many-body localized systems

Juhee Lee, Dongkyu Kim, DH Kim

We investigate the typicality of the growth behavior of the out-of-time-ordered commutator (OTOC) in the many-body localized (MBL) quantum spin chains across random disorder realizations. In the MBL phase of the Heisenberg XXZ chain, we find that the estimate of the OTOC fluctuates significantly with the disorder realizations at the intermediate times of the main growth. Despite the consequent failure of the disorder average in the MBL phase, we argue that the characteristic behavior of the OTOC growth can still be identified by going through individual disorder realizations. We find that a power-law-type growth behavior appears typically after a disorder-dependent relaxation period, which is very close to the ${t}^{2}$ form derived in the effective Hamiltonian of a fully MBL system. The characteristic growth behavior observed at an individual disorder realization is robust in our tests with various state preparations and also verified in another MBL system of the random-transverse-field quantum Ising chain in a uniform longitudinal field.

arXiv:1812.00357•Physical review. B./Physical review. B•May 31, 2019•31 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

Breaking the chains: Extreme value statistics and localization in random spin chains

Jeanne Colbois, Nicolas Laflorencie

Despite a very good understanding of single-particle Anderson localization in one-dimensional (1D) disordered systems, many-body effects are still full of surprises; a famous example is the interaction-driven many-body localization (MBL) problem, about which much has been written, and perhaps the best is yet to come. Interestingly enough the noninteracting limit provides a natural playground to study nontrivial multiparticle physics, offering the possibility to test some general mechanisms with very large-scale exact diagonalization simulations. In this paper, we first revisit the 1D many-body Anderson insulator through the lens of extreme value theory, focusing on the extreme polarizations of the equivalent spin chain model in a random magnetic field. A many-body-induced chain breaking mechanism is explored numerically, and compared to an analytically solvable toy model. A unified description, from weak to large disorder strengths $W$ emerges, where the disorder-dependent average localization length $\ensuremath{\xi}(W)$ governs the extreme events leading to chain breaks. In particular, tails of the local magnetization distributions are controlled by $\ensuremath{\xi}(W)$. Remarkably, we also obtain a quantitative understanding of the full distribution of the extreme polarizations, which is given by a Fr\'echet-type law. In a second part, we explore finite interaction physics and the MBL question. For the available system sizes, we numerically quantify the difference in the extreme value distributions between the interacting problem and the noninteracting Anderson case. Strikingly, we observe a sharp ``extreme-statistics transition'' as $W$ changes, which may coincide with the MBL transition.

arXiv:2305.10574•Physical review. B./Physical review. B•Oct 24, 2023•10 citations

Strongly disordered Anderson insulator chains with generic two-body interaction

Bartosz Krajewski, Lev Vidmar, J. Bonča +1 more

The random-field spin-1/2 XXZ chains, and the corresponding Anderson insulators of spinless fermions with density-density interaction, have been intensively studied in the context of many-body localization. However, we recently argued [B. Krajewski et al., Phys. Rev. Lett. 129, 260601 (2022)] that the two-body density-density interaction in these models is not generic since only a small fraction of this interaction represents a true local perturbation to the Anderson insulator. Here we study ergodicity of strongly disordered Anderson insulator chains, choosing other forms of the two-body interaction for which the strength of the true perturbation is of the same order of magnitude as the bare two-body interaction. Focusing on the strong-interaction regime, numerical results for the level statistics and the eigenstate thermalization hypothesis are consistent with emergence of ergodicity at arbitrary strong disorder.

arXiv:2306.14613•Physical review. B./Physical review. B•Aug 4, 2023•7 citations
We consider the spreading of a local operator A in one-dimensional many-body systems with Hamiltonian H by calculating the k -fold commutator [ H , [ H , [ . . . , [ H , A ] ] ] ] . We derive bounds for the operator norm of this commutator in free and interacting systems with and without disorder thus directly connecting the with questions of localization. We analytically show, in particular, that an almost-factorial growth of the operator norm—as recently proven for the random Ising model and strongly suggested here to also hold for the Heisenberg model with random fields—is inconsistent with an exponential localization of A . Assuming there exists a quasilocal unitary U which maps H onto an effective Hamiltonian H ̃ = U H U † = ∑ n E n τ n z + ∑ i , j J i j τ i z τ j z + ⋯ , we show that A ̃ = U A U † is a quasilocal operator which, in contrast to the Anderson case, indeed in the general many-body case, leading to an almost-factorial growth of the commutator norm. Therefore, either the unitary U in many-body systems with maximal norm growth does not exist, and such systems are always ergodic, or unusual nonergodic phases described by H ̃ exist which violate the operator growth hypothesis and in which local operators spread over the entire lattice, implying that transport will eventually set in. To investigate this issue further, we concentrate on the chain with random magnetic fields. We analytically and symbolically verify our general results for the noninteracting Anderson and Aubry-André models. For the case, the symbolic calculations are consistent with a maximal norm growth. Furthermore, we find no indication of a weakened exponential localization of A , expected for strong disorder and low commutator orders if the unitary U does exist. Finally, we study the differences between the interacting and noninteracting cases when trying to perturbatively construct U by consecutive Schrieffer-Wolff transformations. While it is straightforward to show that this construction converges in the Anderson case, we find no indications for a convergence in the interacting case, suggesting that U does not exist and that many-body localization is absent.
For the disordered Ising chain with transverse and longitudinal fields, where couplings and fields are drawn from strictly positive distributions, Cao~\cite{Cao} has shown that the moments μ2k=∥[H,σ0z](k)∥22μ_{2k} = \|[H,σ^z_0]^{(k)}\|_2^2μ2k​=∥[H,σ0z​](k)∥22​ grow almost factorially, μ2k1/(2k)∼k/ln⁡kμ_{2k}^{1/(2k)}\sim k/\ln kμ2k1/(2k)​∼k/lnk, and thus asymptotically at the maximal allowed rate. We generalize this result by resolving the operator norm in support length and show that the weight at length ℓk∼k/ln⁡k\ell_k \sim k/\ln kℓk​∼k/lnk already exhibits almost factorial growth, ∥[H,σ0z]ℓk(k)∥2≳(k/ln⁡k)k\|[H,σ^z_0]^{(k)}_{\ell_k}\|_2 \gtrsim (k/\ln k)^k∥[H,σ0z​]ℓk​(k)​∥2​≳(k/lnk)k. This implies maximal spatial delocalization of local operators and, in particular, rules out dynamical locality---the strongest form of many-body localization---at any disorder strength. We further establish rigorously a finite-size crossover scale L∼(W/J)2L\sim (W/J)^2L∼(W/J)2, where WWW is the disorder and JJJ the coupling strength. For L≲(W/J)2L\lesssim (W/J)^2L≲(W/J)2 numerical studies only access a pre-asymptotic regime. Finally, we identify a structural path-entropy obstruction to perturbative LIOM constructions, based on the almost factorial branching of operator content and independent of resonance effects; the same mechanism strongly suggests ballistic real-time operator spreading, so sub-ballistic or localized dynamics would require a presently unidentified cancellation principle acting on almost factorially many disorder-dependent paths with random amplitudes.
The presence and character of local integrals of motion---quasilocal operators that commute with the Hamiltonian---encode valuable information about the dynamics of a quantum system. In particular, strongly disordered many-body systems can generically avoid thermalization when there are extensively many such operators. In this work, we explicitly construct local conserved operators in one-dimensional spin chains by directly minimizing their commutator with the Hamiltonian. We demonstrate the existence of an extensively large set of local integrals of motion in the many-body localized phase of the disordered XXZ spin chain. These operators are shown to have exponentially decaying tails, in contrast to the ergodic phase where the decay is (at best) polynomial in the size of the subsystem. We study the algebraic properties of localized operators and confirm that in the many-body localized phase, they are well described by ``dressed'' spin operators.
We investigate the typicality of the growth behavior of the out-of-time-ordered commutator (OTOC) in the many-body localized (MBL) quantum spin chains across random disorder realizations. In the MBL phase of the Heisenberg XXZ chain, we find that the estimate of the OTOC fluctuates significantly with the disorder realizations at the intermediate times of the main growth. Despite the consequent failure of the disorder average in the MBL phase, we argue that the characteristic behavior of the OTOC growth can still be identified by going through individual disorder realizations. We find that a power-law-type growth behavior appears typically after a disorder-dependent relaxation period, which is very close to the t2{t}^{2}t2 form derived in the effective Hamiltonian of a fully MBL system. The characteristic growth behavior observed at an individual disorder realization is robust in our tests with various state preparations and also verified in another MBL system of the random-transverse-field quantum Ising chain in a uniform longitudinal field.
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.
Despite a very good understanding of single-particle Anderson localization in one-dimensional (1D) disordered systems, many-body effects are still full of surprises; a famous example is the interaction-driven many-body localization (MBL) problem, about which much has been written, and perhaps the best is yet to come. Interestingly enough the noninteracting limit provides a natural playground to study nontrivial multiparticle physics, offering the possibility to test some general mechanisms with very large-scale exact diagonalization simulations. In this paper, we first revisit the 1D many-body Anderson insulator through the lens of extreme value theory, focusing on the extreme polarizations of the equivalent spin chain model in a random magnetic field. A many-body-induced chain breaking mechanism is explored numerically, and compared to an analytically solvable toy model. A unified description, from weak to large disorder strengths WWW emerges, where the disorder-dependent average localization length $\ensuremath{\xi}(W)$ governs the extreme events leading to chain breaks. In particular, tails of the local magnetization distributions are controlled by $\ensuremath{\xi}(W)$. Remarkably, we also obtain a quantitative understanding of the full distribution of the extreme polarizations, which is given by a Fr\'echet-type law. In a second part, we explore finite interaction physics and the MBL question. For the available system sizes, we numerically quantify the difference in the extreme value distributions between the interacting problem and the noninteracting Anderson case. Strikingly, we observe a sharp ``extreme-statistics transition'' as WWW changes, which may coincide with the MBL transition.
The random-field spin-1/2 XXZ chains, and the corresponding Anderson insulators of spinless fermions with density-density interaction, have been intensively studied in the context of many-body localization. However, we recently argued [B. Krajewski et al., Phys. Rev. Lett. 129, 260601 (2022)] that the two-body density-density interaction in these models is not generic since only a small fraction of this interaction represents a true local perturbation to the Anderson insulator. Here we study ergodicity of strongly disordered Anderson insulator chains, choosing other forms of the two-body interaction for which the strength of the true perturbation is of the same order of magnitude as the bare two-body interaction. Focusing on the strong-interaction regime, numerical results for the level statistics and the eigenstate thermalization hypothesis are consistent with emergence of ergodicity at arbitrary strong disorder.