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

Many body localization in Disordered One-Dimensional Fermi-Hubbard Model

Harmanjeet Kaur, Vinod Ashokan

Published Jun 3, 2026• 0 citations
arXivHTMLPDF

Abstract

We investigate the non-equilibrium dynamics of the disordered one-dimensional Fermi-Hubbard model with a focus on many-body localization. The system is initialized in a charge-density-wave-state, and its time evolution is analyzed through sublattice imbalance (spin and charge), and bipartite entanglement entropy. A clear crossover from ergodic to non-ergodic behavior is observed with increasing disorder strength. In the weak disorder regime, rapid decay of imbalance and the fast growth of entanglement indicate efficient thermalization. In contrast, a strong disorder leads to persistent imbalance and slow dynamics, signaling the breakdown of ergodicity. The charge and spin sectors exhibit distinct relaxation behavior, providing evidence for partial decoupling between these degrees of freedom. Furthermore, in the interacting regime, the entanglement entropy shows slow logarithmic growth, reflecting the dephasing-driven dynamics characteristic of the many-body localized phase. These results highlight the interplay between disorder and interactions in determining the dynamical properties of the system and establish robust signatures of many-body localization in the Fermi-Hubbard model.

DOI
10.48550/arxiv.2606.04948

Related papers

Connecting entanglement growth with local integrals of motion in the disordered Fermi-Hubbard model

Ahad Nokhostin Helm, Brandon Leipner-Johns, Rachel Wortis

Generically a quantum system initialized in an unentangled state will, under unitary dynamics, rapidly become entangled, a process closely related to information transport and to thermalization. Disorder can suppress the growth of entanglement and result in memory of initial conditions. In non-interacting systems this arises from localization of single-particle states, the occupancy of which is fixed by the initial condition. In interacting systems similar localized conserved quantities persist, but with the added feature that they are coupled, resulting in entanglement growth which is distinct from both non-interacting localized systems and from generic ergodic systems. The Fermi-Hubbard model has two degrees of freedom per site -- charge and spin -- and disorder may be present in both of these. We study the growth of entanglement in two scenarios -- disorder in charge equal and unequal to that in spin, and determine the distinct contributions of charge and spin degrees of freedom by expanding the Hamiltonian in terms of a set of optimally localized conserved quantities with separate charge and spin character. We find that coupling between charge and spin is significantly weaker than charge-charge and spin-spin coupling. While this decoupling is present in all our results, it is only apparent when the strength of the disorder in the two sectors is different such that there is a separation between the characteristic timescales of the contributions to entanglement made by charge and by spin.

quant-pharXiv:2606.15481•Jun 16, 2026•0 citations

Connecting Entanglement Entropy Growth and Local Integrals of Motion in the Disordered Fermi-Hubbard Model

Ahad Nokhostin Helm

We study a one-dimensional Fermi-Hubbard model with disorder in charge and spin degrees of freedom. We calculate the time dependence of the entanglement entropy. While previous research on disordered interacting systems has typically focused on systems with either charge or spin, our model enables us to explore the interplay between charge and spin in shaping the behavior of entanglement. We use a method that identifies optimally local charge- and spin-specific integrals of motion. We ask how the locality level of these integrals of motion influences the capacity of low-order terms in the l-bit Hamiltonian to capture the entanglement entropy. Our results show that increasing the locality level improves the accuracy of low-order terms in capturing entanglement entropy dynamics. With equally strong charge and spin disorder, the behavior of the entanglement entropy closely resembles that observed in single-degree-of-freedom systems, and the l-bit Hamiltonian truncated at second order accurately captures this behavior. Author Keywords: Entanglement Entropy, Fermi-Hubbard Model, Hungarian Algorithm, l-bit Hamiltonian, Local Integrals of Motion, Many-Body Localization

•Jan 1, 2025•0 citations

Many-body localization in a disorder Bose-Hubbard chain

Jie Chen, Xiaoqun Wang

Many-body localization of a disorder interacting boson system in one dimension is studied numerically at the filling factor being one-half, in terms of level statistics, local compressibility, correlation function and entanglement entropies. The von Neumann entanglement entropy is decoupled into a particle number entropy and a configuration entropy. An anomalous volume-law behavior is found for the configuration entanglement entropy to confirm a recent experimental observation [A. Lukin, M. Rispoli, R. Schittko, et al., Science 364, 256 (2019)] for sufficient strong disorder, while the particle number entropy fulfills an area-law corresponding to the total entropy for disordered spin chain. The localization length are extracted from a two-body correlation function for many-body localization states and corresponding time-evolutions states as well. A phase diagrams is established with consisting of an ergodic thermalized region and a many-body-localization region in a parameter space of the disorder strength and the energy density. Two regions are separated by a many-body mobility edge deducted from the standard deviation of the particle-number entanglement entropy, which appears consistent with that based on the localization length. Slow dynamics characterized by a logarithmic time-dependence is explicitly shown for both the particle number entropy and the configuration entropy in an intermediate regime of their time-evolutions, which does not show up in the Anderson localization case, i.e. non-interacting disorder systems.

•Apr 17, 2021•1 citation

