---
title: Krylov Shadow Tomography for QFI Estimation
url: https://www.emergentmind.com/topics/krylov-shadow-tomography-kst
type: topic
---

# Krylov Shadow Tomography for QFI Estimation

Searching arXiv for papers on Krylov shadow tomography and QFI estimation.
Krylov shadow tomography (KST) is a hybrid framework for estimating the quantum Fisher information (QFI) of an unknown quantum state by combining the Krylov subspace method with shadow tomography. It was introduced as a way to make QFI estimation experimentally feasible for large quantum systems, where direct evaluation is difficult because the QFI depends on the density operator in a highly nonlinear manner [2503.01697]. Subsequent work established that KST yields a strict hierarchy of non-polynomial lower bounds on the QFI, that these bounds converge exponentially fast with Krylov order, that they are tighter than state-of-the-art polynomial lower bounds at comparable resource level, and that for certain low-rank states the hierarchy terminates at the exact QFI already at low order [2602.17361].

## 1. QFI estimation problem and the motivation for KST

For a quantum state $\rho=\sum_k p_k |k\rangle\langle k|$ and Hermitian generator $H$, the QFI is given by
$$
F_Q = 2\sum_{p_k+p_l>0}\frac{(p_k-p_l)^2}{p_k+p_l}\, |\langle k|H|l\rangle|^2.
$$
This quantity is central in quantum metrology because it determines the ultimate precision limit in parameter estimation, and it is also used in entanglement detection and related tasks [2503.01697]. The difficulty is that $F_Q$ is highly nonlinear in $\rho$, so it is not directly accessible through the low-order polynomial observables naturally handled by modern randomized-measurement methods.

Earlier work therefore concentrated on experimentally accessible lower bounds on the QFI, including Legendre-transform bounds, sub-QFI bounds, and Taylor-series-based bounds [2602.17361]. These polynomial lower bounds are measurable, but they exhibit an unavoidable systematic gap: no polynomial lower bound can coincide with the QFI for all states [2602.17361]. Previous asymptotic polynomial hierarchies can approach $F_Q$, but only at high order, with rapidly increasing measurement and postprocessing cost [2503.01697].

KST addresses this limitation by replacing the polynomial-approximation viewpoint with a hierarchy of non-polynomial Krylov bounds. This shifts the task from approximating the QFI by polynomial surrogates to approximating an operator associated with the QFI inside a structured Krylov subspace. A plausible implication is that the central advantage of KST is not merely tighter lower bounds, but a qualitatively different approximation principle.

## 2. Operator-theoretic construction

The mathematical setup begins with the weighted inner product
$$
\langle X,Y\rangle_\rho = \operatorname{tr}\!\left[\rho(XY+YX)/2\right],
$$
and the associated norm $\|X\|_\rho=\sqrt{\langle X,X\rangle_\rho}$. The relevant superoperator is the symmetrization map
$$
\mathcal{R}_\rho(X)=\frac{1}{2}(\rho X+X\rho).
$$
In the eigenbasis of $\rho$, this acts as
$$
\mathcal{R}_\rho(X)=\sum_{p_k+p_l>0}\frac{p_k+p_l}{2}X_{kl}|k\rangle\langle l|,
$$
with inverse on the relevant operator subspace
$$
\mathcal{R}_\rho^{-1}(X)=\sum_{p_k+p_l>0}\frac{2}{p_k+p_l}X_{kl}|k\rangle\langle l|.
$$
Defining
$$
L=\mathcal{R}_\rho^{-1}\!\left(i[\rho,H]\right),
$$
one obtains the compact representation
$$
F_Q=\|L\|_\rho^2
$$
[2503.01697].

KST then constructs the nested Krylov subspaces
$$
\mathcal{K}_n = \operatorname{span}\left\{i[\rho,H],\, \mathcal{R}_\rho(i[\rho,H]),\, \ldots,\, \mathcal{R}_\rho^{\,n-1}(i[\rho,H])\right\}.
$$
Within each $\mathcal{K}_n$, one chooses the best approximation $L_n$ to $L$ in the $\rho$-weighted norm, and defines the corresponding Krylov bound
$$
B_n^{(\mathsf{Kry})}=\|L_n\|_\rho^2.
$$
Equivalently,
$$
\|L-L_n\|_\rho=\min_{X\in\mathcal{K}_n}\|L-X\|_\rho.
$$
Because $L_n$ is the orthogonal projection of $L$ onto $\mathcal{K}_n$, the exact decomposition
$$
F_Q=B_n^{(\mathsf{Kry})}+\|L-L_n\|_\rho^2
$$
holds, and therefore
$$
B_n^{(\mathsf{Kry})}\le F_Q
$$
[2503.01697].

The operator space is restricted to
$$
\mathcal{X}=\{X=X^\dagger: X_{kl}=0 \text{ if } p_k=p_l=0\},
$$
so that the weighted inner product is positive definite even when $\rho$ is singular. The nontrivial case assumes $[\rho,H]\neq 0$; otherwise $F_Q=0$ [2503.01697].

