Phase Diagram of Commensurate Two-Dimensional Disordered Bose Hubbard Model

We do quantum Monte-Carlo simulations to calculate the full ground state phase diagram of the two dimensional disordered Bose-Hubbard model at unity filling factor. We observe superfluid regions persisting up to large values of disorder and interaction strength. In the vicinity of the superfluid-insulator transition of the pure system, the system shows almost unmeasurable response to weak disorder.