Integrable Trotterization: Local Conservation Laws and Boundary Driving

Quantum Turnstiles for Robust Measurement of Full Counting Statistics

We present a scalable protocol for measuring full counting statistics (FCS) in experiments or tensor-network simulations. In this method, an ancilla in the middle of the system acts as a turnstile, with its phase keeping track of the time-integrated particle flux. Unlike quantum gas microscopy, the turnstile protocol faithfully captures FCS starting from number-indefinite initial states or in the presence of noisy dynamics. In addition, by mapping the FCS onto a single-body observable, it allows for stable numerical calculations of FCS using approximate tensor-network methods. We demonstrate the wide-ranging utility of this approach by computing the FCS of the transferred magnetization in a Floquet Heisenberg spin chain, as studied in a recent experiment with superconducting qubits, as well as the FCS of charge transfer in random circuits.

Strong Zero Modes in Integrable Quantum Circuits

It is a classic result that certain interacting integrable spin chains host robust edge modes known as strong zero modes (SZMs). In this Letter, we extend this result to the Floquet setting of local quantum circuits, focusing on a prototypical model providing an integrable Trotterization for the evolution of the XXZ Heisenberg spin chain. By exploiting the algebraic structures of integrability, we show that an exact SZM operator can be constructed for these integrable quantum circuits in certain regions of parameter space. Our construction, which recovers a well-known result by Paul Fendley in the continuous-time limit, relies on a set of commuting transfer matrices known from integrability, and allows us to easily prove important properties of the SZM, including normalizabilty. Our approach is different from previous methods and could be of independent interest even in the Hamiltonian setting. Our predictions, which are corroborated by numerical simulations of infinite-temperature autocorrelation functions, are potentially interesting for implementations of the XXZ quantum circuit on available quantum platforms.

Dynamics of magnetization at infinite temperature in a Heisenberg spin chain

Understanding universal aspects of quantum dynamics is an unresolved problem in statistical mechanics. In particular, the spin dynamics of the one-dimensional Heisenberg model were conjectured as to belong to the Kardar-Parisi-Zhang (KPZ) universality class based on the scaling of the infinite-temperature spin-spin correlation function. In a chain of 46 superconducting qubits, we studied the probability distribution of the magnetization transferred across the chain's center, [Formula: see text]. The first two moments of [Formula: see text] show superdiffusive behavior, a hallmark of KPZ universality. However, the third and fourth moments ruled out the KPZ conjecture and allow for evaluating other theories. Our results highlight the importance of studying higher moments in determining dynamic universality classes and provide insights into universal behavior in quantum systems.

Scrambling Is Necessary but Not Sufficient for Chaos

We show that out-of-time-order correlators (OTOCs) constitute a probe for local-operator entanglement (LOE). There is strong evidence that a volumetric growth of LOE is a faithful dynamical indicator of quantum chaos, while OTOC decay corresponds to operator scrambling, often conflated with chaos. We show that rapid OTOC decay is a necessary but not sufficient condition for linear (chaotic) growth of the LOE entropy. We analytically support our results through wide classes of local-circuit models of many-body dynamics, including both integrable and nonintegrable dual-unitary circuits. We show sufficient conditions under which local dynamics leads to an equivalence of scrambling and chaos.

Integrable Digital Quantum Simulation: Generalized Gibbs Ensembles and Trotter Transitions

The Trotter-Suzuki decomposition is a promising avenue for digital quantum simulation (DQS), approximating continuous-time dynamics by discrete Trotter steps of duration τ. Recent work suggested that DQS is typically characterized by a sharp Trotter transition: when τ is increased beyond a threshold value, approximation errors become uncontrolled at large times due to the onset of quantum chaos. Here, we contrast this picture with the case of integrable DQS. We focus on a simple quench from a spin-wave state in the prototypical XXZ Heisenberg spin chain, and study its integrable Trotterized evolution as a function of τ. Because of its exact local conservation laws, the system does not heat up to infinite temperature and the late-time properties of the dynamics are captured by a discrete generalized Gibbs ensemble (dGGE). By means of exact calculations we find that, for small τ, the dGGE depends analytically on the Trotter step, implying that discretization errors remain bounded even at infinite times. Conversely, the dGGE changes abruptly at a threshold value τ_{th}, signaling a novel type of Trotter transition. We show that the latter can be detected locally, as it is associated with the appearance of a nonzero staggered magnetization with a subtle dependence on τ. We highlight the differences between continuous and discrete GGEs, suggesting the latter as novel interesting nonequilibrium states exclusive to digital platforms.

On Two Non-Ergodic Reversible Cellular Automata, One Classical, the Other Quantum

