Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review A

Bose-Einstein condensation in two dimensions: A quantum Monte Carlo study

Sean Pearson, Tao Pang, Changfeng Chen

Published Dec 1, 1998• 15 citations
Publisher

Abstract

A path-integral quantum Monte Carlo method is used to calculate finite temperature properties of up to 1000 hard-core bosons in a two-dimensional isotropic harmonic-oscillator potential. If the interatomic repulsions are sufficiently short range, an abrupt increase in the condensate fraction and a hump in the specific heat occur close to the critical temperature of ideal bosons. The critical temperature and the condensate fraction are in general lowered by an increase in the hard-core radius a. If a is decreased below a certain level, the condensate fraction becomes indistinguishable from the corresponding value of the ideal bosons. For up to 1000 particles, this occurs when ${\mathrm{ln}}^{\ensuremath{-}1}{(1/na}^{2})\ensuremath{\lesssim}0.1,$ where n is the average particle density.

DOI
10.1103/physreva.58.4811

Related papers

Bose-Einstein condensation in one-dimensional power-law traps: A path-integral Monte Carlo simulation

Sean Pearson, Tao Pang, Changfeng Chen

Bose-Einstein condensation is known to occur in one-dimensional power-law potentials $V(x)\ensuremath{\propto}|x{|}^{\ensuremath{\eta}}$ as a true phase transition for \ensuremath{\eta}2, but only if there are no interactions between the particles. We show, via a path-integral quantum Monte Carlo scheme, that for both \ensuremath{\eta}2 and \ensuremath{\eta}=2 the spatial distribution of finite numbers of hard-core bosons suddenly becomes bimodal below a certain temperature, with a central condensate of particles distributed according to the lowest single-particle eigenstate. At still lower temperatures, the hard-core interactions cause the single-particle ground state to become a poor description for the interacting gas. If \ensuremath{\eta}2, the energy per particle undergoes a sudden decrease near the same temperature at which the bimodal distribution appears. It is found that this drop in energy disappears if the range of the interaction potential is sufficiently large.

Physical Review A•Aug 1, 1998•13 citations

Variational Monte Carlo analysis of Bose–Einstein condensation in a two-dimensional trap

Zheng Rong-Jie, Jing Jin, Yi Ping Tang

The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas–Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.

Chinese Physics•Aug 22, 2006•6 citations

Path-Integral Monte Carlo and the Squeezed Trapped Bose-Einstein Gas

Juan Pablo Fernández, William J. Mullin

Bose-Einstein condensation has been experimentally found to take place in finite trapped systems when one of the confining frequencies is increased until the gas becomes effectively two-dimensional (2D). We confirm the plausibility of this result by performing path-integral Monte Carlo (PIMC) simulations of trapped Bose gases of increasing anisotropy and comparing them to the predictions of finite-temperature many-body theory. PIMC simulations provide an essentially exact description of these systems; they yield the density profile directly and provide two different estimates for the condensate fraction. For the ideal gas, we find that the PIMC column density of the squeezed gas corresponds quite accurately to that of the exact analytic solution and, moreover, is well mimicked by the density of a 2D gas at the same temperature; the two estimates for the condensate fraction bracket the exact result. For the interacting case, we find 2D Hartree-Fock solutions whose density profiles coincide quite well with the PIMC column densities and whose predictions for the condensate fraction are again bracketed by the PIMC estimates.

arXiv:cond-mat/0606710•AIP conference proceedings•Jan 1, 2006•1 citation

Bose-Einstein condensation in anisotropic harmonic traps

Tor Haugset, Hårek Haugerud, Jens O. Andersen

We study the thermodynamic behavior of an ideal gas of bosons trapped in a three-dimensional anisotropic harmonic-oscillator potential. The condensate fraction as well as the specific heat is calculated using the Euler-Maclaurin approximation. For a finite number of particles there is no phase transition, but there is a well-defined temperature at which the condensation starts. We also consider condensation in lower dimensions, and, for one-dimensional systems, we discuss the dependence of the condensate fraction and heat capacity on the ensemble used.

arXiv:cond-mat/9611061•Physical Review A•Apr 1, 1997•61 citations

BOSE-EINSTEIN CONDENSATION IN TRAPS: A QUANTUM MONTE CARLO STUDY

Jonathan L. DuBois

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.

•Jan 1, 2002•1 citation

Quantum Monte Carlo simulation of a two-dimensional Bose gas

