These exercises practice the partial trace as a computational rule and as a physical rule for local statistics. The main identity is
Tr B ( ∣ i ⟩ A ⟨ k ∣ A ⊗ ∣ j ⟩ B ⟨ l ∣ B ) = δ j l ∣ i ⟩ A ⟨ k ∣ A . \operatorname{Tr}_B
\bigl(
\lvert i\rangle_A\langle k\rvert_A
\otimes
\lvert j\rangle_B\langle l\rvert_B
\bigr)
=
\delta_{jl}
\lvert i\rangle_A\langle k\rvert_A. Tr B ( ∣ i ⟩ A ⟨ k ∣ A ⊗ ∣ j ⟩ B ⟨ l ∣ B ) = δ j l ∣ i ⟩ A ⟨ k ∣ A .
Equivalently, the traced subsystem contributes an inner product between the ket and bra labels that remain on that subsystem.
Let
∣ Ψ ⟩ A B = ∣ ψ ⟩ A ⊗ ∣ ϕ ⟩ B . \lvert\Psi\rangle_{AB}
=
\lvert\psi\rangle_A\otimes\lvert\phi\rangle_B. ∣ Ψ ⟩ A B = ∣ ψ ⟩ A ⊗ ∣ ϕ ⟩ B .
Compute ρ A = Tr B ∣ Ψ ⟩ ⟨ Ψ ∣ \rho_A=\operatorname{Tr}_B\lvert\Psi\rangle\langle\Psi\rvert ρ A = Tr B ∣ Ψ ⟩ ⟨ Ψ ∣ .
Solution
The joint density operator is
ρ A B = ∣ ψ ⟩ ⟨ ψ ∣ A ⊗ ∣ ϕ ⟩ ⟨ ϕ ∣ B . \rho_{AB}
=
\lvert\psi\rangle\langle\psi\rvert_A
\otimes
\lvert\phi\rangle\langle\phi\rvert_B. ρ A B = ∣ ψ ⟩ ⟨ ψ ∣ A ⊗ ∣ ϕ ⟩ ⟨ ϕ ∣ B .
Taking the trace over B B B gives
ρ A = ∣ ψ ⟩ ⟨ ψ ∣ A Tr ( ∣ ϕ ⟩ ⟨ ϕ ∣ B ) . \rho_A
=
\lvert\psi\rangle\langle\psi\rvert_A
\operatorname{Tr}
\bigl(
\lvert\phi\rangle\langle\phi\rvert_B
\bigr). ρ A = ∣ ψ ⟩ ⟨ ψ ∣ A Tr ( ∣ ϕ ⟩ ⟨ ϕ ∣ B ) .
Since ∣ ϕ ⟩ \lvert\phi\rangle ∣ ϕ ⟩ is normalized,
Tr ( ∣ ϕ ⟩ ⟨ ϕ ∣ ) = 1. \operatorname{Tr}
\bigl(
\lvert\phi\rangle\langle\phi\rvert
\bigr)
=
1. Tr ( ∣ ϕ ⟩ ⟨ ϕ ∣ ) = 1.
Therefore
ρ A = ∣ ψ ⟩ ⟨ ψ ∣ A . \rho_A
=
\lvert\psi\rangle\langle\psi\rvert_A. ρ A = ∣ ψ ⟩ ⟨ ψ ∣ A .
A product pure state has pure reduced states.
Compute the one-qubit reduced state of
∣ Φ + ⟩ = 1 2 ( ∣ 00 ⟩ + ∣ 11 ⟩ ) . \lvert\Phi^+\rangle
=
\frac{1}{\sqrt2}
\bigl(
\lvert00\rangle+\lvert11\rangle
\bigr). ∣ Φ + ⟩ = 2 1 ( ∣ 00 ⟩ + ∣ 11 ⟩ ) .
Solution
The density operator is
ρ A B = 1 2 ( ∣ 00 ⟩ ⟨ 00 ∣ + ∣ 00 ⟩ ⟨ 11 ∣ + ∣ 11 ⟩ ⟨ 00 ∣ + ∣ 11 ⟩ ⟨ 11 ∣ ) . \rho_{AB}
=
\frac12
\bigl(
\lvert00\rangle\langle00\rvert
+\lvert00\rangle\langle11\rvert
+\lvert11\rangle\langle00\rvert
+\lvert11\rangle\langle11\rvert
\bigr). ρ A B = 2 1 ( ∣ 00 ⟩ ⟨ 00 ∣ + ∣ 00 ⟩ ⟨ 11 ∣ + ∣ 11 ⟩ ⟨ 00 ∣ + ∣ 11 ⟩ ⟨ 11 ∣ ) .
Trace over B B B . The two off-diagonal terms vanish because ⟨ 1 ∣ 0 ⟩ = 0 \langle1\vert0\rangle=0 ⟨ 1∣0 ⟩ = 0 and ⟨ 0 ∣ 1 ⟩ = 0 \langle0\vert1\rangle=0 ⟨ 0∣1 ⟩ = 0 :
ρ A = 1 2 ( ∣ 0 ⟩ ⟨ 0 ∣ + ∣ 1 ⟩ ⟨ 1 ∣ ) = 1 2 I . \rho_A
=
\frac12
\bigl(
\lvert0\rangle\langle0\rvert
+\lvert1\rangle\langle1\rvert
\bigr)
=
\frac12 I. ρ A = 2 1 ( ∣ 0 ⟩ ⟨ 0 ∣ + ∣ 1 ⟩ ⟨ 1 ∣ ) = 2 1 I .
The same calculation gives ρ B = I / 2 \rho_B=I/2 ρ B = I /2 .
The classically correlated state
ρ A B = 1 2 ∣ 00 ⟩ ⟨ 00 ∣ + 1 2 ∣ 11 ⟩ ⟨ 11 ∣ \rho_{AB}
=
\frac12\lvert00\rangle\langle00\rvert
+\frac12\lvert11\rangle\langle11\rvert ρ A B = 2 1 ∣ 00 ⟩ ⟨ 00 ∣ + 2 1 ∣ 11 ⟩ ⟨ 11 ∣
has the same one-qubit reduced states as ∣ Φ + ⟩ \lvert\Phi^+\rangle ∣ Φ + ⟩ . Verify this, and state one global difference between the two states.
Solution
Trace over B B B term by term:
Tr B ( ∣ 00 ⟩ ⟨ 00 ∣ ) = ∣ 0 ⟩ ⟨ 0 ∣ , \operatorname{Tr}_B
\bigl(
\lvert00\rangle\langle00\rvert
\bigr)
=
\lvert0\rangle\langle0\rvert, Tr B ( ∣ 00 ⟩ ⟨ 00 ∣ ) = ∣ 0 ⟩ ⟨ 0 ∣ ,
and
Tr B ( ∣ 11 ⟩ ⟨ 11 ∣ ) = ∣ 1 ⟩ ⟨ 1 ∣ . \operatorname{Tr}_B
\bigl(
\lvert11\rangle\langle11\rvert
\bigr)
=
\lvert1\rangle\langle1\rvert. Tr B ( ∣ 11 ⟩ ⟨ 11 ∣ ) = ∣ 1 ⟩ ⟨ 1 ∣ .
Therefore
ρ A = 1 2 I , ρ B = 1 2 I . \rho_A
=
\frac12 I,
\qquad
\rho_B
=
\frac12 I. ρ A = 2 1 I , ρ B = 2 1 I .
The global states are nevertheless different. The Bell state is a pure entangled state with coherence between ∣ 00 ⟩ \lvert00\rangle ∣ 00 ⟩ and ∣ 11 ⟩ \lvert11\rangle ∣ 11 ⟩ . The classically correlated state is mixed and has no such off-diagonal coherence.
Let
∣ G H Z 3 ⟩ = 1 2 ( ∣ 000 ⟩ + ∣ 111 ⟩ ) . \lvert\mathrm{GHZ}_3\rangle
=
\frac{1}{\sqrt2}
\bigl(
\lvert000\rangle+\lvert111\rangle
\bigr). ∣ GHZ 3 ⟩ = 2 1 ( ∣ 000 ⟩ + ∣ 111 ⟩ ) .
Trace out qubit C C C and compute ρ A B \rho_{AB} ρ A B .
Solution
The density operator contains four terms:
ρ A B C = 1 2 ( ∣ 000 ⟩ ⟨ 000 ∣ + ∣ 000 ⟩ ⟨ 111 ∣ + ∣ 111 ⟩ ⟨ 000 ∣ + ∣ 111 ⟩ ⟨ 111 ∣ ) . \rho_{ABC}
=
\frac12
\bigl(
\lvert000\rangle\langle000\rvert
+\lvert000\rangle\langle111\rvert
+\lvert111\rangle\langle000\rvert
+\lvert111\rangle\langle111\rvert
\bigr). ρ A B C = 2 1 ( ∣ 000 ⟩ ⟨ 000 ∣ + ∣ 000 ⟩ ⟨ 111 ∣ + ∣ 111 ⟩ ⟨ 000 ∣ + ∣ 111 ⟩ ⟨ 111 ∣ ) .
Tracing over C C C kills the cross terms because the C C C labels are orthogonal:
⟨ 1 ∣ 0 ⟩ = 0 , ⟨ 0 ∣ 1 ⟩ = 0. \langle1\vert0\rangle=0,
\qquad
\langle0\vert1\rangle=0. ⟨ 1∣0 ⟩ = 0 , ⟨ 0∣1 ⟩ = 0.
The remaining terms are
ρ A B = 1 2 ∣ 00 ⟩ ⟨ 00 ∣ + 1 2 ∣ 11 ⟩ ⟨ 11 ∣ . \rho_{AB}
=
\frac12
\lvert00\rangle\langle00\rvert
+\frac12
\lvert11\rangle\langle11\rvert. ρ A B = 2 1 ∣ 00 ⟩ ⟨ 00 ∣ + 2 1 ∣ 11 ⟩ ⟨ 11 ∣ .
Thus the two-qubit reduction is a separable classical mixture, not a Bell state.
Let
∣ W 3 ⟩ = 1 3 ( ∣ 100 ⟩ + ∣ 010 ⟩ + ∣ 001 ⟩ ) . \lvert W_3\rangle
=
\frac{1}{\sqrt3}
\bigl(
\lvert100\rangle+\lvert010\rangle+\lvert001\rangle
\bigr). ∣ W 3 ⟩ = 3 1 ( ∣ 100 ⟩ + ∣ 010 ⟩ + ∣ 001 ⟩ ) .
Trace out qubit C C C and compute ρ A B \rho_{AB} ρ A B .
Solution
Group the state by the value of qubit C C C :
∣ W 3 ⟩ = 1 3 ( ∣ 10 ⟩ A B + ∣ 01 ⟩ A B ) ∣ 0 ⟩ C + 1 3 ∣ 00 ⟩ A B ∣ 1 ⟩ C . \lvert W_3\rangle
=
\frac{1}{\sqrt3}
\bigl(
\lvert10\rangle_{AB}
+\lvert01\rangle_{AB}
\bigr)
\lvert0\rangle_C
+\frac{1}{\sqrt3}
\lvert00\rangle_{AB}\lvert1\rangle_C. ∣ W 3 ⟩ = 3 1 ( ∣ 10 ⟩ A B + ∣ 01 ⟩ A B ) ∣ 0 ⟩ C + 3 1 ∣ 00 ⟩ A B ∣ 1 ⟩ C .
The two C C C states are orthogonal, so the trace over C C C removes cross terms between the C = 0 C=0 C = 0 and C = 1 C=1 C = 1 branches. The C = 0 C=0 C = 0 branch has norm squared 2 / 3 2/3 2/3 and normalized A B AB A B state
∣ Ψ + ⟩ = 1 2 ( ∣ 01 ⟩ + ∣ 10 ⟩ ) . \lvert\Psi^+\rangle
=
\frac{1}{\sqrt2}
\bigl(
\lvert01\rangle+\lvert10\rangle
\bigr). ∣ Ψ + ⟩ = 2 1 ( ∣ 01 ⟩ + ∣ 10 ⟩ ) .
Therefore
ρ A B = 2 3 ∣ Ψ + ⟩ ⟨ Ψ + ∣ + 1 3 ∣ 00 ⟩ ⟨ 00 ∣ . \rho_{AB}
=
\frac23
\lvert\Psi^+\rangle\langle\Psi^+\rvert
+\frac13
\lvert00\rangle\langle00\rvert. ρ A B = 3 2 ∣ Ψ + ⟩ ⟨ Ψ + ∣ + 3 1 ∣ 00 ⟩ ⟨ 00 ∣ .
This reduction remains partly coherent in the one-excitation subspace.
For ∣ Φ + ⟩ \lvert\Phi^+\rangle ∣ Φ + ⟩ , use the reduced state to compute the probabilities for measuring Z Z Z on subsystem A A A . Then compare with the joint measurement of Z ⊗ Z Z\otimes Z Z ⊗ Z .
Solution
From Exercise 2,
ρ A = 1 2 I . \rho_A=\frac12 I. ρ A = 2 1 I .
The projectors for measuring Z Z Z are
P 0 = ∣ 0 ⟩ ⟨ 0 ∣ , P 1 = ∣ 1 ⟩ ⟨ 1 ∣ . P_0=\lvert0\rangle\langle0\rvert,
\qquad
P_1=\lvert1\rangle\langle1\rvert. P 0 = ∣ 0 ⟩ ⟨ 0 ∣ , P 1 = ∣ 1 ⟩ ⟨ 1 ∣ .
Thus
p ( 0 ) = Tr ( ρ A P 0 ) = 1 2 , p ( 1 ) = Tr ( ρ A P 1 ) = 1 2 . p(0)
=
\operatorname{Tr}(\rho_A P_0)
=
\frac12,
\qquad
p(1)
=
\operatorname{Tr}(\rho_A P_1)
=
\frac12. p ( 0 ) = Tr ( ρ A P 0 ) = 2 1 , p ( 1 ) = Tr ( ρ A P 1 ) = 2 1 .
The local outcome is random. But the joint observable has perfect correlation:
( Z ⊗ Z ) ∣ Φ + ⟩ = ∣ Φ + ⟩ . (Z\otimes Z)\lvert\Phi^+\rangle
=
\lvert\Phi^+\rangle. ( Z ⊗ Z ) ∣ Φ + ⟩ = ∣ Φ + ⟩ .
Therefore ⟨ Z ⊗ Z ⟩ = 1 \langle Z\otimes Z\rangle=1 ⟨ Z ⊗ Z ⟩ = 1 . Local reduced states give local statistics, not all joint correlations.
Construct a purification of
ρ A = p ∣ 0 ⟩ ⟨ 0 ∣ + ( 1 − p ) ∣ 1 ⟩ ⟨ 1 ∣ , 0 ≤ p ≤ 1 , \rho_A
=
p\lvert0\rangle\langle0\rvert
+(1-p)\lvert1\rangle\langle1\rvert,
\qquad
0\le p\le1, ρ A = p ∣ 0 ⟩ ⟨ 0 ∣ + ( 1 − p ) ∣ 1 ⟩ ⟨ 1 ∣ , 0 ≤ p ≤ 1 ,
and verify that tracing out the reference system gives ρ A \rho_A ρ A .
Solution
Introduce a reference qubit R R R and define
∣ Ψ ⟩ A R = p ∣ 0 ⟩ A ∣ 0 ⟩ R + 1 − p ∣ 1 ⟩ A ∣ 1 ⟩ R . \lvert\Psi\rangle_{AR}
=
\sqrt p\,\lvert0\rangle_A\lvert0\rangle_R
+\sqrt{1-p}\,
\lvert1\rangle_A\lvert1\rangle_R. ∣ Ψ ⟩ A R = p ∣ 0 ⟩ A ∣ 0 ⟩ R + 1 − p ∣ 1 ⟩ A ∣ 1 ⟩ R .
The projector is
∣ Ψ ⟩ ⟨ Ψ ∣ = p ∣ 00 ⟩ ⟨ 00 ∣ + p ( 1 − p ) ∣ 00 ⟩ ⟨ 11 ∣ + p ( 1 − p ) ∣ 11 ⟩ ⟨ 00 ∣ + ( 1 − p ) ∣ 11 ⟩ ⟨ 11 ∣ . \begin{aligned}
\lvert\Psi\rangle\langle\Psi\rvert
&=
p\lvert00\rangle\langle00\rvert
+\sqrt{p(1-p)}\,
\lvert00\rangle\langle11\rvert\\
&\quad
+\sqrt{p(1-p)}\,
\lvert11\rangle\langle00\rvert
+(1-p)\lvert11\rangle\langle11\rvert.
\end{aligned} ∣ Ψ ⟩ ⟨ Ψ ∣ = p ∣ 00 ⟩ ⟨ 00 ∣ + p ( 1 − p ) ∣ 00 ⟩ ⟨ 11 ∣ + p ( 1 − p ) ∣ 11 ⟩ ⟨ 00 ∣ + ( 1 − p ) ∣ 11 ⟩ ⟨ 11 ∣ .
Tracing over R R R removes the cross terms because the reference states are orthogonal. The result is
Tr R ( ∣ Ψ ⟩ ⟨ Ψ ∣ ) = p ∣ 0 ⟩ ⟨ 0 ∣ + ( 1 − p ) ∣ 1 ⟩ ⟨ 1 ∣ = ρ A . \operatorname{Tr}_R
\bigl(
\lvert\Psi\rangle\langle\Psi\rvert
\bigr)
=
p\lvert0\rangle\langle0\rvert
+(1-p)\lvert1\rangle\langle1\rvert
=
\rho_A. Tr R ( ∣ Ψ ⟩ ⟨ Ψ ∣ ) = p ∣ 0 ⟩ ⟨ 0 ∣ + ( 1 − p ) ∣ 1 ⟩ ⟨ 1 ∣ = ρ A .
Let
ρ A B = 1 4 ( I ⊗ I + r Z ⊗ I + s I ⊗ Z + t Z ⊗ Z ) . \rho_{AB}
=
\frac14
\bigl(
I\otimes I
+r\,Z\otimes I
+s\,I\otimes Z
+t\,Z\otimes Z
\bigr). ρ A B = 4 1 ( I ⊗ I + r Z ⊗ I + s I ⊗ Z + t Z ⊗ Z ) .
Compute ρ A \rho_A ρ A and identify which coefficient is invisible to measurements on A A A alone.
Solution
Use
Tr ( I ) = 2 , Tr ( Z ) = 0. \operatorname{Tr}(I)=2,
\qquad
\operatorname{Tr}(Z)=0. Tr ( I ) = 2 , Tr ( Z ) = 0.
Then
ρ A = Tr B ρ A B = 1 4 ( 2 I + 2 r Z + s I Tr Z + t Z Tr Z ) = 1 2 ( I + r Z ) . \begin{aligned}
\rho_A
&=
\operatorname{Tr}_B\rho_{AB}\\
&=
\frac14
\bigl(
2I
+2rZ
+sI\,\operatorname{Tr}Z
+tZ\,\operatorname{Tr}Z
\bigr)\\
&=
\frac12
\bigl(
I+rZ
\bigr).
\end{aligned} ρ A = Tr B ρ A B = 4 1 ( 2 I + 2 r Z + s I Tr Z + tZ Tr Z ) = 2 1 ( I + r Z ) .
The coefficients s s s and t t t do not appear in ρ A \rho_A ρ A . In particular, the correlation coefficient t t t is invisible to measurements on A A A alone even though it affects joint statistics.
R. Shankar, Principles of Quantum Mechanics , 2nd ed., Springer, 1994.
L. E. Ballentine, Quantum Mechanics: A Modern Development , 2nd ed., World Scientific, 2014.
M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information , Cambridge University Press, 2010.
H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems , Oxford University Press, 2002.