---
title: Ancilla-Assisted Quantum Process Tomography
url: https://www.emergentmind.com/topics/ancilla-assisted-quantum-process-tomography-aaqpt
type: topic
---

# Ancilla-Assisted Quantum Process Tomography

Ancilla-assisted quantum process tomography (AAQPT), often called ancilla-assisted process tomography (AAPT), is a process-identification framework in which an unknown quantum channel acts on only one subsystem of a bipartite probe state, after which ordinary state tomography is performed on the joint output. In the ideal Choi–Jamiołkowski setting, the probe is a maximally entangled state and the output is the Choi state of the channel; more generally, AAQPT allows non-maximally entangled, mixed, and even separable probe states, provided the input state is sufficiently faithful or useful for reconstruction [2206.05899], [2110.02965], [2605.19182].

## 1. Operational formulation and Choi-state representation

In the basic AAQPT setup, one prepares a bipartite state $\rho_{AB}$, lets the unknown channel $\mathcal E$ act only on system $A$, and obtains
$$
\mathcal E_A(\rho_{AB}) := (\mathcal E \otimes id_B)(\rho_{AB}).
$$
If different channels always produce different bipartite outputs, then tomography of $\mathcal E_A(\rho_{AB})$ determines $\mathcal E$ [2206.05899].

The standard state-channel representation used in AAQPT is the Choi–Jamiołkowski isomorphism. For a channel $E$ on an $n$-qubit system, the Choi matrix is written as
$$
\Lambda_{E} = (\mathds{I} \otimes E)\big(\ket{\phi^{+}}\bra{\phi^{+}}^{\otimes n}\big),
$$
with normalized Choi state
$$
\rho_{\Lambda} = \frac{1}{2^n}\Lambda,
$$
and channel action
$$
E(\rho) = \mathrm{tr}\big[(\rho^T \otimes \mathds{I}) \Lambda_E \big].
$$
This turns process tomography into state tomography on a doubled Hilbert space [2110.02965].

A closely related formulation represents the process by a process matrix $X \in \mathbb C^{d^2\times d^2}$ satisfying
$$
X \ge 0,\qquad \operatorname{Tr}_1(X)\le I_d.
$$
For trace-preserving channels, $\operatorname{Tr}_1(X)=I_d$. In the maximally entangled-input case, the Choi state obeys
$$
\rho_{\mathcal E} = (\mathcal E\otimes I)\bigl(|\Psi\rangle\langle\Psi|\bigr),\qquad X = d\,\rho_{\mathcal E},
$$
so tomography of the bipartite output directly reconstructs the process matrix [2509.05988].

AAQPT is therefore distinguished from standard quantum process tomography by replacing many input preparations with one correlated system–ancilla preparation and a joint-output tomography stage. The literature summarized here also emphasizes that a maximally entangled probe is not mandatory in practice; a known input state plus full joint measurement data can suffice, provided the state satisfies the relevant faithfulness criterion [2410.00892].

## 2. Faithfulness, invertibility, and sensitivity

The central structural notion in AAQPT is faithfulness. A bipartite state $\rho_{AB}$ is faithful on $A$ if the map
$$
\mathcal E \mapsto \mathcal E_A(\rho_{AB})
$$
is injective on the class of channels under consideration [2206.05899].

The key object is the Jamiołkowski map associated with $\rho_{AB}$,
$$
\rho_{A\to B}(\sigma)=\operatorname{Tr}_A\!\left[(\sigma_A^T\otimes \mathbf 1_B)\rho_{AB}\right],
$$
or equivalently
$$
\rho_{B\to A}(\sigma)=\operatorname{Tr}_B\!\left[(\mathbf 1_A\otimes \sigma_B^T)\rho_{AB}\right].
$$
A rigorous characterization states that a bipartite state is faithful for process tomography on $A$ iff its Jamiołkowski map $\rho_{A\to B}$ is left invertible; equivalently, $\rho_{B\to A}$ is surjective [2206.05899].

