Efficient computation of free energy surfaces of chemical reactions using ab initio molecular dynamics with hybrid functionals and plane waves

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 a more accurate description of structural and dynamic properties than the GGA functionals. However, the computational cost for carrying out calculations using hybrid functionals and plane wave basis set is at least two orders of magnitude higher than that using 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 by 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 an alkaline aqueous medium. By actuating our implementation with the well‐sliced metadynamics scheme, we can compute the two‐dimensional free energy surface of this reaction at the level of hybrid‐DFT. This work also investigates the accuracy of the PBE0 (hybrid) and the PBE (GGA) functionals in predicting the free energetics of this chemical reaction.

[1]  M. Tuckerman Ab initio molecular dynamics: basic concepts, current trends and novel applications , 2002 .

[2]  Tamar Schlick,et al.  Molecular Modeling and Simulation: An Interdisciplinary Guide , 2010 .

[3]  Jesús A. Izaguirre,et al.  Verlet-I/R-RESPA/Impulse is Limited by Nonlinear Instabilities , 2003, SIAM J. Sci. Comput..

[4]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. , 1991, Physical review. B, Condensed matter.

[5]  François Gygi,et al.  Performance and accuracy of recursive subspace bisection for hybrid DFT calculations in inhomogeneous systems. , 2015, Journal of chemical theory and computation.

[6]  Bernd Ensing,et al.  Formamide hydrolysis in alkaline aqueous solution: insight from Ab initio metadynamics calculations. , 2006, Angewandte Chemie.

[7]  P. Giannozzi,et al.  Fast hybrid density-functional computations using plane-wave basis sets , 2018, Electronic Structure.

[8]  Francois Gygi,et al.  Electronic Structure of Aqueous Sulfuric Acid from First-Principles Simulations with Hybrid Functionals. , 2014, The journal of physical chemistry letters.

[9]  Jií Kolafa,et al.  Time‐reversible always stable predictor–corrector method for molecular dynamics of polarizable molecules , 2004, J. Comput. Chem..

[10]  Donald G Truhlar,et al.  Density functional theory for transition metals and transition metal chemistry. , 2009, Physical chemistry chemical physics : PCCP.

[11]  K. Burke,et al.  Rationale for mixing exact exchange with density functional approximations , 1996 .

[12]  Thomas E Markland,et al.  Multiple time step integrators in ab initio molecular dynamics. , 2014, The Journal of chemical physics.

[13]  Joost VandeVondele,et al.  Ab initio molecular dynamics using hybrid density functionals. , 2008, The Journal of chemical physics.

[14]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for multiple time scales: Systems with long range forces , 1991 .

[15]  François Gygi,et al.  Efficient Computation of Hartree-Fock Exchange Using Recursive Subspace Bisection. , 2013, Journal of chemical theory and computation.

[16]  Lin Lin,et al.  Adaptively Compressed Exchange Operator. , 2016, Journal of chemical theory and computation.

[17]  M. Klein,et al.  Revisiting the free energy profile for the nucleophilic attack of hydroxide on formamide in aqueous solution , 2006 .

[18]  R. S. Brown,et al.  Aspects of the hydrolysis of formamide: revisitation of the water reaction and determination of the solvent deuterium kinetic isotope effect in base , 2002 .

[19]  Nisanth N. Nair,et al.  Sampling free energy surfaces as slices by combining umbrella sampling and metadynamics , 2015, J. Comput. Chem..

[20]  Teodora Todorova,et al.  Molecular dynamics simulation of liquid water: hybrid density functionals. , 2006, The journal of physical chemistry. B.

[21]  M. Cascella,et al.  Formamide Hydrolysis Investigated by Multiple-Steering ab Initio Molecular Dynamics , 2004 .

[22]  E. Carter,et al.  Time-Reversible Multiple Time Scale ab Initio Molecular Dynamics , 1993 .

[23]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[24]  Stefano de Gironcoli,et al.  Advanced capabilities for materials modelling with Quantum ESPRESSO , 2017, Journal of physics. Condensed matter : an Institute of Physics journal.

[25]  R. Distasio,et al.  Local structure analysis in ab initio liquid water , 2015, 1503.00020.

[26]  François Gygi,et al.  Compact representations of Kohn-Sham invariant subspaces. , 2007, Physical review letters.

[27]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[28]  M. Klein,et al.  Hydroxide diffuses slower than hydronium in water because its solvated structure inhibits correlated proton transfer , 2018, Nature Chemistry.

[29]  Weitao Yang,et al.  Insights into Current Limitations of Density Functional Theory , 2008, Science.

[30]  Sagarmoy Mandal,et al.  Enhanced sampling and free energy calculations with hybrid functionals and plane waves for chemical reactions. , 2018, The Journal of chemical physics.