## 3. Hierarchy, strict improvement, and finite termination

A defining feature of KST is that it yields not a single bound but a strict hierarchy of increasingly tight lower bounds. The Krylov subspaces satisfy
$$
\mathcal{K}_1\subsetneq\mathcal{K}_2\subsetneq\cdots\subsetneq\mathcal{K}_{n^*}=\mathcal{K}_{n^*+j},
$$
for some finite stopping index $n^*$, and the exact operator $L$ belongs to $\mathcal{K}_{n^*}$ but not to any smaller $\mathcal{K}_n$ with $n<n^*$ [2503.01697]. Consequently,
$$
B_1^{(\mathsf{Kry})}<B_2^{(\mathsf{Kry})}<\cdots<B_{n^*}^{(\mathsf{Kry})}=F_Q.
$$
This strict hierarchy is one of the principal structural results of the framework [2602.17361].

The highest Krylov bound is exact because once the Krylov subspace contains $L$, the projection error vanishes. In this sense, KST is not limited to asymptotic approximation: at finite Krylov dimension it can reproduce the exact QFI [2503.01697]. This directly differentiates it from previously studied polynomial lower-bound schemes, which cannot coincide with the QFI for all states [2602.17361].

A common misconception is to regard KST as simply another polynomial lower-bound construction. That characterization is inaccurate. The bound $B_n^{(\mathsf{Kry})}$ is non-polynomial in $\rho$, even though it is assembled from polynomial moments. The distinction is essential because the exactness and finite termination properties of KST lie beyond the reach of polynomial lower bounds in general [2503.01697].

## 4. Convergence properties and comparison with polynomial bounds

The 2026 analysis established that low-order Krylov bounds are already practically effective. Its central quantitative result is an exponential convergence theorem for the relative error
$$
\mathcal{E}_n^{(\mathsf{Kry})}\coloneqq \frac{|B_n^{(\mathsf{Kry})}-F_Q|}{F_Q},
$$
namely
$$
\mathcal{E}_n^{(\mathsf{Kry})}\le 4\left[\frac{\sqrt{\kappa(\rho)}-1}{\sqrt{\kappa(\rho)}+1}\right]^{2n}.
$$
Here $\kappa(\rho)$ is a condition-number-like quantity determined by the spectrum of $\rho$; for full rank it is essentially the ratio of largest to smallest eigenvalue, with an appropriate modification for non-full-rank states [2602.17361]. The theorem implies exponential convergence of the Krylov hierarchy to the exact QFI as the Krylov order increases, with faster convergence for more mixed states.

The same work proved a direct superiority statement against the Taylor hierarchy $B_n^{(\mathsf{Tay})}$:
$$
\mathcal{E}_n^{(\mathsf{Kry})}\le \mathcal{E}_{2n-1}^{(\mathsf{Tay})}.
$$
This comparison is described as resource-fair because $B_n^{(\mathsf{Kry})}$ and $B_{2n-1}^{(\mathsf{Tay})}$ require estimating polynomial functions of $\rho$ of comparable degree [2602.17361]. The result shows that, at essentially the same estimation cost, the Krylov bound is generally tighter than the corresponding Taylor bound.

The significance of this comparison is methodological as well as numerical. Polynomial lower bounds approximate the nonlinear QFI through polynomial surrogates and therefore retain a systematic bias. KST instead constructs a hierarchy of non-polynomial approximants through projection in Krylov space. This suggests that the observed superiority is structurally rooted in the approximation scheme rather than in a favorable choice of constants.

## 5. Exactness for low-rank states

KST has a particularly sharp exactness theory for low-rank states. Let
$$
S=\{p_i+p_j: p_i\neq p_j\}
$$
be the set of sums of distinct eigenvalue pairs. The stopping order satisfies
$$
n^*\le \mathcal{N}[S],
$$
where $\mathcal{N}[S]$ denotes the number of distinct values in $S$ [2602.17361]. Since $\mathcal{N}[S]\le r(r+1)/2$ for a rank-$r$ state,
$$
n^*\le \frac{r(r+1)}{2}.
$$
Thus low-rank structure can force exact recovery of the QFI at very small Krylov order.

The Methods section further refines this by giving the exact formula
$$
n^*=\mathcal{N}[S]-\mathcal{N}[J],
$$
where $J\subseteq S$ consists of those eigenvalue sums $\mu$ for which all relevant matrix elements of $H$ vanish between the corresponding eigenspaces [2602.17361]. This characterizes when the Krylov sequence terminates earlier than the generic rank-based upper bound would suggest.

