PHYSICAL REVIEW B, v. 72, No. 18, Article no. 184507
We consider the two-species bosonic Hubbard model with variable interspecies interaction and hopping strength in the grand canonical ensemble with a common chemical potential. We analyze the superfluid-insulator (SI) transition for the relevant parameter regimes and compute the ground state phase diagram in the vicinity of odd filling Mott states. We find that the superfluid-insulator transition occurs with (a) simultaneous onset of superfluidity of both species or (b) coexistence of Mott insulating state of one species and superfluidity of the other or, in the case of unit filling, (c) complete depopulation of one species. The superfluid-insulator transition can be first order in a large region of the phase diagram. We develop a variational mean-field method which takes into account the effect of second order quantum fluctuations on the superfluid-insulator transition and corroborate the mean-field phase diagram using a quantum Monte Carlo study.