The reversed experimental determination of Al content of aluminous tremolites (Jenkins, 1994) provides powerful constraints on the mixing properties of amphiboles in this system. Jenkins compared several ideal mixing models with his reversals and noted that random mixing on both octahedral and tetrahedral sites was not supported by his experiments, and that a much more ordered situation was implied. Several such models, including those used by Jenkins, have been fitted to his experiments and the one which best fits the data is a partialorder model outlined below. Although the coupled onesite and twosite models of Jenkins yield close fits to his quartzsaturated experiments, they did not fit the silicaundersaturated data so well (Jenkins, 1994). There is also a logical inconsistency in using such models because they imply an ordered tschermakite endmember, yet the fits to the experiments demand a tschermakite entropy which is 4050 J/K higher than estimates for ordered tschermakite. To resolve this, the new partialorder model assumes that Al and Mg mix on the M2 octahedral sites, but that there is less than maximal configurational entropy on the tetrahedral sites. Restricting the mixing of Al and Si to just the 4 T1 sites still yields too high an entropy to account satisfactorily for Jenkins' reversals, as he noted. Thus it will be assumed that the entropy contribution is half that for random mixing on the T1 sites, a situation which can be considered to correspond to chargebalanced coupling between an Al on M2 and one or other of the two nearest T1 site Al cations, rather than strict charge balance coupling to just one single tetrahedral neighbour. The activity of tremolite and tschermakite endmembers in this model are just those of the random mixing model but with the square root of the tetrahedral site terms.
where y = XAlM2 is the mole fraction, or proportion, of the tschermakite endmember. The experimental data of Jenkins (1994) and Hoschek (1995) were used to determine the nonideal mixing enthalpy, using a simple regular solution model
The square root of the activities at 1000K are shown below, with a value for Wtrts = 20 kJ/mol.
Coded up for THERMOCALC:
%========= Tremolitetschermakite amphiboles in CMASH =============== hb 2 tr 1/4 2 1 1 1 1 2 2 1 1 1 2 gamma 1 20 0 0 1 1 1 1 0 1 1 1 2 % W(trts) ts 1 2 0 1 1 1 3 2 1 1 1 1 gamma 1 20 0 0 1 1 1 1 1 1 1 1 2 % W(trts) bulk 1 1 0 1 1 1 y(hb) 0.08
The experimental reversals of Jenkins (1994) involving both silicasaturated and undersaturated assemblages as well as those of Hoschek (1995) are well fitted, as may be seen in the figures below, with this reducedentropy mixing model. The derived enthalpy of the tschermakite endmember depends on this mixing model.
