---
title: Quantum State Tomography
url: https://www.emergentmind.com/topics/quantum-state-tomography
type: topic
---

# Quantum State Tomography

Quantum state tomography (QST) is the task of inferring the state of a quantum system by appropriate measurements, or, equivalently, reconstructing an unknown quantum state \(\rho\) from measured data. In the standard formulation, QST reconstructs the density matrix by estimating measurement probabilities such as \(y_m = \operatorname{Tr}(\rho E_m)\) or \(P(\boldsymbol{a}) = \mathrm{Tr}(M^{\boldsymbol{a}}\rho)\), depending on the measurement model. It is described as the gold standard for verification and benchmarking of quantum devices, but its central difficulty is the exponential growth of parameters and measurement resources with system size: for an \(N\)-qubit system, the density matrix is of size \(2^N \times 2^N\); full QST may require \(4^n-1\) observables, \(3^N\) Pauli measurement settings, or, in the most direct count, \(4^n\) independent density-matrix entries [2506.23560] [1604.02046] [2206.00019] [2401.12864].

## 1. Core formulation and the scalability problem

In its conventional form, QST estimates a density operator from repeated measurements on identically prepared copies. Several equivalent formulations appear across the literature. In compressed-sensing QST, one solves
\[
\min_{\hat{\rho}} \|\hat{\rho}\|_* \quad \text{s.t.} \quad \|\vec{y} - \mathcal{M}(\hat{\rho})\|_2 \leq \epsilon,\; \hat{\rho} \in \mathcal{S},
\]
where \(\hat{\rho}\) is constrained to Hermitian, PSD, and unit-trace matrices [2506.23560]. In informationally complete measurement schemes, the outcome distribution itself determines the state through relations of the form
\[
P(\boldsymbol{a}) \equiv \mathrm{Tr}(M^{\boldsymbol{a}}\rho),
\]
so that learning the measurement distribution is equivalent to learning \(\rho\) [2212.04940].

The principal obstruction is dimensional. Standard QST becomes impractical because the number of quantum measurements and the amount of computation required to process them grows exponentially in the system size, and both memory and computational bottlenecks become severe even for moderately sized systems [1101.4366] [2506.23560]. This is why a large part of the modern literature does not treat tomography as a single protocol, but as a family of protocols adapted to structure: low rank, matrix product structure, tensor-network compressibility, locality, restricted observability, continuous measurement records, or fixed-shot property estimation.

A further aspect of the basic formulation is physicality. Density matrices must be positive semidefinite and trace one, and many practical reconstruction methods are designed to enforce these constraints by construction rather than by a posteriori correction. Cholesky-type parameterizations, lower-triangular factorizations, and tensor-network factorizations are all used for this purpose in different settings [1702.03635] [2506.23560].

## 2. Informational completeness, locality, and identifiability

A central question in QST is not only how to reconstruct a state, but whether the measured data determine it uniquely. Full informational completeness is the strongest setting, but many efficient schemes rely on reduced or restricted data. One prominent example is tomography via local reduced density matrices (RDMs). Measuring all weight-1 and weight-2 Pauli operators can be much more manageable than full tomography, but this only works when the global state is uniquely determined by its local RDMs [1604.02046].

That work distinguishes two notions of uniqueness. A state is **UDA** if it is uniquely determined among all states by its RDMs, and **UDP** if it is uniquely determined among pure states only. It presents explicit 4-qubit examples of states that are UDP but not UDA for their 2-RDMs, thereby disproving the conjectured implication UDP \(\Rightarrow\) UDA. The classification leads to three classes: states that are neither UDP nor UDA, states that are both UDP and UDA, and states that are UDP but not UDA. The practical consequence is precise: for Class B, 2-RDM tomography is sufficient; for Class C, it is sufficient only if global purity is assumed; for Class A, full QST is necessary [1604.02046].