An equivalent criterion is expressed through realignment. For
$$
\check{R}(\rho_{AB}) := (\rho_{AB}^{T_B}E)^{T_A},
$$
with swap operator $E=\sum_{ij}|ij\rangle\langle ji|$, usefulness for AAQPT is equivalent to the existence of $\check{R}(\rho_{AB})^{-1}$. The realignment map
$$
\mathcal{R}(\rho_{AB}) = \sum_{i,j,k,l}\varrho_{ij,kl}\,|i\rangle\langle k|\otimes|j\rangle\langle l|
$$
has the same singular values as $\check R(\rho_{AB})$, so the criterion can be stated as full-rank realignment: a bipartite input is useful for AAQPT iff $\mathcal R(\rho_{AB})^{-1}$ exists [2208.05132].

The same work also introduces a weaker notion, sensitivity, defined by
$$
\mathcal E_A(\rho_{AB})=\rho_{AB} \quad\Rightarrow\quad \mathcal E = id_A.
$$
Faithfulness asks whether the channel can be recovered; sensitivity asks only whether any nontrivial channel can be detected. Faithfulness implies sensitivity in general, but the converse need not hold. For classes of maps that form a group, however, the two notions coincide; unitary operations are the canonical example [2206.05899].

For restricted channel classes, several nontrivial equivalences emerge. In particular,
$$
\text{faithful to quantum channels} \;\Longleftrightarrow\; \text{faithful to unital channels} \;\Longleftrightarrow\; \text{faithful to random unitary channels},
$$
whereas faithfulness to unitary operations alone is strictly weaker [2206.05899]. For sensitivity to unital, random-unitary, and unitary channels, the relevant state-side criterion is the absence of any nontrivial local classical observable on the system side. In the paper’s terminology, a PC-Q state satisfies
$$
\rho_{AB}=\sum_i (\Pi_i\otimes \mathbf 1_B)\rho_{AB}(\Pi_i\otimes \mathbf 1_B),
$$
and the characterization is that only bipartite states that has no local classical observable at all can be used to sense the effect of unital channels [2206.05899].

## 3. Correlations, operator Schmidt rank, and the role of entanglement

AAQPT does not require entanglement per se. A decisive quantity is the operator Schmidt rank (OSR), defined from the operator Schmidt decomposition
$$
\rho_{AB} = \sum_{i=1}^{OSR(\rho)} r_i\, A_i \otimes B_i,
$$
where $\{A_i\}$ and $\{B_i\}$ are orthonormal operator bases and $r_i>0$. In the framework of correlation-assisted process tomography, standard AAQPT is possible when
$$
OSR(\rho_{AB}) = d_A^2.
$$
Thus, even a separable but sufficiently correlated state can be enough for AAPT [1808.10835].

This operator-space viewpoint yields a quantitative interpolation between standard QPT and AAQPT. If $OSR(\rho_{AB}) = k$, then one can generate a faithful set using
$$
\left\lceil \frac{d_A^2}{k} \right\rceil
$$
local operations on the probe. The number of required local input preparations therefore scales inversely proportional to the operator Schmidt rank. For pure probe–ancilla states of Schmidt rank $k$, a faithful set can be obtained with
$$
\left\lceil \frac{d}{k} \right\rceil^2
$$
local unitaries, matching the inverse-OSR law because $OSR = k^2$ for pure states [1808.10835].

Several misconceptions are corrected by later work. Entanglement is not necessary for AAQPT usefulness, but it is also not sufficient. There exist entangled states $\varrho_e$ such that $\mathcal R(\varrho_e)$ is singular, hence they are useless for AAQPT. Explicit examples include a two-qutrit entangled state and a two-qutrit bound entangled state whose realigned matrices have zero singular values [2208.05132].

Conversely, certain PPT entangled and bound entangled states can be faithful. A family
$$
\gamma=\mathbb{I}+\mathbb{F}+\varepsilon\,|v\rangle\langle v|,\qquad \varepsilon>0,
$$
with
$$
|v\rangle=\sum_{i=1}^n |a_i\rangle\otimes |b_i\rangle,
$$
can be PPT yet entangled for suitable parameters, while its realignment remains invertible by continuity. The same work gives an explicit $4\otimes 4$ PPT entangled example $\rho_{\mathrm{CCNR}}$ with all realigned singular values strictly positive, so it is faithful for AAQPT [2605.19182].

This literature also separates entanglement detection from tomography usefulness. The CCNR criterion tests separability through
$$
\|\mathcal{R}(\rho_{\mathrm{sep}})\|_1 \le 1,
$$
but AAQPT requires all singular values of $\mathcal R(\rho)$ to be nonzero. Local filtering operations may improve the trace norm of the realignment criterion, yet rank-reducing subspace filters can destroy faithfulness by creating zero singular values. The resulting state may then be unusable for AAQPT despite improved CCNR value [2605.19182].

