Localized Dynamics in the Floquet Quantum East Model.

We introduce and study the discrete-time version of the quantum East model, an interacting quantum spin chain inspired by simple kinetically constrained models of classical glasses. Previous work has established that its continuous-time counterpart displays a disorder-free localization transition signaled by the appearance of an exponentially large (in the volume) family of nonthermal, localized eigenstates. Here we combine analytical and numerical approaches to show that (i) the transition persists for discrete times, in fact, it is present for any finite value of the time step apart from a zero measure set; (ii) it is directly detected by following the nonequilibrium dynamics of the fully polarized state. Our findings imply that the transition is currently observable in state-of-the-art platforms for digital quantum simulation.

Ergodicity Breaking and Deviation from Eigenstate Thermalization in Relativistic Quantum Field Theory.

The validity of the ergodic hypothesis in quantum systems can be rephrased in the form of the eigenstate thermalization hypothesis (ETH), a set of statistical properties for the matrix elements of local observables in energy eigenstates, which is expected to hold in any ergodic system. We test the ETH in a nonintegrable model of relativistic quantum field theory (QFT) using the numerical method of Hamiltonian truncation in combination with analytical arguments based on Lorentz symmetry and renormalization group theory. We find that there is an infinite sequence of eigenstates with the characteristics of quantum many-body scars-that is, exceptional eigenstates with observable expectation values that lie far from thermal values-and we show that these states are one-quasiparticle states. We argue that in the thermodynamic limit the eigenstates cover the entire area between two diverging lines: the line of one-quasiparticle states, whose direction is dictated by relativistic kinematics, and the thermal average line. Our results suggest that the strong version of the ETH is violated in any relativistic QFT whose spectrum admits a quasiparticle description.

Dynamical Criticality of Magnetization Transfer in Integrable Spin Chains.

Recent studies have found that fluctuations of magnetization transfer in integrable spin chains violate the central limit property. Here, we revisit the problem of anomalous counting statistics in the Landau-Lifshitz field theory by specializing to two distinct anomalous regimes featuring a dynamical critical point. By performing optimized numerical simulations using an integrable space-time discretization, we extract the algebraic growth exponents of time-dependent cumulants which attain their threshold values. The distinctly non-Gaussian statistics of magnetization transfer in the easy-axis regime is found to converge toward the universal distribution of charged single-file systems. At the isotropic point, we infer a weakly non-Gaussian distribution, corroborating the view that superdiffusive spin transport in integrable spin chains does not belong to any known dynamical universality class.

Reply to "Comment on 'Exact large deviation statistics and trajectory phase transition of a deterministic boundary driven cellular automaton' ".

We reply to Whitelam's Comment [Phys. Rev. E 108, 036105 (2023)2470-004510.1103/PhysRevE.108.036105] on our paper [Phys. Rev. E 100, 020103(R) (2019)2470-004510.1103/PhysRevE.100.020103] where we compute the exact large deviation (LD) statistics of a wide class of observables in the rule 54 cellular automaton. Using some heuristic arguments, Whitelam states that despite the fact that the LD functions we compute display singular behavior, this is not indicative of a LD phase transition or of dynamical phase coexistence. Here, we refute this observation and confirm that the (standard) interpretation of our exact results stands.

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.

Spectral form factor in a minimal bosonic model of many-body quantum chaos.

We study spectral form factor in periodically kicked bosonic chains. We consider a family of models where a Hamiltonian with the terms diagonal in the Fock space basis, including random chemical potentials and pairwise interactions, is kicked periodically by another Hamiltonian with nearest-neighbor hopping and pairing terms. We show that, for intermediate-range interactions, the random phase approximation can be used to rewrite the spectral form factor in terms of a bistochastic many-body process generated by an effective bosonic Hamiltonian. In the particle-number conserving case, i.e., when pairing terms are absent, the effective Hamiltonian has a non-Abelian SU(1,1) symmetry, resulting in universal quadratic scaling of the Thouless time with the system size, irrespective of the particle number. This is a consequence of degenerate symmetry multiplets of the subleading eigenvalue of the effective Hamiltonian and is broken by the pairing terms. In such a case, we numerically find a nontrivial systematic system-size dependence of the Thouless time, in contrast to a related recent study for kicked fermionic chains.

