Development of a new variational approach for thermal density matrices

A McLachlan-type variational principle is derived for thermal density matrices. In this approach, the trace of the mean square of the differences between the derivatives of the exact and model density matrices is minimized with respect to the parameters in the model Hamiltonian. Applications to model anharmonic systems in the independent particle model show that the method can provide thermodynamic state functions accurately (within 5% of the converged basis set results) and at the same level of accuracy as the results using Feynman-Gibbs-Bogoliubov variational principle at this level of approximation.