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.13255 · Physical review. B./Physical review. B

The Anderson impurity model from a Krylov perspective: Lanczos coefficients in a quadratic model

Merlin Füllgraf, Jiaozi Wang, Jochen Gemmer, Stefan Kehrein

cond-mat.str-elcond-mat.stat-mechquant-phSubmitted Jan 19, 2026• Announced Jan 20, 2026• 0 citations
v1 · Jan 19, 2026v2 · Apr 13, 2026v3 · May 11, 2026
arXivHTMLPDF

Abstract

We study the Lanczos coefficients in a quadratic model given by an impurity interacting with a multi-mode field of fermions, also known as single impurity Anderson model. We analytically derive closed expressions for the Lanczos coefficients of Majorana fermion operators of the impurity for different structures of the coupling to the hybridisation band at zero temperature. While the model remains quadratic, we find that the growth of the Lanczos coefficients structurally depends strongly on the chosen coupling. Concretely, we find $(i)$ approximately constant, $(ii)$ exactly constant, $(iii)$ square root-like as well $(iv)$ linear growth in the same model. We further argue that in fact through suitably chosen couplings, essentially arbitrary Lanczos coefficients can be obtained in this model. These altogether evince the inadequacy of the Lanczos coefficients as a reliable criterion for classifying the integrability or chaoticity of the systems. Eventually, in the wide-band limit, we find exponential decay of autocorrelation functions in all the settings $(i)-(iv)$, which demonstrates the different structures of the Lanczos coefficients not being indicative of different physical behavior.

Journal reference
Phys. Rev. B 113, 184301 (2026)
DOI
10.1103/l9zp-kg96
Comments
6 pages, 4 figures

Related papers

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

Nonequilibrium quantum many-body physics from perspectives of operator dynamics and the Lanczos picture

Merlin Füllgraf

This cumulative doctoral thesis, based on the works [P1–P5], addresses various topical concepts as well as foundational questions in the nonequilibrium physics of quantum many-body systems. Publication [P1] concerns itself with the question of time’s arrow in quantum mechanics. Based on the recursion method, we attribute autocorrelation functions of few-body observables at infinite temperature (almost) monotonously decreasing functions, constructed from temporal derivatives of the correlation function that indicate a directed approach to equilibrium. We numerically find that in chaotic models, these constructions comprise a notion of simplicity of the underlying dynamics. Addressing the ubiquity of simple relaxation dynamics in otherwise complex chaotic systems, in Ref. [P3] we introduce a scheme based on Lanczos coefficients to approximate correlation functions in the thermodynamic limit. Examining several pertinent condensed-matter type models and benchmarking with state-of-the-art numerics, we find the method to yield good results already for a low number of included Lanczos coefficients, provided that the latter exhibit a smooth and increasing structure. We underpin these findings with analytical arguments on the method’s applicability. Ref. [P2] provides an estimate on the equilibration time of such infinite-temperature correlation functions in the thermodynamic limit. Similarly to Ref. [P3], we find the estimate to capture the equilibration times very reasonably, given that the corresponding Lanczos coefficients exhibit a smooth growth. Substantiated by analytical considerations, the convergence of the estimate implies that equilibration timescales of the relaxation processes are finite and conceivably short. In Ref. [P5] we scrutinize the interpretative power of the Lanczos coefficients concerning the system’s integrability or chaoticity, operator growth and the indicativeness of physical behavior. To this end, we analytically deter- mine the Lanczos coefficients for a Majorana observable in a strictly quadratic model and illustrate that within this model any structure of the Lanczos coefficients may be realised. Different to the other works, Ref. [P4] is not cast in the language of the recursion method. Focusing on the full eigenstate thermalization hypothesis, it explores the connection between thermalization dynamics of multi-point correlation functions and the scaling of thermal free cumulants in setups mimicking a system-bath situation.

•Jul 6, 2026•0 citations

Krylov complexity in Calabi–Yau quantum mechanics

Bao-ning Du, Min-xin Huang

