From the above discussion, we notice that the electronic density matrix ρ el can be diagonalized by unitary matrices, ρ d el = U † ρ el U, consisting of orthonormal eigenvectors, U = (v 1 , v 2 , ..., v N ). Thus the electronic density operator can be written in terms of its eigenvalues and eigenvectors, namely ρel =