I have some doubt about the Bose-Einstein distribution for magnons/spin-waves.
A one-dimensional ferromagnet placed in an external magnetic field $\mathbf{B} = B\, \hat{z}$ obeys the Hamiltonian
$$H = - |J|\sum_i \mathbf{S}_i \cdot \mathbf{S}_{i+1} - B \sum_i {S}^z_i. ,$$ where we set the lattice constant $a = 1$.
By utilizing the Holstein-Primakoff transformation and the Fourier transform we can diagonalize this Hamiltonian and find the dispersion relation for the magnons $$\hbar \omega = 2 J S (1-\cos k) + B$$ where $k$ is the wavenumber of the magnon.
We know that magnons are bosons and should satisfy Bose-Einstein statistics. So the Bose-Einstein distribution for magnons is
$$n_k = \frac{1}{e^{\frac{\hbar \omega - \mu}{k_bT}}-1} = \frac{1}{e^{\frac{2 \hbar J S (1-\cos k) + B - \mu}{k_bT}}-1}.$$
However, shouldn't the distribution be normalized in the sense that if we integrate it over the Brilluin zone we should get unity $$\int_{1. BZ} n_k dk = 1?$$