Spectral scheme for atomic structure calculations in density functional theory
Sayan Bhowmik, John E. Pask, Andrew J. Medford, Phanish Suryanarayana
arXiv:2406.00534·physics.comp-ph·Published 2024-06-01·Updated 2024-06-06
We present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number $Z = 1$ to $83$.
TopicsQuantum Chemistry & Force Fields
Tagsdft
arXiv categoriesphysics.comp-ph, cond-mat.mtrl-sci
arXiv abstract pagePDF