Locality also appears in tensor-network tomography. For one-dimensional systems well approximated by a matrix product state (MPS) ansatz, local reductions on blocks of \(k = \lceil \log_d(D) \rceil + 1\) consecutive sites can determine the global state generically, and only a linear number of experimental operations is required [1101.4366]. This suggests that “efficient tomography” is often inseparable from a structural identifiability statement: without a guarantee that local or compressed data determine the global object, reduced-measurement protocols become ambiguous rather than merely approximate.

Another route to informational completeness replaces many settings with a single generalized measurement. Symmetric informationally complete (SIC) POVMs are described as mathematically optimal measurements for tomography, requiring only \(d^2\) outcomes for a \(d\)-dimensional quantum system—the absolute minimal informationally complete measurement. In that setting, complete tomographic information is contained in every single experimental shot [2206.00019].

## 3. Reconstruction principles and estimators

Reconstruction procedures in QST range from direct linear inversion to constrained optimization and likelihood-based inference. Maximum-likelihood estimation (MLE) remains a standard approach because linear inversion can yield unphysical density matrices. A representative parameterization is
\[
\rho = \frac{R^\dagger R}{\operatorname{Tr}(R^\dagger R)},
\]
with \(R\) lower triangular, so positivity is enforced by construction [1702.03635]. In photonic time-bin tomography, MLE is used together with a likelihood
\[
L(R) = \sum_j \left[ \frac{(n_j^M - n_j^E)^2}{n_j^E} + \ln n_j^E \right],
\]
where \(n_j^M\) are measured counts and \(n_j^E\) are model expectations [1702.03635].

For realistic experiments, the measurement process may not be described by an instantaneous POVM. In that case, trajectory-based MaxLike tomography replaces static POVM elements by trajectory-dependent effective operators \(E^{(n)}\) computed from the adjoint quantum filter. The generalized likelihood takes the form
\[
\mathcal{L}(Y|\rho)=\prod_{n=1}^N c_n\,\mathrm{Tr}\{\rho E^{(n)}\},
\]
and confidence intervals for observables are obtained from an asymptotic expansion of multidimensional Laplace integrals [1510.01726]. This extends QST to continuous-time signals, measurement imperfections, and decoherence, while keeping the inferential object an initial quantum state.

Statistical reliability is a separate issue from point estimation. Confidence-region tomography constructs subsets of state space that contain the true state with probability at least \(1-\epsilon\), independently of any prior assumption on the distribution of possible states. The construction proceeds through a data-dependent distribution
\[
\mu_{B^n}(\sigma)d\sigma=\frac{1}{c_{B^n}}\mathrm{Tr}[B^n\sigma^{\otimes n}]\,d\sigma,
\]
followed by a “fattening” step using fidelity balls of radius \(\delta\), where
\[
\delta^2=\frac{2}{n}\left(\ln\frac{2}{\epsilon}+2\ln c_{2n,d}\right).
\]
The resulting regions are operational, prior-free, and shrink as \(1/\sqrt{n}\) [1108.5329].

Neural-network reconstructions introduce a different estimator design. In Neural-Shadow Quantum State Tomography, the training loss is the infidelity
\[
\mathcal{L}_\lambda = 1 - |\langle \psi_\lambda | \Phi \rangle|^2,
\]
estimated from classical shadows rather than basis-dependent cross-entropy. The stated motivation is that cross-entropy does not correlate well with true distance, especially for phases, whereas infidelity is a natural loss for quantum states [2305.01078].

## 4. Structured-state tomography: low rank, tensor networks, and matrix completion

Many efficient QST protocols are efficient only on restricted state classes. Compressed sensing assumes that the unknown state is low rank or close to pure. In one formulation, only \(\mathcal{O}(RD\log^2 D)\) measurements are needed for rank-\(R\) states, using random Pauli measurements and convex recovery [2506.23560]. Single-observable tomography uses the same low-rank principle but compensates for a lack of informational completeness by introducing an ancilla and random linear mixing, so that the recovery problem becomes low-rank matrix recovery under physical constraints [1704.06695].

