 Density matrix

"Mixed state" redirects here. For the psychiatric condition, see Mixed state (psychiatry).
In quantum mechanics, a density matrix is a selfadjoint (or Hermitian) positivesemidefinite matrix (possibly infinite dimensional) of trace one, that describes the statistical state of a quantum system.^{[1]} The formalism was introduced by John von Neumann^{[2]} (and, less systematically, independently by Lev Landau and Felix Bloch in 1927).^{[3]}^{[4]}
The density matrix is useful for describing and performing calculations with a mixed state, which is a statistical ensemble of several quantum states. (This is in contrast to a pure state, which is a quantum system that is described by a single state vector). The density matrix is the quantummechanical analogue to a phasespace probability measure (probability distribution of position and momentum) in classical statistical mechanics.
Mixed states arise in situations where there is classical uncertainty, i.e. when it is unknown by the experimenter which particular states are being manipulated. (This should not be confused with the unrelated concept of quantum uncertainty, which dictates that even if the experimenter knows which particular states are being manipulated, the results of some measurements cannot be predicted.) For example: A quantum system with an uncertain or randomlyvarying preparation history (so one does not know with certainty which pure quantum state the system is in); or a quantum system in thermal equilibrium (at finite temperatures). Also, if a quantum system has two or more subsystems that are entangled, then each individual subsystem must be treated as a mixed state even if the complete system is in a pure state; the density matrix of the subsystem is calculated as a partial trace of the density matrix of the whole system. Relatedly, the density matrix is a crucial tool in quantum decoherence theory. 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. (When the space of quantum states is infinitedimensional, the term "density operator" is preferred.) The density operator, like the density matrix, is positivesemidefinite, selfadjoint, and has trace one.
Contents
Pure and mixed states
In quantum mechanics, a quantum system is represented by a state vector (or ket) . A quantum system with a state vector is called a pure state. However, it is also possible for a system to be in a statistical ensemble of different state vectors: For example, there may be a 50% probability that the state vector is and a 50% chance that the state vector is . This system would be in a mixed state. The density matrix is especially useful for mixed states, because any state, pure or mixed, can be characterized by a single density matrix.
A mixed state is different from a quantum superposition. In fact, a quantum superposition of pure states is another pure state, for example .
Example: Light polarization
An example of pure and mixed states is light polarization. Photons can have two helicities, corresponding to two orthogonal quantum states, (right circular polarization) and (left circular polarization). A photon can also be in a superposition state, such as (vertical polarization) or (horizontal polarization). More generally, it can be in any state , corresponding to linear, circular, or elliptical polarization.
However, unpolarized light (such as the light from an incandescent light bulb) is different from any of these. Unlike linearly or elliptically polarized light, it passes through a polarizer with 50% intensity loss whatever the orientation of the polarizer; and unlike circularly polarized light, it cannot be made linearly polarized with any wave plate. Indeed, unpolarized light cannot be described as any state of the form . However, unpolarized light can be described perfectly by assuming that each photon is either with 50% probability or with 50% probability. The same behavior would occur if each photon was either vertically polarized with 50% probability or horizontally polarized with 50% probability.
Therefore, unpolarized light cannot be described by any pure state, but can be described as a statistical ensemble of pure states in at least two ways (the ensemble of half left and half right circularly polarized, or the ensemble of half vertically and half horizontally linearly polarized). These two ensembles are completely indistinguishable experimentally, and therefore they are considered the same mixed state. One of the advantages of the density matrix is that there is just one density matrix for each mixed state, whereas there are many statistical ensembles of pure states for each mixed state. Nevertheless, the density matrix contains all the information necessary to calculate any measurable property of the mixed state.
Where do mixed states come from? To answer that, consider how to generate unpolarized light. One way is to use a system in thermal equilibrium, a statistical mixture of enormous numbers of microstates, each with a certain probability (the Boltzmann factor), switching rapidly from one to the next due to thermal fluctuations. Thermal randomness explains why an incandescent light bulb, for example, emits unpolarized light. A second way to generate unpolarized light is to introduce uncertainty in the preparation of the system, for example, passing it through a birefringent crystal with a rough surface, so that slightly different parts of the beam acquire different polarizations. A third way to generate unpolarized light uses an EPR setup: A radioactive decay can emit two photons traveling in opposite directions, in the quantum state . The two photons together are in a pure state, but if you only look at one of the photons and ignore the other, the photon behaves just like unpolarized light.
More generally, mixed states commonly arise from a statistical mixture of the starting state (such as in thermal equilibrium), from uncertainty in the preparation procedure (such as slightly different paths that a photon can travel), or from looking at a subsystem entangled with something else.
Mathematical description
In quantum mechanics, the state vector ^{[5]} of a system completely determines the statistical behavior of a measurement. As an example, take an observable quantity, and let A be the associated observable operator which has a representation on the Hilbert space of the quantum system. For any realvalued function F^{[6]} defined on the real numbers, the expectation value of F(A) ^{[7]} is the quantity
Now consider the example of a "mixed quantum system" prepared by statistically combining two different pure states and , with the associated probabilities p and 1 − p, respectively. The associated probabilities mean that the preparation process for the quantum system ends in the state with probability p and in the state with probability 1 − p.
It is not hard to show that the statistical properties of the observable for the system prepared in such a mixed state are completely determined. However, there is no state vector which determines this statistical behavior in the sense that the expectation value of F(A) is
Nevertheless, there is a unique operator ρ such that the expectation value of F(A) can be written as
where the operator ρ is the density operator of the mixed system. A simple calculation shows that the operator ρ for the above example is given by
Formulation
For a finite dimensional function space, the most general density operator is of the form
where the coefficients p_{j} are nonnegative and add up to one. This represents a statistical mixture of pure states. If the given system is closed, then one can think of a mixed state as representing a single system with an uncertain preparation history, as explicitly detailed above; or we can regard the mixed state as representing an ensemble of systems, i.e. large number of copies of the system in question, where p_{j} is the proportion of the ensemble being in the state . An ensemble is described by a pure state if every copy of the system in that ensemble is in the same state, i.e. it is a pure ensemble.
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 n_{1}, n_{2},...,n_{k} corresponding to the orthonormal states , respectively, where n_{1}+...+n_{k} = N, and, thus, the coefficients p_{j} = n_{j} /N. For a pure ensemble, where all N particles are in state , we have n_{j} = 0, for all j ≠ i, from which we recover the corresponding density operator .
However the density operator of a mixed state does not capture all the information about a mixture; in particular, the coefficients p_{j} and the kets ψ_{j} are not recoverable from the operator ρ without additional information. This nonuniqueness 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
 U ^{*} U = I
