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

Krylov Break Times from an Inhomogeneous Lieb--Robinson Light Cone

Shunji Matsuura, Yoji Kawamura, Joseph Salfi, Satoshi Iso

quant-phcond-mat.dis-nnhep-thSubmitted Aug 27, 2026• Announced Aug 28, 2026• 0 citations
arXivHTMLPDF

Abstract

Krylov and Lanczos approximations are used in quantum dynamics, quantum subspace methods, and Hamiltonian learning. A practical question is how long an $m$-dimensional Krylov truncation can be trusted. We argue that this time is fixed by causal propagation on the associated Jacobi chain. The relevant distance is not the Krylov index itself, but the inhomogeneous transport metric $ρ(m,n) = \sum_{j=\min(m,n)}^{\max(m,n)-1} 1/b_j$, where $b_j$ is the Lanczos hopping across the bond $j \leftrightarrow j+1$. We prove a Lieb--Robinson bound in this metric. Its small-weight limit gives the velocity $v_{\rm LR} = 2$, meaning that propagation is exponentially suppressed outside the cone $ρ(m,n) \simeq 2|t|$. The error of a finite Krylov approximation to the return amplitude is a round-trip effect: information has to travel from the probe to the truncation boundary and back. Combining the Lieb--Robinson bound with Duhamel's formula yields a lower bound on the error of the truncated dynamics. For a fixed tolerance $ε$, let the break time $t_\ast(m;ε)$ denote the longest time for which the $m$-dimensional truncation is guaranteed to reproduce the exact return amplitude within error $ε$. We show that $t_\ast(m;ε) \ge τ_m[1-o(1)]$, where $τ_m = ρ(0,m) = \sum_{j<m} 1/b_j$. When the probe spreads along the chain, this lower bound is also tight, so $t_\ast(m) \simeq τ_m$. The situation is different when the probe excites only a localized part of the spectrum, or a part already resolved by the truncation. In this case, essentially no signal reaches the boundary. Beyond a state-dependent Krylov dimension $m_\ast$, the approximation can therefore remain accurate at all times, and the break time is effectively infinite. Numerical tests on spin chains and random Jacobi matrices support $t_\ast(m) \simeq τ_m$ in the transport-limited regime.

Comments
23 pages, 7 figures

Related papers

On Information Propagation in Krylov Subspaces

Pablo Martinez-Azcona, Maximilian Meyer-Moelleringhof, Adolfo del Campo +1 more

Krylov subspace methods provide an efficient tool for modeling any quantum initial value problem as a one dimensional chain. A time ordered path formulation, in which the dynamical generator successively splits probability between paths, yields a speed limit within such dynamical subspaces. For power law growing Lanczos coefficients $b_n \sim n^δ$ we find super ballistic lightcones, with an exponentially spreading lightcone for maximal growth $δ= 1$. Causality in the Krylov chain constrains the Lanczos coefficients to grow at most linearly. We illustrate the validity of the above in the $SL(2, R)$ and Heisenberg Weyl algebras, and in a many body skew field Ising model. Moreover, we extend the above bound to the family of $Q$ complexities. Lastly, mapping to a spin chain gives a Lieb Robinson bound analogue and opens the problem to standard condensed matter methods.

quant-pharXiv:2610.03883•Oct 6, 2026•0 citations

Tightening the Lieb-Robinson Bound in Locally Interacting Systems

Zhiyuan Wang, Kaden R.A. Hazzard

The Lieb-Robinson (LR) bound rigorously shows that in quantum systems with short-range interactions, the maximum amount of information that travels beyond an effective "light cone" decays exponentially with distance from the light-cone front, which expands at finite velocity.Despite being a fundamental result, existing bounds are often extremely loose, limiting their applications.We introduce a method that dramatically and qualitatively improves LR bounds in models with finite-range interactions.Most prominently, in systems with a large local Hilbert space dimension D, our method gives a LR velocity that grows much slower than previous bounds with D as D → ∞.For example, in the Heisenberg model with spin S, we find v ≤ const.compared to the previous v ∝ S, which diverges at large S, and in multiorbital Hubbard models with N orbitals, we find v ∝ √ N instead of previous v ∝ N , and similarly in the Nstate truncated Bose-Hubbard model and Wen's quantum rotor model.Our bounds also scale qualitatively better in some systems when the spatial dimension or certain model parameters become large, for example in the d-dimensional quantum Ising model and perturbed toric code models.Even in spin-1/2 Ising and Fermi-Hubbard models, our method improves the LR velocity by an order of magnitude with typical model parameters, and significantly improves the LR bound at large distance and early time.

