Skip to content

ρ

The lowercase Greek letter ρ\rho, pronounced “rho,” denotes a density operator by default in quantum-state contexts. The same glyph also commonly denotes a probability density, number or charge density, density of states, radial coordinate, or resistivity. Its arguments, indices, units, and mathematical type are decisive.

For a quantum state, ρ\rho is a positive trace-one operator:

ρ≥0,Tr⁡ρ=1.\rho\ge0, \qquad \operatorname{Tr}\rho=1.

In finite dimensions,

ρ†=ρ,ρ=∑jλj∣j⟩⟨j∣,λj≥0,∑jλj=1.\rho^\dagger=\rho, \qquad \rho=\sum_j\lambda_j \lvert j\rangle\langle j\rvert, \qquad \lambda_j\ge0, \qquad \sum_j\lambda_j=1.

The abstract density operator is basis independent. Its matrix elements in a basis {∣n⟩}\{\lvert n\rangle\} are

ρmn=⟨m∣ρ∣n⟩.\rho_{mn} =\langle m\vert\rho\vert n\rangle.

The phrase density matrix often refers informally to the operator itself, but a matrix is one representation of that operator.

UseMathematical typeTypical normalization or units
ρ\rhoPositive trace-class state operatorDimensionless; Tr⁡ρ=1\operatorname{Tr}\rho=1
ρmn\rho_{mn}Matrix element in a discrete orthonormal basisDimensionless
ρ(x,x′)\rho(x,x')Coordinate kernel ⟨x∣ρ∣x′⟩\langle x\vert\rho\vert x'\rangleL−1L^{-1} in one dimension
ρ(x)\rho(x)Position probability densityL−1L^{-1} in one dimension
n(x)n(\mathbf x) or ρN(x)\rho_N(\mathbf x)Number densityL−3L^{-3}
ρq(x)\rho_q(\mathbf x)Charge densitycharge per volume
ρm(x)\rho_m(\mathbf x)Mass densitymass per volume
ρ(E)\rho(E)Density of statesstates per energy, possibly also per volume
ρ\rho in cylindrical coordinatesRadial distance from an axislength
ρel\rho_{\mathrm{el}}Electrical resistivityΩ m\Omega\,\mathrm m in SI

The coordinate kernel’s units follow from

Tr⁡ρ=∫ρ(x,x) dx=1.\operatorname{Tr}\rho =\int\rho(x,x)\,dx=1.

Thus a dimensionless abstract operator can have dimensionful continuous-basis matrix elements because generalized basis kets and integration measures carry dimensions.

FormMeaning
ρ(t)\rho(t)Time-dependent density operator
ρA\rho_AState reduced to subsystem AA
ρAB\rho_{AB}Joint state of subsystems AA and BB
ρss\rho_{\mathrm{ss}}Steady-state density operator
ρβ\rho_\betaThermal state at inverse temperature β\beta
ρin\rho_{\mathrm{in}}, ρout\rho_{\mathrm{out}}Input and output states
ρ(k)\rho^{(k)}Member of an ensemble or indexed family
ρmn\rho_{mn}Basis matrix element, not usually a subsystem label

Subsystem subscripts and matrix indices play different roles. ρA\rho_A is itself an operator on HA\mathcal H_A, while ρmn\rho_{mn} is a scalar component after choosing a basis.

For an observable AA,

⟨A⟩=Tr⁡(ρA).\langle A\rangle =\operatorname{Tr}(\rho A).

For a measurement effect EkE_k,

p(k)=Tr⁡(ρEk).p(k)=\operatorname{Tr}(\rho E_k).

For a bipartite state,

ρA=Tr⁡B(ρAB).\rho_A=\operatorname{Tr}_B(\rho_{AB}).

Closed-system evolution is

iℏdρdt=[H,ρ],i\hbar\frac{d\rho}{dt} =[H,\rho],

or

ρ(t)=U(t,t0)ρ(t0)U†(t,t0).\rho(t) =U(t,t_0)\rho(t_0)U^\dagger(t,t_0).

These equations identify ρ\rho as an operator. A scalar spatial density does not enter a commutator in this way.

A pure state has

ρ=∣ψ⟩⟨ψ∣,ρ2=ρ,Tr⁡(ρ2)=1.\rho =\lvert\psi\rangle\langle\psi\rvert, \qquad \rho^2=\rho, \qquad \operatorname{Tr}(\rho^2)=1.

A mixed state has Tr⁡(ρ2)<1\operatorname{Tr}(\rho^2)<1. A thermal Gibbs state is often written

ρβ=e−βHZ,Z=Tr⁡(e−βH),\rho_\beta =\frac{e^{-\beta H}}{Z}, \qquad Z=\operatorname{Tr}(e^{-\beta H}),

where

β=1kBT\beta=\frac{1}{k_{\mathrm B}T}

unless thermal units kB=1k_{\mathrm B}=1 have been declared.

The density operator in position representation has kernel

ρ(x,x′)=⟨x∣ρ∣x′⟩.\rho(x,x') =\langle x\vert\rho\vert x'\rangle.

For a pure state,

ρ(x,x′)=ψ(x)ψ∗(x′).\rho(x,x') =\psi(x)\psi^*(x').

The diagonal

ρ(x,x)=∣ψ(x)∣2\rho(x,x)=\lvert\psi(x)\rvert^2

is the position probability density for a scalar pure wavefunction. The off-diagonal kernel contains spatial coherence information.

For a mixed state,

ρ(x,x′)=∑jpjψj(x)ψj∗(x′),\rho(x,x') =\sum_jp_j \psi_j(x)\psi_j^*(x'),

for any chosen ensemble representation. The operator is independent of which ensemble decomposition is used.

In wave mechanics, authors sometimes write

ρprob(x,t)=∣ψ(x,t)∣2,\rho_{\mathrm{prob}}(x,t) =\lvert\psi(x,t)\rvert^2,

with

∫ρprob(x,t) dx=1.\int\rho_{\mathrm{prob}}(x,t)\,dx=1.

For a spinor, the density includes a component sum:

ρprob(x,t)=ψ†(x,t)ψ(x,t).\rho_{\mathrm{prob}}(\mathbf x,t) =\psi^\dagger(\mathbf x,t)\psi(\mathbf x,t).

A particle-number density instead often satisfies

∫n(x) d3x=N.\int n(\mathbf x)\,d^3x=N.

The notation ρ(x)\rho(\mathbf x) does not reveal whether the integral is one, total charge, total mass, or particle number. Units and the defining integral do.

The notation

ρ(E)\rho(E)

often denotes a density of energy levels. In a discrete formal form,

ρ(E)=∑nδ(E−En).\rho(E) =\sum_n\delta(E-E_n).

It has units of inverse energy. A density of states per volume has units of inverse energy per volume. This ρ(E)\rho(E) is a scalar distribution, not a quantum state operator.

In cylindrical coordinates,

ρ=x2+y2\rho=\sqrt{x^2+y^2}

is commonly the distance from the symmetry axis. The volume element is

d3x=ρ dρ dϕ dz.d^3x=\rho\,d\rho\,d\phi\,dz.

Here ρ\rho has units of length and appears without an operator trace. Some texts instead use r⊥r_\perp to avoid collision with density notation.

Density Operators owns the default quantum-state meaning. Density-Matrix Conventions defines subscripts, matrix elements, ensembles, and Bloch-vector notation.

Reduced Density Matrices owns ρA\rho_A and ρAB\rho_{AB}. The Density Matrix Glossary Entry gives the concept-level summary, and the Density-Matrix Expectation Formula Card provides the compact trace rule.

  • ϱ\varrho is sometimes used for a density operator when ρ\rho is needed for a spatial density.
  • ρ^\hat\rho makes operator status explicit; this documentation usually omits the hat when the context is unambiguous.
  • DD or WW can denote density matrices in some older or specialized sources.
  • ρ(1)\rho^{(1)} and ρ(2)\rho^{(2)} can denote reduced one- and two-body density matrices rather than powers.
  • Do not confuse a density operator ρ\rho with a spatial density ρ(x)\rho(x).
  • Matrix entries of ρ\rho depend on the chosen basis.
  • The notation “density matrix” is common even when the abstract object is a density operator.
  • A density operator is dimensionless, but its continuous-basis kernel can carry units.
  • ρ(x,x)\rho(x,x) is a diagonal kernel; ρ(x)\rho(x) may instead be an independently defined scalar density.
  • A subsystem label such as AA is not a matrix index.
  • Superscripts in ρ(k)\rho^{(k)} are not automatically powers.
  • A diagonal density matrix is diagonal only relative to the stated basis.
  • ρ(E)\rho(E) as a density of states is not positive trace one as an operator.
  • The cylindrical radial coordinate ρ\rho is distinguished by its geometric definition and length units.
  • A number, charge, or mass density integrates to the corresponding total, not necessarily to one.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014, chs. 2 and 3.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed., Cambridge University Press, 2010, sec. 2.4.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994, ch. 12.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.