= Solution
For a bosonic degree of freedom, the supplied integral gives
$$
I_B=\int_0^\infty\frac{y^3}{e^y-1}\,dy=\frac{\pi^4}{15}.
$$
For a fermionic degree of freedom, use
$$
\frac1{e^y+1}=\frac1{e^y-1}-\frac2{e^{2y}-1}.
$$
Changing variable $u=2y$ in the second integral gives
$$
I_F=I_B-\frac18I_B=\frac78I_B.
$$
This is the energy-density case of the <fermion-to-boson thermal integral ratio>.
Summing the result of part (a) over all particle species therefore gives
$$
\epsilon
=\frac{4\pi g_*}{(hc)^3}(kT)^4\frac{\pi^4}{15}
=\boxed{\frac{4\pi^5}{15(hc)^3}g_*(kT)^4},
$$
where the <effective number of relativistic energy degrees of freedom> is
$$
\boxed{g_*=\sum_{\rm bosons}g_i+\frac78\sum_{\rm fermions}g_i.}
$$
Back to article page
