Methods for geometry optimization of large molecules.: I.: An O(N2) algorithm for solving systems of linear equations for the transformation of coordinates and forces

被引:58
作者
Farkas, Ö [1 ]
Schlegel, HB
机构
[1] Wayne State Univ, Dept Chem, Detroit, MI 48202 USA
[2] Eotvos Lorand Univ, Dept Organ Chem, H-1518 Budapest, Hungary
关键词
D O I
10.1063/1.477393
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
The most recent methods in quantum chemical geometry optimization use the computed energy and its first derivatives with an approximate second derivative matrix. The performance of the optimization process depends highly on the choice of the coordinate system. In most cases the optimization is carried out in a complete internal coordinate system using the derivatives computed with respect to Cartesian coordinates. The computational bottlenecks for this process are the transformation of the derivatives into the internal coordinate system, the transformation of the resulting step back to Cartesian coordinates, and the evaluation of the Newton-Raphson or rational function optimization (RFO) step. The corresponding systems of linear equations occur as sequences of the form y(i)=M(i)x(i), where M-i can be regarded as a perturbation of the previous symmetric matrix Mi-1. They are normally solved via diagonalization of symmetric real matrices requiring O(N-3) operations. The current study is focused on a special approach to solving these sequential systems of linear equations using a method based on the update of the inverse of the symmetric matrix Mi. For convergence, this algorithm requires a number of O(N-2) operations with an O(N-3) factor for only the first calculation. The method is generalized to include redundant (singular) systems. The application of the algorithm to coordinate transformations in large molecular geometry optimization is discussed. (C) 1998 American Institute of Physics. [S0021-9606(98)30341-4].
引用
收藏
页码:7100 / 7104
页数:5
相关论文
共 37 条
[1]   SEARCH FOR STATIONARY-POINTS ON SURFACE [J].
BANERJEE, A ;
ADAMS, N ;
SIMONS, J ;
SHEPARD, R .
JOURNAL OF PHYSICAL CHEMISTRY, 1985, 89 (01) :52-57
[2]  
Broyden C.G., 1970, J I MATH ITS APPL, V6, P76, DOI DOI 10.1093/IMAMAT/6.1.76
[3]   Analytic energy gradients for the Gaussian very fast multipole method (GvFMM) [J].
Burant, JC ;
Strain, MC ;
Scuseria, GE ;
Frisch, MJ .
CHEMICAL PHYSICS LETTERS, 1996, 248 (1-2) :43-49
[4]  
Crippen G.M., 1988, DISTANCE GEOMETRY MO
[5]   Semiempirical methods with conjugate gradient density matrix search to replace diagonalization for molecular systems containing thousands of atoms [J].
Daniels, AD ;
Millam, JM ;
Scuseria, GE .
JOURNAL OF CHEMICAL PHYSICS, 1997, 107 (02) :425-431
[6]  
DAPPRICH S, IN PRESS THEOR CHEM
[7]  
Dennis, 1996, NUMERICAL METHODS UN
[8]  
FARKAS O, IN PRESS
[9]   THE CALCULATION OF ABINITIO MOLECULAR GEOMETRIES - EFFICIENT OPTIMIZATION BY NATURAL INTERNAL COORDINATES AND EMPIRICAL CORRECTION BY OFFSET FORCES [J].
FOGARASI, G ;
ZHOU, XF ;
TAYLOR, PW ;
PULAY, P .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1992, 114 (21) :8191-8201
[10]  
FRISCH MJ, 1988, GAUSSIAN 98 REVISION