Trace Distance and State Distinguishability

Every quantum protocol eventually reduces to one sharp question: handed a system in one of two known states ρ\rho or σ\sigma, how well can any measurement tell which you were given? Classically the answer is the total variation distance between two probability distributions. The trace distance T(ρ,σ)=12∥ρ−σ∥1T(\rho,\sigma)=\tfrac12\|\rho-\sigma\|_1 is the exact quantum generalization, and its power is that it is not a convenience bolted on afterward — it is the answer. Carl Helstrom's 1967 solution of the two-state detection problem shows the best guess succeeds with probability 12(1+T)\tfrac12(1+T), so TT is the bias of the optimal measurement. Three structural facts organize everything that follows. First, because ρ−σ\rho-\sigma is Hermitian and traceless, its eigenvalues split into a positive bundle and an equal negative bundle; TT is the weight of one bundle, and the optimal measurement projects onto the positive one. Second, TT is a genuine metric — symmetric, nonnegative, zero only when ρ=σ\rho=\sigma, obeying the triangle inequality — which is why error budgets in fault tolerance and cryptography are written in trace distance. Third, TT is contractive under every quantum channel: no processing can make two states easier to tell apart. This article develops all three, ties TT to fidelity through Fuchs–van de Graaf and to relative entropy through Pinsker, and shows where it does real work — from composable key distribution to error-correction thresholds.

Key takeaways
  • The trace distance T(ρ,σ)=12∥ρ−σ∥1T(\rho,\sigma)=\tfrac12\|\rho-\sigma\|_1 is a metric on density operators, equal to the sum of the positive eigenvalues of the traceless operator ρ−σ\rho-\sigma, with range 0≤T≤10\le T\le 1.
  • It is operational: given one of two equiprobable states, the best probability of guessing correctly is psucc=12(1+T)p_{\mathrm{succ}}=\tfrac12(1+T), achieved by the Helstrom measurement projecting onto the positive part of ρ−σ\rho-\sigma.
  • Equivalently TT is the largest classical statistical distance any POVM can produce between the two outcome distributions, and for a single qubit it is half the Euclidean distance between the Bloch vectors — orthogonal states a full diameter apart at T=1T=1.
  • It is contractive under every CPTP map, sandwiched to fidelity by 1−F≤T≤1−F21-F\le T\le\sqrt{1-F^2}, and lower-bounds relative entropy through quantum Pinsker.
  • Key idea: half the trace norm of ρ−σ\rho-\sigma is exactly the bias of the best possible measurement distinguishing the states, and the optimal POVM projects onto the positive eigenspace of ρ−σ\rho-\sigma.
The core insight

Half the trace norm of ρ−σ\rho-\sigma is exactly the bias of the best possible two-outcome measurement distinguishing the states; the optimal POVM projects onto the positive eigenspace of ρ−σ\rho-\sigma, and every other property of TT — metric, contractive, fidelity-bounded — follows from this one operational fact.

Where this lives today

Trace distance is the error bar of quantum information science, doing its most consequential work in security proofs. A composably secure quantum key distribution protocol is defined by a trace-distance statement: the real joint state of key and eavesdropper must lie within ϵ\epsilon of an ideal state whose key is uniformly random and decoupled from the adversary. The Helstrom interpretation makes this the right definition — an adversary distinguishing real from ideal succeeds with advantage at most ϵ\epsilon, and because TT is a metric the guarantee survives composition with any protocol that consumes the key. The same quantity underwrites fault tolerance, where the diamond norm certifying a noisy gate reduces to output trace distances and the triangle inequality gives the linear-in-depth error accumulation threshold theorems need. In the lab, randomized benchmarking and gate-set tomography report infidelities that convert into trace-distance bounds on how far a real processor sits from its target — the number a superconducting or trapped-ion team watches fall as it calibrates. At the frontier, photonic experiments have realized Helstrom-optimal receivers for coherent states, beating the homodyne standard quantum limit and approaching the very bound Helstrom wrote down for radar half a century before the apparatus existed.

The physical picture

Imagine an experimenter who receives a qubit prepared by a fair coin flip: heads it is ρ\rho, tails σ\sigma. She may perform any measurement quantum mechanics allows and must then guess which. Her problem is not noise — she gets one perfect copy — but the irreducible ambiguity of non-orthogonal states, and the trace distance is the exact currency of that ambiguity. For a single qubit the whole situation lives inside the Bloch ball: a state is a point r⃗\vec r, and the trace distance between two states is half the straight-line distance between their points. Pure states sit on the surface, orthogonal pure states at antipodes a full diameter apart, the maximally mixed state dead center.

