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:1910.14048 · WORLD SCIENTIFIC eBooks

Thermalization and Its Absence within Krylov Subspaces of a Constrained Hamiltonian

Sanjay Moudgalya, Abhinav Prem, Rahul Nandkishore, Nicolas Regnault, Bogdan Andrei Bernevig

Published Sep 1, 2021• 92 citations
arXivHTMLPDF

Abstract

We study the quantum dynamics of a simple translation invariant, center-of-mass (CoM) preserving model of interacting fermions in one dimension (1D), which arises in multiple experimentally realizable contexts. We show that this model naturally displays the phenomenology associated with fractonic systems, wherein single charges can only move by emitting dipoles. This allows us to demonstrate the rich Krylov fractured structure of this model, whose Hilbert space shatters into exponentially many dynamically disconnected subspaces. Focusing on exponentially large Krylov subspaces, we show that these can be either be integrable or non-integrable, thereby establishing the notion of Krylov-restricted thermalization. We analytically find a tower of integrable Krylov subspaces of this Hamiltonian, all of which map onto spin-1/2 XX models of various system sizes. We also discuss the physics of the non-integrable subspaces, where we show evidence for weak Eigenstate Thermalization Hypothesis (ETH) restricted to each non-integrable Krylov subspace. Further, we show that constraints in some of the thermal Krylov subspaces cause the long-time expectation values of local operators to deviate from behavior typically expected from translation-invariant systems. Finally, we show using a Schrieffer-Wolff transformation that such models naturally appear as effective Hamiltonians in the large electric field limit of the interacting Wannier-Stark problem, and comment on connections of our work with the phenomenon of Bloch many-body localization.

DOI
10.1142/9789811231711_0009

Related papers

Ergodicity Breaking Arising from Hilbert Space Fragmentation in Dipole-Conserving Hamiltonians

Pablo Sala, Tibor Rakovszky, Ruben Verresen +2 more

We show that the combination of charge and dipole conservation-characteristic of fracton systemsleads to an extensive fragmentation of the Hilbert space, which, in turn, can lead to a breakdown of thermalization.As a concrete example, we investigate the out-of-equilibrium dynamics of one-dimensional spin-1 models that conserve charge (total S z ) and its associated dipole moment.First, we consider a minimal model including only three-site terms and find that the infinite temperature autocorrelation saturates to a finite value-showcasing nonthermal behavior.The absence of thermalization is identified as a consequence of the strong fragmentation of the Hilbert space into exponentially many invariant subspaces in the local S z basis, arising from the interplay of dipole conservation and local interactions.Second, we extend the model by including four-site terms and find that this perturbation leads to a weak fragmentation: The system still has exponentially many invariant subspaces, but they are no longer sufficient to avoid thermalization for typical initial states.More generally, for any finite range of interactions, the system still exhibits nonthermal eigenstates appearing throughout the entire spectrum.We compare our results to charge and dipole moment-conserving random unitary circuit models for which we reach identical conclusions.

arXiv:1904.04266•Physical Review X•Feb 26, 2020•438 citations

Ergodicity Breaking Arising from Hilbert Space Fragmentation in Dipole-Conserving Hamiltonians

Sala, P., Tibor Rakovszky, Ruben Verresen +2 more

We show that the combination of charge and dipole conservation-characteristic of fracton systems-leads to an extensive fragmentation of the Hilbert space, which, in turn, can lead to a breakdown of thermalization. As a concrete example, we investigate the out-of-equilibrium dynamics of one-dimensional spin-1 models that conserve charge (total S-z) and its associated dipole moment. First, we consider a minimal model including only three-site terms and find that the infinite temperature autocorrelation saturates to a finite value-showcasing nonthermal behavior. The absence of thermalization is identified as a consequence of the strong fragmentation of the Hilbert space into exponentially many invariant subspaces in the local S-z basis, arising from the interplay of dipole conservation and local interactions. Second, we extend the model by including four-site terms and find that this perturbation leads to a weak fragmentation: The system still has exponentially many invariant subspaces, but they are no longer sufficient to avoid thermalization for typical initial states. More generally, for any finite range of interactions, the system still exhibits nonthermal eigenstates appearing throughout the entire spectrum. We compare our results to charge and dipole moment-conserving random unitary circuit models for which we reach identical conclusions.

•Jan 1, 2020•362 citations

Thermalization of isolated quantum systems

Sergei Khlebnikov, Martı́n Kruczenski

