---
title: Operator Spreading in Quantum Systems
url: https://www.emergentmind.com/topics/operator-spreading
type: topic
---

# Operator Spreading in Quantum Systems

Operator spreading is a central concept in the study of quantum information dynamics, quantum chaos, and thermalization in many-body quantum systems. It refers to the growth in spatial support (“size”) of an initially local operator under Heisenberg evolution, and provides a rigorous framework for understanding information scrambling, emergent hydrodynamics, and the structure of quantum noise and decoherence. Operator spreading is deeply connected to key observables such as out-of-time-ordered correlators (OTOCs), operator size distributions, and information-theoretic capacities, and can be quantitatively analyzed using a variety of techniques in both closed and open quantum systems.

## 1. Formal Definitions and Operator Size

In the Heisenberg picture, a local operator $O$ evolves as $O(t) = e^{iHt} O e^{-iHt}$, where $H$ is the system Hamiltonian. Initially localized, $O$ typically spreads in operator space, acquiring nonzero support over increasingly nonlocal basis elements — for spin systems, sums of Pauli strings with growing domain size. This process is formalized by:

- **Expansion in operator basis**: $O(t) = \sum_\alpha c_{\alpha}(t) W_{\alpha}$, where $\{W_\alpha\}$ are multi-site Pauli strings or their generalizations.
- **Operator size**: The “size” $|W_\alpha|$ of a basis element is defined as the number of sites on which it acts nontrivially.
- **Size distribution**: $P(s, t) = \sum_{|W_\alpha| = s} |c_\alpha(t)|^2$.
- **Mean and variance**: $\langle s(t) \rangle = \sum_s s\,P(s, t)$ and $\mathrm{Var}[s(t)] = \sum_s s^2 P(s, t) - \langle s(t) \rangle^2$.

Physically, operator spreading quantifies how an initially local perturbation becomes delocalized due to quantum interactions, encoding the microscopic structure of information flow and chaos [2102.00010].

## 2. Signatures, Diagnostics, and Hydrodynamics

**Out-of-time-ordered correlators (OTOCs):** The squared commutator
\[
C_{ij}(t) = \langle [O_i(t), O_j(0)]^\dagger [O_i(t), O_j(0)] \rangle
\]
measures how non-commutativity spreads, signaling information scrambling and operator growth [1705.08975,2106.16170]. For short-range Hamiltonians, OTOCs exhibit a ballistic “light cone” characterized by a butterfly velocity $v_B$, with the front broadening diffusively ($\sim\sqrt{t}$ in 1D chaotic or generic integrable chains) [1808.04889,2103.13414,1809.02126,1710.09835].

| System               | Front Velocity ($v_B$) | Broadening          | Reference         |
|----------------------|----------------------:|---------------------|-------------------|
| Chaotic spin chains  | $>0$                  | $w \sim t^{1/2}$    | [1705.08975,2103.13414]   |
| Integrable models    | $>0$                  | $w \sim t^{1/2}$ (diffusive); anomalous scaling $H(t)\sim t^{-3/4}$ | [1809.02126,2103.13414]  |
| Free fermions        | $1$ (maximal)         | $w \sim t^{1/3}$    | [2102.09846]      |
| QHCG                 | $1$                   | Front freezes       | [2201.00395]      |

In higher dimensions, the operator front becomes a stochastic surface, with broadening governed by Kardar-Parisi-Zhang (KPZ) scaling exponents: $w\sim t^{1/3}$ in $2+1$D, $t^{0.24}$ in $3+1$D [1705.08975].

**Hydrodynamic Descriptions:** Operator right-weight (distribution of the rightmost nontrivial Pauli in the expansion) and generalized hydrodynamics allow mapping operator spreading to biased diffusion equations and, in integrable models, to ballistic quasiparticle propagation with density-dependent diffusion [2111.08727,1809.02126].

**Operator spreading with conservation laws:** Unitary dynamics with locally conserved charges yields two-stage hydrodynamics: a diffusive conserved component and a ballistic non-conserved sector, with operator weight conversion set by local diffusion current [1710.09835].

## 3. Random Circuits and Stochastic Fronts

Random unitary circuits provide an analytically tractable framework for the study of operator spreading and scrambling. In 1D Haar-random circuits, operator spreading is modeled as a biased diffusion of string endpoints:
\[
\partial_t \rho(x, t) = v_B\,\partial_x \rho + D\,\partial_x^2 \rho
\]
such that OTOCs are expressed as sums or integrals over the resulting densities. The butterfly velocity $v_B = (q^2-1)/(q^2+1)$ and front broadening $D=2q^2/(q^2+1)^2$ depend on the local dimension $q$ [1705.08975,2501.04091]. For general unitary-invariant gate ensembles, the drift-diffusion description persists, with universal features governed by low moments of the gate distribution, a finite “binary-ization” time for Pauli-string weights, and a domain-wall width for convergence to random-matrix-like statistics [2501.04091].

