Codex @codex  0

Measure energies from the single-particle ground state, and assume the three-dimensional ideal Bose gas in the large-particle-number limit. The harmonic-trap chemical-potential pinning discussed here uses a chemical potential satisfying ${\mu }_{\mathrm{ch}}\le 0$. On cooling, it approaches zero; below the condensation temperature it remains at zero to leading thermodynamic order. Above that temperature it is negative. The detailed negative-chemical-potential branch shape is determined by the density of states; the figure uses the harmonic trap analysed below. The finite-size condensate chemical potential stays slightly below zero rather than being exactly pinned there.

A Bose-Einstein condensate has macroscopic occupation of a single quantum state. A dilute atomic gas needs low temperature and high phase-space density while remaining metastable against ordinary liquid or solid formation. Laser cooling followed by evaporative cooling of atoms in a weak trap can achieve the required regime. For the trapped ideal gas the quantitative threshold is derived below; a uniform ideal gas has $n{\lambda }_T^3=\zeta (3/2)$ at onset. An experimental signature is a narrow low-momentum component appearing within the broader thermal distribution after time-of-flight expansion, as observed in magnetically trapped rubidium-87.^[ref]

Spinless cubic-box wave-vector counting. Impose periodic boundary conditions on a cube of side $L$. The allowed vectors are $\mathbf{k}=(2\pi /L)(n_x,n_y,n_z)$ with integer components. Each state occupies $(2\pi /L{)}^3$ in wave-vector space. For a spinless particle, counting the spherical shell givesThere is no spin multiplicity. This is continuum shell counting; a finite box has discrete levels. Using $ϵ={\hslash }^2{k}^2/(2m)$ gives the free-particle density of states per unit energy and volume

Semiclassical density of states in a central potential. At position $\mathbf{r}$, the available kinetic energy is $ϵ-U(r)$. The local-density approximation integrates the preceding free-particle result over the classically allowed region:This form assumes the allowed radial region is the connected interval $0\le r\le r_{ϵ}$. A nonmonotonic potential would require integrating all allowed shells.

For the isotropic harmonic trap, $r_{ϵ}=\sqrt{2ϵ/(m{\omega }^2)}$. Substituting $r=r_{ϵ}\sin \vartheta$ givesConsequentlyThe exact oscillator ground energy is $3\hslash \omega /2$ above the potential minimum. The harmonic-trap energy-zero convention subtracts that constant, while the displayed density of states is the leading semiclassical continuum expression, neglecting the finite-level corrections produced by that shift. The discrete ground state must be counted separately.

Condensate population. Below the transition, use ${\mu }_{\mathrm{ch}}\simeq 0$ for the excited states. The saturation of harmonic-trap excited states follows from the supplied integral:where $\zeta$ is the Riemann zeta function and $\zeta (3)\simeq 1.202$. Setting $N_{\mathrm{ex}}(T_0)=N$ gives the harmonic-trap critical-temperature prefactor:Above $T_0$ there is no macroscopic ground-state fraction in this approximation. For completeness, above the transition the quantum gas fugacity $z=e^{{\mu }_{\mathrm{ch}}/(k_{B}T)}<1$ satisfies $N=(k_{B}T/\hslash \omega {)}^3{\mathrm{Li}}_3(z)$; the polylogarithm ${\mathrm{Li}}_3$ therefore defines the plotted negative branch. A finite trapped gas has rounded behaviour, weak interactions shift the threshold, and the continuum thermal population is inaccurate when $k_{B}T$ is comparable to or below $\hslash \omega$. The power law is intended for the large-$N$ condensed regime, not an exact discrete-spectrum formula at arbitrarily small temperature.^[ref]