Studies of integrable systems have profoundly deepened the fundamental understanding of quantum many-body physics. While equilibrium properties such as ground states and thermodynamics can often be characterized efficiently, accurately characterizing nonequilibrium integrable dynamics remains a significant challenge. Here, we address this problem in the "Rule 201" quantum cellular automaton, an integrable Trotterization of the PXP Hamiltonian. Using the tensor-network approach of the influence matrix, we develop local conditions called generalized zipper conditions that allow exact solutions of local dynamics. We also introduce a numerical bootstrap method for solving influence matrices with finite but relatively large bond dimensions. This uncovers a rich landscape of nonequilibrium behavior exhibiting initial-state dependence. As an example, we investigate the fate of persistent oscillating dynamics under local non-integrable perturbations, and present analytical results for non-thermal relaxation constrained by conservation laws. We also obtain numerically exact results for entanglement growth across a broad class of initial states. Furthermore, from an information-theoretic perspective, we identify a refined structure of multitime correlations termed the hidden Markov order: the memory encoded in the dynamics separates into finite-length and long-range distributed components, which becomes transparent in an exact split-index matrix-product-state representation of the influence matrix. Our approach enables unified investigations of nonthermalizing and thermalizing regimes of nonequilibrium dynamics within a single analytically tractable model, and can be tested experimentally in state-of-the-art quantum simulators such as Rydberg atom arrays.
Studies of integrable systems have profoundly deepened the fundamental understanding of quantum many-body physics. While equilibrium properties such as ground states and thermodynamics can often be characterized efficiently, accurately characterizing nonequilibrium integrable dynamics remains a significant challenge. Here, we address this problem in the "Rule 201" quantum cellular automaton, an integrable Trotterization of the PXP Hamiltonian. Using the tensor-network approach of the influence matrix, we develop local conditions called generalized zipper conditions that allow exact solutions of local dynamics. We also introduce a numerical bootstrap method for solving influence matrices with finite but relatively large bond dimensions. This uncovers a rich landscape of nonequilibrium behavior exhibiting initial-state dependence. As an example, we investigate the fate of persistent oscillating dynamics under local non-integrable perturbations, and present analytical results for non-thermal relaxation constrained by conservation laws. We also obtain numerically exact results for entanglement growth across a broad class of initial states. Furthermore, from an information-theoretic perspective, we identify a refined structure of multitime correlations termed the hidden Markov order: the memory encoded in the dynamics separates into finite-length and long-range distributed components, which becomes transparent in an exact split-index matrix-product-state representation of the influence matrix. Our approach enables unified investigations of nonthermalizing and thermalizing regimes of nonequilibrium dynamics within a single analytically tractable model, and can be tested experimentally in state-of-the-art quantum simulators such as Rydberg atom arrays.
Recent experimental and theoretical works have made much progress toward understanding nonequilibrium phenomena in thermalizing systems, which act as thermal baths for their small subsystems, and many-body localized systems, which fail to do so.The description of time evolution in many-body systems is generally challenging due to the dynamical generation of quantum entanglement.In this work, we introduce an approach to study quantum many-body dynamics, inspired by the Feynman-Vernon influence functional.Focusing on a family of interacting, Floquet spin chains, we consider a Keldysh path-integral description of the dynamics.The central object in our approach is the influence matrix, which describes the effect of the system on the dynamics of a local subsystem.For translationally invariant models, we formulate a self-consistency equation for the influence matrix.For certain special values of the model parameters, we obtain an exact solution which represents a perfect dephaser (PD).Physically, a PD corresponds to a many-body system that acts as a perfectly Markovian bath on itself: at each period, it measures every spin.For the models considered here, we establish that PD points include dual-unitary circuits investigated in recent works.In the vicinity of PD points, the system is not perfectly Markovian, but rather acts as a bath with a short memory time.In this case, we demonstrate that the self-consistency equation can be solved using matrix-product states (MPS) methods, as the influence matrix temporal entanglement is low.A combination of analytical insights and MPS computations allows us to characterize the structure of the influence matrix in terms of an effective "statistical-mechanics" description.We finally illustrate the predictive power of this description by analytically computing how quickly an embedded impurity spin thermalizes.The influence matrix approach formulated here provides an intuitive view of the quantum many-body dynamics problem, opening a path to constructing models of thermalizing dynamics that are solvable or can be efficiently treated by MPS-based methods and to further characterizing quantum ergodicity or lack thereof.
The dynamics of local subsystems in a thermodynamically large quantum many-body system can be understood as effectively open as the system produces its own effective bath. The action of this bath can be characterised in terms of the so-called influence matrices. In generic situations, the complexity of these objects grows unfavourably with time, however, there exist solvable cases where influence matrices can be characterised exactly even in the presence of non-trivial interactions. Here we show that Rule 201, a deterministic version of the Floquet-PXP model, is one of these solvable instances. Indeed, it admits influence matrices given by a finite-dimensional matrix-product operator (MPO) that solves a finite set of algebraic conditions. We provide the solution, and characterise multi-time autocorrelation functions.
Exactly solvable models in quantum many-body dynamics provide valuable insights into many interesting physical phenomena and serve as platforms for rigorously investigating fundamental theoretical questions. Nevertheless, they are extremely rare and existing solvable models and solution techniques have serious limitations. In this paper, we present a family of exactly solvable unitary circuits that model quantum many-body dynamics in discrete space and time. Unlike many previous solvable models, one can exactly compute the full quantum dynamics initialized from any matrix product state in this family of models. The time evolution of local observables and correlations, the linear growth of Rényi entanglement entropy, spatiotemporal correlations, and out-of-time-order correlations are all exactly computable. A key property of these models enabling the exact solution is that any time-evolved local operator is an exact matrix product operator with a finite bond dimension, even at arbitrarily long times, which we prove using the underlying C * -(weak) Hopf algebra structure along with tensor network techniques. We lay down the general framework for construction and solution of this family of models, and give several explicit examples. In particular, we study in detail a model constructed out of a C * -weak Hopf algebra that is very close to a Floquet version of the PXP model, and the exact results we obtain may shed light on the phenomenon of quantum many-body scars, and more generally, Floquet quantum dynamics in constrained systems.
Quantum scars are special eigenstates of many-body systems that evade thermalization. They were first discovered in the PXP model, a well-known effective description of Rydberg atom arrays. Despite significant theoretical efforts, the fundamental origin of PXP scars remains elusive. By investigating the discretized dynamics of the PXP model as a function of the Trotter step τ, we uncover a remarkable correspondence between the zero- and two-particle eigenstates of the integrable Floquet-PXP cellular automaton at τ=π/2 and the PXP many-body scars of the time-continuous limit. Specifically, we demonstrate that PXP scars are adiabatically connected to the eigenstates of the τ=π/2 Floquet operator. Building on this result, we propose a protocol for achieving high-fidelity preparation of PXP scars in Rydberg atom experiments.
We study the out-of-equilibrium dynamics of the quantum cellular automaton known as "Rule 54." For a class of low-entangled initial states, we provide an analytic description of the effect of the global evolution on finite subsystems in terms of simple quantum channels, which gives access to the full thermalization dynamics at the microscopic level. As an example, we provide analytic formulas for the evolution of local observables and Rényi entropies. We show that, in contrast to other known examples of exactly solvable quantum circuits, Rule 54 does not behave as a simple Markovian bath on its own parts, and displays typical nonequilibrium features of interacting integrable many-body quantum systems such as finite relaxation rate and interaction-induced dressing effects. Our study provides a rare example where the full thermalization dynamics can be solved exactly at the microscopic level.
Feynman–Vernon influence functional (IF) was originally introduced to describe the effect of a quantum environment on the dynamics of an open quantum system. We apply the IF approach to describe quantum many-body dynamics in isolated spin systems, viewing the system as an environment for its local subsystems. While the IF can be computed exactly only in certain many-body models, it generally satisfies a self-consistency equation, provided the system, or an ensemble of systems, are translationally invariant. We view the IF as a fictitious wavefunction in the temporal domain, and approximate it using matrix-product states (MPS). This approach is efficient provided the temporal entanglement of the IF is sufficiently low. We illustrate the broad applicability of the IF approach by analyzing several models that exhibit a range of dynamical behaviors, from thermalizing to many-body localized. In particular, we study the non-equilibrium dynamics in the quantum Ising model in both Floquet and Hamiltonian settings. We find that temporal entanglement entropy may be significantly lower compared to the spatial entanglement and analyze the IF in the continuous-time limit. We simulate the thermodynamic-limit evolution of local observables in various regimes, including the relaxation of impurities embedded in an infinite-temperature chain, and the long-lived oscillatory dynamics of the magnetization associated with the confinement of excitations. Furthermore, by incorporating disorder-averaging into the formalism, we analyze discrete time-crystalline response using the IF of a bond-disordered kicked Ising chain. In this case, we find that the temporal entanglement entropy scales logarithmically with evolution time. The IF approach offers a new lens on many-body non-equilibrium phenomena, both in ergodic and non-ergodic regimes, connecting the theory of open quantum systems to quantum statistical physics.