|q⟩ Bad Qubits

advanced · Physics · Quantum Channels & Open Systems

Checkpoint: Analyze a Channel

This checkpoint ties the module together. You will take a concrete channel, apply it to a coherent superposition through its Kraus operators, and read off a quantitative property of the output — the kind of end-to-end calculation that underlies every noise analysis.

What to recall

The problem

Prepare the coherent superposition +=12(0+1)|+\rangle = \tfrac{1}{\sqrt2}(|0\rangle+|1\rangle), so ρ=++=12(1111)\rho = |+\rangle\langle+| = \tfrac12\begin{pmatrix}1&1\\1&1\end{pmatrix}. Send it through the amplitude-damping channel with γ=12\gamma = \tfrac12. Then compute the purity of the output state, tr(E(ρ)2)\operatorname{tr}\big(\mathcal{E}(\rho)^2\big).

Working it out by hand

Applying the channel to ρ\rho:

E(ρ)=(ρ00+γρ111γρ011γρ10(1γ)ρ11)=(12+1212121212121212)=(3412212214).\mathcal{E}(\rho) = \begin{pmatrix} \rho_{00}+\gamma\rho_{11} & \sqrt{1-\gamma}\,\rho_{01}\\[4pt] \sqrt{1-\gamma}\,\rho_{10} & (1-\gamma)\rho_{11} \end{pmatrix} = \begin{pmatrix} \tfrac12+\tfrac12\cdot\tfrac12 & \tfrac{1}{\sqrt2}\cdot\tfrac12\\[4pt] \tfrac{1}{\sqrt2}\cdot\tfrac12 & \tfrac12\cdot\tfrac12 \end{pmatrix} = \begin{pmatrix} \tfrac34 & \tfrac{1}{2\sqrt2}\\[4pt] \tfrac{1}{2\sqrt2} & \tfrac14 \end{pmatrix}.

The purity is the sum of the squared magnitudes of all entries (since the matrix is Hermitian, tr(ρ2)=ijρij2\operatorname{tr}(\rho^2)=\sum_{ij}|\rho_{ij}|^2):

tr(E(ρ)2)=(34)2+(14)2+2(122)2=916+116+218=1016+416=1416=0.875.\operatorname{tr}\big(\mathcal{E}(\rho)^2\big) = \Big(\tfrac34\Big)^2 + \Big(\tfrac14\Big)^2 + 2\Big(\tfrac{1}{2\sqrt2}\Big)^2 = \tfrac{9}{16} + \tfrac{1}{16} + 2\cdot\tfrac18 = \tfrac{10}{16}+\tfrac{4}{16} = \tfrac{14}{16} = 0.875 .

The pure input (purity 1\text{purity }1) has been driven to a partially mixed output — damping has injected entropy. Confirm 0.8750.875 by building E(ρ)\mathcal{E}(\rho) from the Kraus operators and squaring.

Your task

Build E(ρ)\mathcal{E}(\rho) for ρ=++\rho=|+\rangle\langle+| and γ=12\gamma=\tfrac12 using E(ρ)=K0ρK0+K1ρK1\mathcal{E}(\rho)=K_0\rho K_0^\dagger + K_1\rho K_1^\dagger, then return tr(E(ρ)2)\operatorname{tr}\big(\mathcal{E}(\rho)^2\big).

Run your code to see the quantum state.

Sign in on the full site to ask questions and join the discussion.