We study higher dimensional models with symmetric correlated hoppings, which generalize a one-dimensional model introduced in the context of dipole-conserving dynamics. We prove rigorously that whenever the local configuration space takes its smallest non-trivial value, these models exhibit localized behavior due to fragmentation, in any dimension. For the same class of models, we then construct a hierarchy of conserved quantities that are power-law localized at the boundary of the system with increasing powers. Combining these with Mazur’s bound, we prove that boundary correlations are infinitely long lived, even when the bulk is not localized. We use our results to construct quantum Hamiltonians that exhibit the analogues of strong zero modes in two and higher dimensions.
We study higher dimensional models with symmetric correlated hoppings, which generalize a one-dimensional model introduced in the context of dipole-conserving dynamics. We prove rigorously that whenever the local configuration space takes its smallest non-trivial value, these models exhibit localized behavior due to fragmentation, in any dimension. For the same class of models, we then construct a hierarchy of conserved quantities that are power-law localized at the boundary of the system with increasing powers. Combining these with Mazur’s bound, we prove that boundary correlations are infinitely long lived, even when the bulk is not localized. We use our results to construct quantum Hamiltonians that exhibit the analogues of strong zero modes in two and higher dimensions.
By mapping the calculation of Mazur bounds to the enumeration of walks on fractal structures, we present exact bounds on the late-time behavior of spin autocorrelation functions in models exhibiting pair-flip dynamics and more general p-flip dynamics. While the pair-flip model is known to exhibit strong Hilbert space fragmentation, the effect of its nontrivial conservation laws on autocorrelation functions has, thus far, only been calculated numerically, which has led to incorrect conclusions about their thermodynamic behavior. Here, using exact results, we prove that infinite-temperature autocorrelation functions exhibit infinite coherence times at the boundary, and that bulk Mazur bounds decay asymptotically as 1/ \sqrt{L} L , rather than 1/L, as had previously been thought. This result implies that the nontrivial conserved operators implied by p-flip dynamics have an important qualitative impact on bulk thermalization properties beyond the constraints imposed by the simple global symmetries of the models.
Hilbert space fragmentation is an ergodicity-breaking phenomenon, in which the Hamiltonian shatters into exponentially many dynamically disconnected sectors. In many fragmented systems, these sectors can be labeled by statistically localized integrals of motion, which are nonlocal operators. We study the paradigmatic nearest-neighbor pair hopping model exhibiting the so-called strong fragmentation. We show that this model hosts local integrals of motion (LIOMs), which correspond to frozen density modes with long wavelengths. The latter modes become subdiffusive when longer-range pair hoppings are allowed. Finally, we make a connection with a tilted (Stark) chain. Contrary to the dipole-conserving effective models, the tilted chain is shown to support either a Hamiltonian or dipole moment as an LIOM. Numerical results are obtained from a numerical algorithm, in which finding LIOMs is reduced to a data compression problem.
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.
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.
Kinetically constrained models were originally introduced to capture slow relaxation in glassy systems, where dynamics are hindered by local constraints instead of energy barriers. Their quantum counterparts have recently drawn attention for exhibiting highly degenerate eigenstates at zero energy—known as zero modes—stemming from chiral symmetry. Yet, the structure and implications of these zero modes remain poorly understood. In this work, we focus on the properties of the zero mode subspace in quantum kinetically constrained models with a U(1) particle-conservation symmetry. We use the U(1) East, which lacks inversion symmetry, and the inversion-symmetric U(1) East-West models to illustrate our two main results. First, we observe that the simultaneous presence of constraints and chiral symmetry generally leads to a parametric increase in the number of zero modes due to the fragmentation of the many-body Hilbert space into disconnected sectors. Second, we generalize the concept of compact localized states from single-particle physics and introduce the notion of collective bound states, a special kind of nonergodic eigenstates that are robust to enlarging the system size. We formulate sufficient criteria for their existence, arguing that the degenerate zero mode subspace plays a central role, and demonstrate bound states in both example models and in a two-dimensional model, the U(1) North-East, and in the pairflip model, a system without particle conservation. Our results motivate a systematic study of bound states and their relation to ergodicity breaking, transport, and other properties of quantum kinetically constrained models.
Particle-hole symmetric localization problems are studied in low dimensions with the emphasis on understanding the origin of the observed singularities in the low-energy density of states. Pure imaginary random hopping and bipartite random hopping problems are considered. In quasi-1D geometry, such systems generically exhibit localized behavior, and are also found to have power-law density of states with nonuniversal dynamical exponents. Griffiths-like mechanism is identified that generates the singularly many low-energy states. Distinct localized phases are separated by delocalized critical points, and at these points the density of states has strong Dyson singularity near the band center. Similar mechanism is argued to be operative in all particle-hole symmetric lo-calized phases in any dimension. Higher-dimensional systems can also exhibit true metallic or critical delocalized phases, and one such system—the 2D bipartite ran-dom hopping problem—is studied in great detail. General such system is known to