Grand canonical Monte Carlo simulations are carried out for the basic Stem model of the electrical double layer formed at the energetically heterogeneous metal oxide/electrolyte interface. The effect of the global (macroscopic) and local (microscopic) adsorption energies correlations as well as the influence of the model parameter on the surface charge density curves were investigated. The linear dependence of point of zero charge (PZC) as a function of H+ ion adsorption energy proves that the acidic/basic properties of the system are mainly governed by proton uptake/release. Two kinds of systems were taken into account: one neglecting lateral interactions and the other one including them. The effect of electrolyte concentrations as well as the surface heterogeneity on the surface charge density curves were shown too. The presented simulation algorithm allows to model two experimentally observed instances of the metal oxide/electrolyte interface: one possessing a common intersection point (CIP) at pH = PZC and the other one with CIP not equal PZC. (c) 2006 Elsevier Inc. All rights reserved.