Spin-charge separation and many-body localization

Jakub Zakrzewski, Dominique Delande

We study many-body localization for a disordered chain of spin-1/2 fermions. When both down and up components are exposed to the same strong disorder, a power-law growth of the entanglement entropy was observed [P. Prelov\ifmmode \check{s}\else \v{s}\fi{}ek et al., Phys. Rev. B 94, 241104 (2016)] that suggests that many-body localization is not complete; the density (charge) degree of freedom is localized, while the spin degree of freedom is apparently delocalized. We show that this powerlike behavior is only a transient effect and that, for longer times, the growth of the entanglement entropy is logarithmic in time. At the same time, the decay of the spin correlation slows down, showing that the spin transport is strongly reduced at long times. The dynamics of the system is quite similar to the one of many-body localized systems. Because of the limited duration of the numerical experiment, it is not possible to determine whether the system is truly many-body localized or diffusing very slowly at long times.

arXiv:1802.02525•Physical review. B./Physical review. B•Jul 16, 2018•45 citations

Extended slow dynamical regime close to the many-body localization transition

David J. Luitz, Nicolas Laflorencie, Fabien Alet

Many-body localization is characterized by a slow logarithmic growth of the entanglement entropy after a global quantum quench while the local memory of an initial density imbalance remains at infinite time. We investigate how much the proximity of a many-body localized phase can influence the dynamics in the delocalized ergodic regime where thermalization is expected. Using an exact Krylov space technique, the out-of-equilibrium dynamics of the random-field Heisenberg chain is studied up to $L=28$ sites, starting from an initially unentangled high-energy product state. Within most of the delocalized phase, we find a sub-ballistic entanglement growth $S(t)\ensuremath{\propto}{t}^{1/z}$ with a disorder-dependent exponent $z\ensuremath{\ge}1$, in contrast with the pure ballistic growth $z=1$ of clean systems. At the same time, anomalous relaxation is also observed for the spin imbalance $\mathcal{I}(t)\ensuremath{\propto}{t}^{\ensuremath{-}\ensuremath{\zeta}}$ with a continuously varying disorder-dependent exponent $\ensuremath{\zeta}$, vanishing at the transition. This provides a clear experimental signature for detecting this nonconventional regime.

arXiv:1511.05141•Physical review. B./Physical review. B•Feb 17, 2016•319 citations

Observation of many-body localization of interacting fermions in a quasirandom optical lattice

Michael Schreiber, Sean S. Hodgman, Pranjal Bordia +6 more

Many-body localization (MBL), the disorder-induced localization of interacting particles, signals a breakdown of conventional thermodynamics because MBL systems do not thermalize and show nonergodic time evolution. We experimentally observed this nonergodic evolution for interacting fermions in a one-dimensional quasirandom optical lattice and identified the MBL transition through the relaxation dynamics of an initially prepared charge density wave. For sufficiently weak disorder, the time evolution appears ergodic and thermalizing, erasing all initial ordering, whereas above a critical disorder strength, a substantial portion of the initial ordering persists. The critical disorder value shows a distinctive dependence on the interaction strength, which is in agreement with numerical simulations. Our experiment paves the way to further detailed studies of MBL, such as in noncorrelated disorder or higher dimensions.

