Role of the surface heterogeneity in adsorption of hydrogen ions on metal oxides: theory and simulations

In this article we study the effect of energetic heterogeneity of a crystalline surface on the adsorption of hydrogen ions (protons) from the liquid phase. In particular, we examine the influence of the shape of the adsorption energy distribution on the equilibrium isotherms of hydrogen ions. To that purpose, a few popular distribution functions, including rectangular, exponential, and asymmetric Gaussian are considered. Additionally, multimodal distribution functions, which may correspond to the adsorption on different crystal planes of the oxide, are also used. Lateral interactions between adsorbed charges are modeled using the potential function proposed by Borkovec et al., which accounts also for polarization of the liquid medium. The results presented here are obtained using both Monte Carlo (MC) simulations and theoretical calculations involving Mean Field Approximation (MFA). They indicate that increased energetic heterogeneity of the adsorbing surface may, in general, considerably change the behavior of the adsorption isotherms, regardless of the assumed distribution function. It is also shown that the predictions of the proposed theory are consistent with the data obtained from the MC simulations.