# Von Neumann entropy

In physics, the von Neumann entropy, named after John von Neumann, is a measure of the statistical uncertainty within a description of a quantum system. It extends the concept of Gibbs entropy from classical statistical mechanics to quantum statistical mechanics, and it is the quantum counterpart of the Shannon entropy from classical information theory. For a quantum-mechanical system described by a density matrix ρ, the von Neumann entropy is $S=-\operatorname {tr} (\rho \ln \rho ),$ where $\operatorname {tr}$ denotes the trace and $\operatorname {ln}$ denotes the matrix version of the natural logarithm. If the density matrix ρ is written in a basis of its eigenvectors $|1\rangle ,|2\rangle ,|3\rangle ,\dots$ as $\rho =\sum _{j}\eta _{j}\left|j\right\rangle \left\langle j\right|,$ then the von Neumann entropy is merely $S=-\sum _{j}\eta _{j}\ln \eta _{j}.$ In this form, S can be seen as the Shannon entropy of the eigenvalues, reinterpreted as probabilities.
The von Neumann entropy and quantities based upon it are widely used in the study of quantum entanglement.

## Fundamentals

In quantum mechanics, probabilities for the outcomes of experiments made upon a system are calculated from the quantum state describing that system. Each physical system is associated with a vector space, or more specifically a Hilbert space. The dimension of the Hilbert space may be infinite, as it is for the space of square-integrable functions on a line, which is used to define the quantum physics of a continuous degree of freedom. Alternatively, the Hilbert space may be finite-dimensional, as occurs for spin degrees of freedom. A density operator, the mathematical representation of a quantum state, is a positive semi-definite, self-adjoint operator of trace one acting on the Hilbert space of the system. A density operator that is a rank-1 projection is known as a pure quantum state, and all quantum states that are not pure are designated mixed. Pure states are also known as wavefunctions. Assigning a pure state to a quantum system implies certainty about the outcome of some measurement on that system (i.e., $P(x)=1$ for some outcome $x$). The state space of a quantum system is the set of all states, pure and mixed, that can be assigned to it. For any system, the state space is a convex set: Any mixed state can be written as a convex combination of pure states, though not in a unique way. The von Neumann entropy quantifies the extent to which a state is mixed.
The prototypical example of a finite-dimensional Hilbert space is a qubit, a quantum system whose Hilbert space is 2-dimensional. An arbitrary state for a qubit can be written as a linear combination of the Pauli matrices, which provide a basis for $2\times 2$ self-adjoint matrices: $\rho ={\tfrac {1}{2}}\left(I+r_{x}\sigma _{x}+r_{y}\sigma _{y}+r_{z}\sigma _{z}\right),$ where the real numbers $(r_{x},r_{y},r_{z})$ are the coordinates of a point within the unit ball and $\sigma _{x}={\begin{pmatrix}0&1\\1&0\end{pmatrix}},\quad \sigma _{y}={\begin{pmatrix}0&-i\\i&0\end{pmatrix}},\quad \sigma _{z}={\begin{pmatrix}1&0\\0&-1\end{pmatrix}}.$ The von Neumann entropy vanishes when $\rho$ is a pure state, i.e., when the point $(r_{x},r_{y},r_{z})$ lies on the surface of the unit ball, and it attains its maximum value when $\rho$ is the maximally mixed state, which is given by $r_{x}=r_{y}=r_{z}=0$.

## Properties

Some properties of the von Neumann entropy:
* S(ρ) is zero if and only if ρ represents a pure state.
* S(ρ) is maximal and equal to $\ln N$ for a maximally mixed state, N being the dimension of the Hilbert space.
* 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. That is, 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.

### 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}).$
The right hand inequality is known as subadditivity, and the left is sometimes known as the triangle inequality. 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. This is expressed by saying that the Shannon entropy is monotonic but the von Neumann entropy is not. For example, take the Bell state of two spin-1/2 particles: $\left|\psi \right\rangle =\left|\uparrow \downarrow \right\rangle +\left|\downarrow \uparrow \right\rangle .$ This is a pure state with zero entropy, but each spin has maximum entropy when considered individually, because its reduced density matrix is the maximally mixed state. This indicates that it is an entangled state; the use of entropy as an entanglement measure is discussed further below.

