We numerically study the ground state phase diagram of the two-dimensional hard-core Bose-Hubbard model with nearest-(V(1)) and next-nearest-neighbour (V(2)) repulsions. In particular, we focus on the quarter-filled phases where one supersolid and two solid phases are observed. Using both canonical and grand canonical quantum Monte Carlo (QMC) methods and a mean-field calculation, we provide evidence for the existence of a commensurate supersolid. Despite the two possible diagonal long-range orderings for the solid phase, only one kind of supersolid phase is found to be energetically stable. The competition between the two solid phases manifests itself as a first-order phase transition around 2V(2) ∼ V(1). The change of order parameters as a function of the chemical potential is also presented.
We numerically study the ground state phase diagram of the two-dimensional hard-core Bose-Hubbard model with nearest-(V(1)) and next-nearest-neighbour (V(2)) repulsions. In particular, we focus on the quarter-filled phases where one supersolid and two solid phases are observed. Using both canonical and grand canonical quantum Monte Carlo (QMC) methods and a mean-field calculation, we provide evidence for the existence of a commensurate supersolid. Despite the two possible diagonal long-range orderings for the solid phase, only one kind of supersolid phase is found to be energetically stable. The competition between the two solid phases manifests itself as a first-order phase transition around 2V(2) ∼ V(1). The change of order parameters as a function of the chemical potential is also presented.
We present a comprehensive numerical study on the ground state phase diagram of the two-dimensional hard core bosonic extended Hubbard model with nearest (V1) and next nearest neighbor (V2) repulsions. In addition to the quantum solid and superfluid phases, we report the existence of a striped supersolid and three-quarter (quarter)-filled supersolid at commensurate density $\ensuremath{\rho}=0.75$ (0.25) due to the interplay of V1 and V2 interactions. The nature of the three-quarter-filled supersolid and the associated quantum solid will be discussed. The observed quantum phase transition between the two supersolids of different symmetries is found to be first order.
We study the ground-state phase diagram of the hard-core extended boson Hubbard model on the square lattice with both nearest- and next-nearest-neighbor hoppings and repulsions, using Gutzwiller mean-field theory and quantum Monte Carlo simulations. We observe the formation of supersolid states with checkerboard, striped, and quarter-filled crystal structures when the system is doped away from commensurate fillings. In the striped supersolid phase, a strong anisotropy in the superfluid density is obtained from the simulations; however, the transverse component remains finite, indicating a true two-dimensional superflow. We find that upon doping, the striped supersolid transitions directly into the supersolid with quarter-filled crystal structure, via a first-order stripe melting transition.
Using an unbiased quantum Monte Carlo method, we obtain convincing evidence of the existence of a checkerboard supersolid at a commensurate filling factor 1/2 (a commensurate supersolid) in the soft-core Bose-Hubbard model with nearest-neighbor repulsions on a cubic lattice. In conventional cases, supersolids are realized at incommensurate filling factors by a doped-defect-condensation mechanism, where particles (holes) doped into a perfect crystal act as interstitials (vacancies) and delocalize in the crystal order. However, in the model, a supersolid state is stabilized even at the commensurate filling factor 1/2 without doping. By performing grand canonical simulations, we obtain a ground-state phase diagram that suggests the existence of a supersolid at a commensurate filling. To obtain direct evidence of the commensurate supersolid, we next perform simulations in canonical ensembles at a particle density ρ=1/2 and exclude the possibility of phase separation. From the obtained snapshots, we discuss its microscopic structure and observe that interstitial-vacancy pairs are unbound in the crystal order.
We study the extended Bose-Hubbard model on a two-dimensional honeycomb lattice by using large-scale quantum Monte Carlo simulations. We present the ground-state phase diagrams for both the hard-core and the soft-core bosons. For the hard-core case, the transition between the $\ensuremath{\rho}=1∕2$ solid and the superfluid is first order, and the supersolid state is unstable toward phase separation. For the soft-core case, due to the presence of multiple occupation, a stable particle-induced supersolid (SS-p) phase emerges when $1∕2<\ensuremath{\rho}<1$. The transition from the solid at $\ensuremath{\rho}=1∕2$ to the SS-p phase is second order with the superfluid density scaling as ${\ensuremath{\rho}}_{s}\ensuremath{\sim}\ensuremath{\rho}\ensuremath{-}1∕2$. The SS-p phase has the same diagonal order as the solid at $\ensuremath{\rho}=1∕2$. As the chemical potential increases further, the SS-p phase turns into a solid where two bosons occupy each site of one sublattice through a first-order transition. We also calculate the critical exponents of the transition between the $\ensuremath{\rho}=1∕2$ solid and superfluid at the Heisenberg point for the hard-core case. We find the dynamical critical exponent z=0.15, which is smaller than results obtained on smaller lattices. This indicates that z approaches zero in the thermodynamic limit, and thus the transition is also first order even at the Heisenberg point.
We study the soft-core extended Bose-Hubbard model in a two-dimensional triangular lattice by using the quantum Monte Carlo methods. The ground-state phase diagram of the system exhibits a very fruitful structure and is much intriguing. Aside from the Mott insulating state, four kinds of solid states with respect to the commensurate filling factors $\ensuremath{\rho}=1∕3,2∕3$ and $\ensuremath{\rho}=1$ are identified. In incommensurate fillings, superfluid, supersolid, as well as phase-separation states are detected. As in the case of the hard-core bosons, a supersolid phase exists in $1∕3<\ensuremath{\rho}<2∕3$ while it is unstable toward the phase separation in $\ensuremath{\rho}<1∕3$. However, this instability is refrained in $2∕3<\ensuremath{\rho}<1$ due to the softening of the bosons, and then a supersolid phase survives.
We investigate the ground-state phase diagram of the soft-core Bose-Hubbard model with the nearest-neighbor repulsion on a square lattice by using an unbiased quantum Monte Carlo method. In contrast to a previous study [P. Sengupta et al., Phys. Rev. Lett. 94, 207202 (2005)], we present the ground-state phase diagrams up to high hopping amplitudes. As a result, in addition to the known supersolid above half-filling, we find a supersolid phase below and at half-filling for high hopping amplitudes. In addition, for a strong nearest-neighbor repulsion, we show that the supersolid phase occupies a remarkably broad region in the phase diagram. These results are in agreement with the results of the Gutzwiller mean-field approximation [M. Iskin, Phys. Rev. A 83, 051606(R) (2011); T. Kimura, Phys. Rev. A 84, 063630 (2011)]. However, it turns out that the regions of the supersolid phases are significantly smaller than the mean-field results.