Motivated by some inconsistencies in the way quantum fluctuations are included beyond the classical treatment of hard-core bosons on a lattice in the recent literature, we revisit the large-S semi-classical approach to hard-core bosons on the square lattice at T=0. First of all, we show that, if one stays at the purely harmonic level, the only correct way to get the 1/S correction to the density is to extract it from the derivative of the ground state energy with respect to the chemical potential, and that to extract it from a calculation of the ground state expectation value of the particle number operator, it is necessary to include 1/\sqrt{S} corrections to the harmonic ground state. Building on this alternative approach to get 1/S corrections, we provide the first semiclassical derivation of the momentum distribution, and we revisit the calculation of the condensate density. The results of these as well as other physically relevant quantities such as the superfluid density are systematically compared to quantum Monte Carlo simulations. This comparison shows that the logarithmic corrections in the dilute Bose gas limit are only captured by the semi-classical approach if the 1/S corrections are properly calculated, and that the semi-classical approach is able to reproduce the 1/k divergence of the momentum distribution at k=0. Finally, the effect of 1/S^2 corrections is briefly discussed.
Motivated by some inconsistencies in the way quantum fluctuations are included beyond the classical treatment of hard-core bosons on a lattice in the recent literature, we revisit the large-S semi-classical approach to hard-core bosons on the square lattice at T=0. First of all, we show that, if one stays at the purely harmonic level, the only correct way to get the 1/S correction to the density is to extract it from the derivative of the ground state energy with respect to the chemical potential, and that to extract it from a calculation of the ground state expectation value of the particle number operator, it is necessary to include 1/\sqrt{S} corrections to the harmonic ground state. Building on this alternative approach to get 1/S corrections, we provide the first semiclassical derivation of the momentum distribution, and we revisit the calculation of the condensate density. The results of these as well as other physically relevant quantities such as the superfluid density are systematically compared to quantum Monte Carlo simulations. This comparison shows that the logarithmic corrections in the dilute Bose gas limit are only captured by the semi-classical approach if the 1/S corrections are properly calculated, and that the semi-classical approach is able to reproduce the 1/k divergence of the momentum distribution at k=0. Finally, the effect of 1/S^2 corrections is briefly discussed.
Motivated by some inconsistencies in the way quantum fluctuations are included beyond the classical treatment of hard-core bosons on a lattice in the recent literature, we revisit the large-S semiclassical approach to hard-core bosons on the square lattice at T=0. First of all, we show that, if one stays at the purely harmonic level, the only correct way to get the 1/S correction to the density is to extract it from the derivative of the ground-state energy with respect to the chemical potential, and that to extract it from a calculation of the ground-state expectation value of the particle number operator, it is necessary to include 1/S corrections to the harmonic ground state. Building on this alternative approach to get 1/S corrections, we provide the first semiclassical derivation of the momentum distribution, and we revisit the calculation of the condensate density. The results of these as well as other physically relevant quantities such as the superfluid density are systematically compared to quantum Monte Carlo simulations. This comparison shows that the logarithmic corrections in the dilute Bose gas limit are only captured by the semiclassical approach if the 1/S corrections are properly calculated, and that the semiclassical approach is able to reproduce the 1/k divergence of the momentum distribution at k=0. Finally, the effect of 1/S2 corrections is briefly discussed.
In this work,we theoretically study hard-core bosons on a two-dimensional square optical superlattice at T = 0.First of all,we present the mean field phase diagram of this model in terms of the chemical potential μ and the alternating potential strength △.Besides a superfluid(SF) phase at △ = 0 and a charge density wave(CDW)phase in the large △ at half filling,we demonstrate that a supersolid(SS) phase emerges in the moderate △.Then,we focus on the μ = 0,e.g.,half filling case,using large-S semi-classical spin-wave approximation to study the SS to CDW quantum phase transition.In particular,we calculate the ground-state energy and the superfluid density at the level of1/S correction.We then compare the spin-wave results with the large scale quantum Monte Carlo(QMC) simulations using the cluster stochastic series expansion(CSSE) algorithm,and find that while the spin wave method is intuitive with clear physical pictures,the quantum critical point is quite different from that of numerical results which is believed to be accurate.We suggest that as simple as it is,this model still exhibits strong quantum fluctuations near the quantum critical point beyond the power of semiclassical spin-wave approach.
Abstract In this work, we theoretically study hard-core bosons on a two-dimensional square optical superlattice at T = 0. First of all, we present the mean field phase diagram of this model in terms of the chemical potential μ and the alternating potential strength Δ. Besides a superfluid (SF) phase at Δ = 0 and a charge density wave (CDW) phase in the large Δ at half filling, we demonstrate that a supersolid (SS) phase emerges in the moderate Δ. Then, we focus on the μ = 0, e.g., half filling case, using large-S semiclassical spin-wave approximation to study the SS to CDW quantum phase transition. In particular, we calculate the ground-state energy and the superfluid density at the level of 1/S correction. We then compare the spin-wave results with the large scale quantum Monte Carlo (QMC) simulations using the cluster stochastic series expansion (CSSE) algorithm, and find that while the spin wave method is intuitive with clear physical pictures, the quantum critical point is quite different from that of numerical results which is believed to be accurate. We suggest that as simple as it is, this model still exhibits strong quantum fluctuations near the quantum critical point beyond the power of semiclassical spin-wave approach.
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.
The one-particle density matrices for hard-core bosons in a one-dimensional harmonic trap are computed numerically for systems with up to 160 bosons. Diagonalization of the density matrix shows that the many-body ground state is not Bose-Einstein condensed. The ground-state occupation, the amplitude of the lowest natural orbital, and the zero momentum peak height scale as powers of the particle number, and the corresponding exponents are related to each other. Close to its diagonal, the density matrix for hard-core bosons is similar to that of noninteracting fermions.
We evaluate the zero temperature properties and Bose-Einstein condensation (BEC) of hard sphere bosons in a harmonic trap using Quantum Monte Carlo methods. The chief purpose of this work is to go beyond the dilute limit, to test the limits of the Gross-Pitaevskii equation and related mean field approximations and to explore the role of interactions in determining the zero temperature properties of the trapped Bose gas. The density is increased by adjusting both the number of trapped bosons, N , and the s-wave scattering length, a, to study systems from the highly dilute region corresponding to early experiments in trapped BEC [2, 3] up to liquid He densities and beyond. Rather than relying on the perturbative Bogoliubov approximation (which assumes a large condensate fraction) to describe the condensate, we use the one-body density matrix (OBDM) formulation of BEC so that the properties of the condensate for systems with arbitrarily large interactions may be studied. In this formulation of BEC, condensate properties are obtained by diagonalizing the OBDM and obtaining the corresponding single particle “natural orbitals” and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single particle orbital(s) with macroscopic occupation (N0 >> 1). Within this framework, we calculate the effects of interactions and increased density on the ground state energy, the density profile, the momentum distribution, the condensate fraction and condensate “wave-function” and several other properties.