Hydrogen adsorption on nickel (100) single-crystal face. A Monte Carlo study of the equilibrium and kinetics

We propose a model of the dissociative adsorption of hydrogen on nickel single-crystal face. In this model, we treat the Ni(100) surface as a strongly correlated energetically heterogeneous surface, because the density functional theory (DFT) studies indicate that hydrogen atoms may adsorb either on hollow sites (energetically more favorable, binding energy 48 kJ/mol H) or bridge sites with the binding energy less by 11 kJ/mol H. The essential assumption of the proposed model is that the dissociation of the hydrogen molecule is possible only over the topmost Ni atom, and the resulting H atoms may adsorb either on two free hollow sites (but the adjacent bridge sites must be free) or two bridge sites (the adjacent hollow sites must be free). If the above condition is not fulfilled, then the dissociation and adsorption are impossible. The second assumption is that the rate (probability) of the associative desorption is limited by the rate of diffusion of H atoms on the surface. This is because the two H atoms desorb, giving an H2 molecule, only when they meet on two adjacent hollow-bridge sites. Our model recovers very well the behavior of the experimental equilibrium adsorption isotherms as well as kinetic isotherms. As a result, we stated that hydrogen atoms are not completely free on the surface, but they cannot also be considered localized at room and elevated temperatures. Additionally, while analyzing the kinetic adsorption isotherms, we stated that the rate-limiting step during the dissociative adsorption of H2 is the disintegration of the activated complex and the subsequent adsorption of hydrogen atoms.