---
title: Quantum Tomography Procedure Overview
url: https://www.emergentmind.com/topics/quantum-tomography-procedure
type: topic
---

# Quantum Tomography Procedure Overview

Quantum tomography procedure denotes the family of experimental and computational workflows used to reconstruct an unknown quantum object from measurement data. In the literature considered here, the reconstructed object may be a density operator, a CPTP map, a POVM or detector, a temporal input–output map with memory, a multi-time process tensor, or a continuous quantum field. Across these variants, the recurring structure is to choose a representation of the target object, design an informationally sufficient measurement scheme, acquire data under repeated preparations or controlled interventions, and reconstruct a physical estimator by linear inversion, maximum-likelihood methods, evidence-based Bayesian regularization, or constrained projection [1404.7830][1803.08245][2103.13973][2106.11722][1406.3632].

## 1. Foundational formalism

Quantum state tomography (QST) reconstructs an unknown density operator $\rho$ of a $d$-dimensional quantum system from measurement data. Quantum process tomography (QPT) reconstructs an unknown quantum process $\mathcal{E}$ acting on density operators. In the standard formulation, a quantum process is modeled as a completely positive, trace-preserving map $\mathcal{E}:\rho \mapsto \mathcal{E}(\rho)$, and standard QPT proceeds by preparing a complete set of $d^2$ linearly independent input states, applying the process to each input, and performing QST on each output [1404.7830].

Two operator representations dominate the procedure. In Kraus form,
$$
\mathcal{E}(\rho)=\sum_k K_k \rho K_k^\dagger,\qquad \sum_k K_k^\dagger K_k = I.
$$
In a Hilbert–Schmidt orthonormal operator basis $\{E_m\}$, the same channel is written as
$$
\mathcal{E}(\rho)=\sum_{m,n}\chi_{mn} E_m \rho E_n^\dagger,
$$
with $\chi$ positive semidefinite. The Choi–Jamiołkowski isomorphism packages the process into a bipartite state
$$
J_{\mathcal{E}}=(I\otimes \mathcal{E})(|\Phi^+\rangle\langle\Phi^+|),
$$
so that reconstructing $J_{\mathcal{E}}$ is equivalent to reconstructing the channel. For state tomography and detector tomography, the basic measurement model is the Born rule, $\langle O_i\rangle_\rho=\operatorname{Tr}(\rho O_i)$ or $p_i=\operatorname{Tr}(E_i\rho)$, with POVM elements constrained by positivity and normalization [1004.0676][2406.11646].

This formalism extends naturally beyond memoryless channels. Temporal tomography targets a time-dependent input–output law for a device that processes a stream of quantum states with memory, rather than a single time-local channel [2103.13973]. Process tensor tomography reconstructs a $k$-step comb $\,\Upsilon_{k:0}\,$ that maps an entire sequence of interventions to a final state, thereby capturing non-Markovian dynamics over multiple times [2106.11722]. In continuous systems, cMPS tomography replaces finite-dimensional operator expansions by a field-theoretic ansatz
$$
|\Psi(Q,R)\rangle = \operatorname{Tr}_A\!\left\{\mathcal{P}\exp\!\left[\int_0^L dx\,(Q\otimes \mathbb{I}+R\otimes \psi^\dagger(x))\right]\right\}|\Omega\rangle,
$$
and reconstructs the state from low-order correlators [1406.3632].

## 2. Standard workflow and resource scaling

The canonical tomography workflow has four stages: state or probe preparation, application of the unknown process or measurement, acquisition of outcome statistics, and inversion or constrained estimation. For standard QPT this means $d^2$ input preparations and a QST subroutine on every output state. Because QST requires measurements of non-commuting observables, the number of experimental configurations grows rapidly with system size [1404.7830].

The scaling is explicit in several platforms. In the $n$-qubit NMR setting, with $d=2^n$ and $N=2^n$, a single QST can be implemented with approximately $\lceil N/n\rceil$ independent NMR measurements, so total QPT measurements scale as
$$
M_{\mathrm{QPT}} \approx N^2\lceil N/n\rceil.
$$
The stated examples are $\approx 8$ measurements for a single qubit, $\approx 32$ for two qubits, and $\approx 192$ for three qubits [1404.7830]. In polarization photonics, an overcomplete six-state tomography uses $6^n$ settings, whereas a dual-output three-base analyzer uses $3^n$ settings for state tomography; process tomography still requires individual input preparations [1809.07521].

