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:2601.07937

Attention in Krylov Space: Transformer-Based Extrapolation of Lanczos Coefficients

Zihao Qi, Christopher Earls

quant-phcond-mat.stat-mechSubmitted Jan 12, 2026• Announced Jan 14, 2026• 0 citations• Revised Jul 31, 2026 (v2)
v1 · Jan 12, 2026v2 · Jul 31, 2026
arXivHTMLPDF

Abstract

The Universal Operator Growth Hypothesis formulates time evolution of operators through Lanczos coefficients. In practice, however, numerical instability and memory cost limit the number of coefficients that can be exactly computed. In response to these challenges, the standard approach relies on fitting early coefficients to asymptotic forms, but such procedures can miss subleading, history-dependent structures in the coefficients that subsequently affect reconstructed observables. In this work, we treat the Lanczos coefficients as a causal time sequence and introduce a transformer-based model to autoregressively predict future Lanczos coefficients from short prefixes. For classical and quantum chaotic systems, our model outperforms asymptotic fits in both coefficient extrapolation and physical observable reconstruction, and achieves an order-of-magnitude reduction in error. The model also accurately extrapolates coefficients in integrable regimes, where no universal asymptotic fit exists. Remarkably, our model transfers across system sizes: it can be trained on smaller systems and then be used to extrapolate coefficients on a larger system \emph{without retraining}. By probing the learned attention patterns and performing targeted attention ablations, we identify portions of the coefficient history that are most influential for accurate forecasts. Our results demonstrate that modern sequence models can serve as practical surrogates for probing operator dynamics deep in Krylov space, where brute-force Lanczos iteration can be computationally prohibitive.

Journal reference
Phys. Rev. B 114, 014319 (2026)
Comments
15 pages, 11 figures

Related papers

Exactly solvable models for universal operator growth

Oleksandr Gamayun, Murtaza Ali Mir, Oleg Valentinovich Lychkovskiy +1 more

A bstract Quantum observables of generic many-body systems exhibit a universal pattern of growth in the Krylov space of operators. This pattern becomes particularly manifest in the Lanczos basis, where the evolution superoperator assumes the tridiagonal form. According to the universal operator growth hypothesis, the nonzero elements of the superoperator, known as Lanczos coefficients, grow asymptotically linearly. We introduce and explore broad families of Lanczos coefficients that are consistent with the universal operator growth and lead to the exactly solvable dynamics. Within these families, the subleading terms of asymptotic expansion of the Lanczos sequence can be controlled and fine-tuned to produce diverse dynamical patterns. For one of the families, the Krylov complexity is computed exactly.

arXiv:2504.03435•Journal of High Energy Physics•Jul 25, 2025•8 citations

Krylov complexity in saddle-dominated scrambling

Budhaditya Bhattacharjee, Xiangyu Cao, Pratik Nandy +1 more

A bstract In semi-classical systems, the exponential growth of the out-of-time-order correlator (OTOC) is believed to be the hallmark of quantum chaos. However, on several occasions, it has been argued that, even in integrable systems, OTOC can grow exponentially due to the presence of unstable saddle points in the phase space. In this work, we probe such an integrable system exhibiting saddle-dominated scrambling through Krylov complexity and the associated Lanczos coefficients. In the realm of the universal operator growth hypothesis, we demonstrate that the Lanczos coefficients follow the linear growth, which ensures the exponential behavior of Krylov complexity at early times. The linear growth arises entirely due to the saddle, which dominates other phase-space points even away from itself. Our results reveal that the exponential growth of Krylov complexity can be observed in integrable systems with saddle-dominated scrambling and thus need not be associated with the presence of chaos.

arXiv:2203.03534•Journal of High Energy Physics•May 25, 2022•138 citations

A Schwinger-Keldysh Formulation of Semiclassical Operator Dynamics

Jeff Murugan, Hendrik J. R. van Zyl

