Skip to main content

Efficient iterative method for solving the Dirac–Kohn–Sham density functional theory

Author(s): Lin, Lin; Shao, Sihong; E, Weinan

To refer to this page use:
Abstract: We present for the first time an efficient iterative method to directly solve the fourcomponent Dirac–Kohn–Sham (DKS) density functional theory. Due to the existence of the negative energy continuum in the DKS operator, the existing iterative techniques for solving the Kohn–Sham systems cannot be efficiently applied to solve the DKS systems. The key component of our method is a novel filtering step (F) which acts as a preconditioner in the framework of the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The resulting method, dubbed the LOBPCG-F method, is able to compute the desired eigenvalues and eigenvectors in the positive energy band without computing any state in the negative energy band. The LOBPCG-F method introduces mild extra cost compared to the standard LOBPCG method and can be easily implemented. We demonstrate our method in the pseudopotential framework with a planewave basis set which naturally satisfies the kinetic balance prescription. Numerical results for Pt2, Au2, TlF, and Bi2Se3 indicate that the LOBPCG-F method is a robust and efficient method for investigating the relativistic effect in systems containing heavy elements
Publication Date: 15-Jul-2013
Electronic Publication Date: 27-Mar-2013
Citation: Lin, Lin, Shao, Sihong, E, Weinan. (2013). Efficient iterative method for solving the Dirac–Kohn–Sham density functional theory. Journal of Computational Physics, 245 (205 - 217. doi:10.1016/
DOI: doi:10.1016/
ISSN: 0021-9991
Pages: 205 - 217
Type of Material: Journal Article
Journal/Proceeding Title: Journal of Computational Physics
Version: Author's manuscript

Items in OAR@Princeton are protected by copyright, with all rights reserved, unless otherwise indicated.