---
title: Sparse Pure-State Tomography Overview
url: https://www.emergentmind.com/topics/sparse-pure-state-tomography
type: topic
---

# Sparse Pure-State Tomography Overview

Searching arXiv for the cited works to ground the article and verify metadata.
Sparse pure-state tomography denotes a family of quantum state tomography protocols that exploit structural promises stronger than generic low-noise reconstruction. In the broadest usage found in the literature, the target is assumed to be pure or nearly pure, so that the relevant model class is low rank; in narrower usages, the state is promised to have only a small number of nonzero computational-basis amplitudes, or the measurement corruption is modeled by a sparse outlier term. For an \(n\)-qubit system with \(d=2^n\), full tomography needs \(4^n-1\) real parameters, whereas a normalized pure state depends on \(2d-2\) real parameters, and this gap motivates measurement-minimal, compressed-sensing, matrix-completion, and circuit-based schemes [1601.05379].

## 1. Conceptual scope and uniqueness notions

A central distinction is between generic mixed-state tomography and tomography restricted to pure states. In the pure-state setting, one measures a collection of Hermitian observables \(\mathbf A=\{A_0=I,A_1,\dots,A_{m-1}\}\), producing outcomes
\[
\boldsymbol{\alpha}=\bigl(\Tr(\rho A_0),\Tr(\rho A_1),\dots,\Tr(\rho A_{m-1})\bigr).
\]
Two uniqueness notions are standard. A pure state \(|\psi\rangle\) is *uniquely determined among all states* (UDA) if no other density operator reproduces the same data, and it is *uniquely determined among pure states* (UDP) if no different pure state does so. If \(\mathcal S(\mathbf A)\) is the real span of the measured observables and \(\mathcal S(\mathbf A)^\perp\) its orthogonal complement in the Hermitian operator space, then \(\mathbf A\) is UDA for every pure state if and only if every nonzero \(H\in\mathcal S(\mathbf A)^\perp\) has at least two positive and two negative eigenvalues [1601.05379].

This literature uses “sparse” in several technically distinct ways. One line treats purity itself as a rank-1 sparsity prior on the density matrix. Another assumes that only \(k\) computational-basis amplitudes are nonzero. A third augments a low-rank density matrix with a sparse outlier matrix in the measurement model. This suggests that “sparse pure-state tomography” is not a single algorithmic paradigm but a cluster of related structural-recovery problems [1401.6533].

| Approach | Structural assumption | Key claim |
|---|---|---|
| ADMM compressive tomography | \(\rho\) low rank, \(S\) sparse outliers | Recover pure or near pure states corrupted by sparse noise [1401.6533] |
| Pauli UDA design | Arbitrary pure \(n\)-qubit states | 11 Pauli measurements for two qubits; 31 for three qubits [1601.05379] |
| MC-QST | Rank-1 density matrix | \(2n+1\) local Pauli settings and direct matrix filling [2111.11071] |
| Sparse-entry circuit tomography | \(k\) nonzero amplitudes | Identity plus at most \(2(k-1)\) further measurements [2407.20298] |
| Phase-estimation tomography | Unknown \(K\)-sparse superposition | Support finding via phase estimation on a designed unitary [2111.04359] |
| Two-basis algebraic tomography | \(s\)-sparse algebraic pure state | Two measurement bases suffice asymptotically [2501.17061] |

## 2. Convex compressed-sensing formulations and robust ADMM

In the low-rank-plus-outlier formulation of Li and Cong, one seeks a density matrix \(\rho\) and a sparse matrix \(S\) that jointly explain the measurements:
\[
\begin{aligned}
&\min_{\rho,S}\;\|\rho\|_* \;+\;\|S\|_1 \\
&\;\text{s.t.}\quad \|\,y - A\,\mathrm{vec}(\rho+S)\|_2^2 \le \varepsilon,
\quad \rho=\rho^*,\;\rho\succeq0~.
\end{aligned}
\]
Here \(\|\rho\|_*\) is the nuclear norm, \(\|S\|_1\) is the entry-wise \(\ell_1\) norm, and \(\mathrm{tr}(\rho)=1\) may be added if desired. The scaled augmented Lagrangian introduces a dual variable \(u\in\mathbb R^m\),
\[
L_{\lambda}(\rho,S,u)
= \|\rho\|_* + I_{\mathcal C}(\rho)
+ \|S\|_1
+ \tfrac{ \lambda}{2}\big\|A\,\mathrm{vec}(\rho+S)-y+u\big\|_2^2,
\]
with \(\mathcal C=\{\rho=\rho^*,\,\rho\succeq0\}\), and then splits the updates over \(\rho\) and \(S\). The \(\rho\)-step uses a least-squares solve followed by projection onto \(\mathcal C\) with low-rank shrinkage, while the \(S\)-step uses a least-squares solve followed by entry-wise soft thresholding. The dual update is
\[
u^{k+1}=u^k+\bigl[A\,\mathrm{vec}(\rho^{k+1}+S^{k+1})-y\bigr].
\]
If the true rank \(r\) is known, the algorithm keeps only the top \(r\) positive eigenvalues in the shrinkage step [1401.6533].

