Multiple bonds and excited states from the Hartree-Fock-Heitler-London method

The recently proposed Hartree-Fock-Heitler-London, HF-HL, method (Corongiu, G. J. Phys. Chem. A 2006, 110, 11584) previously tested for single bond molecules is validated by potential energy computations for open and closed shells, single and multiple bonds, in ground and excited states of homopolar diatomic molecules of the first and second period. The simple HF-HL function, including the configurations for 2s/2p near degeneracy and avoiding state crossing, yields correct dissociation products, qualitatively correct binding, and accounts for non-dynamical correlation. Addition of ionic structures improves the ab initio HF-HL function and yields about 95% of the experimental binding energies on average. Computed excitation energies are also in agreement with laboratory values as verified for the 3 Pi u and 3 Zeta g- excited states of the C2 molecule. Computation of the remaining dynamical correlation using a semiempirical functional yields binding energies with an average deviation of 1.5 kcal/mol from laboratory values, and total energies with an average deviation of 0.7 kcal/mol from exact nonrelativistic dissociation energies.