Efficient Computation of Free Energy Surfaces of Chemical Reactions using Ab Initio Molecular Dynamics with Hybrid Functionals and Plane Waves

Abstract: Ab initio molecular dynamics (AIMD) simulations employing density functional theory (DFT) and plane waves are routinely carried out using density functionals at the level of Generalized Gradient Approximation (GGA). AIMD simulations employing hybrid density functionals are of great interest as it offers more accurate description of structural and dynamic properties than the GGA functionals. However, computational cost for carrying out calculations using hybrid functionals and plane wave basis set is at least two order of magnitude higher than GGA functionals. Recently, we proposed a strategy that combined the adaptively compressed exchange operator formulation and the multiple time step integration scheme to reduce the computational cost about an order of magnitude [J. Chem. Phys. 151, 151102 (2019)]. In this work, we demonstrate the application of this method to study chemical reactions, in particular, formamide hydrolysis in alkaline aqueous medium. By actuating our implementation with the well-sliced metadynamics scheme, we are able to compute the two-dimensional free energy surface of this reaction at the level of hybrid-DFT. Accuracy of PBE0 (hybrid) and PBE (GGA) functionals in predicting the free energetics of the chemical reaction is investigated here.