Understanding the evolution towards thermal equilibrium of an isolated quantum system is at the foundation of statistical mechanics and a subject of interest in such diverse areas as cold atom physics or the quantum mechanics of black holes. Since a pure state can never evolve into a thermal density matrix, the Eigenstate Thermalization Hypothesis (ETH) has been put forward by Deutsch and Srednicki as a way to explain this apparent thermalization, similarly to what the ergodic theorem does in classical mechanics. In this paper this hypothesis is tested numerically. First, it is observed that thermalization happens in a subspace of states (the Krylov subspace) with dimension much smaller than that of the total Hilbert space. We check numerically the validity of ETH in such a subspace, for a system of hard core bosons on a two-dimensional lattice. We then discuss how well the eigenstates of the Hamiltonian projected on the Krylov subspace represent the true eigenstates. This discussion is aided by bringing the projected Hamiltonian to the tridiagonal form and interpreting it as an Anderson localization problem for a finite one-dimensional chain. We also consider thermalization of a subsystem and argue that generation of a large entanglement entropy can lead to a thermal density matrix for the subsystem well before the whole system thermalizes. Finally, we comment on possible implications of ETH in quantum gravity.

arXiv:1312.4612•Dec 17, 2013•10 citations

Thermalization Dynamics in Closed Quantum Many Body Systems: a Precision Large Scale Exact Diagonalization Study

Ivo A. Maceira, Andreas Martin Läuchli

Using a Krylov-subspace time evolution algorithm, we simulate the real-time dynamics of translation invariant non-integrable finite spin rings to quite long times with high accuracy. We systematically study the finite-size deviation between the resulting equilibrium state and the thermal state, and we highlight the importance of the energy variance on the deviations. We find that the deviations are well described by the eigenstate thermalization hypothesis, and that the von Neumann entropy correction scaling is the square of the local operator scaling. We reveal also an area law contribution to the relaxed von Neumann entropy, which we connect to the mutual information between the considered subsystem and its immediate environment. We also find that local observables relax towards equilibrium exponentially with a relaxation time scale that grows linearly with system length and is somewhat independent of the local operator, but depends strongly on the energy of the initial state, with the fastest relaxation times found towards one end of the overall energy spectrum. To contrast this behaviour we also study domain wall initial states, which exhibit clear diffusive behaviour, with a Thouless time scaling quadratically with the systems size, leading to a rather precise estimate for the diffusion constant for states in the vicinity of the middle of the energy spectrum.

arXiv:2409.18863•Sep 27, 2024•0 citations

Aspects of Hilbert space fragmentation in the quantum East model: Fragmentation, subspace-restricted quantum scars, and effects of density-density interactions

Maitri Ganguli, Sreemayee Aditya, Diptiman Sen

We investigate a one-dimensional correlated-hopping model of spinless fermions with an East constraint. We first analytically unravel the complete fragmentation structure of this model by labeling each fragment by a unique root configuration and utilizing the transfer matrix method. We show that the growth of the size of each fragment of the model follows the widely studied Dyck sequence and is therefore analytically tractable with the help of Catalan triangles. While the eigenstate thermalization hypothesis (ETH) does not hold within the full Hilbert space which exhibits Poisson statistics of the energy level spacing, an examination of various quantities restricted to the largest fragments shows that a weaker version of the subspace-restricted thermalization holds. This weaker violation of the ETH within the largest fragments is supported by the presence of subspace-restricted quantum many-body scars due to quantum fragmentation. Our analysis further reveals that the filling fraction at which the model has the largest fragment shifts with increasing system sizes from half-filling due to the distinct fragmentation structure. Next, we show that the inclusion of a nearest-neighbor density-density interaction with strength $V$ induces a spectral transition within the largest fragment from a weakly ETH-violating phase containing scars to a statistical bubble localized phase as $V$ increases. In particular, the $V\ensuremath{\rightarrow}\ensuremath{\infty}$ limit produces an integrable model. We find that the addition of finite-$V$ stabilizes the ground state near half-filling while keeping intact the fragmentation structure of the East model. However, this behavior abruptly changes exactly at $V=\ensuremath{\infty}$ due to the emergence of a distinct fragmentation structure; we analytically explore this in a comprehensive manner. The infinite-$V$ model has many interesting properties, among which the appearance of the ground state and the largest fragment at two different filling fractions is specially noteworthy. Finally, we propose an experimental setup to realize the infinite-$V$ model as a particular limit of a special kind of $t\text{\ensuremath{-}}V$ model with an on-site potential.

arXiv:2409.15943•Physical review. B./Physical review. B•Jan 9, 2025•9 citations

Geometric fragmentation and anomalous thermalization in cubic dimer model

Joel Steinegger, Debasish Banerjee, Emilie Huffman +1 more

