SPARC-atomSFE: Spectral finite-element package for atomic structure calculations in density functional theory

We present SPARC-atomSFE, a spectral finite-element package for accurate and efficient atomic structure calculations within the framework of Kohn-Sham density functional theory. The package supports both all-electron and norm conserving pseudopotential calculations across a comprehensive hierarchy of exchange-correlation approximations, spanning local, semilocal, and nonlocal functionals. The latter includes hybrid functionals and the many-body random phase approximation, for which we implement both the generalized Kohn-Sham approach and the optimized effective potential (OEP) method, with OEP necessary for eigenvalue-dependent functionals. Spatial discretization is based on an adaptive grid with element nodes distributed according to the Legendre–Gauss–Lobatto scheme, high-order $C^0$-continuous Lagrange polynomial basis functions, and Gauss–Legendre quadrature for numerical integration. We present systematic convergence studies and identify the computational parameters required to achieve target accuracies. We validate the accuracy of SPARC-atomSFE through representative calculations spanning the various exchange-correlations approximations, obtaining results that generally agree with values in the literature to within $1\mu \text{Ha}$ or better.