Universal Krylov Complexity Growth Laws

• Universal Krylov Complexity Growth Laws is a framework quantifying operator spreading in many-body quantum systems by projecting Heisenberg-evolved operators onto a Krylov basis derived from the Lanczos algorithm.
• It reveals distinct temporal regimes—exponential, linear, and power-law growth—that reflect the system’s chaotic, integrable, or critical nature through the asymptotic behavior of Lanczos coefficients.
• The universal laws link mathematical structure with physical phenomena by connecting operator entropy, symmetry effects, and holographic dualities to the evolution of complexity.

Universal Krylov Complexity Growth Laws

Krylov complexity quantifies the operator growth in quantum many-body dynamics by projecting Heisenberg-evolved operators onto a specially constructed Krylov basis derived via the Lanczos algorithm. Its evolution—captured by the interplay of the Liouvillian superoperator, the underlying system Hamiltonian, and the associated sequence of Lanczos coefficients—displays a set of universal growth laws. These laws encode essential distinctions between integrable, chaotic, and open quantum systems, as well as structural features such as symmetry, temperature, and quantum criticality.

1. Foundational Framework: Krylov Basis and Operator Growth Hypothesis

The construction of Krylov complexity begins with the Krylov (Lanczos) basis: For an initial operator $O$, the operator evolution is governed by the Liouvillian $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$, generating the Krylov basis $\{|O_n)\}$ recursively via

$b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$

with $b_0=0$ and $|O_0)=O/\|O\|$. In this basis, $\mathcal{L}$ becomes tridiagonal, with off-diagonal elements $b_n$. Krylov complexity is then defined as

$K(t) = \sum_{n=0}^{K-1} n\,|\varphi_n(t)|^2,$

where $|\varphi_n(t)|^2$ is the probability amplitude of the operator wavepacket on the $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$0th Krylov site.

The Universal Operator-Growth Hypothesis (UOGH) posits that in generic many-body quantum systems, the off-diagonal Lanczos coefficients $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$1 grow at most linearly with $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$2 at large index: $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$3 with the "Lanczos slope" $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$4 achieving its maximal value in fully chaotic or nonintegrable (e.g., SYK-type) systems and growing sub-linearly in integrable or free systems (Bhattacharjee et al., 2022, Gamayun et al., 4 Apr 2025, Rabinovici et al., 8 Jul 2025).

2. Universal Growth Laws: Temporal Regimes

Krylov complexity exhibits a sequence of universal regimes connected with the dynamical structure of the $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$5 sequence. The archetypal behavior, valid for a wide class of isolated, closed quantum systems, consists of the following stages (Bhattacharjee et al., 2022, Rabinovici et al., 8 Jul 2025, Tang, 2023):

Regime	Lanczos Asymptotics	Krylov Complexity $\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$6	Typical Models
Early time	$\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$7	$\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$8	Chaotic/fast scrambling systems
Crossover	$\mathcal{L}[\,\cdot\,] = [H,\,\cdot\,]$9 enters plateau $\{|O_n)\}$0	$\{|O_n)\}$1	RMT, post-scrambling in local chaos
Saturation	$\{|O_n)\}$2 (finite chain edge)	$\{|O_n)\}$3	Finite Hilbert space, Heisenberg time

For polynomial $\{|O_n)\}$4 with $\{|O_n)\}$5 (integrable or critical scaling), one finds

$\{|O_n)\}$6

Sublinear $\{|O_n)\}$7 leads to power-law or bounded growth only.

Notably, exact analytical solutions for Krylov complexity exist for numerous solvable instances (e.g., Meixner–Pollaczek, continuous Hahn, alternating odd-even tail families), confirming exponential growth rates $\{|O_n)\}$8 for any $\{|O_n)\}$9 at large $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$0 (Gamayun et al., 4 Apr 2025).

3. Growth in Distinct Physical Contexts

