The improved generator coordinate Hartree± Fock (GCHF) method is extended to molecular systems. The Grin±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 H 2 , Li 2 , 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 eventempered basis sets. For H 2 , Li 2 , and LiH, the dierences between our best energies and the corresponding numerical HF results are about 2´10 )2 , 1, and 4´10 )1 lhartree, respectively. The dipole, quadrupole, and octupole moments at the center of mass and electric ®eld, the electric ®eld gradient, the electrostatic potential, and the electron density at the nuclei were evaluated and compared with results reported in the literature.