Figure
Two qubit states on the Bloch sphere, one at the north pole and one on the equator, are unit Bloch vectors 90° apart, so their Euclidean separation is √2 and the trace distance is exactly half that distance, T = ½·√2 ≈ 0.707, giving Helstrom success ½(1 + 0.707) ≈ 0.854
Trace distance = half the Bloch-vector separation: T = ½ |r_ρ − r_σ| |0⟩ ρ |+⟩ σ √2 Read off the number r_ρ = (0,0,1), r_σ = (1,0,0) |r_ρ − r_σ| = √2 ≈ 1.414 T = ½·√2 = 1/√2 ≈ 0.707 overlap |⟨0|+⟩|² = ½ √(1 − ½) = 1/√2 agrees p_succ = ½(1+T) ≈ 0.854

The Helstrom strategy behind that success probability is almost disappointingly simple. Form Δ=ρ−σ\Delta=\rho-\sigma. Because both states have unit trace, Δ\Delta is traceless, so every positive eigenvalue (where ρ\rho is "heavier") is matched by a negative one (where σ\sigma is heavier). The optimal measurement asks one question — is the state more likely ρ\rho or σ\sigma? — by projecting onto the positive eigenspace of Δ\Delta, guessing ρ\rho there and σ\sigma elsewhere; the advantage over a coin flip is the total positive weight, TT. The second figure makes the split concrete for the states above, where Δ\Delta has eigenvalues ±0.707\pm0.707 and the positive bundle sums to T=0.707T=0.707.

Figure
The difference operator Δ = ρ − σ is traceless, so its eigenvalues split into a positive bundle and an equal negative bundle; for the two states shown the eigenvalues are exactly ±0.707, the trace distance is the positive sum T = 0.707, and the optimal Helstrom measurement projects onto the positive bundle Π₊
Eigenvalues of Δ = ρ − σ split evenly about zero; T = Σ (positive eigenvalues) eigenvalue axis +0.71 −0.71 λ₊ = +0.707 λ₋ = −0.707 Π₊ Π₋ guess ρ here Bookkeeping Tr Δ = 0 T = Σλ₊ = 0.707 ‖Δ‖₁ = 2T = 1.414 p = ½(1+T) ≈ 0.854

Contractivity closes the physical picture. A quantum channel E\mathcal{E} is a noisy, possibly irreversible process — a decohering wire, a lossy fiber, an eavesdropper. Whatever it does, it can only blur the difference between two inputs, never sharpen it, so T(E(ρ),E(σ))≤T(ρ,σ)T(\mathcal{E}(\rho),\mathcal{E}(\sigma))\le T(\rho,\sigma): distinguishability is a resource processing spends but cannot create. Because TT is at once a metric and monotone, "the real device is within ϵ\epsilon of the ideal" composes cleanly across a whole circuit — errors add rather than compound — exactly what makes threshold theorems and composable security provable.

The core equations

The definition is built from the trace norm, the Schatten 1-norm. For the Hermitian ρ−σ\rho-\sigma this is the sum of absolute eigenvalues; the one half pins the range to [0,1][0,1] and matches the classical total variation distance when the states commute.

T(ρ,σ)=12 ∥ρ−σ∥1=12 Tr ∣ρ−σ∣(1)T(\rho,\sigma)=\tfrac12\,\|\rho-\sigma\|_1=\tfrac12\,\mathrm{Tr}\,|\rho-\sigma| \qquad (1)

where: ρ,σ\rho,\sigma — density operators (positive semidefinite, unit trace) on a Hilbert space H\mathcal{H}; ∥A∥1=TrA†A\|A\|_1=\mathrm{Tr}\sqrt{A^\dagger A} — trace norm, the sum of singular values; ∣A∣=A†A|A|=\sqrt{A^\dagger A} — for Hermitian AA, same eigenvectors with absolute-valued eigenvalues. Reading: TT measures how much probability weight separates the states, capped at 11 for perfectly distinguishable ones. Because ρ−σ\rho-\sigma is Hermitian and traceless, the 1-norm collapses to a one-sided sum — the computational workhorse and the key to the measurement.

T(ρ,σ)=∑i:λi>0λi=Tr[Π+(ρ−σ)](2)T(\rho,\sigma)=\sum_{i:\lambda_i>0}\lambda_i = \mathrm{Tr}\big[\Pi_+(\rho-\sigma)\big] \qquad (2)

where: λi\lambda_i — eigenvalues of the traceless Hermitian ρ−σ\rho-\sigma; Π+\Pi_+ — projector onto its strictly positive eigenspace. Reading: diagonalize ρ−σ\rho-\sigma and add only the positive eigenvalues; that same eigenspace defines the winning measurement. For a single qubit this becomes pure geometry. In Bloch form ρ=12(1+r⃗ρ⋅σ⃗)\rho=\tfrac12(\mathbb{1}+\vec r_\rho\cdot\vec\sigma), the difference is 12(r⃗ρ−r⃗σ)⋅σ⃗\tfrac12(\vec r_\rho-\vec r_\sigma)\cdot\vec\sigma, whose eigenvalues are ±12∣r⃗ρ−r⃗σ∣\pm\tfrac12|\vec r_\rho-\vec r_\sigma|.