Boundary chaos.

Scrambling in many-body quantum systems causes initially local observables to spread uniformly over the whole available Hilbert space under unitary dynamics, which in lattice systems causes exponential suppression of dynamical correlation functions with system size. Here, we present a perturbed free quantum circuit model, in which ergodicity is induced by an impurity interaction placed on the system's boundary, that allows for demonstrating the underlying mechanism. This is achieved by mapping dynamical correlation functions of local operators acting at the boundary to a partition function with complex weights defined on a two-dimensional lattice with a helical topology. We evaluate this partition function in terms of transfer matrices, which allow for numerically treating system sizes far beyond what is accessible by exact diagonalization and whose spectral properties determine the asymptotic scaling of correlations. Combining analytical arguments with numerical results, we show that for impurities which remain unitary under partial transpose correlations are exponentially suppressed with system size in a particular scaling limit. In contrast, for generic impurities or generic locations of the local operators, correlations show persistent revivals with a period given by the system size.

Exact Anomalous Current Fluctuations in a Deterministic Interacting Model.

We analytically compute the full counting statistics of charge transfer in a classical automaton of interacting charged particles. Deriving a closed-form expression for the moment generating function with respect to a stationary equilibrium state, we employ asymptotic analysis to infer the structure of charge current fluctuations for a continuous range of timescales. The solution exhibits several unorthodox features. Most prominently, on the timescale of typical fluctuations the probability distribution of the integrated charge current in a stationary ensemble without bias is distinctly non-Gaussian despite diffusive behavior of dynamical charge susceptibility. While inducing a charge imbalance is enough to recover Gaussian fluctuations, we find that higher cumulants grow indefinitely in time with different exponents, implying singular scaled cumulants. We associate this phenomenon with the lack of a regularity condition on moment generating functions and the onset of a dynamical critical point. In effect, the scaled cumulant generating function does not, irrespectively of charge bias, represent a faithful generating function of the scaled cumulants, yet the associated Legendre dual yields the correct large-deviation rate function. Our findings hint at novel types of dynamical universality classes in deterministic many-body systems.

Exact solution of the "Rule 150" reversible cellular automaton.

We study the dynamics and statistics of the Rule 150 reversible cellular automaton (RCA). This is a one-dimensional lattice system of binary variables with synchronous (Floquet) dynamics that corresponds to a bulk deterministic and reversible discretized version of the kinetically constrained "exclusive one-spin facilitated" (XOR) Fredrickson-Andersen (FA) model, where the local dynamics is restricted: A site flips if and only if its adjacent sites are in different states from each other. Similar to other RCA that have been recently studied, such as Rule 54 and Rule 201, the Rule 150 RCA is integrable, however, in contrast is noninteracting: The emergent quasiparticles, which are identified by the domain walls, behave as free fermions. This property allows us to solve the model by means of matrix product ansatz. In particular, we find the exact equilibrium and nonequilibrium stationary states for systems with closed (periodic) and open (stochastic) boundaries, respectively, resolve the full spectrum of the time evolution operator and, therefore, gain access to the relaxation dynamics, and obtain the exact large deviation statistics of dynamical observables in the long-time limit.

Absence of Normal Fluctuations in an Integrable Magnet.

We investigate dynamical fluctuations of transferred magnetization in the one-dimensional lattice Landau-Lifshitz magnet with uniaxial anisotropy, representing an emblematic model of interacting spins. We demonstrate that the structure of fluctuations in thermal equilibrium depends radically on the characteristic dynamical scale. In the ballistic regime, typical fluctuations are found to follow a normal distribution and scaled cumulants are finite. In stark contrast, on the diffusive and superdiffusive timescales, relevant, respectively, for the easy-axis and isotropic magnet at vanishing total magnetization, typical fluctuations are no longer Gaussian and, remarkably, scaled cumulants are divergent. The observed anomalous features disappear upon breaking integrability, suggesting that the absence of normal fluctuations is intimately tied to the presence of soliton modes. In a nonequilibrium setting of the isotropic magnet with weakly polarized step-profile initial state we find a slow drift of dynamical exponent from the superdiffusive towards the diffusive value.

