What is the best alternative to diagonalization of the Hamiltonian in large scale semiempirical calculations?

Journal of Chemical Physics

Volumn 110, Issue 3, 1999, Pages 1321-1328

  Daniels, Andrew D ^a   Scuseria, Gustavo E ^a  

^a RICE UNIVERSITY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

EID: 0000725931     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.478008     Document Type: Article

References (38)

1. T-S. Lee, D. M. York, and W. Yang, J. Chem. Phys. 105, 2744 (1996).
2. S. L. Dixon and K. M. Merz, Jr., J. Chem. Phys. 104, 6643 (1996).
3. X. P. Li, R. W. Nunes, and D. Vanderbilt, Phys. Rev. B 47, 10 891 (1993); R. W. Nunes and D. Vanderbilt, ibid. 50, 17611 (1994).
4. X. P. Li, R. W. Nunes, and D. Vanderbilt, Phys. Rev. B 47, 10 891 (1993); R. W. Nunes and D. Vanderbilt, ibid. 50, 17611 (1994).
5. C. H. Xu and G. E. Scuseria, Chem. Phys. Lett. 262, 219 (1996).
6. A. D. Daniels, J. M. Millam, and G. E. Scuseria, J. Chem. Phys. 107, 425 (1997).
7. J. P. Stewart, P. Császár, and P. Pulay, J. Comput. Chem. 3, 227 (1982).
8. J. P. Stewart, Int. J. Quantum Chem. 58, 133 (1996).
9. S. Goedecker and L. Colombo, Phys. Rev. B 73, 122 (1994).
10. M. J. S. Dewar, E. G. Zoebisch, E. F. Hearly, and J. P. Stewart, J. Am. Chem. Soc. 107, 3902 (1985).
11. GAUSSIAN 95, Development Version (Revision G.3), M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, V. G. Zakrzewski, J. A. Montgomery, R. E. Stratmann, J. C. Burant, S. Dapprich, J. M. Millam, A. D. Daniels, K. N. Kudin, M. C. Strain, O. Farkas, G. A. Petersson, P. Y. Ayala, A. D. Rabuck, K. Raghavachari, Q. Cui, K. Morokuma, J. B. Foresman, J. Cioslowski, J. V. Ortiz, V. Barone, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, W. Chen, M. W. Wong, J. L. Andres, E. S. Replogle, R. Gomperts, R. L. Martin, D. J. Fox, T. Keith, M. A. Al-Laham, A. Nanayakkara, M. Challacombe, C. Y. Peng, P. M. W. Gill, B. Johnson, and J. Pople. Gaussian, Inc., Pittsburgh, PA, 1998.
12. M. J. S. Dewar and W. Thiel, J. Am. Chem. Soc. 99, 4899 (1977).
13. J. J. P. Stewart, J. Comput. Chem. 10, 209, 221 (1989).
14. R. C. Bingham, M. J. S. Dewar, and D. H. Lo, J. Am. Chem. Soc. 97, 1285,1294,1302,1307 (1975).
15. R. McWeeny, Rev. Mod. Phys. 32, 335 (1960).
16. During each SCF iteration between 2 to 4 CG-DMS iterations are needed to produce a sufficiently converged density matrix. During each CG-DMS step 5 matrix multipications are required to form the gradient and 3 to calculate the step size. 2 more matrix multiplications are required to perform each purification transformation. To calculate the form of the matrix to be used in the CG-DMS step, 2 matrix multiples are needed and 2 matrix multiples are needed for the DIIS step. Thus, a calculation with 2 CG-DMS steps per SCF iteration and 2 purification cycles per CG-DMS step requires a total of 28 matrix multiples per SCF cycle.
17. S. Pissanetsky, Sparse Matrix Technology (Academic, London, 1984).
18. ov, and 1 more to rotate the orbitals. (The rotation is similar to a matrix multiplication in CPU cost.) Every SCF iteration 1 more matrix multiplication is required to form estimates for the eigenvalues of the matrix, 1 to form ρ, and 2 for DIIS. Hence, a calculation taking 2 PD steps per SCF cycle requires 10 matrix multiplications per SCF cycle.
19. S. Goedecker and M. Teter, Phys. Rev. B 51, 9455 (1995).
20. B. N. Bartlett, The Symmetric Eigenvalue Problem (Prentice-Hall, Englewood Cliffs, NJ, 1980). R. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users Guide: Solution of Large Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods (SIAM Publications, Philadelphia, 1998); D. C. Sorensen, SIAM J. Matrix Anal. Appl. 13, 357 (1992).
21. B. N. Bartlett, The Symmetric Eigenvalue Problem (Prentice-Hall, Englewood Cliffs, NJ, 1980). R. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users Guide: Solution of Large Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods (SIAM Publications, Philadelphia, 1998); D. C. Sorensen, SIAM J. Matrix Anal. Appl. 13, 357 (1992).
22. B. N. Bartlett, The Symmetric Eigenvalue Problem (Prentice-Hall, Englewood Cliffs, NJ, 1980). R. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users Guide: Solution of Large Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods (SIAM Publications, Philadelphia, 1998); D. C. Sorensen, SIAM J. Matrix Anal. Appl. 13, 357 (1992).
23. K. R. Bates, A. D. Daniels, and G. E. Scuseria, J. Chem. Phys. 109, 3308 (1998).
24. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and W. P. Flannery, Numerical Recipes in Fortran, 2nd ed. (Cambridge University Press, Cambridge, 1992).
25. R. Baer and M. Head-Gordon, J. Chem. Phys. 107, 10003 (1997).
26. i. In addition, it requires 2 matrix multiplications to determine the form of the matrices, and 2 for DIIS.
27. A. H. R. Palser and D. E. Manolopoulos, Phys. Rev. B. 58, 12704 (1998).
28. PDM requires 2 matrix multiplications per purification cycle. It also requires 2 matrix multiplications to determine the form of the matrices, and 2 for DIIS. So a PDM step using 9 purification transformations requires 22 matrix multiplications per SCF cycle.
29. P. Pulay, J. Comput. Chem. 3, 556 (1982).
30. J. Hutter, P. Carloni, and M. Parrinello, J. Am. Chem. Soc. 118, 8710 (1996).
31. The polyglycine chains are similar to those described in Refs. 1, 29, and 30. Coordinates for these molecules are available from us.
32. W. Yang, Phys. Rev. Lett. 66, 1438 (1991).
33. W. Yang, J. Mol. Struct.: THEOCHEM 255, 461 (1992).
34. The water clusters up to 150 water molecules are those used in Ref. 32. We obtained the coordinates from the authors of Ref. 32. Coordinates for the larger water molecule clusters are available from us.
35. E. Schwegler and M. Challacombe, J. Chem. Phys. 105, 2726 (1996).
36. The 368 and 1576 atom RNA molecules are from Refs. 26 and 35, respectively. The 651 atom DNA molecule is from Ref. 34 with sodium atoms added to balance the charge. Hydrogens were added with XPLOR. All larger RNA molecules were made from replicating the 1576 atom RNA molecule.
37. J. D. Baleja, R. T. Pon, and B. D. Sykes, Biochemistry 29, 4828 (1990).
38. S. E. Lietzke, V. E. Carperos, and C. E. Kundrot, Acta Crystallogr., Sect. D: Biol. Crystallogr. 52, 687 (1996).