K. Lejaeghere

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 Tm from first principles. An alternative is to use semi-empirical 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 DFT-PBE level for several semi-empirical approximations of αV and Tm , 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 Gruneisen parameter, provides a similar overall accuracy as exhaustive first-principles phonon calculations. For the melting temperature, expensive ab initio molecular-dynamics simulations still outperform semi-empirical methods.

Predictions of observable properties by density-functional theory calculations (DFT) are used increasingly often by experimental condensed-matter physicists and materials engineers as data. These predictions are used to analyze recent measurements, or to plan future experiments in a rational way. Increasingly more experimental scientists in these fields therefore face the natural question: what is the expected error for such a first-principles prediction? Information and experience about this question is implicitly available in the computational community, scattered over two decades of literature. The present review aims to summarize and quantify this implicit knowledge. This eventually leads to a practical protocol that allows any scientist -- experimental or theoretical -- to determine justifiable error estimates for many basic property predictions, without having to perform additional DFT calculations.

A central role is played by a large and diverse test set of crystalline solids, containing all ground-state elemental crystals (except most lanthanides). For several properties of each crystal, the difference between DFT results and experimental values is assessed. We discuss trends in these deviations and review explanations suggested in the literature.

A prerequisite for such an error analysis is that different implementations of the same first-principles formalism provide the same predictions. Therefore, the reproducibility of predictions across several mainstream methods and codes is discussed too. A quality factor Delta expresses the spread in predictions from two distinct DFT implementations by a single number. To compare the PAW method to the highly accurate APW+lo approach, a code assessment of VASP and GPAW (PAW) with respect to WIEN2k (APW+lo) yields Delta-values of 1.8 and 3.3 meV/atom, respectively. In both cases the PAW potentials recommended by the respective codes have been used. These differences are an order of magnitude smaller than the typical difference with experiment, and therefore predictions by APW+lo and PAW are for practical purposes identical.