A kinetic Monte Carlo/UBI-QEP model of O2 adsorption on fcc (1 1 1) metal surfaces

The previously developed kinetic Monte Carlo model of molecular oxygen adsorption on fcc (1 0 0) metal surfaces has been extended to fcc (1 1 1) surfaces. The model treats uniformly all elementary steps of the process—O₂ adsorption, dissociation, recombination, desorption, and atomic oxygen hopping—at various coverages and temperatures. The model employs the unity bond index—quadratic exponential potential (UBI-QEP) formalism to calculate coverage-dependent energetics (atomic and molecular binding energies and activation barriers of elementary steps) and a Metropolis-type algorithm including the Arrhenius-type reaction rates to calculate coverage- and temperature-dependent features, particularly the adsorbate distribution over the surface. Optimal values of non-energetic model parameters (the spatial constraint, a travel distance of “hot” atoms, attempt frequencies of elementary steps) have been chosen. Proper modifications of the fcc (1 0 0) model have been made to reflect structural differences in the fcc (1 1 1) surface, in particular the presence of two different hollow sites (fcc and hcp). Detailed simulations were performed for molecular oxygen adsorption on Ni(1 1 1). We found that at very low coverages, only O₂ adsorption and dissociation were effective, while O₂ desorption and O₂ and O diffusion practically did not occur. At a certain O + O₂ coverage, the O₂ dissociation becomes the fastest process with a rate one-two orders of magnitude higher than adsorption. Dissociation continuously slows down due to an increase in the activation energy of dissociation and due to the exhaustion of free sites. The binding energies of both molecular and atomic oxygen decrease with coverage, and this leads to greater mobility of atomic oxygen and more pronounced desorption of molecular oxygen. Saturation is observed when the number of adsorbed molecules becomes approximately equal to the number of desorbed molecules. Simulated coverage dependences of the sticking probability and of the atomic binding energy are in reasonable agreement with experimental data. From comparison with the results of the previous work, it appears that the binding energy profiles for Ni(1 1 1) and Ni(1 0 0) have similar shapes, although at any coverage the absolute values of the oxygen binding energy are higher for the (1 0 0) surface. For metals other than Ni, particularly Pt, the model projections were found to be too parameter-dependent and therefore less certain. In such cases further model developments are needed, and we briefly comment on this situation.