Abstract:Sampling high-coverage configurations and predicting adsorbate-adsorbate interactions on surfaces are highly relevant to understand realistic interfaces in heterogeneous catalysis. However, the combinatorial explosion in the number of adsorbate configurations among diverse site environments presents a considerable challenge in accurately estimating these interactions. Here, we propose a strategy combining high-throughput simulation pipelines and a neural network-based model with the MACE architecture to increase sampling efficiency and speed. By training the models on unrelaxed structures and energies, which can be quickly obtained from single-point DFT calculations, we achieve excellent performance for both in-domain and out-of-domain predictions, including extrapolation to different facets, coverage regimes and low-energy configurations. From this systematic understanding of model robustness, we exhaustively sample the configuration phase space of catalytic systems without false positives or active learning. In particular, by evaluating over 14 million structures within the simulated annealing method, we predict coverage-dependent adsorption energies for CO adsorption on six Cu facets (111, 100, 211, 331, 410 and 711) and the co-adsorption of CO and CHOH on Rh(111). When validated by targeted post-sampling relaxations, our results for CO on Cu correctly reproduce experimental interaction energies reported in the literature, and provide atomistic insights on the site occupancy of steps and terraces for the six facets at all coverage regimes. Additionally, the arrangement of CO on the Rh(111) surface is demonstrated to substantially impact the activation barriers for the CHOH bond scission, illustrating the importance of comprehensive sampling on reaction kinetics. Our findings demonstrate that simplified data generation routines and controlled extrapolation of neural networks can be deployed at scale to understand lateral interactions on surfaces, paving the way towards realistic modeling of heterogeneous catalytic processes.

Accelerated calculation of configurational free energy using a combination of reverse Monte Carlo and neural network models: Adsorption isotherm for 2D square and triangular lattices

Thermodynamic calculations using reverse Monte Carlo: Simultaneously tuning multiple short-range order parameters for 2D lattice adsorption problem

An Efficient Multi-Scale Lattice Model Approach To Screening Nano-Porous Adsorbents

Importance Sampling within Configuration Space Integration for Adsorbate Thermophysical Properties: A Case Study for CH3/Ni(111)

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Accelerated structure-stability energy-free calculator

Accessing the Usefulness of Atomic Adsorption Configurations in Predicting the Adsorption Properties of Molecules with Machine Learning

An Artificial Neural Network-based Density Functional Approach for Adiabatic Energy Differences in Transition Metal Complexes

A New Approach to Model Adsorption in Heterogeneous Phase System with Monte Carlo Method

Towards an Ab-initio Description of Adsorbate Surface Dynamics

Predicting Electronic Structure Properties of Transition Metal Complexes with Neural Networks

Symmetry-Adapted High Dimensional Neural Network Representation of Electronic Friction Tensor of Adsorbates on Metals

Transferable Machine Learning Interatomic Potential for Carbon Hydrogen Systems

Controlling neural network extrapolation enables efficient and comprehensive sampling of coverage effects in catalysis

LAMMPS Implementation of Rapid Artificial Neural Network Derived Interatomic Potentials

Reweighting from Molecular Mechanics Force Fields to the ANI-2x Neural Network Potential

Adsorb-Agent: Autonomous Identification of Stable Adsorption Configurations via Large Language Model Agent

AdsorbML: A Leap in Efficiency for Adsorption Energy Calculations using Generalizable Machine Learning Potentials

Transferable prediction of formation energy across lattices of increasing size

Understanding the Origins of Reversible and Hysteretic Pathways of Adsorption Phase Transitions by Mesocanonical Ensemble Monte Carlo Simulations