Hi,
as we know a density operator \rho is defined to be a non-negative definite operator of trace class (with trace 1).
We also know that for a given observable A, which is a (possibly unbounded) self-adjoint operator, the expectation value can be calculated as \operatorname{tr}(A \cdot...