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:2602.19750 · The European Physical Journal C

Krylov Distribution and Universal Convergence of Quantum Fisher Information

Mohsen Alishahiha, Fatemeh Tarighi Tabesh, Mohammad Javad Vasli

quant-phcond-mat.stat-mechhep-thSubmitted Feb 23, 2026• Announced Feb 24, 2026• 0 citations
v1 · Feb 23, 2026v2 · Apr 20, 2026v3 · Aug 19, 2026
arXivHTMLPDF

Abstract

We develop a spectral-resolvent framework for computing the quantum Fisher information (QFI) using Krylov subspace methods, extending the notion of the Krylov distribution. By expressing the QFI as a resolvent moment of the superoperator $\mathcal{K}_ρ$ associated with a density matrix, the Krylov distribution quantifies how the QFI weight is distributed across Krylov levels in operator space and provides a natural measure for controlling the truncation error in Krylov approximations. Leveraging orthogonal polynomial theory, we identify two universal convergence regimes: exponential decay when the Liouville-space spectrum is gapped away from zero, and algebraic decay governed by hard-edge (Bessel) universality when small eigenvalues accumulate near zero. This framework establishes a direct connection between quantum metrology, spectral geometry, and Krylov dynamics, offering both conceptual insight and practical tools for efficient QFI computation in high-dimensional and many-body systems.

Comments
12 pages, 2 figures

Related papers

Efficient Quantum Algorithm for Mixed-State Quantum Geometry

Chenyang Wang, Qi Zhao, Xiao Yuan +1 more

Mixed-state quantum geometry plays a fundamental role in quantum metrology, multiparameter estimation, and the characterization of quantum many-body systems. Yet, unlike its pure-state counterpart, accessing mixed-state geometry requires solving an ill-conditioned Lyapunov equation for the symmetric logarithmic derivative (SLD), posing a major challenge for resource-efficient quantum computation. Here we address the challenge and develop a novel quantum algorithmic framework for efficiently accessing the quantum Fisher information matrix (QFIM) and quantum geometric tensor (QGT) of mixed states. Our key insight is to uncover a previously unexplored connection between the positive semidefinite structure of density matrices and exponentially convergent Gauss-Laguerre quadrature, enabling an efficient integral-based construction of the SLD. A rigorous complexity analysis further establishes the resource efficiency of our algorithm, demonstrating logarithmic precision dependence and double-logarithmic ancillary-qubit overhead at fixed conditioning. Building on this construction, we achieve near-Heisenberg precision scaling for diagonal QFI estimation without quantum phase estimation and construct unified block encodings of the full QFIM and QGT under a bounded spectral-gap assumption. A complementary regularized algorithm extends efficient diagonal estimation to exponentially small spectral gaps under a bounded-response condition. Our results establish a resource-efficient route to coherent mixed-state quantum geometry, enabling its direct integration into quantum linear-algebra algorithms without classical reconstruction and opening new possibilities for quantum-enhanced multiparameter metrology.

quant-pharXiv:2610.24686•Oct 20, 2026•0 citations

Superiority of Krylov shadow tomography in estimating quantum Fisher information: From bounds to exactness

Yuan-Hao Wang, Da-Jian Zhang

Estimating the quantum Fisher information (QFI) is a crucial yet challenging task with widespread applications across quantum science and technologies. The recently proposed Krylov shadow tomography (KST) opens a new avenue for this task by introducing a series of Krylov bounds on the QFI. In this work, we address the practical applicability of the KST, unveiling that the Krylov bounds of low orders already enable efficient and accurate estimation of the QFI. We show that the Krylov bounds converge to the QFI exponentially fast with increasing order and can surpass the state-of-the-art polynomial lower bounds known to date. Moreover, we show that certain low-order Krylov bound can already match the QFI exactly for low-rank states prevalent in practical settings. Such exact match is beyond the reach of polynomial lower bounds proposed previously. These theoretical findings, solidified by extensive numerical simulations, demonstrate practical advantages over existing polynomial approaches, holding promise for fully unlocking the effectiveness of QFI-based applications.

quant-pharXiv:2602.17361•npj Quantum Information•Feb 20, 2026•1 citation

Quantum Fisher information and symmetric logarithmic derivative via anti-commutators

Jing Liu, Jie Chen, Xiao-Xing Jing +1 more

Abstract The symmetric logarithmic derivative (SLD) is a key quantity to obtain quantum Fisher information (QFI) and to construct the corresponding optimal measurements. Here we develop a method to calculate the SLD and QFI via anti-commutators. This method has originated from the Lyapunov representation and would be very useful for cases where the anti-commutators among the state and its partial derivative exhibit periodic properties. As an application, we discuss a class of states whose squares linearly depend on the states themselves, and give the corresponding analytical expressions of SLD and QFI. A noisy scenario of this class of states is also considered and discussed. Finally, we readily apply the method to the block-diagonal states and the multi-parameter estimation scenarios.

arXiv:1501.04290•Journal of Physics A Mathematical and Theoretical•May 24, 2016•65 citations