In this work we develop a real-time Schwinger-Keldysh formulation of Krylov dynamics that treats Krylov complexity as an in-in observable generated by a closed time contour path integral. The resulting generating functional exposes an emergent phase-space description in which the Lanczos coefficients define an effective Hamiltonian governing operator motion along the Krylov chain. In the semiclassical limit, exponential complexity growth arises from hyperbolic trajectories, and asymptotically linear Lanczos growth appears as a universal chaotic fixed point, with sub-leading deformations classified as irrelevant, marginal or relevant. Going beyond the saddle, the Schwinger-Keldysh framework provides controlled access to fluctuations and large deviations of Krylov complexity, revealing sharp signatures of integrability-chaos crossovers that are invisible at the level of the mean. This formulation reorganises Krylov complexity into a dynamical field-theoretic framework and identifies new fluctuation diagnostics of operator growth in closed quantum systems.

quant-pharXiv:2602.02106•Feb 3, 2026•0 citations

Krylov complexity of many-body localization: Operator localization in Krylov basis

Fabian Ballar Trigueros

We study the operator growth problem and its complexity in the many-body localization (MBL) system from the Lanczos algorithm perspective. Using the Krylov basis, the operator growth problem can be viewed as a single-particle hopping problem on a semi-infinite chain with the hopping amplitudes given by the Lanczos coefficients. We find that, in the MBL systems, the Lanczos coefficients scale as ∼ n/ ln(n) asymptotically, same as in the ergodic systems, but with an additional even-odd alteration and an effective randomness. We use a simple linear extrapolation scheme as an attempt to extrapolate the Lanczos coefficients to the thermodynamic limit. With the original and extrapolated Lanczos coefficients, we study the properties of the emergent single-particle hopping problem via its spectral function, integrals of motion, Krylov complexity, wavefunction profile and return probability. Our numerical results of the above quantities suggest that the emergent single-particle hopping problem in the MBL system is localized when initialized on the first site. We also study the operator growth in the MBL phenomenological model, whose Lanczos coefficients also have an even-odd alteration, but approach constants asymptotically. The Krylov complexity grows linearly in time in this case.

•Aug 1, 2022•81 citations

Long-Range Pairing in the Kitaev Model: Krylov Subspace Signatures

Rishabh Jha, Heiko Georg Menzler

Krylov subspace methods quantify operator growth in quantum many-body systems through Lanczos coefficients that encode how operators spread under time evolution. Although these diagnostics were originally motivated by questions of chaos and integrability, quadratic fermionic Hamiltonians are often expected to exhibit trivial Lanczos structure. Here we show that, in the long-range Kitaev chain, Lanczos coefficients generated from local boundary operators sharply diagnose whether the lowest excitation gap is controlled by boundary-localized or bulk-extended modes. We introduce the $Krylov$ $staggering$ $parameter$ for the Lanczos coefficients. In the short-range Kitaev chain with balanced hopping and pairing, we derive analytically for arbitrary system size (valid in the thermodynamic limit) and show that this quantity is exactly constant and its sign cleanly distinguishes the topological phase with Majorana edge modes from the trivial phase. Away from that limit, long-range couplings and pairing-hopping imbalance deform the simple flat structure and analytical control is lost, nevertheless, we show that the sign pattern of the diagnostic still tracks whether the lowest excitation gap is controlled by boundary modes or by bulk excitations. These results are enabled by an exact single-particle operator Lanczos algorithm, as derived in this work, which reduces the recursion from exponentially large operator space to a finite-dimensional linear problem and achieves machine precision for chains of hundreds of sites. Krylov diagnostics thus emerge as practical probes of boundary-versus-bulk low-energy physics in topological superconductors with broken $U(1)$ symmetry and algebraically decaying couplings.

quant-pharXiv:2602.11278•Journal of Physics A Mathematical and Theoretical•Feb 13, 2026•0 citations

Krylov complexity of many-body localization: Operator localization in Krylov basis

Fabian Ballar Trigueros, Cheng-Ju Lin

