Incorporating Linear Synchronous Transit Interpolation into the Growing String Method: Algorithm and Applications

被引:39
作者
Behn, Andrew [1 ]
Zimmerman, Paul M. [2 ]
Bell, Alexis T. [1 ]
Head-Gordon, Martin [2 ]
机构
[1] Univ Calif Berkeley, Dept Chem & Biomol Engn, Berkeley, CA 94720 USA
[2] Univ Calif Berkeley, Dept Chem, Berkeley, CA 94720 USA
基金
美国国家科学基金会;
关键词
NUDGED ELASTIC BAND; DETERMINING REACTION PATHS; POTENTIAL-ENERGY; QUANTUM; POINTS; STATE;
D O I
10.1021/ct200654u
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
The growing string method is a powerful tool in the systematic study of chemical reactions with theoretical methods which allows for the rapid identification of transition states connecting known reactant and product structures. However, the efficiency of this method is heavily influenced by the choice of interpolation scheme when adding new nodes to the string during optimization. In particular, the use of Cartesian coordinates with cubic spline interpolation often produces guess structures which are far from the final reaction path and require many optimization steps (and thus many energy and gradient calculations) to yield a reasonable final structure. In this paper, we present a new method for interpolating and reparameterizing nodes within the growing string method using the linear synchronous transit method of Halgren and Lipscomb. When applied to the alanine dipeptide rearrangement and a simplified cationic alkyl ring condensation reaction, a significant speedup in terms of computational cost is achieved (30-50%).
引用
收藏
页码:4019 / 4025
页数:7
相关论文
共 41 条
[1]   Implementation of an Algorithm Based on the Runge-Kutta-Fehlberg Technique and the Potential Energy as a Reaction Coordinate to Locate Intrinsic Reaction Paths [J].
Aguilar-Mogas, Antoni ;
Gimenez, Xavier ;
Bofill, Josep Maria .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 2010, 31 (13) :2510-2525
[2]   Finding reaction paths using the potential energy as reaction coordinate [J].
Aguilar-Mogas, Antoni ;
Gimenez, Xavier ;
Maria Bofill, Josep .
JOURNAL OF CHEMICAL PHYSICS, 2008, 128 (10)
[3]  
[Anonymous], 2007, J CHEM PHYS
[4]  
[Anonymous], 2002, PHYS REV B
[5]   A combined method for determining reaction paths, minima, and transition state geometries [J].
Ayala, PY ;
Schlegel, HB .
JOURNAL OF CHEMICAL PHYSICS, 1997, 107 (02) :375-384
[6]   SEARCH FOR STATIONARY-POINTS ON SURFACE [J].
BANERJEE, A ;
ADAMS, N ;
SIMONS, J ;
SHEPARD, R .
JOURNAL OF PHYSICAL CHEMISTRY, 1985, 89 (01) :52-57
[7]   Quadratic string method for determining the minimum-energy path based on multiobjective optimization [J].
Burger, SK ;
Yang, WT .
JOURNAL OF CHEMICAL PHYSICS, 2006, 124 (05)
[8]   Methods for finding transition states on reduced potential energy surfaces [J].
Burger, Steven K. ;
Ayers, Paul W. .
JOURNAL OF CHEMICAL PHYSICS, 2010, 132 (23)
[9]   Sequential quadratic programming method for determining the minimum energy path [J].
Burger, Steven K. ;
Yang, Weitao .
JOURNAL OF CHEMICAL PHYSICS, 2007, 127 (16)
[10]  
Cerjan C.J., 1981, J. Chem. Phys, V112, P2129