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

Chinese Physics

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

Zheng Rong-Jie, Jing Jin, Yi Ping Tang

Published Aug 22, 2006• 6 citations
Publisher

Abstract

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.

DOI
10.1088/1009-1963/15/9/010

Related papers

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

Bose-Einstein condensation in trapped bosons: A variational Monte Carlo analysis

Jonathan L. DuBois, H. R. Glyde

Several properties of trapped hard-sphere bosons are evaluated using variational Monte Carlo techniques. A trial wave function composed of a renormalized single-particle Gaussian and a hard-sphere Jastrow function for pair correlations is used to study the sensitivity of condensate and noncondensate properties to the hard-sphere radius and the number of particles. Special attention is given to diagonalizing the one-body density matrix 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 with the highest occupation. The effect of interaction on other quantities such as the ground-state energy, the mean radial displacement, and the momentum distribution is calculated as well. Results are compared with mean-field theory in the dilute limit.

arXiv:cond-mat/0008368•Physical Review A•Jan 5, 2001•83 citations

Natural orbitals and Bose-Einstein condensates in traps: A diffusion Monte Carlo analysis

Jonathan L. DuBois, H. R. Glyde

We investigate the properties of hard-core bosons in harmonic traps over a wide range of densities. Bose-Einstein condensation is formulated using the one-body density matrix (OBDM) which is equally valid at low and high densities. The OBDM is calculated using diffusion Monte Carlo methods and it is diagonalized to obtain the ``natural'' single-particle orbitals and their occupation, including the condensate fraction. At low boson density, ${\mathrm{na}}^{3}<{10}^{\ensuremath{-}5},$ where $n=N/V$ and a is the hard-core diameter, the condensate is localized at the center of the trap. As ${\mathrm{na}}^{3}$ increases, the condensate moves to the edges of the trap. At high density, it is localized at the edges of the trap. At ${\mathrm{na}}^{3}<~{10}^{\ensuremath{-}4},$ the Gross-Pitaevskii theory of the condensate describes the whole system within $1%.$ At ${\mathrm{na}}^{3}\ensuremath{\approx}{10}^{\ensuremath{-}3},$ corrections are $3%$ to the Grass-Pitaevskii energy but $30%$ to the Bogoliubov prediction of the condensate depletion. At ${\mathrm{na}}^{3}\ensuremath{\gtrsim}{10}^{\ensuremath{-}2},$ mean-field theory fails. At ${\mathrm{na}}^{3}\ensuremath{\gtrsim}0.1,$ the bosons behave more like a liquid ${}^{4}\mathrm{He}$ droplet than a trapped boson gas.

arXiv:cond-mat/0303426•Physical Review A•Sep 10, 2003•39 citations

Condensate depletion in two-species Bose gases: A variational quantum Monte Carlo study

Asaad R. Sakhel, Jonathan L. DuBois, H. R. Glyde

We investigate two-species Bose gases in traps with various interactions using variational quantum Monte Carlo (VMC) techniques at zero temperature. The bosons are represented by hard spheres whose diameter is equivalent to the $s$-wave scattering length in the low-energy and long-wavelength approximation. We explore the role of repulsive and attractive interspecies or intraspecies interactions on the condensate properties of the mixtures, particularly the condensate fraction of each species as compared to the case when each species is in a separate trap of its own. We model the repulsive interactions by a hard-core (HC) potential and the attractive interactions by a shallow model potential. The VMC density profiles and energies are evaluated at various interactions and two mass ratios of the species.

arXiv:cond-mat/0703248•Physical Review A•Apr 24, 2008•14 citations

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

Sean Pearson, Tao Pang, Changfeng Chen

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.

Physical Review A•Dec 1, 1998•15 citations

Ground-State Properties of Charged Bosons Confined in a One-Dimensional Harmonic Double-Well Trap: Diffusion Monte Carlo Calculations

Jing Jin, Tang Yi