The stopping criteria monitor the primal residual \(\|A\,\mathrm{vec}(\rho+S)-y\|_2\) and dual residuals \(\|\rho^k-\rho^{k-1}\|_2+\|S^k-S^{k-1}\|_2\) against tolerances \(\varepsilon_p,\varepsilon_d\). Typical \(\varepsilon_p\approx10^{-3}\), \(\varepsilon_d\approx10^{-4}\), and 30–50 iterations suffice. The per-iteration cost is dominated by a least-squares solve and a \(d\times d\) eigendecomposition in \(O(d^3)\), with \(d=2^q\) implying eigendecomposition cost \(O(2^{3q})\). The compressive advantage appears when \(m\approx O(r d\log d)\ll d^2\) [1401.6533].

The numerical findings reported for \(q=5\) qubits and rank \(r=2\) are specific. For relative error \(\|\rho^*-\hat\rho\|_2/\|\rho^*\|_2\), least squares requires \(m/d^2\approx0.9\) for 10% error, the Dantzig method needs \(\approx0.5\), and ADMM needs \(\approx0.35\). With sparse outliers affecting 1% of entries, Dantzig fails with error approaching 1, whereas ADMM still recovers with error \(<0.2\) using \(m/d^2\approx0.5\). The stated interpretation is that the \(\ell_1\) penalty on \(S\) isolates heavily corrupted entries while preserving the low-rank \(\rho\) [1401.6533].

## 3. Minimal Pauli measurement sets and pure-state identifiability

For Pauli tomography, the problem is to find a minimal subset of \(n\)-qubit Pauli operators that is UDA for all pure states. On two qubits, full tomography uses 15 nontrivial Pauli operators, while Ma *et al.* showed that the 11-element set
\[
\mathbf P_2=\{II,IX,IY,IZ,XI,YX,YY,YZ,ZX,ZY,ZZ\}
\]
suffices to UDA any two-qubit pure state, and no smaller Pauli set can do so for all two-qubit pure states. On three qubits, full tomography uses 63 nontrivial Paulis, and a set of 31 Paulis suffices to UDA every three-qubit pure state; again, one cannot use fewer than 31 Paulis for UDA of all pure states [1601.05379].

Known dimension-dependent bounds are sharper in the non-Pauli setting than in the Pauli setting. For a \(d\)-dimensional pure state, UDP requires \(m\ge 4d-4\) observables, whereas UDA requires \(m\ge 5d-6\) observables. For Pauli measurements on general \(n\)-qubit systems, no closed form is known. Compressed sensing shows that \(O(d\log d)=O(2^n n)\) random Paulis suffice to identify almost all pure states, but this fails on a measure-zero set. The conjecture reported in the three-qubit work is that the minimal Pauli count scales linearly in \(d=2^n\), with \(m_{\min}\approx c\,2^n\) and \(c\approx1.7\) for \(n=2,3\); a constructive algorithm for general \(n\) remains open [1601.05379].

The same work also analyzes robustness. Under the depolarizing channel
\[
\rho(\eta)=(1-\eta)|\psi\rangle\langle\psi|+\eta\,\frac{I}{d},
\]
simulations on 100 Haar-random pure states show average fidelity remaining \(\gtrsim0.98\) for \(\eta\lesssim0.05\), with roughly linear decay as \(\eta\) increases. In liquid-state NMR experiments, the two-qubit sparse protocol achieved typical fidelity \(F\approx0.992\) relative to full tomography, and the three-qubit GHZ experiment achieved \(F\approx0.96\). Random Pauli subsets of the same size performed worse: mean fidelity \(\approx0.90\) for 11 random Paulis in the two-qubit case and \(\approx0.87\) for 31 random Paulis in the three-qubit case, which the paper attributes to “failing sets” whose complements admit Hermitians with the wrong eigenvalue signature [1601.05379].