T(ρ,σ)=12 ∣r⃗ρ−r⃗σ∣(3)T(\rho,\sigma)=\tfrac12\,\big|\vec r_\rho-\vec r_\sigma\big| \qquad (3)

where: r⃗ρ,r⃗σ\vec r_\rho,\vec r_\sigma — Bloch vectors, each of length ≤1\le 1; ∣⋅∣|\cdot| — Euclidean length. Reading: on the Bloch ball TT is half the chord between the points, so antipodal (orthogonal) states give T=1T=1 and coincident states T=0T=0. The Helstrom bound is the operational payoff — the maximum probability of a correct guess for an equiprobable pair — and what makes TT physical rather than a mere norm.

psuccmax⁡=12(1+T(ρ,σ))(4)p_{\mathrm{succ}}^{\max}=\tfrac12\big(1+T(\rho,\sigma)\big) \qquad (4)

where: psuccmax⁡p_{\mathrm{succ}}^{\max} — maximal probability of correctly identifying an equiprobable state by any POVM; T(ρ,σ)T(\rho,\sigma) — the discrimination "bias". Reading: with no information the best is a coin flip (T=0T=0); each unit of trace distance buys exactly that advantage over chance. A variational identity re-expresses TT as a maximization over measurement effects — the form from which contractivity falls out in two lines.

T(ρ,σ)=max⁡0≤E≤1  Tr[E(ρ−σ)](5)T(\rho,\sigma)=\max_{0\le E\le \mathbb{1}}\;\mathrm{Tr}\big[E(\rho-\sigma)\big] \qquad (5)

where: EE — a POVM effect, 0≤E≤10\le E\le\mathbb{1}, the "guess ρ\rho" outcome; 1\mathbb{1} — the identity on H\mathcal{H}. Reading: the best effect is E=Π+E=\Pi_+; no admissible effect collects more than the full positive weight. Optimizing instead over the classical distance of a full POVM's outcome distributions returns the same number — the quantum distance is the largest classical distance any measurement can produce.

T(ρ,σ)=max⁡{Em}  12∑m∣ Tr(Emρ)−Tr(Emσ)∣(6)T(\rho,\sigma)=\max_{\{E_m\}}\;\tfrac12\sum_m\big|\,\mathrm{Tr}(E_m\rho)-\mathrm{Tr}(E_m\sigma)\big| \qquad (6)

where: {Em}\{E_m\} — an arbitrary POVM, Em≥0E_m\ge0, ∑mEm=1\sum_m E_m=\mathbb{1}; Tr(Emρ)\mathrm{Tr}(E_m\rho) — probability of outcome mm given ρ\rho. Reading: measure, take the total variation distance of the two outcome distributions, maximize over the measurement — the ceiling is TT. Contractivity states that this ceiling can only descend under processing.

T(E(ρ),E(σ))≤T(ρ,σ)(7)T\big(\mathcal{E}(\rho),\mathcal{E}(\sigma)\big)\le T(\rho,\sigma) \qquad (7)

where: E\mathcal{E} — any CPTP map. Reading: noise, post-processing, and discarding subsystems are all channels, so none can increase distinguishability — data processing in trace-distance form. Two relations tie TT to the fidelity F(ρ,σ)=∥ρσ∥1F(\rho,\sigma)=\|\sqrt{\rho}\sqrt{\sigma}\|_1: an exact identity for pure states, and the two-sided Fuchs–van de Graaf sandwich in general.

T(∣ψ⟩,∣ϕ⟩)=1−∣⟨ψ∣ϕ⟩∣2,1−F(ρ,σ)≤T(ρ,σ)≤1−F(ρ,σ)2(8)T(|\psi\rangle,|\phi\rangle)=\sqrt{1-|\langle\psi|\phi\rangle|^2},\qquad 1-F(\rho,\sigma)\le T(\rho,\sigma)\le\sqrt{1-F(\rho,\sigma)^2} \qquad (8)

where: ∣ψ⟩,∣ϕ⟩|\psi\rangle,|\phi\rangle — normalized pure states with overlap ⟨ψ∣ϕ⟩\langle\psi|\phi\rangle; FF — Uhlmann fidelity. Reading: upper bound tight for pure states, lower bound when one state is pure — so TT and FF track tightly but are never interchangeable, since only TT is a metric.

Derivation AThe Helstrom bound for two-state discrimination

The goal is to prove the single most important fact about trace distance — that the optimal probability of correctly identifying which of two equiprobable states ρ,σ\rho,\sigma you hold is exactly 12(1+T)\tfrac12(1+T) — and to exhibit the measurement that achieves it. We model a decision rule as a two-outcome POVM, effect E0E_0 announcing "ρ\rho" and E1=1−E0E_1=\mathbb{1}-E_0 announcing "σ\sigma", then write the success probability, reduce it to a single trace, and maximize over admissible effects, at which point the spectral form of TT appears on its own.

