The phase-space distribution of semi-infinite nuclear matter is expanded in an ħ series analogous to the low-temperature expansion of the Fermi function. Besides the usual Wigner-Kirkwood expansion, oscillatory terms are derived. In the case of a Woods-Saxon potential, a smallness parameter is defined, which determines the convergence of the series and explains the very rapid convergence of the Wigner-Kirkwood expansion for average (nuclear) binding energies.