While quantum statistical mechanics triumphs in explaining many equilibrium phenomena, there is an increasing focus on going beyond conventional scenarios of thermalization. Traditionally examples of nonthermalizing systems are either integrable or disordered. Recently, examples of translationally invariant physical systems have been discovered whose excited energies avoid thermalization either due to local constraints (whether exact or emergent) or due to higher-form symmetries. In this article, we extend these investigations for the case of 3D U ( 1 ) quantum dimer models, which are lattice gauge theories with finite-dimensional local Hilbert spaces (also generically called quantum link models) with staggered charged static matter. Using a combination of analytical and numerical methods, we uncover a class of athermal states that arise in large winding sectors, when the system is subjected to external electric fields. The polarization of the dynamical fluxes in the direction of applied field traps excitations in 2D planes, while an interplay with the Gauss law constraint in the perpendicular direction causes exotic athermal behavior due to the emergence of new conserved quantities. This causes a geometric fragmentation of the system. We provide analytical arguments showing that the scaling of the number of fragments is exponential in the linear system size, leading to weak fragmentation. Further, we identify sectors which host fractonic excitations with severe mobility restrictions. The unitary evolution of fragments dominated by fractons is qualitatively different from the one dominated by nonfractonic excitations.

arXiv:2508.03802•Physical review. D/Physical review. D.•Dec 4, 2025•1 citation
We study the quantum dynamics of a simple translation invariant, center-of-mass (CoM) preserving model of interacting fermions in one dimension (1D), which arises in multiple experimentally realizable contexts. We show that this model naturally displays the phenomenology associated with fractonic systems, wherein single charges can only move by emitting dipoles. This allows us to demonstrate the rich Krylov fractured structure of this model, whose Hilbert space shatters into exponentially many dynamically disconnected subspaces. Focusing on exponentially large Krylov subspaces, we show that these can be either be integrable or non-integrable, thereby establishing the notion of Krylov-restricted thermalization. We analytically find a tower of integrable Krylov subspaces of this Hamiltonian, all of which map onto spin-1/2 XX models of various system sizes. We also discuss the physics of the non-integrable subspaces, where we show evidence for weak Eigenstate Thermalization Hypothesis (ETH) restricted to each non-integrable Krylov subspace. Further, we show that constraints in some of the thermal Krylov subspaces cause the long-time expectation values of local operators to deviate from behavior typically expected from translation-invariant systems. Finally, we show using a Schrieffer-Wolff transformation that such models naturally appear as effective Hamiltonians in the large electric field limit of the interacting Wannier-Stark problem, and comment on connections of our work with the phenomenon of Bloch many-body localization.
We show that the combination of charge and dipole conservation-characteristic of fracton systemsleads to an extensive fragmentation of the Hilbert space, which, in turn, can lead to a breakdown of thermalization.As a concrete example, we investigate the out-of-equilibrium dynamics of one-dimensional spin-1 models that conserve charge (total S z ) and its associated dipole moment.First, we consider a minimal model including only three-site terms and find that the infinite temperature autocorrelation saturates to a finite value-showcasing nonthermal behavior.The absence of thermalization is identified as a consequence of the strong fragmentation of the Hilbert space into exponentially many invariant subspaces in the local S z basis, arising from the interplay of dipole conservation and local interactions.Second, we extend the model by including four-site terms and find that this perturbation leads to a weak fragmentation: The system still has exponentially many invariant subspaces, but they are no longer sufficient to avoid thermalization for typical initial states.More generally, for any finite range of interactions, the system still exhibits nonthermal eigenstates appearing throughout the entire spectrum.We compare our results to charge and dipole moment-conserving random unitary circuit models for which we reach identical conclusions.
We show that the combination of charge and dipole conservation-characteristic of fracton systems-leads to an extensive fragmentation of the Hilbert space, which, in turn, can lead to a breakdown of thermalization. As a concrete example, we investigate the out-of-equilibrium dynamics of one-dimensional spin-1 models that conserve charge (total S-z) and its associated dipole moment. First, we consider a minimal model including only three-site terms and find that the infinite temperature autocorrelation saturates to a finite value-showcasing nonthermal behavior. The absence of thermalization is identified as a consequence of the strong fragmentation of the Hilbert space into exponentially many invariant subspaces in the local S-z basis, arising from the interplay of dipole conservation and local interactions. Second, we extend the model by including four-site terms and find that this perturbation leads to a weak fragmentation: The system still has exponentially many invariant subspaces, but they are no longer sufficient to avoid thermalization for typical initial states. More generally, for any finite range of interactions, the system still exhibits nonthermal eigenstates appearing throughout the entire spectrum. We compare our results to charge and dipole moment-conserving random unitary circuit models for which we reach identical conclusions.
Understanding the evolution towards thermal equilibrium of an isolated quantum system is at the foundation of statistical mechanics and a subject of interest in such diverse areas as cold atom physics or the quantum mechanics of black holes. Since a pure state can never evolve into a thermal density matrix, the Eigenstate Thermalization Hypothesis (ETH) has been put forward by Deutsch and Srednicki as a way to explain this apparent thermalization, similarly to what the ergodic theorem does in classical mechanics. In this paper this hypothesis is tested numerically. First, it is observed that thermalization happens in a subspace of states (the Krylov subspace) with dimension much smaller than that of the total Hilbert space. We check numerically the validity of ETH in such a subspace, for a system of hard core bosons on a two-dimensional lattice. We then discuss how well the eigenstates of the Hamiltonian projected on the Krylov subspace represent the true eigenstates. This discussion is aided by bringing the projected Hamiltonian to the tridiagonal form and interpreting it as an Anderson localization problem for a finite one-dimensional chain. We also consider thermalization of a subsystem and argue that generation of a large entanglement entropy can lead to a thermal density matrix for the subsystem well before the whole system thermalizes. Finally, we comment on possible implications of ETH in quantum gravity.
Using a Krylov-subspace time evolution algorithm, we simulate the real-time dynamics of translation invariant non-integrable finite spin rings to quite long times with high accuracy. We systematically study the finite-size deviation between the resulting equilibrium state and the thermal state, and we highlight the importance of the energy variance on the deviations. We find that the deviations are well described by the eigenstate thermalization hypothesis, and that the von Neumann entropy correction scaling is the square of the local operator scaling. We reveal also an area law contribution to the relaxed von Neumann entropy, which we connect to the mutual information between the considered subsystem and its immediate environment. We also find that local observables relax towards equilibrium exponentially with a relaxation time scale that grows linearly with system length and is somewhat independent of the local operator, but depends strongly on the energy of the initial state, with the fastest relaxation times found towards one end of the overall energy spectrum. To contrast this behaviour we also study domain wall initial states, which exhibit clear diffusive behaviour, with a Thouless time scaling quadratically with the systems size, leading to a rather precise estimate for the diffusion constant for states in the vicinity of the middle of the energy spectrum.
We investigate a one-dimensional correlated-hopping model of spinless fermions with an East constraint. We first analytically unravel the complete fragmentation structure of this model by labeling each fragment by a unique root configuration and utilizing the transfer matrix method. We show that the growth of the size of each fragment of the model follows the widely studied Dyck sequence and is therefore analytically tractable with the help of Catalan triangles. While the eigenstate thermalization hypothesis (ETH) does not hold within the full Hilbert space which exhibits Poisson statistics of the energy level spacing, an examination of various quantities restricted to the largest fragments shows that a weaker version of the subspace-restricted thermalization holds. This weaker violation of the ETH within the largest fragments is supported by the presence of subspace-restricted quantum many-body scars due to quantum fragmentation. Our analysis further reveals that the filling fraction at which the model has the largest fragment shifts with increasing system sizes from half-filling due to the distinct fragmentation structure. Next, we show that the inclusion of a nearest-neighbor density-density interaction with strength VVV induces a spectral transition within the largest fragment from a weakly ETH-violating phase containing scars to a statistical bubble localized phase as VVV increases. In particular, the $V\ensuremath{\rightarrow}\ensuremath{\infty}$ limit produces an integrable model. We find that the addition of finite-VVV stabilizes the ground state near half-filling while keeping intact the fragmentation structure of the East model. However, this behavior abruptly changes exactly at $V=\ensuremath{\infty}$ due to the emergence of a distinct fragmentation structure; we analytically explore this in a comprehensive manner. The infinite-VVV model has many interesting properties, among which the appearance of the ground state and the largest fragment at two different filling fractions is specially noteworthy. Finally, we propose an experimental setup to realize the infinite-VVV model as a particular limit of a special kind of $t\text{\ensuremath{-}}V$ model with an on-site potential.
While quantum statistical mechanics triumphs in explaining many equilibrium phenomena, there is an increasing focus on going beyond conventional scenarios of thermalization. Traditionally examples of nonthermalizing systems are either integrable or disordered. Recently, examples of translationally invariant physical systems have been discovered whose excited energies avoid thermalization either due to local constraints (whether exact or emergent) or due to higher-form symmetries. In this article, we extend these investigations for the case of 3D U ( 1 ) quantum dimer models, which are lattice gauge theories with finite-dimensional local Hilbert spaces (also generically called quantum link models) with staggered charged static matter. Using a combination of analytical and numerical methods, we uncover a class of athermal states that arise in large winding sectors, when the system is subjected to external electric fields. The polarization of the dynamical fluxes in the direction of applied field traps excitations in 2D planes, while an interplay with the Gauss law constraint in the perpendicular direction causes exotic athermal behavior due to the emergence of new conserved quantities. This causes a geometric fragmentation of the system. We provide analytical arguments showing that the scaling of the number of fragments is exponential in the linear system size, leading to weak fragmentation. Further, we identify sectors which host fractonic excitations with severe mobility restrictions. The unitary evolution of fragments dominated by fractons is qualitatively different from the one dominated by nonfractonic excitations.