Journal
JOURNAL OF CHEMICAL PHYSICS
Volume 113, Issue 16, Pages 6566-6572Publisher
AMER INST PHYSICS
DOI: 10.1063/1.1308551
Keywords
-
Ask authors/readers for more resources
We present an improved version of our recent algorithm [B. Paizs, G. Fogarasi, and P. Pulay, J. Chem. Phys. 109, 6571 (1998)] for optimizing the geometries of large molecules. The approximate Cholesky factorization technique has been generalized to the case of redundant coordinates, and an alternative approach involving use of the B daggerB matrix in the iterative coordinate back transformation is described. The generalized full Cholesky factors of B daggerB are very sparse and the corresponding force and geometry transformations are fast and numerically stable, permitting us to apply this technique for internal coordinate geometry optimization of molecules containing thousands of atoms. As an example we present optimization data on alpha-helical alanine polypeptides, and various globular proteins. Results for the alanine polypeptides indicates that internal coordinate optimization is clearly superior to the first-order Cartesian optimization techniques generally used in force field calculations. The largest system investigated is alpha-helical Ac-(Ala)(999)-NH2 containing 9999 atoms, which was successfully optimized using less than a megaword of memory. Optimization of various globular proteins shows that our procedure can easily deal with highly redundant (including full primitive) coordinate sets. (C) 2000 American Institute of Physics. [S0021-9606(00)30839-X].
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available