i.e., U is unitary and such that
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 a 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 p_{j} 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
 ,
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 selfadjoint 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
Let A be an observable of the system, and suppose the ensemble is in a mixed state such that each of the pure states occurs with probability p_{j}. Then the corresponding density operator is:
The expectation value of the measurement can be calculated by extending from the case of pure states (see Measurement in quantum mechanics):
where tr denotes trace. Moreover, if A has spectral resolution
where , the corresponding density operator after the measurement is given by:
Note that the above density operator describes the full ensemble after measurement. The subensemble for which the measurement result was the particular value a_{i} is described by the different density operator
This is true assuming that is the only eigenket (up to phase) with eigenvalue a_{i}; more generally, P_{i} in this expression would be replaced by the projection operator into the eigenspace corresponding to eigenvalue a_{i}.
Entropy
The von Neumann entropy S of a mixture can be expressed in terms of the eigenvalues of ρ or in terms of the trace and logarithm of the density operator ρ. Since ρ is a positive semidefinite operator, it has a spectral decomposition such that where are orthonormal vectors. Therefore the entropy of a quantum system with density matrix ρ is
Also it can be shown that
where H(p) is the Shannon entropy. This entropy can increase but never decrease^{[8]}^{[9]} 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.
(A subsystem of a larger system can be turned from a mixed to a pure state, but only by increasing the von Neumann entropy elsewhere in the system. This is analogous to how the entropy of an object can be lowered by putting it in a refrigerator: The air outside the refrigerator's heatexchanger warms up, gaining even more entropy than was lost by the object in the refrigerator. See second law of thermodynamics.)
Von Neumann equation
Just as the Schrödinger equation describes how pure states evolve in time, the von Neumann equation (also known as Liouvillevon 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:^{[10]}^{[11]}
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 even though this equation seems at first look to emulate the Heisenberg equation of motion:
where A^{(H)} is some 'Heisenberg' operator. Taking the density operator to be in the Schrödinger picture makes sense since it is composed of 'Schrödinger' kets and bras evolved in time as per the Schrödinger picture. If the Hamiltonian is timeindependent, this differential equation can be solved to get
"Quantum Liouville", Moyal's equation
Under the Wigner map, the density matrix transforms into the Wigner function. The equation for timeevolution of the Wigner function is then the transform of the above von Neumann equation,
where H(q,p) is the Hamiltonian, and { { •,• } } is the Moyal bracket, the transform of the quantum commutator.
The evolution equation for the Wigner function is then analogous to that of its classical limit, the Liouville equation of classical physics. In the limit of vanishing Planck's constant ħ, W(q,p,t) reduces to the classical Liouville probability density function in phase space. The classical Liouville equation can be solved using the method of characteristics for partial differential equations, the characteristic equations being Hamilton's equations. The Moyal equation in quantum mechanics similarly admits formal solutions in terms of quantum characteristics, predicated on the ∗−product of phase space.
Composite Systems
The joint density matrix of a composite system of two systems A and B is described by ρ_{AB}. Then the subsystems are described by their reduced density operator.
is called partial trace over system B. If A and B are two distinct and independent systems then which is a product state.
C*algebraic formulation of states
It is now generally accepted that the description of quantum mechanics in which all selfadjoint operators represent observables is untenable.^{[12]}^{[13]} 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. However, by using the GNS construction, we can recover Hilbert spaces which realize A as a subalgebra 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
 ^ Fano, Ugo (1957), "Description of States in Quantum Mechanics by Density Matrix and Operator Techniques", Reviews of Modern Physics 29: 74–93, Bibcode 1957RvMP...29...74F, doi:10.1103/RevModPhys.29.74.
 ^ von Neumann, John (1927), "Wahrscheinlichkeitstheoretischer Aufbau der Quantenmechanik", Göttinger Nachrichten 1: 245–272.
 ^ Landau, L. D. (1927), "Das Dämpfungsproblem in der Wellenmechanik", Zeitschrift für Physik 45: 430–441, Bibcode 1927ZPhy...45..430L, doi:10.1007/BF01343064
 ^ Landau, L. D., and Lifshitz, E. M. (1977). Quantum Mechanics, NonRelativistic Theory: Volume 3. Oxford: Pergamon Press. pp. 41. ISBN 0080178014.
 ^ We assume is a pure state in this example.
 ^ Technically, F must be a Borel function
 ^ F(A) is defined to be the result of measuring the observable associated with the operator A, and then applying F to the outcome.
 ^ Nielsen, Michael; Chuang, Isaac (2000), Quantum Computation and Quantum Information, Cambridge University Press, ISBN 9780521635035. Chapter 11: Entropy and information, Theorem 11.9, "Projective measurements cannot decrease entropy"
 ^ Everett, Hugh (1973), "The Theory of the Universal Wavefunction (1956) Appendix I. "Monotone decrease of information for stochastic processes"", The ManyWorlds Interpretation of Quantum Mechanics, Princeton Series in Physics, Princeton University Press, pp. 128–129, ISBN 9780691081311
 ^ The theory of open quantum systems, by Breuer and Petruccione, p110.
 ^ Statistical mechanics, by Schwabl, p16.
 ^ See appendix, Mackey, George Whitelaw (1963), Mathematical Foundations of Quantum Mechanics, Dover Books on Mathematics, New York: Dover Publications, ISBN 9780486435176
 ^ Emch, Gerard G. (1972), Algebraic methods in statistical mechanics and quantum field theory, WileyInterscience, ISBN 9780471239000
 S. Kryszewski, [1]
See also
 Fluctuation theorem
 Density functional theory
 Dynamic nuclear polarisation
 Linear response function
 Green's function (manybody theory)
 Lindblad equation
 Green–Kubo relations
 Purification of quantum state
 POVM Generalized measurement of density states
 Quantum logic
 Wigner quasiprobability distribution
Wikimedia Foundation. 2010.