Ab initio calculations often do not straightforwardly yield the thermal properties of a material yet. It requires considerable computational efforts, for example, to predict the volumetric thermal expansion coefficient α V or the melting temperature T m from first principles. An alternative is to use semiempirical approaches. They relate the experimental values to first-principles predictors via fits or approximative models. Before applying such methods, however, it is of paramount importance to be aware of the expected errors. We therefore quantify these errors at the density-functional theory level using the Perdew-Burke-Ernzerhof functional for several semiempirical approximations of α V and T m , and compare them to the errors from fully ab initio methods, which are computationally more intensive. We base our conclusions on a benchmark set of 71 ground-state elemental crystals. For the thermal expansion coefficient, it appears that simple quasiharmonic theory, in combination with different approximations to the Grüneisen parameter, provides a similar overall accuracy as exhaustive firstprinciples phonon calculations. For the melting temperature, expensive ab initio molecular-dynamics simulations still outperform semiempirical methods.