Para fermiones $$\bar{n}_{FD}=\frac{1}{e^{(\epsilon -\mu)/kT}+1}$$ y $\epsilon$ puede ser más grande o más pequeño que el de $\mu$. Sin embargo, para los bosones: $$\bar{n}_{BE}=\frac{1}{e^{(\epsilon -\mu)/kT}-1}$$ lo que implica que $\epsilon >\mu$ de lo contrario, la ocupación es negativo.
Hay otro argumento de por qué $\epsilon >\mu$ debe ser cierto para los bosones, pero no para fermiones? Estoy leyendo Schroeder y dice:
Sería negativo si $\epsilon$ podría ser de menos de $\mu$, pero ya hemos visto que esto no puede suceder.
pero no puedo encontrar donde "ya hemos visto que esto no puede suceder"