Krylov Complexity Meets Confinement

Basic Information

• Paper ID: 2511.03783
• Title: Krylov Complexity Meets Confinement
• Authors: Xuhao Jiang, Jad C. Halimeh, N. S. Srivatsa
• Institutions: Ludwig Maximilian University of Munich, Technical University of Munich, Max Planck Institute of Quantum Optics, Munich Center for Quantum Science and Technology, Kyung Hee University
• Classification: cond-mat.stat-mech, cond-mat.str-el, hep-lat, hep-th, quant-ph
• Publication Date: November 5, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2511.03783

Abstract

This paper investigates the relationship between confinement phenomena and Krylov complexity in quantum many-body systems. Using the Ising model with transverse and longitudinal fields, the authors demonstrate that Krylov state complexity serves as a sensitive and quantitative probe of confinement. Key findings include: quantum quenches within the ferromagnetic phase show significant suppression of Krylov complexity growth upon introducing the longitudinal field, reflecting slowed correlation dynamics due to confinement; quenches within the paramagnetic phase exhibit enhanced complexity with increasing longitudinal field, indicating absence of confinement; quenches across the critical point display orders-of-magnitude larger complexity with weak confinement trends. Importantly, in the confined region, the oscillation frequency of complexity corresponds to meson masses, with power spectrum peaks matching semiclassical predictions precisely.

Research Background and Motivation

Core Problems to Address

1. Quantum Information Characterization of Confinement: Confinement is a fundamental phenomenon in high-energy physics (e.g., quark confinement) and also appears in certain condensed matter systems. How to quantitatively characterize confinement dynamics using quantum information theory tools remains an open question.
2. Krylov Complexity as a Probe: While Krylov operator complexity has been proposed for studying confinement-deconfinement phase transitions in large-N gauge theories, how Krylov state complexity reflects real-time confinement features after quantum quenches has not been systematically studied.

Importance of the Problem

1. Interdisciplinary Significance: Connects high-energy physics, condensed matter physics, and quantum information theory
2. Experimental Relevance: Confinement phenomena in Ising models have been observed in quasi-one-dimensional compounds and cold atom systems, even verified on IBM quantum computers
3. Theoretical Value: Provides new perspective for understanding how quantum information propagation is constrained by confinement

Limitations of Existing Methods

1. Constraints of Traditional Probes: While entanglement entropy and correlation functions reflect confinement, they lack directness and quantitative precision
2. Operator Dependence: Power spectra of traditional observables may miss energy differences due to vanishing matrix elements
3. Lack of Unified Framework: Need a unified tool that both quantitatively characterizes information propagation and encodes bound state spectra

Research Motivation

The authors are inspired by the observation that confinement imposes strict constraints on information and correlation propagation in quantum many-body systems. Therefore, quantum complexity—a measure quantifying quantum information spreading—may be particularly sensitive to confinement phenomena.

Core Contributions

1. Establishing Krylov Complexity as a Confinement Probe: First systematic demonstration that Krylov state complexity serves as a sensitive and quantitative detection tool for confinement phenomena
2. Revealing Complexity Characteristics in Different Quench Regions:
   • Quenches within ferromagnetic phase: Longitudinal field causes significant complexity suppression
   • Quenches within paramagnetic phase: Complexity enhancement with increasing longitudinal field
   • Quenches across critical point: Orders-of-magnitude complexity increase, showing weak confinement trends
3. Krylov Spectroscopy Method: Demonstrates that power spectrum peaks of complexity precisely correspond to meson bound state masses, with high consistency to semiclassical predictions
4. Theoretical Advantage Analysis: Clarifies advantages of Krylov complexity over traditional operator expectation values—it is state-specific, operator-independent, and scans all accessible frequency components

Methods in Detail

Task Definition

Input:

• Initial state $|\Psi_0\rangle$ (fully polarized ferromagnetic or paramagnetic state)
• Post-quench Hamiltonian parameters $(h_x, h_z)$

Output:

• Time evolution of Krylov state complexity $C_k(t)$
• Power spectrum of complexity $S_k(\omega)$

Objective: Identify confinement signals by analyzing growth behavior and oscillation characteristics of $C_k(t)$

Model Architecture

1. Physical Model: Transverse-Longitudinal Field Ising Model

The Hamiltonian is: $\hat{H} = -J\sum_{j=1}^{N}\left[\hat{\sigma}^z_j\hat{\sigma}^z_{j+1} + h_x\hat{\sigma}^x_j + h_z\hat{\sigma}^z_j\right]$

where $J=1$ (energy unit), $h_x$ is the transverse field, and $h_z$ is the longitudinal field.