## 4. Matrix completion from local Pauli data

Matrix-Completion Quantum State Tomography (MC-QST) specializes pure-state tomography to the rank-1 case \(\rho=|\psi\rangle\langle\psi|\). The algebraic basis is that all \(2\times2\) minors of a rank-1 matrix vanish:
\[
\det\!\begin{pmatrix} \rho_{i,r}&\rho_{i,s}\\ \rho_{j,r}&\rho_{j,s} \end{pmatrix} =0
\quad\Longrightarrow\quad
\rho_{i,s}\rho_{j,r}=\rho_{i,r}\rho_{j,s}.
\]
The measurement scheme uses only \(2n+1\) local settings: one computational-basis measurement and, for each qubit \(i\), single-qubit measurements of
\[
\Lambda_1^i=I^{\otimes(i-1)}\otimes \sigma_x\otimes I^{\otimes(n-i)},
\qquad
\Lambda_2^i=I^{\otimes(i-1)}\otimes \sigma_y\otimes I^{\otimes(n-i)}.
\]
The computational basis yields all \(d\) diagonal entries \(\rho_{i,i}\), and the \(X\) and \(Y\) measurements recover off-diagonal entries between basis states that differ in exactly one bit. The directly measured entries therefore total
\[
d+n\frac d2=d\left(1+\frac n2\right)>2d-1,
\]
which matches the spanning-tree style criterion quoted from Strang: any \(2d-1\) entries forming an acyclic row-column graph uniquely determine the full rank-1 matrix [2111.11071].

Reconstruction proceeds in three steps. First, estimate the diagonal entries from the computational basis. Second, recover nearest-neighbor off-diagonals from the \(X\) and \(Y\) data through
\[
\hat\rho_{j,k}
=\tfrac12\bigl(-p_{1,-}^{j,k}+p_{1,+}^{j,k}
-\mathrm i\,p_{2,-}^{j,k}
+\mathrm i\,p_{2,+}^{j,k}\bigr).
\]
Third, fill in all remaining entries by rank-1 completion:
\[
\hat\rho_{i,j}=\hat\rho_{i,r}\hat\rho_{r,j}/\hat\rho_{r,r}.
\]
Finally, renormalize by \(\hat\rho\gets \hat\rho/\mathrm{Tr}\hat\rho\) and extract the leading eigenvector. Because each fill-in is \(O(1)\) and there are \(d^2\) entries, the classical post-processing cost is \(O(d^2)\) [2111.11071].

The stated scaling is \(2n+1=2\log_2 d+1=O(\log d)\) measurement settings, \(O(d\log d)\) estimated probabilities, and \(N=O(d\log d/\epsilon^2)\) total samples to reach mean-square error \(\epsilon^2\) on each probability. The method is exact in the noiseless rank-1 limit but does not provide a fully worked-out finite-sample error bound. Numerically, it outperformed contemporary pure-state tomography baselines on Haar-random states for \(n=2,3,4,5\), and on IBM “ibmq-casablanca” it reconstructed two-qubit pure states with median infidelities below 5%. For a three-qubit GHZ state, after a prerotation \(V=H^{\otimes3}\), the reported median fidelity was 97%. The principal limitation stated in the paper is that the method assumes exact purity; if \(\rho\) is slightly mixed, the minors do not vanish exactly and regularization becomes necessary [2111.11071].

## 5. Computational-basis sparsity: circuit constructions and phase estimation

A different branch of sparse pure-state tomography assumes that only \(k\) computational-basis amplitudes are nonzero. In the circuit-based sparse-entry protocol, one first performs the identity measurement in the computational basis to learn the support set \(S=\{x:a_x\neq0\}\) and the magnitudes \(|a_x|^2\). One then performs \(2m\) additional measurements \(\{U_1|\psi\rangle,\dots,U_{2m}|\psi\rangle\}\) with \(m\le k\), where each \(U_j\) consists of a short CNOT network \(R_j\) mapping a pair of support states \(x,y\in S\) to basis states differing in exactly one bit, followed by either a Hadamard \(H\) or a phase-shifted Hadamard \(V=HD\), \(D=\mathrm{diag}(1,i)\), on the differing qubit. The support graph on \(S\) is weighted by Hamming distance, and a minimum-spanning tree determines the order of pairwise phase recovery. For an edge \((x,y)\), the probabilities \(P_+(x')=|a_x+a_y|^2\), \(P_+(y')=|a_x-a_y|^2\), \(P_{+i}(x')=|a_x+i a_y|^2\), and \(P_{+i}(y')=|a_x-i a_y|^2\) determine the real and imaginary parts of \(a_x\bar a_y\), thereby recovering the unknown amplitude once the adjacent known amplitude is fixed [2407.20298].

