Opening Krylov Space to Access All-Time Dynamics via Dynamical Symmetries

Solving short- and long-time dynamics of closed quantum many-body systems is one of the main challenges of both atomic and condensed matter physics. For locally interacting closed systems, the dynamics of local observables can always be expanded into (pseudolocal) eigenmodes of the Liouvillian, so-called dynamical symmetries. They come in two classes: transient operators, which decay in time and perpetual operators, which either oscillate forever or stay the same (conservation laws). These operators provide a full characterization of the dynamics of the system. Deriving these operators, apart from a very limited class of models, has not been possible. Here, we present a method to numerically and analytically derive some of these dynamical symmetries in infinite closed systems by introducing a naturally emergent open boundary condition on the Krylov chain. This boundary condition defines a partitioning of the Krylov space into system and environment degrees of freedom, where nonlocal operators make up an effective bath for the local operators. We demonstrate the practicality of the method on some numerical examples and derive analytical results in two idealized cases. Our approach lets us directly relate the operator growth hypothesis to thermalization and exponential decay of observables in chaotic systems and provides a powerful approach for computing notably challenging many-body dynamics.