Step 1. Write the success probability. With equal priors, a correct call happens when the state is ρ\rho and E0E_0 fires, or the state is σ\sigma and E1E_1 fires, so psucc=12Tr(E0ρ)+12Tr(E1σ)p_{\mathrm{succ}}=\tfrac12\mathrm{Tr}(E_0\rho)+\tfrac12\mathrm{Tr}(E_1\sigma). Why: the prior weights each branch, and Tr(Eρ)\mathrm{Tr}(E\rho) is the Born-rule probability that EE fires on ρ\rho.

Step 2. Collapse to one operator. Substitute E1=1−E0E_1=\mathbb{1}-E_0 and use Tr σ=1\mathrm{Tr}\,\sigma=1 to get psucc=12+12Tr[E0(ρ−σ)]p_{\mathrm{succ}}=\tfrac12+\tfrac12\mathrm{Tr}[E_0(\rho-\sigma)]. Why: maximizing the success probability now means maximizing the single quantity Tr[E0(ρ−σ)]\mathrm{Tr}[E_0(\rho-\sigma)] over effects.

Step 3. Optimize over effects. With ρ−σ=∑iλi∣i⟩⟨i∣\rho-\sigma=\sum_i\lambda_i|i\rangle\langle i|, Tr[E0(ρ−σ)]=∑iλi⟨i∣E0∣i⟩\mathrm{Tr}[E_0(\rho-\sigma)]=\sum_i\lambda_i\langle i|E_0|i\rangle; since 0≤⟨i∣E0∣i⟩≤10\le\langle i|E_0|i\rangle\le1 the sum is largest when E0E_0 puts full weight on positive λi\lambda_i and none on negative ones, giving ∑λi>0λi\sum_{\lambda_i>0}\lambda_i. Why: each eigenvalue is scaled by a dial in [0,1][0,1] — turn the positive dials up, the negative down, and the optimal measurement is born.

Step 4. Name the measurement. The maximizer is E0=Π+=∑λi>0∣i⟩⟨i∣E_0=\Pi_+=\sum_{\lambda_i>0}|i\rangle\langle i|, the projector onto the positive part of ρ−σ\rho-\sigma — the Helstrom measurement. Why: the derivation hands you the explicit optimal apparatus, so the result is constructive, not just an existence claim.

Step 5. Insert the spectral identity. Since T=∑λi>0λiT=\sum_{\lambda_i>0}\lambda_i, psucc=12+12T=12(1+T)p_{\mathrm{succ}}=\tfrac12+\tfrac12 T=\tfrac12(1+T). Why: this bridges the abstract norm to a laboratory probability, certifying TT as the operational bias. *Verified: orthogonal states give T=1T=1 so psucc=1p_{\mathrm{succ}}=1 (always correct); identical states give T=0T=0 so psucc=12p_{\mathrm{succ}}=\tfrac12 (a coin flip). Matches Nielsen & Chuang Box 9.2 and Helstrom (1976).*

Go deeper · Derivation B · 1 minTrace distance is the largest classical distance any measurement can make
Derivation BTrace distance is the largest classical distance any measurement can make

The goal is to establish that trace distance is the operational ceiling on all quantum distances: for any POVM you apply, the classical total variation distance between the resulting outcome distributions is at most T(ρ,σ)T(\rho,\sigma), and one measurement hits that ceiling exactly — which is what justifies calling TT "the" distinguishability. We fix an arbitrary POVM {Em}\{E_m\}, form the distributions pm=Tr(Emρ)p_m=\mathrm{Tr}(E_m\rho), qm=Tr(Emσ)q_m=\mathrm{Tr}(E_m\sigma), bound their distance using the Jordan decomposition of ρ−σ\rho-\sigma, and show the two-outcome Helstrom POVM saturates it.

Step 1. Split the difference operator. Write Δ=ρ−σ=Δ+−Δ−\Delta=\rho-\sigma=\Delta_+-\Delta_-, its Jordan decomposition into orthogonal positive and negative parts Δ±≥0\Delta_\pm\ge0. For any effect, ∣Tr[EmΔ]∣≤Tr[EmΔ+]+Tr[EmΔ−]|\mathrm{Tr}[E_m\Delta]|\le\mathrm{Tr}[E_m\Delta_+]+\mathrm{Tr}[E_m\Delta_-] since the parts have orthogonal supports and Em≥0E_m\ge0. Why: separating where ρ\rho wins from where σ\sigma wins turns an absolute value into a sum of two nonnegative traces the POVM can close up.

