Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

QCQI – Chapter 5 Detailed Summary: Density Operators & Quantum Operations

Learning Objectives

  • Work with density operators: mixtures, partial trace, Bloch form.

  • Define and use quantum operations (CPTP maps) with Kraus operators.

  • Interpret channels via Stinespring dilation and Choi–Jamiołkowski isomorphism.

  • Compute trace distance and fidelity; understand data-processing monotonicity.

  • Model common noise (depolarizing, dephasing, amplitude damping) and its Bloch-affine action.

Density Operators

A state is a positive semidefinite operator with unit trace: ρ0\rho\ge 0, Trρ=1\operatorname{Tr}\rho=1. For a qubit,

ρ=12 ⁣(I+rxX+ryY+rzZ),rR3, r1.\begin{aligned} \rho = \frac{1}{2}\!\left(I + r_x X + r_y Y + r_z Z\right),\qquad \vec r\in\mathbb{R}^3,\ \|\vec r\|\le 1. \end{aligned}

Purity: Tr(ρ2)=1+r22\operatorname{Tr}(\rho^2)=\tfrac{1+\|\vec r\|^2}{2}.

Eureka! Pure r=1\Leftrightarrow \|\vec r\|=1. Noise shrinks r\|\vec r\| inward.

Quantum Operations (CPTP)

A physical operation E\mathcal{E} is linear, completely positive, and trace-preserving (CPTP). Kraus form:

E(ρ)=iKiρKi,iKiKi=I.\begin{aligned} \mathcal{E}(\rho)=\sum_i K_i \rho K_i^\dagger,\qquad \sum_i K_i^\dagger K_i = I. \end{aligned}

Kraus representations are not unique (related by isometries).

Eureka! “Many realizations, one map”: different {Ki}\{K_i\} sets describe the same channel.

Dilation Pictures

Stinespring. There exists an ancilla EE, unitary UU on S ⁣ES\!E, and fixed 0E|0\rangle_E such that

E(ρ)=TrE[U(ρ0 ⁣0)U].\begin{aligned} \mathcal{E}(\rho)=\operatorname{Tr}_E\left[U\,(\rho\otimes |0\rangle\!\langle 0|)\,U^\dagger\right]. \end{aligned}

Naimark. Any POVM arises from a projective measurement on a larger Hilbert space.

Eureka! Open dynamics = unitary on a bigger space + discard. Measurements = projective on a bigger space.

Choi–Jamiołkowski (CJ) Isomorphism

Define the (unnormalized) maximally entangled state Φ=jjj|\Phi\rangle=\sum_j |j\rangle|j\rangle and the Choi matrix

J(E)=(EI)(Φ ⁣Φ).\begin{aligned} J(\mathcal{E}) = (\mathcal{E}\otimes I)\big(|\Phi\rangle\!\langle \Phi|\big). \end{aligned}

E\mathcal{E} is CP iff J(E)0J(\mathcal{E})\ge 0 and TP iff TroutJ(E)=Iin\operatorname{Tr}_{\text{out}} J(\mathcal{E})=I_{\text{in}}.

Eureka! Checking CP reduces to checking positivity of a single matrix: the Choi matrix.

Distances & Fidelity

Trace distance D(ρ,σ)=12ρσ1D(\rho,\sigma)=\tfrac12\|\rho-\sigma\|_1 contracts under CPTP:

D(E(ρ),E(σ))D(ρ,σ).\begin{aligned} D(\mathcal{E}(\rho),\mathcal{E}(\sigma)) \le D(\rho,\sigma). \end{aligned}

Uhlmann fidelity F(ρ,σ)=(Trρσρ)2F(\rho,\sigma)=\left(\operatorname{Tr}\sqrt{\sqrt{\rho}\sigma\sqrt{\rho}}\right)^2 is non-decreasing under CPTP:

F(E(ρ),E(σ))F(ρ,σ).\begin{aligned} F(\mathcal{E}(\rho),\mathcal{E}(\sigma)) \ge F(\rho,\sigma). \end{aligned}

