Simulation of forsterite protonation by the interatomic potential method