Von Neumann entropy

In quantum statistical mechanics, the von Neumann entropy, named after John von Neumann, is the extension of classical Gibbs entropy concepts to the field of quantum mechanics. For a quantum-mechanical system described by a density matrix ρ, the von Neumann entropy is[1]

$S = - \mathrm{tr}(\rho \ln \rho),$

where tr denotes the trace and ln denotes the (natural) matrix logarithm. If ρ is written in terms of its eigenvectors |1〉, |2〉, |3〉, ... as

$\rho = \sum_j \eta_j \left| j \right\rang \left\lang j \right| ~,$

then the von Neumann entropy is merely[1]

$S = -\sum_j \eta_j \ln \eta_j.$

In this form, S can be seen to amount to the information theoretic Shannon entropy.[1][clarification needed]

Background

John von Neumann established a rigorous mathematical framework for quantum mechanics in his 1932 work Mathematical Foundations of Quantum Mechanics.[2] In it, he provided a theory of measurement, where the usual notion of wave-function collapse is described as an irreversible process (the so-called von Neumann or projective measurement).

The density matrix was introduced, with different motivations, by von Neumann and by Lev Landau. The motivation that inspired Landau was the impossibility of describing a subsystem of a composite quantum system by a state vector.[3] On the other hand, von Neumann introduced the density matrix in order to develop both quantum statistical mechanics and a theory of quantum measurements.

The density matrix formalism was developed to extend the tools of classical statistical mechanics to the quantum domain. In the classical framework we compute the partition function of the system in order to evaluate all possible thermodynamic quantities. Von Neumann introduced the density matrix in the context of states and operators in a Hilbert space. The knowledge of the statistical density matrix operator would allow us to compute all average quantities in a conceptually similar, but mathematically different way. Let us suppose we have a set of wave functions |Ψ〉 that depend parametrically on a set of quantum numbers n1, n2, ..., nN. The natural variable which we have is the amplitude with which a particular wavefunction of the basic set participates in the actual wavefunction of the system. Let us denote the square of this amplitude by p(n1, n2, ..., nN). The goal is to turn this quantity p into the classical density function in phase space. We have to verify that p goes over into the density function in the classical limit, and that it has ergodic properties. After checking that p(n1, n2, ..., nN) is a constant of motion, an ergodic assumption for the probabilities p(n1, n2, ..., nN) makes p a function of the energy only.

After this procedure, one finally arrives at the density matrix formalism when seeking a form where p(n1, n2, ..., nN) is invariant with respect to the representation used. In the form it is written, it will only yield the correct expectation values for quantities which are diagonal with respect to the quantum numbers n1, n2, ..., nN.

Expectation values of operators which are not diagonal involve the phases of the quantum amplitudes. Suppose we encode the quantum numbers n1, n2, ..., nN into the single index i or j. Then our wave function has the form

$\left| \Psi \right\rangle \,=\, \sum_i a_i\, \left| \psi_i \right\rangle .$

The expectation value of an operator B which is not diagonal in these wave functions, so

$\left\langle B \right\rangle \,=\, \sum_{i,j} a_i^{*}a_j\, \left\langle i \right| B \left| j \right\rangle .$

The role which was originally reserved for the quantities $\left| a_i \right| ^2$ is thus taken over by the density matrix of the system S.

$\left\langle j \right| \, \rho \, \left| i \right\rangle \,=\, a_j \, a_i^{*} .$

$\left\langle B \right\rangle \,=\, \mathrm{tr} (\rho \, B) ~.$

The invariance of the above term is described by matrix theory. A mathematical framework was described where the expectation value of quantum operators, as described by matrices, is obtained by taking the trace of the product of the density operator ρ̂ and an operator (Hilbert scalar product between operators). The matrix formalism here is in the statistical mechanics framework, although it applies as well for finite quantum systems, which is usually the case, where the state of the system cannot be described by a pure state, but as a statistical operator ρ̂ of the above form. Mathematically, ρ̂ is a positive-semidefinite hermitian matrix with unit trace.

Definition

Given the density matrix ρ, von Neumann defined the entropy[4][5] as

$S(\rho) \,=\,-\mathrm{tr} (\rho \, {\rm \ln} \rho),$