Tensor-network methods replace the full density matrix by a compact ansatz. In one-dimensional systems, states well approximated by MPS of small bond dimension admit tomography from local measurements with classical postprocessing polynomial in system size, and one of the two main schemes comes with rigorous certification without any a priori assumptions [1101.4366]. A more recent tensor-network development parameterizes the density matrix using a low-rank Block Tensor Train (Block-TT) decomposition,
\[
\rho_{\mathrm{TT}} = \mathcal{A} \bullet_n \mathcal{A}^{\mathrm{H}},
\]
where \(\mathcal{A}\) is a Block-TT tensor with a shared index of size \(K\). When TT-rank \(=R\), the parameter count scales as \(2R^2(\log_2 D - 3) + 2R^2K + 4R\), and for states well approximated by low TT-rank the effective scaling is \(\mathcal{O}(NR^2)\) rather than \(\mathcal{O}(2^{2N})\). The parameterization also inherently preserves PSD [2506.23560].

Pure-state structure allows even more specialized reconstructions. Matrix-Completion QST reconstructs an \(n\)-qubit pure state from only \(2n+1\) local Pauli measurement settings. The method uses the fact that a rank-1 density matrix can be completed from \(2d-1\) entries, together with the rule
\[
\rho_{i,j} = \frac{\rho_{r,i}\rho_{j,r}}{\rho_{r,r}},
\]
derived from vanishing \(2\times 2\) minors. The protocol measures all diagonal elements in the computational basis and selected off-diagonal elements via local \(X\) and \(Y\) basis measurements, then fills the rest algebraically [2111.11071].

Approximate tomography can also exploit sparsity in a chosen basis. Threshold Quantum State Tomography first measures all diagonal elements \(\rho_{ii}\), then measures only those off-diagonal elements \(\rho_{ij}\) for which \(\sqrt{\rho_{ii}\rho_{jj}} \geq t\), using the PSD inequality \(|\rho_{ij}| \le \sqrt{\rho_{ii}\rho_{jj}}\). With \(t=0\), the protocol reduces to standard QST; for \(t>0\), many off-diagonal elements are omitted. The paper gives a lower bound on the fidelity between the true and thresholded reconstructions:
\[
F_\text{bound}(\{\rho_{ii}\}, t) = \left( 1 - \sqrt{r(\bar{\rho}) \sum_{\{ij\} \in \mathcal{C}_t} \rho_{ii}\rho_{jj}} \right)^2.
\]
This suggests a resource-accuracy tradeoff controlled directly at the level of matrix entries [2401.12864].

The following summary organizes the structured-state approaches described in the literature.

| Approach | Structural assumption | Resource statement |
|---|---|---|
| Compressed sensing / low rank | Low rank or nearly pure state | \(\mathcal{O}(RD\log^2 D)\) measurements for rank-\(R\) states |
| MPS / tensor networks | Small bond dimension, local correlations | Linear number of experimental operations; polynomial postprocessing |
| Block-TT QST | Low TT-rank, TT-compressible density matrix | \(\mathcal{O}(NR^2)\) effective scaling |
| Matrix-completion QST | Pure state, rank-1 density matrix | \(2n+1\) local Pauli settings |
| Threshold QST | Small or negligible off-diagonals in a chosen basis | Basis-dependent measurement reduction controlled by threshold \(t\) |

## 5. Measurement design innovations

A major research direction in QST concerns measurement design: reducing the number of settings, adapting to hardware constraints, or extracting full information from unconventional observables. Several proposals replace large setting families by a single setting or a fixed small number of settings.

Single-setting SIC tomography implements a local SIC POVM on trapped ions by mapping each qubit SIC POVM to orthogonal states embedded in a higher-dimensional system, read out using repeated in-sequence detections and without ancilla qubits. Combined with classical shadows, this yields a fixed measurement basis with full tomography information in every shot and enables estimation of arbitrary polynomial functions of the density matrix orders of magnitudes faster than standard methods [2206.00019].