The original KST paper presented explicit examples with $n^*=1$. For the pseudo-pure state
$$
\rho_p=(1-p)|\psi\rangle\langle\psi|+p\frac{\openone}{2^N},
$$
it showed that $n^*=1$ for any $p$, any pure state $|\psi\rangle$, and any $H$, so that
$$
B_1^{(\mathsf{Kry})}=F_Q.
$$
For the bound entangled state
$$
\rho_k=\lambda P_k^++\frac{\lambda}{2}(Q_k^++Q_k^-),
$$
with $H=\frac12\sum_i \sigma_z^{(i)}$, it likewise proved $n^*=1$, hence again $B_1^{(\mathsf{Kry})}=F_Q$ [2503.01697].

The broader relevance of these results follows from the observation that low-rank and nearly low-rank states are common in quantum information processing, for example when ideal pure states are weakly corrupted by noise [2602.17361]. In such settings, low-order Krylov bounds can eliminate the systematic error that affects polynomial approximations.

## 6. Estimation by classical shadows, practical behavior, and limitations

KST becomes experimentally accessible because each bound can be expressed through quantities
$$
T_k=\operatorname{tr}\!\left(i[\rho,H]\mathcal{R}_\rho^k(i[\rho,H])\right),
$$
which are polynomial functionals of $\rho$ [2602.17361]. Writing
$$
A_{ij}=T_{i+j-1},\qquad b_i=T_{i-1},
$$
the bound is computed as
$$
B_n^{(\mathsf{Kry})}=b^T A^{-1} b.
$$
Equivalently, the original derivation formulates $A$ as a positive-definite Hankel matrix for $n\le n^*$, so the matrix inverse is well defined [2503.01697]. The non-polynomial character of $B_n^{(\mathsf{Kry})}$ thus enters through the inverse of a moment matrix built from polynomial data.

The measurement protocol is local. One samples
$$
U=\bigotimes_{j=1}^N u_j,
$$
where each $u_j$ is drawn independently from a suitable single-qubit ensemble such as the single-qubit Clifford group or $U(2)$, applies $U$ to $\rho$, measures in the computational basis, and records a bitstring $s$ [2503.01697]. Repeating this procedure yields classical shadows
$$
\hat{\rho}^{(m)}=\bigotimes_{j=1}^N\left(3u_j^{(m)\dagger}|s_j^{(m)}\rangle\langle s_j^{(m)}|u_j^{(m)}-\openone\right),
$$
with $\mathbb{E}[\hat{\rho}^{(m)}]=\rho$. The moments $T_k$ can be estimated simultaneously from the same shadow data using $U$-statistics, and the batch-shadow method can accelerate post-processing [2602.17361]. The estimator in the original proposal uses a median-of-means construction and finally computes
$$
\hat B_n^{(\mathsf{Kry})}=\hat b^{\,T}\hat A^{-1}\hat b
$$
[2503.01697].

The formal guarantee states that for
$$
I=8\ln\frac{2n}{\delta},
$$
and $L$ chosen according to an explicit bound depending on the operator norms or variances of the relevant $O_l^{(k+2)}$, the final estimate satisfies
$$
\abs{\hat{B}_n^{(\mathsf{Kry})}-B_n^{(\mathsf{Kry})}<2n\epsilon\max_{0\leq k\leq 2n-1}\abs{\frac{\partial}{\partial T_k}B_n^{(\mathsf{Kry})}
$$
with probability at least $1-\delta$ [2503.01697]. The papers emphasize that the protocol uses only single-qubit gates and local measurements, and that its measurement cost is comparable to the best existing randomized-measurement approaches for QFI lower bounds [2503.01697].

At the same time, the resource cost grows exponentially with Krylov order, which is why the practical emphasis falls on low-order bounds [2602.17361]. This is the main limitation of current KST implementations. The method is therefore most attractive precisely in the regimes highlighted by the theory: moderate Krylov order, rapid convergence, or low-rank structure yielding early exactness. Numerical simulations support this picture: low-order Krylov bounds are close to $F_Q$ already for moderate mixed states, outperform Taylor bounds in direct comparisons, and for rank-2 examples the highest relevant Krylov bound matches the QFI exactly [2602.17361]. The original benchmarks further reported that $\hat B_1^{(\mathsf{Kry})}$ for pseudo-pure states obtained from $4.8\times 10^5$ classical shadows tracks $F_Q$ well, with sample complexity for relative error $0.1$ scaling like $2^{0.8N}$, while for the bound entangled example numerical results using $3.2\times 10^6$ classical shadows showed sample cost scaling like $2^{1.2N}$ [2503.01697].

KST is therefore best understood as a framework for QFI estimation that combines three features: a Krylov-subspace approximation to the SLD-related operator $L$, a strict hierarchy of non-polynomial lower bounds $B_n^{(\mathsf{Kry})}$, and shadow-based estimation of the polynomial moments required to evaluate those bounds [2503.01697]. In quantum metrology, entanglement detection, sensing, variational algorithms, and resource quantification, its importance lies in reducing or removing the systematic error inherent to polynomial lower-bound methods while preserving experimental accessibility [2602.17361].

Source: https://www.emergentmind.com/topics/krylov-shadow-tomography-kst