Application of XDM to ionic solids: The importance of dispersion for bulk moduli and crystal geometries

J Chem Phys. 2020 Aug 7;153(5):054121. doi: 10.1063/5.0015133.

Abstract

Dispersion corrections are essential in the description of intermolecular interactions; however, dispersion-corrected functionals must also be transferrable to hard solids. The exchange-hole dipole moment (XDM) model has demonstrated excellent performance for non-covalent interactions. In this article, we examine its ability to describe the relative stability, geometry, and compressibility of simple ionic solids. For the specific cases of the cesium halides, XDM-corrected functionals correctly predict the energy ranking of the B1 and B2 forms, and a dispersion contribution is required to obtain this result. Furthermore, for the lattice constants of the 20 alkali halides, the performance of XDM-corrected functionals is excellent, provided that the base functional's exchange enhancement factor properly captures non-bonded repulsion. The mean absolute errors in lattice constants obtained with B86bPBE-XDM and B86bPBE-25X-XDM are 0.060 Å and 0.039 Å, respectively, suggesting that delocalization error also plays a minor role in these systems. Finally, we considered the calculation of bulk moduli for alkali halides and alkaline-earth oxides. Previous claims in the literature that simple generalized gradient approximations, such as PBE, can reliably predict experimental bulk moduli have benefited from large error cancellations between neglecting both dispersion and vibrational effects. If vibrational effects are taken into account, dispersion-corrected functionals are quite accurate (4 GPa-5 GPa average error), again, if non-bonded repulsion is correctly represented. Careful comparisons of the calculated bulk moduli with experimental data are needed to avoid systematic biases and misleading conclusions.