Resource Theory of Coherence

Superposition is the first thing quantum mechanics teaches and the last it lets you keep. A qubit prepared in 12(∣0⟩+∣1⟩)\tfrac{1}{\sqrt2}(|0\rangle+|1\rangle) carries a definite phase relation between two basis states; leave it exposed to its environment and that phase relation bleeds away, leaving a classical coin. The resource theory of coherence is the accounting system for that phase relation — it treats definite superposition relative to a fixed basis as a commodity that can be measured, spent, transferred, and irreversibly lost, exactly as one treats energy in thermodynamics or entanglement in quantum communication. The construction, due to Baumgratz, Cramer and Plenio in 2014 (building on Åberg's 2006 formalization of superposition), imports the architecture that had already tamed entanglement. Fix a preferred orthonormal basis {∣i⟩}\{|i\rangle\} — computational, energy, or a decoherence-selected pointer basis. Declare the states diagonal in it free (incoherent), because a classical apparatus can prepare them for nothing, and declare the operations that can never manufacture off-diagonal terms free operations. Then any genuine superposition is a resource, and a coherence monotone assigns it a number that free operations can only leave alone or destroy. This article develops the theory to its operational core: the two workhorse measures Cr(ρ)=S(Δρ)−S(ρ)C_r(\rho)=S(\Delta\rho)-S(\rho) and Cℓ1(ρ)=∑i≠j∣ρij∣C_{\ell_1}(\rho)=\sum_{i\neq j}|\rho_{ij}|; the hierarchy of free operation classes and why it is not one class but many; and the crowning result that CrC_r is the physical rate at which pure coherence can be distilled.

Key takeaways
  • The theory fixes a reference basis and calls the diagonal (incoherent) states free; every superposition of basis states is the resource, and the fully dephasing channel Δ is the canonical way to destroy it.
  • Free operations are the incoherent operations (IO): channels with a Kraus decomposition that maps incoherent states to incoherent states branch-by-branch — flanked by looser (MIO) and stricter (SIO, DIO, PIO) classes that are genuinely inequivalent.
  • A coherence monotone vanishes exactly on the incoherent set and never increases under the free operations; the two standard choices are the relative entropy of coherence Cr(ρ)=S(Δρ)−S(ρ)C_r(\rho)=S(\Delta\rho)-S(\rho) and the ℓ1\ell_1-norm Cℓ1(ρ)=∑i≠j∣ρij∣C_{\ell_1}(\rho)=\sum_{i\neq j}|\rho_{ij}|.
  • The maximally coherent state ∣Ψd⟩=1d∑i∣i⟩|\Psi_d\rangle=\tfrac{1}{\sqrt d}\sum_i|i\rangle is the universal unit of the resource, with Cr=log⁡2dC_r=\log_2 d and Cℓ1=d−1C_{\ell_1}=d-1; it can generate any state of the same dimension by IO alone.
  • CrC_r carries an operational meaning: it equals the distillable coherence, the optimal asymptotic rate of extracting maximally coherent qubits under IO — the coherence analogue of the entropy of entanglement.
  • Key idea: pick a basis, call its diagonal states free and the operations that never create off-diagonal terms free; then coherence becomes a bona fide resource whose amount is captured by basis-dependent monotones, the foremost being Cr(ρ)=S(Δρ)−S(ρ)C_r(\rho)=S(\Delta\rho)-S(\rho), equal to the asymptotic distillable coherence.

Where this lives today

Coherence-as-resource has moved from formalism to laboratory bookkeeping. Photonic and NMR experiments now directly witness and quantify Cℓ1C_{\ell_1} and CrC_r by interferometric readout of the off-diagonals. The theory is the natural currency of quantum thermodynamics: with the reference basis taken as the energy eigenbasis, coherence between energy levels is a resource thermal operations cannot create for free, setting limits on quantum heat engines, batteries, and the work extractable from a coherently driven qubit — coherent strokes can beat the incoherent bound, but only by spending a budget the theory tracks exactly. In metrology, a probe's coherence lower-bounds its quantum Fisher information for phase estimation. The live open problems are structural: the theory is asymptotically reversible under MIO but not under IO, where bound coherence — costing more to create than can ever be distilled back — genuinely exists, mirroring bound entanglement. Pinning down the exact gaps between MIO, IO, DIO and SIO remains an active frontier a decade on.

The physical picture

Every resource theory partitions the world into the free and the precious, and the whole content of coherence theory is encoded in the density matrix written in the chosen basis. The diagonal entries ρii\rho_{ii} are populations — the classical probabilities of finding the system in each basis state, an interferometer's which-path information. The off-diagonal entries ρij\rho_{ij} are the coherences — the phase relations that make interference fringes visible. A classical apparatus can dial any diagonal for free, so the diagonal is the free part; the off-diagonals are what quantum mechanics adds and decoherence takes away. The dephasing channel Δ erases the off-diagonals while leaving the diagonal untouched — the natural yardstick, since coherence is what Δ destroys.

Figure
The density matrix in the reference basis splits into a free diagonal (populations, which-path information) and a resourceful off-diagonal (coherences, fringe visibility). The dephasing channel Δ deletes the off-diagonals, projecting ρ onto the incoherent set, and the entropy it releases is the relative entropy of coherence C_r = S(Δρ) − S(ρ)
Diagonal = free populations · off-diagonal = the coherence resource · Δ deletes it state ρ (has coherence) ρ₀₀ ρ₀₁ ρ₁₀ ρ₁₁ populations (free) coherences (resource) Δ dephase Δρ (incoherent) ρ₀₀ 0 0 ρ₁₁ off-diagonals erased entropy rises by C_r C_r(ρ) = S(Δρ) − S(ρ) ≥ 0

The parallel to entanglement theory is a near-isomorphism: incoherent states play the role of separable states, incoherent operations that of LOCC, the maximally coherent state ∣Ψd⟩|\Psi_d\rangle that of the singlet, and distillable coherence mirrors distillable entanglement. The second figure makes the payoff concrete: the two standard measures agree on the qubit unit but diverge as the dimension grows, a reminder that a resource can be scored on more than one honest scale.

Figure
The two standard measures of the maximally coherent state Ψ_d as dimension grows: the relative entropy of coherence C_r = log₂ d in blue rises logarithmically, while the ℓ₁-norm C_ℓ1 = d − 1 in amber rises linearly. They coincide at the qubit unit d = 2 (both equal 1) and split thereafter — at d = 8, C_r = 3 bits but C_ℓ1 = 7
The unit resource |Ψ_d⟩ scored two ways: C_r = log₂ d vs C_ℓ1 = d − 1 0 1 3 5 7 2 4 6 8 dimension d of |Ψ_d⟩ coherence value d = 2: both = 1 C_ℓ1 = 7 C_r = 3 C_r = log₂ d (bits) C_ℓ1 = d − 1

The core equations

The theory is specified by three ingredients — free states, free operations, monotones — closed by an operational theorem, taken here in that order. The first fixes what "free" means for states: a density matrix is incoherent exactly when it is diagonal in the reference basis.

I={ ρ=∑ipi ∣i⟩⟨i∣ : pi≥0, ∑ipi=1 }(1)\mathcal{I}=\Big\{\,\rho=\sum_{i} p_i\,|i\rangle\langle i|\ :\ p_i\geq0,\ \textstyle\sum_i p_i=1\,\Big\} \qquad (1)

where: {∣i⟩}\{|i\rangle\} — the fixed reference orthonormal basis of Cd\mathbb{C}^d, an input to the theory; pip_i — classical probabilities, the eigenvalues of the diagonal state; I\mathcal{I} — the convex set of incoherent states. Reading: I\mathcal{I} is the free set; its convexity matters (mixing free states must stay free), and every resource number below vanishes precisely on it. The dephasing channel projects any state onto I\mathcal{I} by deleting off-diagonals — the canonical free way to destroy coherence and the reference map against which it is measured.

Δ(ρ)=∑i∣i⟩⟨i∣ ρ ∣i⟩⟨i∣=∑iρii ∣i⟩⟨i∣(2)\Delta(\rho)=\sum_i |i\rangle\langle i|\,\rho\,|i\rangle\langle i|=\sum_i \rho_{ii}\,|i\rangle\langle i| \qquad (2)

where: Δ\Delta — the completely dephasing channel; ρii=⟨i∣ρ∣i⟩\rho_{ii}=\langle i|\rho|i\rangle — a diagonal matrix element. Reading: Δ is CPTP, idempotent (Δ∘Δ=Δ\Delta\circ\Delta=\Delta), and fixes exactly the incoherent states (Δρ=ρ  ⟺  ρ∈I\Delta\rho=\rho \iff \rho\in\mathcal{I}) — the coherence-theory analogue of the partial trace, the free operation that annihilates the resource. Free operations are subtler, because "cannot create coherence" can be read at the level of the channel or of its individual outcomes. The original and central choice is the incoherent operations.

E(ρ)=∑nKnρKn†,Kn I Kn†⊆I  ∀n,∑nKn†Kn=I(3)\mathcal{E}(\rho)=\sum_n K_n\rho K_n^\dagger,\qquad K_n\,\mathcal{I}\,K_n^\dagger\subseteq\mathcal{I}\ \ \forall n,\qquad \sum_n K_n^\dagger K_n=I \qquad (3)

where: {Kn}\{K_n\} — a Kraus decomposition of E\mathcal{E}; KnρKn†K_n\rho K_n^\dagger — the unnormalized post-measurement state of outcome nn, required incoherent whenever the input was. Reading: the demand is branch-by-branch — no single selective outcome may create coherence — strictly stronger than merely mapping I→I\mathcal{I}\to\mathcal{I} overall (the weaker demand defining the largest class, MIO). Incoherence is a property of a decomposition, not of the channel alone. The first workhorse measure quantifies coherence as distinguishability from the free set, collapsing to a simple closed form.

Cr(ρ)=min⁡σ∈ID(ρ ∥ σ)=S(Δρ)−S(ρ)(4)C_r(\rho)=\min_{\sigma\in\mathcal{I}}D(\rho\,\|\,\sigma)=S(\Delta\rho)-S(\rho) \qquad (4)

where: D(ρ∥σ)=Tr ρ(log⁡ρ−log⁡σ)D(\rho\|\sigma)=\mathrm{Tr}\,\rho(\log\rho-\log\sigma) — quantum relative entropy; S(ρ)=−Tr ρlog⁡ρS(\rho)=-\mathrm{Tr}\,\rho\log\rho — von Neumann entropy (log base 2, so the unit is bits). Reading: the closest free state to ρ\rho is always its own dephasing Δρ\Delta\rho, so the coherence is exactly the entropy Δ releases (Derivation A). The second workhorse measure is geometric, not entropic: add up the moduli of the off-diagonal entries.

Cℓ1(ρ)=∑i≠j∣ρij∣=(∑i,j∣ρij∣)−1(5)C_{\ell_1}(\rho)=\sum_{i\neq j}|\rho_{ij}|=\Big(\sum_{i,j}|\rho_{ij}|\Big)-1 \qquad (5)

where: ρij=⟨i∣ρ∣j⟩\rho_{ij}=\langle i|\rho|j\rangle — an off-diagonal matrix element in the reference basis. Reading: Cℓ1C_{\ell_1} is zero for diagonal states and trivial to compute, the everyday diagnostic; it is an IO monotone but, unlike CrC_r, is not a monotone under all of MIO — a first warning that measure and operation class must be quoted together. The maximally coherent state is the universal unit — the state from which every other state of the same dimension can be reached for free.

∣Ψd⟩=1d∑i=0d−1∣i⟩,Cr(∣Ψd⟩)=log⁡2d,Cℓ1(∣Ψd⟩)=d−1(6)|\Psi_d\rangle=\frac{1}{\sqrt d}\sum_{i=0}^{d-1}|i\rangle,\qquad C_r(|\Psi_d\rangle)=\log_2 d,\quad C_{\ell_1}(|\Psi_d\rangle)=d-1 \qquad (6)

where: ∣Ψd⟩|\Psi_d\rangle — the uniform superposition, the golden unit of coherence in dimension dd. Reading: because ∣Ψd⟩→∣ψ⟩|\Psi_d\rangle\to|\psi\rangle by IO for any target, its coherence upper-bounds every state's in that dimension; log⁡2d\log_2 d and d−1d-1 are the ceilings, one per measure. Finally, the identity that lifts CrC_r from definition to physics — the rate at which pure coherence can be distilled.

Cd(ρ)=Cr∞(ρ)=Cr(ρ)(distillable coherence under IO)(7)C_d(\rho)=C_r^{\infty}(\rho)=C_r(\rho)\qquad (\text{distillable coherence under IO}) \qquad (7)

where: Cd(ρ)C_d(\rho) — the optimal asymptotic rate of maximally coherent qubits ∣Ψ2⟩|\Psi_2\rangle extractable per copy of ρ\rho under IO; Cr∞C_r^{\infty} — its regularized value. Reading: Winter and Yang proved this single-letter identity — no regularization needed — establishing CrC_r as the coherence analogue of the entropy of entanglement, the resource's true operational currency.

Derivation AThe closed form of the relative entropy of coherence

The goal is to prove the second equality in equation (4): that the variational quantity min⁡σ∈ID(ρ∥σ)\min_{\sigma\in\mathcal{I}}D(\rho\|\sigma), an optimization over the entire convex set of incoherent states, collapses to the explicit S(Δρ)−S(ρ)S(\Delta\rho)-S(\rho), with the minimizer being ρ\rho's own dephasing Δρ\Delta\rho. This is what makes the relative entropy of coherence computable in closed form rather than as a convex program, and it mirrors the classical fact that the closest diagonal distribution to a given one, in relative entropy, is that distribution itself.

Step 1. Expand: D(ρ∥σ)=−S(ρ)−Tr(ρlog⁡σ)D(\rho\|\sigma)=-S(\rho)-\mathrm{Tr}(\rho\log\sigma). Why: the first term is fixed by ρ\rho, so the entire dependence on the free state σ\sigma lives in the cross term −Tr(ρlog⁡σ)-\mathrm{Tr}(\rho\log\sigma), and minimizing DD means minimizing that term alone.

Step 2. For an incoherent σ=∑iqi∣i⟩⟨i∣\sigma=\sum_i q_i|i\rangle\langle i|, log⁡σ\log\sigma is diagonal, so Tr(ρlog⁡σ)=∑iρiilog⁡qi=Tr(Δρ log⁡σ)\mathrm{Tr}(\rho\log\sigma)=\sum_i\rho_{ii}\log q_i=\mathrm{Tr}(\Delta\rho\,\log\sigma). Why: only the diagonal entries ρii\rho_{ii} ever meet the diagonal log⁡σ\log\sigma; the coherences are invisible to any incoherent σ\sigma, so the optimal σ\sigma can only see Δρ\Delta\rho.

Step 3. Minimize the cross term. Writing pi=ρiip_i=\rho_{ii}, Gibbs' inequality −∑ipilog⁡qi≥H(p)-\sum_i p_i\log q_i\geq H(p) gives the optimum at qi=piq_i=p_i, i.e. σ=Δρ\sigma=\Delta\rho, with value S(Δρ)S(\Delta\rho). Why: cross-entropy is minimized by the true distribution, so the best free state is the one whose diagonal matches ρ\rho's exactly.

Step 4. Assemble: Cr(ρ)=S(Δρ)−S(ρ)C_r(\rho)=S(\Delta\rho)-S(\rho), non-negative since dephasing can only raise entropy, and vanishing iff ρ=Δρ\rho=\Delta\rho, i.e. iff ρ∈I\rho\in\mathcal{I}. Why: this delivers the closed form and verifies faithfulness (axiom C1) at once — CrC_r is zero exactly on the free set. *Verified: for ρ=12(I+0.8 σx)\rho=\tfrac12(I+0.8\,\sigma_x), S(ρ)=0.4690S(\rho)=0.4690 bits and S(Δρ)=1S(\Delta\rho)=1 bit give Cr=0.531C_r=0.531 bits; the minimizer Δρ=diag(12,12)\Delta\rho=\mathrm{diag}(\tfrac12,\tfrac12) was checked against Gibbs' inequality.*

Go deeper · Derivation B · 1 minWhy relative-entropy coherence cannot increase under incoherent operations
Derivation BWhy relative-entropy coherence cannot increase under incoherent operations

The goal is to establish the monotonicity axiom C2 for CrC_r: that no incoherent operation E\mathcal{E} can raise it, Cr(E(ρ))≤Cr(ρ)C_r(\mathcal{E}(\rho))\leq C_r(\rho). This is the property that earns CrC_r the title of monotone — a quantity free operations can only waste. The proof rests on two facts in hand: the closed form of Derivation A, and the data-processing inequality, which says relative entropy is contractive under every CPTP map, so no channel can make two states more distinguishable than they already were.

Step 1. Use the closed form as a variational bound. From Derivation A, Cr(ρ)=D(ρ∥Δρ)C_r(\rho)=D(\rho\|\Delta\rho), and for the output, Cr(E(ρ))=min⁡σ∈ID(E(ρ)∥σ)C_r(\mathcal{E}(\rho))=\min_{\sigma\in\mathcal{I}}D(\mathcal{E}(\rho)\|\sigma). Why: writing the input as an exact relative entropy and the output as a minimum lets us bound the output by any convenient incoherent trial state.

Step 2. Choose the trial state σ=E(Δρ)\sigma=\mathcal{E}(\Delta\rho). Because E\mathcal{E} is incoherent it maps the incoherent Δρ\Delta\rho to an incoherent state, so this σ\sigma is admissible, giving Cr(E(ρ))≤D(E(ρ)∥E(Δρ))C_r(\mathcal{E}(\rho))\leq D(\mathcal{E}(\rho)\|\mathcal{E}(\Delta\rho)). Why: the defining property of IO — preserving the free set — is exactly what makes E(Δρ)\mathcal{E}(\Delta\rho) a legal competitor in the minimization.

Step 3. Apply data processing. Relative entropy is contractive under any CPTP map, so D(E(ρ)∥E(Δρ))≤D(ρ∥Δρ)D(\mathcal{E}(\rho)\|\mathcal{E}(\Delta\rho))\leq D(\rho\|\Delta\rho). Why: the same channel on both arguments cannot increase their relative entropy — identical noise can only blur the distinction — the one deep inequality the proof leans on.

Step 4. Chain: Cr(E(ρ))≤D(E(ρ)∥E(Δρ))≤D(ρ∥Δρ)=Cr(ρ)C_r(\mathcal{E}(\rho))\leq D(\mathcal{E}(\rho)\|\mathcal{E}(\Delta\rho))\leq D(\rho\|\Delta\rho)=C_r(\rho). Convexity of CrC_r follows separately from the joint convexity of DD. Why: the combined inequality is axiom C2; with C1 from Derivation A it certifies CrC_r as a full coherence monotone — the abstract requirement is now a theorem. *Verified: the two links use only (i) IO preserving I\mathcal{I} and (ii) the data-processing inequality for quantum relative entropy (Wilde, Ch. 11).*

Go deeper · Derivation C · 1 minCoherence of a pure state equals the entropy of its measurement statistics
Derivation CCoherence of a pure state equals the entropy of its measurement statistics

The goal is to show that for a pure state ∣ψ⟩=∑ici∣i⟩|\psi\rangle=\sum_i c_i|i\rangle the relative entropy of coherence reduces to a classical quantity — the Shannon entropy H({∣ci∣2})H(\{|c_i|^2\}) of the distribution a basis measurement would produce. This closes the conceptual circle: coherence, an intrinsically quantum notion, is measured on pure states by exactly the classical uncertainty of where the state would be found — and the ceiling log⁡2d\log_2 d falls out in one line.

Step 1. The purity term vanishes: a pure state has S(∣ψ⟩⟨ψ∣)=0S(|\psi\rangle\langle\psi|)=0, its density matrix having a single unit eigenvalue. Why: with S(ρ)=0S(\rho)=0, the whole of Cr=S(Δρ)−S(ρ)C_r=S(\Delta\rho)-S(\rho) reduces to the entropy of the dephased state, so only Δ's output matters.

Step 2. Dephase: Δ(∣ψ⟩⟨ψ∣)=∑i∣ci∣2 ∣i⟩⟨i∣\Delta(|\psi\rangle\langle\psi|)=\sum_i|c_i|^2\,|i\rangle\langle i|, since the off-diagonals cicj∗c_ic_j^* are deleted, leaving the diagonal ∣ci∣2|c_i|^2. Why: the moduli-squared amplitudes are the Born-rule probabilities of a basis measurement, so the dephased state is the classical measurement statistics as a density matrix.

Step 3. Take its entropy: S(Δρ)=−∑i∣ci∣2log⁡2∣ci∣2=H({∣ci∣2})S(\Delta\rho)=-\sum_i|c_i|^2\log_2|c_i|^2=H(\{|c_i|^2\}), the entropy of a diagonal state being the Shannon entropy of its diagonal. Why: with the purity term zero, this gives Cr(∣ψ⟩)=H({∣ci∣2})C_r(|\psi\rangle)=H(\{|c_i|^2\}) — pure-state coherence is measurement uncertainty, nothing more.

Step 4. Specialize to the uniform superposition ∣ci∣2=1/d|c_i|^2=1/d: Cr=log⁡2dC_r=\log_2 d, recovering equation (6). Why: the flat distribution maximizes Shannon entropy at log⁡2d\log_2 d, so ∣Ψd⟩|\Psi_d\rangle is literally the maximally coherent pure state, its coherence the ceiling for the whole dimension. *Verified: for ∣Ψ3⟩|\Psi_3\rangle, H({1/3,1/3,1/3})=log⁡23=1.585H(\{1/3,1/3,1/3\})=\log_2 3=1.585 bits, matching equation (6); the pure-state ℓ1\ell_1 value (∑i∣ci∣)2−1=2(\sum_i|c_i|)^2-1=2.*

Limiting cases

The standard measures take instructive values in the regimes that bound the theory:

CaseResultMeaning
Incoherent (diagonal) stateCr(Δρ)=0, Cℓ1(Δρ)=0C_r(\Delta\rho)=0,\ C_{\ell_1}(\Delta\rho)=0Faithfulness — every free state carries zero resource, by construction of the theory.
Maximally mixed state I/dI/dCr=0, Cℓ1=0C_r=0,\ C_{\ell_1}=0Randomness is not coherence: I/dI/d is diagonal in every basis, so it is free despite maximal entropy.
Pure state ∣ψ⟩=∑ici∣i⟩\vert\psi\rangle=\sum_i c_i\vert i\rangleCr=H({∣ci∣2}), Cℓ1=(∑i∣ci∣)2−1C_r=H(\{\vert c_i\vert^2\}),\ C_{\ell_1}=(\sum_i\vert c_i\vert)^2-1CrC_r becomes the Shannon entropy of the measurement statistics; both vanish iff a single ∣ci∣=1\vert c_i\vert=1.
Maximally coherent qudit ∣Ψd⟩\vert\Psi_d\rangleCr=log⁡2d, Cℓ1=d−1C_r=\log_2 d,\ C_{\ell_1}=d-1The unit resource: maximal on both scales and able to generate any dd-dimensional state by IO.
Single qubit ρ=12(I+r⃗⋅σ⃗)\rho=\tfrac12(I+\vec r\cdot\vec\sigma)Cℓ1=rx2+ry2=∣2ρ01∣C_{\ell_1}=\sqrt{r_x^2+r_y^2}=\vert 2\rho_{01}\vertCoherence is the equatorial (Bloch-plane) radius; the zz-component is pure population and contributes nothing.

Worked examples

1 · Coherence of a partially decohered qubit (fully worked)

A qubit is prepared with Bloch vector r⃗=(0.8,0,0)\vec r=(0.8,0,0) — an equatorial state, partway between pure and maximally mixed. Compute Cℓ1C_{\ell_1} and CrC_r in the computational basis, and say what CrC_r means operationally.

Step 1. With ρ=12(I+r⃗⋅σ⃗)\rho=\tfrac12(I+\vec r\cdot\vec\sigma) and r⃗=(0.8,0,0)\vec r=(0.8,0,0), the off-diagonal is ρ01=12(rx−iry)=0.4\rho_{01}=\tfrac12(r_x-ir_y)=0.4, so ρ=(0.50.40.40.5)\rho=\begin{pmatrix}0.5&0.4\\0.4&0.5\end{pmatrix}. Populations equal; all the structure is in the coherence.

Step 2. The ℓ1\ell_1-norm is immediate: Cℓ1=∣ρ01∣+∣ρ10∣=2(0.4)=0.8C_{\ell_1}=|\rho_{01}|+|\rho_{10}|=2(0.4)=\mathbf{0.8}, equivalently rx2+ry2=0.8\sqrt{r_x^2+r_y^2}=0.8 — the fringe visibility a Mach–Zehnder interferometer would read for this state, 80% of the maximum.

Step 3. For CrC_r, the eigenvalues of ρ\rho are 1±∣r⃗∣2={0.9,0.1}\tfrac{1\pm|\vec r|}{2}=\{0.9,0.1\}, so S(ρ)=H(0.9)=−0.9log⁡20.9−0.1log⁡20.1=0.4690S(\rho)=H(0.9)=-0.9\log_2 0.9-0.1\log_2 0.1=0.4690 bits. The dephased state Δρ=diag(0.5,0.5)\Delta\rho=\mathrm{diag}(0.5,0.5) has S(Δρ)=1S(\Delta\rho)=1 bit. Hence Cr=1−0.4690=0.531C_r=1-0.4690=\mathbf{0.531} bits.

Step 4. Interpret. By equation (7), CrC_r is the distillable coherence: from NN copies, incoherent operations distill about 0.531 N0.531\,N perfect maximally coherent qubits ∣Ψ2⟩|\Psi_2\rangle asymptotically. Concretely, 1000 copies yield ≈ 531 units of pure coherence — just over half a "coherbit" apiece. The two measures disagree in absolute value (0.8 versus 0.531): different honest scales, and only CrC_r carries the distillation meaning.

2 · The maximally coherent qutrit as a unit of resource

Evaluate both measures for ∣Ψ3⟩=13(∣0⟩+∣1⟩+∣2⟩)|\Psi_3\rangle=\tfrac{1}{\sqrt3}(|0\rangle+|1\rangle+|2\rangle), and compare its distillation value to that of a qubit unit.

Step 1. Every entry of the density matrix equals 1/31/3. The off-diagonals number 32−3=63^2-3=6, each of modulus 1/31/3, so Cℓ1=6⋅13=2=d−1C_{\ell_1}=6\cdot\tfrac13=\mathbf{2}=d-1, saturating the qutrit ceiling.

Step 2. For CrC_r, the state is pure so S(ρ)=0S(\rho)=0, and Δρ=diag(13,13,13)\Delta\rho=\mathrm{diag}(\tfrac13,\tfrac13,\tfrac13) has S(Δρ)=log⁡23=1.585S(\Delta\rho)=\log_2 3=\mathbf{1.585} bits (Derivation C). So Cr(∣Ψ3⟩)=1.585C_r(|\Psi_3\rangle)=1.585 bits.

Step 3. Compare to the qubit unit ∣Ψ2⟩|\Psi_2\rangle (Cr=1C_r=1 bit). The qutrit carries log⁡23≈1.585\log_2 3\approx1.585 times as much: asymptotically, NN copies of ∣Ψ3⟩|\Psi_3\rangle convert by IO into about 1.585 N1.585\,N copies of ∣Ψ2⟩|\Psi_2\rangle. One maximally coherent qutrit is worth log₂3 ≈ 1.585 maximally coherent qubits — which is why ∣Ψ2⟩|\Psi_2\rangle, the "coherbit," is the standard reference every resource is priced against.

Physics of the portrait

The resource theory of coherence has a clean origin: a single 2014 paper turned a century-old intuition about superposition into a quantitative theory, by an author who had already done the same for entanglement.

Origins · 1 minMartin Plenio and how this was discovered
MP
Martin Plenio
b. 1968, German physicist · quantum information, resource theories of coherence and entanglement

Martin B. Plenio is a German theoretical physicist whose work runs like a spine through modern quantum resource theory. Shaped in the 1990s Imperial College group around Peter Knight, he became one of the architects of entanglement theory: with Vlatko Vedral he introduced the relative entropy of entanglement in 1997, casting a resource as the distance from a state to the free set — the very template later transplanted to coherence. In 2014, with Tillmann Baumgratz and Marcus Cramer, he published "Quantifying Coherence" (Phys. Rev. Lett. 113, 140401), which fixed the axioms C1–C2, proposed the relative-entropy and ℓ1\ell_1-norm measures, and launched coherence as a resource theory in its own right — building on Johan Åberg's 2006 formalization of "quantifying superposition." The move was characteristic Plenio: take the geometric-distance idea proven for entanglement and apply it to the more elementary resource of basis superposition. From his chair at Ulm University he has driven the field into quantum thermodynamics, metrology, and the operational hierarchy of incoherent operations. The resource theory of coherence is, in a real sense, entanglement theory's younger sibling — and Plenio is a parent of both.

History
2006
— Åberg formalizes the idea of quantifying superposition relative to a fixed basis
2014
— Baumgratz, Cramer and Plenio fix the axioms and propose the relative-entropy and ℓ₁ measures
2016
— Winter and Yang prove distillable coherence equals the relative entropy of coherence
2016
— Chitambar, Gour and others map the inequivalent operation classes MIO, IO, DIO, SIO
2017
— Streltsov, Adesso and Plenio consolidate the field in a Reviews of Modern Physics colloquium

Going further

Texts and papers. Baumgratz, Cramer, Plenio, "Quantifying Coherence," Phys. Rev. Lett. 113, 140401 (2014) — the founding axioms · Streltsov, Adesso, Plenio, "Colloquium: Quantum coherence as a resource," Rev. Mod. Phys. 89, 041003 (2017) — the definitive review · Winter, Yang, "Operational Resource Theory of Coherence," Phys. Rev. Lett. 116, 120404 (2016) — distillation and cost · Åberg, "Quantifying Superposition," arXiv:quant-ph/0612146 (2006) — the precursor · Chitambar, Gour, "Quantum resource theories," Rev. Mod. Phys. 91, 025001 (2019) — the general framework and operation hierarchy · Nielsen & Chuang, Quantum Computation and Quantum Information (entropy and channels) · Wilde, Quantum Information Theory, 2nd ed. (relative entropy, data processing).

Provenance. The axioms C1–C2 and both measures are from Baumgratz–Cramer–Plenio (2014); the operation-class hierarchy (MIO, IO, DIO, SIO, PIO) is consolidated in Chitambar–Gour (2019) and the Streltsov–Adesso–Plenio colloquium (2017); the operational identity Cd=CrC_d=C_r is Winter–Yang (2016). The closed forms — Cr=S(Δρ)−S(ρ)C_r=S(\Delta\rho)-S(\rho), Cℓ1=∑i≠j∣ρij∣C_{\ell_1}=\sum_{i\neq j}|\rho_{ij}|, pure-state Cr=H({∣ci∣2})C_r=H(\{|c_i|^2\}), and the qudit unit values log⁡2d\log_2 d and d−1d-1 — are standard and independently verified above (Derivations A and C, worked examples). The monotonicity proof (Derivation B) rests only on IO preserving the incoherent set and the data-processing inequality for quantum relative entropy (Wilde, Ch. 11). A verifier should recompute the Gibbs minimization, the data-processing chain, and the qubit numbers S(ρ)=0.4690S(\rho)=0.4690 bits, Cr=0.531C_r=0.531 bits.

Prerequisites. Von Neumann Entropy · Quantum Relative Entropy · Quantum Channels and CPTP Maps · Kraus Representation · Data-Processing Inequality. Leads to. Holevo Bound · Strong Subadditivity of Quantum Entropy · Choi–Jamiołkowski Isomorphism.

Ask the tutor: Ask me to walk through the incoherent-Kraus structure behind the ℓ1\ell_1 monotone proof, or to show how coherence in one system converts into entanglement with an ancilla.

Why fix a basis at all — isn't singling out a preferred basis unphysical?

The basis is not arbitrary; it is chosen by the physics. Decoherence selects a pointer basis in which a system goes classical, thermodynamics selects the energy eigenbasis because Hamiltonians conserve energy, and quantum computing selects the computational basis in which gates and measurements are defined. Coherence is the resource relative to that physically distinguished basis — the superposition structure the environment, or the task, treats as special. Far from a flaw, basis-dependence is what makes the theory match reality: interference visibility, work extraction, and metrological sensitivity are all defined against a concrete basis, and coherence theory quantifies the resource each task consumes.

How is coherence different from entanglement?

They are cousins built on the same template but partitioned differently. Entanglement is the resource under local operations and classical communication (LOCC), with separable states free; coherence is the resource under incoherent operations, with basis-diagonal states free. Coherence is more elementary and needs only one system — a single qubit can be coherent but not entangled. There is a precise bridge: coherence can be converted into entanglement with an ancilla by incoherent operations, and the maximum entanglement generated equals a coherence measure. So entanglement is, in a sense, coherence distributed across a partition; the two theories share almost every structural theorem.

Why are there so many classes of "free" operations — MIO, IO, DIO, SIO, PIO?

Because "cannot create coherence" can be demanded at different strengths, and the demands are genuinely inequivalent. MIO (maximally incoherent operations) is weakest: any channel mapping the incoherent set to itself. IO strengthens this to require an incoherent Kraus decomposition, so no individual outcome creates coherence. DIO asks the channel to commute with dephasing; SIO asks each Kraus operator to also not use coherence; PIO restricts to what is physically implementable with incoherent ancillas. These are nested and strictly different — a measure can be monotone under one and not another, and reversibility holds under MIO but fails under IO. Unlike LOCC, there is no single obviously correct free class.

Do the two measures ever disagree about which state is more coherent?

Yes, and this is expected. In the worked examples the equatorial qubit had Cℓ1=0.8C_{\ell_1}=0.8 but Cr=0.531C_r=0.531 bits — different absolute values because they measure on different scales. More strikingly, for dimension three or higher CrC_r and Cℓ1C_{\ell_1} can order two states oppositely: A can beat B on CrC_r while B beats A on Cℓ1C_{\ell_1}. This is no contradiction; a resource can be genuinely multidimensional, as two portfolios rank differently by return and by risk. Only CrC_r carries the distillation meaning, so it is the measure of choice for "how much pure coherence can I extract," while Cℓ1C_{\ell_1} is the everyday diagnostic.

What does it mean operationally that C_r is the "distillable coherence"?

Distillation is the many-copy conversion of an imperfect resource into a perfect one. Given NN copies of ρ\rho, incoherent operations produce some number of flawless maximally coherent qubits ∣Ψ2⟩|\Psi_2\rangle; the best output-to-input ratio, as N→∞N\to\infty, is the distillable coherence, and Winter and Yang proved it equals Cr(ρ)C_r(\rho) exactly. So CrC_r is the physical exchange rate between your noisy state and the gold-standard coherbit: a state with Cr=0.531C_r=0.531 bits is worth about 531 pure coherent qubits per thousand copies. This is why the relative entropy of coherence, not the easier ℓ1\ell_1-norm, is regarded as the theory's fundamental measure.

Why does the maximally mixed state have zero coherence despite being maximally random?

Because coherence measures off-diagonal phase order, not statistical spread. The maximally mixed state I/dI/d is proportional to the identity, so it is already diagonal in every basis — including the reference basis — which makes it free: Δ(I/d)=I/d\Delta(I/d)=I/d gives Cr=0C_r=0, and all off-diagonals vanish so Cℓ1=0C_{\ell_1}=0. Its large von Neumann entropy is classical uncertainty about which basis state it is in, which a classical apparatus reproduces for free. Coherence, by contrast, is the quantum phase relation between basis states, which no classical process can create — and I/dI/d has none of it.

Before this

Where this leads