Comparison of Conjugate Gradient Density Matrix Search and Chebyshev Expansion Methods for Avoiding Diagonalization in Large-Scale Electronic Structure Calculations

We report a comparison of two linear-scaling methods which avoid the diagonalization bottleneck of traditional electronic structure algorithms. The Chebyshev expansion method (CEM) is implemented for carbon tight-binding calculations of large systems and its memory and timing requirements compared to those of our previously implemented conjugate gradient density matrix search (CG-DMS). Benchmark calculations are carried out on icosahedral fullerenes from C60 to C8640 and the linear scaling memory and CPU requirements of the CEM demonstrated. We show that the CPU requisites of the CEM and CG-DMS are similar for calculations with comparable accuracy.

[1]  David R. Bowler,et al.  Tight-binding modelling of materials , 1997 .

[2]  Donald J. Kouri,et al.  Direct approach to density functional theory: iterative treatment using a polynomial representation of the Heaviside step function operator , 1995 .

[3]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[4]  Wang.  Calculating the density of states and optical-absorption spectra of large quantum systems by the plane-wave moments method. , 1994, Physical review. B, Condensed matter.

[5]  R. Mcweeny.  Some Recent Advances in Density Matrix Theory , 1960 .

[6]  C. Wang,et al.  Tight-binding molecular dynamics with linear system-size scaling , 1994 .

[7]  Che Ting Chan,et al.  A transferable tight-binding potential for carbon , 1992 .

[8]  Youhong Huang,et al.  Direct Approach to Density Functional Theory: Heaviside−Fermi Level Operator Using a Pseudopotential Treatment , 1996 .

[9]  Michael J. Frisch,et al.  Achieving linear scaling in exchange-correlation density functional quadratures , 1996 .

[10]  Michael J. Frisch,et al.  A linear scaling method for Hartree–Fock exchange calculations of large molecules , 1996 .

[11]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[12]  Gustavo E. Scuseria,et al.  An O(N) tight-binding study of carbon clusters up to C8640: the geometrical shape of the giant icosahedral fullerenes , 1996 .

[13]  Gustavo E. Scuseria,et al.  Semiempirical methods with conjugate gradient density matrix search to replace diagonalization for molecular systems containing thousands of atoms , 1997 .

[14]  Michael J. Frisch,et al.  Achieving Linear Scaling for the Electronic Quantum Coulomb Problem , 1996, Science.

[15]  Sergio Pissanetzky,et al.  Sparse Matrix Technology , 1984 .

[16]  Colombo,et al.  Efficient linear scaling algorithm for tight-binding molecular dynamics. , 1994, Physical review letters.

[17]  R. Silver,et al.  Kernel polynomial method for a nonorthogonal electronic-structure calculation of amorphous diamond , 1997 .

[18]  Daw.  Model for energetics of solids based on the density matrix. , 1993, Physical review. B, Condensed matter.

[19]  R. Silver,et al.  DENSITIES OF STATES OF MEGA-DIMENSIONAL HAMILTONIAN MATRICES , 1994 .

[20]  Gustavo E. Scuseria,et al.  Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations , 1997 .

[21]  Lin-wang Wang,et al.  Solving Schrödinger’s equation around a desired energy: Application to silicon quantum dots , 1994 .

[22]  Martin Head-Gordon,et al.  Chebyshev expansion methods for electronic structure calculations on large molecular systems , 1997 .

[23]  M. Teter,et al.  Tight-binding electronic-structure calculations and tight-binding molecular dynamics with localized orbitals. , 1994, Physical review. B, Condensed matter.