Full tomography of an unknown quantum evolution is resource-intensive and often unnecessary when the goal is only to predict selected properties. This motivates the study of classical shadow estimation of unitary channels (CSEU), a task in which one queries an unknown $d$-dimensional unitary $U$ and stores classical data that can later be used to predict expectation values $\mathrm{tr}[O UρU^\dagger]$ for arbitrary input states $ρ$ and observables $O$. Here we propose a parallel and non-adaptive CSEU protocol that uses $\mathcal{O}(\min\{d\varepsilon^{-2},d^{3/2}\varepsilon^{-1}\})$ queries to $U$ when the observables have constant effective size. This achieves Heisenberg scaling with respect to the additive error $\varepsilon$ and is query-optimal, as we prove a matching lower bound $Ω(\min\{d\varepsilon^{-2},d^{3/2}\varepsilon^{-1}\})$ for the query complexity, which remains valid even when stronger access to the unknown unitary is available. Our CSEU protocol is a versatile primitive with applications to several quantum learning tasks. In particular, it yields a near-optimal algorithm for learning an arbitrary Hamiltonian from its real-time evolution in the high-precision regime and recovers the best-known bounds for several tasks in a unified way.
Full tomography of an unknown quantum evolution is resource-intensive and often unnecessary when the goal is only to predict selected properties. This motivates the study of classical shadow estimation of unitary channels (CSEU), a task in which one queries an unknown d-dimensional unitary U and stores classical data that can later be used to predict expectation values tr[OUρU†] for arbitrary input states ρ and observables O. Here we propose a parallel and non-adaptive CSEU protocol that uses O(min{dε−2,d3/2ε−1}) queries to U when the observables have constant effective size. This achieves Heisenberg scaling with respect to the additive error ε and is query-optimal, as we prove a matching lower bound Ω(min{dε−2,d3/2ε−1}) for the query complexity, which remains valid even when stronger access to the unknown unitary is available. Our CSEU protocol is a versatile primitive with applications to several quantum learning tasks. In particular, it yields a near-optimal algorithm for learning an arbitrary Hamiltonian from its real-time evolution in the high-precision regime and recovers the best-known bounds for several tasks in a unified way.
Classical shadow estimation (CSE) is a powerful tool for learning the properties of quantum states and quantum processes. Here we consider the CSE task for quantum unitary channels. By querying an unknown unitary channel U multiple times in quantum experiments, the goal is to learn a classical description from which one can accurately predict many different linear properties of the channel, i.e., the expectation values of arbitrary observables measured on the output of U upon arbitrary input states. Based on collective measurements on multiple systems, we propose a query-efficient protocol for this task, whose query complexity has a quadratic advantage over the previous best approach for this problem, and almost saturates the information-theoretic lower bound. To further enhance practicality, we also present a variant protocol using only single-copy measurements, which still offers much better query performance than previous protocols that do not use quantum memory, and can serve as a key subroutine for learning an arbitrary unknown Hamiltonian from dynamics. In addition to linear properties of unitary channels, our protocol can also be applied to simultaneously predict many nonlinear properties, such as out-of-time-ordered correlators.
In the study of quantum process tomography, it has remained open how to learn a unitary channel in diamond distance with strictly parallel queries as efficiently as with sequential queries. Sequential queries allow one to exploit adaptivity to refine the estimate and adjust the learning strategy accordingly, whereas such adjustments are impossible for parallel queries. From the perspective of higher-order quantum operations, any parallel learning strategy can be simulated by a sequential one, whereas the converse does not hold in general, suggesting that sequential learning strategies are potentially more powerful. Contrary to this intuition, we present a new protocol for learning an unknown d-dimensional unitary channel to within ε in diamond distance using O(d2/ε) parallel queries, matching the lower bound presented in [Haah, Kothari, O'Donnell, and Tang, FOCS '23]. Our protocol thus closes the gap between parallel and sequential strategies for unitary channel tomography. As applications, it achieves optimal query complexity for boundary-regime quantum channel tomography and improves the query efficiency of the best-known protocol for tomography of fermionic linear optics.
Reversing unknown quantum dynamics is vital for quantum control and learning, yet full unitary inversion requires a resource-intensive O(d2) queries. Because many applications only require the reversed evolution for a given observable, a critical gap remains in understanding the minimal resources for such targeted reversal. Here, we address this by introducing shadow unitary inversion. We establish a lower bound showing the query complexity must scale at least linearly with system dimension for spectrally biased observables, with the constant determined by the observable’s spectral properties. For qubit case, we construct an explicit, deterministic three-query sequential protocol achieving exact shadow inversion and completely characterize all admissible channels, with numerical evidence suggesting optimality. For higher dimensions, we develop a semidefinite-programming formulation and introduce a representation-theoretic symmetry reduction that decomposes the optimization into invariant blocks, substantially reducing the problem size. Shadow unitary inversion thus offers a resource-efficient path to inverse-dynamics estimation for future quantum control, diagnostics, and learning tasks. Reversing unitary operations is crucial in quantum computing and control, yet remains a challenging task. The authors introduce shadow unitary inversion, achieving efficient inversion at the expectation value level with reduced query complexity, offering a promising approach for quantum error correction and information recovery.
How many black-box queries to a quantum channel are needed to learn its full classical description? This question lies at the heart of quantum channel tomography (also known as quantum process tomography), a fundamental task in the characterization and validation of quantum hardware. Despite extensive prior work, the optimal query complexity for quantum channel tomography is far from fully understood.
In this paper, we study tomography of an unknown quantum channel with input dimension d1, output dimension d2, and Kraus rank at most r, to within error ε. We identify the dilation rate τ=rd2/d1 (which always satisfies τ≥1 due to the trace preservation of quantum channels) as a key parameter, and establish that the optimal query complexity of channel tomography exhibits distinct scaling laws across three regimes of τ.
- In the boundary regime (τ=1): we show that the query complexity is Θ(rd1d2/ε) for Choi trace norm error ε, and is upper bounded by O(min{rd11.5d2/ε,rd1d2/ε2}) and lower bounded by Ω(rd1d2/ε) for diamond norm error ε.
- In the away-from-boundary regime (τ≥1+Ω(1)): we show that the query complexity is Θ(rd1d2/ε2) for both Choi trace norm and diamond norm errors ε.
Our results uncover a sharp Heisenberg-to-classical phase transition in the query complexity of quantum channel tomography: at τ=1, the optimal query complexity exhibits Heisenberg scaling 1/ε, whereas for τ≥1+Ω(1), it exhibits classical scaling 1/ε2. In addition, we show that in the near-boundary regime (1<τ<1+o(1)), the query complexity exhibits a mixture of Heisenberg and classical scaling behaviors.
We study the sample complexity of shadow tomography in the high-precision regime under realistic measurement constraints. Given an unknown d-dimensional quantum state ρ and a known set of observables {Oi}i=1m, the goal is to estimate expectation values {tr(Oiρ)}i=1m to accuracy ε in Lp-norm, using possibly adaptive measurements that act on O(polylog(d)) number of copies of ρ at a time. We focus on the regime where ε is below an instance-dependent threshold.
Our main contribution is an instance-optimal characterization of the sample complexity as Θ~(Γp/ε2), where Γp is a function of {Oi}i=1m defined via an optimization formula involving the inverse Fisher information matrix. Previously, tight bounds were known only in special cases, e.g. Pauli shadow tomography with L∞-norm error. Concretely, we first analyze a simpler oblivious variant where the goal is to estimate an observable of the form ∑i=1mαiOi with ∥α∥q=1 (where q is dual to p) revealed after the measurement. For single-copy measurements, we obtain a sample complexity of Θ(Γpob/ε2). We then show Θ~(Γp/ε2) is necessary and sufficient for the original problem, with the lower bound applying to unbiased, bounded estimators. Our upper bounds rely on a two-step algorithm combining coarse tomography with local estimation. Notably, Γ∞ob=Γ∞. In both cases, allowing c-copy measurements improves the sample complexity by at most Ω(1/c).
Our results establish a quantitative correspondence between quantum learning and metrology, unifying asymptotic metrological limits with finite-sample learning guarantees.
We consider process tomography for unitary quantum channels. Given access to an unknown unitary channel acting on a d-dimensional qudit, we aim to output a classical description of a unitary that is ε-close to the unknown unitary in diamond norm. We design an algorithm achieving error ε using O(d2/ε) applications of the unknown channel and only one qudit. This improves over prior results, which use O(d3/ε2) [via standard process tomography] or O(d2.5/ε) [Yang, Renner, and Chiribella, PRL 2020] applications. To show this result, we introduce a simple technique to “bootstrap” an algorithm that can produce constant-error estimates to one that can produce ε-error estimates with the Heisenberg scaling. Finally, we prove a complementary lower bound showing that estimation requires Ω(d2/ε) applications, even with access to the inverse or controlled versions of the unknown unitary. This shows that our algorithm has both optimal query complexity and optimal space complexity.