Abstract: The work discussed algorithm of grand canonical quantum Monte Carlo method for the one-band two-dimensional Hubbard model, and calculated local magnetic moment, magnetic susceptibility, staggered magnetic susceptibility and the internal energy of the model. Results show that grand canonical quantum Monte Carlo method is very effective in dealing with Hubbard model from weak to strong correlation rigim, and also in dealing with many——body problems with strong correlation.