Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism

Abstract

Density functional calculations on atoms are often used for determining accurate initial guesses as well as generating various types of pseudopotential approximations and efficient atomic-orbital basis sets for polyatomic calculations. To reach the best accuracy for these purposes, the atomic calculations should employ the same density functional as the polyatomic calculation. Atomic density functional calculations are typically carried out employing spherically symmetric densities, corresponding to the use of fractional orbital occupations. We have described their implementation for density functional approximations (DFAs) belonging to the local density approximation (LDA) and generalized gradient approximation (GGA) levels of theory as well as Hartree-Fock (HF) and range-separated exact exchange [S. Lehtola, Phys. Rev. A 2020, 101, 012516]. In this work, we describe the extension to meta-GGA functionals using the generalized Kohn-Sham scheme, in which the energy is minimized with respect to the orbitals, which in turn are expanded in the finite element formalism with high-order numerical basis functions. Furnished with the new implementation, we continue our recent work on the numerical well-behavedness of recent meta-GGA functionals [S. Lehtola and M. A. L. Marques, J. Chem. Phys. 2022, 157, 174114]. We pursue complete basis set (CBS) limit energies for recent density functionals, and find many to be ill-behaved for the Li and Na atoms. We report basis set truncation errors (BSTEs) of some commonly used Gaussian basis sets for these density functionals and find the BSTEs to be strongly functional dependent. We also discuss the importance of density thresholding in DFAs and find that all of the functionals studied in this work yield total energies converged to $0.1\mu E_h$ when densities smaller than $10^{-11}{a}_0^{-3}$ are screened out.

arXiv comment30 pages, 2 figures. Fixed some typos

References in corpus (25)

• Climbing the Density Functional Ladder: Non-Empirical Meta-Generalized Gradient Approximation Designed for Molecules and Solids
• QuantumATK: An integrated platform of electronic and atomic-scale modelling tools
• SIESTA: recent developments and applications
• Localized atomic basis set in the projector augmented wave method
• Regularized SCAN functional
• r2SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications
• Issues and Challenges in Orbital-free Density Functional Calculations
• RESCU: a Real Space Electronic Structure Method
• The Elephant in the Room of Density Functional Theory Calculations
• Large scale and linear scaling DFT with the CONQUEST code
• Deorbitalization strategies for meta-GGA exchange-correlation functionals
• 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
• Kohn-Sham Kinetic Energy Density in the Nuclear and Asymptotic Regions: Deviations from the Von Weizsäcker Behavior and Applications to Density Functionals
• 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
• Atomic electronic structure calculations with Hermite interpolating polynomials
• 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
• Ensemble v-representable ab-initio density functional calculation of energy and spin in atoms: atest of exchange-correlation approximations
• Implementation of self-consistent MGGA functionals in augmented plane wave based methods
• How good are recent density functionals for ground and excited states of one-electron systems?

Cited by in corpus (12)

• A call to arms: making the case for more reusable libraries
• Atomic electronic structure calculations with Hermite interpolating polynomials
• 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