Optimal high-precision shadow estimation
Sitan Chen, Jerry Li, Allen Liu
TL;DR
The paper addresses the problem of estimating many linear properties of an unknown $d$-dimensional quantum state from as few copies as possible, focusing on the high-accuracy regime where $\varepsilon = O(d^{-c})$. It introduces a representation-theoretic approach based on Keyl's POVM and Schur-Weyl duality, together with a splitting-based dimensionality-reduction that reduces to balanced states and linearizes around the maximally mixed state $I_d/d$. The main contributions are the first tight, nontrivial sample complexity bounds for shadow tomography and classical shadows in this regime, achieving a dimension-free rate $n = O(\log m / \varepsilon^2)$ for general observables and proving an equivalence between the quantum and classical shadow tasks under these conditions, along with reductions to the medium-accuracy regime. This advances the theoretical understanding of shadow estimation and yields practically efficient protocols for high-precision quantum state profiling with broad implications for quantum data analysis and benchmarking.
Abstract
We give the first tight sample complexity bounds for shadow tomography and classical shadows in the regime where the target error is below some sufficiently small inverse polynomial in the dimension of the Hilbert space. Formally we give a protocol that, given any $m\in\mathbb{N}$ and $ε\le O(d^{-12})$, measures $O(\log(m)/ε^2)$ copies of an unknown mixed state $ρ\in\mathbb{C}^{d\times d}$ and outputs a classical description of $ρ$ which can then be used to estimate any collection of $m$ observables to within additive accuracy $ε$. Previously, even for the simpler task of shadow tomography -- where the $m$ observables are known in advance -- the best known rates either scaled benignly but suboptimally in all of $m, d, ε$, or scaled optimally in $ε, m$ but had additional polynomial factors in $d$ for general observables. Intriguingly, we also show via dimensionality reduction, that we can rescale $ε$ and $d$ to reduce to the regime where $ε\le O(d^{-1/2})$. Our algorithm draws upon representation-theoretic tools recently developed in the context of full state tomography.