Step 2. Sum over outcomes. Because ∑mEm=1\sum_m E_m=\mathbb{1}, summing gives 12∑m∣pm−qm∣≤12(Tr Δ++Tr Δ−)\tfrac12\sum_m|p_m-q_m|\le\tfrac12(\mathrm{Tr}\,\Delta_+ +\mathrm{Tr}\,\Delta_-), and Tr Δ=0\mathrm{Tr}\,\Delta=0 forces Tr Δ+=Tr Δ−=T\mathrm{Tr}\,\Delta_+=\mathrm{Tr}\,\Delta_-=T, so the right side is T(ρ,σ)T(\rho,\sigma). Why: completeness collapses the outcome sum onto Δ±\Delta_\pm, and tracelessness makes both weights TT — the ceiling.

Step 3. Achieve it. The two-outcome POVM {Π+,Π−}\{\Pi_+,\Pi_-\} from the eigenspaces of Δ\Delta gives 12(∣Tr Π+Δ∣+∣Tr Π−Δ∣)=12(T+T)=T\tfrac12(|\mathrm{Tr}\,\Pi_+\Delta|+|\mathrm{Tr}\,\Pi_-\Delta|)=\tfrac12(T+T)=T, meeting the bound with equality. Why: the same Helstrom projectors that solve discrimination also generate the maximally separated classical distributions. *Verified: when [ρ,σ]=0[\rho,\sigma]=0 the states share an eigenbasis and this reduces to the classical total variation distance 12∑i∣pi−qi∣\tfrac12\sum_i|p_i-q_i|, as in Wilde Ch. 9.*

Go deeper · Derivation C · 1 minContractivity of trace distance under quantum channels
Derivation CContractivity of trace distance under quantum channels

The goal is to prove that no quantum channel can increase trace distance, T(Eρ,Eσ)≤T(ρ,σ)T(\mathcal{E}\rho,\mathcal{E}\sigma)\le T(\rho,\sigma) — the trace-distance form of the data-processing inequality and the reason distinguishability behaves like a consumable resource. The route uses the variational form (5) plus two facts about the adjoint map E†\mathcal{E}^\dagger: it is positive and unital. We take the optimal effect for the output states and pull it back through the adjoint to an admissible input effect, where the input maximum must dominate.

Step 1. Start at the output optimum. Let E∗E^\ast achieve the max in (5) for the output pair, so T(Eρ,Eσ)=Tr[E∗ E(ρ−σ)]T(\mathcal{E}\rho,\mathcal{E}\sigma)=\mathrm{Tr}[E^\ast\,\mathcal{E}(\rho-\sigma)]. Why: the variational form lets us move the optimization back to the input space, where the answer is already T(ρ,σ)T(\rho,\sigma).

Step 2. Pass to the Heisenberg picture. By definition of the adjoint, Tr[E∗ E(X)]=Tr[E†(E∗) X]\mathrm{Tr}[E^\ast\,\mathcal{E}(X)]=\mathrm{Tr}[\mathcal{E}^\dagger(E^\ast)\,X], so T(Eρ,Eσ)=Tr[E†(E∗)(ρ−σ)]T(\mathcal{E}\rho,\mathcal{E}\sigma)=\mathrm{Tr}[\mathcal{E}^\dagger(E^\ast)(\rho-\sigma)]. Why: shifting the channel off the states and onto the effect turns an output quantity into an input one without changing its value.

Step 3. Check admissibility. Because E\mathcal{E} is CPTP, its adjoint E†\mathcal{E}^\dagger is positive and unital (E†(1)=1\mathcal{E}^\dagger(\mathbb{1})=\mathbb{1}), so 0≤E∗≤10\le E^\ast\le\mathbb{1} maps to 0≤E†(E∗)≤10\le\mathcal{E}^\dagger(E^\ast)\le\mathbb{1}. Why: this is the crux — trace preservation of E\mathcal{E} is exactly unitality of E†\mathcal{E}^\dagger, guaranteeing the pulled-back operator is a legitimate effect.

Step 4. Invoke the input maximum. Since E†(E∗)\mathcal{E}^\dagger(E^\ast) is one admissible input effect, Tr[E†(E∗)(ρ−σ)]≤max⁡0≤E≤1Tr[E(ρ−σ)]=T(ρ,σ)\mathrm{Tr}[\mathcal{E}^\dagger(E^\ast)(\rho-\sigma)]\le\max_{0\le E\le\mathbb{1}}\mathrm{Tr}[E(\rho-\sigma)]=T(\rho,\sigma). Why: the output distance was realized by an input effect, and no input effect exceeds TT. *Verified: taking E\mathcal{E} to be the partial trace recovers monotonicity under discarding subsystems, and the depolarizing channel gives the explicit contraction T→(1−p)TT\to(1-p)T for qubits (Watrous Ch. 3, Wilde Ch. 9).*

Limiting cases