Physical Picture:

• $h_z=0$: Model is integrable with quantum phase transition at $h_x=1$
  • $h_x<1$: Ferromagnetic phase, domain walls propagate freely
  • $h_x>1$: Paramagnetic phase
• $h_z\neq 0$: Breaks integrability and $Z_2$ symmetry, producing linear confinement potential $\chi = 2Jh_z\bar{\sigma}$, where $\bar{\sigma}=(1-h_x^2)^{1/8}$ is spontaneous magnetization

2. Krylov Complexity Construction

Krylov Subspace Generation: Construct orthonormal basis $\{|K_n\rangle\}$ via Lanczos algorithm: $\mathcal{K}_n(\hat{H}, |\Psi_0\rangle) = \text{span}\{|\Psi_0\rangle, \hat{H}|\Psi_0\rangle, \ldots, \hat{H}^{n-1}|\Psi_0\rangle\}$

Tridiagonalization: $\hat{H}|K_n\rangle = \alpha_n|K_n\rangle + \beta_{n+1}|K_{n+1}\rangle + \beta_n|K_{n-1}\rangle$

where $\alpha_n$ (diagonal) and $\beta_n$ (off-diagonal) are Lanczos coefficients.

Time Evolution: Expand evolved state as $|\Psi(t)\rangle = \sum_n \psi_n(t)|K_n\rangle$, satisfying: $i\partial_t\psi_n(t) = \alpha_n\psi_n(t) + \beta_n\psi_{n-1}(t) + \beta_{n+1}\psi_{n+1}(t)$

This is equivalent to a tight-binding model on a semi-infinite chain.

Complexity Definition: $C_k(t) = \sum_n n|\psi_n(t)|^2$

Physical meaning: Quantifies wave function spreading in Krylov space, analogous to average particle position on a chain.

3. Power Spectrum Analysis

$S_k(\omega) = \left|\int_{-\infty}^{\infty} e^{i\omega t}C_k(t)dt\right|^2$

Computed via discrete Fourier transform (DFT) with time step $\Delta t=0.1$ and total time $t=100$.

Technical Innovations

1. State-Specific Probe

Unlike traditional operator expectation values $\langle\hat{O}(t)\rangle$, Krylov complexity depends only on initial state and Hamiltonian, avoiding information loss from vanishing matrix elements $\langle E_b|\hat{O}|E_a\rangle$.

2. Geometric Interpretation

Maps complexity dynamics to quantum particle motion on a semi-infinite chain, providing geometric picture of information spreading.

3. Full Frequency Scanning

Complexity dynamics scan all frequency components accessible from the initial state, providing complete energy level structure information.

4. Full Orthogonalization Lanczos (FOL)

Employs FOL algorithm for numerical stability, orthogonalizing each new vector against all previously generated basis vectors to avoid precision loss. Termination condition: $\beta_n < 10^{-12}$.

Experimental Setup

System Parameters

• System Size: Main results use $L=14$ lattice sites
• Boundary Conditions: Periodic boundary conditions
• Energy Unit: $J=1$

Quench Protocols

Protocol 1: Ferromagnetic Phase Quench

• Initial State: Fully polarized ferromagnetic state (all spins up)
• Quench Target: $h_x=0.25$ (within ferromagnetic phase)
• Longitudinal Field Range: $h_z \in [0, 0.4]$

Protocol 2: Paramagnetic Phase Quench

• Initial State: Paramagnetic state ($h_x=2$)
• Quench Target: $h_x=1.75$ (within paramagnetic phase)
• Longitudinal Field Range: $h_z \in [0, 0.4]$

Protocol 3: Critical Point Crossing Quench

• Initial State: Paramagnetic state ($h_x=2$)
• Quench Target: $h_x=0.25$ (ferromagnetic phase, crossing $h_x=1$ critical point)
• Longitudinal Field Range: $h_z \in [0, 0.4]$

Meson Mass Calculation

Semiclassical Method: Treat two domain walls as fermions in linear confinement potential, with effective two-body Hamiltonian: $H_{\text{two-body}} = \epsilon(p_1) + \epsilon(p_2) + \chi|x_2-x_1|$

where single-particle dispersion relation: $\epsilon(p) = 2J\sqrt{1+h_x^2-2h_x\cos p}$

Solve bound state energies via Bohr-Sommerfeld quantization condition: $2E_n(P)p_a - \int_{-p_a}^{p_a}dp\,\omega(p;P) = 2\pi\chi\left(n-\frac{1}{4}\right)$

where $\omega(p;P) = \epsilon(p+P/2) + \epsilon(p-P/2)$.

Specific Parameters:

