NURBS-based non-periodic finite element framework for Kohn-Sham density functional theory calculations

被引：6
作者：
Temizer İ. [1 ]
Motamarri P. [2 ]
Gavini V. [2 ,3 ]
机构：
[1] Department of Mechanical Engineering, Bilkent University, Ankara
[2] Department of Mechanical Engineering, University of Michigan, Ann Arbor, 48109, MI
[3] Department of Materials Science and Engineering, University of Michigan, Ann Arbor, 48109, MI
来源：
Journal of Computational Physics | 2020年 / 410卷
关键词：
All-electron; B-splines; Density functional theory; Finite elements; Pseudopotential; Real space;
D O I：
10.1016/j.jcp.2020.109364
中图分类号：
学科分类号：
摘要：
A real-space non-periodic computational framework is developed for Kohn-Sham density functional theory (DFT). The electronic structure calculation framework is based on the finite element method (FEM) where the underlying basis is chosen as non-uniform rational B-splines (NURBS) which display continuous higher-order derivatives. The framework is formulated within a unified presentation that can simultaneously address both all-electron and pseudopotential settings in radial and three-dimensional cases. The canonical Kohn-Sham equation and the Poisson equation are discretized on different meshes in order to ensure that the underlying variational structural of Kohn-Sham DFT is preserved within the weak formulation of FEM. The discrete generalized eigenvalue problem emanating from the Kohn-Sham equation is solved efficiently based on the Chebyshev-filtered subspace iteration method. Numerical investigations in the radial case demonstrate all-electron and local pseudopotential capabilities on single atoms. In the three-dimensional case, all-electron and nonlocal pseudopotential computations on single atoms and small molecules are followed by local and nonlocal pseudopotential studies on larger systems. At all stages, special care is taken to demonstrate optimal convergence rates towards the ground state energy with chemical accuracy. Comparisons with classical Lagrange basis sets indicate the significantly higher per-degree-of-freedom accuracy displayed by NURBS. Specifically, cubic NURBS discretizations can offer a faster route to a prescribed accuracy than even sixth-order Lagrange discretizations on comparable meshes, thereby indicating considerable efficiency gains which are possible with these higher-order basis sets within effective numerical implementations. © 2020 Elsevier Inc.
引用
收藏
相关论文
共 109 条
[1]  
Payne M.C., Teter M.P., Allan D.C., Arias T.A., Joannopoulos J.D., Iterative minimization techniques for ab initio total-energy calculations: molecular dynamics and conjugate gradients, Rev. Mod. Phys., 64, pp. 1045-1097, (1992)
[2]  
Szabo A., Ostlund N.S., Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, (1996)
[3]  
Helgaker T., Jorgensen P., Olsen J., Molecular Electronic-Structure Theory, (2000)
[4]  
Parr R.G., Yang W., Density-Functional Theory of Atoms and Molecules, (1989)
[5]  
Koch W., Holthausen M.C., A Chemist's Guide to Density Functional Theory, (2001)
[6]  
Beck T.L., Real-space mesh techniques in density-functional theory, Rev. Mod. Phys., 72, pp. 1041-1080, (2000)
[7]  
Martin R.M., Electronic Structure: Basis Theory and Practical Methods, (2004)
[8]  
Finnis M., Interatomic Forces in Condensed Matter, (2003)
[9]  
Batcho P.F., Computational method for general multicenter electronic structure calculations, Phys. Rev. E, 61, pp. 7169-7183, (2000)
[10]  
Modine N.A., Zumbach G., Kaxiras E., Adaptive-coordinate real-space electronic-structure calculations for atoms, molecules, and solids, Phys. Rev. B, 55, 10, (1997)
12345678910