Abstract
We investigate the two-dimensional Hubbard model by the variational Monte Carlo method. We consider off-diagonal wave functions where we improve the Gutzwiller wave function by multiplying off-diagonal factors. In the 4×4 system, we compare our results with the exact diagonalization. The variational energies are fairly close to the exact values. We show the results for the 10×10 Hubbard model.