Temperature dependency of equilibrium at metal oxide-aqueous electrolyte solution interface was analyzed by numerical simulation. Derivations of inner surface potential with respect to temperature were performed at constant values of several different parameters. When surface charge density in inner plane was kept constant the reasonable results were obtained, i.e. the electrostatic contribution to enthalpy of protonation of amphotheric surface sites was found to be positive in the pH region below the point of zero potential and negative above this point. All other examined possibilities produced opposite results. Derivation of empirical interfacial equilibrium constant at constant surface potential indicated that electrostatic effect on protonation entropy is negligible and that electrostatic contributions to reaction Gibbs energy and enthalpy are equal and directly related to the surface potential in the inner plane. (c) 2005 Elsevier B.V. All rights reserved.