Bose-Einstein condensation shown by Monte Carlo simulation

A Monte Carlo algorithm was constructed for calculating statistical physics quantities characterizing an ideal bosonic lattice gas. To demonstrate the possibilities of the algorithm the Bose-Einstein condensation was simulated within a three-dimensional isotropic harmonic oscillator where, in the frame of a canonical ensemble, average occupancy of the ground state and specific heat was calculated versus temperature and number of bosons. The algorithm can be further applied for studying both the static properties of ideal bosons within other trapping potentials and the relaxation of the system to the condensate. (C) 1999 Elsevier Science B.V. All rights reserved.