 Solution

In the large-$N$ Gross–Pitaevskii equation, dropping the kinetic energy gives $[\kappa r^2/2+gn(r)-\mu ]\Psi (r)=0$. The physical nonnegative Thomas–Fermi approximation for a condensate is thereforewhere $(z{)}_{+}=\max (z,0)$ and $R$ is the cloud radius. Its particle number isThus the two-dimensional Thomas–Fermi condensate obeysHere $N=\int n{d}^{2}r$ is a total particle number. Keeping the exact pair count in the fixed-$N$ Hartree product state would replace $g$ by $g(1-1/N)$ in these formulas expressed in terms of $n=N|\varphi {|}^2$.

The approximation is best at the centre, $r=0$, and worst near the edge, $r=R$. For example its omitted quantum pressure of a condensate, evaluated on the interior profile, isIt is smallest at the centre and diverges at the artificial sharp edge. The central healing length ${\xi }_0=\hslash /\sqrt{2m\mu }$ is much smaller than $R$ when $\mu \gg \hslash \omega$, since $R/{\xi }_0=2\mu /(\hslash \omega )$. This is the bulk validity condition; the full Gross–Pitaevskii equation smooths the edge and gives a decaying tail outside $R$.