Optimal algorithmic complexity of inference in quantum kernel methods
A new algorithm cuts the query cost of quantum kernel inference to a proven minimum — but fault-tolerant hardware is still years away.

The Thesis
Quantum kernel methods are a family of machine-learning algorithms that use quantum computers to compute similarity scores between data points. The bottleneck has been inference: every prediction on new data requires summing potentially millions of these similarity scores, and the standard approach is painfully slow. This paper proves that encoding the entire sum as a single quantum measurement, combined with quantum amplitude estimation (a technique that extracts probabilities quadratically faster than repeated sampling), achieves the best possible query count — and no future algorithm can do better. The catch is hardware: the approach requires early fault-tolerant quantum devices that don't yet exist at useful scale.
Catalyst
Quantum amplitude estimation has matured from a theoretical curiosity into a concrete near-term subroutine, with several groups demonstrating small-scale implementations on current hardware. Simultaneously, the kernel-method community has produced more realistic estimates of the α coefficients (the trained model weights) that govern inference cost, making it tractable to compare algorithmic strategies on the same footing. The convergence of these two threads made a rigorous, complete complexity landscape possible now.
What's New
The standard approach estimates each of the N kernel values independently via quantum sampling, giving query complexity that grows with both N and 1/ε², where ε is the required precision. Prior improvements addressed one axis at a time — either swapping in amplitude estimation for individual terms, or batching terms differently — without proving whether the result was optimal. This paper jointly optimizes both axes, proves a matching lower bound showing O(‖α‖₁/ε) queries cannot be beaten, and additionally compares gate complexity (total circuit operations) against query complexity, showing the theoretically optimal query strategy is not always the best practical choice depending on hardware.
The Counter
This paper solves an optimization problem inside a bigger problem that may not matter. Quantum kernel methods have yet to show a convincing advantage over classical kernels or modern neural networks on any practically relevant dataset — the field is still searching for a problem where they're the right tool. Even granting that quantum kernels will someday be useful, this result optimizes inference query complexity while explicitly noting that gate complexity — what actually limits real hardware — does not always favor the same strategy. The required subroutine, quantum amplitude estimation, needs coherent quantum memory and low error rates that current and near-term devices cannot sustain for problem sizes where classical methods already struggle. A practitioner building a product today has no hardware on which to run this algorithm. And the lower bound, while elegant, applies within a specific oracle model; a different algorithmic framework or a different classical-quantum hybrid might sidestep the bound entirely. The paper is a clean theoretical contribution, but the distance from 'query-optimal proof' to 'deployed advantage' is enormous.
Longs
- IBM (quantum hardware and Qiskit ecosystem, directly relevant to early-fault-tolerant implementations)
- IonQ (IONQ — trapped-ion hardware targets early fault-tolerant regimes where amplitude estimation is feasible)
Shorts
- Classical kernel-method vendors whose competitive moat assumes inference stays slow enough that approximation shortcuts (like random Fourier features) remain necessary — a quadratic quantum speedup in ε would erode that justification if hardware arrives
Enablers (Picks & Shovels)
- Quantum amplitude estimation libraries (e.g., Qiskit's IterativeAmplitudeEstimation module)
- Open quantum hardware access programs (IBM Quantum Network, IonQ cloud) needed to benchmark these algorithms
- Kernel method benchmarking datasets such as those in the QML benchmark suite used to study realistic α distributions
Private Watchlist
- Quantinuum (trapped-ion quantum computing, early fault-tolerant focus)
- Phasecraft (quantum algorithms for near-term hardware)
- Classiq (quantum circuit compilation tools that would be needed to implement these subroutines)
The Paper
Quantum kernel methods are among the leading candidates for achieving quantum advantage in supervised learning. A key bottleneck is the cost of inference: evaluating a trained model on new data requires estimating a weighted sum $\sum_{i=1}^N α_i k(x,x_i)$ of $N$ kernel values to additive precision $\varepsilon$, where $α$ is the vector of trained coefficients. The standard approach estimates each term independently via sampling, yielding a query complexity of $O(N\lVertα\rVert_2^2/\varepsilon^2)$. In this work, we identify two independent axes for improvement: (1) How individual kernel values are estimated (sampling versus quantum amplitude estimation), and (2) how the sum is approximated (term-by-term versus via a single observable), and systematically analyze all combinations thereof. The query-optimal combination, encoding the full inference sum as the expectation value of a single observable and applying quantum amplitude estimation, achieves a query complexity of $O(\lVertα\rVert_1/\varepsilon)$, removing the dependence on $N$ from the query count and yielding a quadratic improvement in both $\lVertα\rVert_1$ and $\varepsilon$. We prove a matching lower bound of $Ω(\lVertα\rVert_1/\varepsilon)$, establishing query-optimality of our approach up to logarithmic factors. Beyond query complexity, we also analyze how these improvements translate into gate costs and show that the query-optimal strategy is not always optimal in practice from the perspective of gate complexity. Our results provide both a query-optimal algorithm and a practically optimal choice of strategy depending on hardware capabilities, along with a complete landscape of intermediate methods to guide practitioners. All algorithms require only amplitude estimation as a subroutine and are thus natural candidates for early-fault-tolerant implementations.