Occupancy Statistics and Entropy in Bose Systems

In this work, we compare three formulations of thermodynamic entropy for a non-interacting bosonic gas: (i) the grand-canonical Bose-Einstein entropy, (ii) the finite-N canonical entropy obtained from the exact partition function (Ziff, Uhlenbeck, and Kac construction), and (iii) the entropy of a multinomial distribution with Boltzmann categorical probabilities and temperature determined from Clausius' equation. It is well-known that the grand-canonical Bose-Einstein systematically overestimates entropy of canonical systems in regimes where particle-number fluctuations are significant. The exact canonical entropy correctly enforces the particle-number constraint, but recent experimental results suggest that it also overestimates particle-number fluctuations below the crossover temperature. The multinomial distribution is less common in thermodynamics. It addresses in a mathematically exact way the puzzle of the famous -log(N!) term introduced by Gibbs as a deus ex machina in discussions of thermodynamic entropy. One remarkable consequence is that the entropy of the multinomial distribution overcomes the issue of negative entropy at low temperature that affects the Gibbs and the Sackur-Tetrode entropies at low temperature. The analysis presented in the paper shows that the multinomial distribution, equipped with a categorical distribution calibrated in such a way that the resulting multinomial entropy fits Clausius' equation, provides accurate approximations to the canonical entropy in the classical regime, while it is smaller than the canonical entropy below the crossover temperature. One feature of the multinomial distribution is that it predicts lower peak variance of the number of particles in the ground state than the canonical distribution. This is in agreement with recent experimental results; hence, this paper identifies the thermodynamically calibrated multinomial distribution as a candidate alternative to the canonical distribution for thermodynamic bosonic entropy in finite systems.