Recently, a novel measure for the complexity of operator growth is proposed based on Lanczos algorithm and Krylov recursion method. We study this Krylov complexity in quantum mechanical systems derived from some well-known local toric Calabi–Yau geometries, as well as some nonrelativistic models. We find that for the Calabi–Yau models, the Lanczos coefficients grow slower than linearly for small [Formula: see text]’s, consistent with the behavior of integrable models. On the other hand, for the nonrelativistic models, the Lanczos coefficients initially grow linearly for small [Formula: see text]’s, then reach a plateau. Although this looks like the behavior of a chaotic system, it is mostly likely due to saddle-dominated scrambling effects instead, as argued in the literature. In our cases, the slopes of linearly growing Lanczos coefficients almost saturate a bound by the temperature. During our study, we also provide an alternative general derivation of the bound for the slope.

arXiv:2212.02926•International Journal of Modern Physics A•Aug 19, 2023•7 citations

Krylov complexity of fermion chain in double-scaled SYK and power spectrum perspective

Takanori Anegawa, Ryota Watanabe

A bstract We investigate Krylov complexity of the fermion chain operator which consists of multiple Majorana fermions in the double-scaled SYK (DSSYK) model with finite temperature. Using the fact that Krylov complexity is computable from two-point functions, the analysis is performed in the limit where the two-point function becomes simple and we compare the results with those of other previous studies. We confirm the exponential growth of Krylov complexity in the very low temperature regime. In general, Krylov complexity grows at most linearly at very late times in any system with a bounded energy spectrum. Therefore, we have to focus on the initial growth to see differences in the behaviors of systems or operators. Since the DSSYK model is such a bounded system, its chaotic nature can be expected to appear as the initial exponential growth of the Krylov complexity. In particular, the time at which the initial exponential growth of Krylov complexity terminates is independent of the number of degrees of freedom. More generally, and not limited to the DSSYK model, we systematically and specifically study the Lanczos coefficients and Krylov complexity using a toy power spectrum and deepen our understanding of those initial behaviors. In particular, we confirm that the overall sech-like behavior of the power spectrum shows the initial linear growth of the Lanczos coefficient, even when the energy spectrum is bounded.

arXiv:2407.13293•Journal of High Energy Physics•Nov 6, 2024•6 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

Isospectrality and Operator Complexity

Pradip Kattel, Yicheng Tang, Natan Andrei

We study a pair of exactly solvable, isospectral fermion chains, one strongly interacting and one quadratic, that nevertheless display remarkably different phase structures and operator dynamics. A nonlocal nonlinear unitary transformation maps one onto the other while preserving the entire many-body spectrum and converting local fermion operators into extended many-body strings. Thus, operators that evolve within a closed linear subspace in the quadratic model become interacting operators that generate increasingly higher-body terms and exhibit asymptotic Lanczos growth $b_n\propto\sqrt n$. Despite their identical spectra, the two models realize distinct phases and sharply different notions of operator complexity. Our results demonstrate that free many-body spectra and interacting operator dynamics are fundamentally compatible.