Single-observable tomography pursues a related objective from a different angle. A single observable is generally insufficient, so the protocol enlarges the system by an ancilla and applies random linear mixing. Recovery then uses either rank minimization heuristics or constrained least squares under Hermiticity, positivity, and trace-one constraints [1704.06695]. A later auxiliary-system method reduces the number of settings to two by entangling the target system with a quantum auxiliary system or correlating it with a probabilistic classical auxiliary system. The stated sampling complexity is \(O(d^2)\), and the method includes two purity-measurement schemes, one of which achieves measurement precision at the Heisenberg limit [2507.16532].

Measurement design can also be tailored to experimental access restrictions. In a qubit-qutrit system, QST is formulated as a numerical optimization problem over a quorum of rank-3 projectors, with the geometric quality measure
\[
\mathcal{Q} := \left|\det \left( \langle Q_i, Q_j \rangle \right)_{i,j=1}^{n^2-1}\right|^{1/2}
\]
as objective. This framework accommodates subsystem-only measurements and numerically approximates mutually unbiased subspaces in dimension six with a deviation described as irrelevant for practical applications [2012.14494].

Some measurement schemes are architecture-specific. Time-bin qudit tomography uses cascaded Mach-Zehnder interferometers with binary delays, where the number of settings scales as \(2^K \propto d\) for \(K=\lceil \log_2 d\rceil\). For \(d=4\), the implementation used 16 measurement settings and reconstructed a maximally entangled time-bin state with average fidelity \(0.950\) [1702.03635]. In NV centers, a fast protocol measures only a time-independent observable—the electron spin population—and maps off-diagonal density-matrix elements into directly measurable population differences via selective unitary operations, thereby removing the need for measurements at many evolution times [2108.13738].

There are also measurement models in which the raw data are not conventional projective outcomes. Transport-based two-qubit tomography reconstructs the density matrix of an open system from currents, current cross-correlations, and current time-derivatives within a Markovian Lindblad framework. Populations are recovered analytically from currents and noise, while coherences are encoded in current derivatives [2501.16819].

## 6. Learning-based, selective, and unconventional formulations

Recent work increasingly treats QST as a problem of statistical learning rather than explicit inversion of measurement equations. In one formulation, QST is recast as a language-modeling task. The unknown quantum state is treated as an unknown language, quantum correlations as semantic information, and IC-POVM measurement outcomes as token sequences
\[
\mathbf{b} = (\mathrm{sos}, a_1,\ldots,a_N,\mathrm{eos}),
\]
modeled autoregressively by a customized GPT-2-style decoder-only transformer. The method learns \(P(\mathbf{a})\) instead of \(\rho\) directly, and a single trained model can reconstruct a family of similar quantum states simultaneously [2212.04940].

Selective Quantum State Tomography (SQST) changes the objective more radically. Rather than reconstructing the full state first, it estimates arbitrary elements of an unknown quantum state using a fixed measurement record. For a single matrix element, the stated sample complexity is
\[
N = O(\epsilon^{-2}\log \delta^{-1}),
\]
independent of the dimension \(d\), using measurements derived from mutually unbiased bases and estimators built from phase factors \(\eta_{ij}^{(s)}\). The same fixed and dimension-independent data sample can be used to estimate mean values from a continuous class of bounded operators on demand [1909.05880].

Classical-shadow-based neural tomography occupies an intermediate position between full reconstruction and selective estimation. Neural-Shadow QST uses randomized measurements to build shadow objects \(\hat{\rho}_i\), then trains a neural-network ansatz by minimizing estimated infidelity rather than basis cross-entropy. The reported advantage is improved learning of relative phases and robustness against various types of noise without any error mitigation [2305.01078].

