Let $\rho$ a density operator $\rho: H \to H$ representing a mixed state, with $H$ Hilbert space. The operator $\rho$ (together with its corresponding matrix) represents a mixed state; when it additionally satisfies the property $\rho^2 = \rho$, then it is a pure state.
In any case, I was wondering if $\rho^2$ has any physical meaning at all. Clearly $\rho^2$ will still be self-adjoint because $\rho$ is by definition, but does it have any physical meaning at all? Would it make sense to compute the expectation value of $\rho^2$ as an observable?