---
title: Projected Least Squares POVM Tomography
url: https://www.emergentmind.com/topics/projected-least-squares-povm-tomography
type: topic
---

# Projected Least Squares POVM Tomography

Searching arXiv for the cited papers to ground the article in the relevant literature.
arXiv search: 2507.04500 "Fast quantum measurement tomography with dimension-optimal error bounds"
Projected least squares POVM tomography is a two-step protocol for reconstructing an unknown quantum measurement from data obtained with a known informationally complete probe ensemble. In the formulation introduced in "Fast quantum measurement tomography with dimension-optimal error bounds" [2507.04500], one first applies least-squares estimation to obtain an unconstrained approximation of each POVM element and then projects the resulting collection onto the set of valid quantum measurements. For an \(L\)-outcome POVM acting on a \(d\)-dimensional system, the method is designed to retain light classical co-processing while achieving non-asymptotic error guarantees and sample complexity that is sample-optimal in the system dimension \(d\) up to logarithmic factors [2507.04500].

## 1. Definition and estimation framework

The target of inference is an unknown \(L\)-outcome POVM
\[
E=\{E_1,\dots,E_L\}\subset \mathbb{C}^{d\times d},
\]
with \(E_j\ge 0\) and \(\sum_{j=1}^L E_j=I\). The protocol assumes access to a known informationally complete ensemble \(\{\ket{\psi_i}\}_{i=1}^M\). In each of \(N\) independent rounds, one chooses \(i\in\{1,\ldots,M\}\) uniformly at random, prepares \(\ket{\psi_i}\), measures \(E\), and records the observed outcome \(j\). The empirical frequencies are
\[
\hat f_{ij}=\frac{\text{number of times “prepare }i\text{, observe }j”}{N},
\]
and satisfy
\[
\mathbb{E}[\hat f_{ij}]=\frac{1}{M}\bra{\psi_i}E_j\ket{\psi_i}=[\mathcal{M}(E_j)]_i,
\]
where \(\mathcal{M}:\mathbb{C}^{d\times d}\to\mathbb{R}^M\) is the linear measurement map [2507.04500].

The first step estimates each outcome operator independently by unconstrained least squares:
\[
\hat E_j^{\rm raw}=(\mathcal{M}^\dagger\mathcal{M})^{-1}\bigl(\mathcal{M}^\dagger \hat f_j\bigr),
\]
equivalently as the unique \(X\) minimizing \(\|\mathcal{M}(X)-\hat f_j\|_2^2\). The second step enforces physicality jointly across all outcomes by solving
\[
\{E_j^*\}_{j=1}^L=\underset{\{Z_j\ge 0,\;\sum_j Z_j=I\}}{\arg\min}\; d\bigl(\{\hat E_j^{\rm raw}\},\{Z_j\}\bigr),
\]
where \(d(\cdot,\cdot)\) is a chosen metric on POVMs [2507.04500].

A recurrent misconception is to identify POVM projection with separate positivity corrections on the individual \(\hat E_j^{\rm raw}\). The formulation above is stricter: the projection is onto the full feasible set \(\{Z_j\ge 0,\ \sum_j Z_j=I\}\), so the normalization constraint is imposed collectively rather than effect-by-effect. This suggests that the method is best understood as a constrained estimation procedure on the measurement object as a whole, not as \(L\) decoupled matrix denoisings.

## 2. Two-step projected least squares construction

Algorithmically, the protocol is stated in terms of precomputed frame operators \(\{\nu_i\}_{i=1}^M\). The raw least-squares estimate can be written as
\[
\hat E_j^{\rm raw}=\sum_{i=1}^M \hat f_{ij}\,\nu_i,
\]
after which one solves the projection problem over valid POVMs. The resulting estimator is the physical POVM \(\{E_j^*\}_{j=1}^L\) [2507.04500].

The projection step is metric-dependent. The underlying paper allows either the worst-case operational distance or the average-case distance, and the projection can be carried out by an SDP or, for large \(L\), by block-coordinate or first-order methods [2507.04500]. This separates the procedure into a linear inversion stage and a convex feasibility-and-proximity stage.

The construction is closely related to projected least squares in quantum state tomography. In "Fast state tomography with optimal error bounds" [1809.11162], one first computes the least-squares estimator and then projects onto the state space; in "A comparative study of estimation methods in quantum tomography" [1901.07991], least squares and projected least squares are analyzed within a broader family of estimators that project data onto parameter spaces with respect to specific metrics. Projected least squares POVM tomography transfers that paradigm from density operators to measurement operators, but the geometry changes because the feasible set is a product of PSD constraints coupled by the identity-resolution condition.

