Calculations using the density-functional theory (DFT) in combination with the single defect method were carried out to determine the heat of mixing behaviour of mineral solid solution phases. The accuracy of this method was tested on the halite–sylvite (NaCl–KCl) binary, pyrope–grossular garnets (Mg3Al2Si3O12–Ca3Al2Si3O12), MgO–CaO (halite structure) binary, and on Al/Si ordered alkali feldspars (NaAlSi3O8–KAlSi3O8); as members for coupled substitutions, the diopside–jadeite pyroxenes (CaMgSi2O6–NaAlSi2O6) and diopside–CaTs pyroxenes (CaMgSi2O6–CaAlAlSiO6) were chosen for testing and, as an application, the heat of mixing of the tremolite–glaucophane amphiboles (Ca2Mg5Si8O22(OH)2–Na2Mg3Al2Si8O22(OH)2) was computed. Six of these binaries were selected because of their experimentally well-known thermodynamic mixing behaviours. The comparison of the calculated heat of mixing data with calorimetric data showed good agreement for halite–sylvite, pyrope–grossular, and diopside–jadeite binaries and small differences for the Al/Si ordered alkali feldspar solid solution. In the case of the diopside–CaTs binary, the situation is more complex because CaTs is an endmember with disordered cation distributions. Good agreement with the experimental data could be, however, achieved assuming a reasonable disordered state. The calculated data for the Al/Si ordered alkali feldspars were applied to phase equilibrium calculations, i.e. calculating the Al/Si ordered alkali feldspar solvus. This solvus was then compared to the experimentally determined solvus finding good agreement. The solvus of the MgO–CaO binary was also constructed from DFT-based data and compared to the experimentally determined solvus, and the two were also in good agreement. Another application was the determination of the solvus in tremolite–glaucophane amphiboles (Ca2Mg5Si8O22(OH)2–Na2Mg3Al2Si8O22(OH)2). It was compared to solvi based on coexisting amphiboles found in eclogites and phase equilibrium experiments.
Read full abstract