Compressed Sensing Shadow Tomography

Estimating many local expectation values over time is a central measurement bottleneck in quantum simulation and device characterization. We study the task of reconstructing the Pauli-signal matrix S_ij=Tr\(O_i ρ(t_j\)) for a collection of M low-weight Pauli observables \{O_i\}_i=1^M over N timesteps \{t_j\}_j=1^N, while minimizing the total number of device shots. We propose a Compressed Sensing Shadow Tomography (CSST) protocol that combines two complementary reductions. First, local classical shadows reduce the observable dimension by enabling many Pauli expectation values to be estimated from the same randomized snapshots at a fixed time. Second, compressed sensing reduces the time dimension by exploiting the fact that many expectation-value traces are spectrally sparse or compressible in a unitary (e.g., Fourier) transform basis. Operationally, CSST samples mll N timesteps uniformly at random, collects shadows only at those times, and then reconstructs each length-N signal via standard ell₁-based recovery in the unitary transform domain. We provide end-to-end guarantees that explicitly combine shadow estimation error with compressed sensing recovery bounds. For exactly s-sparse signals in a unitary transform basis, we show that m=O left\(slog² s log Nright\) random timesteps suffice (with high probability), leading to total-shot savings scaling as widetildeΘ(N/s) (i.e., up to polylogarithmic factors) relative to collecting shadows at all N timesteps. For approximately sparse signals, the reconstruction error decomposes into a compressibility (tail) term plus a noise term. We present numerical experiments on noisy many-qubit dynamics that support strong Fourier compressibility of Pauli traces and demonstrate substantial shot reductions with accurate reconstruction.