2. Density Matrix Formalism

The density matrix formalism extends quantum mechanics to describe both pure and mixed states, providing a unified framework for statistical mixtures and quantum information theory.

Definition and Motivation

For a pure state $|\psi\rangle$, the density operator is:

$$\hat{\rho} = |\psi\rangle\langle\psi|$$

For a statistical mixture of states $|\psi_i\rangle$ with probabilities $p_i$:

$$\hat{\rho} = \sum_i p_i |\psi_i\rangle\langle\psi_i|$$

Where $\sum_i p_i = 1$ and $p_i \geq 0$

Key insight: Density matrix describes our incomplete knowledge about the quantum state

Properties of Density Matrices

A valid density matrix satisfies three conditions:

  • Hermiticity: $\hat{\rho}^\dagger = \hat{\rho}$ (ensures real eigenvalues)
  • Positive semidefinite: $\langle\phi|\hat{\rho}|\phi\rangle \geq 0$ for all $|\phi\rangle$
  • Unit trace: $\text{Tr}(\hat{\rho}) = 1$ (normalization)

Additional property distinguishing pure vs. mixed states:

$$\text{Tr}(\hat{\rho}^2) \begin{cases} = 1 & \text{pure state} \\ < 1 & \text{mixed state} \end{cases}$$

For pure states: $\hat{\rho}^2 = \hat{\rho}$ (idempotent)

Expectation Values and Measurements

The expectation value of an observable $\hat{A}$ is:

$$\langle A \rangle = \text{Tr}(\hat{\rho}\hat{A})$$

For pure state $\hat{\rho} = |\psi\rangle\langle\psi|$:

$$\langle A \rangle = \text{Tr}(|\psi\rangle\langle\psi|\hat{A}) = \langle\psi|\hat{A}|\psi\rangle$$

Recovers the standard quantum mechanical formula

Probability of measuring eigenvalue $a_n$:

$$P(a_n) = \text{Tr}(\hat{\rho}\hat{P}_n)$$

Where $\hat{P}_n = |a_n\rangle\langle a_n|$ is the projector onto eigenstate $|a_n\rangle$

Time Evolution

The density matrix evolves according to the von Neumann equation:

$$i\hbar\frac{\partial\hat{\rho}}{\partial t} = [\hat{H}, \hat{\rho}]$$

This is the quantum Liouville equation, analogous to classical Liouville theorem

Formal solution for time-independent Hamiltonian:

$$\hat{\rho}(t) = \hat{U}(t)\hat{\rho}(0)\hat{U}^\dagger(t)$$

Where $\hat{U}(t) = e^{-i\hat{H}t/\hbar}$ is the unitary time evolution operator

Conservation: $\text{Tr}(\hat{\rho}^2)$ is conserved under unitary evolution (pure states remain pure)

Reduced Density Matrix and Partial Trace

For a composite system AB with density matrix $\hat{\rho}_{AB}$, the reduced density matrix of subsystem A is:

$$\hat{\rho}_A = \text{Tr}_B(\hat{\rho}_{AB}) = \sum_i \langle b_i|\hat{\rho}_{AB}|b_i\rangle$$

Where $\{|b_i\rangle\}$ is a basis for subsystem B

Key property: Even if $\hat{\rho}_{AB}$ is pure, $\hat{\rho}_A$ can be mixed (signature of entanglement)

Example: Bell state $|\Psi^+\rangle = \frac{1}{\sqrt{2}}(|00\rangle + |11\rangle)$

$$\hat{\rho}_{AB} = |\Psi^+\rangle\langle\Psi^+| \quad \Rightarrow \quad \hat{\rho}_A = \frac{1}{2}(|0\rangle\langle 0| + |1\rangle\langle 1|) = \frac{\mathbb{I}}{2}$$

Subsystem A is maximally mixed despite AB being in a pure state

Purity and Mixedness

The purity of a state is quantified by:

$$\mathcal{P}(\hat{\rho}) = \text{Tr}(\hat{\rho}^2)$$

