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

Spin and energy diffusion versus subdiffusion in disordered spin chains

J. Herbrych, Peter Prelovšek

Published Jun 20, 2025• 2 citations
arXivHTMLPDF

Abstract

While the high-temperature spin diffusion in spin chains with random local fields has been the subject of numerous studies concerning the phenomenon of many-body localization (MBL), the energy diffusion in the same models has been much less explored. We demonstrate that energy diffusion is faster at weak random fields but becomes essentially equal at strong fields; hence, both diffusions determine the slowest relaxation time scale (Thouless time) in the system. Numerically reachable finite-size systems reveal the anomalously large distribution of diffusion constants with respect to actual field configurations. Despite the exponential-like dependence of diffusion on field strength, the results for sensitivity to twisted boundary conditions are incompatible with the Thouless criterion for localization and the presumed transition to MBL, at least for numerically reachable sizes. In contrast, we find indications of the scenario of subdiffusive transport, particularly in the dynamical diffusivity response.

DOI
10.1103/kv2f-m8vk

Related papers

Slow diffusion and Thouless localization criterion in modulated spin chains

Peter Prelovšek, J. Herbrych, Marcin Mierzejewski

In recent years the ergodicity of disordered spin chains has been investigated via extensive numerical studies of the level statistics or the transport properties. However, a clear relationship between these results has yet to be established. We present the relation between the diffusion constant and the energy-level structure, which leads to the Thouless localization criterion. Together with the exponential-like dependence of the diffusion constant on the strength of quasiperiodic or random fields, the Thouless criterion explains the nearly linear drift with the system size of the crossover/transition to the nonergodic regime. Moreover, we show that the Heisenberg spin chain in the presence of the quasiperiodic fields can be well approached via a sequence of simple periodic systems, where diffusion remains finite even at large fields.

arXiv:2302.03325•Physical review. B./Physical review. B•Jul 5, 2023•12 citations

Asymmetry in energy versus spin transport in certain interacting disordered systems

Juan José Mendoza-Arenas, Marko Žnidarič, Vipin Kerala Varma +3 more

We study energy transport in disordered $XXZ$ spin-1/2 chains driven to nonequilibrium configurations by thermal reservoirs of different temperatures at the boundaries, using large-scale matrix product simulations. In particular we discuss the transition between diffusive and subdiffusive transport in sectors of zero and finite magnetization at high temperature. At large anisotropies we find that diffusive energy transport prevails over a large range of disorder strengths, which is in contrast to spin transport that is subdiffusive in the same regime for weak disorder. However, at finite magnetization both energy and spin currents decay as a function of the system size with the same exponent. We conclude that diffusion of energy is much more pervasive than that of magnetization in these disordered spin-1/2 systems, and occurs across a significant range of the interaction-disorder parameter phase space. We support the existence of this asymmetry, reminiscent of that in the clean limit, by an analytical estimation of diffusion constants for weak disorder.

arXiv:1803.11555•Physical review. B./Physical review. B•Mar 22, 2019•40 citations

Energy transport in a disordered spin chain with broken U(1) symmetry: Diffusion, subdiffusion, and many-body localization

Maximilian Schulz, Scott R. Taylor, C. Hooley +1 more

We explore the physics of the disordered XYZ spin chain using two complementary numerical techniques: exact diagonalization (ED) on chains of up to 17 spins, and time-evolving block decimation (TEBD) on chains of up to 400 spins. Our principal findings are as follows. First, we verify that the clean XYZ spin chain shows ballistic energy transport for all parameter values that we investigated. Second, for weak disorder there is a stable diffusive region that persists up to a critical disorder strength that depends on the XY anisotropy. Third, for disorder strengths above this critical value, energy transport becomes increasingly subdiffusive. Fourth, the many-body localization transition moves to significantly higher disorder strengths as the XY anisotropy is increased. We discuss these results, and their relation to our current physical picture of subdiffusion in the approach to many-body localization.

arXiv:1805.01036•Physical review. B./Physical review. B•Nov 5, 2018•39 citations

Subdiffusion in a disordered spin chain with random dephasing: Finite-size corrections, Griffiths effects, and implications for many-body localization

Scott R. Taylor, Antonello Scardicchio

We study transport in the boundary-driven XX spin chain with onsite disorder and randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. We then present an exactly solvable semiclassical model of conductors and insulators, which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The subdiffusion in these models is a consequence of rare insulating regions, and they therefore exhibit, in a more controlled setup, the physics of which have been conjectured to cause the subdiffusive transport observed in interacting many-body localizable systems. For the quantum model we show that finite-size effects come from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. The latter grows logarithmically with system size, and we conjecture that this may be the reason why the heavy-tailed distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.

•Jul 27, 2020•2 citations

Energy transport in disordered classical spin chains

Vadim Oganesyan, Arijeet Pal, David A. Huse

We present a numerical study of the diffusion of energy at high temperature in strongly disordered chains of interacting classical spins evolving deterministically. We find that quenched randomness strongly suppresses transport with the diffusion constant becoming reduced by several orders of magnitude upon the introduction of moderate disorder. We have also looked for but not found signs of a classical many-body localization transition at any nonzero strength of the spin-spin interactions.

arXiv:0905.4112•Physical Review B•Sep 3, 2009•112 citations

Subdiffusion in a one-dimensional Anderson insulator with random dephasing: Finite-size scaling, Griffiths effects, and possible implications for many-body localization

Scott R. Taylor, Antonello Scardicchio

