Prediction of many-electron wavefunctions using atomic potentials: Refinements and extensions to transition metals and large systems

Prediction of many-electron wavefunctions using atomic potentials: extended basis sets and molecular dissociation

A one-electron Schrödinger equation based on special one-electron potentials for atoms is shown to exist that produces orbitals for an arbitrary molecule that are sufficiently accurate to be used without modification to construct single- and multi-determinant wavefunctions. The exact Hamiltonian is used to calculate the energy variationally and to generate configuration interaction expansions. Earlier work on equilibrium geometries is extended to larger basis sets and molecular dissociation. For a test set of molecules representing different bonding environments, a single set of invariant atomic potentials gives wavefunctions with energies that deviate from configuration interaction energies based on SCF orbitals by less than 0.04 eV per bond or valence electron pair. On a single diagonalization of the Fock matrix, the corresponding errors are reduced 0.01 eV. Atomization energies are also in good agreement with CI values based on canonical SCF orbitals. Configuration interaction applications to single bond dissociations of water and glycine, and multiple bond dissociations of ethylene and oxygen produce dissociation energy curves in close agreement with CI calculations based on canonical SCF orbitals for the entire range of internuclear distances.

Assessment of Initial Guesses for Self-Consistent Field Calculations. Superposition of Atomic Potentials: Simple yet Efficient

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 nonrelativistic calculations on a data set 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.