This paper presents a study of the grand canonical Bose-Einstein (BE) statistics for a finite number of particles in an arbitrary quantum system. The thermodynamical quantities that identify BE condensation---namely, the fraction of particles in the ground state and the specific heat---are calculated here exactly in terms of temperature and fugacity. These calculations are complemented by a numerical calculation of fugacity in terms of the number of particles, without taking the thermodynamic limit. The main advantage of this approach is that it does not rely on approximations made in the vicinity of the usually defined critical temperature, rather it makes calculations with arbitrary precision possible, irrespective of temperature. Graphs for the calculated thermodynamical quantities are presented in comparison to the results previously obtained in the thermodynamic limit. In particular, it is observed that for the gas trapped in a three-dimensional box, the derivative of specific heat reaches smaller values than what was expected in the thermodynamic limit---here, this result is also verified with analytical calculations. This is an important result for understanding the role of the thermodynamic limit in phase transitions and makes possible to further study BE statistics without relying neither on the thermodynamic limit nor on approximations near critical temperature.
We have investigated the thermodynamic behaviour of ideal Bose gases with an arbitrary number of particles confined in a harmonic potential. By taking into account the conservation of the total number N of particles and using a saddle-point approximation, we derive analytically the simple explicit expression of mean occupation number in any state of the finite system. The temperature dependence of the chemical potential, specific heat, and condensate fraction for the trapped finite-size Bose system is obtained numerically. We compare our results with the usual treatment which is based on the grand canonical ensemble. It is shown that there exists a considerable difference between them at sufficiently low temperatures, especially for the relative small numbers of Bose atoms. The finite-size scaling at the transition temperature for the harmonically trapped systems is also discussed. We find that the scaled condensate fractions for various system sizes and temperatures collapse onto a single scaled form.
Bose-Einstein condensation (BEC) of a noninteracting Bose gas of N particles in a two-dimensional box with Dirichlet boundary conditions is studied. Confirming previous work, we find that BEC occurs at finite N at low temperatures T without the occurrence of a phase transition. The conventionally-defined transition temperature T(E) for an infinite three-dimensional (3D) system is shown to correspond in a 2D system with finite N to a crossover temperature between a slow and rapid increase in the fractional boson occupation N(0)/N of the ground state with decreasing T. We further show that T(E)∼1/logN at fixed area per boson, so in the thermodynamic limit there is no significant BEC in 2D at finite T. Thus, paradoxically, BEC only occurs in 2D at finite N with no phase transition associated with it. Calculations of thermodynamic properties versus T and area A are presented, including Helmholtz free energy, entropy S, pressure p, ratio of p to the energy density U/A, heat capacity at constant volume (area) C(V) and at constant pressure C(p), isothermal compressibility κ(T) and thermal expansion coefficient α(p), obtained using both the grand-canonical ensemble (GCE) and canonical ensemble (CE) formalisms. The GCE formalism gives acceptable predictions for S, p, p/(U/A), κ(T) and α(p) at large N, T and A but fails for smaller values of these three parameters for which BEC becomes significant, whereas the CE formalism gives accurate results for all thermodynamic properties of finite systems even at low T and/or A where BEC occurs.
Abstract This chapter describes the mechanism of Bose–Einstein condensation in the simplest ideal Bose gas case. An explicit expression for the critical temperature is derived in the analytically soluble case of a gas confined in a box. Various thermodynamic quantities are calculated above and below the critical temperature. Emphasis is given to the isothermal compressibility which diverges below the critical temperature. The behaviour of the off-diagonal one-body density matrix is discussed. Results for the fluctuations of the single-particle occupation numbers are derived in both the canonical and grand canonical ensembles.
We find the universal structure and scaling of the Bose-Einstein condensation (BEC) statistics and thermodynamics (Gibbs free energy, average energy, heat capacity) for a mesoscopic canonical-ensemble ideal gas in a trap with an arbitrary number of atoms, any volume, and any temperature, including the whole critical region. We identify a universal constraint-cutoff mechanism that makes BEC fluctuations strongly non-Gaussian and is responsible for all unusual critical phenomena of the BEC phase transition in the ideal gas. The main result is an analytical solution to the problem of critical phenomena. It is derived by, first, calculating analytically the universal probability distribution of the noncondensate occupation, or a Landau function, and then using it for the analytical calculation of the universal functions for the particular physical quantities via the exact formulas which express the constraint-cutoff mechanism. We find asymptotics of that analytical solution as well as its simple analytical approximations which describe the universal structure of the critical region in terms of the parabolic cylinder or confluent hypergeometric functions. The obtained results for the order parameter, all higher-order moments of BEC fluctuations, and thermodynamic quantities perfectly match the known asymptotics outside the critical region for both low and high temperature limits. We suggest two- and three-level trap models of BEC and find their exact solutions in terms of the cutoff negative binomial distribution (which tends to the cutoff gamma distribution in the continuous limit) and the confluent hypergeometric distribution, respectively. Also, we present an exactly solvable cutoff Gaussian model of BEC in a degenerate interacting gas. All these exact solutions confirm the universality and constraint-cutoff origin of the strongly non-Gaussian BEC statistics. We introduce a regular refinement scheme for the condensate statistics approximations on the basis of the infrared universality of higher-order cumulants and the method of superposition and show how to model BEC statistics in the actual traps. In particular, we find that the three-level trap model with matching the first four or five cumulants is enough to yield remarkably accurate results for all interesting quantities in the whole critical region. We derive an exact multinomial expansion for the noncondensate occupation probability distribution and find its high-temperature asymptotics (Poisson distribution) and corrections to it. Finally, we demonstrate that the critical exponents and a few known terms of the Taylor expansion of the universal functions, which were calculated previously from fitting the finite-size simulations within the phenomenological renormalization-group theory, can be easily obtained from the presented full analytical solutions for the mesoscopic BEC as certain approximations in the close vicinity of the critical point.
We investigate the thermodynamic geometry of an ideal quantum gas obeying inclusion statistics, characterized by a negative statistical parameter
g<0. In this framework the grand-canonical partition function admits a finite maximum fugacity, and the thermodynamic scalar curvature
R is strictly positive for all
g<0, reflecting effective attractive statistical interactions analogous to those of a bosonic system. As the fugacity approaches its maximum value,
R diverges, signaling a phase transition of the Bose-Einstein condensation type. A key distinction from the ordinary ideal Bose gas is that the condensation temperature is elevated relative to the bosonic case, and finite-temperature condensation occurs even in the dimensional regime
1/2<D/σ≤1 where standard bosons do not condense, while for
D/σ≤1/2 the transition temperature vanishes. Three independent criteria; divergence of
R, the maximum fugacity singularity, and the non-analytic cusp in the specific heat, coincide at the same condensation point, confirming the thermodynamic consistency of the transition.