  Using recent advances in auxiliary-field quantum Monte Carlo techniques and
the phaseless approximation to control the sign/phase problem, we determine the
equation of state in the ground state of the two-dimensional repulsive
single-band Hubbard model at intermediate interactions. Shell effects are
eliminated and finite-size effects are greatly reduced by boundary condition
integration. Spin-spin correlation functions and structure factors are also
calculated. In lattice sizes up to $16\times 16$, the results show signal for
phase-separation. Upon doping, the system separates into one phase of density
$n=1$ (hole-free) and the other at density $n_c$ ($\sim 0.9$). The long-range
antiferromagnetic order is coupled to this process, and is lost below $n_c$.
