Semiempirical methods with conjugate gradient density matrix search to replace diagonalization for molecular systems containing thousands of atoms
J. Chem. Phys. 107, 425 (1997); doi:10.1063/1.474404
Issue Date: 8 July 1997
You are not logged in to this journal. Log in
Conventional semiempirical methods using diagonalization are not practical for calculations on molecular systems containing more than a few hundred atoms because of
(N3) time and
(N2) memory requirements, where N is the number of atoms. Currently, the time dominating step is diagonalization of the Fock matrix. This paper demonstrates how
(N3) diagonalization and
(N2) memory requirements are eliminated by using a conjugate gradient search for the density matrix with sparse matrix techniques. Our method makes high accuracy energy calculations on molecules containing thousands of atoms possible on the typical workstation. Benchmark examples are presented on polyglycine chains (20000 atoms), water clusters (up to 1800 atoms), and nucleic acids (up to 6304 atoms). ©1997 American Institute of Physics.
| History: | Received 13 January 1997; accepted 3 April 1997 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/107/425/1 |
REFERENCES (42)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- D. Bakowies and W. Thiel,
J. Am. Chem. Soc. 113, 3704 (1991 ). - D. Bakowies, M. Bühl, and W. Thiel,
J. Am. Chem. Soc. 117, 10 113 (1995 ). - J. P. Stewart,
Int. J. Quantum Chem. 58, 133 (1996 ). - T.-S. Lee, D. M. York, and W. Yang, J. Chem. Phys. 105, 2744 (1996).
- S. L. Dixon and K. M. Merz, Jr., J. Chem. Phys. 104, 6643 (1996).
- X. P. Li, R. W. Nunes, and D. Vanderbilt, Phys. Rev. B 47, 10 891 (1993);
- S.-Y. Qiu, C. Z. Wang, K. M. Ho, and C. T. Chan,
J. Phys. Condensed. Matter 6, 9153 (1994 ). - C. H. Xu and G. E. Scuseria,
Chem. Phys. Lett. 262, 219 (1996 ). - J. M. Millam and G. E. Scuseria, J. Chem. Phys. (in press).
- C. A. White, B. G. Johnson, P. M. W. Gill, and M. Head-Gordon,
Chem. Phys. Lett. 230, 8 (1994 ). - M. C. Strain, G. E. Scuseria, and M. J. Frisch,
Science 271, 51 (1996 ). - M. J. S. Dewar and W. Thiel,
J. Am. Chem. Soc. 99, 4499 (1977 ). - A. D. Becke, J. Chem. Phys. 88, 2547 (1988).
- R. E. Stratmann, G. E. Scuseria, and M. J. Frisch,
Chem. Phys. Lett. 257, 213 (1996 ). - L. P. Davis, R. M. Guidry, J. R. Williams, M. J. S. Dewar, and H. S. Rzepa,
J. Comp. Chem. 2, 433 (1981 ). - M. J. S. Dewar, M. L. McKee, and H. S. Rzepa,
J. Am. Chem. Soc. 100, 3607 (1978 ). - M. J. S. Dewar, E. G. Zoebisch, E. F. Hearly, and J. P. Stewart,
J. Am. Chem. Soc. 107, 3902 (1985 ). - M. J. S. Dewar and C. H. Reynolds, J. Comp. Chem. 2, 140 (1986).
- M. J. S. Dewar and W. Thiel,
J. Am. Chem. Soc. 99, 4899 (1977 ). - M. J. S. Dewar and H. S. Rzepa,
J. Am. Chem. Soc. 100, 777 (1978 ). - M. J. S. Dewar and M. L. McKee,
J. Comp. Chem. 4, 84 (1978 ). - M. J. S. Dewar and E. F. Healy,
J. Comp. Chem. 4, 542 (1983 ). - M. J. S. Dewar, G. L. Grady, and J. J. P. Stewart,
J. Am. Chem. Soc. 106, 6771 (1984 ). - M. S. Dewaret al.,
Organometallics 4, 1964 (1985 ). - J. J. P. Stewart,
J. Comp. Chem. 10, 209 (1989 ). - J. J. P. Stewart,
J. Comp. Chem. 10, 221 (1989 ). - R. C. Bingham, M. Dewar, and Lo,
J. Am. Chem. Soc. 97, 1285, 1294, 1302, 1307 (1975 ). - R. McWeeny,
Rev. Mod. Phys. 32, 335 (1960 ). - W. H. Press, S. A. Teukolsky, W. T. Vetterling, and W. P. Flannery, Numerical Recipes in Fortran, 2nd ed. (Cambridge University Press, Cambridge, 1992).
- S. Pissanetsky, Sparse Matrix Technology (Academic, London, 1984).
- The polyglycine chains are similar to those described in Refs. 4, 37, and 38. Coordinates for these molecules are available from us.
- The water clusters up to 150 water molecules are those used in Ref. 39. We obtained the coordinates from the authors of Ref. 39. Coordinates for the 300 and 600 water molecule cluster are available from us.
- GAUSSIAN 95, Development Version, Revision D.4, M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, M. C. Strain, J. C. Burant, R. E. Stratmann, G. A. Petersson, J. A. Montgomery, V. G. Zakrzewski, T. Keith, K. Raghavachari, M. A. Al-Laham, J. V. Ortiz, J. B. Foresman, J. Cioslowski, B. B. Stefanov, A. Nanayakkara, M. Challacombe, C. Y. Peng, P. Y. Ayala, W. Chen, M. W. Wong, J. L. Andres, E. S. Replogle, R. Gomperts, R. L. Martin, D. J. Fox, J. S. Binkley, D. J. Defrees, J. Baker, J. P. Stewart, C. Gonzalez, M. Head-Gordon, P. M. W. Gill, B. G. Johnson, and J. A. Pople, Gaussian, Inc., Pittsburgh Pennsylvania, 1996.
- J. Almöf, K. Faergi, Jr., and K. Korsell,
J. Comp. Chem. 3, 385 (1982 ). - P. Pulay,
J. Comp. Chem. 3, 556 (1982 ). - The 368 atom RNA molecule is from Ref. 40. The 651 atom DNA molecule is from Ref. 41 with sodium atoms added to balance the charge. The 1576 atom RNA molecule is from Ref. 42. Hydrogens were added with Xplor. All larger RNA molecules were made from replicating the 1576 atom RNA molecule.
- W. Yang, Phys. Rev. Lett. 66, 1438 (1991).
- W. Yang,
J. Mol. Struct. (Theochem) 255, 461 (1992 ). - E. Schwegler and M. Challacombe, J. Chem. Phys. 105, 2726 (1996).
- J. Hutter, P. Carloni, and M. Parrinello,
J. Am. Chem. Soc. 118, 8710 (1996 ). - J. D. Baleja, R. T. Pon, and B. D. Sykes,
Biochemistry 29, 4828 (1990 ). - S. E. Lietzke, V. E. Carperos, and C. E. Kundrot,
Acta Crystallogr. Sect. D 52, 687 (1996 ).








