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.
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 We study operator dynamics in Brownian quantum many-body models with q -local interactions. The operator dynamics are characterized by the time-dependent size distribution, for which we derive an exact master equation in both the Brownian Majorana Sachdev–Ye–Kitaev model and the spin model for general q . This equation can be solved numerically for large systems. Additionally, we obtain the analytical size distribution in the large N limit for arbitrary initial conditions and q . The distributions for both models take the same form, related to the χ -squared distribution by a change of variable, and strongly depend on the initial condition. For small initial sizes, the operator dynamics are characterized by a broad distribution that narrows as the initial size increases. When the initial operator size is below q − 2 for the Majorana model or q − 1 for the spin model, the distribution diverges in the small size limit at all times. The mean size of all operators, which can be directly measured by the out-of-time ordered correlator, grows exponentially during the early time. In the late time regime, the mean size for a single Majorana or Pauli operator for all q decays exponentially as t e − t , much slower than all other operators, which decay as e − t . At finite N , the size distribution exhibits modulo-dependent branching within a symmetry sector for the q ⩾ 8 Majorana model and the q ⩾ 4 spin model. Our results reveal universal features of operator dynamics in q -local quantum many-body systems.
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.
We study quantum many-body systems with a global U(1) conservation law, focusing on a theory of N interacting fermions with charge conservation, or N interacting spins with one conserved component of total spin. We define an effective operator size at finite chemical potential through suitably regularized out-of-time-ordered correlation functions. The growth rate of this density-dependent operator size vanishes algebraically with charge density; hence we obtain new bounds on Lyapunov exponents and butterfly velocities in charged systems at a given density, which are parametrically stronger than any Lieb-Robinson bound. We argue that the density dependence of our bound on the Lyapunov exponent is saturated in the charged Sachdev-Ye-Kitaev model. We also study random automaton quantum circuits and Brownian Sachdev-Ye-Kitaev models, each of which exhibit a different density dependence for the Lyapunov exponent, and explain the discrepancy. We propose that our results are a cartoon for understanding Planckian-limited energy-conserving dynamics at finite temperature.
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.
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.
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.