CaseResultMeaning
Commuting / classical statesT(ρ,σ)=12∑i∣pi−qi∣T(\rho,\sigma)=\tfrac12\sum_i\lvert p_i-q_i\rvertWith [ρ,σ]=0[\rho,\sigma]=0 they share an eigenbasis with eigenvalues pi,qip_i,q_i; TT becomes the classical total variation distance and the quantum problem reduces to a classical hypothesis test.
Pure statesT=1−∣⟨ψ∣ϕ⟩∣2T=\sqrt{1-\lvert\langle\psi\vert\phi\rangle\rvert^2}Orthogonal pure states give T=1T=1 (one measurement identifies them with certainty); identical states give T=0T=0.
Single qubitT=12∣r⃗ρ−r⃗σ∣T=\tfrac12\lvert\vec r_\rho-\vec r_\sigma\rvertHalf the Bloch-ball chord; the maximally mixed center is at distance 12∣r⃗∣\tfrac12\lvert\vec r\rvert from a state of purity ∣r⃗∣\lvert\vec r\rvert.
Fidelity bounds (Fuchs–van de Graaf)1−F≤T≤1−F21-F\le T\le\sqrt{1-F^2}With F=∥ρσ∥1F=\lVert\sqrt{\rho}\sqrt{\sigma}\rVert_1; upper bound tight for pure states, lower bound tight when one state is pure.
Pinsker lower boundD(ρ∥σ)≥12ln⁡2∥ρ−σ∥12=2ln⁡2 T2D(\rho\Vert\sigma)\ge\tfrac{1}{2\ln2}\lVert\rho-\sigma\rVert_1^2=\tfrac{2}{\ln2}\,T^2Quantum Pinsker inequality (relative entropy in bits): a small trace distance forces a small relative entropy, though not conversely.

Worked examples

1 · Distinguishing the two most-used qubit states (fully worked)

You are handed, with equal prior, either ∣0⟩|0\rangle or ∣+⟩=12(∣0⟩+∣1⟩)|+\rangle=\tfrac{1}{\sqrt2}(|0\rangle+|1\rangle). Find the trace distance, the best probability of guessing correctly, and the measurement that achieves it; cross-check against the pure-state formula and Bloch geometry.

Step 1. Form the difference operator. With ∣0⟩⟨0∣=(1000)|0\rangle\langle0|=\begin{pmatrix}1&0\\0&0\end{pmatrix} and ∣+⟩⟨+∣=12(1111)|+\rangle\langle+|=\tfrac12\begin{pmatrix}1&1\\1&1\end{pmatrix}, Δ=ρ−σ=(0.5−0.5−0.5−0.5)\Delta=\rho-\sigma=\begin{pmatrix}0.5&-0.5\\-0.5&-0.5\end{pmatrix} — Hermitian and traceless, as required.

Step 2. Diagonalize. Since Tr Δ=0\mathrm{Tr}\,\Delta=0 and det⁡Δ=(0.5)(−0.5)−(−0.5)2=−0.5\det\Delta=(0.5)(-0.5)-(-0.5)^2=-0.5, the eigenvalues are λ=±0.5=±0.7071\lambda=\pm\sqrt{0.5}=\pm0.7071, so T=∑λ>0λ=0.7071=1/2\mathbf{T=\sum_{\lambda>0}\lambda=0.7071=1/\sqrt2}.

Step 3. Apply the Helstrom bound. psucc=12(1+T)=12(1+0.7071)=0.8536p_{\mathrm{succ}}=\tfrac12(1+T)=\tfrac12(1+0.7071)=\mathbf{0.8536}, achieved by projecting onto the +0.7071+0.7071 eigenvector of Δ\Delta and guessing ∣0⟩|0\rangle there, ∣+⟩|+\rangle otherwise — physically a spin measurement at 22.5∘22.5^\circ, bisecting the two Bloch directions.

Step 4. Cross-check. Pure-state formula: overlap ∣⟨0∣+⟩∣2=1/2|\langle0|+\rangle|^2=1/2, so T=1−1/2=1/2T=\sqrt{1-1/2}=1/\sqrt2. Bloch geometry: r⃗0=(0,0,1)\vec r_0=(0,0,1) and r⃗+=(1,0,0)\vec r_+=(1,0,0) are 90∘90^\circ apart, separation 2\sqrt2, so T=122=1/2T=\tfrac12\sqrt2=1/\sqrt2. All three routes give T=0.707T=0.707, psucc=0.854p_{\mathrm{succ}}=0.854.

2 · How noise erodes distinguishability (fully worked)

The distinguishable states ∣0⟩,∣1⟩|0\rangle,|1\rangle are each sent through a depolarizing channel E(ρ)=(1−p)ρ+p 12\mathcal{E}(\rho)=(1-p)\rho+p\,\tfrac{\mathbb{1}}{2} with p=0.25p=0.25. How far does the best discrimination probability drop?

