Class of variational singlet wave functions for the Hubbard model away from half filling

We present a class of variational wave functions for strong-coupling Heisenberg Hubbard models. These are written in the form of three factors---a pair of determinants and a Jastrow function---and are made out of orbitals, a la Hartree-Fock theory, which solve a fictitious one-body problem. The wave functions respect various constraints known from general principles and appear to be potentially useful in understanding the possible behavior of the models in quantitative terms.