Sorry, I forgot to mention that the summation on k is restricted to first Brillouin zone only, in that case does ## \sum_{\mathbf{k} \in BZ} \rightarrow \frac{N\Omega}{(2\pi)^3} \int_{BZ} d\mathbf{k} ##
apply? I mean directly restricting the integration region to BZ.
I have some doubts on this...