Solution (source code)

= Solution

If the <eigenvalues> of a <density operator> are $p_j\geq0$ with $\sum_jp_j=1$, then
$$
\gamma=\sum_jp_j^2\leq1,
$$
with equality exactly when one <eigenvalue> is one, namely for a pure state. In dimension $d$, Cauchy--Schwarz gives $\gamma\geq1/d$, with equality for the maximally mixed state. On an infinite-dimensional space the universal lower bound is zero as an infimum.