Atomic electronic structure calculations with Hermite interpolating polynomials

Abstract

We have recently described the implementation of atomic electronic structure calculations within the finite element method with numerical radial basis functions of the form ${\chi }_{\mu }(r)=r^{-1}{B}_{\mu }(r)$, where high-order Lagrange interpolating polynomials (LIPs) were used as the shape functions $B_{\mu }(r)$ [S. Lehtola, Int. J. Quantum Chem. 119, e25945 (2019)]. In this work, we discuss how ${\chi }_{\mu }(r)$ can be evaluated in a stable manner at small $r$ and also revisit the choice of the shape functions $B_{\mu }(r)$. Three kinds of shape functions are considered: in addition to the ${\mathcal{C}}^0$ continuous LIPs, we consider the analytical implementation of first-order Hermite interpolating polynomials (HIPs) that are ${\mathcal{C}}^1$ continuous, as well as numerical implementations of $n$-th order (${\mathcal{C}}^n$ 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-$\sigma$ local kinetic energy ${\tau }_{\sigma }$ 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 ${\tau }_{\sigma }$. 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.

arXiv comment26 pages, 9 figures. Fixed typos

References in corpus (21)

• Climbing the Density Functional Ladder: Non-Empirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids
• Nonlocal van der Waals density functional: The simpler the better
• Numerical atomic orbitals for linear scaling
• Systematic generation of finite-range atomic basis sets for linear-scaling calculations
• An assessment of initial guesses for self-consistent field calculations. Superposition of Atomic Potentials: simple yet efficient
• An overview of self-consistent field calculations within finite basis sets
• A review on non-relativistic fully numerical electronic structure calculations on atoms and diatomic molecules
• A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues
• Many recent density functionals are numerically ill-behaved
• Fully numerical calculations on atoms with fractional occupations. Range-separated exchange functionals
• Fully numerical Hartree-Fock and density functional calculations. I. Atoms
• Optimal finite-range atomic basis sets for liquid water and ice
• Fully numerical electronic structure calculations on diatomic molecules in weak to strong magnetic fields
• Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism
• Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets
• Polarized Gaussian basis sets from one-electron ions
• Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules
• On the accurate reproduction of strongly repulsive interatomic potentials
• Meta-local density functionals: a new rung on Jacob's ladder
• How good are recent density functionals for ground and excited states of one-electron systems?
• The importance of basis states: an example using the Hydrogen basis

Cited by in corpus (12)

• Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism
• A call to arms: making the case for more reusable libraries
• Reproducibility of density functional approximations: how new functionals should be reported
• Automatic Generation of Accurate and Cost-efficient Auxiliary Basis Sets
• Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities
• On the accuracy of a recent regularized nuclear potential
• OpenOrbitalOptimizer -- a reusable open source library for self-consistent field calculations
• Ensemble Generalization of the Perdew-Zunger Self-Interaction Correction: a Way Out of Multiple Minima and Symmetry Breaking
• Systematic study of confinement induced effects on atomic electronic structure
• Review of the finite difference Hartree-Fock method for atoms and diatomic molecules, and its implementation in the x2dhf program
• Importance profiles. Visualization of atomic basis set requirements
• Atomic Confinement Potentials and the Generation of Numerical Atomic Orbitals