Predictive Thermodynamic Modells
Predictive Thermodynamic Modells
Predictive Thermodynamic Modells
ABSTRACT
INrnonucrroN
An understanding of the thermodynamics of interphase equilibria is
fundamental to the explanation of observed petrological phase rela-
tions. Moreover, if petrologists are ever to be able to predict phase
relations for unstudied combinations of components, our knowledge
of the thermodynamics of mixtures, of both solid and liquid solutions,
must be brought to a quantitative level. No matter how accurate or
completethe thermodynamic data for pure phasesbecomes,geochemical
calculationswill still be crucially limited by a lack of knowledge about
knowledge about mixtures, for most mineral systems are crystalline
solutions,and few of them are ideal.
If we divide the molar Gibbs energyof an m-componentmineral system
into its ideal and excessparts,
(where r; is the mole fraction of the itr' component in the mixture and
is the chemicalpotential of the pure component)the problem ma1-be
,l.r;0
formulated as a lack of quantitative or predictive knowledge of the
excessGibbs energyof mixing as a function of temperature,pressure,and
composition, G'"(7, P, x;, x.i ,. , , fi,-1).
If the excessfree energy of mixing can be calculatedfrom the proper-
1692
T H ERMODY N A M I C MOD IiLS OF ITA LI T E-SYLV IT ]J 1693
fo(x'o): exp[p""'*(roo)/RTl. Q)
Thirdly, the equilibrium distribution coefficient, which predicts
chemical partitioning between coexisting mineral pairs, ,4 and B, is
obtainable:
(xil.r;)-t - - (/r0"" -
* f (pt'" l.rj''")n lli.")ll
D= ^L
(*01*)" RT I
In addition, geothermometry and geobarometry would be specifiedby
the composition of the minerals present and a theoretical basis rvould
be provided for an understandingof non-equilibrium assemblages.
The presentation hereinafter will consist of the development of a
theoretical model appropriate to the s1-stemunder considerationcon-
taining certain parameters numerically unspecified by the theory;
the determination of the values of these parametersby examination of
the phase relations and distribution coefficientsin well-studied casesl
the correlation of the fitted values of the parameterswith the physico-
chemical properties of the pure components, and the prediction of
parametersfor unstudied systemsl and finally the prediction, by insdr-
tion of these assumedvalues into the theoretical model, of phase rela-
trons and thermodynamic properties of unstudied systemsof geological
interest. This paper will be concernedwith the first two of theseaspects.
A p p l i , c a t i o nt o H e t e r o p o l aS
r o l i d s . I n 1 9 3 5 B e t h e d e v i s e da m e t h o d f o r
deriving the equilibrium properties of a superlattice.Shortly thereafter
Rushbrooke (1938) appliedthe method to the study of mixtures of mole-
culessufficientlyalike in sizeand shapeto be interchangeableon a crystal
Iattice, but in which the configurationalenergy dependsupon the dispo-
sition of the moleculeson the lattice. SubsequentlyFowler and Guggen-
heim (1939) showed that the method was equivalent to a much more
elegantand direct method which has becomeknown as the quasi-chem-
ical (QC) treatment. The basic assumptionsunderlying this treatment
will be subsequentlyoutlined but detailed derivations must be sought
elsewhere.lIn the QC treatment it is assumedthat intermolecularforces
are central, short-ranged,and two-body additive, and that consequently
the internal energy of the system at the absolute zero oI temperature
may be evaluatedby summing the pair-potentialsover all nearestneigh-
bor pairs. It is further assumedthat the configurationaland vibrational
contributions to the partition function (or to the free energy of the mix-
ture) are separable,or equivalently,that the vibrational part of the heat
capacity is a linear function of the composition,as suggestedby -Joule's
(1844) rule.2Of these assumptionsperhaps the most difficult to justifl'
is the neglect of Iong range coulombic interactions when dealing with
essentiallyionic solids.
We believe, however, that although the Madelung term constitutes
about 90 percent of the lattice energy of a simple ionic crystal, the cou-
lombic contribution to the mixing energy may be neglectedin the first
approximation for isomorphoussolid solutions.The reasonsfor this are
twofold. First, becausethe formal valence chargeson substituting iso-
morphousspeciesare identical, the crystal structures and the Madelung
constantsare identical (or, in the caseof nonisometriccrystal systems,
nearll. linearly related). Thus the only differencein Madelung energy is
due to lattice parameter variations with composition.These diflerences
tend to be minimized due to lattice relaxation around substitutins ions.
1 In particular see Guggenheim (1952) in which the essential features are deveioped in
Chapter 4 and the asymmetrical model, subsequentlyused here, in Chapter 11. Throughout
this paper where equations are equivalent to those of Guggenheim, his equation numbers
are shown in square brackets.
'zWhich rule is sometimes invoked in the names of Neumann (1831), Woestyn (1848),
or Kopp (1865).
T H ERMODY N A M I C MO I)II.I-S OIt H A LI T E-S Y LV IT E
That is, the anions around a given substituting foreign cation tend to be
located at distancesfrom it determined b-vthe pair potential minimum
for that given ion pair.
Detailed direct calculationsby Douglas (1966)for very dilute crystal-
line alkali halide solutions,basedon the Born-Mayer (1932,1933)poten-
tial function, show that in theseface-centeredcubic systemsinterchange
of Na+ and K+ is accompaniedby an increasedenergy of about 5 kcal
per mole of replacing ion (although the lattice energiesof the pure salts
are of the order of -170 kcal mole t).1 Ol that 5 kcal mixing energy,by
far the most important contribution is from short range overiap repul-
sion. On the basis of less elaborate calculationsDurham and Hawkins
(1951) had earlier concludedthat the treatment of the specificrepulsion
constant was the factor having the greatesteffect on the heat of mixing
in alkali halide solid solutions.They attributed the successof their model
to allowancesmade for fluctuations in interionic distance. Still earlier
Wallace (1949)attributed the failure of his attempt to calculateheats of
mixing from the Born-Mayer model without lattice relaxation to inad-
equate representation of the replusive pctential. Wasastjerna (1944,
1949) has produced experimental evidencesuggestingthe existenceof
somelocal order in alkali halide solid solutions.He has developeda model
which ailows the ions to be displacedfrom averagelattice positions in
the mixed crystal dependingupon their sizeand upon the degreeof short
range order. With later modifications suggestedby Hovi (1950) the
model has been usedto calculateexcessenthaipy in alkali halide systems
with somedegreeof success.(c/. Lister anclMyers, 1958,for a comparison
with experiment.)The principal point of lailure of this model is the over-
estimationof locai order near the compositionof order (r2:0.5), and the
invariable prediction of negative excessnixing entropy.
In view of the relative importance of lattice configuration,relaxation,
and shorl rangerepulsion,as well as the relative unimportanceof coulom-
bic interaction in determining the mixing energy, the nearestneighbor
interactionmodel may not be as outrageousas it first appears.
p : Ir - 4e,6,[r- exp(2We/zRT)]l+,
(6),Lrr.07.r2l
where 2W6, the molar interchangeenergy, is a parameter of the theory
which measuresthe energy requirement to effect the interchange of a
[1] and a [2] atom betweentwo pure crystals.
An approximateform of the QC mixing expressionmay be derived for
the limiting caseof small interchangeenergy, Wc, ot high temperature.
In this case we note that the expression in square brackets in (6) is
small. Expanding B in binomial seriesyields
The first term on the right will be recognizedas a form of the van Laar
(1906, 1910)-Hildebrand(1929) expressionfor the enthalpy of mixing if
we take the contact factors, h/ tlz:Vrl Vr, to be propbrtional to the
molar volumes of the components.The secondterm on the right can be
identifiedwith Lumsden's (1952,p. 31S)correctionterm for nonrandom-
ness if we take the contact factors to be proportional to the $ power of
Lumsden's atomic radii q1fq2:(rrld+.Thus it is seen that the van
Laar-Hildebrand and Lumsden treatments may be consideredto be
limiting forms of the QC model.
For our purposes,however,it is no more difficult to work with the full
QC expression(4), and the approximate relations (8) will not be treated
further. The enthalpy of mixing may be obtainedfrom (4) using-a version
of the Gibbs-Helmholtzrelation,
d(G^ /RT\
: H'"' (9)
t1t/*)
However, beforedoing so it is usefulto de{rnetwo new quantities:
dWn
Wrt: Wc - T (roa),[4.27
.r]
dT
dll/ c
II/
-
vt s - - (10b)
lT'
.s"* (r2)
+:?r^[1
RT 2
+ o'q-l)-l (13b)
L 6,(0 1- r)J
t 698 E. TEWETT GREI|N
It is well known that the van Laar and Lumsden binary models predict
so it is not surprisingthat the QC model
a miscibiiitl'gap for positive l,tr/6,
does also. The consolute conditions of temperature and composition,
7", xr",r26,could, in principle,be evaluatedfrom the conditions,
({#)' : (1::;"-)"":
o, (14)
-0, ( 1s)
D i f f e r e n t i a t i n g( 1 3 a )w i t h r e s p e c t o r z a n d n o t i n g t h a t r r : 1 - r z in the
binarv system. we obtain
: o:;t- ,.#,::=,h1,
#(*iu)7. (17)
where the first term in the square brackets arisesfrom the ideal contri-
bution and the secondfrom the excesspotentiai of (13a). The spinodal
curve is thus defined by
lqrqrz
- 2(*rq, I xzqz))'
(1e)
This relation involves only qr, {2, and Z;and thus the critical consolute
compositionis definedin terms of theseparameters.The relation is, how-
ever, cubic in 11; and as a consequencethe simplest methods of solution
will be numerical rather than anall'tical. Having found the consolute
compositionit is a straightforward procedureto substitute this value jn
(18) to obtain the consolutetemperature.
I In which case the activitv coefficient of a comoonent at infinite dilution rvill not in
general be unity.
E. JEWETT GREDN
IJt.4 IJtB
(2La)
RT
Zo' l-
ln(rre) * j Inl t
1L
: rn(rlB)tf
+ ^lt (2rb)
lJz't
(2rc)
. Zo,
'-
In(rre) +
2
(21d)
Zo,
: ln(:rr,e) + -'-
2
t/qqz : r, (22)
than a constant. The mean value for this constant and its root-mean-
square(standard)deviation are,
q r / q z : 0 . 6 9 2+ 0 . 0 2 3 . (24)
Guggenheim (1944) takes the contact factors to be proportional to the
surfacearea of a moleculeinvolved in contactswith dissimilarmolecules.
Thus he showsfor simply branchedpolymer chainsof length 11,
2lr{Z-2)
qi: ( 2 s ) ,[ 1 1 . 0 1 . 1 ]
Vr/V2: 0.7199, ( 2 7a )
V12t\fV22t3:0.8033. (27b)
Temp.
References
E. 34
=o
>
32
U
z
U
?n
U
(,
z
I
o zd
E.
UJ
z
o
U
o
l
o
u
t
Frc. 1. Reduced interchange energy for NaCI-KCI crystalline solutions. Vertical bars
connect the two independent calculations. The solid line is the best linear weightedJeast-
squares fit. The dashed lines bound the 9501.,confidence limits (two standard deviations).
o O Borreffond Wolloce(1954)
tr Bunkond Ticheloor(,|953)
* Nocken(i918)
^ 450
C)
o
lrl
E
b 400
E
lrJ
(L
lrj
F --^
25o3
o.2 0.4 0.6 o.8 t.o
NoCl
MOLE FRACTION KCI
Frc. 2. Observed two-phase data and binodal curves calculated from QC model (solid line)
and SR model (dashed line) where the latter diverges from the former.
2RT" 2
'[ve r w h e nq t fq z : 1, (32),[4.12.131
f- z I
ZInl |
LZ _2J
which tends to unity as Z becomeslarge. For Z and qrf qzol the halite-
sylvite system we determinefrom Figure 4:
2RT"
:0.8s38. (33)
w;
Substitution of (30) into (33) yields
T" : 763.0"K. (34)
TH ERMODY NA M I C MO DELS OF II A LI T E-SYLTTI TII, 1705
2
F
(t',
(L
=
o
ul
F
J
o
o.
U)
z
o =6
o.20
Frc. 3. Consolute composition (mole fraction of component two in the binary system
at the maximum temperature at which two separate phases co-exist) as calculated from
the QC model for various coordination numbers.
The critical conditions determined, (31) and (34), are verv closeto those
yieldedby the SR model of Thompson and Waldbaum (1969),r2"--0.448,
T":765.2"K.
It is seen that both the QC and SR models represent the data well
(allhough the QC model does so with one less parameter) and that on
this basis alone there is little to recommend one model over the other.
Such is not the case, however, when the thermodynamic predictions of
the models are compared with the experimental calorimetric data.
I
=
,o
F
(r
N
;
E
f
k
lrJ
(L
trl
F
LU
F
f
<J)
z
()
t!
:l
o
trJ
E
'faer.n
2. Cer,oruMernrcDare ron NaCl,KCl Cnvsr,q.u;Nr Sor,urroNs
Reduced Heat
Mole Fraction Reduced Excess I{educed Excess Capacity between 500"
of KCl Dnthalpy, kcal mole-l Entropy, gibbs mole-t and 630"C,
x2 IIr" f arz S"uf rzxz gibbs moie-r
c;"f ''62
010 5 .1 6 0 .5 6
0.30 418 1.21
050 418 112
070 4.09 114
0 .90 4.47 -0 22
0.04 4 .5 8
0.11 4 .t 6
015 4.35
026 392
034 439
0.+2 42t
0.55 408
068
0.90 3. 7 1
09s 491
0.98 457
O L i s t e r o n d M e y e r s( , l 9 5 8 )
90
E l B o r r e t l o n d W o l l o c e( , | 9 5 4 )
i q)
o
E
380 -..QP
J
\ /L/n
N
x
x-
\
LI 7,O
I(L
J
r 6.0
F
z
Ld
(9
z. tr
i 5.o
= ku. o
o o
l!
O o tr tro tr
3+o tr
IJ
o
E.
80 E B o r r e t l o n d W o l l o c e( , | 9 5 4 )
; q,
o
E
4
o --- sR
o _._44OOr,
- --(
? e.o --
x
N
-\
; \ i
a\
v,
tE +.o
F
z
lrl
o
v,
lrl OC MODEL
9
lrJ
zo
o
trl
o
tr tr
f
o
lrl
tr
Eo
tr
CENTIGRADE T EI\4PERATURE
I
o
(I
O
F
U'
-
@
(D
UJ
X
Ld
R E C I P R O C AA
L E S O L U T ET E M P E R A T U R EI O
: TT
C o"" - (or-
2qgzxrrzlB2 6il,1
D
Z(*rqr+ nrq)(g+ 1)'B
(q#)' (3s)
for the case where dWH/dT:O. This function is shown in Figure 7 for
the equimolar solid as calculatedfrom (35) using (24) and (30). Between
500oand 600oCthe model predicts the equimolar excessheat capacity to
be about 0.44 gibbs mole-l. This is a substantially larger value than the
experimental 0.10 gibbs mole-l, but it should be noted that the heat
capacity is related to the secondtemperature derivative or curvature of
the Gibbs energysurfaceand somewhatlarger errorsin this function are
to be expected.
The values Douglas obtains are 6.04 and 4.09 kcal mole-1,respectively,
for il.2""and ilf". Again the agreementis quite remarkable.
C E N T I G R A D ET E M P E R A T U R E
2o,0
EXPERIMENTAL> -/.,i-,SA
tr' UOOtt
)/,/
o.o 05 10 15 2A 25 3.O
R E C I P R O C AA
L B S O L U T ET E M P E R A T U R EI O
: 3IT
Frc. 8. Reduced excess Gibbs energy of mixing for equimolar NaCl-KCl crystaliine
solutions. The solid line represents the QC model (this work), the dot-dashed line is the
prediction of the SR model (Thompson and Waldbaum, 1969), and the dashed iine repre-
sents the experimental observations (Barrett and Wallace; 1954a, b).
1712 E. IEWETT GREEN
BonN, M., arlo J E. Mnvrn (1932) Zur Gittertheorie der Ionenkristalle. Z. Physik,75,
1-18.
Bumr, A. J. H., lNr G. W. Trcnnee,n (1953) Investigations in the systems NaCI-KCL
Kon. Ned..Akad. Wetensch.Proc. 8.,56, 375-384.
Douores, T. B. (1966) Calculated heats of dilute solid solution among the alkali halides
other than cesium salts. .r. Chern. Ph.ys.4s, 4571 4585.
Dunrrau, G. S., ewo J. A. Hawrrus (1951) Solid solutions of the alkali halides II. The
theoretical calculation of lattice constants, heats of mixing, and distributions between
solid and aqueousphases.J. Chem.Phys.19, 149-156.
Folvr-rn, R. H., lxo E. A. GuccnNnrru (1939) Statistical thermodynamics of superlat-
tices. Proc. Roy. Soc. London Arl74r 189-206.
Grnes, J. W. (1873) A method of geometrical representation of the thermodynamic proper-
ties of substances by means of surfaces. Trans. Conn. Acad,. Arts Sci.21 382404.
TIIDRMODYNAMIC MODELS OI] HALITE-SYLVITE 1713