In a recent Letter  (see also ) the authors presented numerical evidence supporting an idea of a direct transition between the superfluid (SF) and Mott insulating (MI) phases in the disordered Bosonic system, and even studied non-trivial properties of the multicritical line where SF, MI and the Bose Glass (BG) phases meet. The results were obtained from Monte Carlo simulations of the (2+1)-dimensional classical loop-current model  with the lattice action S = 1 2K ÞE ~ J=0 XrƒÑ ~ J2(r, ƒÑ) . 2(ƒÊ + v(r)) ~ JƒÑ (r, ƒÑ) . (1) where r, ƒÑ are spatial and imaginary time coordinates, and ~ J(r, ƒÑ) are integer current vectors with zero divergence. The spatial disorder potential v(r) is uniformly distributed on the interval (.,).
Prokof’ev, N and Svistunov, B, "Comment on “Phase Diagram of a Disordered Boson Hubbard Model in Two Dimensions”" (2003). Physics Department Faculty Publication Series. 1121.