16.7.1 Density matrix for a subsystem
LetHab=Ha⊗Hbbe the Hilbert space organized as the tensor product of two subspacesHa
andHb. Letρabbe the density matrix corresponding to a state of the system. Suppose now that
we observe the system with observables that act only on the subsystema. In all generality, such
observables are of the form,
A=Aa⊗Ib (16.48)The expectation value ofAof this type of operators in the full system reduces in a systematic
manner to the expectation value with respect to the density matrix of the subsystem. To see this,
we compute the expectation values,
〈A〉 = TrHab(
ρabA)
= TrHaTrHb(
ρab(Aa⊗Ib))= TrHa(
Aaρa)
(16.49)where
ρa= TrHb(
ρab)
(16.50)Note that the normalization TrHa(ρa) = 1 as a result of the normalization ofρab. In particular, the
population fraction of a pure state|ψa〉〈ψa|is given by
pa= TrHa(
ρa|ψa〉〈ψa|)
=〈ψa|ρa|ψa〉 (16.51)16.7.2 Example of relations between density matrices of subsystems
It will be helpful to clarify the relations betweenρab,ρa,ρband the tensor productρa⊗ρbwith the
help of an example involving two spin 1/2 systems. We denote the Pauli matrices for these systems
respectively byσiaandσibwherei= 1, 2 ,3. Consider general density matrices for each subsystem,
ρa =1
2
(Ia+~a·~σa) |~a|≤ 1ρb =1
2
(
Ib+~b·~σb)
|~b|≤ 1 (16.52)The tensor product is given by
ρa⊗ρb=1
4
(
Ia⊗Ib+Ia⊗(~b·~σb) + (~a·~σa)⊗Ib+ (~a·~σa)⊗(~b·~σb))
(16.53)Clearly, when we compute the partial traces, we recover the density matrices of the subsystems,
trHa(ρa⊗ρb) = ρb
trHb(ρa⊗ρb) = ρa (16.54)