Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2607.17346

The Information Content of Krylov Observables: A Machine Learning Approach

Ritam Basu

hep-thcond-mat.dis-nncond-mat.stat-mechmath-phquant-phSubmitted Jul 19, 2026• Announced Jul 21, 2026• 0 citations• Revised Oct 12, 2026 (v2)
v1 · Jul 19, 2026v2 · Oct 12, 2026
arXivHTMLPDF

Abstract

We employ machine learning techniques to estimate the amount of information contained in three Krylov-space observables: the spread complexity $C(t)$, the discrete Wigner negativity $N(t)$, and the normalized negativity $χ(t)=N(t)/|S(t)|$, where $S(t)$ is the survival amplitude. We study thermofield-double evolutions on GUE, GOE and Poisson spectra, and on the integrable $SL(2,R)$ Krylov chain of a two-dimensional CFT primary, with Lanczos coefficients $b_n=α\sqrt{n(n+2h-1)}$. To introduce chaos we add a GUE random matrix to the integrable chain, $H(\varepsilon)=H_{SL(2,R)}+\varepsilon R_0 W_{GUE}$, where $R_0$ is its spectral half-width. Either $C(t)$ or $N(t)$ alone determines the thermofield temperature, with $R^2\simeq0.999$. However, neither can reconstruct the fine structure of the spectral form factor. In the integrable sector all three observables carry the same information. When we turn on chaos an asymmetry appears: $χ$ predicts $C$ better than $C$ predicts $χ$. This gap grows from $+0.33$ to $+0.77$, while the same gap for the raw negativity decays to zero. The extra information comes from the survival amplitude, and $χ(t)$ carries it while $C(t)$ can still be recovered from $χ(t)$. We derive a lower bound on this gap, set by the fraction of the variance of $χ$ that comes from $S$.

DOI
10.1103/hqts-lp2n
Comments
54 pages, 34 figures, 4 tables. Accepted for publication in Physical Review D, https://journals.aps.org/prd/accepted/10.1103/hqts-lp2n

Related papers

Krylov complexity in mixed phase space

Kyoung-Bum Huh, Hyun-Sik Jeong, Leopoldo A. Pando Zayas +1 more

We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.

arXiv:2412.04963•Physical review. D/Physical review. D.•May 21, 2025•18 citations

Onset of random-matrix statistics in the Yukawa-Sachdev-Ye-Kitaev model: Spectral correlations and Krylov complexity

Sizheng Cao, Hyun-Sik Jeong, Yi-Li Wang

The Yukawa-Sachdev-Ye-Kitaev (YSYK) model, a solvable model of non-Fermi liquids with a proposed cavity-QED realization, couples $N$ complex fermions to $M$ bosons of frequency $ω_0$ through random Yukawa couplings of strength $g$, and its dynamics depends only on $R\equivω_0/g^{2/3}$. For $R\ll1$ the bosons are slow and act as a static random background in which the fermions are free, while for $R\gg1$ they are only virtually excited and generate a random four-fermion interaction of rank $M$, which is generically chaotic for $M\geq2$. Using exact diagonalization at infinite temperature, we follow the emergence of random-matrix statistics between these limits with unfolded spectral statistics and the Krylov complexity of the thermofield-double state. For complex couplings the level-spacing distribution evolves from Poisson to Gaussian-unitary-ensemble (GUE) statistics, the spectral form factor develops a correlation hole and a ramp, and the Krylov complexity evaluated on the unfolded spectrum develops a peak above its late-time plateau whose height reaches the GUE value. All diagnostics locate the onset around $R\approx0.3$, and adding boson modes moves it to smaller $R$. This delineates the parameter regime in which a cavity-QED realization would display random-matrix spectral correlations. On the other hand, real couplings give the same crossover towards Gaussian-orthogonal-ensemble (GOE) statistics, and along an interpolation between the two the peak height follows the mean gap ratio, so the peak identifies the symmetry class as well as the onset of chaos. Without unfolding the peak stays about 20$\%$ below its random-matrix value in both classes, which we trace to the non-semicircular density of states, and at larger $R$ the boson-number bands slow down the spreading of the state in Krylov space.

nlin.CDCross-listedarXiv:2609.30333•Sep 25, 2026•0 citations

Krylov complexity of density matrix operators

Paweł Caputa, Hyun-Sik Jeong, Sinong Liu +2 more

