Density matrices and mixed preparations
Combine preparation probabilities at the operator level and distinguish statistical mixtures from coherent superpositions in a chosen basis.
On this page
The question for today
What you’ll understand
- Construct a density operator from a probabilistic preparation recipe.
- Check positivity, trace, and purity of a finite density matrix.
- Explain why one density operator can have several ensemble decompositions.
A preparation device may select among pure states using a classical random record that is unavailable to the observer. Predicting outcomes then requires averaging probabilities over that record, rather than adding amplitudes across distinct preparations. The density operator collects these predictions into one positive unit-trace object. Its matrix depends on the chosen basis, but the operator's spectrum, purity, and predictions for specified measurements do not depend on that representation.
Average projectors instead of state vectors
If preparation j occurs with probability wj and produces a normalized ket, average the rank-one projectors weighted by those probabilities. Averaging kets directly is incorrect: independent preparation labels are classical alternatives, so no controlled relative phase links them. The resulting operator is positive because every vector has a nonnegative quadratic form against each projector. Its trace is one because the recipe weights sum to one. These properties characterize finite-dimensional density operators.
For an observable A, expectation is the trace of rho times A. For a projective outcome, replace A by its projector to obtain a Born probability. In infinite dimension, rho must be a positive trace-class unit-trace operator, and unbounded-observable moments need additional finiteness conditions. A matrix that is Hermitian with trace one but has a negative eigenvalue is not a valid density matrix: some possible measurement would receive a negative probability.
Coherence is basis dependent, purity is not
Off-diagonal matrix entries encode phase-sensitive predictions relative to the selected basis. A coherent superposition can have nonzero off-diagonal entries even when its diagonal probabilities match those of an incoherent mixture. Diagonalizing any density operator removes those entries in its eigenbasis, so having a diagonal matrix does not by itself mean a preparation lacks coherence relative to every physically relevant basis. State properties and basis-specific resources must be described separately.
Purity is the trace of rho squared. For a finite d-dimensional normalized state it lies between one divided by d and one. Equality to one means a rank-one projector, which describes a pure state ray; smaller purity means a mixed state. In a qubit, the Bloch-vector length provides the same test. Pure states lie on the unit sphere, mixed states inside it, and the maximally mixed state sits at its centre.
The recipe is not uniquely recoverable from rho
A maximally mixed qubit can be prepared by selecting the two computational states equally, or by selecting the plus and minus X eigenstates equally. Both recipes produce the same operator and therefore identical statistics for every measurement on the qubit alone. Their discarded classical labels are different records, but those records are not recoverable from rho without additional access to a preparation register. This is ensemble nonuniqueness rather than an uncertainty about which matrix entries to use.
The same mathematical density operator can also describe a subsystem entangled with another system. Locally, its measurement predictions match those of an appropriate classical preparation mixture, even though the global correlations differ. Calling every mixed state merely ignorance of a unique pre-existing pure state therefore adds information that the operator does not contain. The next lessons introduce the tensor-product description needed to distinguish a local state from its possible global extensions.
Work it through
A nonorthogonal preparation mixture
A device prepares |0⟩ with probability 3/5 and with probability 2/5, discarding the label. Find ρ, the zero and plus outcome probabilities, and its purity.
Use projectors: and . The weighted sum is , with unit trace and positive determinant 3/25.
The computational zero probability is ρ00=4/5. Projecting onto plus gives . This differs from using the diagonal entries alone because the X measurement is sensitive to off-diagonal coherence.
Square the matrix and take its trace: . Its Bloch vector is (2/5,0,3/5), whose squared length 13/25 gives the same purity.
The mixture predicts p(0)=0.80 and p(+)=0.70, with purity 0.76. It is mixed despite retaining off-diagonal entries in the computational basis.
The living notebook
Explore the model
Interactive lab · quantum
Two pure components do not display their mixture
Reset and set both controls to θ=90° and φ=0°; predict Bloch (1,0,0). For a fresh comparison, Reset and set both controls to θ=90° and φ=180°; predict Bloch (−1,0,0). Both runs have equal computational probabilities. An equally weighted preparation mixture has average Bloch vector zero, which must be calculated separately and is not displayed by either pure-state run.
- Probability |0⟩
- 71.1%
- Probability |1⟩
- 28.9%
- Bloch (x, y, z)
- 0.78, 0.45, 0.42
A normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forThe instrument plots pure two-level geometry only. It has no mixture-weight control, interior mixed-state point, random preparation sampling, or purity readout.
Pause & reason
Check your understanding
Your turn
Try it yourself
Prepare |+⟩ and |−⟩ with equal classical probabilities, where . Find the density matrix, its purity, and probabilities for computational zero and X-plus.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Add the two projectors with weights one half. Their off-diagonal entries have opposite signs.
Compare with the solution
The weighted projectors sum to . Hence . Both p(0) and p(+) equal 1/2. This is different from adding the kets coherently, whose normalized sum would be |0⟩ and would give p(0)=1. The mixed recipe is operationally identical on the qubit to an equal mixture of |0⟩ and |1⟩, although its preparation labels describe different pure components.
Keep the idea
Take it with you
- Classical preparation probabilities weight projectors, rather than coherent amplitudes.
- A valid density operator is positive and unit trace; purity detects a rank-one state.
- Ensemble decompositions are nonunique even though all specified local predictions are fixed.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.