## 4. Reconstruction schemes and estimator families

Because AAQPT maps process tomography to state tomography on the joint output, it supports several estimator families. In projected least-squares QPT, one first computes the least-squares estimator of the Choi matrix and then projects it onto the convex set of Choi matrices,
$$
CPTP = CP \cap TP,
$$
with
$$
CP=\{\Phi:\Phi\ge 0\},\qquad TP=\{\Phi:\operatorname{Tr}_s(\Phi)=\mathbb 1_d/d\}.
$$
For AAQPT with MUB measurements, the least-squares estimator is
$$
\hat{\Phi}_{LS} = (d^2+1)\sum_{i=1}^m f_i\,|v_i\rangle\langle v_i|-\mathbb 1_{d^2},
$$
and the projection stage can be implemented numerically by the hyperplane intersection projection (HIP) algorithm [2107.01060].

A distinct line of work extends a two-stage solution from standard QPT to AAPT. Starting from an operator-Schmidt decomposition
$$
\sigma^{\text{in}}=\sum_{i=1}^{d_A^2} s_i\, A_i\otimes B_i,
$$
the joint output
$$
\sigma^{\text{out}}=(\mathcal E\otimes I)(\sigma^{\text{in}})
$$
yields
$$
\mathcal E(A_j)=\frac{1}{s_j}\operatorname{Tr}_B\!\left[(I_{d_A}\otimes B_j^\dagger)\sigma^{\text{out}}\right].
$$
The method reconstructs a raw process estimate and then projects it onto the physically admissible set. In this framework, the maximally entangled state is the optimal input state for AAPT [2310.20421].

Classical-shadow methods provide another AAQPT realization. In ancilla-assisted ShadowQPT, one prepares a maximally entangled state, applies the channel to one half, and performs randomized measurements on the resulting Choi state. The snapshot estimator is
$$
\hat{\Lambda}_i = 2^n \mathcal{M}^{-1}\!\left(U_i^\dagger |b_i\rangle\langle b_i| U_i\right),
$$
and averaging reconstructs the Choi matrix. This formulation permits arbitrary a posteriori evaluation of input-output quantities through
$$
\mathrm{tr}\!\left(E(\rho^{in})\sigma\right) = \mathrm{tr}\!\left((\rho^{in\,T}\otimes \sigma)\rho_\Lambda\right),
$$
with favorable scaling for reduced-process tasks [2110.02965].

AAQPT can also be combined with ancilla-assisted quantum state tomography to produce single-shot process tomography (SSPT). In that construction, AAPT encodes the process into one enlarged state, AAQST reconstructs that state in a single collective measurement of commuting observables, and the process matrix $\chi$ is then obtained from the linear system
$$
\beta\chi=\lambda.
$$
The method was demonstrated for several single-qubit processes and a twirling process in a three-qubit NMR register [1404.7830].

## 5. Resource scaling, sample complexity, and optimality

AAQPT does not automatically improve all resource measures. For non-adaptive incoherent measurements, the sample complexity of channel tomography in diamond norm is
$$
\tilde{\Theta}\!\left(\frac{d_{\text{in}}^3 d_{\text{out}}^3}{\varepsilon^2}\right).
$$
A lower bound applies even for ancilla-assisted strategies:
$$
N=\Omega\!\left(\frac{d_{\text{in}}^3 d_{\text{out}}^3}{\varepsilon^2}\right),
$$
while a matching ancilla-free upper bound holds up to logarithmic factors. In this regime, ancillas are allowed, but under non-adaptive incoherent measurements they do not improve the sample complexity [2301.12925].

By contrast, adaptivity changes the asymptotic infidelity behavior. A unified formalism for state, detector, and process tomography introduces the error metric $1-F(\hat S,S)$ and proves that
$$
\mathbb E\bigl(1-F(\hat S,S)\bigr)=O(1/N)
$$
iff both the mean-squared error and the total spurious weight in the estimated zero-eigenspace scale as $O(1/N)$. Guided by this criterion, a three-step adaptive AAQPT algorithm first estimates the joint output state, then measures in the estimated eigenbasis, and finally converts the output-state estimate to a process estimate. For both trace-preserving and non-trace-preserving AAQPT, this yields
$$
\mathbb E\bigl(1-F(\hat X,X)\bigr)=O(1/N),
$$
whereas static methods typically give only $O(1/\sqrt N)$ worst-case scaling [2509.05988].