We propose and discuss two variants of kinetic particle models-cellular automata in 1 + 1 dimensions-that have some appeal due to their simplicity and intriguing properties, which could warrant further research and applications. The first model is a deterministic and reversible automaton describing two species of quasiparticles: stable massless matter particles moving with velocity ±1 and unstable standing (zero velocity) field particles. We discuss two distinct continuity equations for three conserved charges of the model. While the first two charges and the corresponding currents have support of three lattice sites and represent a lattice analogue of the conserved energy-momentum tensor, we find an additional conserved charge and current with support of nine sites, implying non-ergodic behaviour and potentially signalling integrability of the model with a highly nested R-matrix structure. The second model represents a quantum (or stochastic) deformation of a recently introduced and studied charged hardpoint lattice gas, where particles of different binary charge (±1) and binary velocity (±1) can nontrivially mix upon elastic collisional scattering. We show that while the unitary evolution rule of this model does not satisfy the full Yang-Baxter equation, it still satisfies an intriguing related identity which gives birth to an infinite set of local conserved operators, the so-called glider operators.

Anomalous transport from hot quasiparticles in interacting spin chains

Many experimentally relevant quantum spin chains are approximately integrable, and support long-lived quasiparticle excitations. A canonical example of integrable model of quantum magnetism is the XXZ spin chain, for which energy spreads ballistically, but, surprisingly, spin transport can be diffusive or superdiffusive. We review the transport properties of this model using an intuitive quasiparticle picture that relies on the recently introduced framework of generalized hydrodynamics. We discuss how anomalous linear response properties emerge from hierarchies of quasiparticles both in integrable and near-integrable limits, with an emphasis on the role of hydrodynamic fluctuations. We also comment on recent developments including non-linear response, full-counting statistics and far-from-equilibrium transport. We provide an overview of recent numerical and experimental results on transport in XXZ spin chains.

Universal Kardar-Parisi-Zhang Dynamics in Integrable Quantum Systems

Although the Bethe ansatz solution of the spin-1/2 Heisenberg model dates back nearly a century, the anomalous nature of its high-temperature transport dynamics has only recently been uncovered. Indeed, numerical and experimental observations have demonstrated that spin transport in this paradigmatic model falls into the Kardar-Parisi-Zhang (KPZ) universality class. This has inspired the significantly stronger conjecture that KPZ dynamics, in fact, occur in all integrable spin chains with non-Abelian symmetry. Here, we provide extensive numerical evidence affirming this conjecture. Moreover, we observe that KPZ transport is even more generic, arising in both supersymmetric and periodically driven models. Motivated by recent advances in the realization of SU(N)-symmetric spin models in alkaline-earth-based optical lattice experiments, we propose and analyze a protocol to directly investigate the KPZ scaling function in such systems.

Formation of robust bound states of interacting microwave photons

Systems of correlated particles appear in many fields of modern science and represent some of the most intractable computational problems in nature. The computational challenge in these systems arises when interactions become comparable to other energy scales, which makes the state of each particle depend on all other particles1. The lack of general solutions for the three-body problem and acceptable theory for strongly correlated electrons shows that our understanding of correlated systems fades when the particle number or the interaction strength increases. One of the hallmarks of interacting systems is the formation of multiparticle bound states2-9. Here we develop a high-fidelity parameterizable fSim gate and implement the periodic quantum circuit of the spin-½ XXZ model in a ring of 24 superconducting qubits. We study the propagation of these excitations and observe their bound nature for up to five photons. We devise a phase-sensitive method for constructing the few-body spectrum of the bound states and extract their pseudo-charge by introducing a synthetic flux. By introducing interactions between the ring and additional qubits, we observe an unexpected resilience of the bound states to integrability breaking. This finding goes against the idea that bound states in non-integrable systems are unstable when their energies overlap with the continuum spectrum. Our work provides experimental evidence for bound states of interacting photons and discovers their stability beyond the integrability limit.

Absence of Superdiffusion in Certain Random Spin Models

The dynamics of spin at finite temperature in the spin-1/2 Heisenberg chain was found to be superdiffusive in numerous recent numerical and experimental studies. Theoretical approaches to this problem have emphasized the role of nonabelian SU(2) symmetry as well as integrability, but the associated methods cannot be readily applied when integrability is broken. We examine spin transport in a spin-1/2 chain in which the exchange couplings fluctuate in space and time around a nonzero mean J, a model introduced by De Nardis et al. [Phys. Rev. Lett. 127, 057201 (2021).PRLTAO0031-900710.1103/PhysRevLett.127.057201]. We show that operator dynamics in the strong noise limit at infinite temperature can be analyzed using conventional perturbation theory as an expansion in J. We find that regular diffusion persists at long times, albeit with an enhanced diffusion constant. The finite time spin dynamics is analyzed and compared with matrix product operator simulations.

