Weber D, Bellinger D, Engels B
Julius-Maximilians-Universität Würzburg, Institut für Physikalische und Theoretische Chemie, Würzburg, Germany.
Julius-Maximilians-Universität Würzburg, Institut für Physikalische und Theoretische Chemie, Würzburg, Germany.
Methods Enzymol. 2016;578:145-67. doi: 10.1016/bs.mie.2016.05.030. Epub 2016 Jun 20.
We present new schemes to improve the convergence of an important global optimization problem and to determine reaction pathways (RPs) between identified minima. Those methods have been implemented into the CAST program (Conformational Analysis and Search Tool). The first part of this chapter shows how to improve convergence of the Monte Carlo with minimization (MCM, also known as Basin Hopping) method when applied to optimize water clusters or aqueous solvation shells using a simple model. Since the random movement on the potential energy surface (PES) is an integral part of MCM, we propose to employ a hydrogen bonding-based algorithm for its improvement. We show comparisons of the results obtained for random dihedral and for the proposed random, rigid-body water molecule movement, giving evidence that a specific adaption of the distortion process greatly improves the convergence of the method. The second part is about the determination of RPs in clusters between conformational arrangements and for reactions. Besides standard approaches like the nudged elastic band method, we want to focus on a new algorithm developed especially for global reaction path search called Pathopt. We started with argon clusters, a typical benchmark system, which possess a flat PES, then stepwise increase the magnitude and directionality of interactions. Therefore, we calculated pathways for a water cluster and characterize them by frequency calculations. Within our calculations, we were able to show that beneath local pathways also additional pathways can be found which possess additional features.
我们提出了新的方案,以改进一个重要的全局优化问题的收敛性,并确定已识别的极小值之间的反应路径(RP)。这些方法已被应用于CAST程序(构象分析与搜索工具)中。本章的第一部分展示了在使用简单模型优化水团簇或水合溶剂壳时,如何改进蒙特卡罗最小化方法(MCM,也称为盆地跳跃法)的收敛性。由于势能面上的随机移动是MCM的一个组成部分,我们建议采用一种基于氢键的算法来改进它。我们展示了随机二面角和所提出的随机刚体水分子移动所获得结果的比较,证明了对畸变过程的特定调整极大地提高了该方法的收敛性。第二部分是关于确定团簇中构象排列之间以及反应的反应路径。除了像推挤弹性带法这样的标准方法外,我们还将重点介绍一种专门为全局反应路径搜索开发的新算法——Pathopt。我们从氩团簇开始,这是一个典型的基准系统,其势能面较为平坦,然后逐步增加相互作用的强度和方向性。因此,我们计算了水团簇的路径,并通过频率计算对其进行表征。在我们的计算中,我们能够表明,除了局部路径外,还能找到具有其他特征的额外路径。