Dear prof. Matthias
Thanks for your reply.
Is there a different number of up and down atoms or is just the
distribution not symmetric? There are both a different number of up and down atoms and a asymmetric distribution. When I increase SWEEPS and MAXSTATES from 30,300 to 50,500, there are slight in the result (< <10^-3).
Restricting Sz=0 on each site means that no singly occupied sites will
be allowed. This is probably not what you want. For a simple Hubbard model with external potential, the operators for up and down atoms in Hamiltonian are symmetry. There will not be singly occupied sites in the result.
Thanks again.
Chen Ahai