Dear Paul, With the new quantum-espresso release, the code cp.x will be able to perform BO molecular dynamics, through conjugate gradient optimization of the electronic states (instead of Verlet). It works also for metals through ensemble DFT. Paolo Umari.