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 consider an interacting, dilute Bose gas trapped in a harmonic potential at a positive temperature. The system is analyzed in a combination of a thermodynamic and a Gross–Pitaevskii (GP) limit where the trap frequency
ω , the temperature T, and the particle number N are related by
N∼(T/ω)3→∞ while the scattering length is so small that the interaction energy per particle around the center of the trap is of the same order of magnitude as the spectral gap in the trap. We prove that the difference between the canonical free energy of the interacting gas and the one of the noninteracting system can be obtained by minimizing the GP energy functional. We also prove Bose–Einstein condensation in the following sense: The one-particle density matrix of any approximate minimizer of the canonical free energy functional is to leading order given by that of the noninteracting gas but with the free condensate wavefunction replaced by the GP minimizer.
Current experiments on atomic gases in highly anisotropic traps present the opportunity to study in detail the low temperature phases of two-dimensional inhomogeneous systems. Although, in an ideal gas, the trapping potential favors Bose-Einstein condensation at finite temperature, interactions tend to destabilize the condensate, leading to a superfluid Kosterlitz-Thouless-Berezinskii phase with a finite superfluid mass density but no long-range order, as in homogeneous fluids. The transition in homogeneous systems is conveniently described in terms of dissociation of topological defects (vortex-antivortex pairs). However, trapped two-dimensional gases are more directly approached by generalizing the microscopic theory of the homogeneous gas. In this paper, we first derive, via a diagrammatic expansion, the scaling structure near the phase transition in a homogeneous system, and then study the effects of a trapping potential in the local density approximation. We find that a weakly interacting trapped gas undergoes a Kosterlitz-Thouless-Berezinskii transition from the normal state at a temperature slightly below the Bose-Einstein transition temperature of the ideal gas. The characteristic finite superfluid mass density of a homogeneous system just below the transition becomes strongly suppressed in a trapped gas.
Bose-Einstein condensation (BEC) of an ideal gas is investigated, beyond the thermodynamic limit, for a finite number
N of particles trapped in a generic three-dimensional power-law potential. We derive an analytical expression for the condensation temperature
Tc in terms of a power series in ${x}_{0}={\ensuremath{\varepsilon}}_{0}/{k}_{B}{T}_{c}$, where ${\ensuremath{\varepsilon}}_{0}$ denotes the zero-point energy of the trapping potential. This expression, which applies in Cartesian, cylindrical, and spherical power-law traps, is given analytically at infinite order. It is also given numerically for specific potential shapes as an expansion in powers of
x0 up to the second order. We show that, for a harmonic trap, the well-known first-order shift of the critical temperature $\ensuremath{\Delta}{T}_{c}/{T}_{c}\ensuremath{\propto}{N}^{\ensuremath{-}1/3}$ is inaccurate when $N\ensuremath{\leqslant}{10}^{5}$, the next order (proportional to ${N}^{\ensuremath{-}1/2}$) being significant. We also show that finite-size effects on the condensation temperature cancel out in a cubic trapping potential, e.g., $V(r)\ensuremath{\propto}{r}^{3}$. Finally, we show that in a generic power-law potential of higher order, e.g., $V(r)\ensuremath{\propto}{r}^{\ensuremath{\alpha}}$ with $\ensuremath{\alpha}>3$, the shift of the critical temperature becomes positive. This effect provides a large increase of
Tc for relatively small atom numbers. For instance, an increase of about $+$40% is expected with
104 atoms in a $V(r)\ensuremath{\propto}{r}^{12}$ trapping potential.