A simplified density matrix minimization for linear scaling self-consistent field theory
J. Chem. Phys. 110, 2332 (1999); doi:10.1063/1.477969
Issue Date: 1 February 1999
You are not logged in to this journal. Log in
A simplified version of the Li, Nunes and Vanderbilt [Phys. Rev. B 47, 10891 (1993)] and Daw [Phys. Rev. B 47, 10895 (1993)] density matrix minimization is introduced that requires four fewer matrix multiplies per minimization step relative to previous formulations. The simplified method also exhibits superior convergence properties, such that the bulk of the work may be shifted to the quadratically convergent McWeeny purification, which brings the density matrix to idempotency. Both orthogonal and nonorthogonal versions are derived. The AINV algorithm of Benzi, Meyer, and T
ma [SIAM J. Sci. Comp. 17, 1135 (1996)] is introduced to linear scaling electronic structure theory, and found to be essential in transformations between orthogonal and nonorthogonal representations. These methods have been developed with an atom-blocked sparse matrix algebra that achieves sustained megafloating point operations per second rates as high as 50% of theoretical, and implemented in the MondoSCF suite of linear scaling SCF programs. For the first time, linear scaling HartreeFock theory is demonstrated with three-dimensional systems, including water clusters and estane polymers. The nonorthogonal minimization is shown to be uncompetitive with minimization in an orthonormal representation. An early onset of linear scaling is found for both minimal and double zeta basis sets, and crossovers with a highly optimized eigensolver are achieved. Calculations with up to 6000 basis functions are reported. The scaling of errors with system size is investigated for various levels of approximation. ©1999 American Institute of Physics.
| History: | Received 1 September 1998; accepted 7 October 1998 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/110/2332/1 |
KEYWORDS and PACS
SELF-CONSISTENT FIELD,
DENSITY MATRIX,
HARTREE-FOCK METHOD,
ELECTRONIC STRUCTURE,
minimisation,
SCF calculations,
convergence of numerical methods
- 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.60.-x
Mathematical methods in physics Numerical approximation and analysis - YEAR: 1999
RELATED DATABASES
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (97)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- D. L. Strout and G. E. Scuseria, J. Chem. Phys. 102, 8448 (1995).
- J. Almlöf, K. Faegri, and K. Korsell,
J. Comput. Chem. 3, 385 (1982) . - M. Häser and R. Ahlrichs,
J. Comput. Chem. 10, 104 (1989) . - C. A. White, B. Johnson, P. Gill, and M. Head-Gordon,
Chem. Phys. Lett. 230, 8 (1994) . - 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).
- M. Challacombe, E. Schwegler, and J. Almlöf, in Computational Chemistry: Review of Current Trends, edited by J. Leczszynski (World Scientific, Singapore, 1996), pp. 53107.
- E. Schwegler and M. Challacombe, J. Chem. Phys. 105, 2726 (1996).
- J. P. Dombroski, S. W. Taylor, and P. M. W. Gill,
J. Phys. Chem. 100, 6272 (1996) . - J. C. Burant, R. E. Stratmann, and M. J. Frisch, J. Chem. Phys. 105, 8969 (1996).
- R. E. Stratmann, G. E. Scuseria, and M. J. Frisch,
Chem. Phys. Lett. 257, 213 (1996) . - M. C. Strain, G. E. Scuseria, and M. J. Frisch,
Science 271, 51 (1996) . - M. Challacombe and E. Schwegler, J. Chem. Phys. 106, 5526 (1997).
- E. Schwelger, M. Challacombe, and M. Head-Gordon, J. Chem. Phys. 106, 9708 (1997).
- C. Ochsenfeld, C. A. White, and M. Head-Gordon, J. Chem. Phys. 109, 1663 (1998).
- E. Schwegler and M. Challacombe (unpublished).
- X. P. Li, R. W. Nunes, and D. Vanderbilt, Phys. Rev. B 47, 10891 (1993).
- M. S. Daw, Phys. Rev. B 47, 10895 (1993).
- R. W. Nunes and D. Vanderbilt, Phys. Rev. B 50, 17611 (1994).
- S. Y. Qiu, C. Z. Wang, K. M. Ho, and C. T. Chan,
J. Phys.: Condens. Matter 6, 9153 (1994) . - E. Hernández and M. J. Gillan, Phys. Rev. B 51, 10157 (1995).
- E. Hernández, M. J. Gillan, and C. Goringe, Phys. Rev. B 53, 7147 (1996).
- J. M. Millam and G. E. Scuseria, J. Chem. Phys. 106, 5569 (1997).
- C. Ochsenfeld and M. Head-Gordon,
Chem. Phys. Lett. 270, 399 (1997) . - F. Mauri, G. Galli, and R. Car, Phys. Rev. B 47, 9973 (1993).
- P. Ordejón, D. A. Drabold, M. P. Grumbach, and R. M. Martin, Phys. Rev. B 48, 14646 (1993).
- J. Kim, F. Mauri, and G. Galli, Phys. Rev. B 51, 1456 (1995).
- D. Sanchezportal, P. Ordejon, E. Artacho, and J. M. Soler,
Int. J. Quantum Chem. 65, 453 (1997) . - P. Ordejon, E. Artacho, and J. M. Soler, Phys. Rev. B 53, 10441 (1996).
- S. Goedecker and L. Colombo, Phys. Rev. Lett. 73, 122 (1994).
- S. Goedecker, J. Comput. Chem. 118, 261 (1995).
- Y. Huang, D. J. Kouri, and D. K. Hoffman,
Chem. Phys. Lett. 243, 367 (1995) . - A. F. Voter, J. D. Kress, and R. N. Silver, Phys. Rev. B 53, 12733 (1996).
- H. Roder, R. N. Silver, D. A. Drablold, and J. J. Dong, Phys. Rev. B 55, 15382 (1997).
- F. Gagel, J. Comput. Chem. 139, 399 (1998).
- G. Galli,
Curr. Opin. Solid State Mater. Sci. 1, 864 (1996) . - D. R. Bowler et al.,
Modell. Simul. Mater. Sci. Eng. 5, 199 (1997) . - S. Goedecker, Rev. Mod. Phys. (submitted).
- W. Kohn,
Int. J. Quantum Chem. 56, 229 (1995) . - P. Maslen et al.,
J. Phys. Chem. A 102, 2215 (1998) . - E. Hernández, M. J. Gillan, and C. M. Goringe, Phys. Rev. B 55, 13485 (1997).
- M. Lepetit, L. Lafon, and X. Lafage,
Int. J. Quantum Chem. 64, 411 (1997) . - S. Q. Wei and M. Y. Chou, Phys. Rev. Lett. 76, 2650 (1996).
- C. J. Tymczak and X. Q. Wang, Phys. Rev. Lett. 78, 3654 (1997).
- R. A. Lippert, T. A. Arias, and A. Edelman, J. Comput. Chem. 140, 278 (1998).
- J. R. Chelikowsky, N. Troullier, and Y. Saad, Phys. Rev. Lett. 72, 1240 (1994).
- J. R. Chelikowsky, N. Troullier, K. Wu, and Y. Saad, Phys. Rev. B 50, 11355 (1994).
- E. L. Briggs, D. J. Sulivan, and J. Bernholc, Phys. Rev. B 52, R5471 (1995).
- R. McWeeny, Rev. Mod. Phys. 126, 1028 (1962).
- M. Benzi and C. D. Meyer,
SIAM J. Sci. Comput. 16, 1159 (1995) . - M. Benzi, C. D. Meyer, and M. T
ma,
SIAM J. Sci. Comput. 17, 1135 (1996) . - M. Challacombe and E. Schwegler, MondoSCF a suite of programs for linear scaling SCF theory (unpublished).
- P. O. Löwdin,
Phys. Rev. 97, 1490 (1955) . - A. Szabo and N. S. Ostlund, Modern Quantum Chemistry, 1st revised ed. (Mc GrawHill. New York, 1989).
- C. C. J. Roothaan,
Rev. Mod. Phys. 23, 69 (1951) . - J. H. Wilkinson, The Algebraic Eigenvalue Problem (Clarendon, Oxford, 1965).
- G. W. Stewart, Introduction to Matrix Computations (Academic, London, 1973).
- P. O. Löwdin, J. Chem. Phys. 18, 365 (1950).
- P. O. Löwdin, Adv. Phys. 5, 3 (1956).
- M. Challacombe, M. Benzi, and M. M. T
ma (unpublished).
- C. A. White, P. Maslen, M. S. Lee, and M. Head-Gordon,
Chem. Phys. Lett. 276, 133 (1997) . - M. Benzi and M. T
ma, Technical Report No. LA-UR-98-2175, Los Alamos National Laborarory, Los Alamos, NM (unpublished).
- M. Benzi, J. Marin, and M. T
ma, in Fourth IMACS International Symposium on Iterative Methods in Scientific Computation (IMACS, Austin, TX, 1998).
- F. Liu,
J. Mol. Struct.: THEOCHEM 230, 47 (1991) . - L. Cohen and C. Frishberg, J. Chem. Phys. 65, 4234 (1976).
- C. Frishberg, L. Cohen, and P. Blumenau, Int. J. Quantum Chem., Symp. 14, 161 (1980).
- L. Massa and L. Cohen, Int. J. Quantum Chem., Symp. 14, 167 (1980).
- A. Redondo, Phys. Rev. A 39, 4366 (1989).
- A. Redondo and J. C. Marshall, J. Chem. Phys. 91, 5492 (1989).
- C. Waggoner and L. L. Combs,
J. Optim. Theory Appl. 76, 225 (1993) . - P. Fantucci and S. Polenzzo,
Int. J. Quantum Chem. 52, 817 (1994) . - E. Anderson et al., LAPACK Users' Guide, release 2.0 ed. (SIAM, 1994).
- W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipies in FORTRAN (Cambridge University Press, Port Chester, NY, 1992).
- L. C. W. Dixon, in Nonlinear Optimization Theory and Algorithms, edited by L. C. W. Dixon, E. Spedicato, and G. P. Szegö (Birkhäuser, Boston, 1980), pp. 124135.
- J. R. Shewchuk, Technical Report No. CMU-CS-94-125, Carnegie Mellon University, Pittsburgh, PA (unpublished).
- L. C. W. Dixon, P. G. Ducksbury, and P. Singh,
J. Optim. Theory Appl. 47, 285 (1985) . - J. F. Annett,
Comput. Mater. Sci. 4, 23 (1995) . - F. Liu,
J. Mol. Struct.: THEOCHEM 226, 197 (1991) . - A. G. Buckley,
Math. Program. 15, 200 (1978) . - P. Pulay,
Chem. Phys. Lett. 73, 393 (1980) . - R. T. McLay, S. Swift, and G. F. Carey,
J. Par. Dist. Comp. 37, 146 (1996) . - S. Goedecker et al., Technical Report No. LA-UR-97-1504, Los Alamos National Laborarory, Los Alamos, NM (unpublished).
- M. R. Field,
SIAM J. Sci. Comput. 19, 27 (1998) . - S. Pissanetzky, Sparse Matrix Technology (Academic, London, 1984).
- Y. Saad, Iterative Methods for Sparse Linear Systems (PWS, Boston, MA, 1996).
- I. S. Duff, A. M. Erisman, and J. K. Reid, Direct Methods for Sparse Matrices (Oxford University Press, London, 1986).
- F. G. Gustavson,
ACM Trans. Math. Softw. 4, 250 (1978) . - G. Golub and C. F. van Loan, Matrix Computations (Johns Hopkins University Press, Baltimore, MD, 1996).
- A. D. Daniels, J. M. Millam, and G. E. Scuseria, J. Chem. Phys. 107, 425 (1997).
- M. Challacombe and E. Schwegler (unpublished).
- P. Pulay,
J. Comput. Chem. 3, 556 (1982) . - E. Polak, Computational Methods in Optimization: A Unified Approach (Academic, London, 1971).
- M. J. D. Powell,
Math. Program. 12, 241 (1977) . - T. R. Gibbs and A. Popolato, LASL Explosive Property Data (University of California Press, Berkeley, CA, 1980), p. 109.
- M. J. Frisch et al., GAUSSIAN 94, Revision E.2, Gaussian Inc., Pittsburgh, PA, 1995.
- M. Smotherman, program mm.c, available by anonymous ftp from ftp.nosc.mil under/pub/aburto/mm. Compiled using cc -DN = 500 -DUNIX -02 -64 -r10000 and executed as mm -w 50, 1997.
- M. J. Daydé and I. S. Duff,
Lect. Notes Comput. Sci. 1215, 108 (1997) .








