An assessment of initial guesses for self-consistent field calculations. Superposition of Atomic Potentials: simple yet efficient

Abstract

Electronic structure calculations, such as in the Hartree-Fock or Kohn-Sham density functional approach, require an initial guess for the molecular orbitals. The quality of the initial guess has a significant impact on the speed of convergence of the self-consistent field (SCF) procedure. Popular choices for the initial guess include the one-electron guess from the core Hamiltonian, the extended Hückel method, and the superposition of atomic densities (SAD). Here, we discuss alternative guesses obtained from the superposition of atomic potentials (SAP), which is easily implementable even in real-space calculations. We also discuss a variant of SAD which produces guess orbitals by purification of the density matrix that could also be used in real-space calculations, as well as a parameter-free variant of the extended Hückel method, which resembles the SAP method and is easy to implement on top of existing SAD infrastructure. The performance of the core Hamiltonian, the SAD and the SAP guesses as well as the extended Hückel variant is assessed in non-relativistic calculations on a dataset of 259 molecules ranging from the first to the fourth periods by projecting the guess orbitals onto precomputed, converged SCF solutions in single- to triple-ζ basis sets. It is shown that the proposed SAP guess is the best guess on average. The extended Hückel guess offers a good alternative, with less scatter in accuracy.

13 pages, 2 figures

References in corpus (7)

• Localized atomic basis set in the projector augmented wave method
• All-electron density functional theory and time-dependent density functional theory with high-order finite elements
• Fully numerical Hartree-Fock and density functional calculations. I. Atoms
• The Python-based Simulations of Chemistry Framework (PySCF)
• Simple and accurate exchange energy for density functional theory
• Prediction of many-electron wavefunctions using atomic potentials
• Prediction of many-electron wavefunctions using atomic potentials: refinements and extensions to transition metals and large systems

Cited by in corpus (33)

• 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
• 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
• Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism
• Communication: Curing basis set overcompleteness with pivoted Cholesky decompositions
• 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
• A call to arms: making the case for more reusable libraries
• Iterative subspace algorithms for finite-temperature solution of Dyson equation
• Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules
• Polarized Gaussian basis sets from one-electron ions
• On the accurate reproduction of strongly repulsive interatomic potentials
• Energy-Based Molecular Orbital Localization in a Specific Spatial Region
• Fermi-Löwdin orbital self-interaction correction using the optimized effective potential method within the Krieger-Li-Iafrate approximation
• SPA${}^{\mathrm{H}}$M: the Spectrum of Approximated Hamiltonian Matrices representations
• A black-box, general purpose quadratic self-consistent field code with and without Cholesky Decomposition of the two-electron integrals
• Repulsive interatomic potentials calculated at three levels of theory
• Understanding electron correlation energy through density functional theory
• Predicting The One-Particle Density Matrix With Machine Learning
• Simple and accurate exchange energy for density functional theory
• How good are recent density functionals for ground and excited states of one-electron systems?
• Economical Quasi-Newton Self Consistent Field Solver
• SPA${}^{\mathrm{H}}$M(a,b): encoding the density information from guess Hamiltonian in quantum machine learning representations
• 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
• Atomic effective potentials for starting molecular electronic structure calculations
• Prediction of many-electron wavefunctions using atomic potentials: extended basis sets and molecular dissociation
• Complexity Reduction in Density Functional Theory: Locality in Space and Energy
• Review of the finite difference Hartree-Fock method for atoms and diatomic molecules, and its implementation in the x2dhf program
• A Reusable Library for Second-Order Orbital Optimization Using the Trust Region Method
• Importance profiles. Visualization of atomic basis set requirements