  The two-dimensional Hubbard model with a bimodal distribution of on-site
interactions, P(U_i) = (1-f)\delta(U_i-U) + f\delta(U_i), is studied using a
finite temperature quantum Monte Carlo technique and dynamical mean-field
theory. We find long range antiferromagnetic order off half-filling is
stabilized by the disorder, due to localization of the dopants on the U=0
sites. Whereas in the clean model there is a single gap at n=1, for nonzero f
we find the compressibility and density of states exhibit gaps at two separate
fillings.
