Physical Review A
Representation of a quantum ensemble as a minimal set of pure states
Thomas E. Skinner, Steffen J. Glaser
Abstract
The density matrix $\ensuremath{\rho}$ for an n-level system is decomposed into the minimum number of pure states necessary to calculate physical observables. The corresponding physical system is first represented by a set $\mathcal{B}$ of n pure states $|{\ensuremath{\beta}}_{i}〉,$ together with their statistical weights. The time evolution of the system is therefore calculated as $\mathcal{B}(t)=U\mathcal{B}{(t}_{0}),$ with the propagator U operating on each member of the set, in contrast to the more laborious $\ensuremath{\rho}(t)=U\ensuremath{\rho}{(t}_{0}{)U}^{\ifmmode\dagger\else\textdagger\fi{}}.$ At least one of the states can be eliminated from the set by reducing its weight to zero via a simple transformation of the density matrix. When there are degenerate weights, the transformation is applied using the weight with the largest degeneracy. Thus, even if none of the original statistical weights are equal to zero, so that $\mathrm{rank}(\ensuremath{\rho})=n,$ the system can be described by a set of m states with $m<n.$ The resulting minimal set provides a concommitant increase in calculational efficiency that depends on the fractional degeneracy. The independent evolution of the states comprising the set also provides an inherent parallelism that can be readily utilized in numerical computations. In addition to efficiency and conceptual simplification, this decomposition can provide physical insight, since the explicit contribution of each pure state to the dynamics of the system is readily apparent.