In the B-E formula for the average number of identical bosons that can occupy a single particle state ψj , out of a set of such states
‹n›j = 1/{exp[β(ϵj – μj)] – 1}
there appears the chemical potential μj . Let me remind the definition of the chemical potential:
The chemical potential of a species in a mixture is the derivative of the free energy of the system with respect to the number of moles of that species, under constant temperature, and given that the concentrations of all the other species in the mixture remain constant.
But, in a gas of identical bosons, there is one single species. Also, since the particles are supposed not to interact, there are no reactions in the gas.
Bottom line, what has to do the number of particles occupying a certain orbital ψj , with chemical potential?
Could it be that the particles occupying different orbitals are considered as different species? But if so it is, how can this number change under a constant temperature? Constant temperature is a state of equilibrium, the average number of particle per orbital doesn't change.
(My question may be naïve, but it's a very long time since I learnt statistical mechanics.)