The resource counts are explicit. The protocol uses the identity measurement plus at most \(2(k-1)\) further measurements, so the total is \(\le 2k\). If the minimum-spanning tree has \(p_j\) edges of weight \(j\), then the total CNOT count is bounded by
\[
2\,(p_2+2p_3+\cdots +(r-1)p_r).
\]
The worst case is \(O(nk)\), but when \(k>2^n-n\) the minimum-spanning tree can be chosen entirely from weight-1 edges, so no CNOTs are needed. The single-qubit gate count is \(O(k)\), the depth of each \(U_j\) is \(O(w_j)\), and the overall depth is \(O(kn)\). Simulations in Qiskit/Aer with the IBM Brisbane noise model showed median fidelity \(>0.99\) when support states were separated by Hamming distance 1, while a randomized backup strategy using \(H^{\otimes n}\) and only single-qubit gates often outperformed the unrandomized MST construction once the MST required at least two CNOTs. The same sparse-state embedding was also applied to process tomography, with reported fidelities \(\gtrsim0.98\) after a polar-decomposition correction [2407.20298].

The phase-estimation approach to \(K\)-sparse tomography assumes no prior knowledge of \(K\), the support bit strings, or the complex coefficients. It constructs a specially designed unitary \(U_\Phi\) on \(n+1\) qubits whose eigenbasis encodes the computational labels, prepares the embedded state
\[
|\Psi_e\rangle
=\sum_{k=1}^K c_k e^{i\theta_k}\bigl(d_{k,1}|E_{s_k,1}\rangle+d_{k,2}|E_{s_k,2}\rangle\bigr),
\]
and uses standard quantum phase estimation with controlled powers \(U_\Phi^{2^j}\). Measuring the phase register yields an estimate of an eigenphase with probability \(|c_k|^2|d_{k,l}|^2\), and after applying \(H^{\otimes n}\) to the data register, a computational-basis measurement reveals the support string \(s_k\). Repeating this \(O(2/|c_{\min}|^2)\) times suffices, with high probability, to observe all \(K\) distinct supports, where \(c_{\min}\) is the least amplitude magnitude appearing in the superposition. Once the support is known, the coefficient vector is reconstructed using conventional compressive-sensing tomography on the \(K\)-dimensional reduced problem, requiring \(O(d\,K(\log K)^c)\) settings with \(c\in[2,6]\) and \(d\) a constant, independent of the ambient dimension \(2^n\). The paper emphasizes open problems concerning the existence of the required phase choices \(\{\phi_i\}\) and the efficient implementation of exponentially large powers \(U_\Phi^{2^j}\) [2111.04359].

## 6. Two-basis algebraic tomography, tensor-network parameterizations, and current open directions

