Pure-State Tomography: Sample-Optimal Methods

• The paper establishes protocols that achieve the theoretical limits of pure-state tomography by saturating the Massar–Popescu bound with collective measurements.
• It details efficient sequential weak measurement strategies and optimal POVM constructions, enabling sample-optimal state reconstruction under various measurement models.
• The work compares methodologies, including classical shadows and oracle-based techniques, demonstrating quadratic improvements and optimal scaling in sample complexity.

Sample-optimal pure-state tomography concerns achieving the minimal possible number of quantum samples (state copies or oracle calls) required to reconstruct or characterize an unknown pure quantum state to a specified accuracy, under various measurement models. This topic encompasses both the fundamental limits—defined precisely by information-theoretic bounds—and explicit, efficiently implementable protocols that saturate these limits. Theoretical and algorithmic developments in this area illuminate the optimal strategies for both full state reconstruction (in ℓ_q-norm or fidelity) and compressed representations such as classical shadows for observable estimation.

1. Theoretical Limits: Fidelity and Sample Complexity

For pure-state qubit tomography with access to $N$ identical copies of an unknown $|\psi\rangle$, the maximal achievable average fidelity of any collective measurement is given by the Massar–Popescu (MP) bound:

$F_N^{\mathrm{opt}} = \frac{N+1}{N+2}$

where the average infidelity $1-F_N^{\mathrm{opt}} = 1/(N+2) = O(1/N)$. This bound is tight and applies to any collective, symmetric measurement on the total spin-$J=N/2$ subspace. Any sample-optimal protocol must saturate this bound (Shojaee et al., 2018).

For general $d$-dimensional systems and observables $O$ with $\mathrm{Tr}(O^2)\leq B$ and $\|O\|\leq 1$, the minimal sample complexity for estimating $\mathrm{Tr}(O\rho)$ of a pure state $\rho$ within error $\epsilon$ is (Grier et al., 2022):

• Joint measurement: $\tilde{\Theta}(\sqrt{B}/\epsilon + 1/\epsilon^2)$
• Independent measurement: $\mathcal{O}(\sqrt{Bd}/\epsilon + 1/\epsilon^2)$

Full state reconstruction in $\ell_2$-norm to error $\epsilon$ requires $\widetilde{\Theta}(d/\epsilon)$ applications of a state-preparation unitary and its inverse when available (Apeldoorn et al., 2022).

2. The Spin-Coherent-State POVM and Optimal Measurements

For ensembles of qubits, the optimal POVM for pure-state tomography is the spin-coherent-state (SCS) POVM:

$$E(\mathbf{n}) = \frac{2J+1}{4\pi} |J, J\rangle_{\mathbf{n}}\langle J, J|_{\mathbf{n}}$$

where $|J, J\rangle_{\mathbf{n}}$ is the spin-coherent state pointing in direction $\mathbf{n}\in S^2$, obtained by rotating the Dicke state $|J, J\rangle_z$. This POVM is overcomplete on the symmetric subspace and is uniquely optimal: any other POVM on the symmetric subspace achieves average fidelity at most equal to the MP bound, with equality if and only if each $E_r\propto|J,J\rangle_{n_r}\langle J,J|_{n_r}$ (Shojaee et al., 2018).

Estimators that maximize fidelity choose the output Bloch vector as

$$\mathbf{n}_\text{est}(r) = \frac{\mathrm{Tr}[E_r J]}{|\mathrm{Tr}[E_r J]|}$$

where $J$ is the total spin operator.

3. Sequential Weak Measurement Protocols

The SCS-POVM, despite its theoretical optimality, is not directly physical. However, Shojaee et al. showed that it is operationally realizable via sequential weak measurements of the collective spin in randomized directions (Shojaee et al., 2018).

Protocol Structure

• Each weak measurement in direction $\mathbf{u}$ applies the Kraus operator:

$$\delta K_{\mathbf{u}}(m) = (\frac{\kappa\delta t}{2\pi})^{1/4} \exp\left[-\frac{\kappa\delta t}{4}(J_{\mathbf{u}} - m)^2\right]$$

• Squeezing (due to the quadratic term) is canceled by averaging blocks of measurements over isotropically sampled directions.
• In the limit of many weak steps, the block evolution is governed by:

$$K_\text{block}(\boldsymbol{\mu}) \propto \exp\left[(\kappa\Delta t/2)\,\boldsymbol{\mu}\cdot J\right]$$

• As the total time $T$ increases, Kraus operators concentrate the measurement outcome onto a single spin-coherent state:

$K(T) = U(T) \exp[\alpha(T)\, \mathbf{n}(T)\cdot J]$

The POVM element after time $T$ is $E(\mu) = K(T)^\dagger K(T)$, which strongly projects onto $|J,\pm J\rangle_{\mathbf{n}(\infty)}$ as $\alpha(T)$ grows.

