Burger Steven K, Liu Yuli, Sarkar Utpal, Ayers Paul W
Department of Chemistry, McMaster University, 1280 Main St. West, Hamilton, Ontario L8S 4M1, Canada.
J Chem Phys. 2009 Jan 14;130(2):024103. doi: 10.1063/1.2996579.
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.
通过使用谢泼德插值法,二次弦法(QSM)和快速行进法(FMM)所需的势能计算次数显著减少,该方法采用移动最小二乘法来拟合势能的高阶导数。势能的导数拟合至五阶。通过对插值值进行误差估计,这种移动最小二乘法增强的谢泼德插值方案大幅减少了FMM中的势能计算次数,通常减少幅度高达80%。对于所有网格间距,拟合至此处测试的最高阶(五阶)能得到最佳结果。对于QSM,使用增强的谢泼德插值法比使用通常的二阶近似、阻尼布罗伊登 - 弗莱彻 - 戈德法布 - 香农更新的海森矩阵来近似表面能得到稍好的结果。为了测试这些方法,我们研究了两种解析势能,丙氨酸二肽的旋转二面角势能以及氯甲烷与氟化物的S(N)2反应。