 Solution

The occupation probability of an individual one-electron state is the Fermi-Dirac distribution,At $T=0$, it is a step: one below $E_F$ and zero above it. Heating rounds the step over an energy range of order a few $k_{B}T$, creating empty states below and occupied states above the chemical potential. At every positive temperature, $f(\mu ,T)=1/2$.

For fixed $N,V$ in the three-dimensional free electron model, the chemical potential is determined by $N=\int g(ϵ)f(ϵ,T)dϵ$; it is not exactly constant. In the degenerate regime $T\ll T_F$, the Sommerfeld expansion givesIt is therefore reasonable in a simple low-temperature sketch to set $\mu \simeq E_F$, provided this approximation is stated. The figure includes the leading small shift. At much higher temperature the fixed-$E_F$ approximation fails and the distribution approaches the classical dilute limit. The plotted occupation is $f$, rather than the energy distribution $gf$.