 Translational and two-level energy variance

For a dilute classical ideal gas whose atoms carry independent nondegenerate internal levels $0,ϵ$, the canonical partition function factorizes into translation and $(1+e^{-\beta ϵ}{)}^N$. The upper-level probability is $p=(1+e^{\beta ϵ}{)}^{-1}$. Taking two derivatives of its logarithm gives the title sum of independent variances. The internal contribution freezes exponentially at low temperature and tends to $N{ϵ}^2/4$ at high temperature. The translational classical result cannot be continued through quantum degeneracy without changing the model.