Characterizing nonequilibrium dynamics in quantum many-body systems is a challenging frontier of physics. In this Letter, we systematically construct solvable nonintegrable quantum circuits that exhibit exact hidden Markovian subsystem dynamics. This feature thus enables accurately calculating local observables for arbitrary evolution time. Utilizing the influence matrix method, we show that the influence of the time-evolved global system on a finite subsystem can be analytically described by sequential, time-local quantum channels acting on the subsystem with an ancilla of finite Hilbert space dimension. The realization of exact hidden Markovian property is facilitated by a solvable condition on the underlying two-site gates in the quantum circuit. We further present several concrete examples with varying local Hilbert space dimensions to demonstrate our approach.
Characterizing nonequilibrium dynamics in quantum many-body systems is a challenging frontier of physics. In this Letter, we systematically construct solvable nonintegrable quantum circuits that exhibit exact hidden Markovian subsystem dynamics. This feature thus enables accurately calculating local observables for arbitrary evolution time. Utilizing the influence matrix method, we show that the influence of the time-evolved global system on a finite subsystem can be analytically described by sequential, time-local quantum channels acting on the subsystem with an ancilla of finite Hilbert space dimension. The realization of exact hidden Markovian property is facilitated by a solvable condition on the underlying two-site gates in the quantum circuit. We further present several concrete examples with varying local Hilbert space dimensions to demonstrate our approach.
Influence matrices encode the action exerted on local subsystems by the rest of an extended quantum many-body system during their evolution. Thus, knowledge of the influence matrix facilitates computationally efficient simulations of local dynamics. Here we propose a new systematic approach to generating quantum circuits with complex dynamics for which the influence matrices can be written down exactly. In contrast to previous frameworks of this kind, such as dual-unitary circuits, the resulting influence matrices are non-Markovian, exhibiting nontrivial temporal correlations. We explicitly construct a broad family of circuits of this kind, based on dressing free-fermion (matchgate) circuits with appropriately chosen interaction terms. We show that, contrary to previous solvable instances, these circuits produce patterns of correlations that closely resemble that of typical many-body systems. Our approach can be directly interpreted in terms of an error correction scheme where the terms breaking the solvability of the influence matrices play the role of errors.
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 develop an influence matrix bootstrap approach for exactly solving nonequilibrium dynamics. This approach begins with an ansatz consisting of a set of local tensor relations, termed generalized zipper conditions, and uses finite-time data to construct exact solutions of local dynamics that extend to arbitrarily long times. We demonstrate this approach in the "Rule 201" quantum cellular automaton, an integrable Trotterization of the PXP Hamiltonian. 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 orde}: 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.
The dynamics of local observables in a quantum many-body system can be formally described in the language of open systems. The problem is that the bath representing the complement of the local subsystem generally does not allow the common simplifications often crucial for such a framework. Leveraging tensor network calculations and optimization tools from machine learning, we extract and characterize the dynamical maps for single- and two-site subsystems embedded in an infinite quantum Ising chain after a global quench. We consider three paradigmatic regimes: integrable critical, integrable non-critical, and chaotic. For each we find the optimal time-local representation of the subsystem dynamics at different times. We explore the properties of the learned time-dependent Liouvillians and whether they can be used to forecast the long-time dynamics of local observables beyond the times accessible through direct quantum many-body numerical simulation. Our procedure naturally suggests a novel measure of non-Markovianity based on the distance between the quasi-exact dynamical map and the closest CP-divisible form and reveals that criticality leads to the closest Markovian representation at large times.