A new method for diagonalising large matrices
- 21 June 1985
- journal article
- Published by IOP Publishing in Journal of Physics A: General Physics
- Vol. 18 (9), 1343-1359
- https://doi.org/10.1088/0305-4470/18/9/018
Abstract
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.Keywords
This publication has 17 references indexed in Scilit:
- New approach for solving the density-functional self-consistent-field problemPhysical Review B, 1982
- Ionicity and the theory of Schottky barriersPhysical Review B, 1977
- Modification of Davidson's method for the calculation of eigenvalues and eigenvectors of large real-symmetric matrices: “root homing procedure”Journal of Computational Physics, 1976
- The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matricesJournal of Computational Physics, 1975
- Algorithm for Diagonalization of Large MatricesThe Journal of Chemical Physics, 1965
- Gaussian-Type Functions for Polyatomic Systems. IThe Journal of Chemical Physics, 1965
- Electronic Spectra of Crystalline Germanium and SiliconPhysical Review B, 1964
- A Note on the Quantum-Mechanical Perturbation TheoryThe Journal of Chemical Physics, 1951
- An iteration method for the solution of the eigenvalue problem of linear differential and integral operatorsJournal of Research of the National Bureau of Standards, 1950
- The solution of natural frequency equations by relaxation methodsQuarterly of Applied Mathematics, 1948