arXiv:1908.03997•PRX Quantum•Sep 3, 2020•62 citations

Lieb-Robinson bounds and out-of-time order correlators in a long-range spin chain

Luis Colmenárez, David J. Luitz

Lieb-Robinson bounds quantify the maximal speed of information spreading in nonrelativistic quantum systems.We discuss the relation of Lieb-Robinson bounds to out-of-time order correlators, which correspond to different norms of commutators C(r, t ) = [A i (t ), B i+r ] of local operators.Using an exact Krylov space-time evolution technique, we calculate these two different norms of such commutators for the spin-1/2 Heisenberg chain with interactions decaying as a power law 1/r α with distance r.Our numerical analysis shows that both norms (operator norm and normalized Frobenius norm) exhibit the same asymptotic behavior, namely, a linear growth in time at short times and a power-law decay in space at long distance, leading asymptotically to power-law light cones for α < 1 and to linear light cones for α > 1.The asymptotic form of the tails of C(r, t ) ∝ t/r α is described by short-time perturbation theory, which is valid at short times and long distances.

arXiv:2005.10257•Physical Review Research•Oct 8, 2020•43 citations

A Converse to Lieb–Robinson Bounds in One Dimension Using Index Theory

Daniel Ranard, Michael Walter, Freek Witteveen

Abstract Unitary dynamics with a strict causal cone (or “light cone”) have been studied extensively, under the name of quantum cellular automata (QCAs). In particular, QCAs in one dimension have been completely classified by an index theory. Physical systems often exhibit only approximate causal cones; Hamiltonian evolutions on the lattice satisfy Lieb–Robinson bounds rather than strict locality. This motivates us to study approximately locality preserving unitaries (ALPUs). We show that the index theory is robust and completely extends to one-dimensional ALPUs. As a consequence, we achieve a converse to the Lieb–Robinson bounds: any ALPU of index zero can be exactly generated by some time-dependent, quasi-local Hamiltonian in constant time. For the special case of finite chains with open boundaries, any unitary satisfying the Lieb–Robinson bound may be generated by such a Hamiltonian. We also discuss some results on the stability of operator algebras which may be of independent interest.

arXiv:2012.00741•Annales Henri Poincaré•Jul 26, 2022•24 citations

Lieb-Robinson Light Cone for Power-Law Interactions

Minh C. Tran, Andrew Y. Guo, Christopher L. Baldwin +3 more

The Lieb-Robinson theorem states that information propagates with a finite velocity in quantum systems on a lattice with nearest-neighbor interactions. What are the speed limits on information propagation in quantum systems with power-law interactions, which decay as 1/r^{α} at distance r? Here, we present a definitive answer to this question for all exponents α>2d and all spatial dimensions d. Schematically, information takes time at least r^{min{1,α-2d}} to propagate a distance r. As recent state transfer protocols saturate this bound, our work closes a decades-long hunt for optimal Lieb-Robinson bounds on quantum information dynamics with power-law interactions.

arXiv:2103.15828•Physical Review Letters•Oct 11, 2021•88 citations

Lieb--Robinson Bounds and Entanglement Limits in Quantum Dynamics with Finite-Dimensional Memory

Jeet Mozumdar

