We substitute this value for I, then find å = 6, A = 0.5092, and B = 0.3283 in Table 3.1, and obtain a value for the activity coefficient of 0.8237, and an activity of 0.308 ×10−3 m. If we did the calculation for other temperatures, we would see that for a dilute solution such as this, the activity coefficient is only a weak function of temperature, decreasing to 0.625 at 300° C.
3.8 IDEAL SOLUTIONS IN CRYSTALLINE SOLIDS AND THEIR ACTIVITIES
When we deal with solid solutions, we are again faced with the inadequacy of the purely macroscopic approach of classical thermodynamics. There is little disadvantage to this approach for gases, where the arrangement of molecules is chaotic. But the crystalline state differs from that of gases in that the arrangement of atoms in the crystal lattice is highly ordered, and the properties of the crystal depend strongly on the nature of the ordering. For this reason, we cannot afford to ignore the arrangement of atoms in solids, particularly with respect to solutions.
Solid solutions differ from those of gases and liquids in several respects. First, solution in the solid state inevitably involves substitution. While we can increase the concentration of HCl in water simply by adding HCl gas, we can only increase the concentration of Fe in biotite solid solution if we simultaneously remove Mg. Second, solid solutions involve substitution at crystallographically distinct sites. Thus, in biotite a solid solution between phlogopite (KMg3AlSi3O10(OH)2) and annite (KFe3AlSi3O10(OH)2) occurs as Fe2+ replaces Mg2+ in the octahedral site; the tetrahedral Si site and the anion (O) sites remain unaffected by this substitution. Third, substitution is often coupled. For example, the solid solution between anorthite (CaAl2Si2O8) and albite (NaAlSi3O8) in plagioclase feldspar involves not only the substitution of Na+ for Ca2+, but also the substitution of Al3+ for Si4+. The anorthite–albite solution problem is clearly simplified if we choose anorthite and albite as our components rather than Na+, Ca2+, Al3+ and Si4+. Such components are known as phase components. Choosing pure phase end members as components is not always satisfactory either because substitution on more than one site is possible, leading to an unreasonably large number of components, or because the pure phase does not exist and hence its thermodynamic properties cannot be measured.
However we choose our components, we need a method of calculating activities that takes account of the ordered nature of the crystalline state. Here we will discuss two ideal solution models of crystalline solids. We tackle the problem of nonideal solid solutions in Chapter 4.
3.8.1 Mixing-on-site model
Many crystalline solids can be successfully treated as ideal solutions. Where this is possible, the thermodynamic treatment and assessment of equilibrium are greatly simplified. A simple and often successful model that assumes ideality but takes account of the ordered nature of the crystalline state is the mixing-on-site model, which considers the substitution of species in sites individually. In this model, the activity of an individual species is calculated as:
(3.78)
where X is the mole fraction of the ith atom and ν is the number of sites per formula unit on which mixing takes place. For example, ν = 2 in the Fe–Mg exchange in olivine, (Mg,Fe)2SiO4. One trick to simplifying this equation is to pick the formula unit such that ν = 1. For example, we would pick (Mg,Fe)Si½O2 as the formula unit for olivine. We must then consistently choose all other thermodynamic parameters to be half those of (Mg,Fe)2SiO4.
The entropy of mixing is given by:
(3.79)
where the subscript j refers to sites and the subscript i refers to components, and n is the number of sites per formula unit. The entropy of mixing is the same as the configurational entropy, residual entropy, or “third law entropy” (i.e., entropy when T = 0 K). For example, in clinopyroxene, there are two exchangeable sites, a sixfold-coordinated M1 site, (Mg, Fe2+, Fe3+, Al3+), and an eightfold-coordinated M2 site (Ca2+, Na+). Here j ranges from 1 to 2 (e.g., 1 = M1, 2 = M2), but n = 1 in both cases (because both sites accept only one atom). i must range over all present ions in each site, so in this example, i ranges from 1 to 4 (1 = Mg, 2 = Fe2+, etc.) when j = 1 and from 1 to 2 when j = 2. Since we have assumed an ideal solution, ΔH = 0 and ΔGideal = −TΔS. In other words, all we need is temperature and eqn. 3.79 to calculate the free energy of solution.
In the mixing-on-site model, the activity of a phase component in a solution, for example, pyrope in garnet, is the product of the activity of the individual species in each site in the phase:
(3.80)
where aφ is the activity of phase component φ, i are the ion components of pure φ, and νi is the stoichiometric proportion of i in pure φ. For example, to calculate the activity of aegirine (NaFe3+Si2O6) in aegirine-augite ([Na,Ca][Fe3+,Fe2+,Mg]Si2O6), we would calculate the product: XNa
A slight complication arises when more than one ion occupies a structural site in the pure phase. For example, suppose we wish to calculate the activity of phlogopite (KMg3Si3AlO10(OH)2) in a biotite of composition K0.8Ca0.2(Mg0.17Fe0.83)3Si2.8Al1.2O10(OH)2. The tetrahedral site is occupied by Si and Al in the ratio of 3:1 in the pure phase end members. If we were to calculate the activity of phlogopite in pure phlogopite using eqn. 3.80, the activities in the tetrahedral site would contribute only