Coulomb-hole-Hartree-Fock functional for molecular systems

We have extended to molecules a density functional previously parametrized for atomic computations. The Coulombhole-Hartree-Fock functional, introduced by Clementi in 1963, estimated the dynamic correlation energy by the computation of a Hartree-Fock type single-determinant wavefunction, where the Hartree-Fock potential was augmented with an effective potential term, related to a hard Coulomb hole enclosing each electron. The method was later revised by S. Chakravorty and E. Clementi, Phys. Rev. A, 38 (1989) 2290, so that a Yukawa-type soft Coulomb hole replaced the previous hard hole. Atomic correlation energies, computed for atoms with Z = 2 to 54, as well as for a number of excited states, validated the method. In this work we have parametrized for molecules a function which controls the width of the soft Coulomb hole by fitting the first and second atomic ionization potentials of atoms with 1 Z 18 and the binding energies of a few diatomic molecules. The parametrization was successfully validated by computing the dissociation energy for a number of molecules. A few-determinant version of the Coulomb-Hartree-Fock method (CHF-N) necessary to account for the non-dynamic correlation correction and to ensure proper dissociation products, is briefly discussed with reference to a previous proposal by G.C. Lie and E. Clementi, J. Chem. Phys., 60 (1974) 1275 and 60 (1974) 1288.