Heisenberg-limited adaptive gradient estimation for multiple observables

TL;DR

This paper introduces an adaptive quantum algorithm that estimates multiple observables' expectation values at the Heisenberg limit with sublinear scaling in the number of observables, optimizing resource use in large quantum systems.

Contribution

The work presents a novel adaptive quantum algorithm achieving Heisenberg-limited precision for multiple observables with efficient resource scaling and stability improvements over existing methods.

Findings

Achieves Heisenberg limit $1/\varepsilon$ in estimation precision.

Scales as $\mathcal{O}(\varepsilon^{-1}\sqrt{M}\log M)$ queries for $M$ observables.

Resource overhead is linear in $M$, independent of precision $\varepsilon$.

Abstract

In quantum mechanics, measuring the expectation value of a general observable has an inherent statistical uncertainty that is quantified by variance or mean squared error of measurement outcome. While the uncertainty can be reduced by averaging several samples, the number of samples should be minimized when each sample is very costly. This is especially the case for fault-tolerant quantum computing that involves measurement of multiple observables of non-trivial states in large quantum systems that exceed the capabilities of classical computers. In this work, we provide an adaptive quantum algorithm for estimating the expectation values of $M$ general observables within root mean squared error $\varepsilon$ simultaneously, using $\mathcal{O}(\varepsilon^{-1}\sqrt{M}\log M)$ queries to a state preparation oracle of a target state. This remarkably achieves the scaling of Heisenberg limit $1/\varepsilon$, a fundamental bound on the estimation precision in terms of mean squared error, together with the sublinear scaling of the number of observables $M$. The proposed method is an adaptive version of the quantum gradient estimation algorithm and has a resource-efficient implementation due to its adaptiveness. Specifically, the space overhead in the proposed method is $\mathcal{O}(M)$ which is independent from the estimation precision $\varepsilon$ unlike non-iterative algorithms. In addition, our method can avoid the numerical instability problem for constructing quantum circuits in a large-scale task (e.g., $\varepsilon\ll 1$ in our case), which appears in the actual implementation of many algorithms relying on quantum signal processing techniques. Our method paves a new way to precisely understand and predict various physical properties in complicated quantum systems using quantum computers.