Density Matrix — Physica

Equivalent forms

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

i\hbar\,\partial_t\hat{\rho} = [\hat{H},\hat{\rho}]

One Hermitian, positive, unit-trace matrix encodes every possible quantum state — pure, mixed, classical, entangled.

Symbol	Name	SI unit	Dimension	Range
$rho$	Density operatoroutput Hermitian, positive, unit-trace operator on Hilbert space	1	1	0 – 1
$p_i$	Classical probability Statistical weight of pure state \|psi_i>	1	1	0 – 1
$purity$	Purity Tr(rho^2)output 1 for pure, 1/d for maximally mixed in dimension d	1	1	0.5 – 1
$T$	Bath temperature Temperature determining thermal mixture weights	K	Theta	1e-9 – 1000

Unit systems

Symbol	SI	CGS	Imperial
$rho$	1	1	1
$p_i$	1	1	1
$purity$	1	1	1
$T$	K	K	R

Where it holds

Universal — pure, mixed, finite or infinite dimension. Required whenever the system is part of a larger system (reduced density matrix via partial trace).

Dimensional analysis

\rho \to

dimensionless

sum of p_i (dimensionless) * |psi><psi| (dimensionless

) \to

dimensionless

rho is dimensionless.

\checkmark

Discovery

John von Neumann / Lev Landau · 1927

Landau introduced the operator in 1927 to describe subsystems; von Neumann formalized it in his 1932 'Mathematical Foundations of Quantum Mechanics'. It is now the universal language of open systems, quantum information, and decoherence.

Try this

How do you describe a quantum system you don't fully know?

Pure states use a wavefunction, but mixed and decohering systems need more. The density matrix is the universal quantum descriptor — and the language of every real experiment.

Research status: stable

Real-world applications

Quantum state tomography
Open quantum systems (Lindblad master equation)
Quantum error correction
Decoherence models in quantum computing and cosmology

Common misconceptions

Mixed states are not superpositions — they describe classical uncertainty, not quantum interference
Two different ensembles can produce the same rho — only rho is observable
The partial trace of an entangled pure state is mixed — even when the global state is pure

Experimental verification

Quantum state tomography (1989) reconstructs rho from many measurements. Entanglement witnesses (concurrence, negativity) use rho directly. Quantum process tomography validates gates in IBM, Google, IonQ devices.

Derivation

Given an ensemble {p_i, |psi_i>} where p_i is the classical probability of preparing |psi_i>, expectation of A is <A>

= sum_i p_i

<psi_i|A|psi_i>

= Tr(\rho

*A) with

\rho = sum_i p_i |\psi _i

><psi_i|. rho is Hermitian, positive semidefinite, with

Tr(\rho ) = 1

.

Time evolution under H follows the von Neumann equation i*hbar*

d_t*\rho = [H

, rho]; with environment coupling, it becomes the Lindblad master equation.

Limiting cases

\rho = |\psi

><psi|⟶ Pure state,

Tr(\rho ^2) = 1

Maximally informed — equivalent to a single wavefunction.

\rho = I/d

⟶ Maximally mixed,

Tr(\rho ^2) = 1/d

Complete ignorance — all outcomes equally likely.

Off-diagonal elements

\to 0

⟶ Decoherence to classical mixtureEnvironment kills coherences but leaves diagonal probabilities.

What if…

What if rho's off-diagonals vanish instantly?

Pure decoherence — system is effectively classical; no interference, no Bell violations.

What if

Tr(\rho ) \neq 1

?

rho is not normalized — probabilities don't sum to 1, the description is ill-defined.

What if rho has a negative eigenvalue?

Not a valid quantum state — violates positivity; basis of entanglement witnesses (PPT criterion).

1

Purity of a 2-level thermal state at T = 300 K

Given ·

level spacing:: 1.6e-21
T:: 300

Find · Tr(rho^2)

Steps

Boltzmann: $p_0 = 1/(1+e^(-\beta *E))$ , $p_1 = e^(-\beta *E)/(1+e^(-\beta *E))$
$p_0 \approx 0.595$ , $p_1 \approx 0.405$
$Tr(\rho ^2) = p_0^2 + p_1^2 \approx 0.518$ — neither pure nor maximally mixed

Result ·

\beta *E = 1.6e-21 / (1.38e-23 * 300) \approx 0.386

;

p_0 \approx 0.595

,

p_1 \approx 0.405

;

Tr(\rho ^2) \approx 0.518

2

Reduced density matrix of Bell state

Given ·

Phi+:: (|00> + |11>)/sqrt(2)

Find · rho_A

Steps

Full $\rho = (|00$ ><00| + |00><11| + |11><00| + |11><11|)/2
Trace out B: $\rho _A = (|0$ ><0| + |1>< $1|)/2 = I/2$
Subsystem is maximally mixed even though global state is pure — hallmark of entanglement

Result · Tr_B(|Phi+><

\Phi +|) = I/2

— maximally mixed

Spin-1/2 Operators→Time-Dependent Schrödinger Equation→Heisenberg Uncertainty Principle→