Article ID Journal Published Year Pages File Type
6892323 Computers & Mathematics with Applications 2017 10 Pages PDF
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).
Related Topics
Physical Sciences and Engineering Computer Science Computer Science (General)
Authors
,