3.1 Integrable and Saddle-Dominated Scrambling

Even in integrable systems, exponential Krylov complexity growth can emerge via saddle-dominated scrambling, where unstable saddles in phase space induce effective linear $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$1 scaling: $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$2 For the integrable Lipkin–Meshkov–Glick model, semiclassical analysis yields $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$3 (Bhattacharjee et al., 2022). This demonstrates that exponential operator growth is not a necessary signature of chaos.

3.2 Random Matrix Theory and Quantum Chaos Benchmarking

In the GUE (random matrix theory) setting, at infinite temperature, $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$4 rapidly approaches a plateau, and the maximal linear Krylov complexity growth $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$5 provides a benchmark: the linear growth speed in arbitrary local chaotic systems is conjectured not to exceed the RMT plateau value, imposing a sharp universal bound (consistent up to $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$6 fluctuations) (Tang, 2023). At low temperatures, $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$7 grows linearly with $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$8 up to a plateau, with slope saturating the chaos bound $b_{n+1}|O_{n+1}) = \mathcal{L}|O_n) - b_n|O_{n-1})$9; Krylov complexity then exhibits $b_0=0$0 until scrambling time $b_0=0$1, before entering linear growth.

3.3 Quantum Field Theory, UV Cutoff, and Holography

Across quantum field theoretic models, the universality of linear $b_0=0$2 (i.e., $b_0=0$3 in continuum theories) holds irrespective of chaos, tying the exponential growth rate to temperature scale: $b_0=0$4 This exponent is ubiquitous in both free and strongly coupled theories and is realized in holographic (AdS/CFT) models as a precise match—the growth of boundary Krylov complexity is dual to the radial momentum of an infalling probe in AdS, yielding $b_0=0$5 in the "Lyapunov" regime (Avdoshkin et al., 2022, Fan, 2024). However, in lattice realizations with finite ultraviolet cutoff, $b_0=0$6 saturates to a constant at large $b_0=0$7, causing a transition to linear $b_0=0$8 growth (or even bounded/oscillatory behavior in finite volume or integrable models).

4. Fundamental Bounds and Generalizations

4.1 Dispersion/Uncertainty Bound

For isolated quantum systems, Krylov complexity growth is universally bounded by a speed-limit derived from a Robertson-type uncertainty relation: $b_0=0$9 where $|O_0)=O/\|O\|$0 is the position variance on the Krylov chain; saturation occurs for coherent-state-like dynamics (e.g., linear $|O_0)=O/\|O\|$1 or SU(2)/Heisenberg-Weyl algebras) (Hörnedal et al., 2022).

In open (Markovian) quantum systems, the generalized non-Hermitian uncertainty relation constrains the normalized Krylov complexity $|O_0)=O/\|O\|$2 so that exponential early-time growth induced by linear $|O_0)=O/\|O\|$3 is always attenuated at late times by decoherence, causing complexity to decay (Bhattacharya et al., 2024).

4.2 Lyapunov Bound

Chaotic operator growth, as measured by the out-of-time-ordered correlator (OTOC) exponent $|O_0)=O/\|O\|$4, is universally bounded above by twice the Lanczos slope: $|O_0)=O/\|O\|$5 the final inequality saturating in holographic models and random matrix theory (Avdoshkin et al., 2022, Tang, 2023, Guo, 2022). Slow scramblers and integrable systems exhibit $|O_0)=O/\|O\|$6.

5. Universality and Extensions

5.1 Symmetry-Resolved Krylov Complexity and Equipartition

Systems with global symmetries (e.g., U(1) charge) admit a block structure in the Krylov basis. The symmetry-resolved Krylov complexity $|O_0)=O/\|O\|$7 for each sector matches (to $|O_0)=O/\|O\|$8) the total complexity averaged over sector weights. Exact equipartition ($|O_0)=O/\|O\|$9 for all $\mathcal{L}$0) occurs when the operator is energy-local or when sector Lanczos chains coincide; generically, sector mixing causes $\mathcal{L}$1, with deviations determined by inter-sector variance (Caputa et al., 2 Jul 2025).

