  Starting from Enskog equation of hard spheres of mass m and diameter D under
the gravity g, we first derive the exact equation of motion for the equilibrium
density profile at a temperature T and examine its solutions via the gradient
expansion. The solutions exist only when \beta\mu \le \mu_o \approx 21.756 in 2
dimensions and \mu_o\approx 15.299 in 3 dimensions, where \mu is the
dimensionless initial layer thickness and \beta=mgD/T. When this inequality
breaks down, a fraction of particles condense from the bottom up to the Fermi
surface.
