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),r⃗∈R3, ∥r⃗∥≤1.\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+∥r⃗∥22\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†,∑iKi†Ki=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 ∣0⟩E|0\rangle_E such that

E(ρ)=Tr⁡E[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 ∣Φ⟩=∑j∣j⟩∣j⟩|\Phi\rangle=\sum_j |j\rangle|j\rangle and the Choi matrix

J(E)=(E⊗I)(∣Φ⟩ ⁣⟨Φ∣).\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 Tr⁡outJ(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): ρ↦(1−p)ρ+p2I\rho \mapsto (1-p)\rho + \frac{p}{2}I; r⃗↦(1−p)r⃗\vec r\mapsto (1-p)\vec r.

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

  • Amplitude damping (γ\gamma): non-unital; roughly rx,y↦1−γ 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/2↦I ⁣ ⁣I/2I\!\!I/2 \mapsto I\!\!I/2 (Bloch shift b⃗=0\vec b=0). Amplitude damping is non-unital (b⃗≠0\vec b\neq 0).

POVMs

POVM elements {Em}\{E_m\} satisfy Em≥0E_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: ∑iKi†Ki≈I\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=0d−1∣j⟩∣j⟩|\Phi\rangle=\sum_{j=0}^{d-1}|j\rangle|j\rangle and Φ=∣Φ⟩ ⁣⟨Φ∣\Phi=|\Phi\rangle\!\langle \Phi| Return J(E)=∑i(Ki⊗I) Φ (Ki†⊗I)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=∑iKi⊗∣i⟩EV=\sum_i K_i\otimes |i\rangle_E; extend VV to a unitary UU on S ⁣ES\!E Realize E(ρ)=Tr⁡E[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.