The same resource issue appears in spin and solid-state platforms. In NV-center spin tomography, standard spin-qubit QST commonly measures transverse magnetization as a function of time and therefore requires many readout points per configuration. The alternative NV protocol based on time-independent observables and unitary pre-rotations reduces the measurements to the minimal informationally complete set, $d^2-1$ for a $d$-dimensional system, with one measurement per operator coefficient [2108.13738]. In CQED microwave tomography, the raw observable is not a quadrature directly but the number $m(\phi)$ of atoms detected in the upper working state in a series of $n$ atoms; quadrature information is extracted statistically from repeated series and then deconvolved from the instrumental function [1510.03155].

A related operational bottleneck is reconfiguration time rather than informational completeness. In photonic tomography, the sequential arrangement of wave-plate settings can be optimized as a traveling salesman problem. For the six-state scheme, the speedup factor grows monotonically with the number of qubits and is approximately $2$ from $n\ge 3$; experimentally, full quantum process tomography of a three-qubit controlled-controlled-phase gate was reduced from approximately $23$ hours to less than $11$ hours without changing the estimator [1809.07521].

## 3. Measurement compression, ancillas, and direct parameter access

Several procedures reduce the number of settings by trading them for ancillas, weak couplings, phase retrieval, or tailored unitary mappings.

| Procedure | Core mechanism | Stated resource property |
|---|---|---|
| Standard QPT | $d^2$ inputs plus QST on each output | $M_{\mathrm{QPT}} \approx N^2\lceil N/n\rceil$ in $n$-qubit NMR |
| AAPT | Encode the channel into a joint system–ancilla state | $M_{\mathrm{AAPT}} \approx \lceil N^2/(2n)\rceil$ |
| SSPT | Combine AAPT with AAQST | $M_{\mathrm{SSPT}}=1$ |
| Direct weak-measurement QPT | Two weak measurements and one post-selection | Each parameter is determined from only five experimental values |
| FQPT | Near-field/far-field distributions plus phase retrieval | 7 measurements total; 5 minimal-in-principle |

The ancilla-assisted route is built around the Choi state. In AAPT, the system and an ancilla of the same dimension are prepared in a maximally entangled state, the process acts on the system only, and a single QST of the joint state yields $J_{\mathcal{E}}$. AAQST further maps non-commuting observables of a target density matrix onto a commuting set in a larger Hilbert space, allowing all required expectations to be acquired in a single collective measurement. SSPT combines the two: one state preparation, one application of $\mathcal{E}$, and one collective measurement of commuting observables on system plus ancillas produce complete process characterization [1404.7830].

The direct weak-measurement scheme replaces full output tomography by a parameter-local relation. For each tuple $(i_1,i_2,i_3,i_4)$, the process parameter $\chi_{i_1 i_2 i_3 i_4}$ is obtained from two weak measurements, one post-selection, four joint pointer correlations $r_1,\dots,r_4$, and the post-selection probability $p_{f|i\hat{A}\hat{B}}$. The key intermediate quantity is
$$
X^{i\hat{A}\hat{B}f}:=\operatorname{tr}\!\left[\Pi_f\,\hat{B}\,\mathcal{E}(\hat{A}\rho_i)\right],
$$
and the paper’s central claim is that each parameter of the process is directly determined from only five experimental values. The number of experimental setups is $d_{\mathrm{in}}$, because one input setting yields all $d_{\mathrm{in}}d_{\mathrm{out}}^2$ parameters associated with that input in parallel [1309.5780].

A different compression mechanism appears in single-qubit reaped tomography. There the system is coupled to a single-qubit pointer initialized in $|+\rangle$, the system is measured in a computational basis $X$, and standard qubit tomography on the pointer yields the data needed for reconstruction. The scheme requires the measurement of only three observables acting jointly on system and pointer, regardless of system size, and the exact pure-state reconstruction reduces to a linear system for amplitude ratios; for finite data and mixed states, the paper gives an iterative maximum-likelihood algorithm that scales via MPS/MPO structure [2206.09562].

Fourier Quantum Process Tomography achieves compression in a different regime. For a space-dependent SU(2) process, it measures probability distributions in two Fourier-conjugate spaces and uses phase retrieval to recover the local unitary. The practically recommended set is 7 measurements total: 3 near-field distributions, 3 far-field distributions, and 1 far-field distribution without polarization projection. The paper reports this as a near-minimal set, benchmarked against standard maximum-likelihood QPT [2312.13458].

## 4. Temporal, dynamical, and non-Markovian generalizations

When the device has memory, the object to be reconstructed is no longer a single channel. Three distinct generalizations appear in the cited literature.

