This view directly correlates with the work of Huang et al. where the authors found similar relations at high potentials, with the addition of constants involved in the formation of OH * . Therefore, the lack of this species in our simulation can be seen as an important limitation. Another limitation is the lack of direct kinetic aspect of our Monte Carlo simulations as experimental studies have shown oxidation to be time dependent. Instead, we are limited to the pure thermodynamic picture which is not likely representative of real conditions. As a result, the sharp oxidation peak visible in our framework might be attenuated or smoothed by kinetic effects under specific experimental conditions.