---
title: 'Pure-State Tomography: Sample-Optimal Methods'
url: https://www.emergentmind.com/topics/sample-optimal-pure-state-tomography
type: topic
---

# Pure-State Tomography: Sample-Optimal Methods

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 [1805.01012].

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 [2211.11810]:
- 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 [2207.08800].

## 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}$ [1805.01012].

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 [1805.01012].

### 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 [1805.01012].

## 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 [2211.11810].

### 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)$ [2211.11810].

### 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 [2207.08800].

- $\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 [1805.01012]. 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 [2211.11810].

## 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)         | [1805.01012]     |
| Observable estimation       | Joint measurement    | $\tilde{\Theta}(\sqrt{B}/\epsilon + 1/\epsilon^2)$ | [2211.11810]     |
| Observable estimation       | Independent          | $O(\sqrt{Bd}/\epsilon + 1/\epsilon^2)$    | [2211.11810]     |
| Full state (ℓ₂-norm)        | Unitary oracle       | $\widetilde{\Theta}(d/\epsilon)$          | [2207.08800]     |

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.

Source: https://www.emergentmind.com/topics/sample-optimal-pure-state-tomography