We systematically investigate whether classical hydrodynamic field theories can predict the long-time dynamics of many-particle quantum systems. We study both numerically and analytically the time evolution of a chain of spins (or qubits) subjected to stroboscopic dynamics. The time evolution is implemented by a sequence of local and nearest-neighbor gates that conserve the total magnetization. The long-time dynamics of such a system is believed to be describable by a hydrodynamic field theory, which, importantly, includes the effect of noise. Based on a field theoretical analysis and symmetry arguments, we map each operator in the spin model to the corresponding fields in hydrodynamics. This allows us to predict which expectation values decay exponentially and which decay with a hydrodynamic long-time tail. We illustrate these findings by studying the time evolution of all 255 Hermitian operators that can be defined on four neighboring sites. All operators not protected by hydrodynamics decay exponentially, while the others show a slow hydrodynamic decay. While most hydrodynamic power laws seem to follow the analytical predictions, we also discuss cases where there is an apparent discrepancy between analytics and the finite-size numerical data.
We systematically investigate whether classical hydrodynamic field theories can predict the long-time dynamics of many-particle quantum systems. We study both numerically and analytically the time evolution of a chain of spins (or qubits) subjected to stroboscopic dynamics. The time evolution is implemented by a sequence of local and nearest-neighbor gates that conserve the total magnetization. The long-time dynamics of such a system is believed to be describable by a hydrodynamic field theory, which, importantly, includes the effect of noise. Based on a field theoretical analysis and symmetry arguments, we map each operator in the spin model to the corresponding fields in hydrodynamics. This allows us to predict which expectation values decay exponentially and which decay with a hydrodynamic long-time tail. We illustrate these findings by studying the time evolution of all 255 Hermitian operators that can be defined on four neighboring sites. All operators not protected by hydrodynamics decay exponentially, while the others show a slow hydrodynamic decay. While most hydrodynamic power laws seem to follow the analytical predictions, we also discuss cases where there is an apparent discrepancy between analytics and the finite-size numerical data.
In thermalizing many-body quantum systems without conservation laws, such as ergodic Floquet and random unitary circuits, local expectation values are predicted to decay to their equilibrium values exponentially quickly. In this work, we derive a relationship between said exponential decay rate g ¯ and the operator spreading properties of a local unitary evolution. A hydrodynamical picture for operator spreading allows us to argue that, for random unitary circuits, g ¯ is encoded by the leading eigenvalue of a dynamical map obtained by enriching unitary dynamics with dissipation, in the limit of weak dissipation. We argue that the size of the eigenvalue does not depend on the details of this weak dissipation (given mild assumptions on properties of the ergodic dynamics), so long as it only suppresses large operators significantly. Our calculations are based on analytical results for random unitary circuits, but we argue that similar results hold for ergodic Floquet systems. These conjectures are in accordance with existing results that numerically obtain quantum many-body analogs of classical Ruelle-Pollicott resonances [Prosen, ; Mori, ] by studying unitary evolutions subject to weak dissipation.
The interplay between symmetry and thermalization governs the late-time dynamics of local quantum and classical many-body systems at nonzero temperature. Recently, two parallel frontiers have emerged: the search for robust anomalous hydrodynamics -- such as superdiffusion -- in generic, non-integrable models, and the formal effort to generalize the fundamental concept of global symmetry. In this paper, we bridge these frontiers by demonstrating that quantum group symmetry provides a novel mechanism for anomalous hydrodynamics in chaotic systems. We study the dynamics of local operators carrying U(1) charge in non-integrable lattice models that also have quantum group symmetry. One example is transverse spin in the XXZ model with integrability breaking deformations. While such excitations are expected to decay very quickly at high temperature because their charge forbids overlap with conventional hydrodynamic densities, we find that protection by the quantum group symmetry makes these modes long-lived, despite the absence of local quantum group charge density or current. Furthermore, the dynamics is superdiffusive across Hamiltonian, Floquet, and classical realizations, and exhibits unusual finite size effects at very late times. We also revisit transverse spin dynamics in the integrable XXZ model.
In ergodic quantum spin chains, locally conserved quantities such as energy or particle number generically evolve according to hydrodynamic equations as they relax to equilibrium. We investigate the complexity of simulating hydrodynamics at infinite temperature with multiple methods: time evolving block decimation (TEBD), TEBD with density matrix truncation (DMT), the recursion method with a universal operator growth hypothesis, and operator-size truncated (OST) dynamics. Density matrix truncation and the OST dynamics give consistent dynamical correlations to t=60/J and diffusion coefficients agreeing within 1. TEBD only converges for $t\ensuremath{\lesssim}20$, but still produces diffusion coefficients accurate within 1%. The universal operator growth hypothesis fails to converge and only matches other methods on short times. We see no evidence of long-time tails in either DMT or OST calculations of the current-current correlator, although we cannot rule out that they appear at longer times. We do observe power-law corrections to the energy density correlator. At finite wavelength, we observe a crossover from purely diffusive, overdamped decay of the energy density to underdamped oscillatory behavior similar to that of cold-atom experiments. We dub this behavior ``hot band second sound'' and offer a microscopically motivated toy model.
Random quantum circuits are paradigmatic models of minimally structured and analytically tractable chaotic dynamics. We study a family of Floquet unitary circuits with Haar random U(1) charge conserving dynamics; the minimal such model has nearest-neighbor gates acting on spin-1/2 qubits, and a single layer of even/odd gates repeated periodically in time. We find that this minimal model is not robustly thermalizing at numerically accessible system sizes, and displays slow subdiffusive dynamics for long times. We map out the thermalization dynamics in a broader parameter space of charge conserving circuits, and understand the origin of the slow dynamics in terms of proximate localized and integrable regimes in parameter space. In contrast, we find that small extensions to the minimal model are sufficient to achieve robust thermalization; these include (i) increasing the interaction range to three-site gates, (ii) increasing the local Hilbert space dimension by appending an additional unconstrained qubit to the conserved charge on each site, or (iii) using a larger Floquet period comprised of two independent layers of gates. Our results should inform future studies of charge conserving circuits, which are relevant for a wide range of topical theoretical questions.
The scrambling of quantum information in closed many-body systems, as measured by out-of-timeordered correlation functions (OTOCs), has received considerable attention lately.Recently, a hydrodynamical description of OTOCs has emerged from considering random local circuits.Numerical work suggests that aspects of this description are universal to ergodic many-body systems, even without randomness; a conjectured explanation for this is that while the random circuits have noise built into them, deterministic quantum systems, much like classically chaotic ones, "generate their own noise" and look effectively random on sufficient length scales and timescales.In this paper, we extend this approach to systems with locally conserved quantities (e.g., energy).We do this by considering local random unitary circuits with a conserved U(1) charge and argue, with numerical and analytical evidence, that the presence of a conservation law slows relaxation in both time-ordered and out-of-time-ordered correlation functions; both can have a diffusively relaxing component or "hydrodynamic tail" at late times.We verify the presence of such tails also in a deterministic, periodically driven system.We show that for OTOCs, the combination of diffusive and ballistic components leads to a wave front with a specific asymmetric shape, decaying as a power law behind the front.These results also explain existing numerical investigations in non-noisy ergodic systems with energy conservation.Moreover, we consider OTOCs in Gibbs states, parametrized by a chemical potential μ, and apply perturbative arguments to show that for μ ≫ 1 the ballistic front of information spreading can only develop at times exponentially large in μ-with the information traveling diffusively at earlier times.We also develop a new formalism for describing OTOCs and operator spreading, which allows us to interpret the saturation of OTOCs as a form of thermalization on the Hilbert space of operators.
Abstract By the hydrodynamic linear response theory, dynamical correlation functions decay as power laws along certain velocities, determined by the flux Jacobian. Such correlations are obtained by hydrodynamic projections, and physically, they are due to propagating ‘sound waves’ or generalisation thereof, transporting conserved quantities between the observables. However, some observables do not emit sound waves, such as order parameters associated to symmetry breaking. In these cases correlation functions decay exponentially everywhere, a behaviour not captured by the hydrodynamic linear response theory. Focussing on spin–spin correlation functions in the XX quantum chain, we first review how hydrodynamic linear response works, emphasising that the necessary fluid cell averaging washes out oscillatory effects. We then show how, beyond linear response, Euler hydrodynamics can still predict the exponential decay of correlation functions of order parameters. This is done by accounting for the large-scale fluctuations of domain walls, via the recently developed ballistic fluctuation theory. We use the framework of generalised hydrodynamics, which is particularly simple in this model due to its free fermion description. In particular, this reproduces, by elementary calculations, the exponential decay in the celebrated formulae by Its et al (1993) and by Jie (1998), which were originally obtained by intricate Fredholm determinant analysis; and gives a new formula in a parameter domain where no result was obtained before. We confirm the results by numerical simulations.