[31]  R. Martin,et al.  Electronic Structure: Basic Theory and Practical Methods , 2004 .

[32]  Fabio Pietrucci,et al.  Formamide reaction network in gas phase and solution via a unified theoretical approach: Toward a reconciliation of different prebiotic scenarios , 2015, Proceedings of the National Academy of Sciences.

[33]  Giulia Galli,et al.  First Principles Simulations of the Infrared Spectrum of Liquid Water Using Hybrid Density Functionals. , 2011, Journal of chemical theory and computation.

[34]  Scott B. Baden,et al.  Parallel implementation of γ‐point pseudopotential plane‐wave DFT with exact exchange , 2011, J. Comput. Chem..

[35]  Ryan P Steele,et al.  Multiple-Time Step Ab Initio Molecular Dynamics Based on Two-Electron Integral Screening. , 2015, Journal of chemical theory and computation.

[36]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[37]  Sagarmoy Mandal,et al.  Speeding-up ab initio molecular dynamics with hybrid functionals using adaptively compressed exchange operator based multiple timestepping. , 2019, The Journal of chemical physics.

[38]  Francois Gygi,et al.  Structural and electronic properties of aqueous NaCl solutions from ab initio molecular dynamics simulations with hybrid density functionals , 2014 .

[39]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[40]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[41]  Efficient mixed-force first-principles molecular dynamics. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[43]  Ryan P Steele,et al.  Communication: Multiple-timestep ab initio molecular dynamics with electron correlation. , 2013, The Journal of chemical physics.

[44]  Weitao Yang,et al.  Challenges for density functional theory. , 2012, Chemical reviews.

[45]  F. Gygi,et al.  Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals. , 2015, The journal of physical chemistry letters.

[46]  B. Leimkuhler,et al.  Stochastic, resonance-free multiple time-step algorithm for molecular dynamics with very large time steps , 2013, 1307.1167.

[47]  A. Becke Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[48]  Annabella Selloni,et al.  Order-N implementation of exact exchange in extended insulating systems , 2008, 0812.1322.

[49]  A. Pasquarello,et al.  Structural, Dynamical, and Electronic Properties of Liquid Water: A Hybrid Functional Study. , 2016, The journal of physical chemistry. B.

[50]  Donald G. Truhlar,et al.  Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods , 2010 .

[51]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for condensed systems with multiple time scales , 1990 .

[52]  Mark E. Tuckerman,et al.  Molecular dynamics in systems with multiple time scales: Systems with stiff and soft degrees of freedom and with short and long range forces , 1991 .

[53]  Wolfram Koch,et al.  A Chemist's Guide to Density Functional Theory , 2000 .

[54]  U. Rothlisberger,et al.  A Versatile Multiple Time Step Scheme for Efficient ab Initio Molecular Dynamics Simulations. , 2018, Journal of chemical theory and computation.

[55]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[56]  Weitao Yang,et al.  Localization and delocalization errors in density functional theory and implications for band-gap prediction. , 2007, Physical review letters.

[57]  Richard M. Martin Electronic Structure: Frontmatter , 2004 .

[58]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[59]  Gregory A. Voth,et al.  Exact exchange in ab initio molecular dynamics: An efficient plane-wave based algorithm , 1998 .

[60]  Chao Yang,et al.  Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene. , 2017, Journal of chemical theory and computation.

[61]  Roberto Car,et al.  The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water. , 2014, The Journal of chemical physics.

[62]  G. Scuseria,et al.  Hybrid functionals based on a screened Coulomb potential , 2003 .

[63]  M. Tuckerman Statistical Mechanics: Theory and Molecular Simulation , 2010 .

[64]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[65]  Thomas M Henderson,et al.  Screened hybrid density functionals for solid-state chemistry and physics. , 2009, Physical chemistry chemical physics : PCCP.

[66]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for multiple time scales: Systems with disparate masses , 1991 .

[67]  M. Klein,et al.  A systematic study of chloride ion solvation in water using van der Waals inclusive hybrid density functional theory , 2015, 1503.07492.

[68]  M. Tuckerman,et al.  Ab initio molecular dynamics: concepts, recent developments, and future trends. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[69]  V. Barone,et al.  Toward reliable density functional methods without adjustable parameters: The PBE0 model , 1999 .

[70]  Ryan P Steele,et al.  Multiple-Timestep ab Initio Molecular Dynamics Using an Atomic Basis Set Partitioning. , 2015, The journal of physical chemistry. A.

[71]  T. Arias,et al.  Iterative minimization techniques for ab initio total energy calculations: molecular dynamics and co , 1992 .