Step 1. Before the channel. ∣0⟩,∣1⟩|0\rangle,|1\rangle are orthogonal, Bloch vectors (0,0,±1)(0,0,\pm1), so T=12(2)=1T=\tfrac12(2)=\mathbf{1} and psucc=1p_{\mathrm{succ}}=\mathbf{1} — a single ZZ-measurement is always right.

Step 2. Push through the channel. Depolarizing shrinks every Bloch vector, r⃗→(1−p)r⃗\vec r\to(1-p)\vec r, so the points move to (0,0,±0.75)(0,0,\pm0.75) and their separation falls from 22 to 1.51.5.

Step 3. Recompute. T′=12(1.5)=0.75T'=\tfrac12(1.5)=\mathbf{0.75}, exactly (1−p)×1(1-p)\times1 as the qubit contraction predicts, and psucc′=12(1+0.75)=0.875p_{\mathrm{succ}}'=\tfrac12(1+0.75)=\mathbf{0.875}. The noise has turned a certainty into a 7/87/8 guess, and no measurement recovers the lost distinguishability — the channel already spent it. T drops 1→0.751\to0.75, psuccp_{\mathrm{succ}} drops 1→0.8751\to0.875.

3 · The classical sanity check

Two commuting states have eigenvalue distributions p=(0.6,0.4)p=(0.6,0.4) and q=(0.5,0.5)q=(0.5,0.5) in a shared basis. What is their trace distance? Because they commute, TT is the classical total variation distance: T=12(∣0.6−0.5∣+∣0.4−0.5∣)=12(0.1+0.1)=0.1T=\tfrac12(|0.6-0.5|+|0.4-0.5|)=\tfrac12(0.1+0.1)=\mathbf{0.1}. The best guess then succeeds with probability 12(1+0.1)=0.55\tfrac12(1+0.1)=0.55 — barely better than chance, as expected for two nearly identical distributions.

Physics of the portrait

The trace distance is the crystallization of a concrete engineering question — how well can a receiver tell two signals apart — asked in the quantum regime for the first time by a radar theorist who saw that at the single-photon level quantum mechanics, not noise, sets the answer.

Origins · 1 minCarl W. Helstrom and how this was discovered
CW
Carl W. Helstrom
1925–2013 · quantum detection and estimation theory

Carl Wilhelm Helstrom was an American electrical engineer and physicist at the University of California, San Diego, whose 1967 paper "Detection Theory and Quantum Mechanics" founded quantum detection theory. Trained in the classical statistics of signal detection — deciding whether a radar echo is a target or just noise — Helstrom asked what changes when the signals are quantum states so weak that the limit is Heisenberg's, not the receiver's electronics. His answer, the Helstrom bound, is that the minimum error probability for distinguishing two states with priors p0,p1p_0,p_1 is 12(1−∥p0ρ0−p1ρ1∥1)\tfrac12(1-\|p_0\rho_0-p_1\rho_1\|_1), which for equal priors is exactly 12(1−T)\tfrac12(1-T) — the trace distance made operational — with optimal measurement the projection onto the positive part of the weighted difference operator. He gathered this and his quantum estimation work into the 1976 monograph Quantum Detection and Estimation Theory, a founding text of quantum information science written years before the field had that name. His perspective — that information-theoretic limits are ultimately statements about state distinguishability — runs straight through the modern theory, from the Holevo bound to composable cryptography.

History
1927
— von Neumann introduces the density operator and the projective theory of quantum measurement
1967
— Helstrom solves the minimum-error discrimination of two quantum states and identifies the optimal measurement
1973
— Holevo bounds the accessible information of a quantum ensemble, tying distinguishability to capacity
1976
— Helstrom's Quantum Detection and Estimation Theory codifies the field
1999
— Fuchs and van de Graaf establish the two-sided inequalities between trace distance and fidelity

Going further

Texts. Nielsen & Chuang, Quantum Computation and Quantum Information, §9.2 (distance measures) · Wilde, Quantum Information Theory, 2nd ed., Ch. 9 · Watrous, The Theory of Quantum Information, Ch. 3 · Helstrom, Quantum Detection and Estimation Theory (1976) · Fuchs & van de Graaf, IEEE Trans. Inf. Theory 45, 1216 (1999).

Provenance. A verifier should recompute: the spectral identity T=∑λi>0λiT=\sum_{\lambda_i>0}\lambda_i from tracelessness (N&C §9.2.1); the Helstrom bound and optimal projector (N&C Box 9.2, Helstrom 1967); the variational and maximal-classical-distance forms and contractivity via the positive unital adjoint (Wilde Ch. 9, Watrous Ch. 3); the qubit relation T=12∣r⃗ρ−r⃗σ∣T=\tfrac12|\vec r_\rho-\vec r_\sigma|; the pure-state formula and Fuchs–van de Graaf bounds (Fuchs & van de Graaf 1999, N&C §9.2.3); and quantum Pinsker D≥2ln⁡2T2D\ge\tfrac{2}{\ln2}T^2 (Wilde). The worked numbers — T=1/2T=1/\sqrt2, psucc=0.854p_{\mathrm{succ}}=0.854; depolarizing T:1→0.75T:1\to0.75; classical T=0.1T=0.1 — follow directly.

