Atomic Confinement Potentials and the Generation of Numerical Atomic Orbitals

Abstract

We aim to develop novel reusable open source infrastructure [Lehtola, J. Chem. Phys. 159, 180901 (2023)] for numerical atomic orbitals (NAOs). Soft confinement potentials are typically used to force the NAO radial basis functions ${\psi }_{nl}(r)$ to vanish smoothly in increasing $r$ and to generate localized unoccupied states; we review such potentials and other commonly-used techniques in NAO generation as a follow-up to our recent study on atoms in hard-wall confinement [Åström and Lehtola, J. Phys. Chem. A 129, 2791 (2025)]. In addition to NAO generation, confinement potentials are also employed to simulate environmental effects in other research areas, such as studies of (i) atoms in solids, (ii) quantum dots, and (iii) high-pressure chemistry. As in our earlier work, we perform fully numerical density functional calculations with spherically averaged densities, as is usual in NAO studies. Our calculations employ the the finite element method (FEM) implemented in the HelFEM program, yielding variational energies and enabling the use of various boundary conditions. We consider four families of potentials to study the Mg and Ca atoms, which are textbook examples of extended electronic structures. We show that the resulting ground-state orbitals are surprisingly insensitive to the employed form of the confinement potential, and that the orbitals decay quickly under confinement. We study increasingly steep potentials and examine how they approach the hard-wall limit. Finally, we assess NAO basis set truncation errors for types of singular potentials that are now broadly used in the NAO literature.

25 pages, 19 figures. Fixed minor issues in a few figures

References in corpus (37)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Numerical atomic orbitals for linear scaling
• Linear-scaling ab-initio calculations for large and complex systems
• Localized atomic basis set in the projector augmented wave method
• Systematic generation of finite-range atomic basis sets for linear-scaling calculations
• GPAW: An open Python package for electronic-structure calculations
• Systematically improvable optimized atomic basis sets for {\it ab inito} calculations
• The Elephant in the Room of Density Functional Theory Calculations
• Large scale and linear scaling DFT with the CONQUEST code
• 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
• Nanoplasmonics simulations at the basis set limit through completeness-optimized, local numerical basis sets
• Fully numerical calculations on atoms with fractional occupations. Range-separated exchange functionals
• Fully numerical Hartree-Fock and density functional calculations. I. Atoms
• Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms
• Optimal finite-range atomic basis sets for liquid water and ice
• Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism
• Fully numerical electronic structure calculations on diatomic molecules in weak to strong magnetic fields
• A call to arms: making the case for more reusable libraries
• Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets
• Atomic electronic structure calculations with Hermite interpolating polynomials
• Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules
• Energy gaps and interaction blockade in confined quantum systems
• On the accurate reproduction of strongly repulsive interatomic potentials
• DFTB modelling of lithium intercalated graphite with machine-learned repulsive potential
• Full correlation single-particle positron potentials for a positron and a positronium interacting with atoms
• A simple Dirac wave function for a Coulomb potential with linear confinement
• First-principles density-functional calculations using localized spherical-wave basis sets
• Simple and efficient LCAO basis sets for the diffuse states in carbon nanostructures
• Meta-local density functionals: a new rung on Jacob's ladder
• Insight on Gaussian basis set truncation errors in weak to intermediate magnetic fields with an approximate Hamiltonian
• OpenOrbitalOptimizer -- a reusable open source library for self-consistent field calculations
• On the accuracy of a recent regularized nuclear potential
• Multiverse cosmological models
• Systematic study of confinement induced effects on atomic electronic structure
• Importance profiles. Visualization of atomic basis set requirements