The first is dynamical tomography with known evolution and a fixed measurement setup. The procedure subjects the unknown state to a known time evolution for different times and performs the same measurement each time. For a suitable unitary dynamics on an $n$-dimensional system, any two states can be discriminated by a measurement with $n$ outcomes at $n+1$ points in time. If prior information restricts the state to a semi-algebraic subset $\mathcal{R}$, then $l$ time steps suffice when $(n-1)l \ge 2\dim \mathcal{R}$; beyond unitary dynamics, feasible CPTP dynamics reduce the measurement setup to the minimal number of two outcomes [1605.06786].

The second is temporal quantum tomography via quantum reservoir computing. Here the target is a temporal map acting on a stream of quantum states with finite memory. The reservoir obeys
$$
\rho_n=\Lambda_{\beta_n}(\rho_{n-1})=\operatorname{Tr}_A\!\left[U(\rho_{n-1}\otimes \beta_n)U^\dagger\right],
$$
local observables define a feature vector $x_n$, and a linear readout $\sigma_n=W^T x_n$ reconstructs the output density matrix. The method is explicitly approximate, but in the reported benchmark tasks it achieved RMSF above $98\%$ with $N_m$ up to $5$ qubits, $M$ up to $5$, training lengths $L\in[500,3000]$, and evaluation horizons $T\in[200,1000]$ [2103.13973].

The third is process tensor tomography. Rather than infer an effective recurrence, it reconstructs the full multi-time object $\Upsilon_{k:0}$ satisfying positivity and recursive causality constraints such as
$$
\operatorname{Tr}_{o_k}[\Upsilon_{k:0}] = \mathbb{I}_{i_k}\otimes \Upsilon_{k-1:0}.
$$
This fully characterizes non-Markovian dynamics over a time frame and predicts outcomes for arbitrary intervention sequences within the control span. Because the setting count scales as $N^k$ for an $N$-element basis over $k$ steps, the paper also develops a finite-Markov-order variant in which the experimental scaling reduces from $O(N^k)$ to $O(k\cdot N^\ell)$ when the effective memory length is $\ell$ [2106.11722].

These extensions shift the role of tomography from static certification to dynamical diagnostics and control. A stated implication is that memory-aware reconstruction can be used to optimize subsequent circuits: on superconducting devices, full process-tensor-based predictive control improved average final-state fidelity by approximately $0.045$, with maxima around $0.10$, and finite-$\ell$ models performed better when neighbor-induced ZZ crosstalk increased effective memory [2106.11722].

## 5. Inference, physicality, and uncertainty quantification

Tomography procedures differ as much in their estimators as in their measurement designs. The supplied literature covers at least six inference strategies.

Evidence-based tomography introduces a prior expectation state $\rho_0$ and an unknown confidence hyperparameter $\alpha$, with prior
$$
p(\rho|\alpha)\propto \exp[-\alpha S(\rho\Vert \rho_0)].
$$
The hyperparameter is selected by maximizing the evidence, and the MAP estimator becomes an entropic interpolation between the MaxEnt state consistent with measured expectations and the anticipated state. The method is efficient when the measured observables are not necessarily informationally complete but are numerous enough and broadly compatible with the prior; it yields error bars and avoids zero-eigenvalue artifacts of small-sample maximum likelihood [1004.0676].

Reliable tomography in the frequentist sense replaces a single estimate by a confidence region. The Christandl–Renner construction defines a likelihood-proportional density $\mu_{B^n}(\sigma)$ over states and forms a region $\Gamma_{\mu_{B^n}}^\delta$ by taking a high-probability set in $\mu_{B^n}$ and thickening it in fidelity by
$$
\delta^2=\frac{2}{n}\left(\ln\frac{2}{\epsilon}+2\ln c_{2n,d}\right),\qquad
c_{N,d}=\binom{N+d^2-1}{d^2-1}.
$$
The resulting confidence region has coverage at least $1-\epsilon$, independently of any prior assumption on the distribution of possible states, and the construction applies to arbitrary measurements including fully coherent collective ones [1108.5329].

A more recent line identifies iterative maximum-likelihood tomography with retrodictive recovery maps. For POVM tomography, the log-likelihood gradient coincides with the Petz recovery map of the measurement channel, and the update
$$
\rho_{t+1}\propto \rho_t^{1/2} R_t \rho_t^{1/2}
$$
is exactly the symmetric R$\rho$R step. In the commuting-output case, repeated Petz updates monotonically increase the likelihood; for channel tomography, the same construction lifts to the Choi matrix and is followed by a TP retraction [2606.23777].

Physicality enforcement is treated algorithmically in boosted projective methods. Starting from a Hermitian linear-inversion estimate $\tilde Z$, the Frobenius projection onto the CPTP set is approximated by alternating exact affine projections with positivity projections and a final Cholesky-based correction. For QPT up to four qubits, the combined Dykstra + CBA method achieved up to approximately three orders of magnitude reduction in Frobenius error to the exact SDP projection relative to earlier projective methods, at comparable runtime [2406.11646].

