advanced · Physics · Quantum Information & von Neumann Entropy
Von Neumann Entropy
The classical Shannon entropy measures the uncertainty in a probability distribution. A quantum
state is described not by a distribution but by a density operatorρ — a positive
semidefinite operator with Trρ=1. John von Neumann's entropy is the direct
generalization, obtained by replacing the probabilities with the operator ρ itself.
Definition
The von Neumann entropy of a density operator ρ is
S(ρ)=−Tr(ρlog2ρ).
Here log2ρ is the matrix logarithm, defined through the spectral decomposition. Writing
ρ in its eigenbasis,
ρ=k∑λk∣k⟩⟨k∣,λk≥0,k∑λk=1,
the operator ρlog2ρ is diagonal in the same basis with entries λklog2λk,
and the trace collects them:
S(ρ)=−k∑λklog2λk=H({λk}).
So the von Neumann entropy is precisely the Shannon entropy of the eigenvalues of ρ, again
with 0log20=0. The eigenvalues form a genuine probability distribution because ρ is a
positive operator of unit trace.
Why the eigenvalues, and nothing else
The entropy is basis-independent: it depends only on the spectrum of ρ, not on how ρ
is written. This is forced by the trace, which is invariant under unitary conjugation. If
ρ′=UρU† for any unitary U, then ρ and ρ′ have the same eigenvalues and
hence S(ρ′)=S(ρ). Physically: applying a reversible (unitary) transformation never changes
the information content of a state.
Basic properties
Non-negativity.S(ρ)≥0, with S(ρ)=0 iff ρ is pure (a single
eigenvalue 1, the rest 0).
Maximum. On a d-dimensional system, S(ρ)≤log2d, with equality iff
ρ=1/d is maximally mixed.
Unitary invariance.S(UρU†)=S(ρ), as shown above.
Additivity on products.S(ρ⊗σ)=S(ρ)+S(σ).
A worked qubit example
Take the density matrix
ρ=(21414121).
It is Hermitian, positive, and has unit trace. Its eigenvalues solve
λ2−λ+detρ=0 with detρ=41−161=163, giving
λ±=21±1−43=21±21={43,41}.
Therefore
S(ρ)=−43log243−41log241=h(41)≈0.8113 bits.
Note that S(ρ)<1: although the off-diagonal coherences make ρ look "spread out", its
eigenvalues are unequal, so it is less mixed than the maximally mixed qubit 1/2 (which has
S=1).
Try it
Compute the von Neumann entropy in bits of
ρ=(0.50.250.250.5)
by diagonalizing it and applying S=−∑kλklog2λk. Return the value.
Run your code to see the quantum state.
Sign in on the full site to ask questions and join the discussion.