which is a proper extension of the Gibbs entropy (up to a factor kB) and the Shannon entropy to the quantum case. To compute S(ρ) it is convenient (see logarithm of a matrix) to compute the Eigendecomposition of $~\rho = \sum_j \eta_j \left| j \right\rangle \left\langle j \right|$. The von Neumann entropy is then given by

$S(\rho) \,=\, - \sum_j \eta_j \ln \eta_j ~.$

Since, for a pure state, the density matrix is idempotent, ρ = ρ2, the entropy S(ρ) for it vanishes. Thus, if the system is finite (finite-dimensional matrix representation), the entropy S(ρ) quantifies the departure of the system from a pure state. In other words, it codifies the degree of mixing of the state describing a given finite system. Measurement decoheres a quantum system into something noninterfering and ostensibly classical; so, e.g., the vanishing entropy of a pure state $\Psi = ( \left| 0 \right\rangle + \left| 1 \right\rangle ) / \sqrt{2}$, corresponding to a density matrix

$\rho = {1\over 2} \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}$

increases to S = ln 2 = 0.69 for the measurement outcome mixture

$\rho = {1\over 2} \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}$

as the quantum interference information is erased.

Properties

Some properties of the von Neumann entropy:

• S(ρ) is zero if and only if ρ represents a pure state.
• S(ρ) is invariant under changes in the basis of ρ, that is, S(ρ) = S(UρU), with U a unitary transformation.
• S(ρ) is concave, that is, given a collection of positive numbers λi which sum to unity ($\Sigma_i \lambda_i = 1$) and density operators ρi, we have
$S\bigg(\sum_{i=1}^k \lambda_i \, \rho_i \bigg) \,\geq\, \sum_{i=1}^k \lambda_i \, S(\rho_i).$
• S(ρ) is additive for independent systems. Given two density matrices ρA , ρB describing independent systems A and B, we have
$S(\rho_A \otimes \rho_B)=S(\rho_A)+S(\rho_B)$.
• S(ρ) is strongly subadditive for any three systems A, B, and C:
$S(\rho_{ABC}) + S(\rho_{B}) \leq S(\rho_{AB}) + S(\rho_{BC}).$.
This automatically means that S(ρ) is subadditive:
$S(\rho_{AC}) \leq S(\rho_{A}) +S(\rho_{C}).$

Below, the concept of subadditivity is discussed, followed by its generalization to strong subadditivity.

If ρA, ρB are the reduced density matrices of the general state ρAB, then

$\left| S(\rho_A)\,-\,S(\rho_B) \right| \,\leq \, S(\rho_{AB}) \, \leq \, S(\rho_A)\,+\,S(\rho_B) ~.$

This right hand inequality is known as subadditivity. The two inequalities together are sometimes known as the triangle inequality. They were proved in 1970 by Huzihiro Araki and Elliott H. Lieb.[6] While in Shannon's theory the entropy of a composite system can never be lower than the entropy of any of its parts, in quantum theory this is not the case, i.e., it is possible that S(ρAB) = 0, while S(ρA) = S(ρB) > 0.

Intuitively, this can be understood as follows: In quantum mechanics, the entropy of the joint system can be less than the sum of the entropy of its components because the components may be entangled. For instance, as seen explicitly, the Bell state of two spin-½s,

$\left| \psi \right\rangle = \left| \uparrow \downarrow \right\rangle + \left| \downarrow \uparrow \right\rangle ,$

is a pure state with zero entropy, but each spin has maximum entropy when considered individually in its reduced density matrix.[7] The entropy in one spin can be "cancelled" by being correlated with the entropy of the other. The left-hand inequality can be roughly interpreted as saying that entropy can only be cancelled by an equal amount of entropy.

If system A and system B have different amounts of entropy, the lesser can only partially cancel the greater, and some entropy must be left over. Likewise, the right-hand inequality can be interpreted as saying that the entropy of a composite system is maximized when its components are uncorrelated, in which case the total entropy is just a sum of the sub-entropies. This may be more intuitive in the phase space formulation, instead of Hilbert space one, where the Von Neumann entropy amounts to minus the expected value of the -logarithm of the Wigner function up to an offset shift.[5] Up to this normalization offset shift, the entropy is majorized by that of its classical limit.

The von Neumann entropy is also strongly subadditive. Given three Hilbert spaces, A, B, C,