We derive an extension of the Lieb--Robinson bound that accounts for memory effects in quantum lattice systems evolving under discrete-time non-Markovian dynamics generated by local interactions with finite-dimensional memory subsystems associated with each bond. The global evolution is unitary on an extended Hilbert space \[ \mathcal{H}_{\mathrm{ext}} = \bigotimes_{i=1}^{N} \mathcal{H}_i \otimes \bigotimes_{\text{bonds}} \mathcal{H}_M^{(i)}, \] while non-Markovian behavior, arising from repeated coupling to the same bond memory at successive time steps, leads to temporal correlations in the system. In this work, “finite memory” refers specifically to memory subsystems of fixed finite Hilbert space dimension whose reduced dynamics satisfy an exponential mixing condition in operator norm, meaning that correlation functions decay exponentially over time, with a characteristic timescale \(\tau\). More precisely, the memory channels satisfy \[ \left\| \mathcal{T}_M^s(X) - \mathrm{Tr}(X)\,\sigma_M \right\| \le c\, e^{-s/\tau}\, \|X\|, \] for all memory observables \(X\), where \(\sigma_M\) is a stationary state, and \(c>0\) depends only on the memory dynamics. For a one-dimensional lattice with nearest-neighbor interactions, we prove that the connected correlation function, defined as \[ C_{i,j}(t) := \sup_{\substack{A_i = A_i^\dagger,\; B_j = B_j^\dagger \\ \|A_i\|\le 1,\; \|B_j\|\le 1}} \left| \mathrm{Tr}\!\left[(A_i \otimes B_j)\rho(t)\right] - \mathrm{Tr}\!\left[A_i \rho_i(t)\right]\, \mathrm{Tr}\!\left[B_j \rho_j(t)\right] \right|, \] satisfies an exponential bound of the form \[ C_{i,j}(t) \le C_0 \exp\!\Big[-\frac{d(i,j) - v_{\mathrm{eff}}(\tau)\, t}{\xi}\Big], \] where \(d(i,j)=|i-j|\) and \(\xi>0\) is independent of system size and time. The effective velocity \(v_{\mathrm{eff}}(\tau)\) is explicitly bounded in terms of the local interaction strength and the memory mixing parameters, ensuring the system exhibits a finite light cone, meaning that information propagation is confined to a region growing linearly with time, with a well-defined maximum speed. In the limit \(\tau \to \infty\), the bound no longer guarantees a finite propagation speed, reflecting the transition to a persistent memory regime where correlations do not decay exponentially. As an operational consequence, the entanglement fidelity between two-site reduced states \(\rho_{i,j}(t)\) and a maximally entangled state satisfies \[ F_{i,j}(t) = \langle \Phi^+ | \rho_{i,j}(t) | \Phi^+ \rangle \le \tfrac{1}{2}\bigl[1 + C_{i,j}(t)\bigr]. \] These results extend Lieb--Robinson theory to a class of discrete-time quantum dynamics with memory subsystems associated with each bond, enabling resummation of multi-step memory effects into an effective correlation propagation kernel that governs the decay of temporal correlations. This establishes a quantitative relationship between the timescale of memory mixing and the effective speed of information propagation in the system.

