Choi–Jamiołkowski Isomorphism

A quantum channel is, in the abstract, a completely positive trace-preserving (CPTP) linear map from operators on one Hilbert space to operators on another — a map between spaces of matrices, unwieldy to reason about or test for physicality. The Choi–Jamiołkowski isomorphism performs a small miracle: it replaces the entire map by a single finite matrix, obtained by feeding one half of a maximally entangled pair through the channel and reading out the state. Every question you might ask of the channel — is it physical, how much information does it carry, does it destroy entanglement, what is its simplest operator-sum form — becomes a question about that one matrix, the Choi matrix. The reason it works is a chain of vector-space identities. Operators on a dd-dimensional space are vectors in a d2d^2-dimensional space; a linear map between operator spaces is therefore a vector in a still larger space; and the maximally entangled state makes this identification physical, probing the channel on a complete basis of inputs simultaneously. The payoff is specific: complete positivity — an infinite family of conditions, one per ancilla — collapses to the single statement J≥0J\ge0; diagonalizing the same matrix hands you a minimal set of Kraus operators; and its normalized version is a density operator whose entropies are the channel's capacities. This is why channel–state duality sits under nearly every proof in quantum Shannon theory. Two roots feed it: Jamiołkowski's 1972 bilinear form and Choi's 1975 characterization of completely positive maps, differing by a partial transpose and a normalization the careless routinely conflate.

Key takeaways
  • The isomorphism is a bijection between linear maps N:L(HA)→L(HB)\mathcal{N}:\mathcal{L}(\mathcal{H}_A)\to\mathcal{L}(\mathcal{H}_B) and Choi matrices J(N)J(\mathcal{N}) on HA⊗HB\mathcal{H}_A\otimes\mathcal{H}_B, built by acting with N\mathcal{N} on half of ∣Ω⟩=∑i∣ii⟩|\Omega\rangle=\sum_i|ii\rangle.
  • Complete positivity of the channel is equivalent to J(N)≥0J(\mathcal{N})\ge 0; trace preservation to the linear condition TrBJ(N)=IA\mathrm{Tr}_B J(\mathcal{N})=I_A. Two hard conditions on a map become two easy ones on a matrix.
  • Diagonalizing JJ yields a Kraus (operator-sum) representation with exactly rank J\mathrm{rank}\,J operators — the minimal number, the Choi rank.
  • The normalized Choi state τ=J/dA\tau=J/d_A is a real density operator whose properties are the channel's: coherent information Ic(N)=−H(A∣B)τI_c(\mathcal{N})=-H(A|B)_\tau, and entanglement-breaking ⇔\Leftrightarrow separable τ\tau.
  • The qubit depolarizing channel has Choi spectrum {2−3p2, p2, p2, p2}\{2-\tfrac{3p}{2},\,\tfrac{p}{2},\,\tfrac{p}{2},\,\tfrac{p}{2}\}, so it is CP exactly for 0≤p≤430\le p\le \tfrac{4}{3}; the transpose map is positive but its Choi matrix is SWAP, spectrum {+1,+1,+1,−1}\{+1,+1,+1,-1\} — the canonical proof that positive ≠\neq CP.
  • Key idea: acting with the channel on half of a maximally entangled state encodes the whole channel into one bipartite operator, N↔J(N)=(I⊗N)(∣Ω⟩⟨Ω∣)\mathcal{N}\leftrightarrow J(\mathcal{N})=(\mathcal{I}\otimes\mathcal{N})(|\Omega\rangle\langle\Omega|); then CP ⇔J≥0\Leftrightarrow J\ge0 and TP ⇔TrBJ=I\Leftrightarrow \mathrm{Tr}_B J=I.
The core insight

Acting with the channel on half of a maximally entangled state encodes the entire channel into one bipartite operator, N↔J(N)=(I⊗N)(∣Ω⟩⟨Ω∣)\mathcal{N}\leftrightarrow J(\mathcal{N})=(\mathcal{I}\otimes\mathcal{N})(|\Omega\rangle\langle\Omega|). Abstract properties of the map become concrete properties of the Choi state: complete positivity is J≥0J\ge0, trace preservation is TrBJ=I\mathrm{Tr}_B J=I, the minimal Kraus number is rank J\mathrm{rank}\,J.

Where this lives today

The isomorphism is the quiet engine under much modern quantum technology. In the laboratory, ancilla-assisted process tomography is the Choi construction made experimental: prepare a maximally entangled pair, send one half through the unknown device, and do full state tomography on the joint output — reconstructing J(N)J(\mathcal{N}), hence the entire channel, with J≥0J\ge0 regularizing the fit. This is how process characterization is benchmarked on superconducting and trapped-ion processors. In certification and design, because CP is J≥0J\ge0 and TP is TrBJ=I\mathrm{Tr}_B J=I, questions like "is this reconstruction a legal channel?" or "what is the diamond-norm distance between two channels?" become semidefinite programs. In quantum Shannon theory, the normalized Choi state converts capacities into state entropies: the LSD quantum-capacity theorem, degradability tests, and the whole apparatus of coherent information rest on manipulating τ(N)\tau(\mathcal{N}). And in fault-tolerant computing, gate teleportation and magic-state injection are channel–state duality in action — a Choi resource state, consumed by Bell measurement and correction, implements the gate on demand, the route to non-transversal gates like the TT gate in surface-code architectures.

