Papers › Atomic electronic structure calculations with Hermite interpolating polynomials
Atomic electronic structure calculations with Hermite interpolating polynomials
Susi Lehtola
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
We have recently described the implementation of atomic electronic structure calculations within the finite element method with numerical radial basis functions of the form χ_μ(r)=r⁻¹B_μ(r), where high-order Lagrange interpolating polynomials (LIPs) were used as the shape functions B_μ(r) [S. Lehtola, Int. J. Quantum Chem. 119, e25945 (2019)]. In this work, we discuss how χ_μ(r) can be evaluated in a stable manner at small r and also revisit the choice of the shape functions B_μ(r). Three kinds of shape functions are considered: in addition to the 𝒞⁰ continuous LIPs, we consider the analytical implementation of first-order Hermite interpolating polynomials (HIPs) that are 𝒞¹ continuous, as well as numerical implementations of n-th order (𝒞ⁿ continuous) HIPs that are expressed in terms of an underlying high-order LIP basis. Furnished with the new implementation, we demonstrate that the first-order HIPs are reliable even with large numbers of nodes and that they also work with non-uniform element grids, affording even better results in atomic electronic structure calculations than LIPs with the same total number of basis functions. We demonstrate that discontinuities can be observed in the spin-σ local kinetic energy τ_σ in small LIP basis sets, while HIP basis sets do not suffer from such issues; however, either set can be used to reach the complete basis set limit with smooth τ_σ. Moreover, we discuss the implications of HIPs on calculations with meta-GGA functionals with a number of recent meta-GGA functionals, and find most Minnesota functionals to be ill-behaved. We also examine the potential usefulness of the explicit control over the derivative in HIPs for forming numerical atomic orbital basis sets, but find that confining potentials are still likely a better option.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections