Finding Free-Energy Landmarks of Chemical Reactions.

We propose a novel approach to search for free-energy landmarks, i.e., minima and the saddle points, of chemical reactions in an automated manner using a combination of steepest descent and gentlest ascent methods. A numerical approach is suggested to improve the sampling efficiency of the second derivatives of the free-energy surface, which is required in the gentlest ascent method. This technique opens a way to identify free-energy landmarks of bond-breaking/creating processes in which the underlying potential energy surface is described using on-the-fly electronic structure calculations. As demonstrations of the approach, we present applications to the ring-opening of cis-1,2-dimethylbenzocyclobutene using the semiempirical PM7 method, focusing on the temperature dependence of the paths and barrier of the reaction, and we study an SN2 reaction in aqueous solution using a semiempirical QM/MM approach combining PM7 with the TIP3P water model.

[1]  E Weinan,et al.  The gentlest ascent dynamics , 2010, 1011.0042.

[2]  G. Henkelman,et al.  A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives , 1999 .

[3]  C. Chipot,et al.  Overcoming free energy barriers using unconstrained molecular dynamics simulations. , 2004, The Journal of chemical physics.

[4]  M. Tuckerman,et al.  Free Energy Reconstruction from Metadynamics or Adiabatic Free Energy Dynamics Simulations. , 2014, Journal of chemical theory and computation.

[5]  M. Tuckerman,et al.  Efficient and direct generation of multidimensional free energy surfaces via adiabatic dynamics without coordinate transformations. , 2008, The journal of physical chemistry. B.

[6]  Frédéric Legoll,et al.  Free energy calculations: an efficient adaptive biasing potential method. , 2009, The journal of physical chemistry. B.

[7]  E. Vanden-Eijnden,et al.  Single-sweep methods for free energy calculations. , 2007, The Journal of chemical physics.

[8]  B. Ruscic,et al.  W4 theory for computational thermochemistry: In pursuit of confident sub-kJ/mol predictions. , 2006, The Journal of chemical physics.

[9]  M. Tuckerman,et al.  An Adiabatic Molecular Dynamics Method for the Calculation of Free Energy Profiles , 2002 .

[10]  E. Vanden-Eijnden,et al.  A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations , 2006 .

[11]  Weinan E,et al.  Atomistic simulations of rare events using gentlest ascent dynamics. , 2011, The Journal of chemical physics.

[12]  M. Tuckerman,et al.  On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , 2002 .

[13]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[14]  Li Li,et al.  Bypassing the Kohn-Sham equations with machine learning , 2016, Nature Communications.

[15]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[16]  D. Truhlar,et al.  The DBH24/08 Database and Its Use to Assess Electronic Structure Model Chemistries for Chemical Reaction Barrier Heights. , 2009, Journal of chemical theory and computation.

[17]  S. Sakai.  Theoretical Studies of the Electrocyclic Reaction Mechanisms of o-Xylylene to Benzocyclobutene , 2000 .

[18]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

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

[20]  D. Marx,et al.  Understanding covalent mechanochemistry. , 2009, Angewandte Chemie.

[21]  H. Schaefer,et al.  Conrotatory and disrotatory stationary points for the electrocyclic isomerization of cyclobutene to cis-butadiene , 1984 .

[22]  Gérald Bernardinelli,et al.  Torquoselectivity in the electrocyclic conversion of benzocyclobutenes to o-xylylenes , 1992 .

[23]  Satoshi Maeda,et al.  Systematic exploration of the mechanism of chemical reactions: the global reaction route mapping (GRRM) strategy using the ADDF and AFIR methods. , 2013, Physical chemistry chemical physics : PCCP.

[24]  Mark E Tuckerman,et al.  Locating landmarks on high-dimensional free energy surfaces , 2015, Proceedings of the National Academy of Sciences.

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

[26]  M. Tuckerman,et al.  Heating and flooding: a unified approach for rapid generation of free energy surfaces. , 2012, The Journal of chemical physics.

[27]  Eric Darve,et al.  Adaptive biasing force method for scalar and vector free energy calculations. , 2008, The Journal of chemical physics.

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