A bstract Quantifying complexity in quantum systems has witnessed a surge of interest in recent years, with Krylov-based measures such as Krylov complexity ( C K ) and Spread complexity ( C S ) gaining prominence. In this study, we investigate their interplay by considering the complexity of states represented by density matrix operators . After setting up the problem, we analyze a handful of analytical and numerical examples spanning generic two-dimensional Hilbert spaces, qubit states, quantum harmonic oscillators, and random matrix theories, uncovering insightful relationships. For generic pure states, our analysis reveals two key findings: (I) a correspondence between moment-generating functions (of Lanczos coefficients) and survival amplitudes, and (II) an early-time equivalence between C K and 2 C S . Furthermore, for maximally entangled pure states, we find that the moment-generating function of C K becomes the Spectral Form Factor and, at late-times, C K is simply related to NC S for N ≥ 2 within the N -dimensional Hilbert space. Notably, we confirm that C K = 2 C S holds across all times when N = 2. Through the lens of random matrix theories, we also discuss deviations between complexities at intermediate times and highlight subtleties in the averaging approach at the level of the survival amplitude.

arXiv:2402.09522•Journal of High Energy Physics•May 31, 2024•59 citations

On Krylov Complexity

Adrián Sánchez-Garrido

This Thesis explores the notion of Krylov complexity as a probe of quantum chaos and as a candidate for holographic complexity. The first Part is devoted to presenting the fundamental notions required to conduct research in this area. Namely, an extensive introduction to the Lanczos algorithm, its properties and associated algebraic structures, as well as technical details related to its practical implementation, is given. Subsequently, an overview of the seminal references and the main debates regarding Krylov complexity and its relation to chaos and holography is provided. The text throughout this first Part combines review material with original analyses which either intend to contextualize, compare and criticize results in the literature, or are the fruit of the investigations leading to the publications on which this Thesis is based. These research projects are the subject of the second Part of the manuscript. In them, methods for the efficient implementation of the Lanczos algorithm in finite many-body systems were developed, allowing to compute numerically the Krylov complexity of models like SYK or the XXZ spin chain up to time scales exponentially large in system size. It was observed that the operator Krylov complexity profile in SYK, a paradigmatic low-dimensional chaotic system with a holographic dual, agrees with holographic expectations, while in the case of integrable models like XXZ complexity is affected by a novel localization effect in the so-called Krylov space which hinders its growth. Finally, an exact, analytical, correspondence between the Krylov complexity of the infinite-temperature thermofield double state in the low-energy regime of the double-scaled SYK model and bulk length in the theory of JT gravity is established.

arXiv:2407.03866•Jul 4, 2024•0 citations

Krylov complexity in the IP matrix model II

Norihiro Iizuka, Mitsuhiro Nishida

We continue the analysis of the Krylov complexity in the IP matrix model. In a previous paper, for a fundamental operator, it was shown that at zero temperature, the Krylov complexity oscillates and does not grow, but in the infinite temperature limit, the Krylov complexity grows exponentially in time as $\sim \exp\left( {\mathcal{O}\left( {\sqrt{t}}\right)} \right)$. We study how the Krylov complexity changes from a zero-temperature oscillation to an infinite-temperature exponential growth. At low temperatures, the spectral density is approximated as collections of infinite Wigner semicircles. We showed that this infinite collection of branch cuts yields linear growth to the Lanczos coefficients and gives exponential growth of the Krylov complexity. Thus the IP model for any nonzero temperature shows exponential growth for the Krylov complexity even though the Green function decays by a power law in time. We also study the Lanczos coefficients and the Krylov complexity in the IOP matrix model taking into account the $1/N^2$ corrections. There, the Lanczos coefficients are constants and the Krylov complexity does not grow exponentially as expected.

arXiv:2308.07567•Aug 15, 2023•0 citations

Complexity Growth and the Krylov-Wigner function

Ritam Basu, Anirban Ganguly, Nath, Souparna +1 more

