La matriz de densidad (operador de densidad) en mecánica cuántica se define como $$\hat{\rho} = \sum_{i} p_i |\psi_i\rangle \langle \psi_i|\, ,$$ donde el $|\psi_i\rangle$ son un sistema ortonormal completo y el $p_i$ son probabilidades, por lo que $$p_i \geq 0, \quad \sum_i p_i =1 \, . $$
¿Cuándo la matriz de densidad no es diagonal? Si expresamos la matriz en la base $|\psi_i\rangle$ ¿No debería ser siempre diagonal?