Density matrix

From Wikipedia, the free encyclopedia

"Mixed state" redirects here. For the psychiatric condition, see Mixed state (psychology).

In quantum mechanics, a density matrix is a self-adjoint (or Hermitian) positive-semidefinite matrix, (possibly infinite dimensional), of trace one, that describes the statistical state of a quantum system. The formalism was introduced by John von Neumann (and according to other sources, independently by Lev Landau and Felix Bloch) in 1927.

It is the quantum-mechanical analogue to a phase-space probability measure (probability distribution of position and momentum) in classical statistical mechanics. The need for a statistical description via density matrices arises when one considers either an ensemble of systems, or one system when its preparation history is uncertain and one does not know with 100% certainty which pure quantum state the system is in.

Situations in which a density matrix is used include the following: a quantum system in thermal equilibrium (at finite temperatures); nonequilibrium time-evolution that starts out of a mixed equilibrium state; and entanglement between two subsystems, where each individual system must be described, via the partial trace operation, by a density matrix even though the complete system may be in a pure state; and in analysis of quantum decoherence. See also quantum statistical mechanics.

The operator that is represented by the density matrix is called the density operator. (The close relationship between matrices and operators is a basic concept in linear algebra; see the article Linear operator for details.) In practice, the terms "density matrix" and "density operator" are often used interchangeably. The density operator, like the density matrix, is positive-semidefinite, self-adjoint, and has trace one.

Contents

The need for a statistical description

In quantum mechanics, the state vector ψ of a system completely determines the statistical behavior of an observable O. This means that if O is represented by an operator A on the Hilbert space H of the system, then for any real-valued function F defined on the real numbers, the expectation value of F(O) is the quantity

 \langle \psi | F(A) \psi \rangle.

or written as

 \langle \psi | F(A) | \psi \rangle.

in Dirac notation.

Now consider the example of a "mixed quantum system" prepared by statistically combining two different pure states φ, ψ each with probability 1/2. The preparation process for such a system consists in tossing an unbiased coin and using the preparation process for φ or for ψ depending on whether the toss outcome is heads or tails.

It is not hard to show that the statistical properties of the observable O for the system prepared in such a mixed state are completely determined. However, there is no vector ξ which determines this statistical behavior in the sense that the expectation value of F(O) is

 \langle \xi | F(A) | \xi \rangle.

Nevertheless: there is a unique operator ρ such that the expectation value can be written as

 \operatorname{Tr}(F(A) \rho).

The operator ρ is the density operator of the mixed system. A simple calculation shows that for the example mentioned above:

 \rho = \frac{1}{2} | \phi\rangle \langle \phi | + \frac{1}{2} | \psi\rangle \langle \psi |.

Formulation

For a finite dimensional Hilbert space, the most general density operator is of the form

 \rho = \sum_j p_j |\psi_j \rang \lang \psi_j|,
|\psi_j \rang

If the system is not closed, however, then it is simply not correct to claim that it has some definite but unknown state vector, as the density operator may record physical entanglements to other systems.

Example Consider a quantum ensemble of size N with occupancy numbers n1, n2,...,nk corresponding to the orthonormal states |1>,...,|k>, respectively, where n1+...+nk = N, and, thus, the coefficients pj = nj /N. For a pure ensemble, where all N particles are in state |i>, we have nj = 0, for all ji, from which we recover the corresponding density operator ρ = |i >< i|.

However the density operator of a mixed state does not capture all the information about a mixture; in particular, the coefficients pj and the kets ψj are not recoverable from the operator ρ without additional information. This non-uniqueness implies that different ensembles or mixtures may correspond to the same density operator. Such equivalent ensembles or mixtures cannot be distinguished by measurement of observables alone. This equivalence can be characterized precisely. Two ensembles ψ, ψ' define the same density operator if and only if there is a matrix U with

 \sum_k u^*_{ik} u_{kj} = \delta_{ij}

i.e, U is unitary and such that

 | \psi_i'\rangle \sqrt p_i' = \sum_{j} u_{ij} | \psi_j\rangle \sqrt {p_j}.

This is simply a restatement of the following fact from linear algebra: for two square matrices M and N, M M* = N N* if and only if M = NU for some unitary U. (See square root of a matrix for more details.) Thus there is an unitary freedom in the ket mixture or ensemble that gives the same density operator. However if the kets in the mixture are orthonormal then the original probabilities pj are recoverable as the eigenvalues of the density matrix.

In operator language, a density operator is a positive semidefinite, hermitian operator acting on the state space of trace 1. A density operator describes a pure state if it is a rank one projection. Equivalently, a density operator ρ is a pure state if and only if

\; \rho = \rho^2

i.e. the state is idempotent. This is true regardless of whether H is finite dimensional or not.

