We question the validity of the grand canonical ensemble for the description of Bose-Einstein condensation of small ideal Bose gas samples in isolated harmonic traps. While the ground state fraction and the specific heat capacity can be well approximated with the help of the conventional grand canonical arguments, the calculation of the fluctuation of the number of particles contained in the condensate requires a microcanonical approach. Resorting to the theory of restricted partitions of integer numbers, we present analytical and numerical results for such fluctuations in one- and three-dimensional traps, and show that their magnitude is essentially independent of the total particle number.
The choice of statistical ensemble has a measurable impact on the statistical description of Bose-Einstein condensates beyond the well-known grand canonical anomaly. We develop a segmented-contour approach for the ideal Bose gas, which refines the saddle-point approximation by including non-Gaussian corrections in the canonical and microcanonical ensembles. The method applies to a broad class of noninteracting spectra, and provides efficient access to condensate fluctuations, higher moments, and correlation functions. For a three-dimensional harmonic trap with fixed anisotropy, we find that the ratio of the maxima of the microcanonical and canonical condensate-number variances approaches the theoretically predicted asymptotic value
1−3ζ(3)2/[4ζ(2)ζ(4)]≈0.391 only extremely slowly with increasing particle number. We also show that ensemble effects are visible in the second-order correlation function
g(2)(x) near the condensation crossover. After ballistic release from a three-dimensional harmonic trap, the time evolution of
g(2)(x) preserves the differences between statistical ensembles.
Grand-canonical fluctuations of Bose-Einstein condensates of light are accessible to state-of-the-art experiments [J. Schmitt et al., Phys. Rev. Lett. 112, 030401 (2014).PRLTAO0031-900710.1103/PhysRevLett.112.030401]. We phenomenologically describe these fluctuations by using the grand-canonical ensemble for a weakly interacting Bose gas at thermal equilibrium. For a two-dimensional harmonic trap, we use two models for which the canonical partition functions of the weakly interacting Bose gas are given by exact recurrence relations. We find that the grand-canonical condensate fluctuations for weakly interacting Bose gases vanish at zero temperature, thus behaving qualitatively similarly to an ideal gas in the canonical ensemble (or microcanonical ensemble) rather than the grand-canonical ensemble. For low but finite temperatures, the fluctuations remain considerably higher than for the canonical ensemble, as predicted by the ideal gas in the grand-canonical ensemble, thus clearly showing that we are not in a regime in which the ensembles are equivalent.
For ideal Bose atoms in an isotropic harmonic trap, we consider thermodynamic variables obtained from microcanonical, canonical, and grand canonical ensembles, each with certain variables specified and other variables fluctuating. For the first two of these ensembles, we derive recursion relations that link partition functions for different dimensions. We discuss fluctuations in general, and obtain expressions for variances of the atom number N, the chemical potential \ensuremath{\mu}, and the temperature T for small, Gaussian fluctuations in the grand canonical ensemble. Then from our recursion relations and others given elsewhere, we obtain probability distributions for N, for ground-state occupation number
n0, for \ensuremath{\mu}, and for T. Below the critical temperature, the shape of the distributions for N,
n0, and \ensuremath{\mu} are definitely not Gaussian for the grand canonical ensemble. For given temperature and small N, we find that the chemical potential values pertaining to the three ensembles differ. We compare the specific-heat function
CN(T) for the three ensembles and propose to use the minimum of
dCN/dT to define the critical temperature to facilitate comparisons with similar configurations of interacting atoms.
Abstract In two-dimensional traps, since the theoretical study of Bose–Einstein condensation will encounter the problem of divergence, the actual contribution of the divergent terms is often estimated in some indirect ways with the accuracy to the leading order. In this paper, by using an analytical continuation method to solve the divergence problem, we obtain the analytical expressions of critical temperature and condensate fraction for Bose gases in a 2D anisotropic box and harmonic trap, respectively. They are consistent with or better than previous studies. Then, we further consider the nonvanishing chemical potential, and obtain the expressions of chemical potential and more precise condensate fraction. These results agree with the numerical calculation well, especially for the case of harmonic traps. The comparison between the grand canonical and canonical ensembles shows that our calculation in the grand canonical ensemble is reliable.
We analyze the general relation between canonical and grand canonical ensembles in the thermodynamic limit. We begin our discussion by deriving, with an alternative approach, some standard results first obtained by Kac and coworkers in the late 1970s. Then, motivated by the Bose-Einstein condensation (BEC) of trapped gases with a fixed number of atoms, which is well described by the canonical ensemble and by the recent groundbreaking experimental realization of BEC with photons in a dye-filled optical microcavity under genuine grand canonical conditions, we apply our formalism to a system of non-interacting Bose particles confined in a two-dimensional harmonic trap. We discuss in detail the mathematical origin of the inequivalence of ensembles observed in the condensed phase, giving place to the so-called grand canonical catastrophe of density fluctuations. We also provide explicit analytical expressions for the internal energy and specific heat and compare them with available experimental data. For these quantities, we show the equivalence of ensembles in the thermodynamic limit.