High-order finite element method for atomic structure calculations

Abstract

We introduce \texttt{featom}, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. To address the slow convergence of the $\kappa =\pm 1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy ($10^{-8}$ Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers ($Z$) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver \texttt{dftatom}. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

arXiv comment30 pages, 3 tables, 6 figures, accepted in Elsevier's Computer Physics Communications

References in corpus (8)

• Libxc: a library of exchange and correlation functionals for density functional theory
• Dual kinetic balance approach to basis set expansions for the Dirac equation
• Periodic Pulay method for robust and efficient convergence acceleration of self-consistent field iterations
• A B-spline Galerkin method for the Dirac equation
• Application of the dual-kinetic-balance sets in the relativistic many-body problem of atomic structure
• dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations
• Variation and Series Approach to the Thomas-Fermi Equation
• Resolving the spurious-state problem in the Dirac equation with finite difference method

Cited by in corpus (4)

• dftatom: A robust and general Schrödinger and Dirac solver for atomic structure calculations
• Relativistic atomic structure calculations in support of spectroscopy
• Systematic study of confinement induced effects on atomic electronic structure
• Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory