Combining graph deep learning and London dispersion interatomic potentials: A case study on pnictogen chalcohalides

TL;DR

This study combines graph neural network-based interatomic potentials with semiempirical dispersion models to improve the accuracy of modeling layered pnictogen chalcohalides, especially in capturing van der Waals interactions.

Contribution

It introduces a straightforward method to incorporate dispersion corrections into graph deep learning potentials without retraining or refitting.

Findings

Dispersion-corrected potentials improve structural predictions.

Enhanced accuracy in van der Waals gap and layer thickness.

Method is broadly applicable to layered polar crystals.

Abstract

Machine-learning interatomic potential models based on graph neural network architectures have the potential to make atomistic materials modeling widely accessible due to their computational efficiency, scalability, and broad applicability. The training datasets for many such models are derived from density-functional theory calculations, typically using a semilocal exchange-correlation functional. As a result, long-range interactions such as London dispersion are often missing in these models. We investigate whether this missing component can be addressed by combining a graph deep learning potential with semiempirical dispersion models. We assess this combination by deriving the equations of state for layered pnictogen chalcohalides BiTeBr and BiTeI and performing crystal structure optimizations for a broader set of V-VI-VII compounds with various stoichiometries, many of which possess van der Waals gaps. We characterize the optimized crystal structures by calculating their X-ray diffraction patterns and radial distribution function histograms, which are also used to compute Earth mover's distances to quantify the dissimilarity between the optimized and corresponding experimental structures. We find that dispersion-corrected graph deep learning potentials generally (though not universally) provide a more realistic description of these compounds due to the inclusion of van der Waals attractions. In particular, their use results in systematic improvements in predicting not only the van der Waals gap but also the layer thickness in layered V-VI-VII compounds. Our results demonstrate that the combined potentials studied here, derived from a straightforward approach that neither requires fine-tuning the training nor refitting the potential parameters, can significantly improve the description of layered polar crystals.