Boulangeot Nathan, Brix Florian, Sur Frédéric, Gaudry Émilie
Univ. de Lorraine, CNRS UMR7198, Institut Jean Lamour, Campus Artem, 2 allée André Guinier, 54000 Nancy, France.
Univ. de Lorraine, INRIA, CNRS UMR7503, Laboratoire Lorrain de Recherche en Informatique et Ses Applications, Campus Scientifique, 615 Rue du Jardin-Botanique, 54506 Vandœuvre-lès-Nancy, France.
J Chem Theory Comput. 2024 Aug 19. doi: 10.1021/acs.jctc.4c00367.
Intermetallic compounds are promising materials in numerous fields, especially those involving surface interactions, such as catalysis. A key factor to investigate their surface properties lies in adsorption energy maps, typically built using first-principles approaches. However, exploring the adsorption energy landscapes of intermetallic compounds can be cumbersome, usually requiring huge computational resources. In this work, we propose an efficient method to predict adsorption energies, based on a Machine Learning (ML) scheme fed by a few Density Functional Theory (DFT) estimates performed on sites selected through the Farthest Point Sampling (FPS) process. We detail its application on the AlCo(100) quasicrystalline approximant surface for several atomic adsorbates (H, O, and Pb). On this specific example, our approach is shown to outperform both simple interpolation strategies and the recent ML force field MACE [arXiv.2206.07697], especially when the number is small, i.e., below 36 sites. The ground-truth DFT adsorption energies are much more correlated with the predicted FPS-ML estimates (Pearson R-factor of 0.71, 0.73, and 0.90 for H, O and Pb, respectively, when = 36) than with interpolation-based or MACE-ML ones (Pearson R-factors of 0.43, 0.39, and 0.56 for H, O, and Pb, in the former case and 0.22, 0.35, and 0.63 in the latter case). The unbiased root-mean-square error (ubRMSE) is lower for FPS-ML than for interpolation-based and MACE-ML predictions (0.15, 0.17, and 0.17 eV, respectively, for hydrogen and 0.17, 0.25, and 0.22 eV for lead), except for oxygen (0.55, 0.47, and 0.46 eV) due to large surface relaxations in this case. We believe that these findings and the corresponding methodology can be extended to a wide range of systems, which will motivate the discovery of novel functional materials.