A new method for diagonalising large matrices

The structure and implementation of a new general iterative method for diagonalising large matrices (the 'residual minimisation/direct inversion in the iterative subspace' method of Bendt and Zunger) are described and contrasted with other more commonly used iterative techniques. The method requires the direct diagonalisation of only a small submatrix, does not require the storage of the large matrix and provides eigensolutions to within a prescribed precision in a rapidly convergent iterative procedure. Numerical results for two rather different matrices (a real 50*50 non-diagonally dominant matrix and a complex Hermitian 181*181 matrix corresponding to the pseudopotential band structure of a semiconductor in a plane wave basis set) are used to compare the new method with the competing methods. the new method converges quickly and should be the most efficient for very large matrices in terms both of computation time and central storage requirements; it is quite insensitive to the properties of the matrices used. This technique makes possible efficient solution of a variety of quantum mechanical matrix problems where large basis set expansions are required.

[1]  J. Cooper.  The solution of natural frequency equations by relaxation methods , 1948 .

[2]  C. Lanczos.  An iteration method for the solution of the eigenvalue problem of linear differential and integral operators , 1950 .

[3]  Per-Olov Löwdin,et al.  A Note on the Quantum‐Mechanical Perturbation Theory , 1951 .

[4]  David Brust,et al.  Electronic Spectra of Crystalline Germanium and Silicon , 1964 .

[5]  R. Nesbet.  Algorithm for Diagonalization of Large Matrices , 1965 .

[6]  I. Shavitt.  Modification of Nesbet's algorithm for the iterative evaluation of eigenvalues and eigenvectors of large matrices , 1970 .

[7]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[8]  W. Butscher,et al.  Modification of Davidson's Method for the Calculation of Eigenvalues and Eigenvectors of Large Real-Symmetric Matrices: , 1976 .

[9]  Steven G. Louie,et al.  Ionicity and the theory of Schottky barriers , 1977 .

[10]  Alex Zunger,et al.  First-principles nonlocal-pseudopotential approach in the density-functional formalism: Development and application to atoms , 1978 .

[11]  A. Zunger.  Contemporary pseudopotentials—Simple reliability criteria , 1979 .

[12]  R. Raffenetti.  A simultaneous coordinate relaxation algorithm for large, sparse matrix eigenvalue problems , 1979 .

[13]  A. Zunger,et al.  First-principles nonlocal-pseudopotential approach in the density-functional formalism. II. Application to electronic and structural properties of solids , 1979 .

[14]  M. T. Hoor.  Gaussian expansions of hydrogenic functions involving few variational parameters , 1980 .

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

[16]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[17]  A. Zunger,et al.  New approach for solving the density-functional self-consistent-field problem , 1982 .

[18]  I. Duff.  Sparse Matrices and Their Uses. , 1983 .