Hard-core and soft-core extended Bose-Hubbard models with fixed boson densities are studied by using the Gutzwiller variational wave function. The phase diagram of the softcore model strongly depends on its transfer integral. Specifically, a supersolid phase can be the ground state even below half-filling for a large transfer integral. The density difference between nearest-neighbor sites, which is an order parameter of the supersolid phase, strongly depends on the boson density, especially for a small transfer integral.
Hard-core and soft-core extended Bose-Hubbard models with fixed boson densities are studied by using the Gutzwiller variational wave function. The phase diagram of the softcore model strongly depends on its transfer integral. Specifically, a supersolid phase can be the ground state even below half-filling for a large transfer integral. The density difference between nearest-neighbor sites, which is an order parameter of the supersolid phase, strongly depends on the boson density, especially for a small transfer integral.
We studied all possible ground states, including supersolid (SS) phases and phase separations of hard-core- and soft-core-extended Bose-Hubbard models with fixed boson densities by using the Gutzwiller variational wave function and the linear programming method. We found that the phase diagram of the soft-core model depends strongly on its transfer integral. Furthermore, for a large transfer integral, we showed that an SS phase can be the ground state even below or at half filling against the phase separation. We also found that the density difference between nearest-neighbor sites, which indicates the density order of the SS phase, depends strongly on the boson density and transfer integral.
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.
We use the Gutzwiller ansatz and analyze the phase diagram of the extended Bose-Hubbard Hamiltonian with on-site (U) and nearest-neighbor (V) repulsions. For d-dimensional hypercubic lattices, when 2dV<U, it is well known that the ground state alternates between the charge-density-wave (CDW) and Mott insulators, and the supersolid (SS) phase occupies small regions around the CDW insulators. However, when 2dV>U, in this Rapid Communication, we show that the ground state has only CDW insulators, and more importantly, the SS phase occupies a much larger region in the phase diagram, existing up to very large hopping values which could be orders of magnitude higher than that of the well-known case. In particular, the SS-superfluid phase boundary increases linearly as a function of hopping when $2\mathit{dV}\ensuremath{\gtrsim}1.5U$, for which the prospects of observing the SS phase with dipolar Bose gases loaded into optical lattices is much higher.
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.
By using the inhomogeneous mean-field theory, the phase diagram of the two-dimensional soft-core extended Bose-Hubbard model with both on-site and nearest-neighbor interactions in the presence of spatial disorder is determined. Rich phases, including Mott-insulator, checkerboard solid, supersolid, superfluid and disordered solid are obtained. Interestingly, due to the presence of disorder, a new disordered solid phase is found between the incompressible lobes, and it is demonstrated that the supersolid phase can survive in two-dimensional soft-core bosons. Furthermore, the system undergoes the phase transition from a disordered solid into a solid order or supersolid landscape when the disorder strength is increased. For weak nearest-neighbor repulsion, the increase of the disorder strength induces both checkerboard solid and supersolid phases shrinking to smaller hopping regions. However, for strong nearest-neighbor repulsion, the disorder stabilizes the supersolid phase and makes it occupying a remarkably broad region in the phase diagram. The analytical expressions for the phase boundaries between the incompressible and compressible phases are also obtained, which are in qualitative agreement with our numerical results.
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.