### Strong subadditivity

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}).$ 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})$ where ρ$_{AB}$, etc. are the reduced density matrices of a density matrix ρ$_{ABC}$. If we apply ordinary subadditivity to the left side of this inequality, we then find $S(\rho _{AC})\leq S(\rho _{AB})+S(\rho _{BC}).$ By symmetry, for any tripartite state ρ$_{ABC}$, each of the three numbers S(ρ$_{AB}$), S(ρ$_{BC}$), S(ρ$_{AC}$) is less than or equal to the sum of the other two.

### Minimum Shannon entropy

Given a quantum state and a specification of a quantum measurement, we can calculate the probabilities for the different possible results of that measurement, and thus we can find the Shannon entropy of that probability distribution. A quantum measurement can be specified mathematically as a positive operator valued measure, or POVM. In the simplest case, a system with a finite-dimensional Hilbert space and measurement with a finite number of outcomes, a POVM is a set of positive semi-definite matrices $\{F_{i}\}$ on the Hilbert space that sum to the identity matrix, $\sum _{i=1}^{n}F_{i}=\operatorname {I} .$ The POVM element $F_{i}$ is associated with the measurement outcome $i$, such that the probability of obtaining it when making a measurement on the quantum state $\rho$ is given by ${\text{Prob}}(i)=\operatorname {tr} (\rho F_{i}).$ A POVM is rank-1 if all of the elements are proportional to rank-1 projection operators. The von Neumann entropy is the minimum achievable Shannon entropy, where the minimization is taken over all rank-1 POVMs.

### Holevo χ quantity

If ρ$_i$ are density operators and λ$_i$ is a collection of positive numbers which sum to unity ($\Sigma _{i}\lambda _{i}=1$), then $\rho =\sum _{i=1}^{k}\lambda _{i}\rho _{i}$ is a valid density operator, and the difference between its von Neumann entropy and the weighted average of the entropies of the ρ$_i$ is bounded by the Shannon entropy of the λ$_i$: $S{\bigg\sum _{i=1}^{k}\lambda _{i}\rho _{i}{\bigg )}-\sum _{i=1}^{k}\lambda _{i}S(\rho _{i})\leq -\sum _{i=1}^{k}\lambda _{i}\log \lambda _{i}.$ Equality is attained when the supports of the ρ$_i$ – the spaces spanned by their eigenvectors corresponding to nonzero eigenvalues – are orthogonal. The difference on the left-hand side of this inequality is known as the Holevo χ quantity and also appears in Holevo's theorem, an important result in quantum information theory.

## Change under time evolution
### Unitary

The time evolution of an isolated system is described by a unitary operator: $\rho \to U\rho U^{\dagger }.$ Unitary evolution takes pure states into pure states, and it leaves the von Neumann entropy unchanged. This follows from the fact that the entropy of $\rho$ is a function of the eigenvalues of $\rho$.

### Measurement