arXiv:1501.05661•Science•Jul 30, 2015•1,722 citations
We investigate the non-equilibrium dynamics of the disordered one-dimensional Fermi-Hubbard model with a focus on many-body localization. The system is initialized in a charge-density-wave-state, and its time evolution is analyzed through sublattice imbalance (spin and charge), and bipartite entanglement entropy. A clear crossover from ergodic to non-ergodic behavior is observed with increasing disorder strength. In the weak disorder regime, rapid decay of imbalance and the fast growth of entanglement indicate efficient thermalization. In contrast, a strong disorder leads to persistent imbalance and slow dynamics, signaling the breakdown of ergodicity. The charge and spin sectors exhibit distinct relaxation behavior, providing evidence for partial decoupling between these degrees of freedom. Furthermore, in the interacting regime, the entanglement entropy shows slow logarithmic growth, reflecting the dephasing-driven dynamics characteristic of the many-body localized phase. These results highlight the interplay between disorder and interactions in determining the dynamical properties of the system and establish robust signatures of many-body localization in the Fermi-Hubbard model.
Generically a quantum system initialized in an unentangled state will, under unitary dynamics, rapidly become entangled, a process closely related to information transport and to thermalization. Disorder can suppress the growth of entanglement and result in memory of initial conditions. In non-interacting systems this arises from localization of single-particle states, the occupancy of which is fixed by the initial condition. In interacting systems similar localized conserved quantities persist, but with the added feature that they are coupled, resulting in entanglement growth which is distinct from both non-interacting localized systems and from generic ergodic systems. The Fermi-Hubbard model has two degrees of freedom per site -- charge and spin -- and disorder may be present in both of these. We study the growth of entanglement in two scenarios -- disorder in charge equal and unequal to that in spin, and determine the distinct contributions of charge and spin degrees of freedom by expanding the Hamiltonian in terms of a set of optimally localized conserved quantities with separate charge and spin character. We find that coupling between charge and spin is significantly weaker than charge-charge and spin-spin coupling. While this decoupling is present in all our results, it is only apparent when the strength of the disorder in the two sectors is different such that there is a separation between the characteristic timescales of the contributions to entanglement made by charge and by spin.
We study a one-dimensional Fermi-Hubbard model with disorder in charge and spin degrees of freedom. We calculate the time dependence of the entanglement entropy. While previous research on disordered interacting systems has typically focused on systems with either charge or spin, our model enables us to explore the interplay between charge and spin in shaping the behavior of entanglement. We use a method that identifies optimally local charge- and spin-specific integrals of motion. We ask how the locality level of these integrals of motion influences the capacity of low-order terms in the l-bit Hamiltonian to capture the entanglement entropy. Our results show that increasing the locality level improves the accuracy of low-order terms in capturing entanglement entropy dynamics. With equally strong charge and spin disorder, the behavior of the entanglement entropy closely resembles that observed in single-degree-of-freedom systems, and the l-bit Hamiltonian truncated at second order accurately captures this behavior. Author Keywords: Entanglement Entropy, Fermi-Hubbard Model, Hungarian Algorithm, l-bit Hamiltonian, Local Integrals of Motion, Many-Body Localization
Many-body localization of a disorder interacting boson system in one dimension is studied numerically at the filling factor being one-half, in terms of level statistics, local compressibility, correlation function and entanglement entropies. The von Neumann entanglement entropy is decoupled into a particle number entropy and a configuration entropy. An anomalous volume-law behavior is found for the configuration entanglement entropy to confirm a recent experimental observation [A. Lukin, M. Rispoli, R. Schittko, et al., Science 364, 256 (2019)] for sufficient strong disorder, while the particle number entropy fulfills an area-law corresponding to the total entropy for disordered spin chain. The localization length are extracted from a two-body correlation function for many-body localization states and corresponding time-evolutions states as well. A phase diagrams is established with consisting of an ergodic thermalized region and a many-body-localization region in a parameter space of the disorder strength and the energy density. Two regions are separated by a many-body mobility edge deducted from the standard deviation of the particle-number entanglement entropy, which appears consistent with that based on the localization length. Slow dynamics characterized by a logarithmic time-dependence is explicitly shown for both the particle number entropy and the configuration entropy in an intermediate regime of their time-evolutions, which does not show up in the Anderson localization case, i.e. non-interacting disorder systems.
We study many-body localization for a disordered chain of spin-1/2 fermions. When both down and up components are exposed to the same strong disorder, a power-law growth of the entanglement entropy was observed [P. Prelov\ifmmode \check{s}\else \v{s}\fi{}ek et al., Phys. Rev. B 94, 241104 (2016)] that suggests that many-body localization is not complete; the density (charge) degree of freedom is localized, while the spin degree of freedom is apparently delocalized. We show that this powerlike behavior is only a transient effect and that, for longer times, the growth of the entanglement entropy is logarithmic in time. At the same time, the decay of the spin correlation slows down, showing that the spin transport is strongly reduced at long times. The dynamics of the system is quite similar to the one of many-body localized systems. Because of the limited duration of the numerical experiment, it is not possible to determine whether the system is truly many-body localized or diffusing very slowly at long times.
Many-body localization is characterized by a slow logarithmic growth of the entanglement entropy after a global quantum quench while the local memory of an initial density imbalance remains at infinite time. We investigate how much the proximity of a many-body localized phase can influence the dynamics in the delocalized ergodic regime where thermalization is expected. Using an exact Krylov space technique, the out-of-equilibrium dynamics of the random-field Heisenberg chain is studied up to L=28L=28L=28 sites, starting from an initially unentangled high-energy product state. Within most of the delocalized phase, we find a sub-ballistic entanglement growth $S(t)\ensuremath{\propto}{t}^{1/z}$ with a disorder-dependent exponent $z\ensuremath{\ge}1$, in contrast with the pure ballistic growth z=1z=1z=1 of clean systems. At the same time, anomalous relaxation is also observed for the spin imbalance $\mathcal{I}(t)\ensuremath{\propto}{t}^{\ensuremath{-}\ensuremath{\zeta}}$ with a continuously varying disorder-dependent exponent $\ensuremath{\zeta}$, vanishing at the transition. This provides a clear experimental signature for detecting this nonconventional regime.
Many-body localization (MBL), the disorder-induced localization of interacting particles, signals a breakdown of conventional thermodynamics because MBL systems do not thermalize and show nonergodic time evolution. We experimentally observed this nonergodic evolution for interacting fermions in a one-dimensional quasirandom optical lattice and identified the MBL transition through the relaxation dynamics of an initially prepared charge density wave. For sufficiently weak disorder, the time evolution appears ergodic and thermalizing, erasing all initial ordering, whereas above a critical disorder strength, a substantial portion of the initial ordering persists. The critical disorder value shows a distinctive dependence on the interaction strength, which is in agreement with numerical simulations. Our experiment paves the way to further detailed studies of MBL, such as in noncorrelated disorder or higher dimensions.