Two other inferential issues recur. First, SPAM bias: standard QPT assumes known state preparation and measurement, and the supplied literature includes an explicit correction procedure based on additional SPAM-calibration data, a gauge-regularized factorization of the SPAM error matrix $\hat E$, and corrected estimators such as
$$
\hat G = \hat E^{-1/2}\hat G_0 \hat E^{-1/2}.
$$
This correction applies to standard and overcomplete linear-inversion QPT and is designed to be extendable to MLE with CPTP constraints [2412.16293]. Second, incomplete identifiability: joint state-and-measurement tomography under imperfect measurements uses alternating maximum-likelihood estimation over unknown states and POVM parameters, and when the engineered measurements are not informationally complete it replaces point identification by SDP bounds on expectation values over the convex set of compatible states [1803.08245].

## 6. Platform realizations and operational limits

The procedures above are not tied to a single hardware class, and the cited papers document implementations or numerical protocols across NMR, trapped ions, superconducting devices, photonics, NV centers, CQED microwave cavities, and ultracold atoms.

In liquid-state NMR, SSPT was demonstrated on a three-qubit register in iodotrifluoroethylene on a $500$ MHz spectrometer at $300$ K. The system qubit and two ancillas enabled one-shot characterization of several single-qubit gates, with gate fidelities approximately $0.99$ for NOP, NOT-X, NOT-Y, and Phase-$\pi$, $0.95$ for Hadamard, and $0.97$ for Phase-$\pi/4$; each full single-qubit SSPT experiment took less than $4$ s from thermal equilibrium [1404.7830]. The same paper used SSPT to characterize a twirling channel and found measured $|\chi_{EE}|$ and $|\chi_{ZZ}|$ curves tracking theory with average gate fidelity greater than $0.96$ and errors below $8\%$.

In NV-center tomography, the time-independent fluorescence observable
$$
O=P_0\otimes E_{\text{nuc}}=\tfrac12(E+Z)\otimes E
$$
is combined with unitary pre-rotations that map desired density-matrix components onto electron-spin population. The paper reports 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$ for prepared electron–nuclear states, using optimized control sequences of duration at most $15\,\mu\text{s}$ [2108.13738].

In photonic platforms, two distinct operational themes appear. One is measurement-order optimization, where full three-qubit process tomography of a controlled-controlled-phase gate was shortened from approximately $23$ hours to less than $11$ hours by reordering settings alone [1809.07521]. The other is FQPT, experimentally tested on space-dependent polarization transformations. In a 1D LCMS example, FQPT achieved average unitary overlap fidelity $\bar F \approx 95.7\%$ with runtime about $1$ minute, versus approximately $30$ minutes for the standard maximum-likelihood comparison; in a 2D cascaded-LCMS example, the reported average overlap fidelity was approximately $91.4\%$ [2312.13458].

Ultracold-atom field tomography uses matter-wave interferometry, even-order phase correlators, and a low-bond-dimension cMPS ansatz. For quenched quasi-condensates, a $D=2$ reconstruction at $3$ ms achieved mean relative deviations $\epsilon(C2)\approx 0.3\%$, $\epsilon(C4)\approx 1.4\%$, and $\epsilon(C6)\approx 2.1\%$; the degradation at $7$ ms and $23$ ms was interpreted as consistent with entanglement growth and the breakdown of a low-$D$ cMPS description [1406.3632]. CQED microwave tomography is formulated differently: it reconstructs the quadrature tomogram from atomic detection statistics, but the measured result is the true quadrature distribution convolved with an instrumental Gaussian of variance $\sigma_s=1/(n v^2)$; numerical simulations with $\lambda\tau=0.04$, $n=300$ or $1000$, and fitted $u$ values $0.76$ or $0.36$ show the protocol’s effectiveness [1510.03155].

Across platforms, the limitations recur with striking regularity. The cited papers name pulse imperfections, RF inhomogeneity, decoherence during the tomography sequence, spectrometer drift, imperfect gradients, finite shots, measurement back-action, ill-conditioned linear inversion, non-identifiability under incomplete measurements, model mismatch, ancilla overhead, spectral crowding, and the need for well-calibrated controls as the dominant obstacles [1404.7830][2103.13973][1803.08245][2406.11646][2412.16293]. This suggests a unifying operational principle: quantum tomography procedures are most effective when the measurement design, estimator, and platform constraints are co-designed rather than treated as separable stages.

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