Quadratic low-temperature energy correction for constant density of states

For constant density of states, weighting the difference between the Fermi-Dirac distribution and the zero-temperature step by $E-E_F$ produces a convergent integral that scales as ${\beta }^{-2}$. Hence the leading low-temperature energy correction is proportional to $T^2$.