In ab initio molecular dynamics, whenever information about the potential energy surface is needed for integrating the equations of motion, it is computed “on the fly” using electronic structure calculations. For Born-Oppenheimer methods, the electronic structure calculations are converged, whereas in the extended Lagrangian approach the electronic structure is propagated along with the nuclei. Some recent advances for both approaches are discussed.