Spectral Statistics of Non-Hermitian Matrices and Dissipative Quantum Chaos.

We propose a measure, which we call the dissipative spectral form factor (DSFF), to characterize the spectral statistics of non-Hermitian (and nonunitary) matrices. We show that DSFF successfully diagnoses dissipative quantum chaos and reveals correlations between real and imaginary parts of the complex eigenvalues up to arbitrary energy scale (and timescale). Specifically, we provide the exact solution of DSFF for the complex Ginibre ensemble (GinUE) and for a Poissonian random spectrum (Poisson) as minimal models of dissipative quantum chaotic and integrable systems, respectively. For dissipative quantum chaotic systems, we show that the DSFF exhibits an exact rotational symmetry in its complex time argument τ. Analogous to the spectral form factor (SFF) behavior for Gaussian unitary ensemble, the DSFF for GinUE shows a "dip-ramp-plateau" behavior in |τ|: the DSFF initially decreases, increases at intermediate timescales, and saturates after a generalized Heisenberg time, which scales as the inverse mean level spacing. Remarkably, for large matrix size, the "ramp" of the DSFF for GinUE increases quadratically in |τ|, in contrast to the linear ramp in the SFF for Hermitian ensembles. For dissipative quantum integrable systems, we show that the DSFF takes a constant value, except for a region in complex time whose size and behavior depend on the eigenvalue density. Numerically, we verify the above claims and additionally show that the DSFF for real and quaternion real Ginibre ensembles coincides with the GinUE behavior, except for a region in the complex time plane of measure zero in the limit of large matrix size. As a physical example, we consider the quantum kicked top model with dissipation and show that it falls under the Ginibre universality class and Poisson as the "kick" is switched on or off. Lastly, we study spectral statistics of ensembles of random classical stochastic matrices or Markov chains and show that these models again fall under the Ginibre universality class.

Many-body quantum chaos and dual-unitarity round-a-face.

We propose a new type of locally interacting quantum circuits-quantum cellular automata-that are generated by unitary interactions round-a-face (IRF). Specifically, we discuss a set (or manifold) of dual-unitary IRFs with local Hilbert space dimension d [DUIRF (d)], which generate unitary evolutions both in space and time directions of an extended 1+1 dimensional lattice. We show how arbitrary dynamical correlation functions of local observables can be evaluated in terms of finite-dimensional completely positive trace preserving unital maps in complete analogy to recently studied circuits made of dual-unitary brick gates (DUBGs). The simplest non-vanishing local correlation functions in dual-unitary IRF circuits are shown to involve observables non-trivially supported on two neighboring sites. We completely characterize the ten-dimensional manifold of DUIRF (2) for qubits ( d=2) and provide, for d=3,4,…,7, empirical estimates of its dimensionality based on numerically determined dimensions of tangent spaces at an ensemble of random instances of dual-unitary IRF gates. In parallel, we apply the same algorithm to determine dimDUBG(d) and show that they are of similar order though systematically larger than dimDUIRF(d) for d=2,3,…,7. It is remarkable that both sets have a rather complex topology for d≥3 in the sense that the dimension of the tangent space varies among different randomly generated points of the set. Finally, we provide additional data on dimensionality of the chiral extension of DUBG circuits with distinct local Hilbert spaces of dimensions d≠d' residing at even/odd lattice sites.

Eigenstate thermalization in dual-unitary quantum circuits: Asymptotics of spectral functions.

The eigenstate thermalization hypothesis provides to date the most successful description of thermalization in isolated quantum systems by conjecturing statistical properties of matrix elements of typical operators in the (quasi)energy eigenbasis. Here we study the distribution of matrix elements for a class of operators in dual-unitary quantum circuits in dependence of the frequency associated with the corresponding eigenstates. We provide an exact asymptotic expression for the spectral function, i.e., the second moment of this frequency resolved distribution. The latter is obtained from the decay of dynamical correlations between local operators which can be computed exactly from the elementary building blocks of the dual-unitary circuits. Comparing the asymptotic expression with results obtained by exact diagonalization we find excellent agreement. Small fluctuations at finite system size are explicitly related to dynamical correlations at intermediate times and the deviations from their asymptotical dynamics. Moreover, we confirm the expected Gaussian distribution of the matrix elements by computing higher moments numerically.

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.

