Using the quantum Monte Carlo simulation technique for the lattice fermion systems, we study the charge susceptibility of the two-dimensional Hubbard near the half-filled band. Lattice sizes up to 6×6 and inverse temperatures up to 6 t ( t : hopping parameter) are studied. We find that the behavior of the charge susceptibility is strongly dependent on the electron concentration at low temperature: the charge susceptibility is suppressed at the half-filled band, while it is considerably enhanced in non-half-filled band cases by the electron correlation. The above behavior of the charge susceptibility is quite similar to the recent Bethe-Ansatz result for the one-dimensional Hubbard model.
A concise general expression for the dynamical charge susceptibility of the Hubbard model is presented. The results are obtained from the spin-rotation-invariant slave-boson theory, whose paramagnetic saddle point corresponds to Gutzwiller's variational approximation. The charge structure factor and charge susceptibility are calculated for a wide range of interactions and temperatures and are compared thoroughly with the available quantum Monte Carlo and variational Monte Carlo results for the one-dimensional Hubbard model. The comparison shows reasonable agreement. All the calculations are done in a totally self-consistent way, without adjustable parameters.