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:2601.12446 · New Journal of Physics

Operator delocalization in disordered spin chains via exact MPO marginals

Jonnathan Pineda, Mario Collura, Gianluca Passarelli, Procolo Lucignaon, Davide Rossini, Angelo Russomanno

quant-phcond-mat.dis-nncond-mat.str-elSubmitted Jan 18, 2026• Announced Jan 20, 2026• 0 citations
v1 · Jan 18, 2026v2 · Mar 2, 2026v3 · Jun 22, 2026
arXivHTMLPDF

Abstract

We investigate operator delocalization in disordered one-dimensional spin chains by introducing -- besides the already known operator mass -- a complementary measure of operator complexity: the operator length. Like the operator nonstabilizerness, both these quantities are defined from the expansion of time-evolved operators in the Pauli basis. They characterize, respectively, the number of sites on which an operator acts nontrivially and the spatial extent of its support. We show that both the operator mass and length can be computed efficiently and exactly within a matrix-product-operator (MPS) framework, providing direct access to their full probability distributions, without resorting to stochastic sampling. Applying this approach to the disordered XXZ spin-1/2 chain, we find sharply distinct behaviors in non-interacting and interacting regimes. In the Anderson-localized case, operator mass, length, and operator entanglement entropy rapidly saturate, signaling the absence of scrambling. By contrast, in the many-body localized (MBL) regime, for arbitrarily weak interactions, all quantities exhibit a robust logarithmic growth in time, consistent with the known logarithmic light cone of quantum-correlation propagation in MBL. We demonstrate that this behavior is quantitatively captured by an effective $\ell$-bit model and persists across system sizes accessible via tensor-network simulations.

Journal reference
New J. Phys. 28 064512 (2026)
DOI
10.1088/1367-2630/ae774b
Comments
5 figures, 31 pages

Related papers

Operator and entanglement growth in nonthermalizing systems: Many-body localization and the random singlet phase

Ian MacCormack, Mao Tian Tan, Jonah Kudler-Flam +1 more

We characterize the growth and spreading of operators and entanglement in two paradigmatic nonthermalizing phases---the many-body localized (MBL) phase and the random singlet phase (RSP)---using the entanglement contour and multipartite operator entanglement measures. The entanglement contour characterizes the spacetime spreading of entanglement and reveals logarithmically growing entanglement light cones in the MBL and RSP phases, sharply contrasting the linear light cones of clean, thermalizing systems. The operator entanglement characterizes scrambling, i.e., the delocalization of information. We find slow scrambling behavior in the MBL phase; the late-time value of the tripartite mutual information scales linearly with system size, but is submaximal. The tripartite logarithmic negativity is also negative and nonzero, but smaller in magnitude, revealing an intriguing distinction between classical and quantum information scrambling in the MBL phase. This is in contrast with the RSP, which, as a noninteracting model, is nonscrambling.

arXiv:2001.08222•Physical review. B./Physical review. B•Dec 2, 2021•50 citations

Many-body localization and delocalization in large quantum chains

Elmer V. H. Doggen, Frank Schindler, Konstantin S. Tikhonov +4 more

We theoretically study the quench dynamics for an isolated Heisenberg spin chain with a random on-site magnetic field, which is one of the paradigmatic models of a many-body localization transition. We use the time-dependent variational principle as applied to matrix product states, which allows us to controllably study chains of a length up to $L=100$ spins, i.e., much larger than $L\ensuremath{\simeq}20$ that can be treated via exact diagonalization. For the analysis of the data, three complementary approaches are used: (i) determination of the exponent $\ensuremath{\beta}$ which characterizes the power-law decay of the antiferromagnetic imbalance with time; (ii) similar determination of the exponent ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ which characterizes the decay of a Schmidt gap in the entanglement spectrum; and (iii) machine learning with the use, as an input, of the time dependence of the spin densities in the whole chain. We find that the consideration of the larger system sizes substantially increases the estimate for the critical disorder ${W}_{c}$ that separates the ergodic and many-body localized regimes, compared to the values of ${W}_{c}$ in the literature. On the ergodic side of the transition, there is a broad interval of the strength of disorder with slow subdiffusive transport. In this regime, the exponents $\ensuremath{\beta}$ and ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ increase, with increasing $L$, for relatively small $L$ but saturate for $L\ensuremath{\simeq}50$, indicating that these slow power laws survive in the thermodynamic limit. From a technical perspective, we develop an adaptation of the ``learning by confusion'' machine-learning approach that can determine ${W}_{c}$.

arXiv:1807.05051•Physical review. B./Physical review. B•Nov 8, 2018•186 citations

Scrambling Dynamics across a Thermalization-Localization Quantum Phase Transition

Subhayan Sahu, Shenglong Xu, Brian G. Swingle