Scaling

The number of required weak measurement steps (or total measurement time $T$) to concentrate onto a rank-1 SCS projector at infidelity $1/(N+2)$ is $O(\log N)$, i.e., the protocol is polylogarithmic in $N$ for sample complexity (Shojaee et al., 2018).

4. Sample-Optimal Classical Shadows for Pure States

Classical shadows are randomized measurement protocols whose classical output suffices for the estimation of arbitrary observables after state preparation. For pure states, the sample-optimal regime displays quadratic improvement over mixed states for joint measurements (Grier et al., 2022).

Joint Measurement Protocol

• Measurement: symmetric subspace POVM on $s$ copies,

$$A_\psi = \kappa_s |\psi\rangle\langle\psi|^{\otimes s}\,d\psi$$

with $$\kappa_s = \text{dim}(\mathrm{sym}^s) = \binom{s+d-1}{d-1}$$.

• Estimator:

$$\hat\rho = \frac{(d+s)|\psi\rangle\langle\psi| - I}{s}$$

• For observable $O$ with $\|O\|\leq 1$, $\mathrm{Tr}(O^2)\leq B$,

$$\mathrm{Var}[\mathrm{Tr}(O \hat\rho)] \leq \frac{\mathrm{Tr}(O^2) + 8s\,\mathrm{Tr}(O^2\rho)}{s^2}$$

Sample Complexity

The minimal sample complexity is

$$s = \tilde{\Theta}(\sqrt{B}/\epsilon + 1/\epsilon^2)$$

for joint measurements, optimal up to logarithmic factors in $B$ and $d$. Independent measurement protocols achieve $s = O(\sqrt{Bd}/\epsilon + 1/\epsilon^2)$ (Grier et al., 2022).

Lower Bounds

Matching lower bounds are established via reductions to the Boolean Hidden Matching problem (for the $\sqrt{B}/\epsilon$ term) and two-state distinguishability (for the $1/\epsilon^2$ term), as well as Holevo’s theorem on information capacity.

5. State-Preparation-Oracular Tomography

In scenarios where a unitary oracle $U$ prepares the pure state $|\psi\rangle = U|0\rangle$ and $U^\dagger$ is available, optimal tomography can be achieved with minimal quantum queries (Apeldoorn et al., 2022).

• $\ell_q$-norm estimation to error $\epsilon$ is achieved with $\widetilde{\Theta}(d/\epsilon)$ calls to $U$ and $U^\dagger$ for $q=2$, saturating the quantum lower bound.
• The algorithm first performs an $\ell_\infty$-approximation of the amplitude vector $\alpha_j$ (via multidimensional phase estimation and block encoding), then recovers phases, and finally converts the estimate to a $\ell_q$ proxy without additional dimension-dependent factors.
• For higher precision ($\epsilon \ll 1/\sqrt{d}$), the query cost is $O(\sqrt{d}/\epsilon)$, with a lower bound also given by reduction to phase-oracle recovery problems.

6. Practical Implementation Considerations

Implementations in atomic-ensemble settings (e.g., Faraday/QND probing in high-cooperativity cavities) achieve $\kappa T\approx 1$ with negligible decoherence, supporting feasibility of the sequential weak measurement protocol (Shojaee et al., 2018). The protocol is robust to deviations from perfect isotropy, with errors only at subleading $1/N$ corrections. Both discrete (randomized axis) and continuous (fixed measurement axis combined with rapid random Euler rotations) protocols realize isotropic POVMs.

For classical shadow protocols, random Clifford measurements (3-design POVMs) suffice for independent measurements, but optimal joint measurement performance requires symmetric subspace measurements. The quadratic improvement in $B$ for pure states is not attainable by the original Huang–Kueng–Preskill approach for mixed states (Grier et al., 2022).

7. Comparative Summary of Sample Complexities

The following table summarizes key sample complexity scalings for pure-state tomography protocols.

Task/Protocol	Model	Sample Complexity	Reference
Qubit pure-state fidelity	Collective	$1/(N+2)$ infidelity ($N$ copies)	(Shojaee et al., 2018)
Observable estimation	Joint measurement	$\tilde{\Theta}(\sqrt{B}/\epsilon + 1/\epsilon^2)$	(Grier et al., 2022)
Observable estimation	Independent	$O(\sqrt{Bd}/\epsilon + 1/\epsilon^2)$	(Grier et al., 2022)
Full state (ℓ₂-norm)	Unitary oracle	$\widetilde{\Theta}(d/\epsilon)$	(Apeldoorn et al., 2022)

These results establish the operational and information-theoretic limits for pure-state tomography across several practical settings, demonstrating that variants of the SCS-POVM, classical-shadow protocols, and oracle-based methods all admit sample-optimal strategies, with explicit constructions and matching lower bounds.