Krylov Distribution

Mohsen Alishahiha, Mohammad Javad Vasli

We introduce the Krylov distribution $\mathcal{D}(ξ)$, a static Krylov-space diagnostic that characterizes how inverse-energy response is organized in Hilbert space. The central object is the resolvent-dressed state $(H-ξ)^{-1}|ψ_0\rangle$, whose decomposition in the Krylov basis generated from a reference state defines a normalized distribution over Krylov levels. Unlike conventional spectral functions, which resolve response solely along the energy axis, the Krylov distribution captures how the resolvent explores the dynamically accessible subspace as the spectral parameter $ξ$ is varied. Using asymptotic analysis, exact results in solvable models, and numerical studies of an interacting spin chain, we identify three universal regimes: saturation outside the spectral support, extensive growth within continuous spectra, and sublinear or logarithmic scaling near spectral edges and quantum critical points. We further show that fidelity susceptibility and the quantum geometric tensor admit natural decompositions in terms of Krylov-resolved resolvent amplitudes.

hep-thCross-listedarXiv:2602.06150•Physical review. D/Physical review. D.•Feb 9, 2026•0 citations

Hall's Exact Uncertainty Relation and the Geometry of Krylov-Space Operator Dynamics

Mohsen Alishahiha, Souvik Banerjee

Krylov representations of operator dynamics encode two inequivalent motions: probability transport along the discrete Krylov chain and motion of the normalized amplitude vector on the associated sphere. We lift Hall's exact uncertainty relation to Liouville space and show that, for the fixed Krylov-basis measurement, its classical Fisher information locally attains the quantum Fisher information. The amplitude vector therefore follows an exact constant-speed spherical trajectory. We construct its moving frame, identifying the first Lanczos coefficient with the speed and the higher coefficients with successive bending. Geodesic curvature controls the leading departure from Mandelstam-Tamm saturation. Separately, we compare cumulative and pointwise bounds on Krylov-index transport, which depends on the full Lanczos sequence.

quant-pharXiv:2512.20359•Dec 24, 2025•0 citations

Quantum Fisher Information under decoherence with explicit wavefunctions

Francesco Musso, Vittorio Vitale, Sara Murciano

