Static lattice energy calculations have been combined with cluster expansion formalism to predict pairwise ordering interactions in the pyrope–grossular solid solution. The ordering interactions, the Js, have been then used to calculate the activity-composition relations over a wide temperature range with the help of the Cluster Variation Method. It is shown that short-range ordering in the system is driven by size mismatch. The prediction of the right signs and magnitudes of the ordering interaction energies requires separation of the mixing enthalpy into the configuration-dependent (chemical) and the configuration-independent (elastic) components. The study predicts the existence of a miscibility gap below 500°C.

The structural and thermodynamic properties of minerals are strongly affected by cation site-ordering processes. We describe methods to determine the main interatomic interactions that drive the ordering process, which are based on parameterizing model Hamiltonians using empirical interatomic potentials and/or ab initio quantum mechanics methods. The methods are illustrated by a number of case study examples, including Al/Si ordering in aluminosilicates, Mg/Ca ordering in garnets, simultaneous Al/Si and Mg/Al ordering in pyroxenes, micas and amphiboles, and Mg/Al non-convergent ordering in spinel using only quantum mechanical methods.

The procedure for deriving interatomic potentials for molecular ionic materials, using empirical fitting procedures, is described. Potentials are obtained for carbonates, phosphates and perchlorates, and used to calculate crystal and lattice properties which are compared with available experimental data.