We investigate the phase diagram of interacting bosons in a tripartite lattice in two dimensions. An analytic computation of the phase diagram in the parameter space defined by the on-site boson–boson repulsion parameter and the chemical potential is done via a second-order strong coupling perturbation theory on the extended Bose–Hubbard model, which is hence supplemented by numerical mean-field calculations. Interesting results emerge in the form of exotic density ordered phases with one-third, two-third and unity filling of bosons per trimer. The pattern repeats as the chemical potential is increased. Further regions of supersolid and superfluid phases are noted in the numerical study of the mean-field model where the phase diagram shows a transition from compressible to insulating phases.
We investigate the phase diagram of interacting bosons in a tripartite lattice in two dimensions. An analytic computation of the phase diagram in the parameter space defined by the on-site boson–boson repulsion parameter and the chemical potential is done via a second-order strong coupling perturbation theory on the extended Bose–Hubbard model, which is hence supplemented by numerical mean-field calculations. Interesting results emerge in the form of exotic density ordered phases with one-third, two-third and unity filling of bosons per trimer. The pattern repeats as the chemical potential is increased. Further regions of supersolid and superfluid phases are noted in the numerical study of the mean-field model where the phase diagram shows a transition from compressible to insulating phases.
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.
Motivated by the possibility of obtaining a trimerized optical lattice, such as a kagome lattice experimentally by employing counter propagating laser beams, we investigate the phase diagram of correlated bosons in a kagome lattice in the presence of a harmonic confinement. To study the phase diagram, we compute the order parameters of an extended Bose–Hubbard model which is solved via the standard mean field technique in the presence of a trapping potential. It is emphasized that such trapping effects (together with the density exchange term) lead to the necessity of considering the detailed geometry of the lattice, and thus the measured quantities appeal to a particular geometry, kagome lattice being the case here. The phase diagram thus obtained yields co-existence of different compressible and incompressible phases. There are distinct signatures of superfluid, insulating, density ordered and supersolid phases as one moves towards the edge of the lattice starting from the centre of the confining potential. A comparison of the results obtained for the order parameters corresponding to other geometries, such as square and honeycomb lattices are presented. Relevant to the time-of-flight experiments, the order parameters have been computed in the momentum space, which corroborate the existence of different phases for a kagome lattice.
We analyse the ground-state quantum phase diagram of hardcore Bosons interacting with repulsive dipolar potentials. The bosons dynamics is described by the extended-Bose-Hubbard Hamiltonian on a two-dimensional lattice. The ground state results from the interplay between the lattice geometry and the long-range interactions, which we account for by means of a classical spin mean-field approach limited by the size of the considered unit cells. This extended classical spin mean-field theory accounts for the long-range density-density interaction without truncation. We consider three different lattice geometries: square, honeycomb, and triangular. In the limit of zero hopping the ground state is always a devil’s staircase of solid (gapped) phases. Such crystalline phases with broken translational symmetry are robust with respect to finite hopping amplitudes. At intermediate hopping amplitudes, these gapped phases melt, giving rise to various lattice supersolid phases, which can have exotic features with multiple sublattice densities. At sufficiently large hoppings the ground state is a superfluid. The stability of phases predicted by our approach is gauged by comparison to the known quantum phase diagrams of the Bose-Hubbard model with nearest-neighbour interactions as well as quantum Monte Carlo simulations for the dipolar case on the square and triangular lattice. Our results are of immediate relevance for experimental realisations of self-organised crystalline ordering patterns in analogue quantum simulators, e.g., with ultracold dipolar atoms in an optical lattice.
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 ground-state properties of repulsively interacting bosons on the honeycomb lattice using large-scale quantum Monte Carlo simulations. In the hard-core limit, the half-filled system develops long-ranged diagonal order for sufficiently strong nearest-neighbor repulsion. This staggered solid melts at a first-order quantum phase transition into the superfluid phase, without the presence of any intermediate supersolid phase. Within the superfluid phase, both the superfluid density and the compressibility exhibit local minima near the particle (hole) density of one quarter, while the density and the condensate fraction show inflection points in this region. Relaxing the hard-core constraint, supersolid phases emerge for soft-core bosons. The suppression of the superfluid density is found to persist for sufficiently large, finite on-site repulsion.
We carry out an extensive study of the phase diagrams of the extended Bose Hubbard model, with a mean filling of one boson per site, in one dimension by using the density matrix renormalization group and show that it contains Superfluid (SF), Mott-insulator (MI), density-wave (DW) and Haldane-insulator (HI) phases. We show that the critical exponents and central charges for the HI-DW, MI-HI and SF-MI transitions are consistent with those for models in the two-dimensional Ising, Gaussian, and Berezinskii-Kosterlitz-Thouless (BKT) universality classes, respectively; and we suggest that the SF-HI transition may be more exotic than a simple BKT transition. We show explicitly that different boundary conditions lead to different phase diagrams.