  Ab initio molecular dynamics (AIMD) simulations using hybrid density
functionals and plane waves are of great interest owing to the accuracy of this
approach in treating condensed matter systems. On the other hand, such AIMD
calculations are not routinely carried out since the computational cost
involved in applying the Hartree Fock exchange operator is very high. In this
work, we make use of a strategy that combines adaptively compressed exchange
operator formulation and multiple time step integration to significantly reduce
the computational cost of these simulations. We demonstrate the efficiency of
this approach for a realistic condensed matter system.
