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.

Runnable Simulations

Density Matrix: Pure vs Mixed States

Python
script.py65 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Partial Trace and Reduced Density Matrix

Fortran
program.f9074 lines

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Rate this chapter: