Heisenberg time evolution under a chaotic many-body Hamiltonian H transforms an initially simple operator into an increasingly complex one, as it spreads over Hilbert space. Krylov complexity, or ‘K-complexity’, quantifies this growth with respect to a special basis, generated by H by successive nested commutators with the operator. In this work we study the evolution of K-complexity in finite-entropy systems for time scales greater than the scrambling time ts> log(S). We prove rigorous bounds on K-complexity as well as the associated Lanczos sequence and, using refined parallelized algorithms, we undertake a detailed numerical study of these quantities in the SYK4 model, which is maximally chaotic, and compare the results with the SYK2 model, which is integrable. While the former saturates the bound, the latter stays exponentially below it. We discuss to what extent this is a generic feature distinguishing between chaotic vs. integrable systems.

Quantum complexity, suitably defined, has been suggested as an important probe of late-time dynamics of black holes, particularly in the context of AdS/CFT. A notion of quantum complexity can be effectively captured by quantifying the spread of an operator in Krylov space as a consequence of time evolution. Complexity is expected to behave differently in chaotic many-body systems, as compared to integrable ones. In this paper we investigate Krylov complexity for the case of interacting integrable models at finite size and find that complexity saturation is suppressed as compared to chaotic systems. We associate this behavior with a novel localization phenomenon on the Krylov chain by mapping the theory of complexity growth and spread to an Anderson localization hopping model with off-diagonal disorder, and find that localization is enhanced in the integrable case due to a stronger disorder in the hopping amplitudes, inducing an effective suppression of Krylov complexity. We demonstrate this behavior for an interacting integrable model, the XXZ spin chain, and show that the same behavior results from a phenomenological model that we define: this model captures the essential features of our analysis and is able to reproduce the behaviors we observe for chaotic and integrable systems via an adjustable disorder parameter.

We apply a notion of quantum complexity, called “Krylov complexity”, to study the evolution of systems from integrability to chaos. For this purpose we investigate the integrable XXZ spin chain, enriched with an integrability breaking deformation that allows one to interpolate between integrable and chaotic behavior. K-complexity can act as a probe of the integrable or chaotic nature of the underlying system via its late-time saturation value that is suppressed in the integrable phase and increases as the system is driven to the chaotic phase. We furthermore ascribe the (under-)saturation of the late-time bound to the amount of disorder present in the Lanczos sequence, by mapping the complexity evolution to an auxiliary off-diagonal Anderson hopping model. We compare the late-time saturation of K-complexity in the chaotic phase with that of random matrix ensembles and find that the chaotic system indeed approaches the RMT behavior in the appropriate symmetry class. We investigate the dependence of the results on the two key ingredients of K-complexity: the dynamics of the Hamiltonian and the character of the operator whose time dependence is followed.