Abstract In this work, we derive an exact hydrodynamical description for the coupled, charge and operator dynamics, in a quantum many-body system with U(1) symmetry. Using an emergent symmetry in the complex Brownian SYK model with charge conservation, we map the operator dynamics in the model to the imaginary-time dynamics of an SU(4) spin-chain. We utilize the emergent SU(4) description to demonstrate that the U(1) symmetry causes quantum-coherence to persist even after disorder-averaging, in sharp contrast to models without symmetries. In line with this property, we write down a ‘restricted’ Fokker-Planck equation for the out-of-time ordered correlator (OTOC) in the large-N limit, which permits a classical probability description strictly in the incoherent sector of the global operator-space. We then exploit this feature to describe the OTOC in terms of a Fisher-Kolmogorov-Petrovsky-Piskun (FKPP)-equation which couples the operator with the charge and is valid at all time-scales and for arbitrary charge-density profiles. The coupled equations obtained belong to a class of models also used to describe the population dynamics of bacteria embedded in a diffusive media. We simulate them to explore operator-dynamics in a background of non-uniform charge configuration, which reveals that the charge transport can strongly affect dynamics of operators, including those that have no overlap with the charge.
Abstract In this work, we derive an exact hydrodynamical description for the coupled, charge and operator dynamics, in a quantum many-body system with U(1) symmetry. Using an emergent symmetry in the complex Brownian SYK model with charge conservation, we map the operator dynamics in the model to the imaginary-time dynamics of an SU(4) spin-chain. We utilize the emergent SU(4) description to demonstrate that the U(1) symmetry causes quantum-coherence to persist even after disorder-averaging, in sharp contrast to models without symmetries. In line with this property, we write down a ‘restricted’ Fokker-Planck equation for the out-of-time ordered correlator (OTOC) in the large-N limit, which permits a classical probability description strictly in the incoherent sector of the global operator-space. We then exploit this feature to describe the OTOC in terms of a Fisher-Kolmogorov-Petrovsky-Piskun (FKPP)-equation which couples the operator with the charge and is valid at all time-scales and for arbitrary charge-density profiles. The coupled equations obtained belong to a class of models also used to describe the population dynamics of bacteria embedded in a diffusive media. We simulate them to explore operator-dynamics in a background of non-uniform charge configuration, which reveals that the charge transport can strongly affect dynamics of operators, including those that have no overlap with the charge.
Abstract In this work, we introduce a symmetry-based approach to study the scrambling and operator dynamics of Brownian SYK models at large finite N and in the infinite N limit. We compute the out-of-time-ordered correlator (OTOC) in the Majorana model without charge conservation and the complex model with charge conservation, and demonstrate that in both models taking the random average of the couplings gives rise to emergent symmetry structures. The random averaging exactly maps the operator dynamics of the Majorana model and the complex model to the imaginary time dynamics of an SU(2) spin and an SU(4) spin respectively, which become solvable in the large N limit. Furthermore, the symmetry structure drastically reduces the size of the Hilbert space required to calculate the OTOC from exponential to linear in N, providing full access to the operator dynamics at all times for large finite N. In the case of the complex model with charge conservation, using this approach, we obtain the OTOC within each charge sector both numerically at finite N and analytically in the large N limit. We find that the time scale of the scrambling dynamics for all times and in each sector is characterized by the charge density. Furthermore, after proper rescaling, the OTOC corresponding to different finite charge densities collapses into a single curve at large finite N. In the large N limit, the rescaled OTOCs at finite density are described by the same hydrodynamic equation as in the Majorana case.
Operator scrambling is commonly diagnosed by the growth of out-of-time-ordered correlators (OTOCs), yet a general symmetry principle underlying their effective dynamics has remained elusive. For Brownian or short-time-correlated large-N Majorana systems, we develop a symmetry-based effective field theory for operator scrambling, organized by a strong-to-weak U(1) symmetry breaking in operator space. The key observation is that, in the noninteracting fermion limit, the four-fold Keldysh contour representation of an OTOC admits an emergent strong U(1) symmetry in a doubled Hilbert-space description, even when the original system has no ordinary conserved quantity. The associated slow mode is the phase of the strong-charge creation operator, whose conjugate density is identified with the local operator size. Generic interactions explicitly break the strong symmetry and generate a mass term at lowest order for the would-be Goldstone mode, thereby converting diffusive operator spreading into chaotic growth. We further show that higher-order symmetry breaking terms are tightly constrained by an emergent duality that combines time reversal with contour permutation. This duality fixes the effective action up to quadratic order in the response field, relates the multiplicative noise strength directly to the Lyapunov exponent, and makes the positivity of the Lyapunov exponent a consequence of convergence of the real-time path integral. The resulting OTOC dynamics is governed by a noisy FKPP equation, which captures within a unified framework the early-time exponential growth, ballistic propagation, nonlinear saturation, and stochastic front broadening of operator scrambling. We verify this construction in a Brownian SYK chain, where a direct saddle-point expansion reproduces the symmetry-based effective action. Our results reveal a symmetry origin of operator-size hydrodynamics and scrambling.
Thermalization of chaotic quantum many-body systems under unitary time evolution is related to the growth in complexity of initially simple Heisenberg operators.Operator growth is a manifestation of information scrambling and can be diagnosed by out-of-time-order correlators (OTOCs).However, the behavior of OTOCs of local operators in generic chaotic local Hamiltonians remains poorly understood, with some semiclassical and large-N models exhibiting exponential growth of OTOCs and a sharp chaos wave front and other random circuit models showing a diffusively broadened wave front.In this paper, we propose a unified physical picture for scrambling in chaotic local Hamiltonians.We construct a random time-dependent Hamiltonian model featuring a large-N limit where the OTOC obeys a Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) type equation and exhibits exponential growth and a sharp wave front.We show that quantum fluctuations manifest as noise (distinct from the randomness of the couplings in the underlying Hamiltonian) in the FKPP equation and that the noise-averaged OTOC exhibits a crossover to a diffusively broadened wave front.At small N, we demonstrate that operator growth dynamics, averaged over the random couplings, can be efficiently simulated for all time using matrix product state techniques.To show that time-dependent randomness is not essential to our conclusions, we push our previous matrix product operator methods to very large size and show that data for a timeindependent Hamiltonian model are also consistent with a diffusively broadened wave front.
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.
We show that operator dynamics in U(1) symmetric systems are constrained by two independent conserved charges and construct a non-Abelian operator size basis that respects both, enabling a symmetry-resolved characterization of operator growth. The non-Abelian operator size depends on the operator's nonlocal structure and is organized by an SU(2) algebra. Operators associated with large total angular momentum are relatively simple, while those with small angular momentum are more complex. Operator growth is thus characterized by a reduction in angular momentum and can be probed using out-of-time-ordered correlators. Using the charge-conserving Brownian Sachdev-Ye-Kitaev model, we derive an exact classical master equation that governs the size distribution, the distribution of an operator expanded in this basis, for arbitrary system sizes. The resulting dynamics reveal that the size distribution follows a chi-squared form, with the two conserved charges jointly determining the overall time scale and the shape of the distribution. In particular, single-particle operators retain a divergent peak at large angular momentum throughout the time evolution.
Abstract Quantum scrambling plays an important role in understanding thermalization in closed quantum systems. By this effect, quantum information spreads throughout the system and becomes hidden in the form of non-local correlations. Alternatively, it can be described in terms of the increase in complexity and spatial support of operators in the Heisenberg picture, a phenomenon known as operator growth. In this work, we study the disordered fully-connected Sachdev-Ye-Kitaev (SYK) model, and we demonstrate that scrambling is absent for disorder-averaged expectation values of observables. In detail, we adopt a formalism typical of open quantum systems to show that, on average and within charge-conserved sectors, operators evolve in a relatively simple way which is governed by their operator size. This feature only affects single-time correlation functions, and in particular it does not hold for out-of-time-order correlators, which are well-known to show scrambling behavior. Making use of these findings, we develop a cumulant expansion approach to approximate the evolution of equal-time observables. We employ this scheme to obtain analytic results that apply to arbitrary system size, and we benchmark its effectiveness by exact numerics. Our findings shed light on the structure of the dynamics of observables in the SYK model, and provide an approximate numerical description that overcomes the limitation to small systems of standard methods.