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 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 [, (H) over cap] = 0 is always satisfied, guaranteeing a very efficient and robust variational minimization algorithm, that can also be extended to nonconventional entropic formulations.
Alfredo Pasquarello, Chunmin Zhang, Assil Bouzid
Giovanni Pizzi, Antimo Marrazzo, Junfeng Qiao, Miki Bonacci, Andrea Ferretti