## 3. Closed-form least squares for 2-design probes

A principal simplification occurs when the probe ensemble is a complex projective 2-design \(\{\ket{\psi_i}\}_{i=1}^M\subset \mathbb{C}^d\), satisfying
\[
\frac1M\sum_i \ket{\psi_i}\!\bra{\psi_i}^{\otimes 2}
=
\int_{\rm Haar}\ket{\phi}\!\bra{\phi}^{\otimes 2}\,d\phi
=
\frac{\mathrm{Sym}_2}{d(d+1)}.
\]
In this case,
\[
\mathcal{M}^\dagger \mathcal{M}(X)=\frac{X+\mathrm{tr}(X)}{d(d+1)M},
\qquad
(\mathcal{M}^\dagger \mathcal{M})^{-1}(Y)=dM\bigl((d+1)Y-\mathrm{tr}(Y)\bigr),
\]
and the raw estimator simplifies to
\[
\hat E_j^{\rm raw}=\sum_{i=1}^M \hat f_{ij}\,\nu_i,
\qquad
\nu_i=d(d+1)\ket{\psi_i}\!\bra{\psi_i}-d\,I.
\]
The paper also states that the method admits an analytic form when using local 2-designs, including tensor-product 2-design probes on \(n\) qubits [2507.04500].

This closed form is the measurement-tomography analogue of the 2-design formulas that make projected least squares especially attractive in state tomography. For state estimation, the same 2-design structure yields explicit inverses for \(M^\dagger M\) and numerically cheap linear inversion [1809.11162]. In the POVM setting, the same phenomenon removes any large-scale matrix inversion beyond the precomputed frame operators, which is central to the method’s intended low co-processing cost.

## 4. Distances and non-asymptotic guarantees

Two distances are used to quantify reconstruction error [2507.04500]. The operational, or worst-case, distance is
\[
d_{\rm op}(E,F)=\max_{\rho}\frac12\sum_j \bigl|\mathrm{tr}(E_j\rho)-\mathrm{tr}(F_j\rho)\bigr|
=
\max_{x\subseteq [L]}\Bigl\|\sum_{j\in x}(E_j-F_j)\Bigr\|.
\]
The average-case distance is
\[
d_{\rm av}^2(E,F)=\frac1{2d}\sum_{j=1}^L
\left(
\|E_j-F_j\|_F^2+\bigl(\mathrm{tr}(E_j-F_j)\bigr)^2
\right).
\]
The operational distance equals the maximum total-variation distance between outcome distributions over all possible input states. The average-case distance corresponds to the expected \(\ell_1\)-error of the outcome distribution, averaged over input states drawn from a spherical 4-design, in the sense that
\[
d_{\rm av}^2(E,F)\approx
\mathbb{E}_{\ket\psi\sim 4\text{-design}}
\left\|
\bigl(\mathrm{tr}(E_j\ket\psi\bra\psi)\bigr)_j
-
\bigl(\mathrm{tr}(F_j\ket\psi\bra\psi)\bigr)_j
\right\|_1^2 .
\]

For global 2-design probes, the non-asymptotic high-probability guarantees are as follows.

| Target metric | Guarantee | Sufficient sample size |
|---|---|---|
| \(d_{\rm op}(E,E^*)\le \epsilon\) | confidence \(1-\delta\) | \(\displaystyle N \ge \frac{8\bigl(d^3+d^2(1+\epsilon/6)\bigr)}{\epsilon^2}\ln\!\bigl(\tfrac{2^{L+1}d}{\delta}\bigr)\) |
| \(d_{\rm av}(E,E^*)\le \epsilon\) | confidence \(1-\delta\) | \(\displaystyle N \ge \frac{8L^2\bigl(d^2+d(1+\epsilon/3L)\bigr)}{\epsilon^2}\ln\!\bigl(\tfrac{4Ld}{\delta}\bigr)\) |

The corresponding asymptotic forms stated in the abstract are \(\mathcal{O}(d^3L\ln(d)/\epsilon^2)\) samples for worst-case distance and \(\mathcal{O}(d^2L^2\ln(dL)/\epsilon^2)\) samples for average-case distance [2507.04500]. For local \(n\)-qubit 2-designs, the theorem replaces \(d^3\) by \(10^n\) and \(d^2\) by \(4^n\) in the worst-case numerator, and replaces \(d^2\) by \(5^n\) and \(d\) by \(2^n\) in the average-case bound [2507.04500].

