Development of a new variational approach for thermal density matrices.