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

Anderson localization: A view from Krylov space

J. Clayton Peacock, Vadim Oganesyan, Dries Sels

Published Feb 9, 2026• 2 citations
arXivHTMLPDF

Abstract

The Krylov subspace expansion is a workhorse method for sparse numerics that has been increasingly explored as source of physical insight into many-body dynamics in recent years. In this work we revisit the venerable Anderson model of localization in dimensions $d=1, 2, 3, 4$ to construct local integrals of motion (LIOM) in Krylov space. These appear as zero eigenvalue edge states of an effective hopping problem in the Krylov superoperator subspace and can be analytically constructed given the Lanczos coefficients. We exploit this idea, focusing on $d=3$, to study the manifestation of the disorder driven Anderson transition in the anatomy of LIOMs. We find that the increasing complexity of the Krylov operators results in a suppression of the fluctuations of the Lanczos coefficients. As such, one can study the phenomenology of the integrals of motion in the disorder averaged Krylov chain. We find edge states localized on vanishing fraction of Krylov space (of dimension $D_K=V^2$ for cubes of volume $V$), both in localized and extended phases. Importantly, in the localized phase, disorder induces powerlaw decaying dimerization in the (Krylov) hopping problem, producing stretched exponential decay of the LIOMs in Krylov space with a stretching exponent $1/2d$. Metallic LIOMs are completely delocalized albeit across only $\propto \sqrt{D_K}$ states. Critical LIOMs exhibit powerlaw decay with an exponent matching the expected value of $0.29$.

DOI
10.1103/y9y9-gfwc

Related papers

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

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

Anderson Localization: A Floquet operator Krylov space perspective

Hsiu-Chung Yeh, Aditi Mitra

The problem of Anderson localization, as well as the single particle localization-delocalizaton transition of the Aubry-André model, is studied employing operator Krylov space methods. It is shown that even when the dynamics is generated by a Hamiltonian, studying the dynamics at stroboscopic rather than continuous times has its advantages. In particular, mapping the dynamics to an effective Floquet problem results in an operator Krylov space description where quantities such as the spectral function can be computed with fewer computational resources, while a moment method exists that allows for the extraction of Krylov parameters directly from the discrete time autocorrelation function. For stroboscopic dynamics, the operator Krylov space corresponds to the dynamics of an edge operator of an inhomogeneous Floquet transverse field Ising model, with the parameters of this effective model generated recursively. The Krylov parameters show disorder-averaged renormalization with their distribution narrowing as the recursion step increases. It is shown that a more physical spectral function is obtained from the Krylov parameters obtained from the disorder-averaged autocorrelation function, rather than the disorder-averaged Krylov parameters. The delocalized (localized) phase is shown to correspond to the appearance (absence) of a Porter-Thomas distribution, a ballistically propagating (localized) wavefront in operator Krylov space, and a smooth (discrete) Berstein-Szegö power-spectrum. The localization-delocalization transition is also demonstrated in operator Krylov space. A Porter-Thomas distribution is also observed at the critical point. The long-time dynamics and the inverse participation ratio at the critical point is shown to exhibit behavior consistent with a multi-fractal scaling with system size.

cond-mat.dis-nnCross-listedarXiv:2605.24115•May 26, 2026•0 citations

Krylov localization and suppression of complexity

Eliezer Rabinovici, Adrián Sánchez-Garrido, Ruth Shir +1 more

A bstract Quantum complexity, suitably defined, has been suggested as an important probe of late-time dynamics of black holes, particularly in the context of AdS/CFT. A notion of quantum complexity can be effectively captured by quantifying the spread of an operator in Krylov space as a consequence of time evolution. Complexity is expected to behave differently in chaotic many-body systems, as compared to integrable ones. In this paper we investigate Krylov complexity for the case of interacting integrable models at finite size and find that complexity saturation is suppressed as compared to chaotic systems. We associate this behavior with a novel localization phenomenon on the Krylov chain by mapping the theory of complexity growth and spread to an Anderson localization hopping model with off-diagonal disorder, and find that localization is enhanced in the integrable case due to a stronger disorder in the hopping amplitudes, inducing an effective suppression of Krylov complexity. We demonstrate this behavior for an interacting integrable model, the XXZ spin chain, and show that the same behavior results from a phenomenological model that we define: this model captures the essential features of our analysis and is able to reproduce the behaviors we observe for chaotic and integrable systems via an adjustable disorder parameter.

arXiv:2112.12128•Journal of High Energy Physics•Mar 1, 2022•123 citations

Krylov delocalization/localization across ergodicity breaking

Heiko Georg Menzler, Rishabh Jha

Krylov complexity has recently gained attention where the growth of operator complexity in time is measured in terms of the off-diagonal operator Lanczos coefficients. The operator Lanczos algorithm reduces the problem of complexity growth to a single-particle semi-infinite tight-binding chain (known as the Krylov chain). Employing the phenomenon of Anderson localization, we propose the phenomenology of inverse localization length on the Krylov chain that undergoes delocalization/localization transition on the Krylov chain while the physical system undergoes ergodicity breaking. On the Krylov chain we find delocalization in an ergodic regime, as we show for the SYK model, and localization in case of a weakly ergodicity-broken regime. Considering the dynamics beyond scrambling, we find a collapse across different operators in the ergodic regime. We test for two settings: (1) the coupled SYK model, and (2) the quantum East model. Our findings open avenues for mapping ergodicity/weak ergodicity-breaking transitions to delocalization/localization phenomenology on the Krylov chain.

