We calculate the superfluid transition temperature of homogeneous interacting Bose gases in three and two spatial dimensions using large-scale path integral Monte Carlo simulations (with up to N=10;{5} particles). In 3D we investigate the limits of the universal critical behavior in terms of the scattering length alone by using different models for the interatomic potential. We find that this type of universality sets in at small values of the gas parameter na3 < or approximately 10(-4). This value is different from the estimate na3 < or approximately 10(-6) for the validity of the asymptotic expansion in the limit of vanishing na3. In 2D we study the Berezinskii-Kosterlitz-Thouless transition of a gas with hard-core interactions. For this system we find good agreement with the classical lattice |psi|4 model up to very large densities. We also explain the origin of the existing discrepancy between previous studies of the same problem.
We calculate the superfluid transition temperature of homogeneous interacting Bose gases in three and two spatial dimensions using large-scale path integral Monte Carlo simulations (with up to N=10;{5} particles). In 3D we investigate the limits of the universal critical behavior in terms of the scattering length alone by using different models for the interatomic potential. We find that this type of universality sets in at small values of the gas parameter na3 < or approximately 10(-4). This value is different from the estimate na3 < or approximately 10(-6) for the validity of the asymptotic expansion in the limit of vanishing na3. In 2D we study the Berezinskii-Kosterlitz-Thouless transition of a gas with hard-core interactions. For this system we find good agreement with the classical lattice |psi|4 model up to very large densities. We also explain the origin of the existing discrepancy between previous studies of the same problem.
We study the interacting homogeneous Bose gas in two spatial dimensions in the thermodynamic limit at fixed density. We shall be concerned with some mathematical aspects of this complicated problem in many-body quantum mechanics. More specifically, we consider the dilute limit where the scattering length of the interaction potential, which is a measure for the effective range of the potential, is small compared to the average distance between the particles. We are interested in a setting with positive (i.e., non-zero) temperature. After giving a survey of the relevant literature in the field, we provide some facts and examples to set expectations for the two-dimensional system. The crucial difference to the three-dimensional system is that there is no Bose–Einstein condensate at positive temperature due to the Hohenberg–Mermin–Wagner theorem. However, it turns out that an asymptotic formula for the free energy holds similarly to the three-dimensional case. We motivate this formula by considering a toy model with δ interaction potential. By restricting this model Hamiltonian to certain trial states with a quasi-condensate we obtain an upper bound for the free energy that still has the quasi-condensate fraction as a free parameter. When minimizing over the quasi-condensate fraction, we obtain the Berezinskii–Kosterlitz–Thouless critical temperature for superfluidity, which plays an important role in our rigorous contribution. The mathematically rigorous result that we prove concerns the specific free energy in the dilute limit. We give upper and lower bounds on the free energy in terms of the free energy of the non-interacting system and a correction term coming from the interaction. Both bounds match and thus we obtain the leading term of an asymptotic approximation in the dilute limit, provided the thermal wavelength of the particles is of the same order (or larger) than the average distance between the particles. The remarkable feature of this result is its generality: the correction term depends on the interaction potential only through its scattering length and it holds for all nonnegative interaction potentials with finite scattering length that are measurable. In particular, this allows to model an interaction of hard disks.
By using exact path-integral Monte Carlo methods we calculate the equation of state of an interacting Bose gas as a function of temperature both below and above the superfluid transition. The universal character of the equation of state for dilute systems and low temperatures is investigated by modeling the interatomic interactions using different repulsive potentials corresponding to the same s-wave scattering length. The results obtained for the energy and the pressure are compared to the virial expansion for temperatures larger than the critical temperature. At very low temperatures we find agreement with the ground-state energy calculated using the diffusion Monte Carlo method.
We study the superfluid transition of weakly interacting Bose gases, classical fields and related systems in two dimensions, where it is a Berezinskii-Kosterlitz-Thouless transition, and in three dimensions, where it is a U(1) transition accompanied by the onset of Bose-Einstein condensation. We show that for small interaction strength U, these systems can be described by a :Ψ:4 model based on classical complex fields and we establish universal behavior among systems of different microscopical details, such as classical and quantum, continuous and lattice. By numerical simulations of the classical lattice field, using Worm Algorithm, we solve all relevant thermodynamical functions of this system, which are density, condensate density (in 3D), superfluid density and quasi-condensate. Through universality considerations, we can translate these results to all systems which can be described by a :Ψ:4 model based on a classical complex field. Our results cover the fluctuation region and extend into mean field and critical region, where we connect them to results of mean field theory and renormalization group theory respectively. We also establish that the quasicondensate, which exists in two and three dimensions, is a basis for a more accurate mean field theory treatment than the condensate, which exists only in three dimensions.
Using renormalization-group arguments we show that the low-temperature thermodynamics of a three- or two-dimensional dilute Bose gas is fully determined by a universal scaling function ${\mathcal{F}}_{d}(\ensuremath{\mu}/{k}_{B}T,\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{g}(T))$ once the mass m and the s-wave scattering length ad of the bosons are known (d is the space dimension). Here $\ensuremath{\mu}$ and T denote the chemical potential and temperature of the gas, and the temperature-dependent dimensionless interaction constant $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{g}(T)$ is a function of $m{a}_{d}^{2}{k}_{B}T/{\ensuremath{\hbar}}^{2}$. We compute the scaling function F2 using a nonperturbative renormalization-group approach and find that both the $\ensuremath{\mu}/{k}_{B}T$ and $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{g}(T)$ dependencies are in very good agreement with recent experimental data obtained for a quasi-two-dimensional Bose gas with or without optical lattice. We also show that the nonperturbative renormalization-group estimate of the Berezinskii-Kosterlitz-Thouless transition temperature compares well with the result obtained from a quantum Monte Carlo simulation of an effective classical field theory.
Using a finite-temperature path integral Monte Carlo simulation (PIMC) method and finite-size scaling, we have investigated the interaction-induced shift of the phase-transition temperature for Bose-Einstein condensation of homogeneous weakly interacting Bose gases in three dimensions, which is given by a proposed analytical expression ${T}_{c}={T}_{c}^{0}{1+{c}_{1}a{n}^{1∕3}+[{c}_{2}^{\ensuremath{'}}\phantom{\rule{0.2em}{0ex}}\mathrm{ln}(a{n}^{1∕3})+{c}_{2}^{\ensuremath{''}}]{a}^{2}{n}^{2∕3}+\mathrm{O}({a}^{3}n)}$, where Tc0 is the critical temperature for an ideal gas, a is the s-wave scattering length, and n is the number density. We have used smaller number densities and more time slices than in the previous PIMC simulations [Gruter et al., Phys. Rev. Lett. 79, 3549 (1997)] in order to understand the difference in the value of the coefficient c1 between their results and the (apparently) other reliable results in the literature. Our results show that ${({T}_{c}\ensuremath{-}{T}_{c}^{0})∕{T}_{c}^{0}}∕(a{n}^{1∕3})$ depends strongly on the interaction strength an1∕3 while the previous PIMC results are considerably flatter and smaller than our results. We obtain ${c}_{1}=1.32\ifmmode\pm\else\textpm\fi{}0.14$, in agreement with results from recent Monte Carlo methods of three-dimensional O(2) scalar ${\ensuremath{\phi}}^{4}$ field theory and variational perturbation theory.
We study the universal critical behavior of two-dimensional (2D) lattice bosonic gases at the Berezinskii-Kosterlitz-Thouless (BKT) transition, which separates the low-temperature superfluid phase from the high-temperature normal phase. For this purpose, we perform quantum Monte Carlo simulations of the hard-core Bose-Hubbard (BH) model at zero chemical potential. We determine the critical temperature by using a matching method that relates finite-size data for the BH model with corresponding data computed in the classical XY model. In this approach, the neglected scaling corrections decay as inverse powers of the lattice size L, and not as powers of 1/lnL, as in more standard approaches, making the estimate of the critical temperature much more reliable. Then, we consider the BH model in the presence of a trapping harmonic potential, and we verify the universality of the trap-size dependence at the BKT critical point. This issue is relevant for experiments with quasi-2D trapped cold atoms.