The packing fraction of the oxygen sublattice: its impact on the heat of mixing

The heat of mixing of some petrological relevant substitutions (i.e., Mg-Al, Si-Al, Mg-Ti, Mg-Ca, and Mg-Fe) was investigated systematically in silicates, titanates, tungstates, carbonates, oxides, hydroxides, and sulphates by density functional theory calculations (e.g., melilite, chlorite, biotite, brucite, cordierite, amphibole, talc, pseudobrookite, pyroxene, olivine, wadsleyite, ilmenite, MgWO₄, ringwoodite (spinel), perovskite, pyrope-grossular, magnesite-calcite, MgO-CaO, anhydrous and different hydrated MgSO₄). A specific substitution is characterised by different microscopic interaction energies in different minerals, e.g., the octahedral Mg-Al exchange on a single crystallographic site in pyroxene has a microscopic interaction energy that is more than twice compared to that in biotite. A comparative investigation of the heat of mixing using microscopic interaction energies on a single crystallographic site has the advantage that they are not influenced by cation ordering. They could be successfully correlated with the stiffnesses of the minerals, which in turn were scaled to the oxygen packing fraction, a parameter that is easily available for poorly investigated minerals. With this information, the interaction energies of a certain substitution can be transferred from minerals where they are well-known to mineral groups where they are less- or unknown. Using the cross-site terms and the microscopic interaction energies, the macroscopic interaction energies of the coupled substitution, e.g., Mg + Si = Al + Al, of biotite and pyroxene were calculated, which are, however, affected by cation ordering and different degrees of local charge balance, for which appropriate models are necessary.