• $h_x=0.25, h_z=0.2$: Two mesons $m_1=4.025J, m_2=4.702J$
• $h_x=0.25, h_z=0.1$: Four mesons $m_1=3.662, m_2=4.127, m_3=4.480, m_4=4.769$

Numerical Implementation Details

• Lanczos Algorithm: FOL variant with repeated Gram-Schmidt orthogonalization for numerical stability
• Time Evolution: Solve Schrödinger equation for tridiagonal Hamiltonian
• Power Spectrum: DFT with time window $t\in[0,100]$ and step $\Delta t=0.1$

Experimental Results

Main Findings

1. Ferromagnetic Phase Quench (Figure 1)

Key Discoveries:

• $h_z=0$: Complexity shows large-amplitude oscillations, amplitude increases with system size (inset a), consistent with light-cone propagation of correlations
• $h_z>0$:
  • Complexity growth is significantly suppressed
  • Suppression degree strengthens with $h_z$
  • Oscillation frequency increases but amplitude decreases
  • System size effects vanish (inset b), indicating blocked correlation propagation

Physical Interpretation: Longitudinal field gives mesons large effective mass; quench energy insufficient for meson motion, only static meson production occurs.

2. Paramagnetic Phase Quench (Figure 2)

Key Discoveries:

• $h_z=0$: Complexity is suppressed with only small-amplitude fluctuations (free fermion behavior)
• $h_z>0$: Complexity amplitude gradually increases

Comparative Conclusion: Completely opposite trend from ferromagnetic phase, indicating paramagnetic phase does not support confinement; longitudinal field introduction enhances quantum chaotic dynamics.

3. Critical Point Crossing Quench (Figure 3)

Key Discoveries:

• Complexity reaches orders-of-magnitude larger values (several orders larger than intra-phase quenches)
• Initially increases with $h_z$, then begins to decrease at larger field strengths
• Superimposed with random fluctuations

Physical Interpretation:

• Large complexity: Quench across critical point excites broad continuum, strong delocalization in Krylov space
• Turnover behavior: May signal emergence of weak confinement, though signal is weaker (possibly due to finite meson propagation velocity)

Krylov Spectroscopy: Meson Mass Extraction (Figure 4)

Experimental Parameters: $h_x=0.25, h_z=0.2$

Power Spectrum Peaks:

• High-frequency peaks precisely correspond to semiclassical predictions of meson masses
  • $m_1 = 4.025J$
  • $m_2 = 4.702J$
• Spectrum also captures meson mass spacing $m_{12} = m_2 - m_1$

Multiple Meson Case (Supplementary Figure S2): At $h_x=0.25, h_z=0.1$, spectrum successfully resolves four meson masses:

• $m_1=3.662, m_2=4.127, m_3=4.480, m_4=4.769$
• Mass spacings $m_{12}, m_{13}$ also captured

Verification Accuracy: Power spectrum peak positions deviate from semiclassical calculations within numerical precision.

Finite Size Effects Analysis

Ferromagnetic Phase ($h_z=0$):

• Significant finite size effects (Figure 1 inset a)
• Reflects ballistic propagation of free domain walls

Confined Region ($h_z=0.4$):

• Finite size effects suppressed (Figure 1 inset b)
• Confirms absence of correlation propagation

Related Work

Confinement Research

High-Energy Physics

1. Confinement in QCD: Quarks and gluons cannot be isolated, forming mesons and baryons1-3
2. Non-Abelian Gauge Theories: Confinement is a core feature

Condensed Matter Physics

1. Confinement in Low-Dimensional Systems:
   • Transverse-longitudinal field Ising chain4-6
   • McCoy-Wu model7
   • Long-range interaction spin chains8,9
2. Experimental Observations:
   • Quasi-one-dimensional compounds10-12: CoNb₂O₆, BaCo₂V₂O₈
   • Cold atom systems13-15
   • IBM quantum computer17
3. Pioneering Work by Kormos et al.16:
   • Demonstrated rapid entanglement entropy saturation from longitudinal field
   • Magnetization oscillations encode meson spectrum information

Krylov Complexity Research

Theoretical Development

1. Foundational Framework18-20:
   • Defined as wave function spreading in Krylov subspace
   • Geometric interpretation: particle motion on semi-infinite chain
2. Quantum Chaos Probe21-30:
   • Captures chaos and thermalization features
   • Order parameter for integrable-chaotic transition23
3. Gauge Theory Applications33:
   • Operator complexity as order parameter for confinement-deconfinement transition in large-N theories
   • But state complexity confinement signals after quenches remain unstudied

Related Complexity Measures