Eureka! Processing information can only make states harder to tell apart (distance down, fidelity up).

Canonical Qubit Channels (Bloch-Affine Forms)

  • Depolarizing (prob. pp): ρ(1p)ρ+p2I\rho \mapsto (1-p)\rho + \frac{p}{2}I; r(1p)r\vec r\mapsto (1-p)\vec r.

  • Dephasing / Phase-damping (pp): rx,ry(12p)rx,yr_x,r_y\mapsto (1-2p)\,r_{x,y}; rzr_z unchanged.

  • Amplitude damping (γ\gamma): non-unital; roughly rx,y1γrx,yr_{x,y}\mapsto \sqrt{1-\gamma}\,r_{x,y}, rz(1γ)rz+γr_z\mapsto (1-\gamma)r_z+\gamma.

Eureka! Unital channels map I ⁣ ⁣I/2I ⁣ ⁣I/2I\!\!I/2 \mapsto I\!\!I/2 (Bloch shift b=0\vec b=0). Amplitude damping is non-unital (b0\vec b\neq 0).

POVMs

POVM elements {Em}\{E_m\} satisfy Em0E_m\ge 0, mEm=I\sum_m E_m=I. Probabilities: p(m)=Tr(Emρ)p(m)=\operatorname{Tr}(E_m\rho).

Realization via Naimark: embed in larger space, apply unitary, measure projectively, and post-process.

Worked Examples

Amplitude damping Kraus set. K0=(1001γ)K_0=\begin{pmatrix}1&0\\0&\sqrt{1-\gamma}\end{pmatrix}, K1=(0γ00)K_1=\begin{pmatrix}0&\sqrt{\gamma}\\0&0\end{pmatrix}.

Choi test. Build J(E)J(\mathcal{E}); check PSD and TP by partial trace.

Monotonicity demo. Pick ρ,σ\rho,\sigma; compute D,FD,F before/after E\mathcal{E}.

Pseudo-code Aids

Pseudocode — Apply a channel in Kraus form

Input: density matrix ρ\rho, list of Kraus ops {Ki}\{K_i\} Check: iKiKiI\sum_i K_i^\dagger K_i \approx I (TP) Return: iKiρKi\sum_i K_i \rho K_i^\dagger

Pseudocode — Build the Choi matrix

Input: channel E\mathcal{E} via {Ki}\{K_i\} on dd-dim space Construct Φ=j=0d1jj|\Phi\rangle=\sum_{j=0}^{d-1}|j\rangle|j\rangle and Φ=Φ ⁣Φ\Phi=|\Phi\rangle\!\langle \Phi| Return J(E)=i(KiI)Φ(KiI)J(\mathcal{E})=\sum_i (K_i\otimes I)\,\Phi\,(K_i^\dagger\otimes I)

Pseudocode — Stinespring dilation (constructive sketch)

Input: Kraus ops {Ki}\{K_i\} Form isometry V=iKiiEV=\sum_i K_i\otimes |i\rangle_E; extend VV to a unitary UU on S ⁣ES\!E Realize E(ρ)=TrE[U(ρ0 ⁣0)U]\mathcal{E}(\rho)=\operatorname{Tr}_E\big[U(\rho\otimes|0\rangle\!\langle 0|)U^\dagger\big]

Schematics

Stinespring Dilation

Choi Test

Hands-on Notebook (Multi-Backend)

Run: L06_QCQI_Ch05_Density_Operators_and_Quantum_Operations

In the first cell set backend to one of: cirq, pennylane, braket, pyquil, or qiskit.

Mini-Lab Ideas

  1. Verify TP/unital conditions for given Kraus sets; compute Bloch-affine (A,b)(A,\vec b).

  2. Build J(E)J(\mathcal{E}) for amplitude damping; check PSD and TP by partial trace.

  3. Numerically confirm DD decreases and FF increases under (i) dephasing, (ii) amplitude damping.