$S(\rho_{ABC}) \, + \, S(\rho_{B}) \, \leq \, S(\rho_{AB}) \,+\, S(\rho_{BC}).$

This is a more difficult theorem and was proved in 1973 by Elliott H. Lieb and Mary Beth Ruskai,[8] using a matrix inequality of Elliott H. Lieb[9] proved in 1973. By using the proof technique that establishes the left side of the triangle inequality above, one can show that the strong subadditivity inequality is equivalent to the following inequality.

$S(\rho_{A}) \, + \, S(\rho_{C}) \, \leq \, S(\rho_{AB}) \,+\, S(\rho_{BC})$

when ρAB, etc. are the reduced density matrices of a density matrix ρABC. If we apply ordinary subadditivity to the left side of this inequality, and consider all permutations of A, B, C, we obtain the triangle inequality for ρABC: Each of the three numbers S(ρAB), S(ρBC), S(ρAC) is less than or equal to the sum of the other two.

Uses

The von Neumann entropy is being extensively used in different forms (conditional entropies, relative entropies, etc.) in the framework of quantum information theory.[10] Entanglement measures are based upon some quantity directly related to the von Neumann entropy. However, there have appeared in the literature several papers dealing with the possible inadequacy of the Shannon information measure, and consequently of the von Neumann entropy as an appropriate quantum generalization of Shannon entropy. The main argument is that in classical measurement the Shannon information measure is a natural measure of our ignorance about the properties of a system, whose existence is independent of measurement.

Conversely, quantum measurement cannot be claimed to reveal the properties of a system that existed before the measurement was made.[11] This controversy has encouraged some authors to introduce the non-additivity property of Tsallis entropy (a generalization of the standard Boltzmann–Gibbs entropy) as the main reason for recovering a true quantal information measure in the quantum context, claiming that non-local correlations ought to be described because of the particularity of Tsallis entropy.

References

1. ^ a b c Bengtsson, Ingemar; Zyczkowski, Karol. Geometry of Quantum States: An Introduction to Quantum Entanglement (1st ed.). p. 301.
2. ^ Von Neumann, John (1932). Mathematische Grundlagen der Quantenmechanik. Berlin: Springer. ISBN 3-540-59207-5.; Von Neumann, John (1955). Mathematical Foundations of Quantum Mechanics. Princeton University Press. ISBN 978-0-691-02893-4.
3. ^ Landau, L. (1927). "Das Daempfungsproblem in der Wellenmechanik". Zeitschrift fuer Physik 45 (5–6): 430–464. doi:10.1007/BF01343064.
4. ^ Geometry of Quantum States: An Introduction to Quantum Entanglement, by Ingemar Bengtsson, Karol Życzkowski, p301
5. ^ a b Zachos, C. K. (2007). "A classical bound on quantum entropy". Journal of Physics A: Mathematical and Theoretical 40 (21): F407. arXiv:hep-th/0609148. Bibcode:2007JPhA...40..407Z. doi:10.1088/1751-8113/40/21/F02.
6. ^ Huzihiro Araki and Elliott H. Lieb, Entropy Inequalities, Communications in Mathematical Physics, vol 18, 160–170 (1970).
7. ^ Zurek, W. H. (2003). "Decoherence, einselection, and the quantum origins of the classical". Reviews of Modern Physics 75 (3): 715. arXiv:quant-ph/0105127. Bibcode:2003RvMP...75..715Z. doi:10.1103/RevModPhys.75.715.
8. ^ Elliott H. Lieb and Mary Beth Ruskai, Proof of the Strong Subadditivity of Quantum-Mechanical Entropy, Journal of Mathematical Physics, vol 14, 1938–1941 (1973).
9. ^ Elliott H. Lieb, Convex Trace Functions and the Wigner–Yanase–Dyson Conjecture, Advances in Mathematics, vol 67, 267–288 (1973).
10. ^ Nielsen, Michael A. and Isaac Chuang (2001). Quantum computation and quantum information (Repr. ed.). Cambridge [u.a.]: Cambridge Univ. Press. p. 700. ISBN 978-0-521-63503-5.
11. ^ Pluch, P. (2006). Theory for Quantum Probability, PhD Thesis, Klagenfurt University.