Slow dynamics from a nested hierarchy of frozen states

We identify the mechanism of slow heterogeneous relaxation in quantum kinetically constrained models (KCMs) in which the potential energy strength is controlled by a coupling parameter. The regime of slow relaxation includes the large-coupling limit. By expanding around that limit, we reveal a of states that remain frozen on time scales determined by powers of the coupling. The classification of such states, together with the evolution of their Krylov complexity, reveals that these time scales are related to the distance between the sites where facilitated dynamics is allowed by the kinetic constraint. While correlations within frozen states relax slowly and exhibit metastable plateaus that persist on time scales set by powers of the coupling parameter, the correlations in the rest of the states decay rapidly. We compute the plateau heights of correlations across all frozen states up to second-order corrections in the inverse coupling. Our results explain slow relaxation in quantum KCMs and elucidate dynamical heterogeneity by relating the relaxation times to the spatial separations between the active regions.