Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
6892323 | Computers & Mathematics with Applications | 2017 | 10 Pages |
Abstract
A method for solving Schrödinger's equation that combines the Rayleigh quotient iteration (RQI) with the conjugate residual (CR) method is described. It is assumed that the eigenfunctions of the Hamiltonian operator HË are expressed as linear combinations of plane waves, which allows the efficient computation of the product Hâ
x via fast Fourier transform, where H is the representation of HË in the plane-wave basis and x an arbitrary vector. The algorithm carries out two kinds of iterations: outer-loop RQI iterations (very few) and inner-loop CR iterations (many more). A very small number of orthogonalizations of the set of eigenvectors is required by the proposed method, just one per RQI iteration, which is one of its key advantages because, for large enough systems, orthogonalization is the costliest component of the calculation. The procedure is much faster than matrix diagonalization, and its cost with respect to the size of the plane-wave basis (NPW) has an effective scaling ranging from linear O(NPW) to quadratic O(NPW2).
Keywords
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science (General)
Authors
José M. Pérez-Jordá,