We introduce a new numerical method for the solution of self-consistent equations in the cluster mean-field theory. The method uses the density matrix renormalization group method to solve the associated cluster problem. We obtain an accurate critical value of the supersolid-superfluid transitions in the hard-core bosonic Hubbard model on a triangular lattice, which is comparable with the recent quantum Monte Carlo results. This algorithm is applicable to more general classes of models with a larger number of degrees of freedom.
The bosonic t-J model is a strong-on-site repulsion limit of the two-component Bose-Hubbard model and is expected to be realized by experiments of cold atoms in an optical lattice. In previous papers, we studied the bosonic t-J model by both analytical methods and numerical Monte - Carlo (MC) simulations. However, in the case of finite
Jz, where
Jz is the
z-component coupling constant of the pseudospin interaction, the phase diagram of the model was investigated by assuming the checkerboard type of boson densities. In this study, we shall continue our previous study of the bosonic t-J model using both the Gross-Pitaevskii (GP) theory and MC simulations without assuming any pattern of boson densities. These two methods complement each other and give reliable results. We show that as
Jz is increased, the superfluid state evolves into a supersolid (SS), and furthermore into a genuine solid with the checkerboard symmetry. In the present study, we propose a method identifying quantum phase transitions in the GP theory. We also study finite-temperature phase transitions of the superfluidity and the diagonal solid order of the SS by MC simulations.