The diffusion Monte Carlo method is applied to study the ground-state properties of charged bosons in one dimension confined in a harmonic double-well trap. The particles interact repulsively through a Coulombic 1/ r potential. Numerical results show that the well separation has significant influence on the ground-state properties of the system. When the interaction of the system is weak, ground-state energy decreases with the increasing well separation and has a minimal value. If the well separation increases continually, the ground-state energy increases and approaches to a constant gradually. This effect will be abatable in the strong interacting system. In addition, by calculating the density of the systems for different interaction strengths with various well separations, we find that the density increases abnormally when the well separation is large at the centre of the system.

Chinese Physics Letters•Aug 23, 2007•1 citation
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.
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.
Several properties of trapped hard-sphere bosons are evaluated using variational Monte Carlo techniques. A trial wave function composed of a renormalized single-particle Gaussian and a hard-sphere Jastrow function for pair correlations is used to study the sensitivity of condensate and noncondensate properties to the hard-sphere radius and the number of particles. Special attention is given to diagonalizing the one-body density matrix 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 with the highest occupation. The effect of interaction on other quantities such as the ground-state energy, the mean radial displacement, and the momentum distribution is calculated as well. Results are compared with mean-field theory in the dilute limit.
We investigate the properties of hard-core bosons in harmonic traps over a wide range of densities. Bose-Einstein condensation is formulated using the one-body density matrix (OBDM) which is equally valid at low and high densities. The OBDM is calculated using diffusion Monte Carlo methods and it is diagonalized to obtain the ``natural'' single-particle orbitals and their occupation, including the condensate fraction. At low boson density, ${\mathrm{na}}^{3}<{10}^{\ensuremath{-}5},$ where n=N/Vn=N/Vn=N/V and a is the hard-core diameter, the condensate is localized at the center of the trap. As na3{\mathrm{na}}^{3}na3 increases, the condensate moves to the edges of the trap. At high density, it is localized at the edges of the trap. At ${\mathrm{na}}^{3}<~{10}^{\ensuremath{-}4},$ the Gross-Pitaevskii theory of the condensate describes the whole system within 11%.1 At ${\mathrm{na}}^{3}\ensuremath{\approx}{10}^{\ensuremath{-}3},$ corrections are 33%3 to the Grass-Pitaevskii energy but 3030%30 to the Bogoliubov prediction of the condensate depletion. At ${\mathrm{na}}^{3}\ensuremath{\gtrsim}{10}^{\ensuremath{-}2},$ mean-field theory fails. At ${\mathrm{na}}^{3}\ensuremath{\gtrsim}0.1,$ the bosons behave more like a liquid 4He{}^{4}\mathrm{He}4He droplet than a trapped boson gas.
We investigate two-species Bose gases in traps with various interactions using variational quantum Monte Carlo (VMC) techniques at zero temperature. The bosons are represented by hard spheres whose diameter is equivalent to the sss-wave scattering length in the low-energy and long-wavelength approximation. We explore the role of repulsive and attractive interspecies or intraspecies interactions on the condensate properties of the mixtures, particularly the condensate fraction of each species as compared to the case when each species is in a separate trap of its own. We model the repulsive interactions by a hard-core (HC) potential and the attractive interactions by a shallow model potential. The VMC density profiles and energies are evaluated at various interactions and two mass ratios of the species.
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.
The diffusion Monte Carlo method is applied to study the ground-state properties of charged bosons in one dimension confined in a harmonic double-well trap. The particles interact repulsively through a Coulombic 1/ r potential. Numerical results show that the well separation has significant influence on the ground-state properties of the system. When the interaction of the system is weak, ground-state energy decreases with the increasing well separation and has a minimal value. If the well separation increases continually, the ground-state energy increases and approaches to a constant gradually. This effect will be abatable in the strong interacting system. In addition, by calculating the density of the systems for different interaction strengths with various well separations, we find that the density increases abnormally when the well separation is large at the centre of the system.