Finite-sample bounds are also available for projected least-squares estimators. For a $k$-qubit channel with Choi rank $r$, the Frobenius- and trace-norm error bounds depend explicitly on the measurement design and rank, and for low ranks the projection step improves the error rates of the least-squares estimator by a factor $d^2$ [2107.01060]. In the two-stage AAPT framework, the overall computational complexity is
$$
O(Md_A^2d_B^2),
$$
and the error upper bound depends on the measurement design, sample size, ancilla dimension, and the operator-Schmidt coefficients of the input state [2310.20421].

These results together imply a nuanced resource picture. AAQPT can reduce the number of input preparations, can support adaptive $O(1/N)$ infidelity scaling, and can be embedded into fast reconstruction pipelines, yet ancilla assistance alone does not remove the fundamental copy complexity imposed by non-adaptive incoherent measurement models [2301.12925], [2509.05988].

## 6. Experimental realizations and related reinterpretations

AAQPT has been implemented in several experimentally distinct settings. A deployed quantum-network demonstration used one photon of a polarization-entangled pair as the ancilla and the other as the system sent through a 1.6 km deployed fiber-optic link. Using 36 joint polarization projections and Bayesian inference with preconditioned Crank–Nicolson MCMC, the experiment reconstructed the input state, output state, Choi matrix, and Pauli-basis process matrix, reporting a steady process fidelity of 95.1(1)% over a 24 h period [2410.00892].

On near-term quantum hardware, ancilla-assisted ShadowQPT was implemented on the IonQ trapped-ion quantum computer for processes up to $n=4$ qubits, using both Pauli and Clifford measurements. The work emphasizes that once the Choi state has been shadow-reconstructed, many different input-output overlaps can be evaluated classically without rerunning the device [2110.02965].

Realignment-based AAQPT has also been verified experimentally on the IBM platform. In that implementation, the channel
$$
\$(\sigma)=\sum_{n=1}^2 K_n \sigma K_n^\dag,\qquad K_1=\frac{I}{\sqrt 2},\quad K_2=\frac{\sigma_x}{\sqrt 2},
$$
was reconstructed from bipartite input-output state tomography via
$$
M = \mathcal R\big((\$\otimes I)(\varrho_{AB})\big)\, \mathcal R(\varrho_{AB})^{-1},
$$
with reported input and output state fidelities $0.974\pm0.011$ and $0.954\pm0.027$, respectively [2208.05132].

In quantum optics, adaptive ancilla-assisted process tomography has been demonstrated on a two-qubit photonic platform, where the principal qubit is photon polarization and the ancilla qubit is photon path. The reported experiments reached, for the first time, the optimal infidelity scaling in ancilla-assisted process tomography [2509.05988]. Earlier NMR work implemented SSPT in a three-qubit register and characterized several single-qubit gates together with a twirling process through a single collective measurement [1404.7830].

AAQPT has also been used as an interpretive framework outside conventional tomography. In “quantum imaging with undetected photons,” the object can be modeled as an unknown process acting on an idler subsystem, while detected signal photons function as ancilla-like degrees of freedom. In the original measurement configuration, the detector probabilities are
$$
P_{\mathbf{h}/\mathbf{g}}=\frac{1\mp T\cos\gamma}{2},
$$
so only the combination $T\cos\gamma$ is accessible; adding a phase shifter
$$
Z_{\phi}=\begin{pmatrix}1&0\\0&e^{i\phi}\end{pmatrix}
$$
upgrades the scheme from partial AAQPT to full tomography by making $T$ and $\gamma$ separately recoverable [1509.02031].

Across these implementations, a common pattern is that AAQPT is best viewed not as a single protocol but as a family of channel-learning methods organized around one principle: correlations convert a dynamical reconstruction problem into a state-reconstruction problem. The technical content of the modern literature concerns exactly which correlations suffice, how reconstruction should be regularized, and which resource measures are or are not improved by the ancilla.

Source: https://www.emergentmind.com/topics/ancilla-assisted-quantum-process-tomography-aaqpt