Adaptive Finite Element Method for Solving the Exact Kohn-Sham Equation of Density Functional Theory
Results of the application of an adaptive finite element (FE) based solution using the FETK library of M. Holst to Density Functional Theory (DFT) approximation to the electronic structure of atoms and molecules are reported. The severe problem associated with the rapid variation of the electronic wave functions in the near singular regions of the atomic centers is treated by implementing completely unstructured simplex meshes that resolve these features around atomic nuclei. This concentrates the computational work in the regions in which the shortest length scales are necessary and provides for low resolution in regions for which there is no electron density. The accuracy of the solutions significantly improved when adaptive mesh refinement was applied, and it was found that the essential difficulties of the Kohn-Sham eigenvalues equation were the result of the singular behavior of the atomic potentials. Even though the matrix representations of the discrete Hamiltonian operator in the adaptive finite element basis are always sparse with a linear complexity in the number of discretization points, the overall memory and computational requirements for the solver implemented were found to be quite high. The number of mesh vertices per atom as a function of the atomic number Z and the required accuracy e (in atomic units) was esitmated to be v (e;Z) = 122:37 * Z2:2346 /1:1173 , and the number of floating point operations per minimization step for a system of NA atoms was found to be 0(N3A*v(e,Z0) (e.g. Z=26, e=0.0015 au, and NA=100, the memory requirement and computational cost would be ~0.2 terabytes and ~25 petaflops). It was found that the high cost of the method could be reduced somewhat by using a geometric based refinement strategy to fix the error near the singularities.
- Research Organization:
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States). Environmental Molecular Sciences Lab. (EMSL)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC05-76RL01830
- OSTI ID:
- 958438
- Report Number(s):
- PNNL-SA-62100; 20904; KJ0101010; TRN: US1000379
- Journal Information:
- Journal of Chemical Theory and Computation, 5(4):937-948, Vol. 5, Issue 4
- Country of Publication:
- United States
- Language:
- English
Similar Records
All-electron Kohn–Sham density functional theory on hierarchic finite element spaces
Parallelization of an Adaptive Multigrid Algorithm for Fast Solution of Finite Element Structural Problems
Related Subjects
ACCURACY
APPROXIMATIONS
ATOMIC NUMBER
ATOMS
EIGENVALUES
ELECTRON DENSITY
ELECTRONIC STRUCTURE
FINITE ELEMENT METHOD
FUNCTIONALS
HAMILTONIANS
MINIMIZATION
NUCLEI
RESOLUTION
WAVE FUNCTIONS
Environmental Molecular Sciences Laboratory