Temporal Entanglement, Quasiparticles, and the Role of Interactions

In quantum many-body dynamics admitting a description in terms of noninteracting quasiparticles, the Feynman-Vernon influence matrix (IM), encoding the effect of the system on the evolution of its local subsystems, can be analyzed exactly. For discrete dynamics, the temporal entanglement (TE) of the corresponding IM satisfies an area law, suggesting the possibility of an efficient representation of the IM in terms of matrix-product states. A natural question is whether integrable interactions, preserving stable quasiparticles, affect the behavior of the TE. While a simple semiclassical picture suggests a sublinear growth in time, one can wonder whether interactions may lead to violations of the area law. We address this problem by analyzing quantum quenches in a family of discrete integrable dynamics corresponding to the real-time Trotterization of the interacting XXZ Heisenberg model. By means of an analytical solution at the dual-unitary point and numerical calculations for generic values of the system parameters, we provide evidence that, away from the noninteracting limit, the TE displays a logarithmic growth in time, thus violating the area law. Our findings highlight the nontrivial role of interactions, and raise interesting questions on the possibility to efficiently simulate the local dynamics of interacting integrable systems.

Macroscopic Effects of Localized Measurements in Jammed States of Quantum Spin Chains

A quantum jammed state can be seen as a state where the phase space available to particles shrinks to zero, an interpretation quite accurate in integrable systems, where stable quasiparticles scatter elastically. We consider the integrable dual folded XXZ model, which is equivalent to the XXZ model in the limit of large anisotropy. We perform a jamming-breaking localized measurement in a jammed state. We find that jamming is locally restored, but local observables exhibit nontrivial time evolution on macroscopic, ballistic scales, without ever relaxing back to their initial values.

Integrable spin chains and cellular automata with medium-range interaction

We study integrable spin chains and quantum and classical cellular automata with interaction range ℓ≥3. This is a family of integrable models for which there was no general theory so far. We develop an algebraic framework for such models, generalizing known methods from nearest-neighbor interacting chains. This leads to a new integrability condition for medium-range Hamiltonians, which can be used to classify such models. A partial classification is performed in specific cases, including U(1)-symmetric three-site interacting models, and Hamiltonians that are relevant for interaction-round-a-face models. We find a number of models which appear to be new. As an application we consider quantum brickwork circuits of various types, including those that can accommodate the classical elementary cellular automata on light cone lattices. In this family we find that the so-called Rule150 and Rule105 models are Yang-Baxter integrable with three-site interactions. We present integrable quantum deformations of these models, and derive a set of local conserved charges for them. For the famous Rule54 model we find that it does not belong to the family of integrable three-site models, but we cannot exclude Yang-Baxter integrability with longer interaction ranges.

Chaos and Ergodicity in Extended Quantum Systems with Noisy Driving

We study the time-evolution operator in a family of local quantum circuits with random fields in a fixed direction. We argue that the presence of quantum chaos implies that at large times the time-evolution operator becomes effectively a random matrix in the many-body Hilbert space. To quantify this phenomenon, we compute analytically the squared magnitude of the trace of the evolution operator-the generalized spectral form factor-and compare it with the prediction of random matrix theory. We show that for the systems under consideration, the generalized spectral form factor can be expressed in terms of dynamical correlation functions of local observables in the infinite temperature state, linking chaotic and ergodic properties of the systems. This also provides a connection between the many-body Thouless time τ_{th}-the time at which the generalized spectral form factor starts following the random matrix theory prediction-and the conservation laws of the system. Moreover, we explain different scalings of τ_{th} with the system size observed for systems with and without the conservation laws.

Ergodic and Nonergodic Dual-Unitary Quantum Circuits with Arbitrary Local Hilbert Space Dimension

Dual-unitary quantum circuits can be used to construct 1+1 dimensional lattice models for which dynamical correlations of local observables can be explicitly calculated. We show how to analytically construct classes of dual-unitary circuits with any desired level of (non-)ergodicity for any dimension of the local Hilbert space, and present analytical results for thermalization to an infinite-temperature Gibbs state (ergodic) and a generalized Gibbs ensemble (nonergodic). It is shown how a tunable ergodicity-inducing perturbation can be added to a nonergodic circuit without breaking dual unitarity, leading to the appearance of prethermalization plateaux for local observables.

