Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets

Abstract

The superposition of atomic potentials (SAP) approach has recently been shown to be a simple and efficient way to initialize electronic structure calculations [S. Lehtola, J. Chem. Theory Comput. 15, 1593 (2019)]. Here, we study the differences between effective potentials from fully numerical density functional and optimized effective potential calculations for fixed configurations. We find that the differences are small, overall, and choose exchange-only potentials at the local density approximation level of theory computed on top of Hartree-Fock densities as a good compromise. The differences between potentials arising from different atomic configurations are also found to be small at this level of theory. Furthermore, we discuss the efficient Gaussian-basis implementation of SAP via error function fits to fully numerical atomic radial potentials. The guess obtained from the fitted potentials can be easily implemented in any Gaussian-basis quantum chemistry code in terms of two-electron integrals. Fits covering the whole periodic table from H to Og are reported for non-relativistic as well as fully relativistic four-component calculations that have been carried out with fully numerical approaches.

arXiv comment12 pages, 8 figures

References in corpus (13)

• Ions in solution: Density Corrected Density Functional Theory (DC-DFT)
• 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
• How is the derivative discontinuity related to steps in the exact Kohn-Sham potential?
• Fully numerical calculations on atoms with fractional occupations. Range-separated exchange functionals
• Fully numerical Hartree-Fock and density functional calculations. I. Atoms
• Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules
• Polarized Gaussian basis sets from one-electron ions
• Atomic effective potentials for starting molecular electronic structure calculations
• Prediction of many-electron wavefunctions using atomic potentials
• Prediction of many-electron wavefunctions using atomic potentials: extended basis sets and molecular dissociation
• Prediction of many-electron wavefunctions using atomic potentials: refinements and extensions to transition metals and large systems

Cited by in corpus (13)

• The DIRAC code for relativistic molecular calculations
• 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
• Atomic electronic structure calculations with Hermite interpolating polynomials
• SPA${}^{\mathrm{H}}$M: the Spectrum of Approximated Hamiltonian Matrices representations
• 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
• Review of the finite difference Hartree-Fock method for atoms and diatomic molecules, and its implementation in the x2dhf program
• Atomic Confinement Potentials and the Generation of Numerical Atomic Orbitals
• A Reusable Library for Second-Order Orbital Optimization Using the Trust Region Method
• A new generation of effective core potentials: Selected lanthanides and heavy elements II
• Implementation of McMurchie-Davidson algorithm for Gaussian AO integrals suited for SIMD processors