Density Operator

A positive semidefinite trace-one operator representing the state of a quantum system, allowing both pure and statistical mixtures.
Density Operator

Let HH be a complex Hilbert space (typically finite-dimensional in basic quantum theory). A density operator (also called a density matrix) is an operator ρ:HH\rho:H\to H such that:

  1. Positivity: ρ0\rho \ge 0, meaning ψ,ρψ0\langle \psi,\rho\psi\rangle \ge 0 for all ψH\psi\in H.
  2. Unit trace: Tr(ρ)=1\operatorname{Tr}(\rho)=1, where Tr\operatorname{Tr} is the operator trace ( ).

In finite dimension, these conditions are equivalent to ρ\rho being a positive semidefinite matrix with trace 11.

Basic structural facts (finite dimension)

  • ρ\rho is automatically self-adjoint: ρ=ρ\rho=\rho^\ast.
  • All eigenvalues of ρ\rho are real and lie in [0,1][0,1].
  • The eigenvalues sum to 11: if ρ\rho has eigenvalues (pi)(p_i), then ipi=1\sum_i p_i=1.

Thus ρ\rho admits a spectral decomposition

ρ=ipiϕiϕi, \rho = \sum_i p_i\,|\phi_i\rangle\langle \phi_i|,

with (ϕi)(\phi_i) orthonormal and pi0p_i\ge 0, ipi=1\sum_i p_i=1.

Pure vs mixed

  • ρ\rho is a pure state iff it has rank 11, equivalently iff ρ2=ρ\rho^2=\rho, equivalently iff Tr(ρ2)=1\operatorname{Tr}(\rho^2)=1. (See .)
  • Otherwise ρ\rho is mixed and can be written as a convex combination ρ=kqkψkψk\rho=\sum_k q_k |\psi_k\rangle\langle\psi_k| with qk0q_k\ge 0, kqk=1\sum_k q_k=1. (See .)

Expectation values (Born rule in operator form)

If AA is an observable (a self-adjoint operator, see ), then the expectation value in state ρ\rho is

Eρ[A]=Tr(ρA). \mathbb{E}_\rho[A] = \operatorname{Tr}(\rho A).

This formula unifies pure and mixed states.

Dynamics and transformations (finite dimension)

  • Unitary evolution: if UU is unitary, then ρ\rho evolves as ρUρU\rho \mapsto U\rho U^\ast.
  • Projective measurement: spectral projectors PiP_i (as in ) define outcome probabilities Pr(i)=Tr(ρPi)\Pr(i)=\operatorname{Tr}(\rho P_i).

Information-theoretic quantities

Two common functionals of ρ\rho are:

  • Von Neumann entropy: S(ρ)=Tr(ρlogρ)S(\rho) = -\operatorname{Tr}(\rho\log\rho) (see ).
  • Relative entropy: D(ρσ)=Tr(ρ(logρlogσ))D(\rho\|\sigma)=\operatorname{Tr}(\rho(\log\rho-\log\sigma)) under suitable support conditions (see ).