PROJECTION OPERATOR METHOD FOR GEOMETRY OPTIMIZATION WITH CONSTRAINTS

被引:32
作者
LU, DH
ZHAO, M
TRUHLAR, DG
机构
[1] UNIV MINNESOTA,DEPT CHEM,CHEM PHYS PROGRAM,MINNEAPOLIS,MN 55455
[2] UNIV MINNESOTA,MINNESOTA SUPERCOMP INST,MINNEAPOLIS,MN 55455
关键词
D O I
10.1002/jcc.540120311
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
A new approach is presented for performing geometry optimization for stationary points on potential energy hypersurfaces with equality constraints on the internal coordinates of a polyatomic system. The working equations are the same as for unconstrained Newton-Raphson optimization in Cartesian coordinates except that projection operators are applied to the gradient and Hessian to enforce the constraints. Two reactive systems with different kinds of constraints are treated as examples: OH + H-2 --> OH3 not-equal --> H2O + H with one constrained OH bond distance and CH3 + H-2 --> CH5 not-equal --> CH4 + H with one constrained H-C-H bond angle in the CH3 group or with one constrained bond distance and one simultaneously constrained bond angle. In each case we optimized all reactants and products as well as the saddle point, all subject to the constraints.
引用
收藏
页码:376 / 384
页数:9
相关论文
共 32 条
[1]  
Arfken G., 1985, MATH METHODS PHYS, V3, P516
[2]   LOCATING TRANSITION-STATES [J].
BELL, S ;
CRIGHTON, JS .
JOURNAL OF CHEMICAL PHYSICS, 1984, 80 (06) :2464-2475
[3]  
BINKLEY JS, 1981, GAUSSIAN 80 AB INITI, V13, P30
[4]  
CALIFANO S, 1976, VIBRATIONAL STATES, pCH4
[5]  
CIARLET PG, 1982, INTRO NUMERICAL LINE
[6]   GLOBAL CONVERGENCE OF A CLASS OF TRUST REGION ALGORITHMS FOR OPTIMIZATION WITH SIMPLE BOUNDS [J].
CONN, AR ;
GOULD, NIM ;
TOINT, PL .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1988, 25 (02) :433-460
[8]  
FLETCHER R, 1982, NONLINEAR OPTIMIZATI, P185
[9]  
Fletcher R., 1987, PRACTICAL METHODS OP
[10]   LOCAL CONVERGENCE OF SECANT METHODS FOR NONLINEAR CONSTRAINED OPTIMIZATION [J].
FONTECILLA, R .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1988, 25 (03) :692-712