Density-matrix-based algorithm for solving eigenvalue problems
Abstract
A fast and stable numerical algorithm for solving the symmetric eigenvalue problem is presented. The technique deviates fundamentally from the traditional Krylov subspace iteration based techniques (Arnoldi and Lanczos algorithms) or other Davidson-Jacobi techniques and takes its inspiration from the contour integration and density-matrix representation in quantum mechanics. It will be shown that this algorithm—named FEAST—exhibits high efficiency, robustness, accuracy, and scalability on parallel architectures. Examples from electronic structure calculations of carbon nanotubes are presented, and numerical performances and capabilities are discussed.
- Publication:
-
Physical Review B
- Pub Date:
- March 2009
- DOI:
- 10.1103/PhysRevB.79.115112
- arXiv:
- arXiv:0901.2665
- Bibcode:
- 2009PhRvB..79k5112P
- Keywords:
-
- 02.60.-x;
- 02.70.Hm;
- 31.15.-p;
- Numerical approximation and analysis;
- Spectral methods;
- Calculations and mathematical techniques in atomic and molecular physics;
- Computer Science - Computational Engineering;
- Finance;
- and Science;
- Computer Science - Mathematical Software
- E-Print:
- 7 pages, 3 figures