arXiv:2403.14384•Physical review. B./Physical review. B•Sep 17, 2024•25 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
The Krylov subspace expansion is a workhorse method for sparse numerics that has been increasingly explored as source of physical insight into many-body dynamics in recent years. In this work we revisit the venerable Anderson model of localization in dimensions d=1,2,3,4d=1, 2, 3, 4d=1,2,3,4 to construct local integrals of motion (LIOM) in Krylov space. These appear as zero eigenvalue edge states of an effective hopping problem in the Krylov superoperator subspace and can be analytically constructed given the Lanczos coefficients. We exploit this idea, focusing on d=3d=3d=3, to study the manifestation of the disorder driven Anderson transition in the anatomy of LIOMs. We find that the increasing complexity of the Krylov operators results in a suppression of the fluctuations of the Lanczos coefficients. As such, one can study the phenomenology of the integrals of motion in the disorder averaged Krylov chain. We find edge states localized on vanishing fraction of Krylov space (of dimension DK=V2D_K=V^2DK​=V2 for cubes of volume VVV), both in localized and extended phases. Importantly, in the localized phase, disorder induces powerlaw decaying dimerization in the (Krylov) hopping problem, producing stretched exponential decay of the LIOMs in Krylov space with a stretching exponent 1/2d1/2d1/2d. Metallic LIOMs are completely delocalized albeit across only ∝DK\propto \sqrt{D_K}∝DK​​ states. Critical LIOMs exhibit powerlaw decay with an exponent matching the expected value of 0.290.290.29.
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.
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.
The problem of Anderson localization, as well as the single particle localization-delocalizaton transition of the Aubry-André model, is studied employing operator Krylov space methods. It is shown that even when the dynamics is generated by a Hamiltonian, studying the dynamics at stroboscopic rather than continuous times has its advantages. In particular, mapping the dynamics to an effective Floquet problem results in an operator Krylov space description where quantities such as the spectral function can be computed with fewer computational resources, while a moment method exists that allows for the extraction of Krylov parameters directly from the discrete time autocorrelation function. For stroboscopic dynamics, the operator Krylov space corresponds to the dynamics of an edge operator of an inhomogeneous Floquet transverse field Ising model, with the parameters of this effective model generated recursively. The Krylov parameters show disorder-averaged renormalization with their distribution narrowing as the recursion step increases. It is shown that a more physical spectral function is obtained from the Krylov parameters obtained from the disorder-averaged autocorrelation function, rather than the disorder-averaged Krylov parameters. The delocalized (localized) phase is shown to correspond to the appearance (absence) of a Porter-Thomas distribution, a ballistically propagating (localized) wavefront in operator Krylov space, and a smooth (discrete) Berstein-Szegö power-spectrum. The localization-delocalization transition is also demonstrated in operator Krylov space. A Porter-Thomas distribution is also observed at the critical point. The long-time dynamics and the inverse participation ratio at the critical point is shown to exhibit behavior consistent with a multi-fractal scaling with system size.
A bstract Quantum complexity, suitably defined, has been suggested as an important probe of late-time dynamics of black holes, particularly in the context of AdS/CFT. A notion of quantum complexity can be effectively captured by quantifying the spread of an operator in Krylov space as a consequence of time evolution. Complexity is expected to behave differently in chaotic many-body systems, as compared to integrable ones. In this paper we investigate Krylov complexity for the case of interacting integrable models at finite size and find that complexity saturation is suppressed as compared to chaotic systems. We associate this behavior with a novel localization phenomenon on the Krylov chain by mapping the theory of complexity growth and spread to an Anderson localization hopping model with off-diagonal disorder, and find that localization is enhanced in the integrable case due to a stronger disorder in the hopping amplitudes, inducing an effective suppression of Krylov complexity. We demonstrate this behavior for an interacting integrable model, the XXZ spin chain, and show that the same behavior results from a phenomenological model that we define: this model captures the essential features of our analysis and is able to reproduce the behaviors we observe for chaotic and integrable systems via an adjustable disorder parameter.
Krylov complexity has recently gained attention where the growth of operator complexity in time is measured in terms of the off-diagonal operator Lanczos coefficients. The operator Lanczos algorithm reduces the problem of complexity growth to a single-particle semi-infinite tight-binding chain (known as the Krylov chain). Employing the phenomenon of Anderson localization, we propose the phenomenology of inverse localization length on the Krylov chain that undergoes delocalization/localization transition on the Krylov chain while the physical system undergoes ergodicity breaking. On the Krylov chain we find delocalization in an ergodic regime, as we show for the SYK model, and localization in case of a weakly ergodicity-broken regime. Considering the dynamics beyond scrambling, we find a collapse across different operators in the ergodic regime. We test for two settings: (1) the coupled SYK model, and (2) the quantum East model. Our findings open avenues for mapping ergodicity/weak ergodicity-breaking transitions to delocalization/localization phenomenology on the Krylov chain.
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.