Abstract. The theory of sublattice solid solution model and optimization methods have been described for modelling the geochemically important multicomponentmultisite silicate solid solution systems. Some new X-ray Mg-Fe z § site occupancy data along with some selection from the existing data on heated orthopyroxene in the temperature range 600 to 1000~ C have been used in thermodynamic modelling of the orthopyroxene (Mg, Fe)2Si206 solid solution using the sublattice solution model. The optimized interaction energy solution parameters are: L~lg, ve = 13 600.7-4.92650 9 T L~2, Ve = 13 308.