Properties:

  • $\frac{1}{d} \leq \mathcal{P} \leq 1$ for $d$-dimensional Hilbert space
  • $\mathcal{P} = 1$: pure state
  • $\mathcal{P} = 1/d$: maximally mixed state $\hat{\rho} = \mathbb{I}/d$

Alternative measure - linear entropy:

$$S_L = 1 - \text{Tr}(\hat{\rho}^2)$$

Measures the degree of mixedness: $S_L = 0$ for pure, $S_L = (d-1)/d$ for maximally mixed

Von Neumann Entropy

The quantum analog of Shannon entropy:

$$S(\hat{\rho}) = -\text{Tr}(\hat{\rho}\ln\hat{\rho}) = -\sum_i \lambda_i \ln\lambda_i$$

Where $\lambda_i$ are the eigenvalues of $\hat{\rho}$

Properties:

  • $S \geq 0$ with equality iff pure state
  • $S \leq \ln d$ with equality iff maximally mixed
  • Invariant under unitary transformations: $S(\hat{U}\hat{\rho}\hat{U}^\dagger) = S(\hat{\rho})$
  • Concave: $S(\sum_i p_i\hat{\rho}_i) \geq \sum_i p_i S(\hat{\rho}_i)$

For composite systems, satisfies subadditivity:

$$S(\hat{\rho}_{AB}) \leq S(\hat{\rho}_A) + S(\hat{\rho}_B)$$

Bloch Sphere Representation (Qubits)

For a qubit, any density matrix can be written as:

$$\hat{\rho} = \frac{1}{2}(\mathbb{I} + \vec{r} \cdot \vec{\sigma})$$

Where $\vec{r} = (r_x, r_y, r_z)$ is the Bloch vector and $\vec{\sigma} = (\sigma_x, \sigma_y, \sigma_z)$ are Pauli matrices

Geometric interpretation:

  • Pure states: $|\vec{r}| = 1$ (on surface of Bloch sphere)
  • Mixed states: $|\vec{r}| < 1$ (inside Bloch sphere)
  • Maximally mixed: $\vec{r} = 0$ (center of sphere)

Purity in terms of Bloch vector:

$$\text{Tr}(\hat{\rho}^2) = \frac{1}{2}(1 + |\vec{r}|^2)$$

Quantum Operations and CPTP Maps

General quantum operations are described by completely positive trace-preserving (CPTP) maps:

$$\mathcal{E}(\hat{\rho}) = \sum_k \hat{K}_k \hat{\rho} \hat{K}_k^\dagger$$

Where Kraus operators $\hat{K}_k$ satisfy $\sum_k \hat{K}_k^\dagger \hat{K}_k = \mathbb{I}$

Examples:

  • Unitary evolution: Single Kraus operator $\hat{K} = \hat{U}$
  • Measurement: $\hat{K}_k = \hat{P}_k$ (projectors)
  • Decoherence: Multiple Kraus operators describing environment interaction

Applications

1. Thermal states

Canonical ensemble at temperature $T$:

$$\hat{\rho} = \frac{e^{-\beta\hat{H}}}{Z}, \quad Z = \text{Tr}(e^{-\beta\hat{H}})$$

Where $\beta = 1/(k_B T)$

2. Open quantum systems

System interacting with environment evolves non-unitarily:

$$\frac{d\hat{\rho}_S}{dt} = -\frac{i}{\hbar}[\hat{H}_S, \hat{\rho}_S] + \mathcal{L}[\hat{\rho}_S]$$

Lindblad master equation, where $\mathcal{L}$ describes dissipation and decoherence

3. Quantum information processing

  • Quantum channels and communication protocols
  • Entanglement quantification via entropy
  • Quantum error correction (tracking mixed states through noisy channels)
  • Quantum tomography (reconstructing $\hat{\rho}$ from measurements)

Physical Insight: The density matrix formalism is essential for describing realistic quantum systems that are never perfectly isolated. It provides a natural bridge between quantum mechanics and statistical mechanics, and is the foundation for quantum information theory and the study of open quantum systems.