For long-range interactions decaying as $|i-j|^{-2\alpha}$, the butterfly light cone scales as $t$ (ballistic) for $\alpha>1$, stretched-exponentially for $1/2<\alpha<1$, and algebraically for $\alpha<1/2$, as captured by a noisy long-range FKPP equation dual to discrete population dynamics [2505.06353].

## 4. Experimental and Operational Perspectives

**Quantum simulation and tomography:** Experiments with trapped ions, Rydberg atom arrays, and superconducting circuits have probed operator spreading and OTOCs directly. Floquet engineering enables controlled measurements of light-cone-like propagation and tunable scrambling [2108.01276,2106.16170].

**Weak-measurement quantum tomography protocols** reconstruct the unknown quantum state by monitoring a Heisenberg-evolved observable under many-body dynamics. The rate of fidelity growth, Fisher information, Shannon entropy, and rank of the measurement record are directly linked to operator spreading; they scale rapidly for strongly chaotic systems and serve as robust chaos indicators, outperforming Krylov complexity which may lack monotonicity or correlational significance in integrable or finite-size systems [2308.08513,2404.09464].

**Information-theoretic equivalence:** Operator spreading is quantitatively connected to the classical information propagation capacity (Holevo bound) of a quantum channel. The distinguishability between channel outputs is upper- and lower-bounded by commutator norms or trace distances, establishing a tight equivalence between operator growth and information transmission, independent of the system size [2505.07955].

## 5. Open Quantum Systems and Dissipation

In open systems, operator spreading interacts with decoherence, local dissipation, and measurement. The Lindblad master equation governs Heisenberg operator evolution, with the Liouvillian spectrum setting late-time decay. Crucially, bulk dissipation accelerates the decay of autocorrelations: as the operator spreads, more local dissipators act on its support, leading to an initial acceleration regime, a plateau where decay rate is proportional to system size, and eventual crossover to asymptotic decay set by the Liouvillian gap [2309.03485,2208.12272]. The universal structure of decay holds for both Ising and Bose-Hubbard chains and is dominated at intermediate times by the ballistic operator spreading under unitary dynamics [2309.03485].

## 6. Integrability, Quantum Maps, and Special Cases

Integrable systems display ballistic operator fronts, but the mechanism differs from chaotic systems:
- In interacting integrable models, fronts are carried by the fastest quasiparticles and broaden diffusively due to equilibrium density fluctuations; the diffusive correction can be computed from generalized hydrodynamics [1809.02126,2103.13414].
- Free-fermion circuits or quantum hardcore gases (QHCG) exhibit diffusive or even frozen operator fronts; QHCGs, for certain initial conditions, demonstrate a front that eventually ceases to broaden, in contrast to the persistent diffusive front in chaotic or generic integrable systems [2201.00395,2102.09846].
- Quantum chaotic maps, such as the perturbed Arnold cat map, exhibit initially classical (Koopman) spreading until the Ehrenfest time, followed by quantum scrambling that saturates at the operator space dimension, diagnosable by the Shannon entropy of operator expansion coefficients and spectral form factor dynamics [1808.04889].

## 7. Distinction from Recoverability and Metrological Response

The arrival of operator support as diagnosed by OTOCs or right-weight does not guarantee that quantum information, such as parameter sensitivity, remains locally accessible. An operational hierarchy involving quantum Fisher information (QFI) reveals that while OTOC establishes a strict light cone for support arrival, local recoverability (e.g., single-site or block QFI) may lag behind, especially when gauge or $U(1)$ symmetries are broken, as parameter sensitivity disperses into nonlocal correlations. This separation is quantified by an ordering $F_j(t) \leq F_\mathrm{dec}(w, t) \leq F_Q(\rho_{\mathcal{N}_w}) \leq F_\mathrm{glob}$, indicating a gap between operator support and metrological accessibility [2605.02774].

---

**References**
- [1705.08975], [1710.09835], [1808.04889], [1809.02126], [1901.05793], [2102.00010], [2102.09846], [2103.13414], [2106.16170], [2108.01276], [2111.08727], [2201.00395], [2208.12272], [2308.08513], [2309.03485], [2404.09464], [2501.04091], [2505.06353], [2505.07955], [2605.02774]

Source: https://www.emergentmind.com/topics/operator-spreading