Geometrically, when the state is not expressible as a convex combination of other states, it is a pure state. The family of mixed states is a convex set and a state is pure if it is an extremal point of that set.

It follows from the spectral theorem for compact self-adjoint operators that every mixed state is an infinite convex combination of pure states. This representation is not unique. Furthermore, a theorem of Andrew Gleason states that certain functions defined on the family of projections and taking values in [0,1] (which can be regarded as quantum analogues of probability measures) are determined by unique mixed states. See quantum logic for more details.

Measurement

|\psi_j\rang
\rho = \sum_j p_j |\psi_j \rang \lang \psi_j| .

The expectation value of the measurement can be calculated by extending from the case of pure states (see Measurement in quantum mechanics):

 \lang A \rang = \sum_j p_j \lang \psi_j|A|\psi_j \rang = \operatorname{tr}[\rho A],

where tr denotes trace. Moreover, if A has spectral resolution

A = \sum_i a_i |a_i \rang \lang a_i| = \sum _i a_i P_i,
P_i = |a_i \rang \lang a_i|
\; \rho ^' = \sum_i P_i \rho P_i.

Note that the above density operator describes the full ensemble after measurement. The sub-ensemble for which the measurement result was the particular value ai is described by the different density operator

\rho_i' = \frac{P_i \rho P_i}{\operatorname{tr}[\rho P_i]}.
|a_i\rang

Entropy

The von Neumann entropy S of a mixture can be expressed in terms of the probabilities pi or in terms of the trace and logarithm of the density operator ρ:

S = -\sum_i p_i \ln \,p_i = -\operatorname{tr}[\rho \ln \rho].

This entropy can increase but never decrease with a measurement. The entropy of a pure state is zero, while that of a proper mixture always greater than zero. Therefore a pure state may be converted into a mixture by a measurement, but a proper mixture can never be converted into a pure state. Thus the act of measurement induces a fundamental irreversible change on the density matrix; this is analogous to the "collapse" of the state vector, or wavefunction collapse.

Von Neumann equation

Just as the Schrödinger equation describes how pure states evolve in time, the von Neumann Equation describes how a density operator evolves in time (in fact, the two equations are equivalent, in the sense that either can be derived from the other.) The von Neumann equation states:

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

where the brackets denote a commutator. (Note that this equation is only true when the density operator is taken to be in the Schrödinger picture, as is usually the case.) This equation is analogous to the Liouville equation in classical physics.

If the Hamiltonian is time-independent, this differential equation can be solved to get

\rho(t) = e^{-i H t/\hbar} \rho(0) e^{i H t/\hbar}.

C*-algebraic formulation of states

It is now generally accepted that the description of quantum mechanics in which all self-adjoint operators represent observables is untenable. For this reason, observables are identified to elements of an abstract C*-algebra A (that is one without a distinguished representation as an algebra of operators) and states are positive linear functionals on A. Note that by using the GNS construction, we can recover Hilbert spaces which realize A as an algebra of operators.

Geometrically, a pure state on a C*-algebra A is a state which is an extreme point of the set of all states on A. By properties of the GNS construction these states correspond to irreducible representations of A.

The states of the C*-algebra of compact operators K(H) correspond exactly to the density operators and therefore the pure states of K(H) are exactly the pure states in the sense of quantum mechanics.

The C*-algebraic formulation can be seen to include both classical and quantum systems. When the system is classical, the algebra of observables become an abelian C*-algebra. In that case the states become probability measures, as noted in the introduction.

Notes and references

  1. ^ We assume ψ is a pure state in this example.
  2. ^ Technically, F must be a Borel function
  3. ^ F(O) is defined to be the result of measuring O and then applying F to the outcome.
  4. ^ Nielsen, Michael & Chuang, Isaac (2000), Quantum Computation and Quantum Information, Cambridge University Press, ISBN 978-0-521-63503-5 . Chapter 11: Entropy and information, Theorem 11.9, "Projective measurements cannot decrease entropy"
  5. ^ Everett, Hugh (1973), "The Theory of the Universal Wavefunction (1956) Appendix I. "Monotone decrease of information for stochastic processes"", The Many-Worlds Interpretation of Quantum Mechanics, Princeton Series in Physics, Princeton University Press, pp. 128–129, ISBN 978-0-691-08131-1 
  6. ^ See appendix, Mackey, George Whitelaw (1963), Mathematical Foundations of Quantum Mechanics, Dover Books on Mathematics, New York: Dover Publications, ISBN 978-0-486-43517-6 
  7. ^ Emch, Gerard G. (1972), Algebraic methods in statistical mechanics and quantum field theory, Wiley-Interscience, ISBN 978-0-471-23900-0