For any state in a $D$-dimensional Hilbert space with a choice of basis, one can define a discrete version of the Wigner function -- a quasi-probability distribution which represents the state on a discrete phase space. The Wigner function can, in general, take on negative values, and the amount of negativity in the Wigner function has an operational meaning as a resource for quantum computation. In this note, we study the growth of Wigner negativity for a generic initial state under time evolution with chaotic Hamiltonians. We introduce the Krylov-Wigner function, i.e., the Wigner function defined with respect to the Krylov basis (with appropriate phases), and show that this choice of basis minimizes the early time growth of Wigner negativity in the large $D$ limit. We take this as evidence that the Krylov basis (with appropriate phases) is ideally suited for a dual, semi-classical description of chaotic quantum dynamics at large $D$. We also numerically study the time evolution of the Krylov-Wigner function and its negativity in random matrix theory for an initial pure state. We observe that the negativity rises gradually for a time of $O(D)$ and then saturates close to its upper bound of $\sqrt{D}$.

arXiv:2402.13694•Feb 21, 2024•0 citations
We employ machine learning techniques to estimate the amount of information contained in three Krylov-space observables: the spread complexity C(t)C(t)C(t), the discrete Wigner negativity N(t)N(t)N(t), and the normalized negativity χ(t)=N(t)/∣S(t)∣χ(t)=N(t)/|S(t)|χ(t)=N(t)/∣S(t)∣, where S(t)S(t)S(t) is the survival amplitude. We study thermofield-double evolutions on GUE, GOE and Poisson spectra, and on the integrable SL(2,R)SL(2,R)SL(2,R) Krylov chain of a two-dimensional CFT primary, with Lanczos coefficients bn=αn(n+2h−1)b_n=α\sqrt{n(n+2h-1)}bn​=αn(n+2h−1)​. To introduce chaos we add a GUE random matrix to the integrable chain, H(ε)=HSL(2,R)+εR0WGUEH(\varepsilon)=H_{SL(2,R)}+\varepsilon R_0 W_{GUE}H(ε)=HSL(2,R)​+εR0​WGUE​, where R0R_0R0​ is its spectral half-width. Either C(t)C(t)C(t) or N(t)N(t)N(t) alone determines the thermofield temperature, with R2≃0.999R^2\simeq0.999R2≃0.999. However, neither can reconstruct the fine structure of the spectral form factor. In the integrable sector all three observables carry the same information. When we turn on chaos an asymmetry appears: χχχ predicts CCC better than CCC predicts χχχ. This gap grows from +0.33+0.33+0.33 to +0.77+0.77+0.77, while the same gap for the raw negativity decays to zero. The extra information comes from the survival amplitude, and χ(t)χ(t)χ(t) carries it while C(t)C(t)C(t) can still be recovered from χ(t)χ(t)χ(t). We derive a lower bound on this gap, set by the fraction of the variance of χχχ that comes from SSS.
We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.
The Yukawa-Sachdev-Ye-Kitaev (YSYK) model, a solvable model of non-Fermi liquids with a proposed cavity-QED realization, couples NNN complex fermions to MMM bosons of frequency ω0ω_0ω0​ through random Yukawa couplings of strength ggg, and its dynamics depends only on R≡ω0/g2/3R\equivω_0/g^{2/3}R≡ω0​/g2/3. For R≪1R\ll1R≪1 the bosons are slow and act as a static random background in which the fermions are free, while for R≫1R\gg1R≫1 they are only virtually excited and generate a random four-fermion interaction of rank MMM, which is generically chaotic for M≥2M\geq2M≥2. Using exact diagonalization at infinite temperature, we follow the emergence of random-matrix statistics between these limits with unfolded spectral statistics and the Krylov complexity of the thermofield-double state. For complex couplings the level-spacing distribution evolves from Poisson to Gaussian-unitary-ensemble (GUE) statistics, the spectral form factor develops a correlation hole and a ramp, and the Krylov complexity evaluated on the unfolded spectrum develops a peak above its late-time plateau whose height reaches the GUE value. All diagnostics locate the onset around R≈0.3R\approx0.3R≈0.3, and adding boson modes moves it to smaller RRR. This delineates the parameter regime in which a cavity-QED realization would display random-matrix spectral correlations. On the other hand, real couplings give the same crossover towards Gaussian-orthogonal-ensemble (GOE) statistics, and along an interpolation between the two the peak height follows the mean gap ratio, so the peak identifies the symmetry class as well as the onset of chaos. Without unfolding the peak stays about 20%\%% below its random-matrix value in both classes, which we trace to the non-semicircular density of states, and at larger RRR the boson-number bands slow down the spreading of the state in Krylov space.
A bstract Quantifying complexity in quantum systems has witnessed a surge of interest in recent years, with Krylov-based measures such as Krylov complexity ( C K ) and Spread complexity ( C S ) gaining prominence. In this study, we investigate their interplay by considering the complexity of states represented by density matrix operators . After setting up the problem, we analyze a handful of analytical and numerical examples spanning generic two-dimensional Hilbert spaces, qubit states, quantum harmonic oscillators, and random matrix theories, uncovering insightful relationships. For generic pure states, our analysis reveals two key findings: (I) a correspondence between moment-generating functions (of Lanczos coefficients) and survival amplitudes, and (II) an early-time equivalence between C K and 2 C S . Furthermore, for maximally entangled pure states, we find that the moment-generating function of C K becomes the Spectral Form Factor and, at late-times, C K is simply related to NC S for N ≥ 2 within the N -dimensional Hilbert space. Notably, we confirm that C K = 2 C S holds across all times when N = 2. Through the lens of random matrix theories, we also discuss deviations between complexities at intermediate times and highlight subtleties in the averaging approach at the level of the survival amplitude.
This Thesis explores the notion of Krylov complexity as a probe of quantum chaos and as a candidate for holographic complexity. The first Part is devoted to presenting the fundamental notions required to conduct research in this area. Namely, an extensive introduction to the Lanczos algorithm, its properties and associated algebraic structures, as well as technical details related to its practical implementation, is given. Subsequently, an overview of the seminal references and the main debates regarding Krylov complexity and its relation to chaos and holography is provided. The text throughout this first Part combines review material with original analyses which either intend to contextualize, compare and criticize results in the literature, or are the fruit of the investigations leading to the publications on which this Thesis is based. These research projects are the subject of the second Part of the manuscript. In them, methods for the efficient implementation of the Lanczos algorithm in finite many-body systems were developed, allowing to compute numerically the Krylov complexity of models like SYK or the XXZ spin chain up to time scales exponentially large in system size. It was observed that the operator Krylov complexity profile in SYK, a paradigmatic low-dimensional chaotic system with a holographic dual, agrees with holographic expectations, while in the case of integrable models like XXZ complexity is affected by a novel localization effect in the so-called Krylov space which hinders its growth. Finally, an exact, analytical, correspondence between the Krylov complexity of the infinite-temperature thermofield double state in the low-energy regime of the double-scaled SYK model and bulk length in the theory of JT gravity is established.
We continue the analysis of the Krylov complexity in the IP matrix model. In a previous paper, for a fundamental operator, it was shown that at zero temperature, the Krylov complexity oscillates and does not grow, but in the infinite temperature limit, the Krylov complexity grows exponentially in time as ∼exp⁡(O(t))\sim \exp\left( {\mathcal{O}\left( {\sqrt{t}}\right)} \right)∼exp(O(t​)). We study how the Krylov complexity changes from a zero-temperature oscillation to an infinite-temperature exponential growth. At low temperatures, the spectral density is approximated as collections of infinite Wigner semicircles. We showed that this infinite collection of branch cuts yields linear growth to the Lanczos coefficients and gives exponential growth of the Krylov complexity. Thus the IP model for any nonzero temperature shows exponential growth for the Krylov complexity even though the Green function decays by a power law in time. We also study the Lanczos coefficients and the Krylov complexity in the IOP matrix model taking into account the 1/N21/N^21/N2 corrections. There, the Lanczos coefficients are constants and the Krylov complexity does not grow exponentially as expected.
For any state in a DDD-dimensional Hilbert space with a choice of basis, one can define a discrete version of the Wigner function -- a quasi-probability distribution which represents the state on a discrete phase space. The Wigner function can, in general, take on negative values, and the amount of negativity in the Wigner function has an operational meaning as a resource for quantum computation. In this note, we study the growth of Wigner negativity for a generic initial state under time evolution with chaotic Hamiltonians. We introduce the Krylov-Wigner function, i.e., the Wigner function defined with respect to the Krylov basis (with appropriate phases), and show that this choice of basis minimizes the early time growth of Wigner negativity in the large DDD limit. We take this as evidence that the Krylov basis (with appropriate phases) is ideally suited for a dual, semi-classical description of chaotic quantum dynamics at large DDD. We also numerically study the time evolution of the Krylov-Wigner function and its negativity in random matrix theory for an initial pure state. We observe that the negativity rises gradually for a time of O(D)O(D)O(D) and then saturates close to its upper bound of D\sqrt{D}D​.