We present a method to estimate the quantum Fisher information (QFI) of many-body quantum states in the presence of decoherence, where its direct evaluation requires the full spectral resolution of the density matrix. We show that, for many-body wave functions known analytically in the occupation-number basis, systematic lower bounds to the QFI can be mapped onto expectation values over a classical probability distribution defined by the wave function amplitudes. This mapping enables efficient estimation via Markov-chain Monte Carlo sampling, with a computational cost that scales as a `slow' exponential ($e^{b L}$ with $b \lesssim 0.6$) and remains manageable for system sizes well beyond exact diagonalization. We specify this framework to Jastrow-Gutzwiller wave functions. We characterize their metrological content by identifying the observables that maximize the QFI and the corresponding scaling with $L$. Then, we analyze the QFI under three physically motivated noise channels: local dephasing, local amplitude damping, and global depolarizing. We compare polynomial and Krylov-based lower bounds across these channels, relating their behavior to the effective rank of the noisy density matrix and to the structure of the operator generating the parameter encoding. The framework extends naturally to other analytically known wave functions and to a broader class of information-theoretic quantities beyond the QFI.

quant-pharXiv:2605.22917•May 25, 2026•0 citations
We develop a spectral-resolvent framework for computing the quantum Fisher information (QFI) using Krylov subspace methods, extending the notion of the Krylov distribution. By expressing the QFI as a resolvent moment of the superoperator Kρ\mathcal{K}_ρKρ​ associated with a density matrix, the Krylov distribution quantifies how the QFI weight is distributed across Krylov levels in operator space and provides a natural measure for controlling the truncation error in Krylov approximations. Leveraging orthogonal polynomial theory, we identify two universal convergence regimes: exponential decay when the Liouville-space spectrum is gapped away from zero, and algebraic decay governed by hard-edge (Bessel) universality when small eigenvalues accumulate near zero. This framework establishes a direct connection between quantum metrology, spectral geometry, and Krylov dynamics, offering both conceptual insight and practical tools for efficient QFI computation in high-dimensional and many-body systems.
Mixed-state quantum geometry plays a fundamental role in quantum metrology, multiparameter estimation, and the characterization of quantum many-body systems. Yet, unlike its pure-state counterpart, accessing mixed-state geometry requires solving an ill-conditioned Lyapunov equation for the symmetric logarithmic derivative (SLD), posing a major challenge for resource-efficient quantum computation. Here we address the challenge and develop a novel quantum algorithmic framework for efficiently accessing the quantum Fisher information matrix (QFIM) and quantum geometric tensor (QGT) of mixed states. Our key insight is to uncover a previously unexplored connection between the positive semidefinite structure of density matrices and exponentially convergent Gauss-Laguerre quadrature, enabling an efficient integral-based construction of the SLD. A rigorous complexity analysis further establishes the resource efficiency of our algorithm, demonstrating logarithmic precision dependence and double-logarithmic ancillary-qubit overhead at fixed conditioning. Building on this construction, we achieve near-Heisenberg precision scaling for diagonal QFI estimation without quantum phase estimation and construct unified block encodings of the full QFIM and QGT under a bounded spectral-gap assumption. A complementary regularized algorithm extends efficient diagonal estimation to exponentially small spectral gaps under a bounded-response condition. Our results establish a resource-efficient route to coherent mixed-state quantum geometry, enabling its direct integration into quantum linear-algebra algorithms without classical reconstruction and opening new possibilities for quantum-enhanced multiparameter metrology.
Estimating the quantum Fisher information (QFI) is a crucial yet challenging task with widespread applications across quantum science and technologies. The recently proposed Krylov shadow tomography (KST) opens a new avenue for this task by introducing a series of Krylov bounds on the QFI. In this work, we address the practical applicability of the KST, unveiling that the Krylov bounds of low orders already enable efficient and accurate estimation of the QFI. We show that the Krylov bounds converge to the QFI exponentially fast with increasing order and can surpass the state-of-the-art polynomial lower bounds known to date. Moreover, we show that certain low-order Krylov bound can already match the QFI exactly for low-rank states prevalent in practical settings. Such exact match is beyond the reach of polynomial lower bounds proposed previously. These theoretical findings, solidified by extensive numerical simulations, demonstrate practical advantages over existing polynomial approaches, holding promise for fully unlocking the effectiveness of QFI-based applications.
Abstract The symmetric logarithmic derivative (SLD) is a key quantity to obtain quantum Fisher information (QFI) and to construct the corresponding optimal measurements. Here we develop a method to calculate the SLD and QFI via anti-commutators. This method has originated from the Lyapunov representation and would be very useful for cases where the anti-commutators among the state and its partial derivative exhibit periodic properties. As an application, we discuss a class of states whose squares linearly depend on the states themselves, and give the corresponding analytical expressions of SLD and QFI. A noisy scenario of this class of states is also considered and discussed. Finally, we readily apply the method to the block-diagonal states and the multi-parameter estimation scenarios.
We introduce the Krylov distribution D(ξ)\mathcal{D}(ξ)D(ξ), a static Krylov-space diagnostic that characterizes how inverse-energy response is organized in Hilbert space. The central object is the resolvent-dressed state (H−ξ)−1∣ψ0⟩(H-ξ)^{-1}|ψ_0\rangle(H−ξ)−1∣ψ0​⟩, whose decomposition in the Krylov basis generated from a reference state defines a normalized distribution over Krylov levels. Unlike conventional spectral functions, which resolve response solely along the energy axis, the Krylov distribution captures how the resolvent explores the dynamically accessible subspace as the spectral parameter ξξξ is varied. Using asymptotic analysis, exact results in solvable models, and numerical studies of an interacting spin chain, we identify three universal regimes: saturation outside the spectral support, extensive growth within continuous spectra, and sublinear or logarithmic scaling near spectral edges and quantum critical points. We further show that fidelity susceptibility and the quantum geometric tensor admit natural decompositions in terms of Krylov-resolved resolvent amplitudes.
Krylov representations of operator dynamics encode two inequivalent motions: probability transport along the discrete Krylov chain and motion of the normalized amplitude vector on the associated sphere. We lift Hall's exact uncertainty relation to Liouville space and show that, for the fixed Krylov-basis measurement, its classical Fisher information locally attains the quantum Fisher information. The amplitude vector therefore follows an exact constant-speed spherical trajectory. We construct its moving frame, identifying the first Lanczos coefficient with the speed and the higher coefficients with successive bending. Geodesic curvature controls the leading departure from Mandelstam-Tamm saturation. Separately, we compare cumulative and pointwise bounds on Krylov-index transport, which depends on the full Lanczos sequence.
We present a method to estimate the quantum Fisher information (QFI) of many-body quantum states in the presence of decoherence, where its direct evaluation requires the full spectral resolution of the density matrix. We show that, for many-body wave functions known analytically in the occupation-number basis, systematic lower bounds to the QFI can be mapped onto expectation values over a classical probability distribution defined by the wave function amplitudes. This mapping enables efficient estimation via Markov-chain Monte Carlo sampling, with a computational cost that scales as a `slow' exponential (ebLe^{b L}ebL with b≲0.6b \lesssim 0.6b≲0.6) and remains manageable for system sizes well beyond exact diagonalization. We specify this framework to Jastrow-Gutzwiller wave functions. We characterize their metrological content by identifying the observables that maximize the QFI and the corresponding scaling with LLL. Then, we analyze the QFI under three physically motivated noise channels: local dephasing, local amplitude damping, and global depolarizing. We compare polynomial and Krylov-based lower bounds across these channels, relating their behavior to the effective rank of the noisy density matrix and to the structure of the operator generating the parameter encoding. The framework extends naturally to other analytically known wave functions and to a broader class of information-theoretic quantities beyond the QFI.