Blackbody ionizing photon function ($F(y)={\int }_y^{\infty }{x}^2/(e^x-1)dx$)

A spherical blackbody star emits $Q(E_0)=8{\pi }^2{R}_s^2(k{T}_s{)}^{3}F(E_0/k{T}_s)/(h^3{c}^2)$ photons per time above threshold $E_0$. Integrating specific intensity over the outward hemisphere gives flux $\pi B_{\nu }$; dividing the resulting spectral luminosity by photon energy and substituting $x=h\nu /k{T}_s$ proves the formula. The convergent series $F(y)={\sum }_{m\ge 1}{e}^{-my}(y^2/m+2y/m^2+2/m^3)$ follows by integrating the geometric-series expansion termwise.