Adsorption and Stability of Monoatomic Adsorbate Adlayers on FCC and HCP Metals Using the Sphere-in-Contact Model

In this paper, we show that the sphere-in-contact model can predict long-range surface adsorption phenomena based on adsorbate-adsorbate repulsions and their geometric distance, assuming that their negative surface-induced charge is smeared on the surface of the adsorbate atoms. Additionally, it can be used to model collective surface diffusion mechanisms such as the domino-type surface diffusion of adsorbate rows on close-packed metal HCP and FCC surfaces. We have recently shown that the sphere-in-contact model can be used as an educational and research tool in various contexts, such as the visualization of carbon structures (e.g., graphene, carbon nanotubes, carbon nanocones, and graphite), heterogeneous catalysts, metal nanoparticles, and organic molecules. Here we present how it can be used to model the adsorbate structure of monoatomic elements on the hexagonal close-packed surface of HCP and FCC metals to study long-range ordering phenomena of monoatomic adsorbates on metals. We have used atoms of varying radius and color to represent the metal surface atoms and the adsorbate atoms. The study reveals that many surface configurations are possible for a fixed adsorbate coverage (θ) by the movement of the adsorbate atoms in response to surface adsorbate-adsorbate repulsions. The movement of the particles (e.g., particle diffusion) can be seen directly in the model, and this is caused by the user intervention. This has great educational and research value, as one can directly see how the adsorbate atoms reorder on the surface of a metal and therefore study diffusion mechanisms. We calculate the repulsive interaction energy of adsorbates using the sphere-in-contact model and can identify which surface-adsorbed configuration is the lowest energy. We find that at a surface coverage of 1/3 (0.333 ML), the most stable adsorbate configuration places adsorbates at the third nearest neighbor 3-fold hollow sites, forming a hexagonal pattern. We find that this model will be useful in the rational design of catalytic materials and material coatings with new technological applications where long-range ordering of surface adsorbates is essential and adsorbate interactions are mainly repulsive interatomic interactions.