We study quantum information scrambling, specifically the growth of Heisenberg operators, in large disordered spin chains using matrix product operator dynamics to scan across the thermalization-localization quantum phase transition. We observe ballistic operator growth for weak disorder, and a sharp transition to a phase with subballistic operator spreading. The critical disorder strength for the ballistic to subballistic transition is well below the many body localization phase transition, as determined from finite size scaling of energy eigenstate entanglement entropy in small chains. In contrast, we find that the transition from subballistic to logarithmic behavior at the actual eigenstate localization transition is not resolved in our finite numerics. These data are discussed in the context of a universal form for the growing operator shape and substantiated with a simple phenomenological model of rare regions.

arXiv:1807.06086•Physical Review Letters•Oct 18, 2019•54 citations

Effect of long-range hopping and interactions on entanglement dynamics and many-body localization

Rajeev Singh, Roderich Moessner, Dibyendu Roy

We numerically investigate the dynamics of entanglement in a chain of spinless fermions with nonrandom but long-range hopping and interactions, and with random on-site energies. For moderate disorder in the absence of interactions, the chain hosts delocalized states at the top of the band which undergo a delocalization-localization transition with increasing disorder. We find an interesting regime in this noninteracting disordered chain where the long-time entanglement entropy scales as $S(t)\ensuremath{\sim}lnt$ and the saturated entanglement entropy scales with system size $L$ as $S(L,t\ensuremath{\rightarrow}\ensuremath{\infty})\ensuremath{\sim}lnL$. We further study the interplay of long-range hopping and interactions on the growth of entanglement and the many-body localization (MBL) transition in this system. We develop an analogy to higher-dimensional short-range systems to compare and contrast such behavior with the physics of MBL in a higher dimension.

arXiv:1606.04542•Physical review. B./Physical review. B•Mar 21, 2017•47 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

Large deviations in the many-body localization transition: The case of the random-field XXZ chain

Greivin Alfaro Miranda, Fabien Alet, Giulio Biroli +3 more

The effect of rare system-wide resonances in the many-body localization (MBL) transition has recently attracted significant attention. They are expected to play a prominent role in the stability of the MBL phase, prompting the development of new theoretical frameworks to properly account for their statistical weight. We employ a method based on an analogy with mean-field disordered glassy systems to characterize the statistics of transmission amplitudes between distant many-body configurations in Hilbert space, and apply it to the random-field XXZ spin chain. By introducing a Lagrange multiplier, which formally plays the role of an effective temperature controlling the influence of extreme outliers in the heavy-tailed distribution of propagators, we identify three distinct regimes: (i) an ergodic phase with uniform spreading in Hilbert space, (ii) an intermediate regime where delocalization is driven by rare, disorder-dependent long-range resonances, and (iii) a robust MBL phase where such resonances cannot destabilize localization. We derive a finite-size phase diagram in the disorder--interaction plane both in the spin and in the Anderson basis that quantitatively agrees with recent numerical results based on real-space spin-spin correlation functions. We further demonstrate that even infinitesimal interactions can destroy the Anderson insulator at finite disorder, with the critical disorder remaining finite down to small interaction strengths. By visualizing resonant transmission pathways on the Hilbert space graph, we provide a complementary perspective to real-space and spectral probes, revealing how the destabilization of the MBL phase at finite sizes stems from the emergence of resonant paths that become progressively rarer and shorter-ranged deep in the localized phase.