A measurement upon a quantum system will generally bring about a change of the quantum state of that system. Writing a POVM does not provide the complete information necessary to describe this state-change process. To remedy this, further information is specified by decomposing each POVM element into a product: $E_{i}=A_{i}^{\dagger }A_{i}.$ The Kraus operators $A_{i}$, named for Karl Kraus, provide a specification of the state-change process. They are not necessarily self-adjoint, but the products $A_{i}^{\dagger }A_{i}$ are. If upon performing the measurement the outcome $E_{i}$ is obtained, then the initial state $\rho$ is updated to $\rho \to \rho '={\frac {A_{i}\rho A_{i}^{\dagger }}{\mathrm {Prob} (i)}}={\frac {A_{i}\rho A_{i}^{\dagger }}{\operatorname {tr} (\rho E_{i})}}.$ An important special case is the Lüders rule, named for Gerhart Lüders. If the POVM elements are projection operators, then the Kraus operators can be taken to be the projectors themselves: $\rho \to \rho '={\frac {\Pi _{i}\rho \Pi _{i}}{\operatorname {tr} (\rho \Pi _{i})}}.$ If the initial state $\rho$ is pure, and the projectors $\Pi _{i}$ have rank 1, they can be written as projectors onto the vectors $|\psi \rangle$ and $|i\rangle$, respectively. The formula simplifies thus to $\rho =|\psi \rangle \langle \psi |\to \rho '={\frac {|i\rangle \langle i|\psi \rangle \langle \psi |i\rangle \langle i|}{|\langle i|\psi \rangle |^{2}}}=|i\rangle \langle i|.$ We can define a linear, trace-preserving, completely positive map, by summing over all the possible post-measurement states of a POVM without the normalisation: $\rho \to \sum _{i}A_{i}\rho A_{i}^{\dagger }.$ It is an example of a quantum channel, and can be interpreted as expressing how a quantum state changes if a measurement is performed but the result of that measurement is lost. Channels defined by projective measurements can never decrease the von Neumann entropy; they leave the entropy unchanged only if they do not change the density matrix. A quantum channel will increase or leave constant the von Neumann entropy of every input state if and only if the channel is unital, i.e., if it leaves fixed the maximally mixed state. An example of a channel that decreases the von Neumann entropy is the amplitude damping channel for a qubit, which sends all mixed states towards a pure state.

## Thermodynamic meaning

The quantum version of the canonical distribution, the Gibbs states, are found by maximizing the von Neumann entropy under the constraint that the expected value of the Hamiltonian is fixed. A Gibbs state is a density operator with the same eigenvectors as the Hamiltonian, and its eigenvalues are $\lambda _{i}={\frac {1}{Z}}\exp \left(-{\frac {E_{i}}{k_{B}T}}\right),$ where T is the temperature, $k_{B}$ is the Boltzmann constant, and Z is the partition function. The von Neumann entropy of a Gibbs state is, up to a factor $k_{B}$, the thermodynamic entropy.

## Generalizations and derived quantities
### Conditional entropy

Let $\rho _{AB}$ be a joint state for the bipartite quantum system AB. Then the conditional von Neumann entropy $S(A|B)$ is the difference between the entropy of $\rho _{AB}$ and the entropy of the marginal state for subsystem B alone: $S(A|B)=S(\rho _{AB})-S(\rho _{B}).$ This is bounded above by $S(\rho _{A})$. In other words, conditioning the description of subsystem A upon subsystem B cannot increase the entropy associated with A.
Quantum mutual information can be defined as the difference between the entropy of the joint state and the total entropy of the marginals: $S(A:B)=S(\rho _{A})+S(\rho _{B})-S(\rho _{AB}),$ which can also be expressed in terms of conditional entropy: $S(A:B)=S(A)-S(A|B)=S(B)-S(B|A).$

### Relative entropy

Let $\rho$ and $\sigma$ be two density operators in the same state space. The relative entropy is defined to be $S(\sigma |\rho )=\operatorname {tr} [\rho (\log \rho -\log \sigma )].$ The relative entropy is always greater than or equal to zero; it equals zero if and only if $\rho =\sigma$. Unlike the von Neumann entropy itself, the relative entropy is monotonic, in that it decreases (or remains constant) when part of a system is traced over: $S(\sigma _{A}|\rho _{A})\leq S(\sigma _{AB}|\rho _{AB}).$

### Entanglement measures

