Improved generator coordinate Hartree–Fock method for molecular systems: application to H2, Li2 and LiH

The improved generator coordinate Hartree–Fock (GCHF) method is extended to molecular systems. The Griffin–Hill–Wheeler–HF equations were solved by an integral discretization technique. The method is then implemented with the use of the GAMESS program and applied to the H2, Li2, and LiH molecules. For these molecules, sequences of basis sets of atom-centred Gaussian-type functions are employed to explore the accuracy achieved with our approach. For all systems studied, our ground-state HF total energies are better than those obtained with basis sets generated with the original GCHF method for molecules and larger even-tempered basis sets. For H2, Li2, and LiH, the differences between our best energies and the corresponding numerical HF results are about 2 × 10−2, 1, and 4 × 10−1 μhartree, respectively. The dipole, quadrupole, and octupole moments at the center of mass and electric field, the electric field gradient, the electrostatic potential, and the electron density at the nuclei were evaluated and compared with results reported in the literature.