Purification of quantum state

In quantum mechanics, especially quantum information, purification refers to the fact that every mixed state acting on finite dimensional Hilbert spaces can be viewed as the reduced state of some pure state.

In purely linear algebraic terms, it can be viewed as a statement about positive-semidefinite matrices.

Statement

Let ρ be a density matrix acting on a Hilbert space $H\_A$ of finite dimension "n". Then there exist a Hilbert space $H\_B$ and a pure state $|\; psi\; angle\; in\; H\_A\; otimes\; H\_B$ such that the partial trace of $|\; psi\; angle\; langle\; psi\; |$ with respect to $H\_B$

:$operatorname\{Tr\}\_B\; |\; psi\; angle\; langle\; psi\; |\; =\; ho.$

Proof

A density matrix is by definition positive semidefinite. So ρ has square root factorization $ho\; =\; A\; A^*\; =\; sum\_\{i\; =1\}\; ^n\; |\; i\; angle\; langle\; i\; |$. Let $H\_B$ be another copy of the "n"-dimensional Hilbert space with any orthonormal basis $\{\; |\; i\text{'}\; angle\; \}$. Define $|\; psi\; angle\; in\; H\_A\; otimes\; H\_B$ by

:$|\; psi\; angle\; =\; sum\_\{i\}\; |i\; angle\; otimes\; |\; i\text{'}\; angle.$

Direct calculation gives

:$operatorname\{Tr\}\_B\; |\; psi\; angle\; langle\; psi\; |\; =\; operatorname\{Tr\}\_B\; sum\_\{i,\; j\}\; |i\; angle\; langle\; j\; |\; otimes\; |\; i\text{'}\; angle\; langle\; j\text{'}|\; =\; ho.$

This proves the claim.

Note

* The vectorial pure state $|\; psi\; angle$ is in the form specified by the Schmidt decomposition.

* Since square root decompositions of a positive semidefinite matrix are not unique, neither are purifications.

* In linear algebraic terms, a square matrix is positive semidefinite if and only if it can be purified in the above sense. The "if" part of the implication follows immediately from the fact that the partial trace is a positive map.

An application: Stinespring's theorem

By combining Choi's theorem on completely positive maps and purification of a mixed state, we can recover the Stinespring dilation theorem for the finite dimensional case.