Signatures of Chaos in Nonintegrable Models of Quantum Field Theories.

We study signatures of quantum chaos in (1+1)D quantum field theory (QFT) models. Our analysis is based on the method of Hamiltonian truncation, a numerical approach for the construction of low-energy spectra and eigenstates of QFTs that can be considered as perturbations of exactly solvable models. We focus on the double sine-Gordon, also studying the massive sine-Gordon and ϕ^{4} model, all of which are nonintegrable and can be studied by this method with sufficiently high precision from small to intermediate perturbation strength. We analyze the statistics of level spacings and of eigenvector components, which are expected to follow random matrix theory predictions. While level spacing statistics are close to the Gaussian orthogonal ensemble (GOE) as expected, on the contrary, the eigenvector components follow a distribution markedly different from the expected Gaussian. Unlike in the typical quantum chaos scenario, the transition of level spacing statistics to chaotic behavior takes place already in the perturbative regime. Moreover, the distribution of eigenvector components does not appear to change or approach Gaussian behavior, even for relatively large perturbations. Our results suggest that these features are independent of the choice of model and basis.

Inhomogeneous matrix product ansatz and exact steady states of boundary-driven spin chains at large dissipation.

We find novel site-dependent Lax operators in terms of which we demonstrate exact solvability of a dissipatively driven XYZ spin-1/2 chain in the Zeno limit of strong dissipation, with jump operators polarizing the boundary spins in arbitrary directions. We write the corresponding nonequilibrium steady state using an inhomogeneous matrix product ansatz, where the constituent matrices satisfy a simple set of linear recurrence relations. Although these matrices can be embedded into an infinite-dimensional auxiliary space, we have verified that they cannot be simultaneously put into a tridiagonal form, not even in the case of axially symmetric (XXZ) bulk interactions and general nonlongitudinal boundary dissipation. We expect our results to have further fundamental applications for the construction of nonlocal integrals of motion for the open XYZ model with arbitrary boundary fields, or the eight-vertex model.

Exact Nonequilibrium Steady State of Open XXZ/XYZ Spin-1/2 Chain with Dirichlet Boundary Conditions.

We investigate a dissipatively driven XYZ spin-1/2 chain in the Zeno limit of strong dissipation, described by the Lindblad master equation. The nonequilibrium steady state is expressed in terms of a matrix product ansatz using novel site-dependent Lax operators. The components of Lax operators satisfy a simple set of linear recurrence equations that generalize the defining algebraic relations of the quantum group U_{q}(sl_{2}). We reveal connection between the nonequilibrium steady state of the nonunitary dynamics and the respective integrable model with edge magnetic fields, described by coherent unitary dynamics.

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.

Quantum chaos challenges many-body localization.

Characterizing states of matter through the lens of their ergodic properties is a fascinating new direction of research. In the quantum realm, the many-body localization (MBL) was proposed to be the paradigmatic ergodicity breaking phenomenon, which extends the concept of Anderson localization to interacting systems. At the same time, random matrix theory has established a powerful framework for characterizing the onset of quantum chaos and ergodicity (or the absence thereof) in quantum many-body systems. Here we numerically study the spectral statistics of disordered interacting spin chains, which represent prototype models expected to exhibit MBL. We study the ergodicity indicator g=log_{10}(t_{H}/t_{Th}), which is defined through the ratio of two characteristic many-body time scales, the Thouless time t_{Th} and the Heisenberg time t_{H}, and hence resembles the logarithm of the dimensionless conductance introduced in the context of Anderson localization. We argue that the ergodicity breaking transition in interacting spin chains occurs when both time scales are of the same order, t_{Th}≈t_{H}, and g becomes a system-size independent constant. Hence, the ergodicity breaking transition in many-body systems carries certain analogies with the Anderson localization transition. Intriguingly, using a Berezinskii-Kosterlitz-Thouless correlation length we observe a scaling solution of g across the transition, which allows for detection of the crossing point in finite systems. We discuss the observation that scaled results in finite systems by increasing the system size exhibit a flow towards the quantum chaotic regime.