quant-pharXiv:2606.05294•Jun 5, 2026•0 citations
We study the Lanczos coefficients in a quadratic model given by an impurity interacting with a multi-mode field of fermions, also known as single impurity Anderson model. We analytically derive closed expressions for the Lanczos coefficients of Majorana fermion operators of the impurity for different structures of the coupling to the hybridisation band at zero temperature. While the model remains quadratic, we find that the growth of the Lanczos coefficients structurally depends strongly on the chosen coupling. Concretely, we find (i)(i)(i) approximately constant, (ii)(ii)(ii) exactly constant, (iii)(iii)(iii) square root-like as well (iv)(iv)(iv) linear growth in the same model. We further argue that in fact through suitably chosen couplings, essentially arbitrary Lanczos coefficients can be obtained in this model. These altogether evince the inadequacy of the Lanczos coefficients as a reliable criterion for classifying the integrability or chaoticity of the systems. Eventually, in the wide-band limit, we find exponential decay of autocorrelation functions in all the settings (i)−(iv)(i)-(iv)(i)−(iv), which demonstrates the different structures of the Lanczos coefficients not being indicative of different physical behavior.
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.
This cumulative doctoral thesis, based on the works [P1–P5], addresses various topical concepts as well as foundational questions in the nonequilibrium physics of quantum many-body systems. Publication [P1] concerns itself with the question of time’s arrow in quantum mechanics. Based on the recursion method, we attribute autocorrelation functions of few-body observables at infinite temperature (almost) monotonously decreasing functions, constructed from temporal derivatives of the correlation function that indicate a directed approach to equilibrium. We numerically find that in chaotic models, these constructions comprise a notion of simplicity of the underlying dynamics. Addressing the ubiquity of simple relaxation dynamics in otherwise complex chaotic systems, in Ref. [P3] we introduce a scheme based on Lanczos coefficients to approximate correlation functions in the thermodynamic limit. Examining several pertinent condensed-matter type models and benchmarking with state-of-the-art numerics, we find the method to yield good results already for a low number of included Lanczos coefficients, provided that the latter exhibit a smooth and increasing structure. We underpin these findings with analytical arguments on the method’s applicability. Ref. [P2] provides an estimate on the equilibration time of such infinite-temperature correlation functions in the thermodynamic limit. Similarly to Ref. [P3], we find the estimate to capture the equilibration times very reasonably, given that the corresponding Lanczos coefficients exhibit a smooth growth. Substantiated by analytical considerations, the convergence of the estimate implies that equilibration timescales of the relaxation processes are finite and conceivably short. In Ref. [P5] we scrutinize the interpretative power of the Lanczos coefficients concerning the system’s integrability or chaoticity, operator growth and the indicativeness of physical behavior. To this end, we analytically deter- mine the Lanczos coefficients for a Majorana observable in a strictly quadratic model and illustrate that within this model any structure of the Lanczos coefficients may be realised. Different to the other works, Ref. [P4] is not cast in the language of the recursion method. Focusing on the full eigenstate thermalization hypothesis, it explores the connection between thermalization dynamics of multi-point correlation functions and the scaling of thermal free cumulants in setups mimicking a system-bath situation.
Recently, a novel measure for the complexity of operator growth is proposed based on Lanczos algorithm and Krylov recursion method. We study this Krylov complexity in quantum mechanical systems derived from some well-known local toric Calabi–Yau geometries, as well as some nonrelativistic models. We find that for the Calabi–Yau models, the Lanczos coefficients grow slower than linearly for small [Formula: see text]’s, consistent with the behavior of integrable models. On the other hand, for the nonrelativistic models, the Lanczos coefficients initially grow linearly for small [Formula: see text]’s, then reach a plateau. Although this looks like the behavior of a chaotic system, it is mostly likely due to saddle-dominated scrambling effects instead, as argued in the literature. In our cases, the slopes of linearly growing Lanczos coefficients almost saturate a bound by the temperature. During our study, we also provide an alternative general derivation of the bound for the slope.
A bstract We investigate Krylov complexity of the fermion chain operator which consists of multiple Majorana fermions in the double-scaled SYK (DSSYK) model with finite temperature. Using the fact that Krylov complexity is computable from two-point functions, the analysis is performed in the limit where the two-point function becomes simple and we compare the results with those of other previous studies. We confirm the exponential growth of Krylov complexity in the very low temperature regime. In general, Krylov complexity grows at most linearly at very late times in any system with a bounded energy spectrum. Therefore, we have to focus on the initial growth to see differences in the behaviors of systems or operators. Since the DSSYK model is such a bounded system, its chaotic nature can be expected to appear as the initial exponential growth of the Krylov complexity. In particular, the time at which the initial exponential growth of Krylov complexity terminates is independent of the number of degrees of freedom. More generally, and not limited to the DSSYK model, we systematically and specifically study the Lanczos coefficients and Krylov complexity using a toy power spectrum and deepen our understanding of those initial behaviors. In particular, we confirm that the overall sech-like behavior of the power spectrum shows the initial linear growth of the Lanczos coefficient, even when the energy spectrum is bounded.
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.
We study a pair of exactly solvable, isospectral fermion chains, one strongly interacting and one quadratic, that nevertheless display remarkably different phase structures and operator dynamics. A nonlocal nonlinear unitary transformation maps one onto the other while preserving the entire many-body spectrum and converting local fermion operators into extended many-body strings. Thus, operators that evolve within a closed linear subspace in the quadratic model become interacting operators that generate increasingly higher-body terms and exhibit asymptotic Lanczos growth bn∝nb_n\propto\sqrt nbn​∝n​. Despite their identical spectra, the two models realize distinct phases and sharply different notions of operator complexity. Our results demonstrate that free many-body spectra and interacting operator dynamics are fundamentally compatible.