Another unconventional formulation is single-qubit reaped tomography. The system is coupled to a single ancillary qubit pointer via
\[
U_\mathrm{int} = \exp\left(i\theta P \otimes |1\rangle\langle 1|\right),
\]
and after measuring the system in the computational basis, the pointer state contains the reaped wavefunction information. Standard two-state tomography on the pointer then yields the system state using only three joint observables,
\[
X\otimes \sigma^z,\qquad X\otimes \sigma^x,\qquad X\otimes \sigma^y,
\]
regardless of system size. An iterative maximum-likelihood algorithm is provided for statistically incomplete data [2206.09562].

Finally, successive-measurement tomography reconstructs finite-dimensional states from sequences of two measurements: projectors \(P_k\) and \(P_\mu\) onto two bases with no mutually orthogonal vectors. The protocol measures position-position and momentum-position correlations of two von Neumann meters, which encode the real and imaginary parts of a joint quasi-probability \(W^\rho(P_\mu + P_k)(\epsilon_1)\). The density matrix is then recovered by an explicit inversion formula [1204.6094].

## 7. Experimental validation, performance criteria, and scope

Experimental studies emphasize that “tomography” is not a single benchmark but a combination of reconstruction accuracy, sample complexity, classical cost, and suitability to the underlying platform. Reported performance criteria include fidelity, trace distance, purity, Rényi entropies, and convergence of observable estimates.

Several platforms are represented in the literature. A 4-qubit NMR processor was used to test tomography via RDMs, including robustness under artificially added Gaussian noise, with reconstructed states maintaining fidelities \(>0.9\) for the reconstructible classes [1604.02046]. A trapped-ion implementation of SIC tomography demonstrated online QST in real time on an 8-qubit entangled state and entanglement analysis on a 5-qubit absolutely maximally entangled state [2206.00019]. Time-bin photonic tomography reconstructed a 4-dimensional maximally entangled state with average fidelity \(0.950 \pm 0.003\) over 15 trials [1702.03635]. An NV-center protocol reported single-qubit fidelities \(0.994\), \(0.985\), \(0.995\), and \(0.986\) for representative states, and two-qubit fidelities \(0.98\), \(0.97\), \(0.97\), and \(0.97\) [2108.13738]. Matrix-completion tomography on IBMQ-Casablanca achieved \(97\%\) fidelity for 3-qubit GHZ states [2111.11071]. Threshold QST was demonstrated up to 7 qubits, with accurate reconstruction using only 170 measurements for a 7-qubit W state, compared with \(4^7=16{,}384\) in the conventional count [2401.12864]. Transport tomography is presented as analytically complete for two qubits in open solid-state systems when currents, noise, and derivative data are available [2501.16819]. Block-TT tomography remained tractable up to \(N=12\) qubits and reported computational time for the measurement contraction step growing linearly with the number of measurements [2506.23560].

A recurrent theme across these results is that efficiency always depends on what is being asked of the protocol. Full, assumption-free, universally valid tomography remains constrained by exponential scaling. The practical literature therefore separates at least four regimes. One regime seeks full reconstruction with rigorous statistical guarantees, as in confidence-region tomography [1108.5329]. A second exploits physical structure such as low rank, locality, or tensor-network compressibility [1101.4366] [2506.23560]. A third targets restricted tasks—selected matrix elements, purities, entropies, or entanglement witnesses—using universal data samples or classical shadows [1909.05880] [2206.00019]. A fourth adapts tomography to nonstandard observables or hardware limitations, including subsystem-only access, continuous measurement records, and transport data [2012.14494] [1510.01726] [2501.16819].

This suggests that, in contemporary usage, “quantum state tomography” denotes a broader inferential program rather than a single reconstruction recipe. The common object is still the quantum state, but the operational form of the problem—full or selective, exact or approximate, direct or learned, isolated or open-system—depends on the measurement model, the accessible controls, and the structural assumptions that can be justified for the state class under study.

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