 Classical grand potential in a harmonic trap

For noninteracting particles in a three-dimensional isotropic quantum harmonic oscillator, measure the one-particle ground energy as zero. The single-particle partition function is $Z_1=(1-e^{-\beta \hslash \omega }{)}^{-3}$, from the product of three geometric series. In the dilute Maxwell–Boltzmann statistics regime, $\ln \Xi =e^{\beta \mu }{Z}_1$, so the grand potential is the displayed expression and $N=e^{\beta \mu }{Z}_1$. For $k_{B}T\gg \hslash \omega$,The same approximation follows from integrating the density of states $g(E)\simeq E^2/[2(\hslash \omega {)}^3]$ against $e^{-\beta E}$. Extrapolating the classical formula until the fugacity is one gives the scale of Bose condensation in a three-dimensional harmonic trap, but the quantitative transition temperature requires the Bose occupation law rather than Maxwell-Boltzmann statistics.