Prerequisites. Von Neumann Entropy · Shannon Entropy and Classical Information Measures · Quantum Channels and Completely Positive Trace-Preserving Maps. Leads to. Quantum Relative Entropy · Data-Processing Inequality · Holevo Bound.

Ask the tutor: ask me to work any two density matrices you like all the way to their Helstrom success probability, or to show why post-selection is the one operation that can pull two states apart.

Why is there a factor of one half in the definition?

The one half normalizes the range to [0,1][0,1] and makes the quantum distance agree with the classical one. The trace norm ∥ρ−σ∥1\|\rho-\sigma\|_1 equals 2P2P, twice the total positive eigenvalue weight, because the traceless operator has equal positive and negative parts; halving recovers PP, the quantity with the clean meaning psucc=12(1+T)p_{\mathrm{succ}}=\tfrac12(1+T). It also matches the classical total variation distance 12∑i∣pi−qi∣\tfrac12\sum_i|p_i-q_i|. Some references drop the factor and use ∥ρ−σ∥1\|\rho-\sigma\|_1 directly, so a formula that looks "off by two" is usually a convention clash.

How is trace distance different from fidelity?

Both measure closeness of states but are not interchangeable. Trace distance is a genuine metric — symmetric, zero only when states are equal, obeying the triangle inequality — so errors accumulate additively, indispensable for composing error bounds. Fidelity F=∥ρσ∥1F=\|\sqrt\rho\sqrt\sigma\|_1 runs the other way (F=1F=1 for equal states) and is not a metric. The two are tied by 1−F≤T≤1−F21-F\le T\le\sqrt{1-F^2}, so knowing one pins the other to a range but never a single value. Rule of thumb: use fidelity for a smooth, easily computed overlap; use trace distance when you need a true metric with operational meaning.

What exactly is the Helstrom measurement, physically?

It is the two-outcome measurement whose "yes" projector is Π+\Pi_+, the projection onto the subspace where ρ−σ\rho-\sigma is positive. It asks the most informative binary question — is the state where ρ\rho outweighs σ\sigma, or the reverse? — and you guess accordingly. For qubits this is a projective spin measurement along the Bloch axis bisecting the two states; for ∣0⟩|0\rangle versus ∣+⟩|+\rangle it sits at 22.5∘22.5^\circ between the ZZ and XX axes. Crucially it is a fixed measurement you can build, and for two states it is provably optimal on a single copy — no POVM beats it.

Can a quantum channel ever increase trace distance?

No — not a genuine channel. Every CPTP map is contractive, T(Eρ,Eσ)≤T(ρ,σ)T(\mathcal{E}\rho,\mathcal{E}\sigma)\le T(\rho,\sigma): processing cannot create distinguishability that was not there, and unitaries preserve TT exactly. What people mistake for an increase is post-selection: conditioning on a rare outcome and renormalizing the survivor can grow the distance. But conditioning throws away the failed runs, so it is not trace-preserving and contractivity does not apply. Account for the discarded outcomes and monotonicity is restored.

How does trace distance connect to relative entropy?

Through the quantum Pinsker inequality, D(ρ∥σ)≥2ln⁡2T(ρ,σ)2D(\rho\|\sigma)\ge\tfrac{2}{\ln2}T(\rho,\sigma)^2 (relative entropy in bits). This is one-way: a small relative entropy forces a small trace distance. The converse fails badly — relative entropy can be infinite (when the support of ρ\rho is not contained in that of σ\sigma) while TT stays below one. That asymmetry is why both are needed: TT is the bounded, symmetric, operational metric, while relative entropy is the unbounded, asymmetric quantity governing asymptotic hypothesis testing through quantum Stein's lemma. Pinsker is the bridge between them.

Why do cryptographers insist on trace distance rather than fidelity for security?

Because security has to compose. A key feeds encryption, authentication, and further protocols, so a definition is meaningful only if a key certified "good" stays good in every larger context. Trace distance delivers this: if the real key is within ϵ\epsilon in trace distance of an ideal key, then by the Helstrom bound any distinguisher, including the eavesdropper's best strategy, tells real from ideal with advantage at most ϵ\epsilon, and the triangle inequality lets these add across composed steps. Fidelity has no triangle inequality and no clean distinguishing-advantage reading, so it would not compose — which is why the modern "ϵ\epsilon-security" framework for QKD is written entirely in trace distance.

Before this

Where this leads