Moving least-squares enhanced Shepard interpolation for the fast marching and string methods.

The number of the potential energy calculations required by the quadratic string method (QSM), and the fast marching method (FMM) is significantly reduced by using Shepard interpolation, with a moving least squares to fit the higher-order derivatives of the potential. The derivatives of the potential are fitted up to fifth order. With an error estimate for the interpolated values, this moving least squares enhanced Shepard interpolation scheme drastically reduces the number of potential energy calculations in FMM, often by up 80%. Fitting up through the highest order tested here (fifth order) gave the best results for all grid spacings. For QSM, using enhanced Shepard interpolation gave slightly better results than using the usual second order approximate, damped Broyden-Fletcher-Goldfarb-Shanno updated Hessian to approximate the surface. To test these methods we examined two analytic potentials, the rotational dihedral potential of alanine dipeptide and the S(N)2 reaction of methyl chloride with fluoride.

[1]  A Fast Method to Simulate Travelling Waves in Nonhomogeneous Chemical or Biological Media , 2001, nlin/0112020.

[2]  M. A. Collins,et al.  The utility of higher order derivatives in constructing molecular potential energy surfaces by interpolation , 1995 .

[3]  M. A. Collins,et al.  Molecular potential energy surfaces by interpolation , 1994 .

[4]  D. Wales,et al.  A doubly nudged elastic band method for finding transition states. , 2004, The Journal of chemical physics.

[5]  Michael A. Collins,et al.  Molecular potential-energy surfaces for chemical reaction dynamics , 2002 .

[6]  Weitao Yang,et al.  Ab initio QM/MM study shows there is no general acid in the reaction catalyzed by 4-oxalocrotonate tautomerase. , 2003, Journal of the American Chemical Society.

[7]  Michael A. Collins,et al.  Learning to interpolate molecular potential energy surfaces with confidence: A Bayesian approach , 1999 .

[8]  David J Wales,et al.  Folding of the GB1 hairpin peptide from discrete path sampling. , 2004, The Journal of chemical physics.

[9]  W. E,et al.  Finite temperature string method for the study of rare events. , 2002, Journal of Physical Chemistry B.

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

[11]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.

[12]  George C. Schatz,et al.  The analytical representation of electronic potential-energy surfaces , 1989 .

[13]  H. Huber,et al.  Non-Hamiltonian molecular dynamics implementation of the Gibbs ensemble method. I. Algorithm. , 2007, The Journal of chemical physics.

[14]  Jerry M. Parks,et al.  Quantum mechanics/molecular mechanics minimum free-energy path for accurate reaction energetics in solution and enzymes: sequential sampling and optimization on the potential of mean force surface. , 2008, The Journal of chemical physics.

[15]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[16]  D. Truhlar.  Potential Energy Surfaces and Dynamics Calculations , 1981 .

[17]  Paul W. Ayers,et al.  Computing tunneling paths with the Hamilton–Jacobi equation and the fast marching method , 2007 .

[18]  David J Wales,et al.  Effect of salt bridges on the energy landscape of a model protein. , 2004, The Journal of chemical physics.

[19]  Richard Dawes,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: computing high-density potential energy surface data from low-density ab initio data points. , 2007, The Journal of chemical physics.

[20]  H. Jónsson,et al.  Nudged elastic band method for finding minimum energy paths of transitions , 1998 .

[21]  Steven K. Burger,et al.  Quadratic string method for determining the minimum-energy path based on multiobjective optimization. , 2006, The Journal of chemical physics.

[22]  George C. Schatz,et al.  A local interpolation scheme using no derivatives in quantum-chemical calculations , 1999 .

[23]  George C. Schatz,et al.  A local interpolation scheme using no derivatives in potential sampling: Application to O(1D) + H2 system , 2003, J. Comput. Chem..

[24]  William H. Miller.  Dynamics of Molecular Collisions , 1976 .

[25]  E. Rouy,et al.  A viscosity solutions approach to shape-from-shading , 1992 .

[26]  Bernard R. Brooks,et al.  Exploring the quantum mechanical/molecular mechanical replica path method: a pathway optimization of the chorismate to prephenate Claisen rearrangement catalyzed by chorismate mutase , 2003 .

[27]  Steven K. Burger,et al.  Sequential quadratic programming method for determining the minimum energy path. , 2007, The Journal of chemical physics.

[28]  A CFL-like constraint for the fast marching method in inhomogeneous chemical kinetics , 2008 .

[29]  Paul W. Ayers,et al.  Fast Marching Method for Calculating Reactive Trajectories for Chemical Reactions , 2007 .

[30]  G. Ciccotti,et al.  String method in collective variables: minimum free energy paths and isocommittor surfaces. , 2006, The Journal of chemical physics.

[31]  P. Pechukas.  Statistical Approximations in Collision Theory , 1976 .

[32]  Reinhard Farwig,et al.  Rate of convergence of Shepard's global interpolation formula , 1986 .

[33]  H. Schlegel,et al.  A combined method for determining reaction paths, minima, and transition state geometries , 1997 .

[34]  Weitao Yang,et al.  Reaction path determination for quantum mechanical/molecular mechanical modeling of enzyme reactions by combining first order and second order "chain-of-replicas" methods. , 2005, The Journal of chemical physics.

[35]  Wolfgang Quapp,et al.  A growing string method for the reaction pathway defined by a Newton trajectory. , 2005, The Journal of chemical physics.

[36]  H. Werner,et al.  Local interpolation of ab initio potential energy surfaces for direct dynamics studies of chemical reactions , 1999 .

[37]  A. Chakraborty,et al.  A growing string method for determining transition states: comparison to the nudged elastic band and string methods. , 2004, The Journal of chemical physics.

[38]  Pathways for conformational change in nitrogen regulatory protein C from discrete path sampling. , 2008, The journal of physical chemistry. B.

[39]  David J. Wales,et al.  Global optimization and folding pathways of selected α-helical proteins , 2005 .

[40]  Paul W Ayers,et al.  Hamilton-Jacobi equation for the least-action/least-time dynamical path based on fast marching method. , 2004, The Journal of chemical physics.

[41]  Weitao Yang,et al.  QM/MM Minimum Free Energy Path: Methodology and Application to Triosephosphate Isomerase. , 2007, Journal of chemical theory and computation.

[42]  P. Ayers,et al.  A Hamilton–Jacobi type equation for computing minimum potential energy paths , 2006 .

[43]  M. Tuckerman,et al.  IN CLASSICAL AND QUANTUM DYNAMICS IN CONDENSED PHASE SIMULATIONS , 1998 .

[44]  Akio Kawano,et al.  Improving the accuracy of interpolated potential energy surfaces by using an analytical zeroth-order potential function. , 2004, The Journal of chemical physics.