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 auto-correlation saturates to a finite value-showcasing non-thermal 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 non-thermal 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.04266v2 [cond-mat.str-el]

Certain disorder-free Hamiltonians can be non-ergodic due to a strong fragmentation of the Hilbert space into disconnected sectors. Here, we characterize such systems by introducing the notion of 'statistically localized integrals of motion' (SLIOM), whose eigenvalues label the connected components of the Hilbert space. SLIOMs are not spatially localized in the operator sense, but appear localized to sub-extensive regions when their expectation value is taken in typical states with a finite density of particles. We illustrate this general concept on several Hamiltonians, both with and without dipole conservation. Furthermore, we demonstrate that there exist perturbations which destroy these integrals of motion in the bulk of the system, while keeping them on the boundary. This results in statistically localized strong zero modes, leading to infinitely long-lived edge magnetizations along with a thermalizing bulk, constituting the first example of such strong edge modes in a non-integrable model. We also show that in a particular example, these edge modes lead to the appearance of topological string order in a certain subset of highly excited eigenstates. Some of our suggested models can be realized in Rydberg quantum simulators.CONTENTS * These authors contributed equally to this work.

The thermalization of isolated quantum many-body systems is deeply related to fundamental questions of quantum information theory. While integrable or many-body localized systems display non-ergodic behavior due to extensively many conserved quantities, recent theoretical studies have identified a rich variety of more exotic phenomena in between these two extreme limits. The tilted one-dimensional Fermi-Hubbard model, which is readily accessible in experiments with ultracold atoms, emerged as an intriguing playground to study non-ergodic behavior in a clean disorder-free system. While non-ergodic behavior was established theoretically in certain limiting cases, there is no complete understanding of the complex thermalization properties of this model. In this work, we experimentally study the relaxation of an initial charge-density wave and find a remarkably long-lived initial-state memory over a wide range of parameters. Our observations are well reproduced by numerical simulations of a clean system. Using analytical calculations we further provide a detailed microscopic understanding of this behavior, which can be attributed to emergent kinetic constraints.