Ensemble density-functional theory for ab-initio molecular dynamics of metals and finite-temperature insulators

A new method is presented for performing first-principles molecular-dynamics simulations of systems with variable occupancies. We adopt a matrix representation for the one-particle statistical operator Gamma, to introduce a ``projected'' free energy functional G that depends on the Kohn-Sham orbitals only and that is invariant under their unitary transformations. The Liouville equation [ Gamma , H ] = 0 is always satisfied, guaranteeing a very efficient and stable variational minimization algorithm that can be extended to non-conventional entropic formulations or fictitious thermal distributions.