We investigate augmenting conventional Gaussian-type one-electron orbital basis sets with two-electron functions that have a Gaussian dependence on the interelectronic distance. We observe substantial improvements in calculated correlation energies for helium and neon atoms and for the water molecule. A feature of our approach is that there is no nonlinear optimization of the two-electron basis function parameters at all.