We study transport in a one-dimensional boundary-driven Anderson insulator (the XX spin chain with onsite disorder) with randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. This model is intended to mimic the passage of an excitation through (many-body) insulating regions or ergodic bubbles, therefore providing a toy model for the diffusion-subdiffusion transition observed in the disordered Heisenberg model by \ifmmode \check{Z}\else \v{Z}\fi{}nidari\ifmmode \check{c}\else \v{c}\fi{} et al. [Phys. Rev. Lett. 117, 040601 (2016)]. We also present the exact solution of a semiclassical model of conductors and insulators introduced by Agarwal et al. [Phys. Rev. Lett. 114, 160401 (2015)], which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The critical properties of both models, when passing from diffusion to subdiffusion, are interpreted in terms of ``Griffiths effects.'' We show that the finite-size scaling comes from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. We conjecture that the latter, which grows logarithmically with system size, may potentially be responsible for the fact that heavy-tailed resistance distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.

arXiv:2007.13783•Physical review. B./Physical review. B•May 7, 2021•35 citations
While the high-temperature spin diffusion in spin chains with random local fields has been the subject of numerous studies concerning the phenomenon of many-body localization (MBL), the energy diffusion in the same models has been much less explored. We demonstrate that energy diffusion is faster at weak random fields but becomes essentially equal at strong fields; hence, both diffusions determine the slowest relaxation time scale (Thouless time) in the system. Numerically reachable finite-size systems reveal the anomalously large distribution of diffusion constants with respect to actual field configurations. Despite the exponential-like dependence of diffusion on field strength, the results for sensitivity to twisted boundary conditions are incompatible with the Thouless criterion for localization and the presumed transition to MBL, at least for numerically reachable sizes. In contrast, we find indications of the scenario of subdiffusive transport, particularly in the dynamical diffusivity response.
In recent years the ergodicity of disordered spin chains has been investigated via extensive numerical studies of the level statistics or the transport properties. However, a clear relationship between these results has yet to be established. We present the relation between the diffusion constant and the energy-level structure, which leads to the Thouless localization criterion. Together with the exponential-like dependence of the diffusion constant on the strength of quasiperiodic or random fields, the Thouless criterion explains the nearly linear drift with the system size of the crossover/transition to the nonergodic regime. Moreover, we show that the Heisenberg spin chain in the presence of the quasiperiodic fields can be well approached via a sequence of simple periodic systems, where diffusion remains finite even at large fields.
We study energy transport in disordered XXZXXZXXZ spin-1/2 chains driven to nonequilibrium configurations by thermal reservoirs of different temperatures at the boundaries, using large-scale matrix product simulations. In particular we discuss the transition between diffusive and subdiffusive transport in sectors of zero and finite magnetization at high temperature. At large anisotropies we find that diffusive energy transport prevails over a large range of disorder strengths, which is in contrast to spin transport that is subdiffusive in the same regime for weak disorder. However, at finite magnetization both energy and spin currents decay as a function of the system size with the same exponent. We conclude that diffusion of energy is much more pervasive than that of magnetization in these disordered spin-1/2 systems, and occurs across a significant range of the interaction-disorder parameter phase space. We support the existence of this asymmetry, reminiscent of that in the clean limit, by an analytical estimation of diffusion constants for weak disorder.
We explore the physics of the disordered XYZ spin chain using two complementary numerical techniques: exact diagonalization (ED) on chains of up to 17 spins, and time-evolving block decimation (TEBD) on chains of up to 400 spins. Our principal findings are as follows. First, we verify that the clean XYZ spin chain shows ballistic energy transport for all parameter values that we investigated. Second, for weak disorder there is a stable diffusive region that persists up to a critical disorder strength that depends on the XY anisotropy. Third, for disorder strengths above this critical value, energy transport becomes increasingly subdiffusive. Fourth, the many-body localization transition moves to significantly higher disorder strengths as the XY anisotropy is increased. We discuss these results, and their relation to our current physical picture of subdiffusion in the approach to many-body localization.
We study transport in the boundary-driven XX spin chain with onsite disorder and randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. We then present an exactly solvable semiclassical model of conductors and insulators, which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The subdiffusion in these models is a consequence of rare insulating regions, and they therefore exhibit, in a more controlled setup, the physics of which have been conjectured to cause the subdiffusive transport observed in interacting many-body localizable systems. For the quantum model we show that finite-size effects come from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. The latter grows logarithmically with system size, and we conjecture that this may be the reason why the heavy-tailed distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.
We study transport in a one-dimensional boundary-driven Anderson insulator (the XX spin chain with onsite disorder) with randomly positioned onsite dephasing, observing a transition from diffusive to subdiffusive spin transport below a critical density of sites with dephasing. This model is intended to mimic the passage of an excitation through (many-body) insulating regions or ergodic bubbles, therefore providing a toy model for the diffusion-subdiffusion transition observed in the disordered Heisenberg model by \ifmmode \check{Z}\else \v{Z}\fi{}nidari\ifmmode \check{c}\else \v{c}\fi{} et al. [Phys. Rev. Lett. 117, 040601 (2016)]. We also present the exact solution of a semiclassical model of conductors and insulators introduced by Agarwal et al. [Phys. Rev. Lett. 114, 160401 (2015)], which exhibits both diffusive and subdiffusive phases, and qualitatively reproduces the results of the quantum system. The critical properties of both models, when passing from diffusion to subdiffusion, are interpreted in terms of ``Griffiths effects.'' We show that the finite-size scaling comes from the interplay of three characteristic lengths: one associated with disorder (the localization length), one with dephasing, and the third with the percolation problem defining large, rare, insulating regions. We conjecture that the latter, which grows logarithmically with system size, may potentially be responsible for the fact that heavy-tailed resistance distributions typical of Griffiths effects have not been observed in subdiffusive interacting systems.