Just as energy is a resource that facilitates mechanical operations, entanglement is a resource that facilitates performing tasks that involve communication and computation. The mathematical definition of entanglement can be paraphrased as saying that maximal knowledge about the whole of a system does not imply maximal knowledge about the individual parts of that system. If the quantum state that describes a pair of particles is entangled, then the results of measurements upon one half of the pair can be strongly correlated with the results of measurements upon the other. However, entanglement is not the same as "correlation" as understood in classical probability theory and in daily life. Instead, entanglement can be thought of as potential correlation that can be used to generate actual correlation in an appropriate experiment. The state of a composite system is always expressible as a sum, or superposition, of products of states of local constituents; it is entangled if this sum cannot be written as a single product term. Entropy provides one tool that can be used to quantify entanglement. If the overall system is described by a pure state, the entropy of one subsystem can be used to measure its degree of entanglement with the other subsystems. For bipartite pure states, the von Neumann entropy of reduced states is the unique measure of entanglement in the sense that it is the only function on the family of states that satisfies certain axioms required of an entanglement measure. It is thus known as the entanglement entropy.
It is a classical result that the Shannon entropy achieves its maximum at, and only at, the uniform probability distribution {1/n, ..., 1/n}. Therefore, a bipartite pure state ρ ∈ H$_A$ ⊗ H$_B$ is said to be a maximally entangled state if the reduced state of each subsystem of ρ is the diagonal matrix ${\begin{pmatrix}{\frac {1}{n}}&&\\&\ddots &\\&&{\frac {1}{n}}\end{pmatrix}}.$
For mixed states, the reduced von Neumann entropy is not the only reasonable entanglement measure. Some of the other measures are also entropic in character. For example, the relative entropy of entanglement is given by minimizing the relative entropy between a given state $\rho$ and the set of nonentangled, or separable, states. The entanglement of formation is defined by minimizing, over all possible ways of writing of $\rho$ as a convex combination of pure states, the average entanglement entropy of those pure states. The squashed entanglement is based on the idea of extending a bipartite state $\rho _{AB}$ to a state describing a larger system, $\rho _{ABE}$, such that the partial trace of $\rho _{ABE}$ over E yields $\rho _{AB}$. One then finds the infimum of the quantity ${\frac {1}{2}}[S(\rho _{AE})+S(\rho _{BE})-S(\rho _{E})-S(\rho _{ABE})],$ over all possible choices of $\rho _{ABE}$.

### Quantum Rényi entropies

Just as the Shannon entropy function is one member of the broader family of classical Rényi entropies, so too can the von Neumann entropy be generalized to the quantum Rényi entropies: $S_{\alpha }(\rho )={\frac {1}{1-\alpha }}\ln[\operatorname {tr} \rho ^{\alpha }]={\frac {1}{1-\alpha }}\ln \sum _{i=1}^{N}\lambda _{i}^{\alpha }.$ In the limit that $\alpha \to 1$, this recovers the von Neumann entropy. The quantum Rényi entropies are all additive for product states, and for any $\alpha$, the Rényi entropy $S_{\alpha }$ vanishes for pure states and is maximized by the maximally mixed state. For any given state $\rho$, $S_{\alpha }(\rho )$ is a continuous, nonincreasing function of the parameter $\alpha$. A weak version of subadditivity can be proven: $S_{\alpha }(\rho _{A})-S_{0}(\rho _{B})\leq S_{\alpha }(\rho _{AB})\leq S_{\alpha }(\rho _{A})+S_{0}(\rho _{B}).$ Here, $S_{0}$ is the quantum version of the Hartley entropy, i.e., the logarithm of the rank of the density matrix.

## History

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. On the other hand, von Neumann introduced the density matrix in order to develop both quantum statistical mechanics and a theory of quantum measurements. He introduced the expression now known as von Neumann entropy by arguing that a probabilistic combination of pure states is analogous to a mixture of ideal gases. Von Neumann first published on the topic in 1927. His argument was built upon earlier work by Albert Einstein and Leo Szilard.
Max Delbrück and Gert Molière proved the concavity and subadditivity properties of the von Neumann entropy in 1936. Quantum relative entropy was introduced by Hisaharu Umegaki in 1962. The subadditivity and triangle inequalities were proved in 1970 by Huzihiro Araki and Elliott H. Lieb. Strong subadditivity is a more difficult theorem. It was conjectured by Oscar Lanford and Derek Robinson in 1968. Lieb and Mary Beth Ruskai proved the theorem in 1973, using a matrix inequality proved earlier by Lieb.