Sebastiano Pilati, Jordi Boronat, J. Casulleras +1 more

The equation of state of a homogeneous two-dimensional Bose gas is calculated using quantum Monte Carlo methods. The low-density universal behavior is investigated using different interatomic model potentials, both finite ranged and strictly repulsive and zero ranged, supporting a bound state. The condensate fraction and the pair distribution function are calculated as a function of the gas parameter, ranging from the dilute to the strongly correlated regime. In the case of the zero-range pseudopotential we discuss the stability of the gaslike state for large values of the two-dimensional scattering length, and we calculate the critical density where the system becomes unstable against cluster formation.

arXiv:cond-mat/0411049•Physical Review A•Feb 11, 2005•51 citations
A path-integral quantum Monte Carlo method is used to calculate finite temperature properties of up to 1000 hard-core bosons in a two-dimensional isotropic harmonic-oscillator potential. If the interatomic repulsions are sufficiently short range, an abrupt increase in the condensate fraction and a hump in the specific heat occur close to the critical temperature of ideal bosons. The critical temperature and the condensate fraction are in general lowered by an increase in the hard-core radius a. If a is decreased below a certain level, the condensate fraction becomes indistinguishable from the corresponding value of the ideal bosons. For up to 1000 particles, this occurs when ln\ensuremath−1(1/na2)\ensuremath≲0.1,{\mathrm{ln}}^{\ensuremath{-}1}{(1/na}^{2})\ensuremath{\lesssim}0.1,ln\ensuremath−1(1/na2)\ensuremath≲0.1, where n is the average particle density.
Bose-Einstein condensation is known to occur in one-dimensional power-law potentials $V(x)\ensuremath{\propto}|x{|}^{\ensuremath{\eta}}$ as a true phase transition for \ensuremath{\eta}2, but only if there are no interactions between the particles. We show, via a path-integral quantum Monte Carlo scheme, that for both \ensuremath{\eta}2 and \ensuremath{\eta}=2 the spatial distribution of finite numbers of hard-core bosons suddenly becomes bimodal below a certain temperature, with a central condensate of particles distributed according to the lowest single-particle eigenstate. At still lower temperatures, the hard-core interactions cause the single-particle ground state to become a poor description for the interacting gas. If \ensuremath{\eta}2, the energy per particle undergoes a sudden decrease near the same temperature at which the bimodal distribution appears. It is found that this drop in energy disappears if the range of the interaction potential is sufficiently large.
The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas–Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.
Bose-Einstein condensation has been experimentally found to take place in finite trapped systems when one of the confining frequencies is increased until the gas becomes effectively two-dimensional (2D). We confirm the plausibility of this result by performing path-integral Monte Carlo (PIMC) simulations of trapped Bose gases of increasing anisotropy and comparing them to the predictions of finite-temperature many-body theory. PIMC simulations provide an essentially exact description of these systems; they yield the density profile directly and provide two different estimates for the condensate fraction. For the ideal gas, we find that the PIMC column density of the squeezed gas corresponds quite accurately to that of the exact analytic solution and, moreover, is well mimicked by the density of a 2D gas at the same temperature; the two estimates for the condensate fraction bracket the exact result. For the interacting case, we find 2D Hartree-Fock solutions whose density profiles coincide quite well with the PIMC column densities and whose predictions for the condensate fraction are again bracketed by the PIMC estimates.
We study the thermodynamic behavior of an ideal gas of bosons trapped in a three-dimensional anisotropic harmonic-oscillator potential. The condensate fraction as well as the specific heat is calculated using the Euler-Maclaurin approximation. For a finite number of particles there is no phase transition, but there is a well-defined temperature at which the condensation starts. We also consider condensation in lower dimensions, and, for one-dimensional systems, we discuss the dependence of the condensate fraction and heat capacity on the ensemble used.
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.
The equation of state of a homogeneous two-dimensional Bose gas is calculated using quantum Monte Carlo methods. The low-density universal behavior is investigated using different interatomic model potentials, both finite ranged and strictly repulsive and zero ranged, supporting a bound state. The condensate fraction and the pair distribution function are calculated as a function of the gas parameter, ranging from the dilute to the strongly correlated regime. In the case of the zero-range pseudopotential we discuss the stability of the gaslike state for large values of the two-dimensional scattering length, and we calculate the critical density where the system becomes unstable against cluster formation.