•Apr 26, 2026•0 citations
Krylov and Lanczos approximations are used in quantum dynamics, quantum subspace methods, and Hamiltonian learning. A practical question is how long an mmm-dimensional Krylov truncation can be trusted. We argue that this time is fixed by causal propagation on the associated Jacobi chain. The relevant distance is not the Krylov index itself, but the inhomogeneous transport metric ρ(m,n)=∑j=min⁡(m,n)max⁡(m,n)−11/bjρ(m,n) = \sum_{j=\min(m,n)}^{\max(m,n)-1} 1/b_jρ(m,n)=∑j=min(m,n)max(m,n)−1​1/bj​, where bjb_jbj​ is the Lanczos hopping across the bond j↔j+1j \leftrightarrow j+1j↔j+1. We prove a Lieb--Robinson bound in this metric. Its small-weight limit gives the velocity vLR=2v_{\rm LR} = 2vLR​=2, meaning that propagation is exponentially suppressed outside the cone ρ(m,n)≃2∣t∣ρ(m,n) \simeq 2|t|ρ(m,n)≃2∣t∣. The error of a finite Krylov approximation to the return amplitude is a round-trip effect: information has to travel from the probe to the truncation boundary and back. Combining the Lieb--Robinson bound with Duhamel's formula yields a lower bound on the error of the truncated dynamics. For a fixed tolerance εεε, let the break time t∗(m;ε)t_\ast(m;ε)t∗​(m;ε) denote the longest time for which the mmm-dimensional truncation is guaranteed to reproduce the exact return amplitude within error εεε. We show that t∗(m;ε)≥τm[1−o(1)]t_\ast(m;ε) \ge τ_m[1-o(1)]t∗​(m;ε)≥τm​[1−o(1)], where τm=ρ(0,m)=∑j<m1/bjτ_m = ρ(0,m) = \sum_{j<m} 1/b_jτm​=ρ(0,m)=∑j<m​1/bj​. When the probe spreads along the chain, this lower bound is also tight, so t∗(m)≃τmt_\ast(m) \simeq τ_mt∗​(m)≃τm​. The situation is different when the probe excites only a localized part of the spectrum, or a part already resolved by the truncation. In this case, essentially no signal reaches the boundary. Beyond a state-dependent Krylov dimension m∗m_\astm∗​, the approximation can therefore remain accurate at all times, and the break time is effectively infinite. Numerical tests on spin chains and random Jacobi matrices support t∗(m)≃τmt_\ast(m) \simeq τ_mt∗​(m)≃τm​ in the transport-limited regime.
Krylov subspace methods provide an efficient tool for modeling any quantum initial value problem as a one dimensional chain. A time ordered path formulation, in which the dynamical generator successively splits probability between paths, yields a speed limit within such dynamical subspaces. For power law growing Lanczos coefficients bn∼nδb_n \sim n^δbn​∼nδ we find super ballistic lightcones, with an exponentially spreading lightcone for maximal growth δ=1δ= 1δ=1. Causality in the Krylov chain constrains the Lanczos coefficients to grow at most linearly. We illustrate the validity of the above in the SL(2,R)SL(2, R)SL(2,R) and Heisenberg Weyl algebras, and in a many body skew field Ising model. Moreover, we extend the above bound to the family of QQQ complexities. Lastly, mapping to a spin chain gives a Lieb Robinson bound analogue and opens the problem to standard condensed matter methods.
The Lieb-Robinson (LR) bound rigorously shows that in quantum systems with short-range interactions, the maximum amount of information that travels beyond an effective "light cone" decays exponentially with distance from the light-cone front, which expands at finite velocity.Despite being a fundamental result, existing bounds are often extremely loose, limiting their applications.We introduce a method that dramatically and qualitatively improves LR bounds in models with finite-range interactions.Most prominently, in systems with a large local Hilbert space dimension D, our method gives a LR velocity that grows much slower than previous bounds with D as D → ∞.For example, in the Heisenberg model with spin S, we find v ≤ const.compared to the previous v ∝ S, which diverges at large S, and in multiorbital Hubbard models with N orbitals, we find v ∝ √ N instead of previous v ∝ N , and similarly in the Nstate truncated Bose-Hubbard model and Wen's quantum rotor model.Our bounds also scale qualitatively better in some systems when the spatial dimension or certain model parameters become large, for example in the d-dimensional quantum Ising model and perturbed toric code models.Even in spin-1/2 Ising and Fermi-Hubbard models, our method improves the LR velocity by an order of magnitude with typical model parameters, and significantly improves the LR bound at large distance and early time.
Lieb-Robinson bounds quantify the maximal speed of information spreading in nonrelativistic quantum systems.We discuss the relation of Lieb-Robinson bounds to out-of-time order correlators, which correspond to different norms of commutators C(r, t ) = [A i (t ), B i+r ] of local operators.Using an exact Krylov space-time evolution technique, we calculate these two different norms of such commutators for the spin-1/2 Heisenberg chain with interactions decaying as a power law 1/r α with distance r.Our numerical analysis shows that both norms (operator norm and normalized Frobenius norm) exhibit the same asymptotic behavior, namely, a linear growth in time at short times and a power-law decay in space at long distance, leading asymptotically to power-law light cones for α < 1 and to linear light cones for α > 1.The asymptotic form of the tails of C(r, t ) ∝ t/r α is described by short-time perturbation theory, which is valid at short times and long distances.
Abstract Unitary dynamics with a strict causal cone (or “light cone”) have been studied extensively, under the name of quantum cellular automata (QCAs). In particular, QCAs in one dimension have been completely classified by an index theory. Physical systems often exhibit only approximate causal cones; Hamiltonian evolutions on the lattice satisfy Lieb–Robinson bounds rather than strict locality. This motivates us to study approximately locality preserving unitaries (ALPUs). We show that the index theory is robust and completely extends to one-dimensional ALPUs. As a consequence, we achieve a converse to the Lieb–Robinson bounds: any ALPU of index zero can be exactly generated by some time-dependent, quasi-local Hamiltonian in constant time. For the special case of finite chains with open boundaries, any unitary satisfying the Lieb–Robinson bound may be generated by such a Hamiltonian. We also discuss some results on the stability of operator algebras which may be of independent interest.
The Lieb-Robinson theorem states that information propagates with a finite velocity in quantum systems on a lattice with nearest-neighbor interactions. What are the speed limits on information propagation in quantum systems with power-law interactions, which decay as 1/r^{α} at distance r? Here, we present a definitive answer to this question for all exponents α>2d and all spatial dimensions d. Schematically, information takes time at least r^{min{1,α-2d}} to propagate a distance r. As recent state transfer protocols saturate this bound, our work closes a decades-long hunt for optimal Lieb-Robinson bounds on quantum information dynamics with power-law interactions.
We derive an extension of the Lieb--Robinson bound that accounts for memory effects in quantum lattice systems evolving under discrete-time non-Markovian dynamics generated by local interactions with finite-dimensional memory subsystems associated with each bond. The global evolution is unitary on an extended Hilbert space Hext=⨂i=1NHi⊗⨂bondsHM(i), \mathcal{H}_{\mathrm{ext}} = \bigotimes_{i=1}^{N} \mathcal{H}_i \otimes \bigotimes_{\text{bonds}} \mathcal{H}_M^{(i)}, Hext​=⨂i=1N​Hi​⊗⨂bonds​HM(i)​, while non-Markovian behavior, arising from repeated coupling to the same bond memory at successive time steps, leads to temporal correlations in the system. In this work, “finite memory” refers specifically to memory subsystems of fixed finite Hilbert space dimension whose reduced dynamics satisfy an exponential mixing condition in operator norm, meaning that correlation functions decay exponentially over time, with a characteristic timescale τ\tauτ. More precisely, the memory channels satisfy ∥TMs(X)−Tr(X) σM∥≤c e−s/τ ∥X∥, \left\| \mathcal{T}_M^s(X) - \mathrm{Tr}(X)\,\sigma_M \right\| \le c\, e^{-s/\tau}\, \|X\|, ∥TMs​(X)−Tr(X)σM​∥≤ce−s/τ∥X∥, for all memory observables XXX, where σM\sigma_MσM​ is a stationary state, and c>0c>0c>0 depends only on the memory dynamics. For a one-dimensional lattice with nearest-neighbor interactions, we prove that the connected correlation function, defined as Ci,j(t):=sup⁡Ai=Ai†,Bj=Bj†∥Ai∥≤1,∥Bj∥≤1∣Tr ⁣[(Ai⊗Bj)ρ(t)]−Tr ⁣[Aiρi(t)] Tr ⁣[Bjρj(t)]∣, C_{i,j}(t) := \sup_{\substack{A_i = A_i^\dagger,\; B_j = B_j^\dagger \\ \|A_i\|\le 1,\; \|B_j\|\le 1}} \left| \mathrm{Tr}\!\left[(A_i \otimes B_j)\rho(t)\right] - \mathrm{Tr}\!\left[A_i \rho_i(t)\right]\, \mathrm{Tr}\!\left[B_j \rho_j(t)\right] \right|, Ci,j​(t):=supAi​=Ai†​,Bj​=Bj†​∥Ai​∥≤1,∥Bj​∥≤1​​∣Tr[(Ai​⊗Bj​)ρ(t)]−Tr[Ai​ρi​(t)]Tr[Bj​ρj​(t)]∣, satisfies an exponential bound of the form Ci,j(t)≤C0exp⁡ ⁣[−d(i,j)−veff(τ) tξ], C_{i,j}(t) \le C_0 \exp\!\Big[-\frac{d(i,j) - v_{\mathrm{eff}}(\tau)\, t}{\xi}\Big], Ci,j​(t)≤C0​exp[−ξd(i,j)−veff​(τ)t​], where d(i,j)=∣i−j∣d(i,j)=|i-j|d(i,j)=∣i−j∣ and ξ>0\xi>0ξ>0 is independent of system size and time. The effective velocity veff(τ)v_{\mathrm{eff}}(\tau)veff​(τ) is explicitly bounded in terms of the local interaction strength and the memory mixing parameters, ensuring the system exhibits a finite light cone, meaning that information propagation is confined to a region growing linearly with time, with a well-defined maximum speed. In the limit τ→∞\tau \to \inftyτ→∞, the bound no longer guarantees a finite propagation speed, reflecting the transition to a persistent memory regime where correlations do not decay exponentially. As an operational consequence, the entanglement fidelity between two-site reduced states ρi,j(t)\rho_{i,j}(t)ρi,j​(t) and a maximally entangled state satisfies Fi,j(t)=⟨Φ+∣ρi,j(t)∣Φ+⟩≤12[1+Ci,j(t)]. F_{i,j}(t) = \langle \Phi^+ | \rho_{i,j}(t) | \Phi^+ \rangle \le \tfrac{1}{2}\bigl[1 + C_{i,j}(t)\bigr]. Fi,j​(t)=⟨Φ+∣ρi,j​(t)∣Φ+⟩≤21​[1+Ci,j​(t)]. These results extend Lieb--Robinson theory to a class of discrete-time quantum dynamics with memory subsystems associated with each bond, enabling resummation of multi-step memory effects into an effective correlation propagation kernel that governs the decay of temporal correlations. This establishes a quantitative relationship between the timescale of memory mixing and the effective speed of information propagation in the system.