5.2 Growth Laws at Quantum Phase Transitions

When a system is driven across a quantum critical point, Krylov complexity exhibits power-law growth controlled by Kibble–Zurek scaling. For a linear quench,

$\mathcal{L}$2

where $\mathcal{L}$3 and $\mathcal{L}$4 are critical exponents, $\mathcal{L}$5 spatial dimension, $\mathcal{L}$6 defect dimensionality; higher cumulants of the Krylov probability distribution share the same scaling, converging to a Gaussian distribution in the adiabatic limit (Grabarits et al., 15 Oct 2025).

5.3 Entropic Growth and Irreversibility

At long times, whenever the operator wavefunction irreversibly spreads over an infinite Krylov chain (i.e., in systems exhibiting true scrambling), the operator entropy $\mathcal{L}$7 universally satisfies

$\mathcal{L}$8

with $\mathcal{L}$9 for fully chaotic (linear $b_n$0) systems, $b_n$1 for integrable or sublinear growth (Fan, 2022). This relation is a hallmark of quantum irreversibility and fails for finite or periodic systems.

5.4 Generalized Krylov Complexities and Master Inequalities

Extensions to generalized Krylov complexities (moments of order $b_n$2, $b_n$3) reveal that (at long times) $b_n$4 holds universally for "slow scramblers," with model-dependent prefactors in fast scramblers (e.g., SYK-type). The relative variances satisfy the master inequality

$b_n$5

with the direction of stringency depending on the scrambler class (Fan, 2023).

6. Key Model Realizations and Exact Solutions

A range of exactly solvable models provide analytic access to the universal regime:

• Linear $b_n$6 (e.g., $b_n$7): archetype for maximal exponential growth, $b_n$8.
• Meixner–Pollaczek type ($b_n$9): retains universal exponential rate with tunable subleading corrections.
• Alternating families, continuous-Hahn, or double-slope $K(t) = \sum_{n=0}^{K-1} n\,|\varphi_n(t)|^2,$0: realize bounded/oscillatory $K(t) = \sum_{n=0}^{K-1} n\,|\varphi_n(t)|^2,$1 or two-slope behavior as seen in finite volume QFT or thermal AdS backgrounds (Gamayun et al., 4 Apr 2025, Avdoshkin et al., 2022).

These families confirm that subleading or staggered corrections in $K(t) = \sum_{n=0}^{K-1} n\,|\varphi_n(t)|^2,$2 affect only transients, not the late-time universal law.

7. Synthesis: Universality Classes and Physical Implications

The universal growth laws for Krylov complexity offer a precise characterization of operator spreading and information scrambling in diverse quantum many-body settings. The growth rate and crossover structure of $K(t) = \sum_{n=0}^{K-1} n\,|\varphi_n(t)|^2,$3 are direct fingerprints of the asymptotics of the Lanczos coefficients, serving as both a diagnostic for chaos/integrability and a kinematic "speed limit" for quantum dynamics. They uniformly govern systems as diverse as integrable chains, random matrix models, quantum field theories, and black hole duals, with explicit dynamical signatures in the structure of $K(t) = \sum_{n=0}^{K-1} n\,|\varphi_n(t)|^2,$4.

Universal Krylov complexity growth laws thus provide a robust organizing framework to quantify and classify complex operator dynamics in quantum systems, independent of model-specific control parameters, and establish sharp links between physical universality, mathematical structure, and information-theoretic limits (Bhattacharjee et al., 2022, Avdoshkin et al., 2022, Rabinovici et al., 8 Jul 2025, Gamayun et al., 4 Apr 2025, Caputa et al., 2 Jul 2025, Tang, 2023, Grabarits et al., 15 Oct 2025).