TY - JOUR
T1 - Efficient multiscale algorithms for solution of self-consistent eigenvalue problems in real space
AU - Wijesekera, Nimal R.
AU - Feng, Guogang
AU - Beck, Thomas L.
PY - 2007/3/2
Y1 - 2007/3/2
N2 - Real-space multiscale methods provide efficient algorithms for large-scale electronic structure calculations. In this paper, we present multigrid strategies for solving self-consistent problems in density functional theory. The full approximation scheme (FAS) formulation of the multigrid method allows for transfer of the expensive orthogonalization and Ritz projection operations to coarse levels. In addition, the effective potential may be updated on coarse levels during multiscale processing of the eigenfunctions. We investigate modifications of a previously proposed algorithm which are necessary to yield robust convergence rates. With these modifications, rapid convergence is observed without orthonormalization or Ritz projection for the full occupied subspace on the fine level. Calculations comparing the various algorithms are performed on three many-electron examples: benzene, benzenedithiol, and the amino acid glycine. The modified algorithm is also illustrated on several larger test cases. Recently developed relativistic separable dual-space Gaussian pseudopotentials are utilized to remove the core electrons.
AB - Real-space multiscale methods provide efficient algorithms for large-scale electronic structure calculations. In this paper, we present multigrid strategies for solving self-consistent problems in density functional theory. The full approximation scheme (FAS) formulation of the multigrid method allows for transfer of the expensive orthogonalization and Ritz projection operations to coarse levels. In addition, the effective potential may be updated on coarse levels during multiscale processing of the eigenfunctions. We investigate modifications of a previously proposed algorithm which are necessary to yield robust convergence rates. With these modifications, rapid convergence is observed without orthonormalization or Ritz projection for the full occupied subspace on the fine level. Calculations comparing the various algorithms are performed on three many-electron examples: benzene, benzenedithiol, and the amino acid glycine. The modified algorithm is also illustrated on several larger test cases. Recently developed relativistic separable dual-space Gaussian pseudopotentials are utilized to remove the core electrons.
UR - http://www.scopus.com/inward/record.url?scp=33847395478&partnerID=8YFLogxK
U2 - 10.1103/PhysRevB.75.115101
DO - 10.1103/PhysRevB.75.115101
M3 - Article
AN - SCOPUS:33847395478
SN - 1098-0121
VL - 75
JO - Physical Review B - Condensed Matter and Materials Physics
JF - Physical Review B - Condensed Matter and Materials Physics
IS - 11
M1 - 115101
ER -