We study the operator growth problem and its complexity in the many-body localization (MBL) system from the Lanczos algorithm perspective. Using the Krylov basis, the operator growth problem can be viewed as a single-particle hopping problem on a semi-infinite chain with the hopping amplitudes given by the Lanczos coefficients. We find that, in the MBL systems, the Lanczos coefficients scale as \sim n/\ln(n) ∼ n / ln ( n ) asymptotically, same as in the ergodic systems, but with an additional even-odd alteration and an effective randomness. We use a simple linear extrapolation scheme as an attempt to extrapolate the Lanczos coefficients to the thermodynamic limit. With the original and extrapolated Lanczos coefficients, we study the properties of the emergent single-particle hopping problem via its spectral function, integrals of motion, Krylov complexity, wavefunction profile and return probability. Our numerical results of the above quantities suggest that the emergent single-particle hopping problem in the MBL system is localized when initialized on the first site. We also study the operator growth in the MBL phenomenological model, whose Lanczos coefficients also have an even-odd alteration, but approach constants asymptotically. The Krylov complexity grows linearly in time in this case.

arXiv:2112.04722•SciPost Physics•Aug 29, 2022•0 citations
The Universal Operator Growth Hypothesis formulates time evolution of operators through Lanczos coefficients. In practice, however, numerical instability and memory cost limit the number of coefficients that can be exactly computed. In response to these challenges, the standard approach relies on fitting early coefficients to asymptotic forms, but such procedures can miss subleading, history-dependent structures in the coefficients that subsequently affect reconstructed observables. In this work, we treat the Lanczos coefficients as a causal time sequence and introduce a transformer-based model to autoregressively predict future Lanczos coefficients from short prefixes. For classical and quantum chaotic systems, our model outperforms asymptotic fits in both coefficient extrapolation and physical observable reconstruction, and achieves an order-of-magnitude reduction in error. The model also accurately extrapolates coefficients in integrable regimes, where no universal asymptotic fit exists. Remarkably, our model transfers across system sizes: it can be trained on smaller systems and then be used to extrapolate coefficients on a larger system \emph{without retraining}. By probing the learned attention patterns and performing targeted attention ablations, we identify portions of the coefficient history that are most influential for accurate forecasts. Our results demonstrate that modern sequence models can serve as practical surrogates for probing operator dynamics deep in Krylov space, where brute-force Lanczos iteration can be computationally prohibitive.
A bstract Quantum observables of generic many-body systems exhibit a universal pattern of growth in the Krylov space of operators. This pattern becomes particularly manifest in the Lanczos basis, where the evolution superoperator assumes the tridiagonal form. According to the universal operator growth hypothesis, the nonzero elements of the superoperator, known as Lanczos coefficients, grow asymptotically linearly. We introduce and explore broad families of Lanczos coefficients that are consistent with the universal operator growth and lead to the exactly solvable dynamics. Within these families, the subleading terms of asymptotic expansion of the Lanczos sequence can be controlled and fine-tuned to produce diverse dynamical patterns. For one of the families, the Krylov complexity is computed exactly.
A bstract In semi-classical systems, the exponential growth of the out-of-time-order correlator (OTOC) is believed to be the hallmark of quantum chaos. However, on several occasions, it has been argued that, even in integrable systems, OTOC can grow exponentially due to the presence of unstable saddle points in the phase space. In this work, we probe such an integrable system exhibiting saddle-dominated scrambling through Krylov complexity and the associated Lanczos coefficients. In the realm of the universal operator growth hypothesis, we demonstrate that the Lanczos coefficients follow the linear growth, which ensures the exponential behavior of Krylov complexity at early times. The linear growth arises entirely due to the saddle, which dominates other phase-space points even away from itself. Our results reveal that the exponential growth of Krylov complexity can be observed in integrable systems with saddle-dominated scrambling and thus need not be associated with the presence of chaos.
In this work we develop a real-time Schwinger-Keldysh formulation of Krylov dynamics that treats Krylov complexity as an in-in observable generated by a closed time contour path integral. The resulting generating functional exposes an emergent phase-space description in which the Lanczos coefficients define an effective Hamiltonian governing operator motion along the Krylov chain. In the semiclassical limit, exponential complexity growth arises from hyperbolic trajectories, and asymptotically linear Lanczos growth appears as a universal chaotic fixed point, with sub-leading deformations classified as irrelevant, marginal or relevant. Going beyond the saddle, the Schwinger-Keldysh framework provides controlled access to fluctuations and large deviations of Krylov complexity, revealing sharp signatures of integrability-chaos crossovers that are invisible at the level of the mean. This formulation reorganises Krylov complexity into a dynamical field-theoretic framework and identifies new fluctuation diagnostics of operator growth in closed quantum systems.
We study the operator growth problem and its complexity in the many-body localization (MBL) system from the Lanczos algorithm perspective. Using the Krylov basis, the operator growth problem can be viewed as a single-particle hopping problem on a semi-infinite chain with the hopping amplitudes given by the Lanczos coefficients. We find that, in the MBL systems, the Lanczos coefficients scale as ∼ n/ ln(n) asymptotically, same as in the ergodic systems, but with an additional even-odd alteration and an effective randomness. We use a simple linear extrapolation scheme as an attempt to extrapolate the Lanczos coefficients to the thermodynamic limit. With the original and extrapolated Lanczos coefficients, we study the properties of the emergent single-particle hopping problem via its spectral function, integrals of motion, Krylov complexity, wavefunction profile and return probability. Our numerical results of the above quantities suggest that the emergent single-particle hopping problem in the MBL system is localized when initialized on the first site. We also study the operator growth in the MBL phenomenological model, whose Lanczos coefficients also have an even-odd alteration, but approach constants asymptotically. The Krylov complexity grows linearly in time in this case.
Krylov subspace methods quantify operator growth in quantum many-body systems through Lanczos coefficients that encode how operators spread under time evolution. Although these diagnostics were originally motivated by questions of chaos and integrability, quadratic fermionic Hamiltonians are often expected to exhibit trivial Lanczos structure. Here we show that, in the long-range Kitaev chain, Lanczos coefficients generated from local boundary operators sharply diagnose whether the lowest excitation gap is controlled by boundary-localized or bulk-extended modes. We introduce the KrylovKrylovKrylov staggeringstaggeringstaggering parameterparameterparameter for the Lanczos coefficients. In the short-range Kitaev chain with balanced hopping and pairing, we derive analytically for arbitrary system size (valid in the thermodynamic limit) and show that this quantity is exactly constant and its sign cleanly distinguishes the topological phase with Majorana edge modes from the trivial phase. Away from that limit, long-range couplings and pairing-hopping imbalance deform the simple flat structure and analytical control is lost, nevertheless, we show that the sign pattern of the diagnostic still tracks whether the lowest excitation gap is controlled by boundary modes or by bulk excitations. These results are enabled by an exact single-particle operator Lanczos algorithm, as derived in this work, which reduces the recursion from exponentially large operator space to a finite-dimensional linear problem and achieves machine precision for chains of hundreds of sites. Krylov diagnostics thus emerge as practical probes of boundary-versus-bulk low-energy physics in topological superconductors with broken U(1)U(1)U(1) symmetry and algebraically decaying couplings.
We study the operator growth problem and its complexity in the many-body localization (MBL) system from the Lanczos algorithm perspective. Using the Krylov basis, the operator growth problem can be viewed as a single-particle hopping problem on a semi-infinite chain with the hopping amplitudes given by the Lanczos coefficients. We find that, in the MBL systems, the Lanczos coefficients scale as \sim n/\ln(n) ∼ n / ln ( n ) asymptotically, same as in the ergodic systems, but with an additional even-odd alteration and an effective randomness. We use a simple linear extrapolation scheme as an attempt to extrapolate the Lanczos coefficients to the thermodynamic limit. With the original and extrapolated Lanczos coefficients, we study the properties of the emergent single-particle hopping problem via its spectral function, integrals of motion, Krylov complexity, wavefunction profile and return probability. Our numerical results of the above quantities suggest that the emergent single-particle hopping problem in the MBL system is localized when initialized on the first site. We also study the operator growth in the MBL phenomenological model, whose Lanczos coefficients also have an even-odd alteration, but approach constants asymptotically. The Krylov complexity grows linearly in time in this case.