A parameterized classical potential model for the interaction of water and hydroxide with iron oxide was used to calculate the optimal proton arrangement and proton binding energies on the (012) surface of hematite. Energy minimization calculations with the parameterized potential model indicate that approximately 75% of adsorbed water molecules are dissociated on this surface, in agreement with recent TPD and HREELS measurements. Surface protonation/deprotonation energies were calculated from the predicted optimal arrangement of protons on the neutral (012) surface. A supercell geometry with translational symmetry in two dimensions and finite in the third dimension (2-D PBC) was assumed. The calculated surface protonation energies were then used to model the experimentally observed surface-charging curve of hematite in aqueous solution. Excellent agreement was found between the calculated and measured surface charge for ionic strengths ranging from 0.001 to 0.1M. Our calculations favor the value of 8.5 for the pH of zero charge of hematite over the more recent result of 6.7.