Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations
J. Chem. Phys. 106, 5569 (1997); doi:10.1063/1.473579
Issue Date: 1 April 1997
You are not logged in to this journal. Log in
Advances in the computation of the Coulomb, exchange, and correlation contributions to Gaussian-based Hartree–Fock and density functional theory Hamiltonians have demonstrated near-linear scaling with molecular size for these steps. These advances leave the
(N3) diagonalization bottleneck as the rate determining step for very large systems. In this work, a conjugate gradient density matrix search (CG-DMS) method has been successfully extended and computationally implemented for use with first principles calculations. A Cholesky decomposition of the overlap matrix and its inverse is used to transform to and back from an orthonormal basis, which can be formed in near-linear time for sparse systems. Linear scaling of CPU time for the density matrix search and crossover of CPU time with diagonalization is demonstrated for polyglycine chains containing up to 493 atoms and water clusters up to 900 atoms. ©1997 American Institute of Physics.
| History: | Received 18 October 1996; accepted 31 December 1996 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/106/5569/1 |
KEYWORDS and PACS
- 31.15.Ew
Electronic structure of atoms, molecules and their ions: theory Calculations and mathematical techniques in atomic and molecular physics (excluding electron correlation calculations) Density-functional theory - 31.15.Ne
Electronic structure of atoms, molecules and their ions: theory Calculations and mathematical techniques in atomic and molecular physics (excluding electron correlation calculations) Self-consistent-field methods - 02.10.Sp
Mathematical methods in physics Logic, set theory, and algebra Linear and multilinear algebra; matrix theory (finite and infinite) - YEAR: 1996-97
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (46)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- J. Almöf, K. Faergi, Jr., and K. Korsell,
J. Comput. Chem. 3, 385 (1982 ). - M. Häser and R. Ahlrichs,
J. Comput. Chem. 10, 104 (1989 ). - D. L. Strout and G. E. Scuseria, J. Chem. Phys. 102, 8448 (1995).
- M. C. Strain, G. E. Scuseria, and M. J. Frisch,
Science 271, 51 (1996 ). - C. A. White, B. G. Johnson, P. M. W. Gill, and M. Head-Gordon,
Chem. Phys. Lett. 253, 268 (1996 ). - M. Challacombe, E. Schwegler, and J. Almlöf, J. Chem. Phys. 104, 4685 (1996).
- 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 ). - E. Schwegler and M. Challacombe, J. Chem. Phys. 105, 2726 (1996).
- W. Kohn,
Int. J. Quantum Chem. 56, 229 (1995 ). - J. C. Burant, G. E. Scuseria, and M. J. Frisch, J. Chem. Phys. 105, 8969 (1996).
- C. Lanczos, J. Res. Nat. Bur. Stand. 45, 255 (1950);
- D. C. Sorensen, SIAM J. Mater. Anal. Appl. 13, 357 (1992).
- R. B. Lehoucq and D. C. Sorensen,
SIAM J. Matr. Anal. Appl. 17, 789 (1996 ). - M. P. Teter, M. C. Payne, D. C. Allan, Phys. Rev. B 40, 12 255 (1989).
- G. Galli and M. Parrinello, Phys. Rev. Lett. 69, 3547 (1992).
- P. Ordejón, D. A. Drabold, M. P. Grumbach, and R. M. Martin, Phys. Rev. B 48, 14 646 (1993).
- F. Mauri, G. Galli, and R. Car, Phys. Rev. B 47, 9973 (1993).
- F. Mauri and G. Galli, Phys. Rev. B 50, 4316 (1994).
- E. B. Stechel, A. R. Williams, and P. J. Feibelman, Phys. Rev. B 49, 10 088 (1994).
- J. P. Stewart,
Int. J. Quantum Chem. 58, 133 (1995 ). - P. Ordejón, D. A. Drabold, R. M. Martin, and M. P. Grumbach, Phys. Rev. B 51, 1456 (1995).
- X. Chen, J.-M. Langlois, and W. A. Goddard III, Phys. Rev. B 52, 2348 (1995).
- X.-P. Li, R. W. Nunes, and D. Vanderbilt, Phys. Rev. B 47, 10 891 (1993).
- M. S. Daw, Phys. Rev. B 47, 10 895 (1993).
- S.-Y. Qiu, C. Z. Wang, K. M. Ho, and C. T. Chan,
J. Phys. Condens. Matter 6, 9153 (1994 ). - E. Hernńdez and M. J. Gillan, Phys. Rev. B 51, 10 157 (1994).
- W. Kohn, Phys. Rev. Lett. 76, 3168 (1996).
- C. H. Xu and G. E. Scuseria,
Chem. Phys. Lett. 262, 219 (1996 ). - A. D. Livelsberger, J. M. Millam, and G. E. Scuseria, J. Chem. Phys. (submitted).
- W. Yang, Phys. Rev. Lett. 66, 1438 (1991).
- W. Yang,
J. Mol. Struc. (Theochem) 255, 461 (1992 ). - T.-S. Lee, D. M. York, and W. Yang, J. Chem. Phys. 105, 2744 (1996).
- R. T. Gallant and A. St-Amant,
Chem. Phys. Lett. 256, 569 (1996 ). - S. L. Dixon and K. M. Merz, Jr., J. Chem. Phys. 104, 6643 (1996).
- R. McWeeny,
Rev. Mod. Phys. 32, 335 (1960 ). - C. Z. Wang, C. T. Chan, and K. M. Ho, Phys. Rev. B 39, 8586 (1989).
- R. W. Nunes and D. Vanderbilt, Phys. Rev. B 50, 17 611 (1994).
- P.-O. Löwdin, J. Chem. Phys. 18, 365 (1950).
- W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in FORTRAN, 2nd ed. (Cambridge University Press, New York, 1992).
- S. Pissanetsky, Sparse Matrix Technology (Academic, New York, 1984).
- The polyglycine chains are similar to those described in Refs. 31–33. Coordinates for these molecules are available from us.
- The water clusters up to 150 water molecules are those used in Ref. 9. We obtained the coordinates from the authors of Ref. 9. Coordinates for the 300 water molecule cluster are available from us.
- 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 95, Development Version (Revision D.4), Gaussian, Inc., Pittsburgh, PA, 1996.
- P. Pulay,
Chem. Phys. Lett. 73, 393 (1980 ). - P. Pulay,
J. Comp. Chem. 3, 556 (1982 ).