• Spread complexity19
• Spectral complexity28
• Bi-Lanczos dynamics for open systems24

Relationship to Related Work

Main Distinctions:

1. First Systematic Study: Krylov state complexity as real-time confinement probe
2. Multiple Quench Protocols: Reveals different behaviors in different regions
3. Spectroscopy Method: New tool for extracting meson masses
4. Advantage Analysis: Theoretical advantages over traditional operator methods

Conclusions and Discussion

Main Conclusions

1. Krylov Complexity is a Sensitive Confinement Probe:
   • Ferromagnetic phase quench: Longitudinal field causes significant suppression
   • Paramagnetic phase quench: Complexity enhancement, no confinement
   • Critical point crossing quench: Extreme complexity, weak confinement trends
2. Krylov Spectroscopy Method:
   • Power spectrum peaks precisely correspond to meson masses
   • High consistency with semiclassical predictions
   • Provides complete energy level structure information
3. Theoretical Advantages:
   • State-specific, operator-independent
   • Scans all accessible frequency components
   • No matrix element vanishing problems

Limitations

1. System Size Constraints:
   • Main results use $L=14$
   • Larger systems computationally expensive (Krylov space dimension grows exponentially)
2. Weak Confinement Signals:
   • Confinement signals in critical point crossing quenches are weaker
   • May require more refined analysis
3. Model Specificity:
   • Only Ising model studied
   • Other confinement models (e.g., gauge theories) await exploration
4. Temperature Effects:
   • Only zero-temperature quenches considered
   • Finite temperature algorithms55 await application

Future Directions

Authors explicitly propose the following open questions:

1. Lattice Gauge Theories:
   • How does Krylov complexity reflect confinement from gauge constraints39-48?
   • Can it capture non-perturbative phenomena like string breaking51,52?
2. Other Non-Perturbative Phenomena:
   • False vacuum decay49,50
   • Quantum many-body scars53,54
3. Operator Growth:
   • Study how confinement affects operator growth in Krylov space
   • Apply finite temperature algorithms to probe near-ground-state phenomena
4. Experimental Implementation:
   • Direct measurement in cold atom systems
   • Verification on quantum simulators
5. Theoretical Extensions:
   • Higher-dimensional systems
   • Long-range interactions
   • Topological effects

In-Depth Evaluation

Strengths

1. Conceptual Innovation (★★★★★)

• Interdisciplinary Bridge: Cleverly connects quantum information (Krylov complexity) with high-energy/condensed matter physics (confinement)
• New Probe Development: First systematic application of Krylov state complexity to confinement research
• Theoretical Depth: Provides geometric and physical intuition for complexity-confinement relationship

2. Methodological Rigor (★★★★★)

• Multi-Protocol Verification: Three different quench protocols comprehensively test complexity behavior
• Semiclassical Comparison: Meson mass calculations provide independent verification
• Numerical Stability: FOL algorithm ensures high precision
• System Size Verification: Finite size effects analysis confirms physical picture

3. Result Convincingness (★★★★★)

• Quantitative Consistency: Exact match between power spectrum peaks and semiclassical predictions (Figure 4)
• Qualitative Contrast Clarity: Opposite trends in ferromagnetic vs. paramagnetic phases clearly demonstrate confinement/non-confinement
• Order-of-Magnitude Differences: Extreme complexity in critical point crossing provides strong signal

4. Physical Insight (★★★★★)

• Unified Understanding: Directly connects complexity suppression to restricted correlation propagation
• Spectroscopy Advantages: Clarifies essential advantages of Krylov method over traditional operators
• Mechanism Explanation: Clear picture of static meson production from increased effective mass

5. Writing Quality (★★★★☆)

• Clear Logic: Well-structured progression from motivation through methods to results
• Quality Figures: Figures 1-4 intuitively present core findings
• Supplementary Materials: Detailed algorithms and additional data enhance reproducibility

Weaknesses

1. System Size Limitations (★★★☆☆)

• Computational Bottleneck: $L=14$ relatively small, possible finite size effects
• Missing Scaling Analysis: Lacks systematic study of complexity scaling with $L$
• Thermodynamic Limit: Cannot confirm behavior as $L\to\infty$

2. Insufficient Critical Point Analysis (★★★☆☆)

• Weak Confinement Signal: Only observes "trends," lacks quantitative characterization
• Random Fluctuations: Physical origin of fluctuations not deeply analyzed
• Incomplete Mechanism: Detailed mechanism for complexity increase-then-decrease insufficiently discussed

3. Limited Model Coverage (★★★☆☆)

• Single Model: Only Ising model studied
• Missing Gauge Theories: Though proposed as future direction, not addressed in this work
• Unknown Universality: Method applicability to other confinement systems requires verification