The two-basis protocol of Feng *et al.* works with *\(s\)-sparse algebraic* pure states
\[
|\psi\rangle=\sum_{k=0}^{d-1}\gamma_k|k\rangle,\qquad \gamma_k=r_k e^{i\varphi_k},
\]
where at most \(s\) amplitudes are nonzero and each nonzero amplitude is an algebraic number. The first basis is the computational basis \(\mathcal B_1=\{|k\rangle\}\), from which one reads off \(P_k=|\langle k|\psi\rangle|^2=r_k^2\), the magnitudes \(r_k=\sqrt{P_k}\), and the support \(V=\{k:r_k\neq0\}\). The second basis is a shifted-Fourier basis \(\mathcal B_2=\{S'|k\rangle\}\) with \(S'=DF\), where \(F\) is the quantum Fourier transform and \(D=\mathrm{diag}(e^{i\theta_0},\dots,e^{i\theta_{d-1}})\) is chosen so that the differences \(\theta_s-\theta_\ell\) are distinct algebraic numbers in absolute value. The observed probabilities are
\[
Q_j
=\frac1d
+\frac2d
\sum_{\substack{s<\ell\\ s,\ell\in V}}
r_s r_\ell
\cos\!\Bigl[(\varphi_s-\varphi_\ell)-(\theta_s-\theta_\ell)-\tfrac{2\pi j(s-\ell)}d\Bigr].
\]
Reconstruction fixes \(\varphi_0=0\) and minimizes \(\|\tilde Q(\varphi)-Q_{\rm meas}\|_2\) over the remaining phases, for example by simulated annealing or another nonlinear solver. The uniqueness proof uses the Lindemann–Weierstrass theorem to show that the distinct algebraic phase differences force all coefficients in the induced trigonometric relations to vanish, so the relative phases are uniquely determined. The scheme is asymptotically informationally complete for pure states because every pure state can be approximated arbitrarily well by algebraic states, and for GHZ-like and W-like states it admits two local product bases using \((RH)^{\otimes n}\). Numerical tests with up to 20 qubits showed reconstructed-state fidelity \(>99.9\%\) at small residual probability error and \(\gtrsim95\%\) fidelity under depolarizing noise up to \(p=0.1\) [2501.17061].

A complementary large-scale direction uses tensor-train parameterizations. In Tensor Train Quantum State Tomography using Compressed Sensing, the pure state \(|\psi\rangle\in\mathbb C^{d^N}\) is written as a TT tensor with cores \(G_n\in\mathbb C^{r_{n-1}\times d\times r_n}\),
\[
\psi_{i_1\ldots i_N}
=
G_1(i_1)\cdots G_N(i_N),
\]
and one solves
\[
\min_{G_1,\dots,G_N}\;\frac12\bigl\|y-\mathcal A(|\psi(G)\rangle\langle\psi(G)|)\bigr\|_2^2
\quad\text{s.t.}\quad \|\psi(G)\|_2=1.
\]
The proposed algorithm is projected gradient descent in TT form with TT-SVD truncation and reorthogonalization. For pure-state tomography, the stated per-measurement complexity is \(O(N\,r^3\log d)\), and storage is \(O(N\,d\,r^2)\). The paper points to existing TT iterative hard-thresholding guarantees of \(O(r^2 d N\log N)\) random Gaussian measurements for low-TT-rank tensors and conjectures comparable Pauli-measurement analogues. This is not a sparsity model in the computational basis; it is a low-TT-rank model that includes pure states, nearly pure states, and ground states of Hamiltonians [2506.23560].

Recent optimization work has also shifted attention from reconstruction itself to certifying when a reduced measurement set uniquely specifies a pure state. The ALM framework for investigating pure-state uniqueness formulates UDP as a constrained nonconvex optimization over pure states and UDA as either an SDP or a low-rank constrained problem. A key theorem states that, given \(m\) non-identity measurements, any compatible density operator can be replaced by another compatible operator of rank \(<\sqrt{m+2}\); with \(|\psi\rangle\langle\psi|\) added as an extra observable, the rank bound becomes \(<\sqrt{m+3}\). Using the ensemble parameterization
\[
\rho=\sum_{i=1}^r p_i|\phi_i\rangle\langle\phi_i|,
\qquad r=\lfloor\sqrt{m+3}\rfloor,
\]
and ALM updates
\[
\lambda^{(k+1)}=\lambda^{(k)}+\mu^{(k)}g^{(k)},\qquad
\mu^{(k+1)}=\min\{\gamma\mu^{(k)},\mu_{\max}\},
\]
the method classified qutrit and four-qubit symmetric states into three categories: UDA, UDP but not UDA, and neither. This suggests that the design of sparse measurement sets is inseparable from the geometry of uniqueness, not merely from sample complexity [2501.00327].

Several open questions recur across these approaches. No closed-form minimal Pauli count is known for arbitrary \(n\)-qubit pure states, and a scalable constructive algorithm remains open. MC-QST assumes exact rank 1 and must be stabilized for mixedness. The phase-estimation protocol leaves unresolved the structured implementation of \(U_\Phi^{2^j}\). The two-basis algebraic scheme is asymptotic for generic pure states rather than exact in finite dimension. TT-based methods trade exact convexity for nonconvex low-rank optimization. Taken together, these lines of work define sparse pure-state tomography as a technically diverse program: reducing measurement count, circuit depth, and classical complexity by matching the tomography model to the structure actually present in the quantum state.

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