The lower bounds establish that any non-adaptive, single-copy POVM tomography protocol must use at least
\[
\Omega\!\left(\frac{d^3}{\epsilon^2}\right)
\quad\text{shots}
\]
to achieve accuracy \(\epsilon\) in \(d_{\rm op}\), and at least
\[
\Omega\!\left(\frac{d^2L}{\epsilon^2}\right)
\quad\text{shots}
\]
to achieve accuracy \(\epsilon\) in \(d_{\rm av}\) [2507.04500]. Accordingly, the protocol is sample-optimal in dimension \(d\) up to logarithmic factors, but the lower-bound comparison is explicitly restricted to non-adaptive, single-copy schemes. This suggests that the optimality claim is dimension-theoretic rather than universal across all conceivable tomography models.

## 5. Proof architecture and computational profile

The proof strategy follows the same general template as projected least squares state tomography, but now at the level of measurement effects. For the worst-case bound, one shows by matrix-Bernstein that
\[
\|\hat E_j^{\rm raw}-E_j\|\lesssim \frac{\epsilon}{2}\qquad \forall j
\]
once \(N\sim d^3L/\epsilon^2\ln(dL/\delta)\). A union bound over \(2^L\) subsets then controls \(d_{\rm op}(E,\hat E^{\rm raw})\le \epsilon/2\), and the projection onto POVMs can only improve the worst-case distance by triangle inequality and the definition of projection [2507.04500]. For the average-case theorem, the Frobenius part and trace part in \(d_{\rm av}\) are bounded separately via matrix-Bernstein and Hoeffding’s inequality and then combined [2507.04500].

The lower bounds are obtained by packing-and-Fano arguments: one constructs an exponentially large family of well-separated POVMs, encodes a random label, upper-bounds the information gained per measurement via a \(\chi^2\)-divergence argument, and applies Fano’s inequality [2507.04500].

The classical post-processing cost is polynomial in \((d,L)\). Given \(\{\hat f_{ij}\}\) and precomputed \(\{\nu_i\}\), each raw \(E_j^{\rm raw}\) is a weighted sum of \(M\) rank-1 terms, with cost \(O(Md^2)\) per outcome or \(O(Md^2L)\) overall; if \(M\sim d^2\) for a global 2-design, this becomes \(O(d^4L)\). The projection step is an SDP of size \(d^2L\). Off-the-shelf interior-point solvers run in \(O(d^6L^3)\) worst-case, while first-order or alternating-projection schemes reduce it in practice to roughly \(O(d^4L)\) per iteration [2507.04500]. Thus the method avoids any large-scale matrix inversion beyond the closed-form 2-design formulas, but its overall runtime still depends materially on the projection stage.

## 6. Empirical demonstration and related developments

The protocol has been demonstrated on two flux-tunable transmon qubits with chip \(T_1\approx 30\,\mu{\rm s}\), \(T_2\approx 10\,\mu{\rm s}\), readout fidelities \(\approx 97\)–\(98\%\), and single-qubit gate errors \(\approx 10^{-3}\) [2507.04500]. In one experiment, a single-qubit SIC-POVM was implemented via an ancilla qubit and controlled rotations. Using \(N=1.66\times 10^5\) random input states from a single-qubit 2-design, each POVM element was recovered with worst-case distance \(\sim 10^{-2}\) from the ideal. The same study reconstructed the corresponding half-sided measurement channel in the Pauli basis, which was needed for classical-shadow error-mitigation, and reported that the experimental inversion closely matches the ideal depolarizing-channel inversion. Additional two-qubit POVMs—identity, Hadamard, Bell, and Haar-random—were studied with \(N=10^6\) shots, showing high fidelity to the targets and behavior consistent with the \(O(1/\sqrt{N})\) scaling predicted by the non-asymptotic bounds [2507.04500].

Within the broader projected least-squares literature, this work extends a pattern already visible in state tomography. The state-tomography version provides rigorous non-asymptotic confidence regions in trace distance and, for the uniform POVM, saturates known lower bounds up to constants [1809.11162]. Comparative analysis of quantum tomography estimators further shows that projected linear estimators can combine strong statistical behavior with high computational efficiency, especially in low-rank regimes [1901.07991]. A distinct line of work applies constrained least squares to structured state classes such as matrix product operators, using informationally complete POVMs including SIC-POVMs and spherical \(t\)-designs and obtaining sample complexity proportional, up to logarithms, to the number of independent parameters [2410.02583]. Taken together, these developments place projected least squares POVM tomography within a larger family of projection-based quantum estimation methods in which closed-form linear inversion, concentration inequalities, and convex projection are the defining ingredients.

Source: https://www.emergentmind.com/topics/projected-least-squares-povm-tomography