arXiv:2510.18545•Physical review. B./Physical review. B•Jul 10, 2026•1 citation
We investigate operator delocalization in disordered one-dimensional spin chains by introducing -- besides the already known operator mass -- a complementary measure of operator complexity: the operator length. Like the operator nonstabilizerness, both these quantities are defined from the expansion of time-evolved operators in the Pauli basis. They characterize, respectively, the number of sites on which an operator acts nontrivially and the spatial extent of its support. We show that both the operator mass and length can be computed efficiently and exactly within a matrix-product-operator (MPS) framework, providing direct access to their full probability distributions, without resorting to stochastic sampling. Applying this approach to the disordered XXZ spin-1/2 chain, we find sharply distinct behaviors in non-interacting and interacting regimes. In the Anderson-localized case, operator mass, length, and operator entanglement entropy rapidly saturate, signaling the absence of scrambling. By contrast, in the many-body localized (MBL) regime, for arbitrarily weak interactions, all quantities exhibit a robust logarithmic growth in time, consistent with the known logarithmic light cone of quantum-correlation propagation in MBL. We demonstrate that this behavior is quantitatively captured by an effective ℓ\ellℓ-bit model and persists across system sizes accessible via tensor-network simulations.
We characterize the growth and spreading of operators and entanglement in two paradigmatic nonthermalizing phases---the many-body localized (MBL) phase and the random singlet phase (RSP)---using the entanglement contour and multipartite operator entanglement measures. The entanglement contour characterizes the spacetime spreading of entanglement and reveals logarithmically growing entanglement light cones in the MBL and RSP phases, sharply contrasting the linear light cones of clean, thermalizing systems. The operator entanglement characterizes scrambling, i.e., the delocalization of information. We find slow scrambling behavior in the MBL phase; the late-time value of the tripartite mutual information scales linearly with system size, but is submaximal. The tripartite logarithmic negativity is also negative and nonzero, but smaller in magnitude, revealing an intriguing distinction between classical and quantum information scrambling in the MBL phase. This is in contrast with the RSP, which, as a noninteracting model, is nonscrambling.
We theoretically study the quench dynamics for an isolated Heisenberg spin chain with a random on-site magnetic field, which is one of the paradigmatic models of a many-body localization transition. We use the time-dependent variational principle as applied to matrix product states, which allows us to controllably study chains of a length up to L=100L=100L=100 spins, i.e., much larger than $L\ensuremath{\simeq}20$ that can be treated via exact diagonalization. For the analysis of the data, three complementary approaches are used: (i) determination of the exponent $\ensuremath{\beta}$ which characterizes the power-law decay of the antiferromagnetic imbalance with time; (ii) similar determination of the exponent ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ which characterizes the decay of a Schmidt gap in the entanglement spectrum; and (iii) machine learning with the use, as an input, of the time dependence of the spin densities in the whole chain. We find that the consideration of the larger system sizes substantially increases the estimate for the critical disorder Wc{W}_{c}Wc​ that separates the ergodic and many-body localized regimes, compared to the values of Wc{W}_{c}Wc​ in the literature. On the ergodic side of the transition, there is a broad interval of the strength of disorder with slow subdiffusive transport. In this regime, the exponents $\ensuremath{\beta}$ and ${\ensuremath{\beta}}_{\mathrm{\ensuremath{\Lambda}}}$ increase, with increasing LLL, for relatively small LLL but saturate for $L\ensuremath{\simeq}50$, indicating that these slow power laws survive in the thermodynamic limit. From a technical perspective, we develop an adaptation of the ``learning by confusion'' machine-learning approach that can determine Wc{W}_{c}Wc​.
We study quantum information scrambling, specifically the growth of Heisenberg operators, in large disordered spin chains using matrix product operator dynamics to scan across the thermalization-localization quantum phase transition. We observe ballistic operator growth for weak disorder, and a sharp transition to a phase with subballistic operator spreading. The critical disorder strength for the ballistic to subballistic transition is well below the many body localization phase transition, as determined from finite size scaling of energy eigenstate entanglement entropy in small chains. In contrast, we find that the transition from subballistic to logarithmic behavior at the actual eigenstate localization transition is not resolved in our finite numerics. These data are discussed in the context of a universal form for the growing operator shape and substantiated with a simple phenomenological model of rare regions.
We numerically investigate the dynamics of entanglement in a chain of spinless fermions with nonrandom but long-range hopping and interactions, and with random on-site energies. For moderate disorder in the absence of interactions, the chain hosts delocalized states at the top of the band which undergo a delocalization-localization transition with increasing disorder. We find an interesting regime in this noninteracting disordered chain where the long-time entanglement entropy scales as $S(t)\ensuremath{\sim}lnt$ and the saturated entanglement entropy scales with system size LLL as $S(L,t\ensuremath{\rightarrow}\ensuremath{\infty})\ensuremath{\sim}lnL$. We further study the interplay of long-range hopping and interactions on the growth of entanglement and the many-body localization (MBL) transition in this system. We develop an analogy to higher-dimensional short-range systems to compare and contrast such behavior with the physics of MBL in a higher dimension.
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 effect of rare system-wide resonances in the many-body localization (MBL) transition has recently attracted significant attention. They are expected to play a prominent role in the stability of the MBL phase, prompting the development of new theoretical frameworks to properly account for their statistical weight. We employ a method based on an analogy with mean-field disordered glassy systems to characterize the statistics of transmission amplitudes between distant many-body configurations in Hilbert space, and apply it to the random-field XXZ spin chain. By introducing a Lagrange multiplier, which formally plays the role of an effective temperature controlling the influence of extreme outliers in the heavy-tailed distribution of propagators, we identify three distinct regimes: (i) an ergodic phase with uniform spreading in Hilbert space, (ii) an intermediate regime where delocalization is driven by rare, disorder-dependent long-range resonances, and (iii) a robust MBL phase where such resonances cannot destabilize localization. We derive a finite-size phase diagram in the disorder--interaction plane both in the spin and in the Anderson basis that quantitatively agrees with recent numerical results based on real-space spin-spin correlation functions. We further demonstrate that even infinitesimal interactions can destroy the Anderson insulator at finite disorder, with the critical disorder remaining finite down to small interaction strengths. By visualizing resonant transmission pathways on the Hilbert space graph, we provide a complementary perspective to real-space and spectral probes, revealing how the destabilization of the MBL phase at finite sizes stems from the emergence of resonant paths that become progressively rarer and shorter-ranged deep in the localized phase.