4. Weak Experimental Connection (★★☆☆☆)

• Missing Measurement Schemes: No discussion of experimental Krylov complexity measurement
• Noise Effects: Experimental noise impact on complexity not considered
• Feasibility Analysis: Specific implementation schemes for cold atoms or quantum simulators lacking

5. Limited Theoretical Depth (★★★☆☆)

• Limited Analytical Results: Primarily numerical, lacks analytical theory
• Causal Relationship: Causal chain between complexity suppression and confinement could be more rigorously established
• Critical Behavior: Scaling behavior near critical point $h_x\to 1$ not studied

Impact Assessment

Contribution to Field (★★★★★)

1. Methodological Contribution:
   • Provides new tool for confinement research
   • Extends Krylov complexity application scope
   • Establishes new quantum information-high-energy physics connection
2. Theoretical Significance:
   • Deepens understanding of restricted quantum information propagation
   • Provides quantum information perspective on non-perturbative phenomena
3. Inspirational Value:
   • Opens multiple follow-up research directions (gauge theories, scars, etc.)
   • May inspire experimental groups to develop new measurement schemes

Practical Value (★★★★☆)

1. Quantum Simulation:
   • Provides theoretical guidance for cold atom experiment design
   • Confinement verification on quantum computers
2. Diagnostic Tool:
   • Can identify confinement in unknown systems
   • New method for meson spectrum extraction
3. Limitations:
   • Currently limited to small-system numerical simulation
   • Experimental implementation pathway unclear

Reproducibility (★★★★☆)

Strengths:

• Simple, clearly specified model
• Standard Lanczos algorithm with detailed description (supplementary materials)
• Complete parameter specification

Improvement Space:

• Code not publicly released (though algorithm is clear)
• Some numerical details (e.g., time evolution method) could be more detailed

Applicable Scenarios

Ideal Application Scenarios

1. Low-Dimensional Quantum Spin Chains:
   • Systems with well-defined confinement potential
   • Sizes amenable to exact numerical simulation
2. Cold Atom Quantum Simulation:
   • Controllable quench protocols
   • Long coherence time evolution
3. Quantum Computer Verification:
   • Small-scale confinement models
   • Algorithm benchmarking

Limiting Scenarios

1. Large Systems/High Dimensions: Exponential computational cost growth
2. Strongly Dissipative Systems: Requires open system generalization24
3. Finite Temperature: Requires new algorithms55
4. Complex Gauge Theories: Method effectiveness requires verification

Technical Highlights

1. Full Orthogonalization Lanczos

• Innovation: Employs FOL rather than standard Lanczos
• Advantage: Avoids numerical orthogonality loss
• Implementation: Repeated Gram-Schmidt, termination at $\beta_n<10^{-12}$

2. Power Spectrum Computation

• Parameter Choice: $t\in[0,100], \Delta t=0.1$ balances resolution and computational cost
• Peak Identification: Clearly resolves meson masses and mass spacings

3. Semiclassical Theory

• Bohr-Sommerfeld Quantization: Provides independent theoretical prediction
• Numerical Solution: Converts to nonlinear equation root-finding

Selected References

Core Related Literature

1. 16 Kormos et al., Nat. Phys. 13, 246 (2017): Pioneering confinement study in Ising model
2. 18 Rabinovici et al., arXiv:2507.06286: Krylov complexity review
3. 33 Anegawa et al., JHEP 2024, 119: Operator complexity in large-N gauge theories
4. 38 Rutkevich, J. Stat. Phys. 131, 917 (2008): Semiclassical theory of meson spectrum

Experimental Related

5. 17 Vovrosh & Knolle, Sci. Rep. 11, 11577 (2021): Confinement on IBM quantum computer
6. 10-12: Experimental observations in quasi-one-dimensional compounds
7. 47 Zhang et al., Nat. Phys. (2024): Topological θ-angle control in cold atoms

Methodology

8. 19,20: Geometric interpretation of Krylov complexity
9. 35,36: Classical Lanczos algorithm literature

Overall Evaluation

Recommended Publication Level

• Suitable Journals: Physical Review Letters, Nature Physics, PRX Quantum
• Rationale: Novel concepts, significant interdisciplinary significance, clear results, opens new research directions

Recommendations for Researchers

1. Follow-up Work: Promptly extend to gauge theories to verify universality
2. Experimental Collaboration: Partner with cold atom groups to develop measurement schemes
3. Theoretical Deepening: Develop analytical theory to understand deep connections between complexity and confinement
4. Open Source Code: Enhance reproducibility and impact