Undular Diffusion in Nonlinear Sigma Models

We discuss general features of charge transport in nonrelativistic classical field theories invariant under non-Abelian unitary Lie groups by examining the full structure of two-point dynamical correlation functions in grand-canonical ensembles at finite charge densities (polarized ensembles). Upon explicit breaking of non-Abelian symmetry, two distinct transport laws characterized by dynamical exponent z=2 arise. While in the unbroken symmetry sector, the Cartan fields exhibit normal diffusion, the transversal sectors governed by the nonlinear analogs of Goldstone modes disclose an unconventional law of diffusion, characterized by a complex diffusion constant and undulating patterns in the spatiotemporal correlation profiles. In the limit of strong polarization, one retrieves the imaginary-time diffusion for uncoupled linear Goldstone modes, whereas for weak polarizations the imaginary component of the diffusion constant becomes small. In models of higher rank symmetry, we prove absence of dynamical correlations among distinct transversal sectors.

Spectral Statistics and Many-Body Quantum Chaos with Conserved Charge

We investigate spectral statistics in spatially extended, chaotic many-body quantum systems with a conserved charge. We compute the spectral form factor K(t) analytically for a minimal Floquet circuit model that has a U(1) symmetry encoded via spin-1/2 degrees of freedom. Averaging over an ensemble of realizations, we relate K(t) to a partition function for the spins, given by a Trotterization of the spin-1/2 Heisenberg ferromagnet. Using Bethe ansatz techniques, we extract the "Thouless time" t_{Th} demarcating the extent of random matrix behavior, and find scaling behavior governed by diffusion for K(t) at t≲t_{Th}. We also report numerical results for K(t) in a generic Floquet spin model, which are consistent with these analytic predictions.

Exact Correlation Functions for Dual-Unitary Lattice Models in 1+1 Dimensions

We consider a class of quantum lattice models in 1+1 dimensions represented as local quantum circuits that enjoy a particular dual-unitarity property. In essence, this property ensures that both the evolution in time and that in space are given in terms of unitary transfer matrices. We show that for this class of circuits, generically nonintegrable, one can compute explicitly all dynamical correlations of local observables. Our result is exact, nonpertubative, and holds for any dimension d of the local Hilbert space. In the minimal case of qubits (d=2) we also present a classification of all dual-unitary circuits which allows us to single out a number of distinct classes for the behavior of the dynamical correlations. We find noninteracting classes, where all correlations are preserved, the ergodic and mixing one, where all correlations decay, and, interestingly, also classes that are both interacting and nonergodic.

Integrable Many-Body Quantum Floquet-Thouless Pumps

We construct an interacting integrable Floquet model featuring quasiparticle excitations with topologically nontrivial chiral dispersion. This model is a fully quantum generalization of an integrable classical cellular automaton. We write down and solve the Bethe equations for the generalized quantum model and show that these take on a particularly simple form that allows for an exact solution: essentially, the quasiparticles behave like interacting hard rods. The generalized thermodynamics and hydrodynamics of this model follow directly, providing an exact description of interacting chiral particles in the thermodynamic limit. Although the model is interacting, its unusually simple structure allows us to construct operators that spread with no butterfly effect; this construction does not seem possible in other interacting integrable systems. This model exemplifies a new class of exactly solvable, interacting quantum systems specific to the Floquet setting.

Kardar-Parisi-Zhang Physics in the Quantum Heisenberg Magnet

Equilibrium spatiotemporal correlation functions are central to understanding weak nonequilibrium physics. In certain local one-dimensional classical systems with three conservation laws they show universal features. Namely, fluctuations around ballistically propagating sound modes can be described by the celebrated Kardar-Parisi-Zhang (KPZ) universality class. Can such a universality class be found also in quantum systems? By unambiguously demonstrating that the KPZ scaling function describes magnetization dynamics in the SU(2) symmetric Heisenberg spin chain we show, for the first time, that this is so. We achieve that by introducing new theoretical and numerical tools, and make a puzzling observation that the conservation of energy does not seem to matter for the KPZ physics.

Ballistic Spin Transport in a Periodically Driven Integrable Quantum System

We demonstrate ballistic spin transport of an integrable unitary quantum circuit, which can be understood either as a paradigm of an integrable periodically driven (Floquet) spin chain, or as a Trotterized anisotropic (XXZ) Heisenberg spin-1/2 model. We construct an analytic family of quasilocal conservation laws that break the spin-reversal symmetry and compute a lower bound on the spin Drude weight, which is found to be a fractal function of the anisotropy parameter. Extensive numerical simulations of spin transport suggest that this fractal lower bound is in fact tight.