The physical picture

Picture the maximally entangled state ∣Ω⟩=∑i∣i⟩A′∣i⟩A|\Omega\rangle=\sum_i|i\rangle_{A'}|i\rangle_A as a perfectly correlated twin pair: whatever basis you choose, measuring A′A' tells you AA. Push only the AA half through N\mathcal{N}, leaving A′A' untouched. Because ∣Ω⟩|\Omega\rangle is the vectorized identity — every input basis element ∣i⟩⟨j∣|i\rangle\langle j| present coherently at once — the joint output learns how N\mathcal{N} acts on the whole operator basis in one shot. That output, as an operator, is the Choi matrix J(N)J(\mathcal{N}): a single bipartite object standing in for the whole map, as the first figure shows.

Figure
The Choi construction sends the A half of a maximally entangled pair Ω through the channel N while the reference A′ is left alone; the untouched A′ and output B jointly carry the Choi operator J, with Tr J = d_A = 2 for a qubit channel, so complete positivity becomes the single check J ≥ 0 and trace preservation becomes Tr over B of J = I
One entangled pair, pushed through N, encodes the whole channel as J(N) prepare |Ω⟩⟨Ω| reference A′ — untouched A N B J(N) on A′ ⊗ B A′ B Tr J = d_A = 2 CP ⟺ J ≥ 0 TP ⟺ Tr_B J = I_A

Why the maximally entangled state, and not some other probe? Because it is the only input whose correlations, after the channel acts, retain full information about N\mathcal{N} on every basis element, coherences included; a less-entangled probe blurs some inputs. This also makes the positivity test powerful: complete positivity demands IR⊗N\mathcal{I}_R\otimes\mathcal{N} preserve positivity for every reference dimension RR, yet the single reference of dimension dAd_A carried by ∣Ω⟩|\Omega\rangle is the worst case — if J≥0J\ge0 here, all larger ancillas pass automatically. The second figure turns this certificate into a quantitative dial, tracking the Choi eigenvalues of a real channel as its noise rises.

Figure
Choi spectrum of the qubit depolarizing channel N(ρ) = (1−p)ρ + p I/2 versus mixing p: the nondegenerate eigenvalue 2 − 3p/2 in blue falls to zero at p = 4/3, the exact complete-positivity boundary, while the triple eigenvalue p/2 in amber rises; the two meet at p = 1 where all four equal 1/2 and J = I/2
Reading CP off the eigenvalues: the depolarizing channel is CP iff 0 ≤ p ≤ 4/3 0 0.5 1 1.5 2 0 0.5 1.0 4/3 1.5 depolarizing parameter p Choi eigenvalue CP fails p > 4/3 2 − 3p/2 → 0 p = 1: all four = 1/2 J = I/2, fully depolarizing λ₁ = 2 − 3p/2 (rank-1 eigenvector) λ₂,₃,₄ = p/2 (triple)

The core equations

Everything begins with a fixed reference state. The unnormalized maximally entangled vector is the vectorized identity — every input basis element with unit weight — and its norm is dAd_A, not 11, the source of the two competing normalization conventions.

∣Ω⟩=∑i=0dA−1∣i⟩A′⊗∣i⟩A(1)|\Omega\rangle = \sum_{i=0}^{d_A-1} |i\rangle_{A'} \otimes |i\rangle_{A} \qquad (1)

where: ∣Ω⟩|\Omega\rangle — unnormalized maximally entangled vector on HA′⊗HA\mathcal{H}_{A'}\otimes\mathcal{H}_A with ⟨Ω∣Ω⟩=dA\langle\Omega|\Omega\rangle=d_A; dAd_A — dimension of the input space; {∣i⟩}\{|i\rangle\} — a fixed orthonormal computational basis. Reading: fixing this basis is what makes the isomorphism basis-dependent, so entanglement statements about JJ must name it. The Choi matrix does all the work: act with N\mathcal{N} on the AA half of ∣Ω⟩⟨Ω∣=∑ij∣i⟩⟨j∣⊗∣i⟩⟨j∣|\Omega\rangle\langle\Omega|=\sum_{ij}|i\rangle\langle j|\otimes|i\rangle\langle j|. The result is a table of the channel's action on every matrix unit.

J(N)=(IA′⊗N)(∣Ω⟩⟨Ω∣)=∑i,j=0dA−1∣i⟩⟨j∣A′⊗N(∣i⟩⟨j∣)B(2)J(\mathcal{N}) = (\mathcal{I}_{A'} \otimes \mathcal{N})\big(|\Omega\rangle\langle\Omega|\big) = \sum_{i,j=0}^{d_A-1} |i\rangle\langle j|_{A'} \otimes \mathcal{N}(|i\rangle\langle j|)_B \qquad (2)

where: J(N)J(\mathcal{N}) — Choi matrix, a dAdB×dAdBd_Ad_B\times d_Ad_B operator on HA′⊗HB\mathcal{H}_{A'}\otimes\mathcal{H}_B with Tr J=dA\mathrm{Tr}\,J=d_A for TP maps; IA′\mathcal{I}_{A'} — identity on the reference; N\mathcal{N} — the map represented; ∣i⟩⟨j∣|i\rangle\langle j| — matrix units, a basis of operators on HA\mathcal{H}_A. Reading: the A′BA'B block indexed by (i,j)(i,j) is exactly N(∣i⟩⟨j∣)\mathcal{N}(|i\rangle\langle j|), so JJ stores the whole map. The construction is invertible: given JJ, the channel is recovered by contracting against the transposed input — the inverse that makes "isomorphism" literal, no information lost.

N(X)=TrA′ ⁣[(XT⊗IB) J(N)](3)\mathcal{N}(X) = \mathrm{Tr}_{A'}\!\big[(X^{T}\otimes I_B)\,J(\mathcal{N})\big] \qquad (3)

where: XX — arbitrary input operator on HA\mathcal{H}_A; XTX^{T} — transpose in the fixed basis; TrA′\mathrm{Tr}_{A'} — partial trace over the reference. Reading: contracting the reference index of JJ against XTX^T superposes the stored blocks with the coefficients of XX. Now the two theorems that make the isomorphism indispensable. Choi's theorem collapses complete-positivity to one semidefinite condition.

N is completely positive  ⟺  J(N)≥0(4)\mathcal{N}\ \text{is completely positive} \iff J(\mathcal{N}) \ge 0 \qquad (4)

where: J(N)≥0J(\mathcal{N})\ge 0 — the Choi matrix is positive semidefinite on HA′⊗HB\mathcal{H}_{A'}\otimes\mathcal{H}_B. Reading: one eigenvalue computation certifies physicality; a negative eigenvalue disqualifies the map. Trace preservation — probabilities summing to one — becomes an equally clean linear condition.

N is trace-preserving  ⟺  TrB J(N)=IA(5)\mathcal{N}\ \text{is trace-preserving} \iff \mathrm{Tr}_B\,J(\mathcal{N}) = I_A \qquad (5)

where: TrB\mathrm{Tr}_B — partial trace over the output; IAI_A — identity on the input space. Reading: (4) and (5) together say a legal CPTP channel is exactly a positive JJ whose BB-marginal is the identity — the constraints of a semidefinite program. Diagonalizing JJ produces Kraus operators directly: each de-vectorized eigenvector is a Kraus operator, and their number is rank J\mathrm{rank}\,J.

J(N)=∑k∣Kk⟩ ⁣⟩⟨ ⁣⟨Kk∣,∣Kk⟩ ⁣⟩=∑ij(Kk)ji ∣i⟩∣j⟩(6)J(\mathcal{N}) = \sum_k |K_k\rangle\!\rangle\langle\!\langle K_k|, \qquad |K_k\rangle\!\rangle = \sum_{ij}(K_k)_{ji}\,|i\rangle|j\rangle \qquad (6)

where: KkK_k — Kraus operators with N(ρ)=∑kKkρKk†\mathcal{N}(\rho)=\sum_k K_k\rho K_k^\dagger; ∣Kk⟩ ⁣⟩|K_k\rangle\!\rangle — vectorization of KkK_k; scaled eigenvectors of JJ furnish the set. Reading: the Choi picture is the most direct route from a channel to its operator-sum form, delivering the minimal set automatically. Finally, normalizing JJ makes a density operator whose entropies are channel information quantities — the bridge into quantum Shannon theory.

Ic(N)=H(B)τ−H(AB)τ=−H(A∣B)τ,τ=1dA J(N)(7)I_c(\mathcal{N}) = H(B)_{\tau} - H(AB)_{\tau} = -H(A|B)_{\tau}, \qquad \tau = \tfrac{1}{d_A}\,J(\mathcal{N}) \qquad (7)

where: τ\tau — normalized Choi state on A′BA'B; H(B)τH(B)_\tau — entropy of the output marginal; H(A∣B)τH(A|B)_\tau — conditional entropy of the Choi state; Ic(N)I_c(\mathcal{N}) — coherent information for maximally mixed input. Reading: "how much quantum information the channel preserves" becomes "how negative the conditional entropy of τ\tau is," where strong subadditivity and data processing apply.

Derivation AChoi's theorem — complete positivity equals one positivity check

The goal is to prove, in both directions, that a map N\mathcal{N} is completely positive if and only if its Choi matrix J(N)=(I⊗N)(∣Ω⟩⟨Ω∣)J(\mathcal{N})=(\mathcal{I}\otimes\mathcal{N})(|\Omega\rangle\langle\Omega|) is positive semidefinite, and to see that the hard direction is constructive: positivity of JJ literally manufactures Kraus operators, from which the operator-sum form and hence complete positivity follow. Recall N\mathcal{N} is completely positive precisely when In⊗N\mathcal{I}_n\otimes\mathcal{N} preserves positivity for every ancilla dimension nn — an infinite family of requirements this theorem reduces to one.

Step 1. (⇒\Rightarrow) Suppose N\mathcal{N} is completely positive. The reference ∣Ω⟩⟨Ω∣|\Omega\rangle\langle\Omega| is a rank-one projector, hence positive; CP with n=dAn=d_A says IdA⊗N\mathcal{I}_{d_A}\otimes\mathcal{N} preserves positivity, so J≥0J\ge 0. Why: the definition of CP applied to the single entangled input forces J≥0J\ge0 with no further work.

Step 2. (⇐\Leftarrow) Suppose J≥0J\ge0. Diagonalize J=∑kλk∣vk⟩⟨vk∣J=\sum_k\lambda_k|v_k\rangle\langle v_k| with λk≥0\lambda_k\ge0 and absorb eigenvalues into ∣Kk⟩ ⁣⟩=λk ∣vk⟩|K_k\rangle\!\rangle=\sqrt{\lambda_k}\,|v_k\rangle, so J=∑k∣Kk⟩ ⁣⟩⟨ ⁣⟨Kk∣J=\sum_k|K_k\rangle\!\rangle\langle\!\langle K_k|. Why: a positive matrix always factorizes as a sum of ∣w⟩⟨w∣|w\rangle\langle w| — the only place positivity is used; the rest is bookkeeping.

Step 3. De-vectorize. Every ∣Kk⟩ ⁣⟩∈HA′⊗HB|K_k\rangle\!\rangle\in\mathcal{H}_{A'}\otimes\mathcal{H}_B corresponds to a unique operator Kk:HA→HBK_k:\mathcal{H}_A\to\mathcal{H}_B via ∣Kk⟩ ⁣⟩=(I⊗Kk)∣Ω⟩|K_k\rangle\!\rangle=(I\otimes K_k)|\Omega\rangle, i.e. (Kk)ji=⟨j∣Kk∣i⟩(K_k)_{ji}=\langle j|K_k|i\rangle. Why: this vector–operator dictionary is a bijection, so the eigenvectors name a genuine operator set.

Step 4. Reconstruct. Feeding this JJ into the inverse formula (3), the algebra collapses to N(X)=∑kKkXKk†\mathcal{N}(X)=\sum_k K_k X K_k^{\dagger}. Why: contracting the factorized JJ reproduces the sandwich form, tying abstract JJ to concrete operators.

Step 5. Conclude. Any operator-sum map is manifestly CP, since In⊗N=∑k(In⊗Kk)(⋅)(In⊗Kk)†\mathcal{I}_n\otimes\mathcal{N}=\sum_k(I_n\otimes K_k)(\cdot)(I_n\otimes K_k)^\dagger preserves positivity for every nn. Hence J≥0⇒NJ\ge0\Rightarrow\mathcal{N} CP, closing the equivalence. Why: the single positivity check is sufficient as well as necessary, exhibiting the Kraus operators explicitly. *Verified: cross-checked against Watrous Ch. 2 (Thm 2.22) and Nielsen–Chuang §8.2.3; the de-vectorization (I⊗K)∣Ω⟩(I\otimes K)|\Omega\rangle matches (K)ji=⟨j∣K∣i⟩(K)_{ji}=\langle j|K|i\rangle.*

Go deeper · Derivation B · 1 minThe inverse map and the trace-preservation condition
Derivation BThe inverse map and the trace-preservation condition

The goal is to establish two facts that together promote the Choi map to a full isomorphism with a clean physicality criterion: first, that equation (3) genuinely inverts equation (2), recovering N(X)\mathcal{N}(X) from JJ for any input XX; second, that conservation of probability Tr N(X)=Tr X\mathrm{Tr}\,\mathcal{N}(X)=\mathrm{Tr}\,X is equivalent to the marginal condition TrBJ=IA\mathrm{Tr}_B J=I_A. Both follow from expanding in matrix units and pushing indices, using nothing beyond linearity and the partial trace.

Step 1. Expand X=∑ijXij∣i⟩⟨j∣X=\sum_{ij}X_{ij}|i\rangle\langle j|; by linearity it suffices to verify the inverse on each matrix unit and reassemble. Why: channel and trace are both linear, so tracking a basis of operators controls the general case.

Step 2. Insert the definition of JJ. Since (XT)ij=Xji(X^T)_{ij}=X_{ji}, the reference-block product (XT⊗IB) J=∑ijklXji ∣i⟩⟨j∣k⟩⟨l∣⊗N(∣k⟩⟨l∣)(X^{T}\otimes I_B)\,J=\sum_{ijkl}X_{ji}\,|i\rangle\langle j|k\rangle\langle l|\otimes\mathcal{N}(|k\rangle\langle l|) forces j=kj=k via ⟨j∣k⟩=δjk\langle j|k\rangle=\delta_{jk}. Why: the transpose is the twist that pairs the coefficients of XX with the correct stored blocks of JJ.

Step 3. Partial-trace over A′A'. The surviving deltas collapse the sums to TrA′[⋯ ]=∑jlXjl N(∣j⟩⟨l∣)=N(X)\mathrm{Tr}_{A'}[\cdots]=\sum_{jl}X_{jl}\,\mathcal{N}(|j\rangle\langle l|)=\mathcal{N}(X). Why: the contraction reassembles precisely N\mathcal{N} acting on the original XX, so (3) inverts (2).

Step 4. Now impose trace preservation. On a matrix unit, Tr N(∣i⟩⟨j∣)=δij\mathrm{Tr}\,\mathcal{N}(|i\rangle\langle j|)=\delta_{ij}; the reference marginal is ⟨i∣(TrBJ)∣j⟩=TrB N(∣i⟩⟨j∣)=δij\langle i|(\mathrm{Tr}_B J)|j\rangle=\mathrm{Tr}_B\,\mathcal{N}(|i\rangle\langle j|)=\delta_{ij}, exactly IAI_A. Why: probability conservation on every input is captured by the single equation TrBJ=IA\mathrm{Tr}_B J=I_A. *Verified: direct index manipulation reproduces N(X)\mathcal{N}(X) and TrBJ=IA\mathrm{Tr}_B J=I_A; consistency Tr J=Tr IA=dA\mathrm{Tr}\,J=\mathrm{Tr}\,I_A=d_A matches equation (2).*

Go deeper · Derivation C · 1 minFrom the Choi matrix to a minimal Kraus representation
Derivation CFrom the Choi matrix to a minimal Kraus representation

The goal is to show that the eigendecomposition of J(N)J(\mathcal{N}) yields a Kraus set of size rank J\mathrm{rank}\,J, that this is the smallest possible number (the Choi rank), and that the operators form a canonical orthogonal set with trace preservation doubling as the Kraus completeness relation. This sharpens the previous derivation: not only does J≥0J\ge0 give some Kraus operators, its spectrum gives the minimal and canonical ones.

Step 1. Spectrally decompose J=∑k=1rλk∣vk⟩⟨vk∣J=\sum_{k=1}^{r}\lambda_k|v_k\rangle\langle v_k|, keeping only the r=rank Jr=\mathrm{rank}\,J strictly positive eigenvalues. Why: zero eigenvalues contribute nothing, so the sum runs over exactly rr terms and fixes the operator count.

Step 2. Define Kk=λk vec−1(∣vk⟩)K_k=\sqrt{\lambda_k}\,\mathrm{vec}^{-1}(|v_k\rangle), i.e. (Kk)ji=λk ⟨ji∣vk⟩(K_k)_{ji}=\sqrt{\lambda_k}\,\langle ji|v_k\rangle. Why: this converts the eigenbasis of JJ into Kraus operators, one per positive eigenvalue.

Step 3. Read off completeness. Term by term, TrBJ=IA\mathrm{Tr}_B J=I_A is exactly ∑kKk†Kk=IA\sum_k K_k^{\dagger}K_k=I_A, the Kraus completeness relation; and eigenvector orthonormality gives Tr(Kk†Kl)=λk δkl\mathrm{Tr}(K_k^{\dagger}K_l)=\lambda_k\,\delta_{kl}, so the eigen-Kraus operators are Hilbert–Schmidt orthogonal — a canonical set. Why: the JJ-conditions that certify physicality reappear as the standard Kraus constraints.

Step 4. Minimality. Any representation {Lm}\{L_m\} satisfies J=∑m∣Lm⟩ ⁣⟩⟨ ⁣⟨Lm∣J=\sum_m|L_m\rangle\!\rangle\langle\!\langle L_m|, so the ∣Lm⟩ ⁣⟩|L_m\rangle\!\rangle span the range of JJ, of dimension rr; no representation uses fewer than rr operators, and the eigen-construction attains the bound. Why: the range of JJ is a representation-independent invariant pinning the minimal Kraus number at rank J\mathrm{rank}\,J. *Verified: Choi rank == minimal Kraus number is standard (Wilde Ch. 4); orthogonality via ⟨ ⁣⟨Kk∣Kl⟩ ⁣⟩=Tr(Kk†Kl)=λkδkl\langle\!\langle K_k|K_l\rangle\!\rangle=\mathrm{Tr}(K_k^\dagger K_l)=\lambda_k\delta_{kl}; for the depolarizing channel below, r=4r=4 recovers the four Pauli Kraus operators.*

Limiting cases

Special channels wear especially transparent Choi matrices — each a check the formalism must reproduce.

CaseResultMeaning
Identity channelJ(I)=∣Ω⟩⟨Ω∣J(\mathcal{I}) = \vert\Omega\rangle\langle\Omega\vertThe Choi matrix is the unnormalized maximally entangled state — rank one; 12J\tfrac12 J is the Bell state Φ+\Phi^+, entanglement entropy 11 bit.
Completely depolarizingJ(Ndep)=1d IA⊗IBJ(\mathcal{N}_{\mathrm{dep}}) = \tfrac{1}{d}\,I_A\otimes I_BOutputting I/dI/d gives a product (uncorrelated) Choi operator; for a qubit 12I4\tfrac12 I_4, the p=1p=1 point of the figure.
Unitary channelJ(U)=(I⊗U)∣Ω⟩⟨Ω∣(I⊗U†)J(\mathcal{U}) = (I\otimes U)\vert\Omega\rangle\langle\Omega\vert(I\otimes U^{\dagger})Rank one and maximally entangled for any UU; Choi entanglement entropy log⁡d\log d, a unitary preserving all information.
Entanglement-breakingJ(N)=∑kpk σk⊗ρkJ(\mathcal{N}) = \sum_k p_k\,\sigma_k\otimes\rho_kSeparable Choi state — the defining characterization: a channel breaks entanglement iff its Choi state is separable.

Worked examples

1 · The qubit depolarizing channel — reading CP and rank off the spectrum (fully worked)

The qubit depolarizing channel is N(ρ)=(1−p)ρ+p Tr(ρ) I2\mathcal{N}(\rho)=(1-p)\rho+p\,\mathrm{Tr}(\rho)\,\tfrac{I}{2}. Build its Choi matrix, find the spectrum, and determine for which pp it is a legal channel and how many Kraus operators it needs.

Step 1. Compute N\mathcal{N} on the operator basis (using Tr(∣i⟩⟨j∣)=δij\mathrm{Tr}(|i\rangle\langle j|)=\delta_{ij}): N(∣0⟩⟨0∣)=diag(1−p2, p2)\mathcal{N}(|0\rangle\langle0|)=\mathrm{diag}(1-\tfrac{p}{2},\,\tfrac{p}{2}), N(∣1⟩⟨1∣)=diag(p2, 1−p2)\mathcal{N}(|1\rangle\langle1|)=\mathrm{diag}(\tfrac{p}{2},\,1-\tfrac{p}{2}), off-diagonals shrinking: N(∣0⟩⟨1∣)=(1−p)∣0⟩⟨1∣\mathcal{N}(|0\rangle\langle1|)=(1-p)|0\rangle\langle1|, N(∣1⟩⟨0∣)=(1−p)∣1⟩⟨0∣\mathcal{N}(|1\rangle\langle0|)=(1-p)|1\rangle\langle0|. These four blocks are the A′BA'B blocks of JJ.

Step 2. Assemble J=∑ij∣i⟩⟨j∣⊗N(∣i⟩⟨j∣)J=\sum_{ij}|i\rangle\langle j|\otimes\mathcal{N}(|i\rangle\langle j|) in the basis {00,01,10,11}\{00,01,10,11\}:

J=(1−p2001−p0p20000p201−p001−p2)J = \begin{pmatrix} 1-\tfrac{p}{2} & 0 & 0 & 1-p \\ 0 & \tfrac{p}{2} & 0 & 0 \\ 0 & 0 & \tfrac{p}{2} & 0 \\ 1-p & 0 & 0 & 1-\tfrac{p}{2} \end{pmatrix}

The {01,10}\{01,10\} sector is diagonal with entries p2\tfrac{p}{2}; the {00,11}\{00,11\} sector is (1−p/21−p1−p1−p/2)\begin{pmatrix}1-p/2 & 1-p\\ 1-p & 1-p/2\end{pmatrix}.

Step 3. Diagonalize. The 2×22\times2 block has eigenvalues (1−p2)±(1−p)={ 2−3p2, p2 }(1-\tfrac{p}{2})\pm(1-p)=\{\,2-\tfrac{3p}{2},\ \tfrac{p}{2}\,\}, so the full spectrum is { 2−3p2, p2, p2, p2 }\boxed{\{\,2-\tfrac{3p}{2},\ \tfrac{p}{2},\ \tfrac{p}{2},\ \tfrac{p}{2}\,\}} with Tr J=2=dA\mathrm{Tr}\,J=2=d_A. This is the spectrum plotted in the second figure.

Step 4. Read off complete positivity. All eigenvalues ≥0\ge0: the triple needs p≥0p\ge0, the singleton 2−3p2≥02-\tfrac{3p}{2}\ge0, i.e. p≤43p\le\tfrac{4}{3}. Hence N\mathcal{N} is CP exactly for 0≤p≤430\le p\le\tfrac{4}{3} — a nontrivial physical bound from Choi eigenvalues alone. At p=1.4>43p=1.4>\tfrac43 the singleton is −0.1<0-0.1<0 and the "channel" is unphysical.

Step 5. Count Kraus operators. For generic p∈(0,43)p\in(0,\tfrac43) the rank is 44, so the minimal set is four operators — the familiar {1−3p4 I, p4 X, p4 Y, p4 Z}\{\sqrt{1-\tfrac{3p}{4}}\,I,\ \sqrt{\tfrac{p}{4}}\,X,\ \sqrt{\tfrac{p}{4}}\,Y,\ \sqrt{\tfrac{p}{4}}\,Z\}. At the endpoints the rank drops: p=0p=0 gives rank 11, p=43p=\tfrac43 gives rank 33. CP range p∈[0,43]p\in[0,\tfrac43], Choi rank 44 generically.

2 · The transpose map — positive but not completely positive

The transpose T(ρ)=ρT\mathcal{T}(\rho)=\rho^{T} maps density operators to density operators, so it is positive. Is it a channel?

Step 1. Apply the definition: J(T)=∑ij∣i⟩⟨j∣⊗(∣i⟩⟨j∣)T=∑ij∣i⟩⟨j∣⊗∣j⟩⟨i∣J(\mathcal{T})=\sum_{ij}|i\rangle\langle j|\otimes(|i\rangle\langle j|)^{T}=\sum_{ij}|i\rangle\langle j|\otimes|j\rangle\langle i| — the SWAP operator on two qubits.

Step 2. SWAP acts as +1+1 on the symmetric subspace (dimension 3) and −1-1 on the antisymmetric singlet, so its spectrum is {+1,+1,+1,−1}\{+1,+1,+1,-1\}.

Step 3. The eigenvalue −1-1 means J(T)≱0J(\mathcal{T})\not\ge0, so by Choi's theorem the transpose is not completely positive, though it is positive. Concretely, (I⊗T)(\mathcal{I}\otimes\mathcal{T}) on a Bell state gives an operator with a negative eigenvalue — the partial transpose behind the PPT entanglement criterion. Positive ≠\neq CP; the witness is the eigenvalue −1-1 of the SWAP.

3 · Coherent information from the Choi state

Use the normalized Choi state to compute the coherent information of the qubit depolarizing channel and find where it vanishes.

Step 1. Normalize: τ=J/2\tau=J/2 has spectrum {1−3p4, p4, p4, p4}\{1-\tfrac{3p}{4},\,\tfrac{p}{4},\,\tfrac{p}{4},\,\tfrac{p}{4}\} and BB-marginal I/2I/2, so H(B)τ=1H(B)_\tau=1 bit for all pp. Then Ic=H(B)τ−H(AB)τ=1−H(τ)I_c=H(B)_\tau-H(AB)_\tau=1-H(\tau).

Step 2. At p=0p=0, H(τ)=0H(\tau)=0 and Ic=1I_c=1 bit. Increasing pp raises H(τ)H(\tau): at p=0.25p=0.25, H(τ)≈0.993H(\tau)\approx0.993 so Ic≈0.007I_c\approx0.007 bit, barely positive; by p=13p=\tfrac13, Ic≈−0.21I_c\approx-0.21 bit. The single-letter coherent information crosses zero near p≈0.25p\approx0.25 — a capacity statement from a state entropy alone.

Physics of the portrait

The isomorphism is where operator-algebra theorems from the 1950s–1970s became the daily tools of quantum information; two names for two conventions differing by a partial transpose.

Origins · 1 minChoi–Jamiołkowski isomorphism and how this was discovered
CI
Choi–Jamiołkowski isomorphism
channel–state duality J(N) = (I⊗N)( · Ω⟩⟨Ω

The result grew from operator-algebra theory, decades before "quantum channel" was a phrase anyone used. In 1955 W. Forrest Stinespring proved his dilation theorem, characterizing completely positive maps as compressions of *-homomorphisms — the ancestor of every channel representation. In 1972 the Polish mathematical physicist Andrzej Jamiołkowski associated to each linear map a bipartite operator (a bilinear form), the first version of the correspondence now bearing his name; his convention is basis-independent and sits one partial transpose from the modern matrix. In 1975 Man-Duen Choi, a Hong Kong-born operator theorist then at Toronto, proved the sharp statement: a linear map on complex matrices is completely positive iff the single Choi matrix ∑ijEij⊗N(Eij)\sum_{ij}E_{ij}\otimes\mathcal{N}(E_{ij}) is positive semidefinite, its eigen-decomposition yielding the operator-sum form. For two decades this lived in mathematics; then quantum information found its central objects — channels, capacities, degradability — most naturally studied as properties of the Choi state. The 2003 characterization of entanglement-breaking channels by Horodecki, Shor and Ruskai (separable Choi state) is emblematic, and reading τ=1dJ\tau=\tfrac1d J as a state whose entropies are capacities completes the arc from abstract theorem to workhorse.

History
1955
— Stinespring proves the dilation theorem for completely positive maps
1972
— Jamiołkowski associates a bipartite operator to each linear map
1975
— Choi characterizes complete positivity by positivity of the Choi matrix
2003
— Horodecki, Shor and Ruskai tie entanglement-breaking to separable Choi states

Going further

References. Nielsen & Chuang, Quantum Computation and Quantum Information, §8.2 · Watrous, The Theory of Quantum Information, Ch. 2 and 4 · Wilde, Quantum Information Theory (2nd ed.), Ch. 4–5, §13 · Preskill, Lecture Notes for Physics 219, Ch. 3 · Kitaev, Shen, Vyalyi, Classical and Quantum Computation, Ch. 11 · M.-D. Choi, "Completely positive linear maps on complex matrices," Lin. Alg. Appl. 10 (1975) 285 · A. Jamiołkowski, Rep. Math. Phys. 3 (1972) 275 · Horodecki, Shor, Ruskai, "Entanglement breaking channels," Rev. Math. Phys. 15 (2003) 629.

Provenance. Definitions and theorems follow Watrous (Ch. 2), Wilde (Ch. 4–5), and Nielsen & Chuang (§8.2). A verifier should: (1) confirm the inverse-map and TP derivations by direct index manipulation; (2) check the transpose map's Choi matrix is SWAP with spectrum {+1,+1,+1,−1}\{+1,+1,+1,-1\}; (3) verify the depolarizing Choi spectrum {2−3p2,p2,p2,p2}\{2-\tfrac{3p}{2},\tfrac{p}{2},\tfrac{p}{2},\tfrac{p}{2}\} and CP window p∈[0,43]p\in[0,\tfrac43] (numerically confirmed); (4) verify Ic=−H(A∣B)τI_c=-H(A|B)_\tau. Anchor numbers (Bell-state entropy 11 bit, Tr J=dA\mathrm{Tr}\,J=d_A, Choi rank == minimal Kraus number, depolarizing bound 43\tfrac43) are standard. Attributions: Stinespring (1955), Jamiołkowski (1972), Choi (1975), Horodecki–Shor–Ruskai (2003).

Prerequisites. Quantum Channels and CPTP Maps · Kraus Representation of Quantum Channels · Von Neumann Entropy. Leads to. Stinespring Dilation Theorem · Complementary Channels and Degradability · Data-Processing Inequality · Holevo Bound.

Ask the tutor: Want to see the Choi matrix, its eigenvalues, and a minimal Kraus set worked out for any channel you name — amplitude damping, dephasing, or your own noisy gate?

Why does one entangled input certify complete positivity, when the definition quantifies over all ancilla dimensions?

Because the maximally entangled reference of dimension dAd_A is the worst case. CP asks that In⊗N\mathcal{I}_n\otimes\mathcal{N} preserve positivity for every nn, but if it passes on the maximally entangled state of the input dimension, it passes for every larger ancilla automatically: ∣Ω⟩|\Omega\rangle already carries maximal correlations across the full input space, so no bigger ancilla builds a harder input. This is the whole miracle — an infinite family of checks reduces to one eigenvalue computation on JJ.

What is the difference between the Choi matrix and the Choi state?

Only normalization, but the distinction matters. The Choi matrix J=(I⊗N)(∣Ω⟩⟨Ω∣)J=(\mathcal{I}\otimes\mathcal{N})(|\Omega\rangle\langle\Omega|) uses the unnormalized ∣Ω⟩|\Omega\rangle, has Tr J=dA\mathrm{Tr}\,J=d_A, and obeys the clean marginal law TrBJ=IA\mathrm{Tr}_B J=I_A. The Choi state τ=J/dA\tau=J/d_A is a genuine density operator with Tr τ=1\mathrm{Tr}\,\tau=1. Use JJ for CP/TP certification and Kraus extraction; use τ\tau for information theory, where you need a normalized state for von Neumann entropies. Mixing the two normalizations is a frequent slip.

How is the Choi–Jamiołkowski isomorphism related to Kraus and Stinespring?

They are three faces of one channel. Diagonalizing the Choi matrix and de-vectorizing its eigenvectors gives Kraus operators directly, the count equalling the minimal (Choi-rank) number. Purifying the Choi state and tracing out the environment recovers the Stinespring dilation, in which the channel is an isometry followed by a partial trace. So Choi →\to Kraus is an eigen-decomposition and Choi →\to Stinespring is a purification; the isomorphism lets you move freely between representations, choosing whichever suits the proof at hand.

Why does the transpose map show that positive is weaker than completely positive?

The transpose ρ↦ρT\rho\mapsto\rho^T sends every density operator to a valid density operator, so it is a positive map. But its Choi matrix is the SWAP operator, whose spectrum {+1,+1,+1,−1}\{+1,+1,+1,-1\} has a negative eigenvalue, so it fails J≥0J\ge0 and is not completely positive. Applied to half of an entangled pair, the transpose produces a joint operator with a negative eigenvalue — impossible for a physical process, and exactly the signal exploited by the PPT entanglement test. No accident that the textbook "positive but not CP" map also underlies a famous entanglement criterion.

Does the Choi matrix depend on the choice of basis, and does that ever cause trouble?

Yes — it is defined relative to the fixed basis named by ∣Ω⟩=∑i∣ii⟩|\Omega\rangle=\sum_i|ii\rangle. For the CP and TP criteria this is harmless: whether J≥0J\ge0 or TrBJ=I\mathrm{Tr}_B J=I is unchanged under a coordinated basis rotation. Where it bites is entanglement and separability statements about JJ, and the relation to Jamiołkowski's convention, which differs by a partial transpose. Since a partial transpose swaps entanglement structure with PPT structure, pin down both the basis and the transpose convention before claiming anything about JJ's entanglement — the practical content of the "two Choi matrices" warning.

Can the isomorphism be used experimentally, or is it only a proof technique?

Very much experimentally. Ancilla-assisted process tomography is the Choi construction run in a lab: prepare a maximally entangled pair, send one half through the device, and do full state tomography on the joint output. The reconstructed bipartite state is a normalization of the Choi matrix, from which you obtain the complete channel — every Kraus operator, the process fidelity, the diamond-norm error — with J≥0J\ge0 imposed during the fit. This is often more robust and sample-efficient than preparing a tomographically complete set of separate inputs, and it is standard in benchmarking real processors.

Before this

Where this leads