Low-temperature thermal phonon number ($N_{\mathrm{ph}}\sim 18\zeta (3)N(T/{\Theta }_D{)}^3$)

Integrating the Bose-Einstein distribution against the three-dimensional Debye model phonon density of states gives $N_{\mathrm{ph}}=9N(T/{\Theta }_D{)}^3{\int }_0^{{\Theta }_D/T}{x}^2/(e^x-1)dx$. For $T\ll {\Theta }_D$, expand $(e^x-1{)}^{-1}={\sum }_{j\ge 1}{e}^{-jx}$ and integrate termwise: the limiting integral is $2{\sum }_{j\ge 1}{j}^{-3}=2\zeta (3)$. This counts thermal quanta, not available modes or half quanta of zero-point energy.